Independent Mesh Realizations Introduce Percent-Level Variability in Temporal Interference Simulations

Computational models of temporal interference stimulation (TIS) commonly report a single electric-field estimate for a given anatomy and electrode montage. Because non-deterministic tetrahedral mesh generation does not produce a unique discretisation of a fixed tissue-label image, a single mesh realisation may introduce numerical variability. We quantified variation across independent mesh realisations and contrasted it with repeated downstream simulation execution on a single selected mesh. Ten head models were evaluated for stimulation of the left hippocampus and right primary motor cortex (M1). For every model and target, we generated 40 independent meshes and performed one complete simulation on each. Separately, we selected the mesh whose parcel-level field estimate was closest to the median and repeated downstream operations 40 times while holding that geometry fixed, yielding 1,600 TIS simulations in total. The primary outcome was the spatial median of the TIS envelope field within a spherical target region. Across independently remeshed runs, within-participant coefficients of variation were 1.81-3.65% for the hippocampus and 1.62-2.79% for M1. Repeated execution on a fixed mesh reduced run-to-run standard deviation by more than 99%, demonstrating that workflow variability is driven almost entirely by non-deterministic mesh generation rather than solver instability, numerical rounding, or post-processing. Single-run mesh realisations preserved overall cohort ordering (median Kendall's τ of 0.867 for the hippocampus and 0.911 for M1) but frequently inverted the rank order of participant pairs with similar predicted fields. Furthermore, a bootstrap analysis demonstrated that averaging five to ten independent remesh runs effectively suppressed this stochastic noise. These results quantify single-workflow repeatability rather than absolute error. Stochastic mesh variation should therefore be controlled or mitigated through multi-run averaging whenever experimental conclusions depend on subtle field differences or fixed neuromodulation thresholds.