{"id":"a372be12-430d-4cf0-ae59-b26ede69b3ca","arxiv_id":"2412.09343","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a driven three-state Ising model, the order-disorder transition splits into two distinct points with different scaling, one continuous (β=1) and one discontinuous.","lead":"A driven three-state Ising model coupled to two heat baths at different temperatures and pushed by opposite driving forces develops two different order-disorder transition points, one for each dominant spin state. The transitions have different characters and thermodynamic properties, including a phase that runs as a low-dissipation heat engine near ideal efficiency.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (2) itself makes m=0 linearly stable for ε<ε1c, so ε1c is a spinodal/metastability limit rather than a continuous critical point; the claimed β=1 exponent is not a true critical exponent.","rationale":"The reader's weakest_assumption already identified the possibility that ε1c is a spinodal and β=1 is not a true critical exponent. My analysis shows this is not merely a possibility but is forced by Eq. (2) combined with the stated sign a>0: for ε<ε1c the symmetric state is linearly stable. This makes the identification of ε1c as a continuous phase transition internally inconsistent with the model's own deterministic dynamics. The strongest claim rests on this identification: the headline novelty is that one transition is continuous with a non-standard exponent while the other is discontinuous. If ε1c is instead a spinodal, the two transitions differ only in the location and shape of metastable branches, and the 'unique scaling laws' for phase A (δ=2) are scaling of a metastable branch rather than of a thermodynamic phase transition. The all-to-all exact solution and the square-lattice simulations may still show two transition points, but the qualitative claim of a continuous nonequilibrium critical point with β=1 is unsupported. This is a decisive issue, not a cosmetic one, because the title and abstract promise a splitting of 'phase transitions' and a critical exponent distinct from β=1/2. I therefore recommend REJECT as submitted, with the possibility of a major revision that reinterprets ε1c as a spinodal and removes the critical-exponent claims. Agreement with reader: the reader's weakest assumption is the same concern, though the reader treated it as addressable; I consider it fatal to the central claim as stated.","tokens_in":16363,"tokens_out":9745,"duration_ms":99512,"concrete_test":"For the parameters of Fig. 1 (β1=2, β2=1, F=2), solve the full all-to-all steady-state mean-field equations (without truncating in m) and compute the Jacobian eigenvalue at the symmetric fixed point (m=0, q=2/3) for a range of ε below ε1c. If the eigenvalue is negative, the disordered state is stable on the ordered side, confirming ε1c as a spinodal. Complement with square-lattice Gillespie simulations (e.g., N=64 or 128) at ε just below ε1c starting from uniform m=0 and from m>0 initial conditions; if both initial conditions persist on simulation timescales (or the steady-state order-parameter distribution is bimodal), the transition is not continuous and β=1 is not a critical exponent.","verdict_should_be":"REJECT","load_bearing_attack":"In 'Emergence of phase transitions...' the expansion dm/dt ≈ a(ε−ε1c)m + b m^2 + c m^3 (Eq. 2) is used to identify ε1c as a continuous transition with m ∼ (ε1c−ε)/b, i.e. β=1. But the text states a>0. Therefore for every ε<ε1c (the ordered side), the linear coefficient a(ε−ε1c) is negative and m=0 is a linearly stable fixed point. The ordered branch with m>0 must then coexist with a stable disordered state, creating a bistable/spinodal region of exactly the kind the paper reserves for ε2c ('Unlike a critical transition, the latter case features a spinodal region...'). At a genuine continuous order-disorder transition the disordered state should be unstable on the ordered side, so the order parameter grows from zero as the stable state is selected. The transcritical branch endpoint at ε1c is at best a limit of metastability of the m>0 phase; the exponent β=1 describes the vanishing of a metastable branch, not a critical exponent. Since the claimed δ=2 scaling of entropy production, power fluctuations, and efficiency is derived from this same expansion (δ=2β), those scaling laws inherit the same mischaracterization. Thus the central distinction between a 'continuous critical' transition at ε1c and a 'discontinuous' transition at ε2c is not established; the splitting may instead be between two spinodal points.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a three-state Ising model in which spins are coupled to two thermal baths with opposite driving forces, and it claims that the order-disorder transition splits into two distinct transition points, one for each ordered phase. For the all-to-all interaction case, the authors derive an effective Landau-type equation for the magnetization near the symmetric state and identify the m>0 transition as a continuous critical point with exponent β=1, while the m<0 transition is said to be discontinuous and to feature a spinodal region. From the same expansion they obtain scaling laws δ=2 for entropy production, power fluctuations, and efficiency, and they argue that the m>0 phase is less dissipative and supports efficient heat-engine operation. The same scenario is asserted to hold for a square lattice on the basis of Gillespie simulations with N=102.","tokens_in":16694,"tokens_out":13698,"duration_ms":127739,"significance":"The claimed splitting of a single order-disorder transition into two phase-transition points with different classifications and distinct thermodynamic scaling would be a novel nonequilibrium phenomenon and would be of interest to the statistical-physics and stochastic-thermodynamics communities. The all-to-all analysis is a genuine strength: the master equation is stated, the expansion coefficients are provided in the Supplemental Material, and no parameter is fitted to data, making the mean-field results internally checkable. The prediction that phase A is less dissipative and has lower power fluctuations than phase B is falsifiable by simulation and experiment. However, the central conceptual claim—that ε1c is a continuous critical point with a true exponent β=1—is undermined by the linear stability of Eq. (2), as detailed below.","major_comments":[{"comment":"Equation (2) is not consistent with the claim that ε1c is a continuous critical point. The text states a>0 and the ordered phase corresponds to ε<ε1c; hence for all ε<ε1c the linear coefficient a(ε−ε1c) is negative, so m=0 is a linearly stable fixed point on the ordered side of the transition. In a Landau equation of the form dm/dt = L m + b m^2 + c m^3 with L<0, the nonzero branch that vanishes linearly at L=0, m* ≈ −L/b, has linearized growth rate f'(m*) ≈ −L > 0, so it is an unstable (separatrix) branch, not the stable order parameter of phase A. The stable ordered branch, if present, must terminate instead at a saddle-node spinodal point. Thus ε1c is at best a transcritical bifurcation or limit of metastability of the m>0 branch, and the exponent β=1 is not a critical exponent for the stable order parameter. The subsequent scaling laws for entropy production, power fluctuations, and efficiency (δ=2β=2) are derived from this same expansion and inherit the same mischaracterization. This also contradicts the paper's own distinction between a critical transition and a spinodal region: the coexistence of a stable disordered state and an ordered branch on the same side of ε1c is exactly the kind of spinodal behavior the text reserves for ε2c. The central claim of a continuous critical transition at ε1c therefore needs to be reanalyzed; if the stability calculation is correct, the two transitions are both spinodal-like.","section":"Emergence of phase transitions..., Eq. (2)"},{"comment":"The square-lattice results, which are used to claim robustness beyond the all-to-all case, are based on Gillespie simulations for N=102 with no error bars and no finite-size scaling analysis. A single small lattice size cannot support the quantitative claims that the two-transition scenario, the scaling exponents, and the thermodynamic asymmetries persist for short-range interactions. In particular, the caption of Fig. 1(e) itself notes a finite-size jump from m≈1 to −1, which signals strong finite-size effects. Please provide error bars, at least two or three system sizes, and a finite-size-scaling or extrapolation analysis, or explicitly restrict the claims to the all-to-all model.","section":"Emergence of phase transitions..., Figs. 1e, 2a-b, 3"}],"minor_comments":[{"comment":"In the sentence defining ⟨σc⟩, the phrase \"for ϵ1≥ ϵ1c\" appears to be a typo for \"ε ≥ ε1c\".","section":"Thermodynamics section"},{"comment":"The expression \"γ(A)_P −γ2(c)_P\" should read \"γ(A)_P − γ_P^{(c)}\" for consistency with the notation introduced for the critical value of the variance.","section":"Scaling of power fluctuations"},{"comment":"The displayed formula for ε1c is typeset with a leading negative sign and a denominator that are difficult to parse; please check the formatting so that readers can verify the sign and value.","section":"Eq. for ε1c in the main text"},{"comment":"The expression for the coefficient c in Eq. (B2) contains several unbalanced parentheses and repeated opening brackets, which make the formula effectively unreadable and prevent verification; please rewrite it cleanly, possibly with intermediate definitions.","section":"Supplemental Material, coefficient c"},{"comment":"The finite-size jump from m≈1 to −1 is not explained; a brief statement about the simulation protocol (initial conditions, number of runs) would help the reader assess the reliability of the square-lattice data.","section":"Fig. 1(e) caption"}],"recommendation":"major_revision","confidential_remarks":"The main issue is the internal inconsistency between the stated sign of a>0 in Eq. (2) and the claim that ε1c is a continuous critical point. This goes beyond a presentation problem: it affects the interpretation of the central scaling laws. If the authors can show that the stable order parameter actually vanishes at ε1c with β=1 (e.g., because the projection onto the m dynamics changes the effective sign of a, or because the relevant branch is not the one considered here), the paper could be acceptable. Alternatively, reframing both transitions as spinodal/metastability limits would remove the critical-exponent claims and weaken the novelty but might still present a valid phenomenon. The all-to-all exact solution and the absence of fitted parameters are positive features that deserve credit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core finding here is worth taking seriously: in a three-state driven Ising model, the order-disorder transition splits into two separate points, one for each ordered phase, with different thermodynamic properties. The all-to-all solution is clean and the expansion coefficients are supplied; the square-lattice simulations, though small, qualitatively support the split. That is a solid step beyond the two-state predecessors.\n\nThe soft spot is the central classification. Equation (2) has a > 0, so on the ordered side (ε < ε1c) the linear coefficient a(ε−ε1c) is negative. That makes m = 0 a linearly stable fixed point there. The ordered branch must then coexist with a stable disordered state — bistability, or a spinodal limit — not a continuous critical point at which the disordered state becomes unstable. The authors reserve the word \"spinodal\" for the ε2c transition, but the same logic applies to ε1c. Consequently, the claimed β = 1 is not a true critical exponent, and the derived δ = 2 scaling for entropy production, power fluctuations, and efficiency inherits that mischaracterization. This is not a minor quibble; it undercuts the \"two different scalings\" narrative.\n\nThat said, the phenomenon of two transition points appears robust. Whether both transitions are first-order-like or one is a spinodal, the asymmetry between the phases is real and worth understanding. The square-lattice data lack error bars and finite-size scaling, which is a minor issue for a Letter but worth noting.\n\nThis paper deserves a serious referee — the idea is fresh and the analytic machinery is reproducible — but it needs major revision before publication. The authors should re-examine the nature of ε1c and either reframe the claims or find a diagnostic that distinguishes a genuine critical point from a metastability limit in this nonequilibrium setting. I'd send it to review with that explicit request.","headline":"The splitting of the transition into two distinct points is a real and interesting effect, but the paper's claim that one of them is a continuous critical point doesn't survive a look at the paper's own linear stability analysis.","tokens_in":17206,"tokens_out":4095,"would_cite":false,"duration_ms":41206,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.70.Ln","05.50.+q","05.70.Fh"],"model":"deepseek-v4-flash","headline":"In a three-state Ising model driven by two opposite thermal baths, the order-disorder transition splits into two distinct transition points.","keywords":["nonequilibrium phase transition","driven Ising model","two thermal baths","entropy production","power fluctuations","critical exponent","collective heat engine"],"falsifier":"Solve the exact all-to-all steady-state equations at F=2, β1=2, β2=1, and plot the stable order parameter m against ε just below ε1c; if the data follow (ε1c−ε)^1/2 rather than (ε1c−ε), or if a stable m=0 solution coexists with the m>0 branch over a finite interval, the transition at ε1c is not the continuous β=1 transition claimed. A Binder cumulant crossing would also distinguish the two cases.","tokens_in":16168,"feed_emoji":"🔥","tokens_out":7129,"duration_ms":60745,"temperature":0.7,"pith_summary":"The paper introduces a three-state Ising model in which each spin is coupled to a hot and a cold bath that drive transitions in opposite directions. It claims that, unlike equilibrium order-disorder transitions, the critical point splits: the phase dominated by −1 spins orders continuously at ε1c, while the phase dominated by +1 spins orders discontinuously at a different point, ε2c. The continuous transition has an unusual mean-field exponent β=1, and the entropy production, power, and power fluctuations all scale with the same exponent δ=2. The two ordered phases also differ thermodynamically, with the less dissipative phase supporting heat-engine operation near maximum power and efficiency. If correct, the work shows that coupling to multiple driving baths can create new universality classes of nonequilibrium phase transitions.","feed_headline":"Opposing baths split one phase transition into two","feed_subtitle":"The two ordered phases have different critical exponents, dissipation, and heat-engine performance.","key_machinery":"The load-bearing object is the Landau-type expansion of the order-parameter dynamics, dm/dt ≈ a(ε−ε1c)m + b m² + c m³, obtained from the exact master equations for the densities n−, n0, n+ in the all-to-all limit. The coefficient b ∝ sinh[F(β1−β2)/4] vanishes only when the two baths are identical; its presence makes the transition at ε1c continuous with β=1 and shifts the critical point for the other ordered phase, producing ε2c≠ε1c. The same expansion is used to derive the scaling of entropy production, power fluctuations, and efficiency near ε1c.","core_discovery":"The central discovery is that the interplay of two thermal baths with opposite nonconservative driving forces splits the order-disorder transition of a minimal Ising model into two separate transition points, one for each ordered phase. In the all-to-all mean-field limit, which is solved exactly, the order parameter m obeys dm/dt ≈ a(ε−ε1c)m + b m² + c m³ near the first transition. Because b is nonzero whenever the baths are asymmetric (F≠0 and β1≠β2), the m² term breaks the up-down symmetry and forces m ∼ (ε1c−ε), i.e., β=1 instead of the standard β=1/2. The transition at the second point, ε2c, is discontinuous and has a spinodal region where final states depend on initial conditions. Near ε1c, entropy production, power variance, and efficiency all exhibit the same scaling exponent δ=2β=2. The same qualitative behavior is found numerically for nearest-neighbor couplings on a square lattice.","pith_inferences":["If ε1c turns out to be a spinodal rather than a true critical point, the claimed exponent β=1 would characterize the limit of metastability rather than a critical point, which would change the universality statement; this distinction is testable with a Binder-cumulant analysis.","The mechanism of a nonzero m² term from opposing baths is generic, so similar split transitions and phase-dependent thermodynamic advantages might appear in driven chemical reaction networks or active matter models.","The convergence of optimal power and efficiency conditions to the discontinuous transition point suggests that collective heat engines could generally exploit proximity to a nonequilibrium transition for performance, without requiring engineered control protocols."],"forward_implications":["The order-disorder transition in this nonequilibrium Ising model is characterized by two distinct transition points, ε1c and ε2c, rather than a single critical point.","The continuous transition at ε1c has mean-field exponent β=1 and thermodynamic scaling exponent δ=2, both distinct from the equilibrium values β=1/2 and δ=1.","The two ordered phases have different thermodynamic signatures: phase A has lower entropy production and power fluctuations than phase B.","Heat-engine operation occurs only in phase A, whose maximum-power and maximum-efficiency operating points converge to the discontinuous transition point as the temperature asymmetry increases.","These results hold beyond the all-to-all mean-field limit, as shown by square-lattice simulations."],"supporting_citations":[{"why":"Supplies the driven two-bath master-equation framework and the thermodynamic definitions of power, heat, and entropy production.","marker":"[5]"},{"why":"Provides the phenomenological steady-state approximation used for |m|≈1 and reports hints of phase-dependent transitions in a simpler two-state model.","marker":"[7]"},{"why":"Introduces the collective heat-engine/dud transition that this work extends to two distinct ordered phases.","marker":"[15]"},{"why":"Supplies the equilibrium three-state Ising model whose continuous, discontinuous, and tricritical transitions are the baseline being modified.","marker":"[34]"},{"why":"Provides the spanning-tree method used to derive approximate steady-state densities in the ordered phases.","marker":"[35]"},{"why":"Supplies the large-deviation method used to compute power fluctuations γP.","marker":"[36]"}],"fun_headline_variants":["Competing baths split Ising criticality into two","Driven Ising model shows two distinct transitions","Asymmetric drives yield split phase boundaries","Dual baths produce dual critical points in Ising","Two ordered phases from competing thermal baths"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper relies on the assumption that the expansion dm/dt ≈ a(ε−ε1c)m + b m² + c m³ faithfully describes the full dynamics near the 'up' transition, meaning that ε1c is a genuine continuous critical point rather than a spinodal limit of metastability.","fun_headline_variants_meta":{"raw":{"variants":["Competing baths split Ising criticality into two","Driven Ising model shows two distinct transitions","Asymmetric drives yield split phase boundaries","Dual baths produce dual critical points in Ising","Two ordered phases from competing thermal baths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000415,"raw_usage":{"total_tokens":2144,"prompt_tokens":950,"completion_tokens":1194,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":1139}},"tokens_in":566,"tokens_out":1194,"duration_ms":9078,"temperature":1.0,"reasoning_tokens":1139,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:07:13.186770+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the exact all-to-all steady-state equations at F=2, β1=2, β2=1, and plot the stable order parameter m against ε just below ε1c; if the data follow (ε1c−ε)^1/2 rather than (ε1c−ε), or if a stable m=0 solution coexists with the m>0 branch over a finite interval, the transition at ε1c is not the continuous β=1 transition claimed. A Binder cumulant crossing would also distinguish the two cases.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the phenomenological steady-state approximation used for |m|≈1 and reports hints of phase-dependent transitions in a simpler two-state model."},{"cited_title":"Vroylandt, M","cited_arxiv_id":null,"evidence_quote":"Introduces the collective heat-engine/dud transition that this work extends to two distinct ordered phases."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the equilibrium three-state Ising model whose continuous, discontinuous, and tricritical transitions are the baseline being modified."},{"cited_title":"Schnakenberg, Reviews of Modern physics 48, 571 (1976)","cited_arxiv_id":null,"evidence_quote":"Provides the spanning-tree method used to derive approximate steady-state densities in the ordered phases."}],"review_version":1}