{"id":"a7e5a56b-8bdc-4d93-ab23-762fe5961b37","arxiv_id":"2602.12787","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A closed-form quantum Fisher information for the multilevel quantum Rabi model in the adiabatic regime shows dark-manifold saturation approaches the ideal thermometer limit and bright-manifold saturation gives broadband, self-averaging thermometric response.","lead":"This paper derives a fast approximate formula for the thermometric sensitivity of a multilevel quantum Rabi model, where many nearly degenerate atomic levels couple to one cavity mode. It shows that loading the model with many 'dark' levels creates a sharp sensitivity peak near the ideal quantum thermometer, while loading it with many 'bright' levels creates a broad, disorder-robust response.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dark-sector fine structure is ignored in the large-D regime; at weak coupling the intraband width epsilon can exceed the peak thermal scale T* ~ E_BD/ln D, so the approach-to-ideal claim is unverified.","rationale":"The reader's weakest assumption correctly identifies the most load-bearing concern: the large-D dark-manifold results assume that O(epsilon) intraband structure in the dark sector is thermometrically irrelevant, and this is only validated for D_g=2, D_e=4. My analysis strengthens this concern by noting that the relevant thermal scale at the bright–dark peak, T* ~ E_BD/ln D, decreases with D, while the dark bandwidth epsilon is fixed. For weak coupling and parameters used in the paper, epsilon can exceed T* by an order of magnitude at D=1000, so the dark states are not effectively degenerate at the peak. The effective degeneracy then grows only as D/ln D, which would slow the approach to the ideal limit. Since Fig. 5 uses epsilon=0 and Fig. 6 (right) relies on the AA formula that already ignores the fine structure, neither provides a check of the assumption at large D. The proposed exact-diagonalization test for moderate D would directly probe whether the fine-structure insensitivity holds and whether the deviation worsens with D. If the test passes, the conditional verdict could be upgraded; if it fails, the central claim is significantly weakened. The paper is otherwise clear and technically sound, with small-system validation and public code, but the large-D claim remains conditional on this unresolved point.","tokens_in":19832,"tokens_out":7232,"duration_ms":64532,"concrete_test":"Brute-force diagonalize the MQRM for D_g=1 (or 2), D=10, 20, 50, with omega_a=0.2 omega_f, epsilon=0.1 omega_a and 0.25 omega_a, at weak coupling g=0.1 omega_f (and intermediate g=1.2 omega_f). Use a photon cutoff with N_max=50 or convergence checks. Compute the exact thermal QFI and the bright–dark peak F*_BD, and compare with the AA formula Eq. (30) evaluated with a degenerate dark manifold. If the exact peak differs from the AA prediction by more than 10% at D=20, or if the discrepancy grows with D, the dark fine-structure insensitivity assumption fails. Alternatively, tune epsilon from 0 to 0.5 omega_a and check whether the peak ratio F*_BD/F*_ideal degrades as epsilon/T* increases.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that dark-manifold saturation makes the bright–dark QFI peak approach the ideal thermometer limit relies on the assertion (Sec. II B) that “thermometric precision is insensitive to [dark-sector] fine structure, so we simply ignore it.” This is validated only for the D_g=2, D_e=4 system in Fig. 3. For the headline D=1000 cases (Figs. 5–6), the dark manifold is treated as exactly degenerate within the AA, and no exact-diagonalization check is provided. In the full model, intraband detunings of order epsilon spread the dark band over a width ~epsilon. The bright–dark peak occurs at T* ~ E_BD/ln D (from Eq. B2). For weak coupling, e.g., g=0.1 omega_f and omega_a=0.2 omega_f, E_BD ~ 0.01 omega_f, while epsilon=0.1 omega_a = 0.02 omega_f; for D=1000, T* ~ 0.0014 omega_f, so epsilon >> T*. Only a fraction ~T*/epsilon of the dark states are thermally active, so the effective degeneracy is D_eff ~ D * T*/epsilon, which grows as D/ln D rather than D. The peak QFI may thus be suppressed relative to the ideal benchmark. Fig. 5 sets epsilon=0, so it cannot test this; Fig. 6 (right) includes epsilon=0.25 omega_a but evaluates the AA formula, which ignores the fine structure, so it cannot validate the approximation. If the fine-structure insensitivity fails at large D, the central promise of near-ideal dark-saturated thermometry is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies equilibrium thermometry in a multilevel quantum Rabi model (MQRM) with two near-degenerate atomic manifolds coupled to a single cavity mode. Working in the adiabatic regime (ω_a/ω_f ≪ 1), the authors perform an SVD of the coupling matrix and derive an approximate spectrum: bright doublets with energies given by Eq. (21) and a dark sector treated as exactly degenerate and decoupled, Eqs. (15)–(16). From this spectrum they obtain a closed-form thermal QFI, Eqs. (30)–(31), decomposed into intra-doublet, bright–bright (Eq. (39)), and bright–dark (Eq. (37)) contributions, together with an energy-ordered truncation scheme for numerical evaluation. The approximation is validated against exact diagonalization for a D_g=2, D_e=4 system in Figs. 2–3. The paper then analyzes two large-manifold limits: dark-manifold saturation (D_e ≫ D_g), where the bright–dark QFI peak is claimed to approach the ideal thermometer benchmark as D grows (Fig. 5), and bright-manifold saturation (D_g = D_e = M), where the response is claimed to be broadband and self-averaging (Fig. 7). The central assertion is that these regimes survive for large manifolds and finite intraband detunings, making the MQRM a versatile and sensitive equilibrium thermometer.","tokens_in":20310,"tokens_out":15266,"duration_ms":134128,"significance":"The paper's analytical QFI formula and its decomposition into physically distinct transition classes are valuable and may provide design principles for cavity-QED thermometry. Concrete strengths include the closed-form expression, the validation against exact diagonalization in Figs. 2–3, the Laguerre–Wishart modal-spectrum validation in Fig. 9, and the public data/code availability. The ideal-thermometer comparison is meaningful because matching the peak temperature does not fix the peak height, so the approach-to-ideal claim is not forced. However, the headline large-manifold claims (D=1000, M=250) rest on two uncontrolled approximations—the neglect of dark-sector fine structure and the energy-ordered truncation—neither of which is directly validated at those system sizes. The bright-manifold saturation claim also relies on the SVD-detuning approximation at non-negligible intraband spreads. These gaps need to be addressed before the large-D conclusions can be regarded as established.","major_comments":[{"comment":"The statement that dark-sector fine structure is thermometrically irrelevant is load-bearing for the dark-manifold saturation claim, but it is validated only for D_g=2, D_e=4 (Fig. 3). The D=1000 cases are not checked against exact diagonalization: Fig. 5 uses ε=0, and Fig. 6(right), although it includes ε=0.25ω_a, evaluates the AA formula that ignores the dark fine structure. This matters because the bright–dark peak temperature scales as T* ~ E_BD/ln D (Eq. B2). For weak coupling with ω_a=0.2ω_f and ε=0.25ω_a=0.05ω_f, D=1000 gives T*≈0.03ω_f, so ε exceeds T*. Then only a fraction ~T*/ε of the dark band is thermally active, reducing the effective degeneracy below D and potentially suppressing the peak QFI relative to the ε=0 benchmark. Please provide exact-diagonalization checks for intermediate D (e.g., D=10–50) with nonzero ε, or an explicit bound showing that the fine-structure corre","section":"Sec. II B, Eq. (15), Figs. 5–6"},{"comment":"The energy-ordered truncation cutoff Θ is introduced to make large-manifold computations feasible, and the text states that convergence was tested 'across all parameter sets tested.' However, the only documented convergence check is the small D_g=2, D_e=4 system in Fig. 3. For the D=1000 and M=250 results in Figs. 5–7 no Θ-dependence is reported. This is not a technicality: the bright-manifold saturation claim relies on many high-lying bright ladders contributing at intermediate/high temperatures, and an insufficient cutoff could artificially broaden or flatten the QFI profile. Please report the values of Θ used for each figure and demonstrate convergence with respect to Θ for representative large-manifold parameters.","section":"Appendix A, Eq. (A1)"},{"comment":"The bright-manifold saturation self-averaging result for M=250 is obtained entirely from the AA formula. The SVD-detuning approximation in Sec. II B retains only diagonal intraband terms in the bright sector (Eq. (14)) and neglects off-diagonal couplings between different bright doublets. At ε=0.25ω_a these neglected terms are not obviously small for the dense spectrum at M=250, and the small-system validation does not cover this regime. Please provide an intermediate-size exact-diagonalization check (e.g., M=6–10) with nonzero ε, or estimate the magnitude of the neglected inter-doublet terms and their effect on the QFI.","section":"Sec. IV.B, Fig. 7"}],"minor_comments":[{"comment":"The caption says the AA Formula retains terms up to the 5-th dark manifold; please state the corresponding truncation parameter Θ explicitly, and do the same for the large-manifold figures.","section":"Fig. 3 caption"},{"comment":"The notation S_2 and S_3 is overloaded: Eq. (31b) defines S_2(n) with both bright and dark contributions, while Eq. (37) uses S_2^B, S_3^B, etc. without defining these block-resolved moments. Please define the block-resolved quantities before Eq. (37).","section":"Eq. (31), Eq. (37)"},{"comment":"The effective ideal gap E_ideal_eff is fixed by requiring the ideal probe to peak at the same temperature as the MQRM. The text describes this procedure, but it should state more explicitly that E_ideal_eff is a calibration output, not an independently measured or free parameter of the model.","section":"Sec. IV.A.1"},{"comment":"The parameter matching for the Laguerre zeros should be checked carefully: for complex Wishart matrices (β=2) the joint density exponent is n−m, and the relation α'=n−m−2/β should be reconciled with the standard Laguerre parameter. This does not affect the numerical validation in Fig. 9, but a short clarification would help.","section":"Appendix C, Eq. (C7)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the analytical QFI formula is a useful contribution. My main reservation is that the headline large-manifold claims are supported by an approximation whose key assumptions—dark-sector fine-structure irrelevance and truncation convergence—are only validated on a D=2 system. The requested intermediate-size exact checks are feasible and should settle whether the near-ideal dark-saturated thermometry claim survives finite intraband structure. I do not think this requires rejection, but the additional validation should be required before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me cut to the part you'll care about. This paper is a legitimate step forward for equilibrium thermometry in multilevel cavity-QED. The authors derive a closed-form thermal QFI for the multilevel quantum Rabi model within the adiabatic approximation, and it separates cleanly into intra-doublet, inter-doublet bright-bright, and bright-dark contributions. That decomposition is new, and it lets them identify two design regimes: dark-manifold saturation, where a large nearly-degenerate dark sector boosts the peak QFI toward the ideal probe limit, and bright-manifold saturation, where many bright channels produce a broad, self-averaging response. I believe this is correct as far as the small-system checks go: they validate the AA spectrum and QFI against exact diagonalization for D_g=2, D_e=4 at weak and strong coupling, and the agreement is good. The energy-ordered truncation scheme is sensible, and the use of Laguerre-Wishart typical spectra is neat. They also post code and data.\n\nThe soft spot is exactly where the reader and stress-test point. The claim that dark-sector fine structure is thermometrically irrelevant (Sec. II B) is asserted and only checked for a system with two dark states. For the headline D=1000 dark-saturation runs, the dark band is taken as exactly degenerate, and no exact or numerically converged cross-check is provided. The stress-test's estimate is not nonsense: at g=0.1 omega_f, omega_a=0.2 omega_f, the primary bright-dark gap is about 0.01 omega_f, while the intraband width is epsilon=0.02 omega_f. The ideal peak sits at T* ~ E_BD/ln D ~ 0.0014 omega_f, so epsilon >> T*. In that situation only a fraction T*/epsilon of the dark states are thermally populated, and the effective degeneracy that sets the peak height is closer to D*T*/epsilon than to D. If that is right, the approach-to-ideal scaling is slowed or broken in the weak-coupling regime with non-zero intraband spread. The paper's Fig. 6 (right) includes epsilon=0.25 omega_a but at intermediate coupling, where E_BD is much larger than epsilon, so it cannot rule out the weak-coupling problem.\n\nI don't think this sinks the paper. The QFI formula and the two saturation mechanisms are still useful, and the intermediate-coupling dark-saturation claims appear robust. But the statement \"thermometric precision is insensitive to this fine structure\" is load-bearing and under-supported. The authors should either provide exact diagonalization checks at intermediate D (say D=50–200) with non-zero epsilon, or give an analytic argument for why a dark band of width epsilon does not reduce the effective degeneracy at the peak. If they can do that, this is a solid PRA-level contribution. As is, it deserves a serious referee, but the referee should push on this point.\n\nI'd bring it to a reading group if the group works on quantum thermometry or ultrastrong coupling; otherwise it's a specialist's paper. I'd probably cite it once the fine-structure question is settled, with a caveat.","headline":"Genuinely useful closed-form QFI decomposition for multilevel Rabi thermometry, but the claim that dark fine structure can be ignored at large D is under-supported and needs a targeted check.","tokens_in":20723,"tokens_out":5625,"would_cite":true,"duration_ms":46247,"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 temperature sensitivity of a multilevel quantum Rabi model is captured by an analytic formula that decomposes into bright and dark contributions, and saturating either manifold yields a distinct thermometric advantage.","keywords":["quantum thermometry","quantum Fisher information","multilevel quantum Rabi model","adiabatic approximation","bright states","dark states","cavity quantum electrodynamics","random coupling matrices"],"falsifier":"Take a D_g=10, D_e=1010 multilevel quantum Rabi model with random intraband detunings in the dark manifold, compute the exact thermal QFI by numerical diagonalization at high photon cutoff and intermediate coupling, and compare with Eq. (30). If the bright–dark peak is suppressed, shifted, or broadened relative to the analytic formula, the dark-sector decoupling assumption fails.","tokens_in":1300,"feed_emoji":"🌡️","tokens_out":1348,"duration_ms":62657,"temperature":0.7,"pith_summary":"This paper claims that the temperature sensitivity of a multilevel quantum Rabi model—two nearly degenerate atomic manifolds of many levels coupled to a single cavity mode—can be described by an analytic formula in the adiabatic regime, where the atomic frequency is much smaller than the cavity frequency. The formula splits the thermal quantum Fisher information (QFI) into three pieces: intra-doublet bright excitations, inter-doublet bright–bright transitions, and bright–dark population transfer. When the dark manifold is large, the bright–dark piece produces a peak sensitivity that approaches the ultimate limit of an ideal thermometer as the dark degeneracy grows. When instead the bright manifold is large, the response is broadband and becomes self-averaging under random couplings. This matters because it offers a route to designing sensitive or broad-range equilibrium thermometers in cavity quantum electrodynamics without full numerical diagonalization of large systems.","feed_headline":"Dark-manifold saturation nears ideal quantum thermometry","feed_subtitle":"A closed-form expression splits thermal sensitivity into bright and dark channels, guiding large-manifold design.","key_machinery":"The load-bearing construction is the multilevel adiabatic approximation in a superradiant basis: a singular-value decomposition of the coupling matrix separates atomic transitions into bright doublets (collective superpositions that couple strongly to the cavity) and dark states (superpositions that do not couple at all). Assuming a small atomic frequency, each bright doublet becomes a displaced-oscillator ladder with overlap factors set by Laguerre polynomials, while the dark manifold is an equally spaced, highly degenerate ladder; intraband detunings enter as first-order shifts. Combining these spectra in the partition function gives the closed-form QFI, whose log-sum-exp structure separat","core_discovery":"Within the multilevel adiabatic approximation, the thermal QFI of the multilevel quantum Rabi model reduces to a closed-form expression, Eq. (30), built from a partition function separated into bright ladders and a dark manifold. In the dark-manifold-saturated regime, the bright–dark contribution dominates and its peak approaches the ideal-thermometer benchmark as the dark degeneracy grows at fixed bright structure; the approach is robust to disorder in couplings and detunings, and can be improved at intermediate light–matter coupling. In the bright-saturated regime, the bright–bright contribution yields a broadband response whose sample-to-sample fluctuations shrink as the manifold size gro","pith_inferences":["A direct experimental test could tune the degeneracy of an excitonic band in a cavity quantum electrodynamics device and search for the predicted bright–dark peak; the paper does not propose a specific platform.","The same bright/dark decomposition may apply to other bosonic multilevel coupling models with an SVD structure, such as multi-mode or collective-coupling settings, although the paper does not claim this.","Engineering the singular-value spectrum deliberately, rather than drawing couplings randomly, could allow the intermediate-coupling enhancement to be optimized or suppressed; the paper notes that the enhancement is not universal.","A stress test of the main assumption would be to increase intraband disorder inside the dark manifold for large manifolds and check whether the bright–dark peak remains stable; the paper only performs this check for a small system."],"forward_implications":["A probe with a small bright sector and a large dark manifold can reach peak sensitivities close to the fundamental ideal-thermometer bound for the same effective degeneracy and gap.","The peak-QFI ratio relative to the ideal thermometer rises monotonically with dark-manifold degeneracy in both weak and intermediate coupling.","Intermediate light–matter coupling can improve the relative peak QFI for some small bright sectors because spectral reshuffling isolates a single dominant bright–dark gap.","A fully bright manifold with 25 to 250 levels produces a broadband thermometric window that becomes smoother and less sample-dependent as the number of levels grows.","The closed-form QFI and energy-ordered truncation allow ensemble studies over thousands of disordered realisations, which would be prohibitive with brute-force diagonalization."],"fun_headline_variants":["Dark-manifold saturation nears ideal quantum thermometer","Closed-form thermal QFI for multilevel quantum Rabi model","Bright-dard manifolds tune quantum thermometry sensitivity","Large dark manifolds approach ideal sensitivity in thermometry","Quantum thermometer sensitivity peaks at dark-manifold saturation"],"cache_read_input_tokens":22016,"weakest_assumption_plain":"The load-bearing premise is that the small energy splittings inside the dark manifold are thermometrically irrelevant; the paper checks this numerically only for a small system and assumes it for the large-dark-manifold saturation results.","fun_headline_variants_meta":{"raw":{"variants":["Dark-manifold saturation nears ideal quantum thermometer","Closed-form thermal QFI for multilevel quantum Rabi model","Bright-dard manifolds tune quantum thermometry sensitivity","Large dark manifolds approach ideal sensitivity in thermometry","Quantum thermometer sensitivity peaks at dark-manifold saturation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000257,"raw_usage":{"total_tokens":1410,"prompt_tokens":735,"completion_tokens":675,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":597}},"tokens_in":479,"tokens_out":675,"duration_ms":6560,"temperature":1.0,"reasoning_tokens":597,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T23:43:16.200309+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a D_g=10, D_e=1010 multilevel quantum Rabi model with random intraband detunings in the dark manifold, compute the exact thermal QFI by numerical diagonalization at high photon cutoff and intermediate coupling, and compare with Eq. (30). If the bright–dark peak is suppressed, shifted, or broadened relative to the analytic formula, the dark-sector decoupling assumption fails.","supporting_citations":[],"review_version":1}