{"id":"59efbee3-6645-46fa-8648-aa66b0ccdc62","arxiv_id":"2608.12625","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"In a two-atom Dicke model, atom-atom interaction tunes the superradiance threshold, and at the degeneracy point J=2δ the threshold vanishes and the photon-number scaling becomes quadratic.","lead":"A minimal model of two atoms inside a cavity shows that a tunable atom-atom interaction can move the point where light-matter coupling triggers a superradiant transition, and at a special interaction strength the transition threshold disappears entirely. The result gives experimentalists a new control knob for cavity quantum electrodynamics and changes the expected scaling of photon emission near the transition.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dissipative phase boundary for J/δ>2 rests on a sign error in the linearized spin equations; Eq. (S74) fails the γ↓→0 limit.","rationale":"I read the paper in good faith and found the closed-system analysis and the zero-threshold result at J=2δ to be internally consistent. The reader's weakest assumption, the classical oscillator limit ω≪δ, is standard for this problem and explicitly stated; it is not the most load-bearing weakness. The load-bearing concern is a concrete sign error in the dissipative stability analysis for the J>2δ side. The linearized equations around the normal state |t−⟩ are only valid for Δ>0; for Δ<0 the normal state is |t0⟩, reversing the sign of the 2Gx term. The paper carries the same characteristic equation into the Δ<0 regime and selects a finite-frequency root, producing a phase boundary formula that fails the physically required γ↓→0 limit: it gives g_c=0 while the photon-loss-only model has a nonzero threshold. The corrected reduced-model analysis yields a symmetric formula g_c² = (ω²+κ²/4)(Δ²+γ↓²/4)/(8ω|Δ|), which reduces properly. This error does not invalidate the central zero-threshold claim at J=2δ, since that result is obtained from the closed set of equations (9)–(11) and is supported by numerical agreement in Fig. 4(c). It does, however, invalidate the claimed analytic dissipative phase boundary for J/δ>2 and the agreement shown in Fig. 4(a) in that region. The manuscript should be accepted only conditionally, after correcting the reduced-model stability analysis and re-plotting Fig. 4(a) with the correct boundary.","tokens_in":127,"tokens_out":52315,"duration_ms":964871,"concrete_test":"Test the γ↓→0 limit: set γ↓=0 in Eq. (S74) and observe that it gives g_c=0, contradicting the pure-photon-loss boundary of SM Eq. (S37), which is nonzero for J>2δ. Then perform exact-diagonalization of the full Liouvillian (ω=1, δ=50, κ=γ↓=0.1) at a representative point J/δ=3 and locate the steady-state photon-number onset in g̃. The corrected reduced-model boundary is g̃_c ≈ sqrt[(ω²+κ²/4)(Δ²+γ↓²/4)/(8ω|Δ|)/ωδ] ≈ 0.33, whereas Eq. (S74) predicts g̃_c ≈ 8.8. If the numerical onset appears near g̃≈0.33, the published analytic boundary for J>2δ is wrong.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section III.B of the SM linearizes the reduced master equation around the normal state |t−⟩, writing ṡ_y = 2Gx + Δs_x − Rs_y (SM Eq. S65). This is the correct normal state only for Δ>0 (J<2δ). For Δ<0 (J>2δ), the normal spin state is |t0⟩, which has τ_z = +1, so the cavity-spin coupling in the τ_y equation changes sign: ṡ_y = −2Gx + Δs_x − Rs_y. Using the same characteristic equation (SM Eq. S67) with −4ΔωG² for both signs is therefore invalid for Δ<0. The resulting boundary, main-text Eq. (S74), does not reduce to the pure-photon-loss boundary when γ↓→0: Eq. (S74) contains a factor κγ↓/(κ+γ↓)² and gives g_c→0 as γ↓→0, whereas for γ↓=0 the problem is the photon-loss-only model with a nonzero threshold (SM Eq. S37). With the correct sign choice, the characteristic equation becomes [(λ+K)²+ω²][(λ+R)²+Δ²] − 4|Δ|ωG² = 0, whose zero-frequency marginal mode gives g_c² = (ω²+κ²/4)(Δ²+γ↓²/4)/(8ω|Δ|), symmetric in the sign of Δ and matching Eq. (S37) for γ↓→0. The paper's asymmetric formula, which scales as |Δ|³ for large |Δ| and vanishes when γ↓→0, is therefore incorrect. This affects the analytic phase boundary in Fig. 4(a) for J/δ>2, although the zero-threshold result at J=2δ, derived separately from the closed equations (9)–(11), remains valid.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a two-atom Dicke model with flip-flop interatomic interaction in the classical oscillator limit. It derives analytic expressions for the closed-system superradiance boundary, finding that the critical coupling vanishes at the interaction-induced degeneracy point J=2δ and that the photon-number onset becomes quadratic there. The paper extends the analysis to open systems with photon loss and local spin relaxation, deriving phase boundaries, showing that the zero threshold survives, and studying metastable dynamics and population distributions.","tokens_in":18906,"tokens_out":31779,"duration_ms":299320,"significance":"If correct, the results establish atom-atom interaction as a qualitative tuning knob for superradiant criticality, including the complete suppression of the threshold and a change of the photon-number onset exponent. The strength of the paper is its analytic tractability: the phase boundaries and exponents are derived from the Hamiltonian without fitted parameters and are confirmed by exact diagonalization. The open-system analysis, including the metastable manifold and the reduced two-level model, is well executed. The paper should be of interest to the quantum-optics and circuit-QED communities.","major_comments":[{"comment":"The Δ<0 branch of Eq. (9) vanishes in the limit γ↓→0, whereas the pure-photon-loss phase boundary (SM Eq. S37) is finite for J/δ>2. This discontinuity arises because the reduced model is derived by projecting out the |t+⟩ and singlet states, which is justified only for finite spin relaxation; the γ↓→0 limit lies outside the reduced model's validity. The authors should state this validity condition and, ideally, verify the Δ<0 branch against full-Liouvillian numerics for finite γ↓.","section":"Main text Eq. (9) and SM Sec. III.B, Eq. (S74)"},{"comment":"The claim that the transition at J=2δ belongs to a 'different universality class' is based solely on the photon-number exponent γ=2 measured at a single fine-tuned point in a two-atom model. Since no other critical exponents are computed and the limit is the classical-oscillator limit, the term 'universality class' is stronger than the evidence supports; please replace it with 'critical exponent' or provide additional scaling information.","section":"Main text near Eq. (3) and Conclusion"}],"minor_comments":[{"comment":"The definition of the critical exponent uses n(g̃)∝(g̃−g̃_c)^γ; near a zero-threshold point (g̃_c=0) the standard scaling form with a multiplicative factor is not given, which could confuse the γ=2 result. Consider stating the scaling more explicitly.","section":"Main text Eq. (3)"},{"comment":"The labels NP1 and NP2 are explained in the text but do not appear in the figure; adding the labels directly to the figure would improve readability.","section":"Fig. 1(b)"},{"comment":"The rate-equation derivation of the triplet populations assumes the cavity is in the vacuum state; this is a valid leading-order choice but the reader should be reminded that it neglects the classical displacement of the cavity in the superradiant phase.","section":"SM Sec. II.C"},{"comment":"After introducing the flip-flop interaction, the singlet state is mentioned as a dark state and then ignored. In the presence of spin relaxation the singlet is no longer dark; a short comment on its role in the reduced model would help the reader.","section":"Main text around Eq. (6)"}],"recommendation":"minor_revision","confidential_remarks":"The reader's take and the skeptic's note differ on whether the Δ<0 branch of Eq. (9) is erroneous. On reading the paper, I find that the sign used in the linearization around |t−⟩ is consistent with the reduced model's steady state for γ↓>0, so the skeptic's sign-error accusation does not land directly. The genuine residual concern is the singular γ↓→0 limit, which I have raised as a major comment. If the authors can numerically confirm the Δ<0 boundary for representative finite γ↓, I would be satisfied."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The closed-system part of this paper is a clean, legitimate extension of the Rabi/Dicke criticality story, and the zero-threshold point at J=2δ with photon number scaling n∝g² is a genuinely nice result. The interaction-tunable degeneracy that erases the threshold is real physics, and the analytic Landau expansions match the exact diagonalization points. If the paper were only about the closed system, I would be confident in it.\n\nMy concern is with the dissipative analysis. The pure-photon-loss boundary (SM Eq. S37) is derived correctly and reduces properly to the closed-system result. But the boundary for J>2δ in the presence of local spin relaxation (SM Eq. S74) is not trustworthy. The linearized equations in SM Section III.B are written around the normal state |t−,0⟩, which is the correct normal state only for Δ>0 (J<2δ). For Δ<0, the closed-system normal state is |t0⟩, and with spin relaxation the actual steady state is a mixture of |t−⟩ and |t0⟩ whose composition depends on γ↓/|Δ|. The paper's own rate equations in S78–S79 make that clear. So deriving a boundary by linearizing around pure |t−⟩ for all Δ is wrong, and Eq. (S74) indeed fails the γ↓→0 limit: it vanishes, whereas the photon-loss-only problem has a finite nonzero threshold. The stress-test note is right that something is off; I would not defend that expression as written.\n\nThat said, the fix is not just flipping a sign. The correct procedure for Δ<0 and finite γ↓ is to linearize around the mixed steady state, not around a pure spin state. The resulting boundary will interpolate between the two limits and may not be as simple as the symmetric formula in the stress-test. The important point is that the analytic dashed curve in Fig. 4(a) for J/δ>2 is not reliable as presented. The claim that the formula covers that regime needs to be either derived properly or removed.\n\nThe zero-threshold result at J=2δ (Eq. 11) is derived from closed equations and survives. The Δ>0 boundary is fine, and the metastability analysis in Fig. 3 is a nice addition. The 'different universality class' language is a bit strong for a two-atom toy model with a single exponent at one point, but the authors do call it a toy model, so I would treat that as a wording issue rather than a fatal one. The citation pattern looks fine; the one self-citation is not load-bearing.\n\nRecommendation: send to peer review, but expect a major revision. The closed-system result is worth publishing, and the dissipative section can be fixed if the authors redo the Δ<0 boundary properly. A referee who checks S74 against the γ↓→0 limit will catch this, so it needs to be addressed before publication.","headline":"The closed-system analysis is solid and the J=2δ zero-threshold result is genuinely new; but the dissipative phase boundary for J>2δ with spin relaxation rests on a flawed linearization and needs to be rederived.","tokens_in":19457,"tokens_out":11750,"would_cite":true,"duration_ms":102948,"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":"This paper shows that a flip-flop interaction between two atoms in a cavity can completely remove the superradiance threshold and change its critical scaling.","keywords":["Dicke model","superradiant phase transition","atom-atom interaction","zero-threshold superradiance","critical exponent","dissipative phase transition","quantum Rabi model","cavity QED"],"falsifier":"A direct experimental check: fix the two-atom interaction at $J=2\\delta$ in a simulator with $\\omega\\ll\\delta$ (for example, two superconducting qubits coupled to a resonator) and measure the mean photon number as $g$ is swept from zero upward. The paper predicts $n\\propto\\tilde g^2$ for every nonzero $g$; a linear onset $n\\propto(\\tilde g-\\tilde g_c)$ with a finite $\\tilde g_c>0$, or an absence of photons below a threshold, would contradict the central claim. A numerical check beyond the $\\omega\\ll\\delta$ limit would also settle it: if exact diagonalization at comparable $\\omega$ and $\\delta$ shows a finite threshold at $J=2\\delta$, the zero-threshold and $\\gamma=2$ results are artifacts of the classical-oscillator limit.","tokens_in":18334,"feed_emoji":"⚛️","tokens_out":10156,"duration_ms":97658,"temperature":0.7,"pith_summary":"This paper asks what a direct interaction between emitters does to superradiance and answers it with the simplest model that still allows exact analytics: two interacting atoms coupled to a single cavity mode in the classical oscillator limit. It shows that the superradiant threshold $\\tilde g_c$ depends non-monotonically on the flip-flop interaction strength $J$ and vanishes exactly at $J=2\\delta$, where the two lowest spin states become degenerate. At that point, any nonzero atom-photon coupling produces photons, and the photon number grows as $\\tilde g^2$ rather than linearly, placing the transition in a different universality class. The zero threshold and the quadratic onset persist under both photon loss and local spin relaxation, and near the degeneracy the dissipative dynamics shows a metastable two-stage relaxation. The result identifies atom-atom interaction, not just light-matter coupling, as a control knob for superradiant criticality in a system realizable with two qubits in circuit QED or a trapped ion.","feed_headline":"A tuned atom-atom interaction erases the superradiance threshold","feed_subtitle":"At one interaction strength, any nonzero atom-photon coupling produces photons, and the photon number rises quadratically.","key_machinery":"The central object is the triplet manifold $|t_-\\rangle$, $|t_0\\rangle$, $|t_+\\rangle$ of the two-atom spin Hamiltonian, with the singlet inert because $\\Sigma_x|s\\rangle=0$. In the classical oscillator limit $\\omega\\ll\\delta$, the photon operators are replaced by a c-number $\\alpha$, producing an effective transverse field $\\lambda=2g\\alpha$ on the spins; the ground state follows from the mean-field energy functional $E(\\alpha)=\\omega\\alpha^2+\\epsilon_0(\\lambda)$. The argument turns on the degeneracy between $|t_-\\rangle$ and $|t_0\\rangle$ at $J=2\\delta$: nondegenerate perturbation theory fails there, and projecting onto the degenerate subspace gives $\\epsilon_0(\\lambda)=-2\\delta-\\sqrt{2}|\\lambda|$, a linear cusp that makes $\\alpha=0$ unstable for any $g>0$. That cusp is what produces the quadratic photon-number onset. For the open system, the load-bearing reduction is the spin-projected effective Liouvillian, which maps the dissipative problem onto a damped Rabi model with detuning $\\Delta=2\\delta-J$, so all threshold calculations reduce to linear stability of that effective model.","core_discovery":"Within the triplet sector of two coupled spins, the cavity couples only through $\\Sigma_x$, and the spin ground state changes from $|t_-\\rangle$ to $|t_0\\rangle$ at $J=2\\delta$. Perturbation theory in the effective transverse field $\\lambda=2g\\alpha$ yields the closed-form critical coupling $\\tilde g_c^2=(2-J/\\delta)/8$ for $J<2\\delta$ and $\\tilde g_c^2=(J/\\delta-4\\delta/J)/16$ for $J>2\\delta$, with $\\tilde g_c=0$ at the degeneracy. At degeneracy the mean-field energy acquires a linear cusp in $\\alpha$, so the origin is unstable for every $g>0$; minimizing gives $n\\propto\\tilde g^2$, i.e. the exponent $\\gamma=2$, whereas away from degeneracy $\\gamma=1$. Cavity decay only rescales the threshold by the factor $1+(\\kappa/2\\omega)^2$, and with local spin relaxation the system reduces to an effective Rabi model with detuning $\\Delta=2\\delta-J$, whose threshold still vanishes at $\\Delta=0$ while the photon number at degeneracy obeys $n_{\\rm ss}=2(1+\\gamma_\\downarrow/\\kappa)g^2/((\\kappa+\\gamma_\\downarrow)^2/4+\\omega^2)$. The normal phase on the $J>2\\delta$ side is nearly maximally entangled, and crossing from it into the superradiant phase decreases rather than increases atomic excitation.","pith_inferences":["Extension not in the paper: the mechanism looks generic — any interaction that degenerates two low-energy atomic states connected by the cavity-coupling operator should produce a zero-threshold soft mode, so similar threshold cusps should appear in multi-atom and multimode Dicke generalizations.","Extension not in the paper: the predicted crossover from $\\gamma=1$ to $\\gamma=2$ as $J/\\delta$ approaches 2 from either side is directly measurable by fitting $n(\\tilde g)$ at several fixed interaction strengths; this would test the universality-class change without needing exact degeneracy.","Extension not in the paper: the effective Rabi detuning $\\Delta=2\\delta-J$ means the interaction can be used as a knob to reach effective strong coupling at small bare $g$, with practical implications for ultrastrong-coupling experiments."],"forward_implications":["At $J=2\\delta$, the cavity is populated for every nonzero $g$; a measurement of $n(\\tilde g)$ should show $n\\propto\\tilde g^2$ ($\\gamma=2$), and this holds under cavity photon loss and under local spin relaxation.","The phase boundary $\\tilde g_c(J/\\delta)$ is non-monotonic and has a cusp at zero at $J/\\delta=2$; detuning the interaction either side restores the conventional linear onset $\\gamma=1$, so the degeneracy is the control parameter that changes the universality class.","In the presence of photon loss alone, the steady state near degeneracy has all three triplet states almost equally populated, and relaxation from $|t_-,0\\rangle$ proceeds through a metastable plateau on the cavity-loss time scale before the slow Liouvillian relaxation sets in.","With local spin relaxation, the two-atom problem becomes an effective damped Rabi model with interaction-tunable detuning $\\Delta=2\\delta-J$; at $\\Delta=0$ the threshold remains zero and the $|t_0\\rangle$ population rises toward $1/2$, so dissipation combined with interaction increases atom-atom entanglement relative to the equal-population triplet steady state."],"supporting_citations":[{"why":"Supplies the classical-oscillator Rabi-model quantum phase transition that the two-atom model extends to interacting atoms.","marker":"[5]"},{"why":"Supplies the mean-field construction for finite-atom Dicke superradiance used to derive the phase boundary.","marker":"[4]"},{"why":"Supplies the open quantum Rabi dissipative phase transition and Liouvillian-spectrum framework used for photon loss and spin relaxation.","marker":"[8]"},{"why":"Contains the derivations of the closed-system critical couplings, photon-number scaling, reduced Liouvillian, and rate equations that the main text relies on.","marker":"[44]"},{"why":"Supplies the spectral theory of Liouvillians used to define dissipative phase transitions and the Liouvillian gap.","marker":"[46]"},{"why":"Supplies the metastability framework used to interpret the slow eigenvalue cluster as a metastable manifold.","marker":"[47]"},{"why":"Supplies the concurrence used to quantify atom-atom entanglement in the two normal phases and across the phase diagram.","marker":"[45]"}],"fun_headline_variants":["Atom-atom interaction erases superradiance threshold","Two-atom coupling eliminates photon threshold","Interactions abolish Dicke threshold entirely","Spin-spin coupling nullifies photon threshold","Dicke threshold vanishes with atom interactions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The sharp phase boundary and the zero threshold assume the photon frequency $\\omega$ is much smaller than the atomic transition frequency $\\delta$, so the photon field can be replaced by a classical amplitude; if this separation of scales fails, the claimed threshold suppression and quadratic onset are not guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["Atom-atom interaction erases superradiance threshold","Two-atom coupling eliminates photon threshold","Interactions abolish Dicke threshold entirely","Spin-spin coupling nullifies photon threshold","Dicke threshold vanishes with atom interactions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000332,"raw_usage":{"total_tokens":1867,"prompt_tokens":986,"completion_tokens":881,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":816}},"tokens_in":602,"tokens_out":881,"duration_ms":8193,"temperature":1.0,"reasoning_tokens":816,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:03:57.823047+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct experimental check: fix the two-atom interaction at $J=2\\delta$ in a simulator with $\\omega\\ll\\delta$ (for example, two superconducting qubits coupled to a resonator) and measure the mean photon number as $g$ is swept from zero upward. The paper predicts $n\\propto\\tilde g^2$ for every nonzero $g$; a linear onset $n\\propto(\\tilde g-\\tilde g_c)$ with a finite $\\tilde g_c>0$, or an absence of photons below a threshold, would contradict the central claim. A numerical check beyond the $\\omega\\ll\\delta$ limit would also settle it: if exact diagonalization at comparable $\\omega$ and $\\delta$ shows a finite threshold at $J=2\\delta$, the zero-threshold and $\\gamma=2$ results are artifacts of the classical-oscillator limit.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the derivations of the closed-system critical couplings, photon-number scaling, reduced Liouvillian, and rate equations that the main text relies on."}],"review_version":1}