{"id":"bd35d327-3b35-4371-aac4-a7f8474f3ff0","arxiv_id":"2502.08480","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In a simulated 1D cardiac fiber, the magnetic field computed from the action potential reproduces spatially discordant alternans node positions, supporting magnetic field monitoring as a potential non-invasive arrhythmia biomarker.","lead":"The authors simulate a 1D cardiac fiber and compute the tiny magnetic field its electrical activity produces, showing that this field mirrors beat-to-beat changes in action potential duration, including spatially discordant alternans. A smart generalist might read this because it suggests magnetocardiography could become a non-invasive way to spot arrhythmia precursors that today are seen only with electrical recordings.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"SDA detection claim relies on a 50 um standoff and neglect of extracellular return currents; at realistic sensor distances the B-field node may not survive, so 'reliably detect SDA' is not yet established.","rationale":"The reader's weakest_assumption already identifies the load-bearing issue: the field is computed from the same voltage that defines SDA, the evaluation distance is 50 um, volume currents are neglected, and no noise or sensor standoff is included. My stress-test focus is narrower and slightly more pointed: Eq. (19) evaluated at 50 um is a spatially local (~100 um-wide) filtering of dVm/dz, so the node correspondence in Fig. 3 is essentially a consequence of the local relation between Vm and B in the tuned model, rather than an empirical or independent observable. This does not make the paper internally inconsistent; the forward calculation is a valid mathematical exercise. But the phrase 'reliably detect SDA' implies a measurable signal under realistic conditions, and none of the three obstacles to that claim (volume-current cancellation, standoff, sensor noise) is addressed. I therefore agree with the reader's conditional verdict: the computational demonstration is sound and interesting, but the clinical/diagnostic conclusion should be conditional on a more realistic sensor/volume-conductor model. No verdict change is needed relative to the reader.","tokens_in":11625,"tokens_out":6510,"duration_ms":77838,"concrete_test":"Extend the forward model to include the extracellular volume conductor: simulate the same 1D fiber surrounded by saline with interstitial conductivity sigma_e, solve the bidomain or finite-element volume-conductor problem, and compute B from the total current (intracellular axial current plus extracellular return currents) at radial standoffs of 50 um, 250 um, 500 um, and 1 mm, adding white sensor noise at state-of-the-art SQUID or NV magnetometer levels. Compare Z_B to Z_APD and the beat-to-beat |B|max alternans amplitude at each standoff. If Z_B deviates from Z_APD by more than a few millimeters, or if the alternans difference falls below the noise floor, at any standoff >250 um, then the paper's central claim does not survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that magnetic field measurements reliably detect spatially discordant alternans because Z_APD and Z_B agree within millimeters in Fig. 3. However, B is computed from the same simulated Vm via Eqs. (18)-(19) and evaluated only 50 um from the fiber, where the Biot-Savart kernel has width ~50 um and effectively acts as a local spatial derivative of Vm. The agreement between APD and |B|max nodes is therefore a property of this near-field, no-volume-conductor limit, not an independent validation. In a real measurement, the sensor is outside the tissue/saline bath (standoff of at least hundreds of microns to millimeters), extracellular volume currents produce a partially cancelling return field, and noise limits the smallest detectable alternans difference. The paper excludes all three by construction. The observation that at 29 C the node positions already differ by 1.0 mm (1.11 vs 1.21 cm) while at 37 C they differ by 0.3 mm (1.08 vs 1.11 cm) underscores that the correspondence is not exact and could degrade further at realistic standoff. Without an end-to-end test that includes a volume conductor and a sensor model, the claim that magnetic measurements reliably detect SDA is not supported by the present evidence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper simulates a one-dimensional cardiac fiber with a four-variable phenomenological model under a pacing-down restitution protocol at temperatures between 29 °C and 40 °C, computes the magnetic field from the simulated transmembrane potential using J = −σ∂Vm/∂z and Biot–Savart integration (Eqs. (18)–(19)), and compares APD80-based and |B|max-based markers of spatially discordant alternans (SDA). The authors report electric and magnetic restitution curves, spatial maps of APD80 and |B|max, and node positions for SDA at selected pacing cycle lengths, and they conclude that magnetic field measurements reliably detect SDA and temperature-dependent changes in cardiac action potentials.","tokens_in":11795,"tokens_out":6599,"duration_ms":75955,"significance":"As an in silico proof of concept, the paper offers a clear and internally consistent forward-modeling pipeline: the electrical model, current-density estimate, and magnetic-field integral are explicit and reproducible in structure, and the observation that a simple scalar magnetic marker inherits the spatial organization of SDA is a useful methodological pointer for biomagnetic monitoring. The connection between |B|max and the upstroke dynamics via Eqs. (18)–(20) is clearly laid out. However, the significance as a measurement claim is limited: the magnetic field is computed from the same simulated voltage that defines the APD-based ground truth, so the node agreement in Fig. 3 is an in-model consistency result rather than an independent validation. The absence of a volume conductor, sensor standoff, and noise means the paper supports a qualitative prediction, not the stated conclusion that magnetic field measurements 'reliably detect' SDA in a practical setting.","major_comments":[{"comment":"The abstract and the final paragraph of the SDA results state that magnetic field measurements 'effectively detect' and 'reliably detect' SDA. This is the central claim, but the only supporting evidence is the spatial agreement between APD80 and |B|max nodes in Fig. 3, where both quantities are computed from the same simulated Vm through Eqs. (18)–(19). The agreement is therefore an in-model consistency result, not an independent measurement. To support the detection wording, the paper should either add an end-to-end model with a realistic volume conductor and a sensor model (standoff, noise, finite sensitivity) or systematically rephrase all detection conclusions as statements about the computed magnetic field in the idealized simulation. The second option is the minimum needed for the current scope.","section":"Abstract and Results, 'Electric and magnetic SDA'"},{"comment":"The neglect of extracellular volume currents is justified by the statement that the action potential wavelength is much larger than the 50 µm evaluation distance. That condition alone is not sufficient: magnetic-field cancellation by extracellular return currents depends on the conductivity and geometry of the bath and on the sensor distance, as in the cited Swinney–Wikswo calculation [47]. Additionally, a 50 µm distance from the source axis is only 10 µm from the surface of the 40 µm radius fiber, which is far below any practical sensor standoff in tissue. Please provide a quantitative estimate of the extracellular contribution at the evaluation distance, or include a bidomain/saline-bath simulation, and state over what standoff range the computed node positions remain stable. Without this, the claim that the computed field approximates a measurable signal is not fully established.","section":"Methods, 'Magnetic Field Modeling'"},{"comment":"The model parameters are described as 'opportunely fine-tuned to reproduce SDA' in the one-dimensional cable. Because the SDA regime is partly constructed by parameter choice, the subsequent conclusion that the magnetic field 'demonstrates' SDA detection overstates the generality of the result. The paper should explicitly frame the SDA conditions as a model scenario and, ideally, include a sensitivity analysis over parameters and pacing cycle lengths to show that the node correspondence between APD and |B|max is robust rather than a consequence of the chosen tuning.","section":"Methods, 'Cardiac Electrophysiology Modeling'"}],"minor_comments":[{"comment":"Table I reports initial conditions as 'u =, v = 1, w = 0, s = 1', but the numerical value for u is missing; please supply it.","section":"Methods, 'Cardiac Electrophysiology Modeling'"},{"comment":"Equation (19) uses I(t, z′) but the current I is not defined in terms of the current density J from Eq. (18); please specify I = J · A (or the equivalent cross-sectional integral) and state the units.","section":"Methods, 'Magnetic Field Modeling'"},{"comment":"The pacing-down restitution protocol is described only qualitatively. Please report the PCL range, decrement step, number of beats at each PCL, and the criteria used to define 'alternans onset' and 'fully developed alternans' so that the results are reproducible.","section":"Methods, 'Cardiac Electrophysiology Modeling' and Results"},{"comment":"Fig. 3 reports node positions for a single PCL per temperature and regime. Please state whether the node positions are stable across neighboring PCLs or provide a measure of variability, since the claim of 'a few millimeters' agreement is based on point values.","section":"Results, 'Electric and magnetic SDA'"},{"comment":"The statement that the magnetic restitution curves show the alternans onset 'more easily identifiable' than APD curves is not quantified. Please define a quantitative detection criterion (e.g., beat-to-beat difference threshold) or soften the claim.","section":"Results, 'Electric and magnetic alternans'"},{"comment":"There are minor typographical issues, including 'elctrical' in the Methods heading and 'Moore (η(T))' where the intended name may be something else; please proofread.","section":"Methods and Discussion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of physics.med-ph and builds naturally on the authors' prior work [44]. The main barrier is the gap between the simulation-based evidence and the abstract's 'reliably detect' wording; this is fixable by a careful revision that either adds modeling of the measurement chain or consistently qualifies the detection claims. I do not see a need for new experiments at this stage, but the authors should also consider citing experimental magnetocardiography work on alternans to anchor the practical relevance of the proposed biomarker."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a competent simulation study with a genuinely new observation—SDA node locations are visible in the |B|max spatial map and track APD nodes within a millimeter or so. The magnetic restitution curves following upstroke velocity is a nice consistency check, not a surprise given Eqs. (18)–(20). The paper knows what it is: a 1D, homogeneous cable, parameters fine-tuned to produce SDA, field evaluated 50 μm from the source with volume currents neglected.\n\nWhat it does well: methods are specified in enough detail to reproduce (no code, but equations and parameters are complete), the forward calculation from Vm to B is internally consistent, and the node-position agreement is a real result in that model. The discussion honestly names the main limiting factors: geometry and homogeneity.\n\nSoft spots, in order: the detection claim outruns the evidence. B is a deterministic function of the same Vm that defines SDA, so correspondence is built in; and at 50 μm standoff the Biot-Savart kernel is basically a local derivative operator. A real sensor sits further away and sees a volume conductor with return currents that partially cancel the field—the paper explicitly relies on [47] to wave this away, but the wavelength argument doesn't save the node map if the standoff is a millimeter or more. The 29 °C case already shows 1.0 mm discrepancy, so the exactness is fragile. The 'opportunely fine-tuned' parameters are another sign this is a proof-of-concept, not a predictive claim.\n\nIf the abstract and title were softened to 'a possible non-invasive indicator' and a feasibility framing, I'd be comfortable. As it stands, 'reliably detect SDA' is not supported beyond the model's own near-field limit. This is a moderate issue, not a fatal one; the core idea is worth a serious referee.\n\nMy verdict: send to review, but with a request for a volume-conductor/sensor-standoff test or a careful rewrite of the claims. The paper is for someone working on biomagnetic forward problems, not for clinical cardiology yet.","headline":"A clean computational feasibility study showing SDA nodes are reproduced in the near-field magnetic norm, but the 'reliably detect' claim needs a volume-conductor and sensor model before it is earned.","tokens_in":12450,"tokens_out":1461,"would_cite":false,"duration_ms":16058,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using a 1D cardiac fiber model, this paper shows that the magnetic field produced by heart tissue exposes spatially discordant alternans, with magnetic and electrical node locations matching within a few millimeters.","keywords":["spatially discordant alternans","cardiac magnetic field","biomagnetism","cardiac alternans","restitution curve","action potential upstroke","temperature dependence","Biot-Savart law"],"falsifier":"Compute the same SDA protocol with a bidomain model that includes extracellular volume currents and evaluate the magnetic field at a realistic sensor standoff of 1 mm to 1 cm; if the $|B|_{\\max}$ node shifts by more than the few-millimeter agreement reported here, or if the field amplitude falls below sensor noise, the central claim fails. An experimental check would be a paced cardiac fiber or monolayer measured with a high-sensitivity magnetometer at about 50 µm standoff, testing whether the magnetic node coincides with the APD node within a few millimeters.","tokens_in":11312,"feed_emoji":"🧲","tokens_out":7228,"duration_ms":64690,"temperature":0.7,"pith_summary":"This paper uses a one-dimensional computer model of a cardiac fiber to argue that the magnetic field generated by heart tissue can expose spatially discordant alternans (SDA), a rhythm disturbance in which neighboring regions alternate out of phase, creating steep repolarization gradients that can trigger reentry. The authors show that the spatial map of the peak magnetic field magnitude $|B|_{\\max}$ reproduces the out-of-phase pattern seen in action potential duration, with the node position—the location where alternans changes phase—matching between electrical and magnetic indicators to within a few millimeters ($Z_{\\mathrm{APD}}=1.08$ cm versus $Z_B=1.11$ cm at 37 °C). They also find that magnetic restitution curves invert relative to classic action-potential-duration curves but closely track the maximum upstroke velocity $\\partial V/\\partial t_{\\max}$, because the magnetic field is tied to the action potential's spatial derivative. If these results carry to real tissue, magnetic field measurements could become a non-invasive way to locate SDA nodes and assess arrhythmic risk without inserting electrodes.","feed_headline":"Magnetic fields expose hidden nodes of risky heart rhythm","feed_subtitle":"If it holds, biomagnetic sensors could spot arrhythmia precursors without touching the heart.","key_machinery":"The machinery that carries the argument is the derivative relation between the action potential and the magnetic field. For a straight cylindrical fiber the current density is $J(t,z')=-\\sigma\\,\\partial V_m/\\partial z'$, and Biot-Savart integration gives a field whose magnitude peaks at the depolarization upstroke, where $\\partial V_m/\\partial t$ is largest. Because alternans changes the upstroke timing and amplitude in adjacent regions, the peak field magnitude inherits the alternating pattern, and the position where the two beats have equal $|B|_{\\max}$ defines the magnetic node. The simulations use a temperature-dependent four-variable phenomenological model with Moore and Arrhenius factors, and the analysis uses pointwise restitution curves rather than spatial averages to avoid washing out SDA.","core_discovery":"The central claim is that spatially discordant alternans have a magnetic signature that a simple scalar indicator can capture. On a 3 cm fiber, the paper demonstrates that for two consecutive beats (n−1 and n), the spatial profile of the peak magnetic field norm $|B|_{\\max}$ alternates out of phase in the same way as $\\mathrm{APD}_{80}$, and the nodes of the two indicators coincide within a few millimeters at 37 °C and 33 °C. The paper further claims that the magnetic restitution curves—$|B|_{\\max}$ versus pacing cycle length—are ordered inversely to APD restitution curves but mirror the upstroke velocity curves, because the current density $J=-\\sigma\\,\\partial V_m/\\partial z'$ and hence the Biot-Savart field are proportional to the spatial derivative of the transmembrane potential. This establishes, in the authors' words, a 'strong spatial correlation between electrical and magnetic indicators' and supports the capability of magnetic field measurements to reliably detect SDA.","pith_inferences":["In a real tissue bath, extracellular volume currents are likely to partially cancel the intracellular magnetic field; at standoff distances of millimeters to centimeters, the node correspondence seen at 50 µm may blur or shift, a risk the paper does not quantify.","The magnetic field's lock to the upstroke suggests a two-channel diagnostic: combining magnetic maps (sodium-driven upstroke) with optical or electrical APD maps (repolarization) could separate channelopathy subtypes, a testable hypothesis for Brugada or Long-QT 3.","If single-cell or fiber-level NV-diamond magnetometers reach the nT sensitivity this model predicts at 50 µm, the predicted node coincidence could be checked directly in a paced cardiac fiber preparation.","The model parameters were fine-tuned to reproduce SDA; a natural next step is to test whether the magnetic signature persists with parameter sets derived from experimental restitution data rather than tuned to the phenomenon."],"forward_implications":["Magnetic field mapping can identify the location of SDA nodes without electrode contact, potentially guiding ablation or pacing strategies aimed at arrhythmia prevention.","Because magnetic restitution curves track the upstroke rather than APD, magnetic measurements are a more direct readout of fast sodium-channel dynamics than standard electrical indicators.","The temperature dependence of node position in the model implies that magnetic field maps could track how hypothermia shifts arrhythmic risk in clinical cooling protocols.","The magnetic bifurcation point in the restitution protocol is easier to detect than the APD bifurcation point, so magnetic biomarkers may detect alternans onset earlier.","The same 1D thermo-electro-magnetic framework can be extended to more complex geometries to test whether the node correspondence persists in tissue and whole hearts."],"supporting_citations":[{"why":"Provides the thermo-electro-magnetic modeling framework that this paper extends to alternans and SDA.","marker":"[44]"},{"why":"Supplies the temperature-dependent Moore and Arrhenius factors used to modulate currents and gating variables.","marker":"[43]"},{"why":"Source of the four-variable phenomenological cardiac model and alternans dynamics.","marker":"[4]"},{"why":"Supplies the minimal human ventricular action potential model and initial conditions used in the simulations.","marker":"[45]"},{"why":"Justifies neglecting extracellular volume currents at the 50 µm evaluation distance because the action potential wavelength is much larger.","marker":"[47]"},{"why":"Supports the feasibility of measuring cellular-level magnetic signals with diamond quantum sensors, motivating the detection claim.","marker":"[39]"},{"why":"Provides the experimental temperature-dependent alternans onset data that the model is compared with.","marker":"[11]"}],"fun_headline_variants":["Magnetic fields unmask hidden nodes of arrhythmic risk","Heart magnetism exposes subtle rhythm danger spots","Magnetic field scans reveal risky beat alternation nodes","Biomagnetic sensing spots cardiac alternans nodes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that magnetic field measurements can detect SDA rests on the assumption that the field computed from an isolated 1D fiber in vacuum—neglecting extracellular volume currents and sensor noise, and evaluated at 50 µm—faithfully represents a measurable signal; if volume-current cancellation or practical sensor distance distorts the field, the node correspondence would not survive.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic fields unmask hidden nodes of arrhythmic risk","Heart magnetism exposes subtle rhythm danger spots","Magnetic field scans reveal risky beat alternation nodes","Biomagnetic sensing spots cardiac alternans nodes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1319,"prompt_tokens":895,"completion_tokens":424,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":363}},"tokens_in":511,"tokens_out":424,"duration_ms":5256,"temperature":1.0,"reasoning_tokens":363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T04:57:04.429496+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same SDA protocol with a bidomain model that includes extracellular volume currents and evaluate the magnetic field at a realistic sensor standoff of 1 mm to 1 cm; if the $|B|_{\\max}$ node shifts by more than the few-millimeter agreement reported here, or if the field amplitude falls below sensor noise, the central claim fails. An experimental check would be a paced cardiac fiber or monolayer measured with a high-sensitivity magnetometer at about 50 µm standoff, testing whether the magnetic node coincides with the APD node within a few millimeters.","supporting_citations":[{"cited_title":"Crispino, M","cited_arxiv_id":null,"evidence_quote":"Provides the thermo-electro-magnetic modeling framework that this paper extends to alternans and SDA."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the temperature-dependent Moore and Arrhenius factors used to modulate currents and gating variables."},{"cited_title":"Gizzi, E","cited_arxiv_id":null,"evidence_quote":"Source of the four-variable phenomenological cardiac model and alternans dynamics."},{"cited_title":"Bueno-Orovio, E","cited_arxiv_id":null,"evidence_quote":"Supplies the minimal human ventricular action potential model and initial conditions used in the simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies neglecting extracellular volume currents at the 50 µm evaluation distance because the action potential wavelength is much larger."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the feasibility of measuring cellular-level magnetic signals with diamond quantum sensors, motivating the detection claim."},{"cited_title":"Crispino, A","cited_arxiv_id":null,"evidence_quote":"Provides the experimental temperature-dependent alternans onset data that the model is compared with."}],"review_version":1}