{"id":"e26a5af8-884b-44e8-8f4f-ede9468d06da","arxiv_id":"2411.18310","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive an exact propagator for the open dispersive Jaynes-Cummings model and a second-order perturbative solution for the open quantum Rabi model, revealing a non-thermal qubit steady state that is unique to the Rabi model.","lead":"This paper derives exact and perturbative analytical solutions for two fundamental models of light-matter interaction, the dispersive Jaynes-Cummings model and the quantum Rabi model, when the light mode is coupled to a thermal bath. A new steady-state formula shows the qubit's late-time population in the Rabi model depends on both qubit and cavity frequencies, unlike in the Jaynes-Cummings model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The non-thermal Rabi steady state Eq. (57) is a property of the local Lindblad equation (10); the paper's own Sec. II B 1 concedes the steady state of local master equations is debated, so the central physical claim may be an artifact of the local approximation.","rationale":"The reader's weakest assumption is exactly the local-vs-global master equation concern, and I agree that this is the most load-bearing issue. The internal perturbative calculation appears consistent: Eq. (57) follows from the multiscale solution displayed in Eq. (C6a), the numerics reproduce the local-equation steady state, and the exact dispersive-JC propagator is independently supported by comparison with QuTiP and with prior literature. I do not find an internal algebraic or logical error that would require rejection. However, the paper's headline physical claim—that the Rabi and Jaynes-Cummings models have non-perturbatively different steady states regardless of g—is not a theorem about the physical system unless the local master equation is the correct reduced description. The paper itself concedes in Sec. II B 1 that the steady-state issue is open and that the true Gibbs state may be reached on longer timescales. Because the paper explicitly frames its results as solutions of the local Lindblad equation, the mathematical core is secure, but the physical interpretation and scope of the claims should be conditional. The reader's CONDITIONAL verdict already captures this, so no change to the verdict is needed. The proposed global-master-equation test would settle whether Eq. (57) is physically realized or an artifact of the local approximation.","tokens_in":30505,"tokens_out":7352,"duration_ms":79893,"concrete_test":"Construct the global master equation for the same system-bath setup: truncate the bosonic Hilbert space, diagonalize H_Rabi, express the bath coupling a⊗B in the eigenbasis, and build Lindblad jump operators from the eigenoperator decomposition of a at temperature T with the same spectral density (or solve the Redfield equation without the secular approximation). Compute the steady-state ⟨σz⟩ for the parameters of Fig. 5 (ω=1, Ω=1.5, γ=T=0.1, g=0.1 and g=0.001). If the global steady state tends to the thermal value -1/(1+2n̄) as g→0 instead of Eq. (57), the non-perturbative difference is an artifact of the local master equation; if it reproduces Eq. (57), the central claim survives the local-vs-global objection.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the O(1) non-thermal steady state of the open Rabi model, Eq. (57), which is derived from the local Lindblad master equation (10) with bare jump operators a and a†. The paper explicitly acknowledges in Sec. II B 1 that the steady state of local master equations is a subject of ongoing debate and that a global master equation with dressed jump operators may give different late-time predictions. For the present system, the local dissipator treats all frequency components of a(t) under H_Rabi as a single decay channel, whereas a global (or Redfield) treatment resolves components at frequencies Ω and Ω±ω with bath-dependent rates. The balance of these rates is what selects the qubit steady state; changing the master equation can therefore change the O(1) steady-state polarization. The numerical checks in Figs. 5 and 6 validate Eq. (57) only against the same local Lindblad equation, not against a microscopically derived global equation. Thus the paper's strongest physical conclusion—a non-perturbative, g-independent difference between Rabi and Jaynes-Cummings late-time predictions—is load-bearing on the validity of the local approximation, and the manuscript itself flags this as unsettled.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript solves two open qubit-boson models within the local Lindblad master equation (10). For the dispersive Jaynes-Cummings model it obtains exact Bargmann-space propagation kernels, Eqs. (30)-(35), and derives qubit and boson dynamics for thermal and coherent initial states, recovering earlier zero-temperature and known-limit results. For the quantum Rabi model it develops perturbation theory in the qubit-boson coupling g up to second order, using multiscale perturbation theory to remove secular terms, and extracts a zeroth-order steady-state formula for the qubit, Eq. (57), that differs from the Jaynes-Cummings thermal value, Eq. (58). The analytical results are benchmarked against numerical solutions of the same local Lindblad equation in Figs. 2-6.","tokens_in":30688,"tokens_out":8836,"duration_ms":88218,"significance":"The exact dispersive-JC propagator is a useful technical advance: it removes earlier restrictions to coherent or thermal boson initial states and provides an analytical benchmark for truncation in numerical simulations. The Rabi steady-state formula Eq. (57) is an explicit, falsifiable prediction within the local-Lindblad model, and the numerical agreement in Figs. 4-6 supports the internal consistency of the perturbative calculation. The strengths of the paper are the detailed appendix derivations, the recovery of known limits (Refs. [21,24]), the absence of fitted parameters in the steady-state selection, and the clear statement of the assumptions (separable initial state; thermal boson state at second order). The main unresolved question is scope: the central late-time claim is a property of the local master equation with bare jump operators, and the manuscript itself flags the local/global master-equation debate in Sec. II B 1.","major_comments":[{"comment":"The O(1) non-thermal Rabi steady state, Eq. (57), is derived from the local Lindblad equation (10) with bare jump operators a and a†. The manuscript explicitly concedes in Sec. II B 1 that the steady state of local master equations is debated and that a global master equation with dressed jump operators can give different predictions. The numerical checks in Figs. 5 and 6 integrate the same local equation, so they do not test whether Eq. (57) survives a microscopically derived global/Redfield treatment, where the dissipator resolves bath correlations at frequencies Ω and Ω±ω instead of a single decay channel. Since the stated highlight of the paper (a non-perturbative difference between the late-time predictions of the Jaynes-Cummings and Rabi models) is load-bearing on this choice, I ask the authors either to compare Eq. (57) with the steady state of a global master equation for H_Rabi at leading order in g, or to reframe Eqs. (57)-(59) and the corresponding conclusions as statements about the local Lindblad dynamics only, removing the implication that they are generic late-time predictions of the open Rabi model.","section":"Sec. IV B 2, after Eq. (55)"},{"comment":"The manuscript states that 'we cannot rule out the presence of other terms of order g^2 that might arise from solvability conditions involved in higher orders of the multiscale perturbation theory.' This means the second-order dynamical solution is not proven to be complete at O(g^2). The numerical agreement in Fig. 4 supports the result but does not by itself close the analytic gap. Please either supply the missing solvability analysis or a bound on the omitted terms, or explicitly characterize Eq. (55) as a partial second-order result whose completeness is verified numerically rather than analytically.","section":"Section IV B 2, after Eq. (55)"}],"minor_comments":[{"comment":"The heading 'Pertubative solution' contains a typo; it should read 'Perturbative solution'.","section":"Sec. IV B 2, heading"},{"comment":"The sentence 'This will yield' is ambiguous because Eq. (59) contains both the zeroth-order term and a g^2 correction, and the preceding text contains the typo 'inifinite' for 'infinite'; please rephrase.","section":"Before Eq. (59)"},{"comment":"The display of the 4x4 matrix H^(B)(t) is difficult to parse as typeset; standard aligned matrix formatting or explicit row/column labels would make the entries easier to check.","section":"Eq. (35)"},{"comment":"The phrase 'non-perturbative difference' is potentially confusing, since Eq. (57) is the O(1) term of a perturbative expansion and receives O(g^2) corrections in Eq. (59); consider wording such as 'a difference that does not vanish as g → 0' instead.","section":"Sec. IV B 3, after Eq. (58)"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically competent and the algebraic core seems sound. The main risk is not an internal inconsistency but an overreach in the physical interpretation of the local-Lindblad steady state. Since the authors themselves acknowledge the local/global debate, the appropriate revision is to temper the abstract and conclusions or to add a concrete global-master-equation comparison. The incomplete second-order multiscale result, also admitted in the text, should be either closed or explicitly labeled as partial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth reading, but I'd advise not taking the flashiest claim at face value. The exact propagation kernel for the open dispersive Jaynes-Cummings model, Eqs. (30)-(35), is the real meat: it extends previous coherent- and thermal-initial-state results to arbitrary separable initial states, recovers known limits, and matches numerics. That part is solid and should be cited. The perturbative open-Rabi solution up to second order is also a substantive technical achievement, with the multiscale treatment of secular terms and the explicit coefficients in Appendices B and C. The paper gives credit where the work is hard and the derivations are detailed.\n\nNow the soft spots, in proportion. The headline physical claim, Eq. (57), is a property of the local Lindblad equation with bare jump operators, not of the Rabi model as a physical system. The authors themselves concede in Sec. II B 1 that the steady state of local master equations is debated and that a longer timescale may drive the system to the Gibbs state. The stress-test note gets this right: the numerical checks in Figs. 5 and 6 validate Eq. (57) only against the same local equation, so they cannot settle whether the Rabi-versus-JC steady-state difference is physical or an artifact of the local approximation. The paper should phrase Eq. (57) as a prediction of the local Lindblad dynamics, not as a non-perturbative statement about the Rabi model. Calling the difference \"non-perturbative\" is also loose: it is a zeroth-order-in-g quantity selected by a second-order consistency condition.\n\nThere are two smaller issues worth flagging for any referee. First, the paper says Eqs. (57) and (58) are \"valid also at resonance,\" but the perturbative solution is derived away from resonance and the secular-term treatment relies on the off-resonance structure. That claim needs an actual argument, not an assertion. Second, the multiscale treatment explicitly says it cannot rule out missing O(g^2) terms from higher-order solvability conditions, which is honest but means the improved second-order expression in Eq. (55) is not fully controlled. Neither issue kills the paper, but both should be tightened in revision.\n\nThe exact dispersive-JC result alone justifies engaging with this paper. I'd send it to a serious referee, with instructions to focus on the local-versus-global steady-state question and the resonance claim. If those are addressed honestly, it is a publishable contribution.\n\nFor my own work: I would cite the dispersive-JC propagator, and I'd bring it to a reading group if the discussion is about local master equations. Eq. (57) I would cite only with a caveat.","headline":"A genuinely useful exact solution for the dispersive JC model, plus a Rabi steady-state formula whose physical significance depends on how much you trust the local Lindblad equation.","tokens_in":31290,"tokens_out":1718,"would_cite":true,"duration_ms":20720,"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":"The open quantum Rabi model has a unique qubit steady state that is not thermal even at zero coupling, while the Jaynes-Cummings model's steady state is thermal.","keywords":["open quantum Rabi model","dispersive Jaynes-Cummings model","local Lindblad master equation","Bargmann space","holomorphic formalism","multiscale perturbation theory","steady state","qubit-boson coupling"],"falsifier":"Compute the long-time $\\langle\\sigma_z\\rangle$ for the same Hamiltonian and bath using a global master equation (dressed jump operators, or a Redfield or mean-force Gibbs treatment) in the limit $g\\to 0$ with fixed $\\omega$, $\\Omega$, $\\gamma$, and $T$; if the result approaches the Jaynes-Cummings thermal value $-1/(1+2\\bar n)$ instead of $-8\\omega\\Omega/[(1+2\\bar n)(\\gamma^2+4(\\omega^2+\\Omega^2))]$, the claimed non-perturbative difference collapses. Equivalently, a circuit-QED experiment measuring qubit polarization after a long thermalization should see a deviation from thermal equilibrium that depends on $\\Omega$ if the paper's claim is correct.","tokens_in":30248,"feed_emoji":"⚛️","tokens_out":6104,"duration_ms":67302,"temperature":0.7,"pith_summary":"The paper derives exact and approximate analytical solutions for two foundational open light-matter models: a qubit coupled to a single bosonic mode that is damped by a thermal bath. In the dispersive Jaynes-Cummings model, the authors obtain the most general propagator for any separable initial state, assuming only that the qubit-boson system starts factorized. For the quantum Rabi model, they construct a second-order perturbative solution in the qubit-boson coupling $g$, regularized by multiscale perturbation theory to eliminate secular terms. The central result is that the two models, usually regarded as equivalent for small $g$, have genuinely different steady states: the Rabi qubit's polarization at zeroth order depends on both qubit and boson frequencies, while the Jaynes-Cummings qubit relaxes to the thermal Gibbs value. This difference is non-perturbative in $g$, since it survives in the limit $g\\to 0$.","feed_headline":"Open Rabi qubit's late-time steady state is not thermal","feed_subtitle":"At vanishing coupling, the Rabi polarization depends on both frequencies; the Jaynes-Cummings model stays thermal.","key_machinery":"The mathematical engine is the holomorphic (Bargmann) representation of bosonic states, in which creation and annihilation operators become multiplication by $z$ and differentiation with respect to $z$. The master equation becomes a set of complex partial differential equations for $2\\times 2$ propagator kernels. For the dispersive Jaynes-Cummings model these PDEs are uncoupled and are solved exactly by a Gaussian ansatz, yielding closed-form kernels (Eqs. (30)-(35)). For the Rabi model the PDEs couple, and the solution is built as a perturbation series in $g$; the secular (time-divergent) terms at second order are removed by promoting the qubit's initial populations to slow functions of $\\tau = g^2 t$, a multiscale perturbation technique. The steady state is then obtained by solving the static second-order equations, which yields a consistency condition that fixes the zeroth-order populations of the qubit.","core_discovery":"The paper claims that the steady state of the local Lindblad master equation for the open quantum Rabi model is unique and non-thermal already at zeroth order in the qubit-boson coupling. Specifically, Eq. (57) gives $\\lim_{t\\to\\infty}\\langle\\sigma_z\\rangle^{(\\text{Rabi})} = -8\\omega\\Omega/[(1+2\\bar n)(\\gamma^2+4(\\omega^2+\\Omega^2))] + O(g^2)$, whereas the non-dispersive Jaynes-Cummings model yields the thermal value $-1/(1+2\\bar n)$ of Eq. (58). Because both expressions are evaluated at $g=0$, the difference between the two models is independent of $g$ and therefore cannot be removed by making the coupling arbitrarily small. The paper argues that this sharp contrast arises because the Rabi Hamiltonian does not commute with $\\sigma_z$, so the bath-mediated dissipation acts on the qubit even though the bath couples only to the bosonic mode. The result is stated within the local Lindblad approximation, whose validity is the subject of an open debate that the paper itself acknowledges.","pith_inferences":["If the local master equation's steady state is later shown to be an artifact of the local approximation in certain parameter regimes, the same Bargmann-kernel machinery can be rerun with global, dressed jump operators; the sharpest testable signature would be whether the $\\Omega$-dependence of the qubit polarization survives in a Redfield or dressed-operator treatment.","The difference may have practical consequences for dispersive qubit readout in circuit QED, where the frequency-dependent steady-state shift could serve as a calibration signal or as a probe of which master-equation description is physically realized.","The same holomorphic-kernel method could be extended to multiple qubits or to additional qubit dissipation channels; one would expect the non-thermal correction to appear whenever $[H,\\sigma_z]\\neq 0$, that is, whenever counter-rotating or transverse terms open a bath-mediated decay path for the qubit.","An experiment with a tunable resonator frequency $\\Omega$ at fixed qubit frequency $\\omega$ could map the predicted curve of steady-state $\\langle\\sigma_z\\rangle$ versus $\\Omega$ and distinguish the Rabi prediction from the flat Jaynes-Cummings prediction even when $g$ is too small to detect dynamically."],"forward_implications":["The open quantum Rabi model, even at infinitesimal coupling, has a qubit steady state that is diagonal but not Gibbs, with populations set by $\\omega$, $\\Omega$, $\\gamma$, and the bath temperature.","The Jaynes-Cummings and Rabi models cannot both be used interchangeably at late times, even in the rotating-wave regime where their early-time dynamics agree.","The exact dispersive Jaynes-Cummings propagator provides a closed form for qubit coherence decay, giving an analytic expression for the dephasing rate $\\Gamma_2$ and a benchmark for truncated numerical simulations.","The second-order multiscale solution supplies a finite-time expression for $\\langle\\sigma_z\\rangle$ that remains accurate at intermediate and long times, not only asymptotically.","Since the difference is independent of $g$, any experiment probing the steady state at weak coupling can distinguish the two models without entering the ultra-strong coupling regime."],"supporting_citations":[{"why":"Supplies the standard weak-coupling derivation of the Lindblad master equation and the complete-positivity framework on which the paper's local master equation rests.","marker":"[17]"},{"why":"Justifies the perturbative treatment of inter-site couplings in the local description of open quantum networks, the approximation the paper uses for bare jump operators.","marker":"[18]"},{"why":"Compares local versus global master equations and is cited for the open debate on steady states, which is the context for the paper's steady-state claim.","marker":"[19]"},{"why":"Provides the numerical benchmark and parameter set for the dispersive Jaynes-Cummings comparison shown in Fig. 2, against which the exact analytical solution is validated.","marker":"[15]"},{"why":"Contains the earlier analytical solution of the dispersive Jaynes-Cummings model at zero temperature with coherent initial states, which the paper generalizes to arbitrary separable initial states.","marker":"[21]"},{"why":"Gives the finite-temperature Wigner-function solution for a thermal boson initial state; the paper's coherence-decay rate matches their Eq. (44).","marker":"[24]"},{"why":"Shows that the local master equation for the Jaynes-Cummings/Rabi model has a non-Gibbs steady state, the phenomenon the paper quantifies explicitly.","marker":"[49]"},{"why":"Discusses open-system dynamics and the mean-force Gibbs state, cited for the ongoing debate about steady states and the timescales on which the local equation is valid.","marker":"[51]"}],"fun_headline_variants":["Rabi steady state stays non-thermal even at zero coupling","Zero-coupling Rabi still avoids thermal equilibrium","Jaynes-Cummings thermal but Rabi isn't even at g=0","Rabi qubit's steady state defies thermal bath at g=0","Open Rabi model: non-thermal steady state even without coupling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All results rest on treating the open dynamics with a local Lindblad master equation whose jump operators are the bare bosonic operators $a$ and $a^\\dagger$; if a global master equation with dressed jump operators is the more faithful description for weak qubit-boson coupling, the predicted Rabi steady state of Eq. (57) could be an artifact rather than a physical property.","fun_headline_variants_meta":{"raw":{"variants":["Rabi steady state stays non-thermal even at zero coupling","Zero-coupling Rabi still avoids thermal equilibrium","Jaynes-Cummings thermal but Rabi isn't even at g=0","Rabi qubit's steady state defies thermal bath at g=0","Open Rabi model: non-thermal steady state even without coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1386,"prompt_tokens":1065,"completion_tokens":321,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":230}},"tokens_in":681,"tokens_out":321,"duration_ms":3212,"temperature":1.0,"reasoning_tokens":230,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:19:27.531166+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the long-time $\\langle\\sigma_z\\rangle$ for the same Hamiltonian and bath using a global master equation (dressed jump operators, or a Redfield or mean-force Gibbs treatment) in the limit $g\\to 0$ with fixed $\\omega$, $\\Omega$, $\\gamma$, and $T$; if the result approaches the Jaynes-Cummings thermal value $-1/(1+2\\bar n)$ instead of $-8\\omega\\Omega/[(1+2\\bar n)(\\gamma^2+4(\\omega^2+\\Omega^2))]$, the claimed non-perturbative difference collapses. Equivalently, a circuit-QED experiment measuring qubit polarization after a long thermalization should see a deviation from thermal equilibrium that depends on $\\Omega$ if the paper's claim is correct.","supporting_citations":[{"cited_title":"Jaynes and F","cited_arxiv_id":null,"evidence_quote":"Supplies the standard weak-coupling derivation of the Lindblad master equation and the complete-positivity framework on which the paper's local master equation rests."},{"cited_title":"Larson and T","cited_arxiv_id":null,"evidence_quote":"Compares local versus global master equations and is cited for the open debate on steady states, which is the context for the paper's steady-state claim."},{"cited_title":"Larson, T","cited_arxiv_id":null,"evidence_quote":"Contains the earlier analytical solution of the dispersive Jaynes-Cummings model at zero temperature with coherent initial states, which the paper generalizes to arbitrary separable initial states."},{"cited_title":"Haroche and J.-M","cited_arxiv_id":null,"evidence_quote":"Gives the finite-temperature Wigner-function solution for a thermal boson initial state; the paper's coherence-decay rate matches their Eq. (44)."},{"cited_title":"Bargmann, On a hilbert space of analytic functions and an associated integral transform part i, Commun","cited_arxiv_id":null,"evidence_quote":"Shows that the local master equation for the Jaynes-Cummings/Rabi model has a non-Gibbs steady state, the phenomenon the paper quantifies explicitly."},{"cited_title":"Verhulst, Methods and Applications of Singular Perturba- tions: Boundary Layers and Multiple Timescale Dynamics , Texts in Applied Mathematics (Springer New York, 2006)","cited_arxiv_id":null,"evidence_quote":"Discusses open-system dynamics and the mean-force Gibbs state, cited for the ongoing debate about steady states and the timescales on which the local equation is valid."}],"review_version":1}