{"id":"4eb323f0-7d01-4b43-8106-64fc3c19f2f9","arxiv_id":"2608.05414","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Chaotic field-line transport in the conservative Tokamap follows fitted power-law scalings, including a recurrence-time law τ_c ∝ x_q^{-0.213}.","lead":"This paper uses a simplified computer model of a fusion reactor to study how magnetic field lines wander and spread. It finds a mathematical pattern that links a magnetic-shear parameter to the time scale of that spreading.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scaling claim τ_c ∝ x_q^{-0.213} rests on an exponential fit to the initial segment S(τ)>0.08 of single-orbit survival probabilities; stickiness lives in the discarded tail, so the exponent may be a fit-window artifact rather than a transport law.","rationale":"The paper's headline claim is the algebraic scaling of the characteristic transport time. The only support for this claim is the decay coefficient B extracted from S(τ)>0.08 on single orbits. If the survival distribution is not exponential, B is not the characteristic transport time and the exponent 0.213 does not quantify long-time transport. This is not an outside-consensus disagreement; it is an internal mismatch between the observable fitted (short-time exponential decay) and the mechanism claimed (long-time stickiness). The reader's weakest assumption names the same issue, and I agree. A threshold-sensitivity and tail-exponent test would settle it. Other defects—lack of code/data, radial-coordinate inconsistency, undefined δ in a caption—are real but secondary; they would support a conditional or unverdictable status but are not the central logical flaw. Since the reader already returned CONDITIONAL with this concern identified, my independent read does not change the verdict; it sharpens the required test.","tokens_in":13208,"tokens_out":8484,"duration_ms":75814,"concrete_test":"Compute for each x_q ∈ {10^-4,...,10^-8} 100 independent 10^9-iteration orbits, record all recurrences to the same box, and fit |B| in three windows S>0.02, S>0.08, S>0.2, with bootstrap errors; separately fit the tail S<0.08 to S∼τ^{-γ}. If the inferred exponent in Eq. (44) shifts by more than its bootstrap error across windows or disagrees with the tail scaling, then τ_c∝x_q^{-0.213} is a fitting artifact, not a robust transport law.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result, Eq. (46), rewrites the fitted decay coefficient of Eq. (43) as an inverse transport time. But Eq. (43) is fit only over the region S(τ)>0.08 (Section 2.2, Fig. 4). In a mixed Hamiltonian phase space the survival probability is generically expected to develop a power-law tail, S(τ)∼τ^{-γ}, controlled by sticky dynamics near islands and cantori. Fitting only the initial exponential portion measures the transient escape rate from the chaotic component, not the asymptotic trapping that the paper invokes to explain the slowdown. Thus the statement that 'long-time dynamics are increasingly dominated by stickiness' is not supported by the observable actually fitted; the discarded tail is precisely where stickiness would appear. The 'characteristic transport time' τ_c=1/|B| is only meaningful for an exponential law; for a power-law-tailed distribution the mean recurrence time can be infinite or dominated by rare events, so an exponential rate is not a robust transport characteristic. Moreover, each survival curve comes from a single orbit (Section 2.2), no error bars are reported on B, and the five decade points in Fig. 4(b) carry no uncertainties; the fitted exponent 0.213 is therefore not protected against threshold choice, orbit length, or initial-condition dependence. Secondary inconsistencies—Ψ stated to lie in [0,1] but Figures 1 and 2 and the recurrence box (T,Ψ)=(0.5,5.0) using Ψ≈5–40, and a caption introducing an undefined δ=10^{-6}—further weaken confidence that the reported parameter regime is well controlled.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the conservative Tokamap, an exactly symplectic map for magnetic field-line dynamics, and combines phase-space portraits, Lyapunov-exponent calculations, ensemble-averaged radial transport, and Poincaré recurrence statistics. The central claims are: (i) ensemble-averaged transport exhibits dynamic scaling with growth, saturation, and crossover exponents satisfying z = δ/β and producing a universal data collapse; and (ii) the characteristic recurrence time grows algebraically as τ_c ∝ x_q^{-0.213}, which the authors attribute to increasing stickiness near KAM islands and cantori as the magnetic-shear parameter x_q is reduced. The dynamic-scaling part is internally consistent and the collapse is visually convincing. The recurrence part, however, is the weakest link: it rests on a single orbit per parameter value, an exponential fit restricted to S(τ) > 0.08, and no reported uncertainties, so the headline exponent 0.213 is not yet robustly established.","tokens_in":13550,"tokens_out":4013,"duration_ms":40218,"significance":"If the results are confirmed, the paper would provide a useful quantitative bridge between phase-space geometry and long-time transport in a widely studied plasma-physics map. The dynamic-scaling framework is a genuine strength: the exponents β, δ, and z are measured independently from the ensemble transport, the scaling relation z = δ/β is satisfied within the reported errors, and the collapse in Fig. 2(b) is a nontrivial consistency test. The recurrence analysis aims at an important and timely question—how partial barriers and stickiness control asymptotic transport—and a validated power law for τ_c would be a valuable addition. However, the current evidence for the headline recurrence law is incomplete, so the paper's significance depends on additional numerical work rather than on the present manuscript alone.","major_comments":[{"comment":"The central result τ_c ∝ x_q^{-0.213} is not yet supported with adequate statistical rigor. For each x_q, the survival probability S(τ) is computed from a single orbit after discarding only 10^3 iterations, and the decay coefficient B is obtained by fitting S(τ) = P0 e^{Bτ} over the restricted range S(τ) > 0.08. No error bars are reported for |B|, no check of convergence with respect to orbit length or number of recurrences is given, and the sensitivity of the fitted exponent to the fit window and to the choice of recurrence box is not examined. With five points spanning five decades, a claimed exponent of 0.213 should be accompanied by a stability analysis; otherwise the value could be an artifact of the fit window rather than a transport law.","section":"Section 2.2, Eq. (43) and Fig. 4"},{"comment":"The interpretation of the exponential fit is in tension with the paper's own discussion of stickiness. The text states that the long-time tail of S(τ) carries information about sticky motion near KAM islands and cantori, but the fit deliberately discards precisely that tail (S(τ) ≤ 0.08). If the true recurrence-time distribution develops a power-law tail, the exponential decay constant is not a robust characteristic transport time and may not reflect the asymptotic trapping that the paper invokes to explain the slowdown. The manuscript should either justify the exponential model over the fitted range, quantify the tail separately, or present a complementary analysis of the long-time behavior before attributing τ_c ∝ x_q^{-0.213} to stickiness.","section":"Section 2.2, Eqs. (40)-(46)"},{"comment":"There is a direct inconsistency in the definition of the phase-space domain. The text repeatedly states that the radial variable satisfies 0 ≤ Ψ ≤ 1, with Ψ = 1 at the plasma boundary, yet the recurrence box is centered at Ψ = 5.0 and Fig. 1(a) shows Ψ ranging up to about 40, while Fig. 1's caption introduces an undefined parameter δ = 10^{-6}. The recurrence statistics are therefore obtained in a region that appears to lie outside the stated physical domain. This must be resolved: either the normalization of Ψ is different from what is written, or the recurrence box is misplaced, in which case the reported scaling may be specific to a particular (possibly non-physical) phase-space region.","section":"Section 2, Eq. (17) and Section 2.2, Eq. (40)"},{"comment":"The dynamic-scaling exponents are presented with standard errors from least-squares fits, but the universal collapse in Fig. 2(b) is assessed only visually. Given that the collapse is a central supporting result, the authors should provide a quantitative measure of collapse quality, for example the maximum or rms deviation of the rescaled curves from a master curve. This would also serve as a check on the claim that the scaling function is universal over the full range of x_q studied.","section":"Section 2.1, Fig. 3"}],"minor_comments":[{"comment":"The caption introduces the parameter δ = 10^{-6} without defining it anywhere in the text; if it is a leftover from a previous version, it should be removed or defined.","section":"Fig. 1 caption"},{"comment":"Several axis labels and tick labels are garbled or hard to read, for example '10°1' in Fig. 2 and the exponent notation in Fig. 4. The figures should be regenerated with clear mathematical notation, including negative exponents.","section":"Figs. 2 and 4"},{"comment":"The notation S(τ) = P0 e^{Bτ} with B negative is unconventional; using S(τ) = P0 e^{-τ/τ_c} would make the interpretation of τ_c as a decay time immediate and would avoid the absolute-value notation |B| introduced later.","section":"Section 2.2, Eq. (43)"},{"comment":"The ensemble average in Eq. (23) averages over the first N iterations of each trajectory, but the text does not state whether the initial transient is discarded for the ensemble transport analysis; for consistency with the recurrence analysis, this should be specified.","section":"Section 2.1, Eq. (23)"},{"comment":"Reference [9] is cited as 'Rechester and Rosenbluth 2020' in the bibliography, but the original work is much older; the citation should point to the original publication or clarify that it is a reprint.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is legitimate and is the main reason for the major-revision recommendation. The recurrence exponent is the paper's headline contribution, but it is currently based on a single-orbit, single-exponential fit over a restricted survival-probability range with no error bars or convergence checks. This is fixable within the manuscript's scope by adding multi-orbit statistics, reporting uncertainties, testing fit-window sensitivity, and examining the long-time tail. The dynamic-scaling part is considerably stronger and should be preserved. I would not recommend rejection, because the central idea is defensible and the requested evidence is numerical rather than conceptual."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper reports new scaling exponents for the conservative Tokamap, in particular the claim τ_c ∝ x_q^{-0.213} from recurrence statistics. The transport scaling part is solid: the three-regime behavior is clear, the exponents β=0.847, δ=-0.198, and z=-0.234 are measured, and the scaling law z=δ/β holds with excellent data collapse. That is a legitimate numerical result for a standard model.\n\nThe weakness is in the recurrence analysis that produces the headline exponent. The survival probability is fitted only over S(τ)>0.08 with a single exponential. The paper justifies this by saying the tail has statistical fluctuations, but in a mixed Hamiltonian system stickiness shows up precisely as a power-law tail. Discarding that tail means the fit measures the initial escape rate from the chaotic component, not the asymptotic trapping that the paper invokes to explain the slowdown. So the statement that 'long-time dynamics are increasingly dominated by stickiness' is not actually supported by the fitted observable. On top of that, each survival curve comes from a single orbit, with no error bars on B, and the 0.213 exponent has no uncertainty. A few checks - multiple orbits, different fit windows, different box locations - would either confirm it or expose it as a window artifact.\n\nThere are also internal inconsistencies in the manuscript. The text says Ψ lies in [0,1] with Ψ=1 at the plasma boundary, but Figure 1(a) shows Ψ up to 40 and the recurrence box is centered at Ψ=5.0. Either the normalization is wrong or the reported parameter regime is mis-specified. The Figure 1 caption also introduces an undefined δ=10^-6. These are fixable but need to be corrected.\n\nThe paper should be sent to peer review: the scaling transport results are worth checking, and the recurrence claim needs referee scrutiny. The authors should be asked for error bars, multiple orbits, and robustness tests, plus a clean definition of the radial coordinate. As it stands, the central claim is plausible but not established.","headline":"Plausible scaling exponents for the conservative Tokamap, but the headline recurrence scaling rests on a single exponential fit that excludes the sticky tail.","tokens_in":14127,"tokens_out":3307,"would_cite":false,"duration_ms":29939,"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":"In the conservative Tokamap, lowering magnetic shear slows chaotic transport as $\\tau_c\\propto x_q^{-0.213}$.","keywords":["Tokamap","symplectic mapping","magnetic field-line transport","Hamiltonian chaos","Poincaré recurrence statistics","stickiness","dynamic scaling","magnetic shear"],"falsifier":"Repeat the recurrence measurement with many independent chaotic orbits at $x_q=10^{-4}$ and $x_q=10^{-8}$, fitting $|B|$ to each survival curve over $S(\\tau)>0.08$ after the same $10^3$-iteration transient; if the orbit-to-orbit spread in $|B|$ is comparable to the factor of $10^4$ separating the two parameter values, or if the early tail is visibly non-exponential, then $\\tau_c\\propto x_q^{-0.213}$ would not survive.","tokens_in":12989,"feed_emoji":"🧲","tokens_out":12659,"duration_ms":106365,"temperature":0.7,"pith_summary":"The paper aims to establish that long-time magnetic field-line transport in a tokamak is governed by the geometry of invariant structures, not just by local chaos. Working with the conservative Tokamap, an exact area-preserving map of field-line motion, the authors show that lowering the magnetic-shear parameter $x_q$ systematically slows transport. The quantitative core is the algebraic law $\\tau_c \\propto x_q^{-0.213}$: the characteristic time between returns to a fixed phase-space region grows as the shear decreases. Because the dynamics stay chaotic throughout, the slower transport is attributed to stickiness near KAM islands, resonance chains, and cantori. If correct, the result connects microscopic phase-space structure to macroscopic confinement-relevant transport times and implies that Lyapunov instability alone cannot predict transport efficiency.","feed_headline":"Chaotic transport in the Tokamap obeys a fixed power law","feed_subtitle":"Lower magnetic shear traps trajectories longer near islands and cantori, slowing return times.","key_machinery":"The argument is carried by the conservative Tokamap itself, the exact symplectic map defined by $\\Psi_{n+1}=\\tfrac{1}{2}(P_n+\\sqrt{P_n^2+4\\Psi_n})$ and $T_{n+1}=T_n+1/q(\\Psi_{n+1})-\\tfrac{x_L}{4\\pi^2}\\cos(2\\pi T_n)/(1+\\Psi_{n+1})^2$ (mod 1), with $P_n=\\Psi_n-1-\\tfrac{x_L}{2\\pi}\\sin(2\\pi T_n)$ and the safety-factor profile $q(\\Psi)=x_q(1+\\Psi^2)^2$. This map preserves phase-space area exactly, so the invariant tori, islands, and cantori are genuine features of the dynamics rather than numerical artifacts. The transport analysis rests on two further pieces of machinery: the homogeneous scaling hypothesis $\\Psi(n,x_q)=\\ell F(\\ell^a n,\\ell^b x_q)$, which yields the exponent relation $z=\\delta/\\beta$, and the survival probability $S(\\tau)=P(T>\\tau)$, whose exponential decay constant $B$ in the interval $S(\\tau)>0.08$ defines the characteristic transport time $\\tau_c=1/|B|$.","core_discovery":"The central discovery is that the conservative Tokamap's mixed phase space produces a scale-invariant transport law with a well-defined long-time exponent. For ensembles of $10^4$ trajectories, the mean radial coordinate grows as $\\Psi\\propto n^{\\beta}$, saturates as $\\Psi_{\\mathrm{sat}}\\propto x_q^{\\delta}$, and crosses over at $n_x\\propto x_q^{z}$, with $\\beta=0.847(9)$, $\\delta=-0.198(1)$, and $z=-0.234(1)$; the homogeneous scaling assumption predicts $z=\\delta/\\beta$, and the measured curves collapse onto one universal function under the rescalings $n\\to n/x_q^z$ and $\\Psi\\to\\Psi/x_q^{\\delta}$. At long times, the survival probabilities $S(\\tau)$ of returns to a fixed box shift monotonically toward longer recurrence times as $x_q$ decreases from $10^{-4}$ to $10^{-8}$. Fitting the early portion $S(\\tau)>0.08$ to $S(\\tau)=P_0 e^{B\\tau}$ gives $|B|\\propto x_q^{0.213}$, hence $\\tau_c=1/|B|\\propto x_q^{-0.213}$; the paper interprets this as progressively dominant stickiness generated by KAM islands, resonance chains, and cantori.","pith_inferences":["Inference: The paper's exponent $0.213$ is measured with a single long orbit per parameter value; an ensemble of independent orbits would show whether the survival probability and its decay constant are orbit-independent or only typical.","Inference: Varying the recurrence box's location and size at fixed $x_q$ would test whether $\\tau_c\\propto x_q^{-0.213}$ is a global property of the chaotic sea or a local diagnostic tied to the chosen window.","Inference: A quantitative link to cantorus flux formulas is left implicit; computing the flux across the relevant partial barriers and comparing it with $1/\\tau_c$ would place the exponent on firmer mechanistic ground.","Inference: For reversed-shear or dissipative Tokamaps, the same analysis could reveal a different exponent or a non-exponential survival tail, since separatrix reconnection and crisis-induced intermittency change the sticky structures; the paper flags these systems as future work but does not predict the outcome."],"forward_implications":["The scaling law $z=\\delta/\\beta$ means only two of the three transport exponents are independent; the observed collapse of the transport curves onto one master function is the direct test of this prediction.","Because $\\tau_c\\propto x_q^{-0.213}$, reducing magnetic shear over the investigated range slows the exploration of phase space by a power law rather than by an abrupt transition.","Positive largest Lyapunov exponents coexist with progressively slower transport, so local exponential instability is not a sufficient measure of how efficiently field lines explore the plasma.","The same combination of phase-space portraits, Lyapunov exponents, scaling collapse, and recurrence statistics can be applied to other Tokamap variants and to generic Hamiltonian systems with mixed phase space."],"supporting_citations":[{"why":"This reference defines the Tokamap and its symplectic construction, the model whose transport is analyzed throughout the paper.","marker":"[34]"},{"why":"It supplies the survival-probability formulation used to quantify recurrence times and long-time transport.","marker":"[49]"},{"why":"It provides the cantorus-as-partial-barrier picture that the paper uses to explain stickiness and slow recurrence.","marker":"[22]"},{"why":"It documents sticky motion near island boundaries, the mechanism invoked for the algebraic growth of transport times.","marker":"[47]"},{"why":"It reviews anomalous diffusion and stickiness in Hamiltonian systems, supporting the interpretation of the slowly decaying survival curves.","marker":"[23]"},{"why":"It provides the dynamic-scaling hypothesis and homogeneous-function framework the paper adapts to derive the exponent relation and the collapse.","marker":"[44]"},{"why":"It characterizes transport in mixed Hamiltonian phase space and justifies the need to go beyond Lyapunov exponents for long-time predictions.","marker":"[19]"}],"fun_headline_variants":["Tokamap transport follows one universal scaling law","Magnetic shear sets the pace of chaotic transport","A single exponent controls tokamak field-line escape","Sticky islands slow magnetic field-line transport","Universal collapse: tokamap chaos obeys one law"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the early tail of the return-time distribution from one long chaotic orbit, after a $10^3$-iteration transient, is a converged single exponential; if the survival probability is not converged or is genuinely non-exponential on $S(\\tau)>0.08$, the fitted $B$ and the exponent $0.213$ would not be a robust transport law.","fun_headline_variants_meta":{"raw":{"variants":["Tokamap transport follows one universal scaling law","Magnetic shear sets the pace of chaotic transport","A single exponent controls tokamak field-line escape","Sticky islands slow magnetic field-line transport","Universal collapse: tokamap chaos obeys one law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000169,"raw_usage":{"total_tokens":1294,"prompt_tokens":1004,"completion_tokens":290,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":620,"completion_tokens_details":{"reasoning_tokens":217}},"tokens_in":620,"tokens_out":290,"duration_ms":2965,"temperature":1.0,"reasoning_tokens":217,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T13:36:30.581948+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the recurrence measurement with many independent chaotic orbits at $x_q=10^{-4}$ and $x_q=10^{-8}$, fitting $|B|$ to each survival curve over $S(\\tau)>0.08$ after the same $10^3$-iteration transient; if the orbit-to-orbit spread in $|B|$ is comparable to the factor of $10^4$ separating the two parameter values, or if the early tail is visibly non-exponential, then $\\tau_c\\propto x_q^{-0.213}$ would not survive.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This reference defines the Tokamap and its symplectic construction, the model whose transport is analyzed throughout the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the survival-probability formulation used to quantify recurrence times and long-time transport."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the cantorus-as-partial-barrier picture that the paper uses to explain stickiness and slow recurrence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It documents sticky motion near island boundaries, the mechanism invoked for the algebraic growth of transport times."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It reviews anomalous diffusion and stickiness in Hamiltonian systems, supporting the interpretation of the slowly decaying survival curves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the dynamic-scaling hypothesis and homogeneous-function framework the paper adapts to derive the exponent relation and the collapse."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It characterizes transport in mixed Hamiltonian phase space and justifies the need to go beyond Lyapunov exponents for long-time predictions."}],"review_version":1}