{"id":"45b24b01-ded3-4525-ba68-5ea2c56207b9","arxiv_id":"2607.08425","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"2D cubically driven Bose-Hubbard lattices with single-photon losses show 2D classical three-state Potts criticality; 1D chains with three-photon losses show 1D quantum three-state Potts criticality.","lead":"This paper simulates lattices of three-photon-driven optical resonators and finds that their ordered steady states show the same critical behavior as the three-state Potts model. If correct, it is the first demonstration of Potts universality in driven-dissipative quantum matter, extending the known Ising cases to a richer symmetry class.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Potts-class claim rests on fixed-exponent collapses and a confounded two-parameter comparison; neither uniquely tests the universality class.","rationale":"The paper is a serious numerical study, and the VMG method is benchmarked at the single-mode level, but the many-body evidence for the headline claim is thinner than the text suggests. The load-bearing concern is not that the results are wrong, but that the finite-size scaling collapses are presented as demonstrations of Potts universality while using the Potts exponents as fixed inputs. A genuine test would extract exponents from the data and compare statistically against alternatives. The 1D 'quantum Potts' identification is especially fragile because the exponents used are the same as for 2D classical Potts; the only thing distinguishing 'quantum' from 'classical' in the 1D setting is the dimension and the presence of three-photon losses, but the two parameters are changed together. The reader's listed weakest assumption (VMG many-body convergence) is a real secondary concern, and the reader's rationale does mention the fixed-exponent and confounded-comparison issues, so there is partial agreement. However, I do not think this rises to the level of rejection; the central claim is plausible and could be established with a more rigorous scaling analysis and additional control simulations. Therefore the CONDITIONAL verdict remains appropriate, and no change to the verdict is needed.","tokens_in":10564,"tokens_out":4725,"duration_ms":48648,"concrete_test":"Re-analyze the raw steady-state order-parameter data of Figs. 2 and 3 with Gc, β, and ν as free parameters, using a collapse-optimization metric (e.g., minimize the integrated squared deviation between curves after rescaling, with bootstrap resampling over the data). Compare the best-fit exponents and their confidence intervals to the Potts values (β=1/9, ν=5/6), and also test alternative fixed exponent pairs (e.g., Ising β=1/8,ν=1; mean-field β=1/2,ν=1/2). If the Potts exponents lie outside the bootstrap errors, or if an alternative pair yields a comparable or better collapse, then the universality-class identification is not established by the present data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the transition belongs to the 3-state Potts universality class. The evidence consists of finite-size scaling collapses in Figs. 2(b) and 3(b), using Eq. (9) with exponents β=1/9 and ν=5/6 fixed to the target class. This is a consistency check, not a test: the exponents are inputs, not outputs. With only four lattice sizes (2D: 5×5–8×8; 1D: L=12–32) and no reported collapse residual, a range of exponent pairs can produce acceptable visual collapses, so the data do not uniquely discriminate Potts criticality from, e.g., Ising or mean-field Z3 criticality. Moreover, in the 1D case the quantum 3-state Potts exponents coincide with the 2D classical Potts exponents via the quantum-to-classical correspondence, so the Fig. 3 collapse cannot by itself certify 'quantum' criticality. The classical-to-quantum interpretation is inferred by comparing 2D with η3=0 (Fig. 2) against 1D with η3=γ (Fig. 3), varying both dimension and dissipation channel simultaneously; no data are given for 1D with η3=0 or 2D with η3>0. Thus the mechanism 'multiphoton losses promote quantum universality' (Table I) is not isolated. Separately, the single-mode benchmark (Fig. 1) does not validate VMG convergence in the many-body regime, so size-dependent variational errors could also distort the raw curves and their collapses.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies a Bose-Hubbard lattice with three-photon parametric driving, whose Lindblad dynamics has an exact Z3 symmetry. In the strong-driving ordered phase, the authors derive an effective energy functional (Eq. (6)) equal to the three-state clock/Potts ferromagnet. They simulate the steady state with the Variational Multi-Gaussian (VMG) phase-space ansatz, benchmark it against exact single-mode dynamics, and perform finite-size scaling (FSS) on four system sizes per dimension. For 2D square lattices with single-photon losses only, they report a collapse with the 2D classical 3-state Potts exponents β=1/9, ν=5/6 and fitted Gc=0.661. For 1D chains with additional three-photon losses, they report a collapse with the same exponents (which also characterize the 1D quantum Potts transition via quantum-to-classical mapping) and fitted Gc=1.141. They conclude that driven-dissipative bosonic lattices realize classical and quantum Potts criticality and propose a general correspondence: n-photon drive fixes the Z_n symmetry, while n-photon losses promote classical to quantum criticality.","tokens_in":10957,"tokens_out":5670,"duration_ms":56876,"significance":"If established, this would be a substantial extension of the emergent-equilibrium-universality program from the Z2 Ising paradigm to the three-state Potts universality class. The VMG method is a nontrivial methodological step: stochastic phase-space approaches fail when three-photon losses are present, and the exact single-mode benchmark (Fig. 1) is a genuine validation of the ansatz in a strongly non-Gaussian regime. The proposed classical-to-quantum knob (presence of the n-photon loss channel) is elegant and experimentally relevant. However, as presented, the universality-class identification is a consistency check with fixed exponents, not a model-independent extraction; the 1D claim is not distinguishable from the 2D classical class by the collapse data; and no control simulation separates the change in dimension from the change in dissipation channel. The paper therefore currently establishes a plausible and interesting scenario rather than a definitive classification.","major_comments":[{"comment":"The FSS collapses use β=1/9 and ν=5/6 as fixed inputs, not as outputs. With only four lattice sizes (2D: 5×5–8×8; 1D: L=12–32), no collapse residuals, no error bars on ⟨Z3⟩ or Gc, and no correction-to-scaling analysis, Eq. (9) cannot discriminate the claimed Potts class from nearby universality classes or from a weakly first-order transition with finite-size rounding. Please extract β/ν and 1/ν from the data (or at least report a collapse-quality metric over a grid of exponents), include bootstrap confidence intervals on Gc, and test stability with correction-to-scaling terms. Without this, Figs. 2(b) and 3(b) demonstrate compatibility, not identification.","section":"Universality class; Eq. (9), Figs. 2 and 3"},{"comment":"The 1D quantum 3-state Potts exponents coincide with the 2D classical 3-state Potts exponents through the standard quantum-to-classical correspondence. Therefore the collapse in Fig. 3(b) cannot certify quantum criticality; it is equally consistent with 2D classical Potts criticality. The comparison also changes dimension and dissipation channel simultaneously (2D with η3=0 versus 1D with η3=γ), so the mechanism 'three-photon losses promote quantum universality' summarized in Table I is not isolated. Please provide control data (e.g., 1D with η3=0 and 2D with η3>0) or use an intrinsically quantum diagnostic (entanglement structure, central charge, or dynamical exponent) that distinguishes the two interpretations.","section":"Universality class; Fig. 3(b), Table I"},{"comment":"The VMG accuracy is benchmarked only against exact single-mode dynamics. No many-body convergence check in NG is shown for the 5×5–8×8 or L=12–32 lattices used in the FSS. Since NG=24 and NG=30 are used in different settings, it is not evident that variational error is controlled at these sizes; size-dependent variational bias could be mistaken for genuine finite-size scaling. Please include convergence of ⟨Z3⟩ with NG for representative system sizes, or error bars derived from the variational uncertainty, before interpreting the collapse.","section":"Numerical method; Fig. 1, Fig. 2(a), Fig. 3(a)"},{"comment":"Equation (6) is derived in the strong-driving product-state coherent-state limit and is a valid effective classical Potts Hamiltonian there. Its use at the critical point rests on the assertion that the exact Z3 symmetry of the Lindbladian determines universality. That assertion is plausible but not demonstrated; Z3 symmetry alone does not fix the universality class, as different three-state models can have different critical behavior. The finite-size data are the only evidence and, as noted above, are currently insufficient. Please clarify what part of the Potts-class claim is derived and what part is numerical hypothesis-testing.","section":"Model and Potts structure / End Matter; Eq. (6)"}],"minor_comments":[{"comment":"The caption labels panels (b) and (c), but the text refers to 'Figure 1(a) shows a snapshot...' with the VMG Wigner function; the panel lettering should be made consistent.","section":"Fig. 1"},{"comment":"The choices NG=24 (2D) and NG=30 (1D) are stated but not justified. Please state how the number of Gaussian components was chosen and how sensitive the order-parameter curves are to this choice.","section":"Figs. 2 and 3"},{"comment":"The steady-state averaging procedure is described as 'a time window of 5γ^{-1} (500 time steps)' but the time step is not given. Please report the integration time step and the thermalization/relaxation criterion.","section":"Fig. 3(a)"},{"comment":"The column headers are terse. Define 'd' and clarify entries such as 'absent'/'present' (also, 'n-loss' meaning n-photon loss). A short caption explaining the correspondence would help.","section":"Table I"},{"comment":"The footnote claims the collapse quality is comparable to Ref. [14] despite smaller sizes, but no quantitative metric is given. Please provide the measure used for this comparison.","section":"Footnote 42"},{"comment":"There are accented-character artifacts in the author affiliation (e.g., 'Universit´ e Paris Cit´ e'). These are cosmetic but should be cleaned before final submission.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the VMG approach is interesting. My main concerns are evidentiary, not ethical or methodological-in-principle: the fixed-exponent collapses and the lack of control simulations leave the central claim underdetermined. If the authors can extract exponents from the data, supply control runs separating dimension from dissipation channel, and quantify variational convergence, a revised version could meet the bar. No concerns about novelty or citation behavior."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a genuine new direction — the first concrete attempt I know to put a driven-dissipative bosonic lattice into the Z3/Potts universality class — but the evidence is thinner than the conclusions. The 2D classical-Potts identification is plausible; the 1D quantum-Potts claim is not actually isolated by their data. I'd referee it, but I'd send it back for serious work.\n\nWhat's good: the effective Hamiltonian derivation (Eq. 6 and End Matter) is clean. In the strong-driving limit, the hopping term maps exactly onto the 3-state Potts ferromagnet with coupling 2J|α|^2/z. That's a nice argument and it's not circular. The Weyl-symbol calculation for the Z3 parity operator is formal and useful. The VMG method is interesting — the single-mode benchmark in Fig. 1 shows real accuracy with Wigner negativity at η3t=10. The idea of using multiphoton losses to promote classical to quantum criticality is attractive and worth testing.\n\nWhere it's soft — and this is the load-bearing issue — the FSS collapses in Figs. 2 and 3 use the target exponents β=1/9, ν=5/6 as fixed inputs. That's a consistency check, not a measurement. With four small lattice sizes (5–8 in 2D, 12–32 in 1D), no error bars, no collapse residual, and no correction-to-scaling, a range of exponent pairs will give acceptable visual collapses. I'm not saying the Potts class is wrong; I'm saying the data don't uniquely pick it out.\n\nThe 1D quantum claim is confounded. They compare 2D with η3=0 (classical Potts) against 1D with η3=γ (quantum Potts). Both dimension and loss channel change at once. The quantum exponents for the 1D chain are the same numbers as the 2D classical exponents via quantum-to-classical correspondence, so the collapse in Fig. 3 cannot distinguish classical from quantum criticality. To make the 'multiphoton losses promote quantum universality' mechanism stick, they'd need 1D with η3=0 and 2D with η3>0, and ideally error-quantified exponent extraction.\n\nAlso, the VMG benchmark is single-mode only. The many-body variational ansatz is the authors' own; without a convergence check at these sizes and NG values, method error could distort the raw curves. That is a premise issue, not a claim about the physics.\n\nWho should read this: people working on driven-dissipative lattice criticality and variational phase-space methods. It's a useful opening, but treat the central claims as hypotheses, not demonstrations.\n\nRecommendation: send it to peer review — it's a significant claim with a real method — but with referees who will push for controlled runs, error bars, and a decomposition of the dimension/loss-channel confound.","headline":"A genuine first step toward Z3/Potts criticality in driven-dissipative lattices, but the evidence is a consistency check with assumed exponents and the 1D quantum claim is confounded — needs controlled runs before the claims hold.","tokens_in":11419,"tokens_out":2244,"would_cite":false,"duration_ms":21286,"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":"Driven-dissipative bosonic lattices host emergent three-state Potts criticality, classical in 2D and quantum in 1D.","keywords":["driven-dissipative bosonic lattices","Potts universality class","Z3 symmetry breaking","finite-size scaling","multiphoton dissipation","open quantum systems","variational multi-Gaussian ansatz","nonequilibrium steady state"],"falsifier":"Compute the same 2D transition with an independent method—tensor network, exact diagonalization, or stochastic wavefunction—for lattice sizes 4x4 through 8x8, and check whether the ⟨Z3⟩ curves still collapse with beta=1/9 and nu=5/6 at the same critical drive. If the collapse requires different exponents or breaks down as sizes grow, the Potts-universality claim is false.","tokens_in":10464,"feed_emoji":"⚛️","tokens_out":5501,"duration_ms":47208,"temperature":0.7,"pith_summary":"This paper asks whether a genuinely nonequilibrium system—a lattice of photon-lossy nonlinear resonators with three-photon driving—can display the same critical behavior as the three-state Potts model, a three-state cousin of the Ising model from equilibrium statistical mechanics. The authors argue that the steady state spontaneously breaks a Z3 symmetry and that the transition falls in the universality class of the 2D classical three-state Potts model when only single-photon losses are present, and in the 1D quantum three-state Potts class when three-photon losses are included. If right, this extends the known Z2/Ising correspondence to higher discrete symmetries and establishes multiphoton losses as the switch between classical and quantum criticality.","feed_headline":"3-state Potts criticality emerges in driven bosonic lattices","feed_subtitle":"Driven-dissipative photon arrays reproduce an equilibrium transition; multiphoton loss flips the universality class.","key_machinery":"The central objects are the Z3 parity operator (a three-valued discrete rotation of the photon-number phase) as order parameter, and the Variational Multi-Gaussian (VMG) ansatz, which represents the many-body Wigner function as a sum of rotating Gaussian components. The linking identity is the effective lattice energy: in the ordered phase the hopping term alone distinguishes configurations and reduces to a ferromagnetic cosine coupling that is exactly the three-state Potts interaction, so the universality class is fixed by the Liouvillian's exact Z3 symmetry rather than by microscopic details.","core_discovery":"For a lattice of Kerr-nonlinear cavities pumped by three-photon parametric driving, the nonequilibrium steady state breaks the Z3 symmetry of the Liouvillian. The authors show that the effective interaction between sites, derived in the strong-driving limit, is exactly the three-state Potts Hamiltonian, and they provide finite-size scaling collapses of the Z3 parity order parameter with the Potts exponents beta=1/9 and nu=5/6. The claim is that in two dimensions with single-photon losses the transition belongs to the 2D classical three-state Potts universality class, while in one dimension with three-photon losses it belongs to the 1D quantum three-state Potts class, with the baths' multipho","pith_inferences":["The paper's quantum-classical distinction rests on dimension and loss channels; a natural extension is to test whether four-photon drives realize Ashkin-Teller criticality, as the authors conjecture.","The strong-driving derivation suggests the Potts identification may hold even outside the semiclassical limit because the exact Z3 symmetry of the Lindbladian, not microscopic parameters, controls universality—if so, the exponents should be robust to moderate changes in detuning and Kerr nonlinearity.","If the 1D quantum Potts universality is realized, a lattice with boundaries might exhibit parafermionic edge behavior, linking this steady-state transition to topological proposals—something the paper does not address.","A direct experimental signature: measuring ⟨Z3⟩ versus drive amplitude in a 3-photon-driven resonator array should show the predicted critical drift and shift exponents."],"forward_implications":["The Z3 parity drop sharpens with lattice size and collapses onto one curve using beta=1/9 and nu=5/6, so the 2D transition is continuous and Potts-like.","Turning on three-photon losses in 1D yields the same exponents, signaling the 1D quantum Potts class, not the classical one.","The pattern suggests a general law: n-photon drive sets the Zn symmetry and universality class, while n-photon losses promote classical to quantum criticality.","Circuit-QED and photonic platforms that already realize tri-squeezed states become candidate experimental settings for these transitions."],"fun_headline_variants":["Potts criticality from driving: photon lattices mirror equilibrium physics","Driven photonic lattices show classical and quantum Potts criticality","Multiphoton loss tips photon lattice into quantum Potts class","Z3 symmetry breaking in driven Bose-Hubbard mimics Potts model","Nonequilibrium photon arrays hit 2D and 1D Potts universal behavior"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The variational multi-Gaussian wavefunction with only 24–30 Gaussian components is accurate enough that the finite-size scaling collapses are real Potts criticality rather than artifacts of the ansatz.","fun_headline_variants_meta":{"raw":{"variants":["Potts criticality from driving: photon lattices mirror equilibrium physics","Driven photonic lattices show classical and quantum Potts criticality","Multiphoton loss tips photon lattice into quantum Potts class","Z3 symmetry breaking in driven Bose-Hubbard mimics Potts model","Nonequilibrium photon arrays hit 2D and 1D Potts universal behavior"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000514,"raw_usage":{"total_tokens":2345,"prompt_tokens":771,"completion_tokens":1574,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":1477}},"tokens_in":515,"tokens_out":1574,"duration_ms":10799,"temperature":1.0,"reasoning_tokens":1477,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T07:50:10.585744+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same 2D transition with an independent method—tensor network, exact diagonalization, or stochastic wavefunction—for lattice sizes 4x4 through 8x8, and check whether the ⟨Z3⟩ curves still collapse with beta=1/9 and nu=5/6 at the same critical drive. If the collapse requires different exponents or breaks down as sizes grow, the Potts-universality claim is false.","supporting_citations":[],"review_version":3}