{"id":"1ca65cd2-f58e-4b83-9a72-e7e6e183fe6b","arxiv_id":"1908.01646","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Persistent atrial fibrillation in the CMP cellular automaton model is caused by complex re-entrant circuits with an asymmetry between activation and termination rates.","lead":"Atrial fibrillation in a simple computer model shows two patterns: short self-stopping episodes or long persistent ones, even when the simulated heart tissue is set up identically. The paper finds this difference comes from special re-entrant circuits that are easier to switch on than off, producing episode lengths from seconds to months.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cCMP control removes not only complex critical structures but also most conduction-blocking nodes; the CMP–cCMP gap therefore conflates complex re-entry with changes in total block density and wavefront interactions.","rationale":"The central claim is that the entire spectrum of AF persistence in the CMP model stems from complex critical structures with asymmetric activation/deactivation probabilities. The cleanest support for that claim is the CMP-vs-cCMP comparison (Figs. 8 and 12). The load-bearing assumption is that removing the random conduction-blocking nodes outside simple critical segments changes nothing except the presence of those complex structures. That assumption is insecure. The construction in Section IV A changes not only the spatial distribution of block nodes but also their total number: cCMP contains roughly δ×M block nodes, while CMP contains δ×L². A system with fewer conduction-blocking nodes has different wavefront fragmentation, different refractory interactions, and different suppression of nearby circuits even if the number of simple critical structures is matched. The paper's own admission that no systematic detection of complex critical structures was performed (Section IV A) further weakens the causal attribution; the two mechanisms in Figs. 9–10 were identified by inspection of selected configurations. The collapse with decreasing δ in Fig. 12 is suggestive but not decisive, since the difference in absolute block-node count also decreases with δ. The concrete ablation test proposed here would settle the issue: if removing only multi-block clusters from an otherwise intact CMP lattice reproduces the full gap, the mechanism is confirmed; if not, the gap is a density artefact. Because the existing reader verdict is already CONDITIONAL and identifies essentially the same weakest assumption, my read does not move the verdict; it sharpens the check that would move it to ACCEPT or REJECT.","tokens_in":30650,"tokens_out":9579,"duration_ms":105537,"concrete_test":"Perform a surgical ablation control: for each CMP realization at ν⊥ = 0.11 and δ = 0.01 (plus ν⊥ in {0.09, 0.10, 0.12} and δ in {0.05, 0.005}), identify every cluster containing two or more block nodes on the same isolated run of consecutive nodes of length ≥ τ/2, and delete all but one block node per such cluster. This removes the multi-block configurations needed for complex critical structures while leaving total block-node density almost unchanged. Simulate the ablated lattice for S = 10⁶ time steps, 200 realizations, and compare P(AF) and time in AF against the original CMP and cCMP. If the ablated CMP collapses onto cCMP, the gap is due to multi-block complex structures; if it remains near CMP, the original comparison is confounded by total block density and non-structure wavefront interactions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV A constructs cCMP by placing block-susceptible nodes only on isolated segments of length ≥ τ/2, at fraction δ of those eligible nodes, while the CMP comparison places block nodes at probability δ on all L² nodes. Thus CMP has approximately δL² block nodes whereas cCMP has only approximately δM with M (eligible segment nodes) much smaller than L². The two models therefore differ not only in the presence of multi-block complex structures but also in total block-node density, in the refractory landscape, and in all wavefront-collision effects caused by the deleted random nodes. The inference in Section IV B that 'same number of simple critical structures, therefore higher-order critical structures must exist' is a non sequitur unless total density and non-structure effects are controlled. The collapse of CMP and cCMP in Fig. 12 at low δ does not settle this, because the absolute density difference also shrinks. If added block density alone reproduces the CMP–cCMP gap, then the central claim that asymmetric complex critical structures drive persistent AF is not established by this comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes the Christensen-Manani-Peters (CMP) cellular-automaton model of atrial fibrillation to explain why different simulations at identical parameters produce AF episodes ranging from seconds to months. The authors first derive a mean-field (MF) birth-death model in which the number of 'simple' re-entrant circuits, counted from the lattice, determines AF statistics, and they show analytically that this MF model underestimates time in AF and cannot produce persistent AF for realistic N. They then introduce a 'controlled' CMP (cCMP) model in which conduction-block-susceptible nodes are placed only on isolated segments that can support simple re-entrant circuits, and they compare this with a random-placement CMP copy. The gap in AF probability and time in AF between CMP and cCMP is attributed to complex critical structures whose activation requires fewer conduction-block failures than their deactivation, including self-contained asymmetric structures and coupled structures that re-initiate on termination. The paper concludes that persistent AF arises from these asymmetric re-entrant circuits, so that parameter variation or remodelling is not needed to explain diverse AF persistence.","tokens_in":30924,"tokens_out":4847,"duration_ms":49519,"significance":"If the central claim holds, the paper provides a mechanism for the paroxysmal-to-persistent transition in a minimal model, with potential implications for how local fibrosis architecture, rather than global burden, determines AF persistence. The work has genuine strengths: the continuous-time mean-field solution in Appendix A is an exact derivation; the enhanced MF model in Appendix B is a useful robustness check; the simulations extend to 10^9 timesteps, much longer than typical biophysical models; and the structural examples in Figs. 9 and 10 are clearly presented. The paper also explicitly acknowledges several limitations of the CMP model and the difficulty of verifying complete detection of critical structures. However, the central evidence for the mechanism rests on a controlled comparison that is confounded by changes in total block-node density, and the inference to higher-order structures is not quantitatively validated. These issues are load-bearing for the main claim, so the paper requires major revision rather than acceptance in its current form.","major_comments":[{"comment":"The cCMP control does not isolate complex critical structures because it changes the total density of conduction-blocking nodes. In cCMP, block-susceptible nodes are placed only on eligible nodes within isolated segments of length at least tau/2, at fraction delta of those eligible nodes, while the CMP comparison places block nodes at probability delta on all L^2 nodes. The two models therefore differ not only in the presence of multi-block complex structures but also in total block-node count, refractory landscape, and wavefront-collision effects. The conclusion in Section IV B that the CMP-cCMP gap is caused by complex asymmetric structures requires a control in which total block-node density and, ideally, the number of simple critical structures are matched while only the complex multi-block configurations are suppressed.","section":"Section IV A"},{"comment":"The statement that 'by construction, the cCMP and CMP models contain the same number of simple critical structures' is not supported by the construction described in Section IV A. The cCMP model places block nodes only on selected segment nodes, whereas the CMP copy places them randomly across the whole lattice; these two procedures do not generally yield the same number of simple critical structures. If the counts differ, the subsequent inference that 'some higher order critical structures must exist' is a non sequitur, because the observed gap could be due to a different number of simple structures or to different total density rather than to complex structures.","section":"Section IV B"},{"comment":"The collapse of the CMP and cCMP curves at low delta in Fig. 12 is presented as evidence that complex structures dominate at high delta, but this collapse is also exactly what would be expected if the gap were driven by the total density of conduction-blocking nodes, since the absolute density difference between the two models shrinks as delta decreases. The figure therefore does not discriminate between the proposed asymmetric-structure mechanism and a generic density effect. A matched-density control, or a direct quantitative count of complex structures, is needed to make the mechanism claim.","section":"Section IV B and Fig. 12"},{"comment":"The paper itself states in Section IV A that there is 'no easy method to verify that all circuits have been detected,' and the central mechanism is supported by hand-selected examples (Figs. 9-11) and by elimination through the cCMP comparison. Given the confound identified above, the existence and causal role of asymmetric complex structures is not quantitatively established. The authors should either provide a validated detection algorithm for complex structures or design a cleaner intervention that affects only multi-block complex configurations.","section":"Section IV A and Section VI"}],"minor_comments":[{"comment":"There is a typo in the first sentence: 'arrhytmia' should be 'arrhythmia', and later 'assymetry' should be 'asymmetry'.","section":"Abstract"},{"comment":"The phase diagrams show averages over 200 simulations but no error bars or confidence intervals, making it difficult to assess whether the reported differences, especially near the sharp transitions, are statistically meaningful. Adding confidence bands would strengthen the quantitative claims.","section":"Figs. 7, 8, 12"},{"comment":"The AF threshold is defined as 1.1 times L, which equals 220 for L=200, but the text states '1.1×L (220) nodes' without noting that this is specific to the chosen lattice size. Please clarify whether the threshold is 220 nodes or 1.1L generally.","section":"Section II C, Eq. (5)"},{"comment":"The statement that the MF model 'spends significantly less time in AF' than the CMP model is not accompanied by any statistical test or uncertainty estimate; a quantitative comparison of distributions or confidence intervals would be appropriate.","section":"Section III"},{"comment":"The eMF model is described as 'perfectly compatible' with the MF model, but again no error bars or statistical measures are given; a quantitative statement of compatibility would be more informative.","section":"Appendix B"},{"comment":"The error bars in Fig. 19 are stated to be 95% confidence intervals over 50 simulations, but the main-text figures (e.g., Figs. 7 and 8) do not report the number of simulations beyond '200 simulations'; please make the simulation counts and uncertainty measures consistent throughout.","section":"Appendix E, Fig. 19"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a clinically relevant question and contains a genuinely useful analytic mean-field solution. However, the central controlled comparison (cCMP vs. CMP) is confounded by total block-node density, and the claim of equal simple-structure counts by construction appears unsupported. If the authors can provide a density-matched control or a direct validated detection of complex structures, the paper could become acceptable. The lack of error bars on the main phase diagrams is also a concern given the emphasis on quantitative differences."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWhat you should know: this is a genuinely useful paper. It takes the Christensen-Manani-Peters model and, instead of just running it, pulls it apart layer by layer. The mean-field model with analytic steady-state solution, the enhanced MF, and the controlled CMP are each introduced to isolate what makes persistent AF possible in the model. The central claim is that asymmetric critical structures - circuits that require fewer failures to activate than to terminate - explain why the same parameters can yield episodes from seconds to months. That claim is actually supported by the structure of the paper: the MF and cCMP models both get the phase boundary roughly right but miss the long tails, and only the full CMP has the long-lived asymmetric structures. The authors are honest that results are qualitative and that event durations need not map clinically. Good credit where due: the cMF formula in Eq. A7 is a real derivation, and the eMF model rules out the obvious alternative explanation that non-spatial details of simple circuits account for the gap.\n\nThe soft spots are real but not fatal. The cCMP control is the main one. The stress-test note is right: cCMP removes both the complex structures and most of the random block nodes, so the CMP-cCMP gap conflates total block density with structure complexity. The authors notice the density effect and show that at low delta the curves collapse, but that does not establish that the remaining gap at high delta is caused by the proposed asymmetric circuits rather than by generic density effects on wavefront collisions and refractoriness. They also identify the complex structures by inspection - hand-picked examples from figs 9 and 10 - and admit there is no systematic detection algorithm. The phase diagrams have no error bars or statistical tests, so the reader cannot judge whether the CMP-cCMP differences at intermediate delta are robust. None of this refutes the core idea, but it means the paper's central comparison is suggestive rather than conclusive.\n\nWho this is for: anyone modeling AF with discrete cellular automata, or interested in how simple rules produce heterogeneous arrhythmia persistence. It deserves a serious referee and, with the control fixed (matching total block density, adding a detection algorithm, or at least variance estimates), it would be a solid contribution. I would cite it for the MF derivation and the asymmetric-structure hypothesis.\n\nRecommendation: send it to review, with the clear expectation that the cCMP comparison needs a density-matched control before the mechanistic claim is accepted.","headline":"A careful dissection of the CMP model that identifies asymmetric re-entrant circuits as the driver of persistent AF, with a real but fixable weakness in the cCMP control.","tokens_in":31397,"tokens_out":608,"would_cite":true,"duration_ms":9511,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that in a simple model of atrial fibrillation, persistent episodes arise from re-entrant circuits that are easier to activate than to deactivate, with no change in model parameters.","keywords":["atrial fibrillation","paroxysmal to persistent transition","re-entrant circuits","micro-anatomical re-entry","conduction block","cellular automaton model","mean-field birth-death model","critical structures"],"falsifier":"Run the cCMP construction but instead of deleting the 'extra' conduction-blocking nodes, keep them in place with failure probability set to zero (inert nodes), so the lattice geometry and connectivity are identical to the CMP model while no conduction block can occur outside simple segments. If time in AF in this inert-node control matches the cCMP value, the gap is due to the presence of blocking-prone nodes in complex configurations; if it rises toward the CMP value, the paper's attribution to complex circuits is not the whole story.","tokens_in":30462,"feed_emoji":"🫀","tokens_out":6429,"duration_ms":65460,"temperature":0.7,"pith_summary":"This paper is a dissection of the Christensen-Manani-Peters (CMP) model of atrial fibrillation, a cellular automaton on a square lattice in which randomly placed conduction-blocking nodes occasionally fail to fire. The authors want to know why the model, at fixed parameter values, can show every clinical AF pattern: permanent sinus rhythm, short self-terminating paroxysmal episodes, occasional long episodes, and persistent AF lasting a simulated month. They establish that the difference is made by 'complex' re-entrant circuits whose activation probability exceeds their deactivation probability, because they require fewer consecutive conduction-block failures to start than to stop. A mean-field birth-death version and a controlled spatial version that contain only simple circuits both underestimate time in AF, and only the full model with complex circuits reproduces the observed spectrum. The clinical significance, if the claim transfers, is that a patient's paroxysmal versus persistent status could depend on the local arrangement of fibrosis-related conduction blocks rather than on total fibrosis burden or progressive remodelling.","feed_headline":"AF can last seconds or months from one local circuit asymmetry","feed_subtitle":"A lattice model shows re-entrant circuits that are easier to start than to stop explain both paroxysmal and persistent AF.","key_machinery":"The load-bearing object is the CMP lattice model's 'critical structure', a local arrangement of a conduction-blocking node and a transverse-connection-free segment of length at least $\\tau/2$ that can host a re-entrant circuit when the block fails. The paper classifies structures as simple or complex depending on how many successive conduction-block failures are needed to activate and to terminate them. Two auxiliary constructions carry the argument: a mean-field model in which each simple critical structure is a particle undergoing birth-death dynamics with rates set by $\\epsilon/T$ and $\\epsilon/\\langle \\ell \\rangle$, and a controlled CMP model (cCMP) in which conduction-blocking nodes are placed only on isolated simple-critical segments so that complex circuits cannot form. The difference in AF probability and time in AF between the full model and these simplifications is what isolates the contribution of asymmetric complex circuits. The analytic risk formula $R=1-[1-(1-\\nu_\\perp)^\\tau]^{\\delta L^2}$ provides the baseline against which the model's excess AF is measured.","core_discovery":"The paper's central claim is that persistent AF in the CMP model is caused by re-entrant circuits with an asymmetry in the probability of activation relative to deactivation. A simple circuit turns on when one susceptible node fails and turns off when one susceptible node fails, so its on and off rates are balanced, and a mean-field birth-death description captures its behaviour; the continuous mean-field solution shows that a finite set of such circuits cannot keep the system in AF almost all the time. Complex circuits break this balance: some need two successive cell failures to initiate but four to terminate, and others are coupled so that the termination of one circuit immediately re-ignites a neighbour. With these circuits present, the same lattice parameters can generate fibrillatory episodes from a few seconds to more than a month, and the gap between the full model and the controlled model collapses as the density of conduction-blocking nodes is reduced and complex circuits become rare.","pith_inferences":["A testable extension is to build an 'inert-node control': in a cCMP lattice, keep the deleted conduction-blocking nodes in place but set their failure probability to zero; if the cCMP-CMP gap persists, connectivity changes, not circuit asymmetry, are responsible.","If the asymmetry mechanism generalizes, the distribution of AF episode durations in a fixed patient substrate should be heavy-tailed over many orders of magnitude, much wider than a single exponential; clinical event-duration histograms could be inspected for such spread before attributing persistent AF to remodelling.","The number of required simultaneous failures suggests a quantitative link between local fibrosis density and persistence: scanning a realistic fibre map for regions where multiple vulnerable nodes lie within one activation cycle of a re-entry path should predict where persistent drivers are anchored."],"forward_implications":["If the claim is correct, the full spectrum of AF persistence in the model emerges without changing any parameter; a single fixed coupling value $\\nu_\\perp=0.11$ can yield sinus rhythm, paroxysmal AF, and persistent AF depending only on the local positions of vulnerable nodes.","Mean-field birth-death descriptions with a finite number of independent drivers cannot explain persistent AF, so any model that treats AF persistence as a simple on/off balance of fixed-rate drivers will miss the mechanism.","Reducing the fraction of conduction-blocking nodes $\\delta$ suppresses complex circuits and makes the full model's behaviour coincide with the controlled model, identifying high local density of vulnerable nodes as proarrhythmic in a way that simple-circuit count alone does not capture.","Ablation strategies would preferentially need to target complex or coupled structures whose termination requires many simultaneous failures, rather than merely any region that can host a simple re-entrant circuit."],"supporting_citations":[{"why":"Introduces the CMP percolation model and the theoretical risk formula that this paper dissects and extends.","marker":"[29]"},{"why":"Documents that identical parameters produce heterogeneous AF persistence in the CMP model, the phenomenon the paper attributes to complex circuits.","marker":"[38]"},{"why":"Provides the parameter-space results on how delta, tau, and epsilon affect time in AF, used to interpret the cCMP versus CMP phase diagrams.","marker":"[43]"},{"why":"Serves as the historical cellular-automaton model of AF against which the current micro-anatomical re-entry account is distinguished.","marker":"[54]"},{"why":"A comparable discrete fibrosis model whose percolation-threshold result connects fibrosis clusters to re-entry, used as a neighbouring result to situate the CMP mechanism.","marker":"[56]"}],"fun_headline_variants":["Re-entrant circuit asymmetry flips paroxysmal to persistent AF","One asymmetry causes AF episodes from seconds to months","Asymmetric re-entry circuits decide AF duration in model","From seconds to months: AF duration set by circuit asymmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument that the cCMP-CMP gap is caused by complex asymmetric circuits assumes that removing all conduction-blocking nodes that are not part of an isolated simple critical segment changes nothing else about wavefront collisions, refractoriness, or suppression of neighbouring circuits; if deleting those nodes alters wave dynamics for reasons unrelated to circuit complexity, the measured gap could be produced without the proposed mechanism.","fun_headline_variants_meta":{"raw":{"variants":["Re-entrant circuit asymmetry flips paroxysmal to persistent AF","One asymmetry causes AF episodes from seconds to months","Asymmetric re-entry circuits decide AF duration in model","From seconds to months: AF duration set by circuit asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000241,"raw_usage":{"total_tokens":1533,"prompt_tokens":970,"completion_tokens":563,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":586,"completion_tokens_details":{"reasoning_tokens":497}},"tokens_in":586,"tokens_out":563,"duration_ms":6487,"temperature":1.0,"reasoning_tokens":497,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:06:16.674249+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the cCMP construction but instead of deleting the 'extra' conduction-blocking nodes, keep them in place with failure probability set to zero (inert nodes), so the lattice geometry and connectivity are identical to the CMP model while no conduction block can occur outside simple segments. If time in AF in this inert-node control matches the cCMP value, the gap is due to the presence of blocking-prone nodes in complex configurations; if it rises toward the CMP value, the paper's attribution to complex circuits is not the whole story.","supporting_citations":[{"cited_title":"Christensen, K","cited_arxiv_id":null,"evidence_quote":"Introduces the CMP percolation model and the theoretical risk formula that this paper dissects and extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents that identical parameters produce heterogeneous AF persistence in the CMP model, the phenomenon the paper attributes to complex circuits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the parameter-space results on how delta, tau, and epsilon affect time in AF, used to interpret the cCMP versus CMP phase diagrams."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Serves as the historical cellular-automaton model of AF against which the current micro-anatomical re-entry account is distinguished."},{"cited_title":"AF begets AF","cited_arxiv_id":null,"evidence_quote":"A comparable discrete fibrosis model whose percolation-threshold result connects fibrosis clusters to re-entry, used as a neighbouring result to situate the CMP mechanism."}],"review_version":1}