{"id":"172f1a80-ec6e-4c32-9a9e-c456ccc42069","arxiv_id":"2511.11903","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A polaron-frame truncation to a few cavity levels reproduces exact Rabi spectra and Dicke–Heisenberg phase diagrams with far fewer photon states than bare truncation.","lead":"An entangling transformation plus a few-level truncation reduces strongly coupled light–matter problems to small, closed-form Hamiltonians that match exact spectra. The approach could make ultrastrong-coupling simulations in polaritonic chemistry and cavity materials much cheaper.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CDH convergence for observables is governed by M_P, not M; SI documents non-monotonic M-convergence and M_P=20, undermining the 'systematically convergent' central claim.","rationale":"The reader's weakest assumption correctly identified the convergence of the CDH series for arbitrary strongly coupled systems and observables as the load-bearing point. My analysis of the SI confirms this: Sec. I C reports non-monotonic convergence with M and the need for M_P=20 for ground-state observables at resonance, while the main text claims systematically convergent observables. The spectral results are credible and well-supported by exact Rabi comparisons, so the paper merits conditional acceptance pending clarification or error bounds, exactly as the reader concluded. No new independent objection found; the conditional verdict stands.","tokens_in":23812,"tokens_out":13559,"duration_ms":115370,"concrete_test":"Test: for the resonant quantum Rabi model (Ω=2Δ), compute the ground-state magnetization ⟨σz⟩_GS at λ/Ω=2 using CDH with M=1,2,3,4 and M_P=M, and compare to exact and to M_P=20. If the M_P=M results do not improve monotonically with M and deviate by >5% while M_P=20 matches, then the M-series does not converge to the exact observable limit on its own, confirming the SI's caveat and invalidating the abstract's unqualified 'systematically convergent' claim for ground-state observables.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The letter's central claim (abstract, Introduction) is that truncating the polaron-frame Hamiltonian to M excitation sectors yields a systematically convergent series of compact CDHs that 'outpace conventional approaches' for spectra and observables. The spectral benchmarks (Fig. 2) support convergence for the first six Rabi eigenenergies with M=4. However, the supplementary material (Sec. I C, Figs. 4-5) explicitly states that for ground-state observables in 'difficult' parameter ranges (low temperature, resonance) convergence with M is non-monotonic and one must use M_P=20 for the operator rotation—the same dimension needed for bare convergence. Thus the accuracy of matter observables is not controlled by the compact M parameter alone; it depends on a separately chosen M_P that is not part of the CDH construction and is fixed by external bare-convergence knowledge. The abstract claims 'converging ground-state and thermal observables' and 'systematically convergent'; the SI's own evidence shows this is false for ground-state observables unless M_P is chosen much larger than M. No proof or error bound is given to show that the M-series converges for arbitrary couplings or observables. This is the load-bearing gap: if the user cannot rely on low-M CDH alone for observables, the advertised compactness and systematic improvability do not hold in the claimed generality.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a cavity-dressed Hamiltonian (CDH) framework for strongly coupled light-matter systems. The method applies a multimode polaron transformation and then truncates the transformed Hamiltonian to the M lowest Fock states, yielding closed-form effective Hamiltonians for matter observables and spectra. The approach is applied to the quantum Rabi model, where closed-form M=2 and M=4 CDHs are given, and to a Dicke-Heisenberg spin chain, where M=1 and M=3 CDHs reproduce phase diagrams and structure factors. The central claim is that the finite-M CDHs form a systematically convergent series that outperforms conventional bare-truncation approaches from weak through deep-strong coupling.","tokens_in":24169,"tokens_out":3897,"duration_ms":37365,"significance":"If the convergence claim can be substantiated, the CDH framework is a useful non-perturbative tool for ultrastrong-coupling cavity QED: it is analytically transparent, computationally cheap, and benchmarked against the exact Braak solution and brute-force diagonalization. The closed-form expressions and the explicit comparison with bare truncation are valuable. However, the manuscript's advertised scope—'systematically convergent' and 'arbitrarily strong coupling'—is broader than what the numerical benchmarks and the supplementary material actually establish, particularly for ground-state observables.","major_comments":[{"comment":"The abstract and Introduction claim that the CDH series is 'systematically convergent' for spectra and ground-state/thermal observables. The supplementary material itself states that for ground-state observables in 'difficult' parameter ranges (low temperature, resonance) convergence in M is non-monotonic and one must use M_P=20 for the operator rotation—the same dimension required for bare convergence. This directly contradicts the advertised compactness and systematic improvability for observables, since accuracy is not controlled by M alone. Please either restrict the central claim to spectral properties, or provide a rigorous convergence criterion/protocol for choosing M_P, or state the observable-convergence claim as a conjecture with explicit numerical support.","section":"SI Sec. I C; Fig. 5"},{"comment":"The CDH mapping is a unitary transformation, but the computational scheme uses a Hamiltonian truncated to M levels while the rotated operator in the expectation value is truncated separately at M_P levels. There is no proof or error bound that the sequence of such hybrid truncations converges to the exact limit for arbitrary couplings and observables. The benchmarks cover two models. The title 'arbitrarily strong light-matter coupling' and the statement of exactness in the deep-strong limit require either a proof for a defined class of Hamiltonians/observables or a clear qualification that these are empirical convergence properties, not established general theorems.","section":"Eq. (2), Eq. (34), SI Sec. I C"},{"comment":"The main-text ground-state phase diagrams in Fig. 3 are computed with M_P=M, as stated in the caption, while SI Sec. I C reports that at resonance ground-state observables require M_P=20 to obtain physical behavior in the deep-strong limit. This raises a concrete concern about whether the Dicke-Heisenberg phase diagrams, which extend to strong coupling, are quantitatively reliable in the deep-strong region. Please clarify whether Fig. 3 was produced with M_P=M or M_P=20, and discuss the sensitivity of the reported phase boundaries to M_P.","section":"Fig. 3 caption; SI Sec. I C"}],"minor_comments":[{"comment":"The notation E_{±,±′} is not defined precisely; the expression has two nested square roots, so the reader has to infer how the four sign combinations are assigned. Please spell this out or label the states explicitly.","section":"Eq. (6)"},{"comment":"The phrase 'arbitrarily strong light-matter coupling' is stronger than the quantitative benchmarks: Fig. 2 covers λ/Δ up to 5 for spectra, and the deep-strong limit is only discussed qualitatively. Please qualify the scope in the abstract or cite the specific benchmark ranges.","section":"Abstract/Introduction"},{"comment":"There is a typo in the supplementary title: 'Dicke-Hiesenberg' should be 'Dicke-Heisenberg'.","section":"SI title"},{"comment":"The main text says 'For ground state observables, M_P > M was required to achieve convergence' but does not reconcile this with the earlier assertion that the method is systematically convergent for observables. This sentence should be expanded or moved to qualify the claim.","section":"Main text, paragraph after Eq. (9)"}],"recommendation":"major_revision","confidential_remarks":"The derivation and benchmarks are solid, and the method is potentially publishable. The key issue is scope: the observable-convergence claim is contradicted by the authors' own supplementary material, so the revision needs to either tighten the claims or add a theoretical error bound/protocol. The paper is also heavily dependent on the companion 'Supplemental material' reference; I would advise the editor to ensure the SI is self-contained enough for referees."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague —\n\nWorth a look, with one important caveat. The paper builds a 'cavity-dressed Hamiltonian' from a polaron transform plus Fock-block truncation and shows that low-order truncations reproduce the Rabi spectrum well beyond the ultrastrong regime. The benchmarks against Braak's exact solution are genuine, and the closed-form M=2, M=3, M=4 Hamiltonians for Rabi and Dicke-Heisenberg models are useful results. The general block formula for noncommuting operators and the dressing functions are new enough to matter. This is not just a restatement of De Bernardis' polaron-frame work; the systematic truncation viewpoint and the lattice-model formulas are the contribution.\n\nThe soft spot is the 'systematically convergent' language. It is true for the eigenenergies shown in Fig. 2 — the M=4 CDH tracks the first six levels far better than bare N=10 truncation. But for ground-state observables, the supplementary material shows the operator rotation needs M_P up to 20 — the same dimension as bare convergence — and that M-convergence is non-monotonic at resonance. The paper discloses this in the SI and mentions M_P>M in the text, but the abstract and the title's 'arbitrarily strong light-matter coupling' overstate the method's compactness for observable calculations. The method is still convergent; it just loses some of its advertised advantage when you want ground-state properties rather than spectra.\n\nOther soft spots: no code or data, no error bounds for general couplings, and the claimed extensions to multimode, noncommuting, and leaky cavities are asserted in the discussion but not benchmarked. The phase diagrams for the Dicke-Heisenberg chain look plausible and consistent with the bare results, and the identification of the antiferromagnetic-superradiant phase via the CDH is a nice bonus.\n\nWho should read this: anyone building effective Hamiltonians for ultrastrong-coupling cavity QED, polaritonic chemistry, or spin-cavity models. The spectral method deserves referee time. A referee should ask for a more nuanced abstract, a convergence test with noncommuting coupling operators, and ideally a release of the code used for the benchmarks.","headline":"Polaron-truncation method that shines for spectra, overclaims for ground-state observables.","tokens_in":24595,"tokens_out":3806,"would_cite":true,"duration_ms":31959,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.50.Pq"],"model":"deepseek-v4-flash","headline":"The paper claims a polaron-style displacement followed by a small photon-number truncation yields compact Hamiltonians that remain accurate from weak to deep-strong light-matter coupling.","keywords":["cavity-dressed Hamiltonian","ultrastrong coupling","deep strong coupling","quantum Rabi model","Dicke-Heisenberg model","polaron transformation","non-perturbative mapping","strong light-matter coupling"],"falsifier":"Compute the exact ground-state photon number ⟨a†a⟩ for the single-mode Rabi model at resonance (Ω=2Δ) and λ/Δ=3, where the paper's M=4 CDH is accurate for the first six energies; if the CDH prediction with M_P=20 does not approach the exact value as M grows from 1 to 6, the claim that the CDH framework captures strongly coupled light-matter systems at arbitrary coupling would be refuted for cavity observables. Alternatively, a simple two-mode model with noncommuting coupling operators at strong coupling would test the general multimode claim: if increasing M fails to reduce the ground-state en","tokens_in":23732,"feed_emoji":"⚛️","tokens_out":4187,"duration_ms":35953,"temperature":0.7,"pith_summary":"The paper claims that strongly coupled light-matter systems can be described by a compact, closed-form Hamiltonian obtained by applying an entangling (polaron) displacement to the photons and then keeping only the lowest M photon-number sectors. It shows for the quantum Rabi model that this cavity-dressed Hamiltonian with M=4 reproduces the first six eigenenergies across weak, ultrastrong, and deep-strong coupling, whereas a bare Fock-truncated Hamiltonian with N=10 diverges. For the Dicke-Heisenberg spin chain, low-degree CDHs reproduce magnetization phase diagrams, entanglement entropy, and structure factors with a small fraction of the Hilbert-space dimension used in brute-force simulations. The method is claimed to be systematically convergent and to extend to multimode and leaky cavities, providing analytical forms that reveal the mechanisms of cavity-mediated interactions.","feed_headline":"Four dressed levels beat ten bare levels at ultrastrong coupling","feed_subtitle":"A compact truncation of a polaron-dressed Rabi Hamiltonian stays accurate where bare Fock truncation diverges.","key_machinery":"The central object is the cavity-dressed Hamiltonian (CDH): after the polaron displacement U_P is applied, the transformed Hamiltonian is decomposed into blocks labeled by photon occupation numbers, and all but the lowest M sectors are discarded. The key identities are closed-form matrix elements of rotated operators in the number-state basis, obtained either by momentum-space integrals over Hermite polynomials or by spectral decomposition of the coupling operator into displacement operators. These give analytical expressions for every block, so the CDH is a finite-dimensional, systematically improvable effective Hamiltonian rather than a numerically constructed one.","core_discovery":"On its own terms, the paper's central claim is that the unitary transformation U_P = exp[Σ (λ_n/Ω_n) S_n(a_n† − a_n)] dresses the matter operators with photon displacements, and that truncating the transformed Hamiltonian to M photon-number sectors yields a series of cavity-dressed Hamiltonians that converge to the exact spectrum and equilibrium observables at arbitrarily strong coupling. The evidence is construction plus numerical benchmarking: closed-form M=1,2,3,4 CDHs for the Rabi model, a closed-form M=3 CDH for the Dicke-Heisenberg model, and the demonstration that resonant and ultrastrong regimes, where bare truncation fails badly, are captured with small M. The paper does not claim a","pith_inferences":["The strongest untested implication is that convergence extends to arbitrary observables; the paper itself computes no cavity observables, and ground-state photon statistics are the natural stress test. If CDH fails there, the 'arbitrarily strong coupling' claim should be read as limited to matter observables and spectra.","Because the dressed Hamiltonian is an effective spin model with all-to-all interactions, the same closed forms could be fed into tensor-network or mean-field treatments of longer chains and higher dimensions; the paper's L=8 numerics suggest the advantage grows with L.","The non-monotonic convergence reported for difficult parameters hints that a rigorous error bound in terms of M, λ/Ω, and the operator norm of S_n would be needed to turn the recipe into a method with guaranteed accuracy; without it, users must benchmark each new model.","For open systems, combining the leaky-cavity CDH with standard dissipative evolution of the dressed Hamiltonian may give a practical route to ultrastrong-coupling dynamics, which the paper lists as future work."],"forward_implications":["Rabi-model spectra in the resonant ultrastrong regime can be computed with M=4 dressed levels where N=10 bare levels diverge, cutting the Hilbert-space dimension needed for comparable accuracy.","Thermal and ground-state observables such as magnetization converge with increasing M; ground-state calculations may require a larger rotation space M_P (up to 20) even when M=4 suffices for energies.","For the Dicke-Heisenberg chain, M=1 already captures the ferromagnetic-paramagnetic crossover and M=3 reproduces structure factors, so phase diagrams can be mapped with a small effective dimension.","The mapping exposes the physics: an effective spin splitting suppressed by exp(−2λ²/Ω²), a cavity-mediated all-to-all spin coupling in the x direction, and renormalization and mixing of Heisenberg couplings with closed-form dressing functions.","The same construction applies to multimode cavities with noncommuting coupling operators and to leaky cavities, extending the method beyond the benchmarked single-mode models."],"fun_headline_variants":["Dressed photon sectors converge at ultrastrong coupling","Compact dressed Hamiltonians tame strong light-matter coupling","New mapping beats bare truncation in strong coupling","Cavity-dressed truncation converges where bare fails","Dressed levels beat bare at ultrastrong coupling"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that truncating the dressed Hamiltonian to the M lowest photon sectors converges to the exact limit for every strongly coupled system and every observable of interest; this is asserted and shown numerically for two models, but not proven, and the paper's own supplementary material documents non-monotonic convergence that forces M_P=20 for ground-state observables.","fun_headline_variants_meta":{"raw":{"variants":["Dressed photon sectors converge at ultrastrong coupling","Compact dressed Hamiltonians tame strong light-matter coupling","New mapping beats bare truncation in strong coupling","Cavity-dressed truncation converges where bare fails","Dressed levels beat bare at ultrastrong coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000143,"raw_usage":{"total_tokens":1013,"prompt_tokens":752,"completion_tokens":261,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":185}},"tokens_in":496,"tokens_out":261,"duration_ms":3177,"temperature":1.0,"reasoning_tokens":185,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:07:45.507207+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact ground-state photon number ⟨a†a⟩ for the single-mode Rabi model at resonance (Ω=2Δ) and λ/Δ=3, where the paper's M=4 CDH is accurate for the first six energies; if the CDH prediction with M_P=20 does not approach the exact value as M grows from 1 to 6, the claim that the CDH framework captures strongly coupled light-matter systems at arbitrary coupling would be refuted for cavity observables. Alternatively, a simple two-mode model with noncommuting coupling operators at strong coupling would test the general multimode claim: if increasing M fails to reduce the ground-state en","supporting_citations":[],"review_version":1}