{"id":"9b2833c8-de95-40a0-ae19-0bb5a624d026","arxiv_id":"2608.05936","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A 2D Ising system driven by a periodic magnetic field at its critical point obeys universal dynamic scaling with variables tau = t/P and sigma = A P^kappa.","lead":"This paper shows that a 2D Ising magnet driven by a periodic magnetic field at its critical temperature responds according to a universal scaling law involving the field amplitude, period, and elapsed time. The result provides a general framework for periodically driven systems near continuous phase transitions, which could help experiments measure dynamic critical exponents.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scaling collapse is tested only over P=100-300 along the line A=sigma P^{-kappa}, so the load-bearing assumption that P enters only through sigma remains under-tested.","rationale":"The reader's weakest-assumption analysis identifies the same point: the scaling theory in Sec. III B combines W_a and W_p into sigma by assuming that the amplitude scales like the static field and the period scales like time. My stress-test agrees with this identification as the most load-bearing element of the paper. The paper's own language -- 'a reasonable hypothesis' -- is the only support for the crucial step, and the numerical evidence is a collapse along the constraint A=sigma P^{-kappa} over a modest P range. That evidence is consistent with the hypothesis but is not a strong falsification test. A separate relevant variable, if present, would destroy the universal two-variable scaling claim, which is the core of the paper. I do not see an internal inconsistency in the scaling argument as stated; the issue is under-determination of the hypothesis by the data. The peripheral exponent u fitted from four points, the lack of error bars, and the absence of a test with a different dynamics are secondary: they weaken confidence but do not target the logical hinge as directly. The proposed wider-P collapse test would settle the main concern in a way that is both computationally feasible and directly tied to the claim. Since the reader already issued a CONDITIONAL verdict and my concern does not move the verdict to a different category, UNCHANGED is the appropriate recommendation. If the wider-P test fails, the verdict should move toward REJECT; if it passes, the paper could be accepted with the current evidence supplemented by quantitative error estimates.","tokens_in":16739,"tokens_out":8816,"duration_ms":111638,"concrete_test":"Run the periodic-driving protocol at sigma=1 for P=50,100,200,400,800, setting A=P^{-kappa} with kappa=0.8653, on lattices large enough that finite-size effects are negligible (e.g., L=400 for P=50 and L=800 for P=800), with at least 500 independent trajectories. At fixed tau=1, 3, and in the stationary regime, compute the rescaled magnetization A^{-1/15} M(t) and the corresponding bootstrap errors. Quantify the collapse residual as the maximum or RMS spread of the rescaled curves across the five P values. If the spread grows systematically with P and exceeds about 2-3 statistical errors, P is not fully absorbed by sigma and the central scaling hypothesis is falsified. If the spread remains at the noise level over this factor of 16 in P, the separate-relevant-variable concern is refuted for practical purposes.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is anchored to the hypothesis in Sec. III B, Eqs. (8)-(15), that the only relevant combination of the driving parameters is sigma = A P^{y_h/z} and tau = t/P. This is explicitly introduced as a 'reasonable hypothesis' rather than derived from a dynamic RG calculation. The numerical verification, however, is restricted to the one-dimensional curve A = sigma P^{-kappa}: for sigma=1 the periods are P=100,200,300; for sigma=1/2 the periods are P=100,200; for sigma=10 the periods are P=200,400. This tests the predicted scaling form along a single line in the (A,P) plane, which is necessary but not sufficient. If the periodic drive generated an additional relevant scaling variable not captured by sigma -- for example a Floquet frequency variable, a stroboscopic phase variable, or an independent dependence on P/L^z -- then two runs with the same sigma but very different P should fail to collapse. With the current factor of 2-3 in P, and with the paper's own estimate that scaling corrections decay only as P^{-0.92}, a weak extra P-dependence could be hidden in the visually assessed collapse. The large-time relaxation time tau_s ~ sigma^{-u} with u about 3.9 is a separate, non-standard exponent; its presence shows that the scaling functions contain a slow variable, but it does not by itself invalidate the two-variable form. Still, the central claim would fail if a third scaling variable exists, and the paper does not provide a test that could detect it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the 2D Ising model at T_c under purely relaxational (Metropolis) dynamics and a periodic magnetic field h(t)=-A cos(2πt/P). It proposes a dynamic scaling theory in which the only relevant combinations are σ=A P^{y_h/z} and τ=t/P, leading to M ≈ A^{1/15} M(σ,τ) and E_s ≈ A^{8/15} E_s(σ,τ). The authors support the theory with Monte Carlo simulations for several values of σ and P, checks against the σ→∞ and σ→0 limits, a square-wave protocol, and a periodically varying temperature protocol. They also study the transient approach to the stationary state and report a non-standard divergence exponent u≈3.9 for the transient time scale.","tokens_in":17083,"tokens_out":10938,"duration_ms":121890,"significance":"If established, the result is a clean extension of dynamic scaling to periodic driving, with the central exponents y_h and z taken from independent equilibrium results, so the collapse is essentially parameter-free apart from the peripheral transient exponent u. The paper includes useful consistency checks, including the σ→∞ and σ→0 limits and a square-wave protocol that are not obtained by fitting the collapse exponents. The main limitations are that the numerical test of the two-variable scaling form covers only a narrow range of P at fixed σ, and the collapse claims are not accompanied by displayed error bars; the universality claim also rests on a single model and a single dynamics. These issues are addressable with additional analysis rather than requiring a change in the theoretical framework.","major_comments":[{"comment":"The central claim that the driving enters only through σ=A P^{y_h/z} and τ=t/P is tested only along the one-parameter family A=σ P^{-κ} with P spanning a factor of 2-3 for each σ (for example, σ=1 uses P=100,200,300; σ=1/2 uses P=100,200; σ=10 uses P=200,400). Given the paper's own estimate that scaling corrections decay as P^{-0.92}, a weak additional P-dependence could be hidden in the visually assessed collapse. I request a more stringent test: at least one σ with P varied by an order of magnitude, or a quantitative correction-to-scaling analysis (for instance, fitting the data with an additional P^{-0.92} term and showing its amplitude is consistent with zero). Without this, the two-variable scaling form is plausible but not established.","section":"Sec. III B, Eqs. (8)-(15), and Figs. 2, 4, 5, 11, 12"},{"comment":"The collapse evidence is presented without error bars, and statements such as \"they collapse onto single curves\" and \"corrections are smaller than the statistical errors\" cannot be independently checked. Please display representative error bars, or provide a residual plot with errors, for at least the main collapse figures. This is important because the visual collapse is the principal numerical support for Eqs. (14) and (15).","section":"Sec. IV A, Figs. 1, 2, 4, 5, 10, 11, 12"}],"minor_comments":[{"comment":"The exponent u=3.9(2) is inferred from only four τ_s estimates with relative uncertainties of 10-20%, each obtained from a fit over a limited range of ln|M_a|. The reported uncertainty likely underestimates systematic errors; if this result is retained, a more complete analysis including fit-range dependence and model selection is needed.","section":"Sec. IV A, Eq. (23) and Fig. 6"},{"comment":"There are several typographical errors: \"the evolution of a a single system\" in Sec. IV B, \"broken be the starting condition\" in the Conclusions, and \"non-exaustive\" in the Introduction.","section":"Sec. IV B and Conclusions"},{"comment":"The oscillation amplitude is quoted as A_m≈0.17 for the rescaled magnetization and as A_m≈0.1254 for the raw magnetization at P=200, σ=1; the notation should state explicitly which quantity is being reported in each case.","section":"Sec. IV A, around Fig. 3 and Sec. IV B"},{"comment":"The notation E_s(t) for the subtracted bond-energy density in the magnetic-driving case and E_se(t) in the temperature-driving case is similar; a more distinct symbol would help the reader avoid confusion.","section":"Sec. V, Eqs. (6) and (30)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the central idea is interesting. My main concern is that the numerical evidence for the two-variable scaling form is narrower than the claim, and the absence of error bars makes the collapse assessment difficult. These concerns are addressable in revision; I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is an exploratory study, but a serious and largely sound one. The genuinely new result is the prediction that a periodic magnetic drive at Tc produces dynamic scaling in terms of sigma = A P^{y_h/z} and tau = t/P, with universal scaling functions for the magnetization and energy density. The numerical collapse across P = 100-300 (up to 400 for sigma = 10) looks convincing, and the sigma -> 0 and sigma -> infinity limits match the predicted forms. Two extras strengthen the case: the square-wave protocol gives the same scaling, and the temperature-driving protocol uses a different combination, B P^{y_t/z}, so the RG form is not tailored to one specific drive. Exponents are all taken from independent published results, so there is no circularity in the central check.\n\nSoft spots are real but, I think, not disqualifying. The hypothesis in Sec. III B is introduced as a 'reasonable hypothesis' and is not derived from a dynamic RG calculation. The numerical test is confined to the curve A = sigma P^{-kappa}; that is necessary but not sufficient. As the stress-test note says, a factor of 2-3 in P, with corrections decaying as P^{-0.92}, could hide a weak extra P-dependence. I do not find that likely: the square-wave and temperature results make a conspiracy less plausible, and the sigma -> infinity limit is quantitatively consistent with equilibrium. Still, a referee should ask for a wider P range or a two-parameter check before the two-variable form is fully established. The lack of error bars on the main figures is annoying, though the multi-size checks provide some context. The fitted exponent u for the transient time scale is peripheral, and u = 3.9(2) comes from four points; the authors are appropriately cautious about its uncertainty.\n\nWho is this for? People working on out-of-equilibrium critical dynamics and dynamic scaling. It is a useful extension of a known framework, not a major reorganization. I would send it to a serious referee; the central claim is plausible, the evidence is decent, and the potential gap is well-defined. This deserves refereeing, not desk rejection.","headline":"A plausible and mostly convincing extension of dynamic scaling to periodically driven Ising systems; the main gap is that the load-bearing scaling variable is only tested along one line in parameter space.","tokens_in":17582,"tokens_out":2387,"would_cite":true,"duration_ms":27274,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B20","82B27","82C05"],"pacs":["64.60.Ht","75.10.Hk","05.10.Ln"],"model":"deepseek-v4-flash","headline":"Periodic magnetic sweeps across a 2D Ising critical point obey a universal dynamic scaling law.","keywords":["dynamic scaling","periodic driving","Ising model","critical dynamics","renormalization group","model A dynamics","Z2 symmetry breaking","Monte Carlo simulation"],"falsifier":"A run at fixed $\\sigma=1$ with $P=200,300,\\ldots,1000$ in the thermodynamic limit should show no residual $P$ dependence of $A^{-\\zeta}M(\\tau)$; a systematic drift with $P$ would falsify the scaling form. A sharper check is to measure $\\kappa$ from the collapse and compare it with $15/8$ divided by an independent determination of $z$; because the paper's exponents are tied to the equilibrium dynamic universality class, any effective $\\kappa$ that varies with $\\sigma$ or $P$ would break the central identity.","tokens_in":16557,"feed_emoji":"🧲","tokens_out":13005,"duration_ms":118725,"temperature":0.7,"pith_summary":"This paper asks what happens when a critical ferromagnet is pushed back and forth through its transition by a periodic magnetic field, and answers that the response is not chaotic but collapses onto a universal curve. At the critical temperature of the two-dimensional Ising model with purely relaxational dynamics, the magnetization and the subtracted bond-energy density depend on the drive only through the combination $\\sigma=AP^\\kappa$, with $\\kappa=y_h/z\\approx 0.865$, and on time only through $\\tau=t/P$. Consequently runs with very different amplitudes and periods fall on the same scaling functions, which oscillate in sync with the field, lag behind it by a phase that shrinks as $\\sigma$ grows, and settle into a stationary regime at large $\\tau$. The paper argues by renormalization-group scaling and confirms the collapse in local-update Monte Carlo simulations, including square-wave drives and periodic temperature drives.","feed_headline":"Periodic magnetic sweeps obey universal scaling at Ising criticality","feed_subtitle":"Magnetization and bond energy collapse onto curves controlled by A P^κ and t/P, with κ≈0.865.","key_machinery":"The load-bearing object is the renormalization-group scaling variable $\\sigma=A P^\\kappa$, obtained by treating the field amplitude $A$ as a static magnetic field with RG dimension $y_h$ and the period $P$ as a time scale with dynamic dimension $z$. In finite size this means $W_a=A L^{y_h}$ and $W_p=P L^{-z}$; in the thermodynamic limit these two variables combine into $\\sigma=W_a W_p^{y_h/z}=A P^\\kappa$, while $\\tau=t/P$ is the ratio of the two time scales. This construction carries the argument because it turns a two-parameter family of protocols into one-parameter families of universal functions, and it fixes the exponents $\\zeta=y_\\phi/y_h=1/15$ and $\\varepsilon=y_e/y_h=8/15$ that make data from different $A$ and $P$ collapse.","core_discovery":"For the 2D Ising model at $T_c$ evolving under purely relaxational (model-A) dynamics in a field $h(t)=-A\\cos(2\\pi t/P)$, the infinite-volume magnetization obeys $M(t,A,P)\\approx A^{1/15}\\mathcal{M}(\\sigma,\\tau)$ and the subtracted bond-energy density obeys $E_s(t,A,P)\\approx A^{8/15}\\mathcal{E}_s(\\sigma,\\tau)$, with $\\tau=t/P$ and $\\sigma=A P^\\kappa$, $\\kappa=y_h/z=0.8653(4)$. The scaling functions are universal for the model-A equilibrium dynamic universality class; they are synchronized with the drive (period one in $\\tau$), show a phase delay of about 0.2 to 0.25 periods that vanishes as $\\sigma\\to\\infty$, and at large $\\tau$ reach a stationary state whose period-averaged magnetization is zero, so the $Z_2$ symmetry broken by the initial condition is recovered. The same scaling theory, with $\\rho=B P^{y_t/z}$ and $y_t/z\\approx 0.4615$, describes periodic temperature variation at zero field, where the susceptibility scales as $B^{-7/4}$ times a universal function. The argument is expected to carry to higher-dimensional Ising systems and, with caveats about energy injection, to quantum transitions.","pith_inferences":["A direct test of universality would be to run the same protocol in a three-dimensional Ising or $O(N)$ model and verify collapse at $\\kappa=y_h/z\\approx 1.226$; the paper predicts this but does not simulate it.","The divergent relaxation time $\\tau_s\\sim\\sigma^{-4}$ resembles the slow dynamics seen near the dynamic phase transition of driven Ising magnets, and connecting the two could give an independent route to measure $u$.","Near a zero of the field the out-of-equilibrium interval shrinks as $\\Delta\\tau\\sim\\sigma^{-1/(1+\\kappa)}$, so measuring the phase delay versus $\\sigma$ would probe whether quench-scaling physics controls the synchronization.","In a thin ferromagnetic film, sweeping frequency and amplitude while recording AC susceptibility should show the predicted data collapse, making the phase delay and the stationary large-$\\tau$ state directly observable."],"forward_implications":["At fixed $\\sigma$, changing $P$ from 100 to 300 leaves the rescaled magnetization and bond energy unchanged in the thermodynamic limit; the collapse is the direct signature of the scaling law.","In the large-$\\sigma$ limit the magnetization approaches the equilibrium curve $\\mathrm{sgn}[h(t)]\\,c\\,|\\cos(2\\pi\\tau)|^{1/15}$ with $c=1.00268751$, so the periodic protocol becomes effectively adiabatic.","For $\\sigma\\lesssim 2$ the system takes many periods to forget its initial magnetization; the rescaled relaxation time grows as $\\tau_s\\approx a\\,\\sigma^{-u}$ with $u\\approx 3.9$, so $Z_2$ symmetry is recovered only after a number of cycles that diverges as $\\sigma\\to 0$.","Square-wave magnetic drives and periodic temperature drives obey the same dynamic scaling structure with their own universal functions; in the temperature case the susceptibility scales as $B^{-7/4}$ and the properly subtracted energy as $B^{-1}$.","In the three-dimensional Ising universality class the predicted exponent is $\\kappa\\approx 1.226$, and for quantum Ising transitions $z=1$ would set a different scaling regime; these numbers give concrete targets for future tests."],"supporting_citations":[{"why":"Classifies purely relaxational (model-A) critical dynamics, defining the dynamic universality class the scaling theory targets.","marker":"[8]"},{"why":"Supplies the RG dimensions $y_h$, $y_\\phi$, $y_e$ and $z$, plus the finite-size scaling theory this paper extends to periodic driving.","marker":"[11]"},{"why":"Provides the exact critical temperature and the static exponents $\\nu=1$, $\\eta=1/4$ that enter $y_h=15/8$.","marker":"[45]"},{"why":"Defines the local spin-flip update rule used in the Monte Carlo evolution.","marker":"[47]"},{"why":"Gives one of the precise dynamic-exponent estimates used to set $z=2.167(1)$ and hence $\\kappa\\approx 0.865$.","marker":"[48]"},{"why":"Supplies the equilibrium magnetization amplitude $c=1.00268751$ used in the $\\sigma\\to\\infty$ asymptotic prediction.","marker":"[58]"}],"fun_headline_variants":["Universal scaling for Ising magnet under periodic field drive","Periodic magnetic drive: Ising dynamics obey universal scaling","Ising criticality: periodic field yields scaling with A P^0.865","2D Ising under periodic field: universal dynamic scaling revealed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire scaling structure rests on the hypothesis that the amplitude $A$ enters only through the static-field scaling variable and the period only through the time scaling variable, so their combined effect is captured by $\\sigma=A P^{y_h/z}$; a periodic drive that creates an independent relevant perturbation would destroy the collapse.","fun_headline_variants_meta":{"raw":{"variants":["Universal scaling for Ising magnet under periodic field drive","Periodic magnetic drive: Ising dynamics obey universal scaling","Ising criticality: periodic field yields scaling with A P^0.865","2D Ising under periodic field: universal dynamic scaling revealed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000242,"raw_usage":{"total_tokens":1593,"prompt_tokens":1079,"completion_tokens":514,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":695,"completion_tokens_details":{"reasoning_tokens":442}},"tokens_in":695,"tokens_out":514,"duration_ms":5733,"temperature":1.0,"reasoning_tokens":442,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T20:44:52.073182+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A run at fixed $\\sigma=1$ with $P=200,300,\\ldots,1000$ in the thermodynamic limit should show no residual $P$ dependence of $A^{-\\zeta}M(\\tau)$; a systematic drift with $P$ would falsify the scaling form. A sharper check is to measure $\\kappa$ from the collapse and compare it with $15/8$ divided by an independent determination of $z$; because the paper's exponents are tied to the equilibrium dynamic universality class, any effective $\\kappa$ that varies with $\\sigma$ or $P$ would break the central identity.","supporting_citations":[{"cited_title":"Metropolis, A","cited_arxiv_id":null,"evidence_quote":"Defines the local spin-flip update rule used in the Monte Carlo evolution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives one of the precise dynamic-exponent estimates used to set $z=2.167(1)$ and hence $\\kappa\\approx 0.865$."},{"cited_title":"Caselle and M","cited_arxiv_id":null,"evidence_quote":"Supplies the equilibrium magnetization amplitude $c=1.00268751$ used in the $\\sigma\\to\\infty$ asymptotic prediction."}],"review_version":1}