{"id":"a82016f8-a28d-4f11-b669-ac6f05ffe582","arxiv_id":"2607.21296","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Fluctuation-induced bistability in cavity-coupled many-body systems is shown to be generic across bosons, fermions, and spins, with a new dressed-state rate-equation method capturing its dynamics.","lead":"This paper shows that 'fluctuation-induced bistability'—where a cavity-driven quantum system can settle into two very different steady states—occurs in bosonic, fermionic, and spin models coupled to a single lossy cavity mode, not just in the fermion case found earlier. It introduces a dressed-state rate-equation method for large systems, checks it against tensor-network simulations, and argues the effect is a genuine many-body phenomenon visible even in small systems.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The thermodynamic-limit 'true stabilization' rests on a three-point DSRE exponential fit, and the DSRE closure discards long-lived same-Θ coherences that exact tMPS shows are present; the scaling may be a closure artifact.","rationale":"The reader's weakest assumption was the three-point thermodynamic extrapolation. I agree that this is the immediate vulnerable step, but I locate the deeper problem one level earlier: the scaling is computed from DSRE, whose no-coherence ansatz is known from the paper's own tMPS comparison to be violated (coherences C2/Ekin persist and reduce ⟨θ²⟩). This makes the γ∝a^{-L} fit a fit to an uncontrolled closure, not to the exact dynamics. The proposed check—either extending DSRE to larger L with alternative fits, or better, including the discarded same-Θ coherence sectors—would settle whether the exponential stabilization is physical. The paper otherwise gives real support: three model classes, two independent approximate methods that agree qualitatively, and exact small-system tMPS showing non-monotonic ⟨θ²⟩ and slow dynamics near the predicted bistable region. So the correct response is to keep the CONDITIONAL verdict: the core finite-time metastable physics is credible, but 'true steady-state bistability in the thermodynamic limit' should not be accepted until the scaling is either anchored by exact/controlled calculations or softened in the claims.","tokens_in":22832,"tokens_out":13796,"duration_ms":138303,"concrete_test":"Re-derive the DSRE slow manifold without the diagonal truncation: keep the same-Θ coherence blocks that Eq. (6) discards (a generalized master equation for populations plus those coherences) and recompute γ(L) for L=112, 224, 336 at the Fig. 2 parameters (U=40J, δ=5J, Γ=3J, g=18J, N=L/2). If γ changes by more than a factor of two or the exponential-in-L trend is not preserved, the thermodynamic stabilization is not supported. As a cheaper cross-check, run tMPS at L=16 with the same parameters and compare the DSRE peak lifetime and ⟨θ²⟩; if they disagree beyond the stated stochastic uncertainty, the closure error is present at the sizes that anchor the extrapolation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the low-|θ| peak becomes truly stable in the thermodynamic limit is carried by Fig. 2(f): an exponential rate γ ∝ a^{-L} fitted to three system sizes (L=112, 224, 336). This is doubly fragile. First, the quantity being extrapolated is produced by DSRE, a projected rate equation built on the diagonal ansatz Eq. (6). That ansatz explicitly drops coherences between states with equal Θ_n, which the authors note are long-lived under the approximate strong symmetry; tMPS at L=14/16 finds substantial coherences and finite-tunneling corrections that reduce ⟨θ²⟩ relative to DSRE. If those coherences participate in the slow manifold, the DSRE rates—and hence γ(L)—need not be the true Lindblad rates. Second, even accepting DSRE, three points without error bars cannot distinguish exponential decay from a power law or stretched exponential; a small curvature at larger L would change 'stabilization' into merely long metastability. The exact tMPS checks reach only L≤16 and times ~500–1000 J/ℏ, so they do not constrain the scaling regime where the claim lives.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the dissipative dynamics of one-dimensional interacting many-body systems (spin-1, spinful fermions, and spinless bosons) globally coupled to a lossy cavity. It claims that fluctuation-induced bistability, previously found for fermionic atoms, is a robust and generic phenomenon across particle statistics, arising from a common resonance mechanism between photonic transitions and many-body energy scales. The authors use two approximate methods (mean-field plus fluctuations, MF+fluct., and a newly introduced dressed-state rate equation, DSRE) and compare with numerically exact time-dependent matrix product state (tMPS) simulations for small systems. The paper reports bistable steady-state solutions in all three models, shows dynamical signatures of bistability at finite times, and proposes an exponential system-size slowdown of the metastable relaxation, suggesting true stabilization in the thermodynamic limit.","tokens_in":23110,"tokens_out":7403,"duration_ms":70145,"significance":"If the claims hold, the paper establishes a generic beyond-mean-field phenomenon in cavity-matter systems and provides a tractable large-scale method (DSRE) for studying dissipative many-body dynamics. The main strengths are the breadth of evidence: three qualitatively different microscopic models, two independent approximate methods with a detailed derivation of DSRE in the Supplemental Material, documented convergence parameters for tMPS, and a data-availability statement. The paper is significant because it suggests that the previously reported fermionic fluctuation-induced bistability is not an artifact of a particular approximation or of fermionic statistics, and it identifies experimentally relevant finite-time signatures in small systems.","major_comments":[{"comment":"The claim that the low-|θ| peak becomes truly stable in the thermodynamic limit rests entirely on the exponential rate law γ ∝ a^{-L} extracted from a linear fit to log(ℏγ/J) at only three system sizes (L=112, 224, 336), with no error bars. With three points one cannot distinguish exponential from power-law or stretched-exponential decay, and a small curvature at larger L would convert 'true stabilization' into merely long-lived metastability. The exact tMPS results are restricted to L≤16 and times Jt/ℏ≤1000 and thus do not constrain this scaling regime. Please provide additional system sizes (e.g., L=448, 560) or an analytical argument for the exponential scaling; otherwise, the abstract and conclusion should be rephrased to claim exponentially long-lived metastability rather than steady-state bistability in the thermodynamic limit.","section":"Scaling with system size, Fig. 2(e)-(f)"},{"comment":"The DSRE ansatz Eq. (6) explicitly neglects coherences between eigenstates with equal Θ_n, which the paper itself notes are long-lived under the approximate strong symmetry (Ref. [64]). The tMPS results in the SM (Figs. 3-4) show substantial off-diagonal atomic coherences (quantified by E_kin and C_2) and finite-tunneling corrections that reduce ⟨θ²⟩ relative to DSRE. Since the rates R_{n→n′} in Eq. (8) are derived under this closure, the system-size dependence γ(L) in Fig. 2(f) may be a closure artifact rather than the true Lindblad relaxation rate. Please assess quantitatively the effect of the omitted same-Θ coherences on the relaxation rate, or at least add an explicit caveat in the main text where the exponential scaling is presented.","section":"Dressed-State Rate Equation, Eqs. (6)-(8), and SM §'Dissipative dynamics results'"},{"comment":"The 'common microscopic mechanism' is summarized by the resonance condition 2ℏgλ≈pU±ℏδ (p=1 for spins and fermions, integer p≥1 for bosons). This condition is not derived in the main text or the SM for the spin-1 and boson models. For the spin-1 model, H∥ is a sum of on-site (S^z_j)^2 terms, and it is not transparent why photon-assisted transitions couple to energy differences pU. Since this resonance condition is presented as the unifying mechanism, please provide the derivation (using Eq. (3) in the small-J limit) or explicitly state it as a heuristic conjecture. This is important for the abstract's claim of a common microscopic mechanism.","section":"Models and Lindblad equation; text near 'The fluctuation-induced bistability emerges...'"}],"minor_comments":[{"comment":"The blue lines are plotted as −|θ_MF|, but neither the text nor the caption explains why the mean-field solution is shown with a negative sign while the fluctuation-induced solutions are plotted as positive |θ_th|. Please clarify the convention.","section":"Fig. 1"},{"comment":"The notation |n, χΘ_n⟩ is used for the coherent-state displaced state, but the displacement parameter χΘ_n is not explicitly defined at first use. Define it near Eq. (4) to improve readability.","section":"Eq. (6)"},{"comment":"The expression for the energy-balance condition is written as '0 = ∂/∂t ⟨H_MF_eff⟩ ∝ ...'. The proportionality is not exact and the various prefactors are suppressed. A brief derivation or an explicit statement that this is an approximation (with the proportionality constant positive) would make the MF+fluct. method more transparent.","section":"Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the journal's readership, but the thermodynamic-limit claim is the main concern. The exponential stabilization rests on a three-point fit, and the DSRE closure problem is acknowledged by the authors. If the authors provide more system sizes or weaken the claim to exponentially long-lived metastability, the paper would be much stronger. The detailed SM and data availability are commendable. No concerns about citation ethics; the references are appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious read. The paper shows fluctuation-induced bistability appears in bosons and spin-1 systems, not just fermions, and develops a dressed-state rate equation (DSRE) that gives real-time dynamics for large systems. The DSRE is derived from a polaron transformation plus projector, no fitted parameters, and the authors check it against exact tMPS on small systems. That is real work, and the qualitative agreement across three models plus two methods makes a strong case that the effect is not an artifact of mean-field-plus-fluctuations.\n\nThe soft spots are where the stress test lands. The claim that the low-|theta| peak becomes truly stable in the thermodynamic limit is carried by Fig. 2(f) — an exponential rate gamma ~ a^{-L} fitted to three system sizes. Three points without error bars cannot distinguish exponential from power-law decay, and the quantity being extrapolated comes from DSRE, which deliberately drops coherences between degenerate Theta_n states. The tMPS data at L=14-16 show these coherences are present and reduce <theta^2> relative to DSRE. So thermodynamic-limit stabilization is plausible but not proven; the paper does say \"suggests\" in the body, but the abstract's \"universal\" framing overstates it.\n\nThat said, the core message does not depend on that extrapolation. The bistability itself is demonstrated in small systems: tMPS shows slow dynamics and non-monotonic <theta^2> vs g, matching DSRE qualitatively. That is sufficient to establish fluctuation-induced bistability as a genuine beyond-mean-field phenomenon. The MF+fluct. thermalization closure is an approximation, but DSRE bypasses it and the two agree, so it is not load-bearing. The \"generic/universal\" claim is broader than the 1D, single-mode, uniform-coupling geometry tested, but that is a standard overreach and easy to fix.\n\nWho is this for? Anyone working on cavity QED with interacting atoms or on open many-body systems. The DSRE is a potentially useful tool for large-scale dynamics. I would send it to review. The authors should be asked to soften the thermodynamic-limit language, add error bars or a fourth system size, and say explicitly that the three-point fit is suggestive. This is a solid paper with one over-extrapolated claim, not a flawed one.","headline":"Solid evidence that fluctuation-induced bistability is a genuine beyond-mean-field effect across fermions, bosons, and spins, but the thermodynamic-limit stabilization claim rests on a fragile three-point extrapolation.","tokens_in":23642,"tokens_out":2194,"would_cite":true,"duration_ms":21584,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Fluctuation-induced bistability—two coexisting steady states generated by quantum fluctuations—is, the paper claims, a generic feature of strongly interacting cavity-coupled matter, arising from photon-assisted resonances in spins, fermions","keywords":["fluctuation-induced bistability","cavity quantum electrodynamics","dissipative many-body dynamics","polaron transformation","dressed-state rate equation","metastability","Fermi-Hubbard model","Bose-Hubbard model"],"falsifier":"Run the same DSRE simulation for L=448 and L=672 and check whether ln(ħγ/J) still falls on the straight line fitted to the three smaller sizes; alternatively, push exact tMPS (or a better method) to L=28 or larger and see whether the low-|θ| peak persists in the true steady state. If the rate at larger L deviates upward from the exponential law, or if any finite-L system eventually fully transfers probability to the high-|θ| peak with a rate that does not vanish exponentially, the claimed thermodynamic-limit bistability is falsified.","tokens_in":22595,"feed_emoji":"⚛️","tokens_out":7347,"duration_ms":72809,"temperature":0.7,"pith_summary":"The authors set out to prove that fluctuation-induced bistability is not a special-case artifact but a robust consequence of strong interactions, strong light-matter coupling, and dissipation. They find the same bistable behavior in spin-1 chains, Fermi-Hubbard, and Bose-Hubbard models globally coupled to a lossy cavity, and trace it to a common mechanism: photon-assisted resonances that match photonic transitions to many-body energy gaps, making energy transfer efficient and stabilizing multiple macroscopic states. To study its dynamics they build a dressed-state rate equation based on a polaron transformation, which lets them simulate hundreds of sites over long times. The relaxation between bistable peaks slows down exponentially with system size, suggesting true steady-state bistability in the thermodynamic limit, and even small systems show precursor signatures in exact tensor-network simulations. A sympathetic reader would care because this establishes a universal phenomenon and gives concrete, experimentally accessible predictions across a wide class of hybrid quantum platforms.","feed_headline":"One mechanism drives bistability in cavity-coupled matter","feed_subtitle":"Spins, fermions, and bosons all show two coexisting steady states when quantum fluctuations and photon-assisted resonances align.","key_machinery":"The central object is the polaron transformation D = exp[(χ↠− χ*â)Θ], which displaces the cavity field by an amount proportional to the matter imbalance operator Θ. In the transformed frame, the dressed states |n, χΘ_n⟩ (matter eigenstates of Θ times a cavity coherent state) are stationary in the absence of tunneling; the rate equation for the probabilities P_n on these dressed states, derived by projecting out fast photonic degrees of freedom, carries all the slow dynamics. Its rates R_{n→n'} contain an integral over ΔE_{n',n} and ΔΘ_{n',n}, and it is the resonance condition encoded in this integral—not any particular statistics—that produces the bistable peaks and their exponential size","core_discovery":"The central claim: fluctuation-induced bistability, first reported for fermionic atoms in a cavity, appears equally in bosonic and spin systems, unified by one mechanism—resonances between photonic transitions and many-body energy scales (in the small-tunneling limit, 2ħgλ ≈ pU ± ħδ). The bistability appears only with quantum fluctuations of the light-matter coupling; mean-field theory misses it. The dressed-state rate equation, derived via a polaron transformation, treats light-matter coupling exactly and tunneling perturbatively, showing the low-imbalance peak is metastable at finite size with an escape rate decaying exponentially with system size—implying stabilization in the thermodynami","pith_inferences":["Editorial: if the exponential slowing persists, the bistability functions as a dissipative two-state memory whose switching time is tunable by system size—an engineering consequence the paper does not discuss.","Editorial: because the mechanism is purely about resonance between a photonic mode and many-body energy scales, it should transfer to other platforms with a globally coupled lossy bosonic mode and strong interactions, such as optomechanical or photonic arrays.","Editorial: a clean experimental probe would be the cavity field statistics: in the bistable regime the Glauber-Sudarshan P function, which the DSRE maps directly to p(θ), should become bimodal.","Editorial: a testable extension is to go beyond the perturbative-tunneling DSRE for moderate J, where finite-tunneling corrections (already visible in tMPS as reduced ⟨θ²⟩) may shrink or shift the bistable region; comparing the two approaches as a function of J/U would settle how much of the effect survives at larger tunneling."],"forward_implications":["Fluctuation-induced bistability is independent of particle statistics: the same phenomenon appears for interacting spins, fermions, and bosons coupled to a dissipative cavity.","In the small-tunneling limit the bistable regions sit at the resonance condition 2ħgλ ≈ pU ± ħδ, giving a quantitative guide for where to search experimentally in each model.","For finite systems the low-imbalance peak is only metastable, but its lifetime grows exponentially with system size (γ ∝ a^{-L}), so in the thermodynamic limit it becomes a true steady state.","Signatures of the bistability—non-monotonic imbalance versus coupling and strongly slowed dynamics—appear already for systems as small as L=14 in exact tMPS simulations, meaning small-scale experiments can see them.","The dressed-state rate equation opens access to parameter regimes (large L, long times) that are out of reach for exact tensor-network methods, enabling systematic studies of metastability across models."],"fun_headline_variants":["Quantum fluctuations drive universal bistability in cavity-coupled matter","Photon-matter resonances unify bistability across spin, fermion, boson models","Bistability emerges from quantum fluctuations in any cavity-matter system","Cavity-coupled systems show universal fluctuation-induced bistability","Metastable two-state behavior links spins, fermions, and bosons in cavities"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The thermodynamic-limit claim that the low-imbalance peak becomes truly stable rests on an exponential fit of escape rate versus system size taken from only three sizes (L=112, 224, 336); if that three-point extrapolation fails at larger L, the phenomenon is long-lived metastability rather than genuine steady-state bistability.","fun_headline_variants_meta":{"raw":{"variants":["Quantum fluctuations drive universal bistability in cavity-coupled matter","Photon-matter resonances unify bistability across spin, fermion, boson models","Bistability emerges from quantum fluctuations in any cavity-matter system","Cavity-coupled systems show universal fluctuation-induced bistability","Metastable two-state behavior links spins, fermions, and bosons in cavities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000182,"raw_usage":{"total_tokens":1131,"prompt_tokens":709,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":453,"completion_tokens_details":{"reasoning_tokens":326}},"tokens_in":453,"tokens_out":422,"duration_ms":4429,"temperature":1.0,"reasoning_tokens":326,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T07:52:29.254764+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same DSRE simulation for L=448 and L=672 and check whether ln(ħγ/J) still falls on the straight line fitted to the three smaller sizes; alternatively, push exact tMPS (or a better method) to L=28 or larger and see whether the low-|θ| peak persists in the true steady state. If the rate at larger L deviates upward from the exponential law, or if any finite-L system eventually fully transfers probability to the high-|θ| peak with a rate that does not vanish exponentially, the claimed thermodynamic-limit bistability is falsified.","supporting_citations":[],"review_version":1}