After localized ablation of a spiral-wave source, atrial fibrillation (AF) terminates stochastically: in planar 2D simulations with an idealized disk-shaped heterogeneity, and in a small clinical cohort of early (paroxysmal) AF patients, the number of spiral-wave tips behaves as a birth-death process whose measured transition rates W_r(n), quasi-stationary tip-number distribution P_qs(n), and mean termination time tau characterize termination (Rappel et al., 2022). Whether this framework survives anatomically realistic, patient-specific 3D atrial geometry with wall thickness, fiber architecture, and complex lesion shapes; how atrial remodeling and dynamical anchoring to inexcitable obstacles modify the birth and death rates and tau; and whether the predicted termination-time statistics hold in larger clinical cohorts including persistent AF remain unresolved. The problem is to establish the domain of validity and quantitative form of the stochastic termination framework in realistic atrial substrate and across the clinical spectrum of AF, and to test its clinical predictions.

Background

Atrial fibrillation (AF), the most common clinical cardiac arrhythmia, is widely modeled as spatially extended chaos of electrical spiral waves in excitable cardiac tissue. A spiral wave is a rotating pattern of electrical activation whose tip is a phase singularity, i.e. a point where the activation phase is undefined; in AF, spiral waves continuously break up, so tips are repeatedly created (birth) and pairwise annihilated (death), a state called spiral defect chaos. One influential clinical hypothesis holds that AF is maintained by localized spiral-wave sources anchored to tissue heterogeneities, motivating radiofrequency catheter ablation (deliberate creation of inexcitable lesion scars) targeted at those sources. Rappel, Krummen, Baykaner, Zaman, Donsky, Swarup, Miller and Narayan (Frontiers in Network Physiology, 2022) combined computer simulations of two standard excitable-media models, the Luo-Rudy and Fenton-Karma ionic models (systems of coupled equations for the cardiac cell's membrane voltage and ion-channel currents), with intracardiac mapping in AF patients. In their 2D simulations, a stable spiral wave was trapped by a disk-shaped heterogeneity and surrounded by breakup; after virtual ablation made the source region inexcitable, the residual spiral-wave activity terminated stochastically, with Poisson-distributed termination times. The dynamics were accurately captured by a master equation for the probability P(n,t)P(n,t) of having nn spiral-wave tips, driven by empirical birth-death transition rates Wr(n)W_r(n) with steps r=±1r=\pm1 and r=±2r=\pm2; the resulting quasi-stationary distribution Pqs(n)P_{qs}(n) (the long-lived metastable distribution of tip number before extinction, in the sense of extinction of metastable stochastic populations) and the mean termination time τ\tau agreed with direct simulation. Notably, the ±1\pm1 rates collapsed onto a single curve when rescaled by the total length of non-conducting boundary, and the ±2\pm2 rates scaled with domain area. Clinically, targeted source ablation terminated AF but only after variable delays up to minutes, without preceding spatiotemporal organization, consistent with the stochastic picture. Related stochastic-population descriptions had earlier been established: Vidmar and Rappel (Physical Review E, 2019) showed with a birth-death 'extinction dynamics' approach that mean spiral-defect-chaos lifetimes can be computed for arbitrary geometries and depend exponentially on tissue size, and Dharmaprani et al. (Circulation: Arrhythmia and Electrophysiology, 2019) showed that phase-singularity regeneration times in simulations, human AF, and animal models are exponentially distributed, consistent with renewal-type statistics. The 2022 study was, however, restricted to flat 2D domains with idealized heterogeneity geometry, and to a modest cohort of early (paroxysmal, i.e. self-terminating within about a week) AF patients with limited lesion sets; it also left open how atrial remodeling (the progressive electrophysiological change, fibrosis deposition, and conduction slowing that accompany AF) and dynamical anchoring of waves to obstacles modify the birth-death rates and τ\tau. These limitations define the open problem addressed in this task.

Problem Statement

Establish whether the stochastic birth-death description of spiral-wave termination after source removal, demonstrated by Rappel et al. (2022) in planar 2D domains with an idealized disk-shaped heterogeneity, extends to realistic atrial substrate and across the clinical spectrum of atrial fibrillation. Concretely: (1) Geometry -- determine whether, and in what modified form, the measured transition rates Wr(n)W_r(n) (steps r=±1r=\pm1, ±2\pm2), the observed scalings of the ±1\pm1 rates with total non-conducting boundary length and of the ±2\pm2 rates with tissue area, the quasi-stationary tip-number distribution Pqs(n)P_{qs}(n), and the mean termination time τ\tau persist in anatomically realistic, patient-specific 3D atrial geometries that include wall thickness, fiber architecture, and complex ablation-lesion shapes, reading the tip number nn and the boundary-length and area scalings through their explicit three-dimensional counterparts (for example, tip filaments or their surface intersections; non-conducting boundary surface area or lesion perimeter; tissue volume or mapped surface area). (2) Substrate -- determine whether and how atrial remodeling (electrophysiological remodeling, fibrosis progression, conduction slowing) and dynamical anchoring of waves to inexcitable obstacles change the birth and death rates and hence τ\tau. (3) Clinical -- test, in cohorts larger and more heterogeneous than the original paroxysmal-AF (self-terminating) series and including persistent AF (AF sustained beyond about seven days), whether termination after targeted source ablation follows the predicted stochastic (exponential) termination-time statistics, and whether the model's inferred rates predict termination behavior across patient subgroups. These three components share one scientific objective: fixing the domain of validity and the quantitative form of the stochastic termination framework for atrial fibrillation in realistic substrate. A resolution requires evidence bearing on all three components: a demonstrated failure of the framework in any one of them is itself a refutation of the claimed generalization, whereas a positive result restricted to a single component does not by itself establish the generalization.

References

  1. Wouter-Jan Rappel, David E. Krummen, Tina Baykaner, Junaid Zaman, Alan Donsky, Vijay Swarup, John M. Miller, Sanjiv M. Narayan (2022). Stochastic Termination of Spiral Wave Dynamics in Cardiac Tissue. Frontiers in Network Physiology 2, 809532. DOI: 10.3389/fnetp.2022.809532. https://doi.org/10.3389/fnetp.2022.809532
  2. David Vidmar, Wouter-Jan Rappel (2019). Extinction dynamics of spiral defect chaos. Physical Review E 99, 012407. DOI: 10.1103/PhysRevE.99.012407. https://doi.org/10.1103/PhysRevE.99.012407 (preprint: arXiv:1808.01991, titled 'Extinction Dynamics of Cardiac Fibrillation')
  3. Dhani Dharmaprani, Madeline Schopp, Pawel Kuklik, Darius Chapman, Anandaroop Lahiri, Lukah Dykes, Feng Xiong, Martin Aguilar, Benjamin Strauss, Lewis Mitchell, Kenneth Pope, Christian Meyer, Stephan Willems, Fadi G. Akar, Stanley Nattel, Andrew D. McGavigan, Anand N. Ganesan (2019). Renewal Theory as a Universal Quantitative Framework to Characterize Phase Singularity Regeneration in Mammalian Cardiac Fibrillation. Circulation: Arrhythmia and Electrophysiology 12, e007569. DOI: 10.1161/CIRCEP.119.007569. https://doi.org/10.1161/CIRCEP.119.007569
  4. Timothy J. Tyree, Patrick Murphy, Wouter-Jan Rappel (2024). Annihilation dynamics during spiral defect chaos revealed by particle models. Chaos 34, 053131. DOI: 10.1063/5.0203319. https://doi.org/10.1063/5.0203319 (arXiv:2402.10308)
  5. Timothy J. Tyree, Michael Reiss, Wouter-Jan Rappel (2025). Particle model of creation and annihilation captures termination dynamics of spiral defect chaos. Chaos 35, 051103. DOI: 10.1063/5.0277113. https://doi.org/10.1063/5.0277113
  6. Mahesh Kumar Mulimani, Wouter-Jan Rappel (2025). Spiral defect chaos with intermittency increases mean termination time. arXiv:2505.06427. https://arxiv.org/abs/2505.06427
  7. Melvin Dix, Thomas Lilienkamp, Stefan Luther, Ulrich Parlitz (2024). Influence of conduction heterogeneities on transient spatiotemporal chaos in cardiac excitable media. Physical Review E 110, 044207. DOI: 10.1103/PhysRevE.110.044207. https://doi.org/10.1103/PhysRevE.110.044207
  8. Wouter-Jan Rappel, Tina Baykaner, Junaid Zaman, Prasanth Ganesan, Albert J. Rogers, Sanjiv M. Narayan (2024). Spatially Conserved Spiral Wave Activity During Human Atrial Fibrillation. Circulation: Arrhythmia and Electrophysiology 17, e012041. DOI: 10.1161/CIRCEP.123.012041. https://doi.org/10.1161/CIRCEP.123.012041
  9. Junaid A. B. Zaman, Abdulhaseeb Khan, Jan Nielsen, Steen B. Kristiansen, Mads B. Kronborg, Christoffer T. Witt, Christian Gerdes, Jens Kristensen, Henrik K. Jensen, Peter Lukac, Sharad C. Agarwal (2024). Repeat ablation of atrial fibrillation using electrogram dispersion to identify additional areas of mechanistic significance. Heart Rhythm O2 5, 543-550. DOI: 10.1016/j.hroo.2024.07.007. https://doi.org/10.1016/j.hroo.2024.07.007
  10. Evan V. Jenkins, Dhani Dharmaprani, Madeline Schopp, Jing Xian Quah, Kathryn Tiver, Lewis Mitchell, Martyn P. Nash, Richard H. Clayton, Kenneth Pope, Anand N. Ganesan (2023). Markov modeling of phase singularity interaction effects in human atrial and ventricular fibrillation. Chaos 33, 063136. DOI: 10.1063/5.0141890. https://doi.org/10.1063/5.0141890
  11. Thomas Lilienkamp, Ulrich Parlitz (2020). Terminating transient chaos in spatially extended systems. Chaos 30, 051108. DOI: 10.1063/5.0011506. https://doi.org/10.1063/5.0011506