Characterizing the arrhythmogenic substrate in personalized models of atrial fibrillation: sensitivity to mesh...

Computationally guided persistent atrial fibrillation (PsAF) ablation has emerged as an alternative to conventional treatment planning. To make this approach scalable, computational cost and the time required to conduct simulations must be minimized while maintaining predictive accuracy. Here, we assess the sensitivity of the process to finite-element mesh resolution. We also compare methods for pacing site distribution used to evaluate inducibility arrhythmia sustained by re-entrant drivers (RDs).

Simulations were conducted in low- and high-resolution models (average edge lengths: 400/350 µm) reconstructed from PsAF patients' late gadolinium enhancement magnetic resonance imaging scans. Pacing was simulated from 80 sites to assess RD inducibility. When pacing from the same site led to different outcomes in low-/high-resolution models, we characterized divergence dynamics by analysing dissimilarity index over time. Pacing site selection schemes prioritizing even spatial distribution and proximity to fibrotic tissue were evaluated. There were no RD sites observed in low-resolution models but not high-resolution models, or vice versa. Dissimilarity index analysis suggested that differences in simulation outcome arising from differences in discretization were the result of isolated conduction block incidents in one model but not the other; this never led to RD sites unique to one mesh resolution. Pacing site selection based on fibrosis proximity led to the best observed trade-off between number of stimulation locations and predictive accuracy.

Simulations conducted in meshes with 400 µm average edge length and ∼40 pacing sites proximal to fibrosis are sufficient to reveal the most comprehensive possible list of RD sites, given feasibility constraints.