{"id":"e0fbcbdb-37e5-4c90-972e-6645d3971a77","arxiv_id":"2608.05082","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"The time-evolution operator for any Hamiltonian quadratic in creation and annihilation operators coupled to classical sources factorizes exactly into displacement, squeeze, and rotation operators, whose parameters obey linear differential equations, and the same formalism applies to thermal…","lead":"This paper shows how quantum particles produced by a classical source, such as axions converting to photons in a microwave cavity or gravitons from a black-hole merger, can be described by simple products of displacement, squeezing, and rotation operators. It gives the equations for these operators and extends the framework to thermal states, which matters for designing future axion-search experiments that use squeezed and single-photon-counting detectors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified. The generalized Glauber theorem derivation is sound; minor typos in Eq. (4.11) and a motivational graviton application do not affect the central claim.","rationale":"The paper is a careful, well-structured derivation of a known theorem, with a useful linear formulation of the parameter equations and a thermal-state extension. My independent check of Sections 2–3 found no mathematical error: the Bogoliubov coefficient equations, the symplectic consistency conditions, and the uniqueness argument are all sound. The comparison with the naive Dyson-series evaluation is also correct and illustrates the role of time ordering. The reader's weakest assumption (classical sources, up-to-quadratic Hamiltonian) is indeed the natural boundary of the theorem's applicability, and the paper states it explicitly; this is a scope condition, not a flaw. The only concrete defects I identified are typographical errors in the multi-mode moment formula (4.11), which do not affect the central theorem or the single-mode variance used in the discussion of axion detection. The graviton application is presented as motivation and is supported by citations to the literature, so its brevity does not undermine the central mathematical claim. For these reasons, the reader's ACCEPT verdict remains appropriate.","tokens_in":25660,"tokens_out":27476,"duration_ms":322217,"concrete_test":"Recompute Eq. (4.11) by direct Gaussian P-function integration for a two-mode thermal coherent state with n_1 ≠ n_2, and compare the coefficient of δ_i1 δ_j1 δ_k2 δ_l2 (i.e., ⟨N_1 N_2⟩ cross-correlation) against the published formula; the correct term is n_1 n_2, not n_1^2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation in Section 3.2 is internally consistent: the linear ODEs (3.19) for μ, ν, and α follow from matching the Heisenberg evolution of b'(t) and b(t), the symplectic condition (3.20) preserves constraints (3.13), and uniqueness of the linear initial-value problem gives U'_I(t) = U_I(t) up to a phase. I find no gap in the proof that a quadratic Hamiltonian with c-number sources admits the exact displacement-squeeze-rotation factorization. The reader's weakest assumption concerning classical sources and quadratic interactions is accurate, and this is a stated condition of the theorem, not an unacknowledged limitation. Two minor issues do not rise to load-bearing status: Eq. (4.11) contains subscript and phase typos (the first δ_ij δ_kl term should multiply n_i n_k, and the final phase should be ω_i + ω_k − ω_j − ω_l), though the single-mode variance discussion is unaffected; and the black-hole-merger graviton application is asserted as motivation rather than derived, but it is not the core theorem.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper revisits the generalized Glauber theorem for Hamiltonians that are quadratic in annihilation and creation operators with c-number, time-dependent coefficients. The central result is that the interaction-picture time-evolution operator factorizes exactly into displacement, squeeze, and rotation operators, with parameters satisfying the linear ODEs in Eq. (3.19). The derivation is based on matching the Heisenberg-picture operators b(t) from the exact evolution with the operators b'(t) constructed from the factorized ansatz; uniqueness of the linear initial-value problem then identifies the two evolutions up to a phase. Special cases reproduce the Glauber displacement-only result and the parametric-amplifier squeeze-only result, and a comparison with a naive, non-time-ordered Dyson exponent illustrates the role of time ordering. The formalism is then applied to thermal initial states, yielding explicit first- and second-moment formulas and a thermal-coherent-state probability representation.","tokens_in":25844,"tokens_out":23872,"duration_ms":271816,"significance":"The main theorem is significant because it reduces quantum particle production from classical backgrounds to solving linear c-number equations, with no perturbative approximation. The proof in Section 3.2 is internally consistent: the symplectic condition (3.20) preserves the constraints (3.13), and the same-evolution/same-initial-condition argument gives U'_I = U_I up to a phase. The classical-source and up-to-quadratic restrictions are explicitly stated in Eq. (3.4), so the scope is not overstated; the paper does not claim validity for quantum-fluctuating sources or cubic/higher-order interactions. The thermal-state formulas in Section 4 provide a concrete and useful description of microwave-cavity axion detectors, including nontrivial fourth-order correlations. The numerical convergence check in Appendix C and the explicit special cases are additional strengths. I find no circularity in the derivation, and I agree with the reader that the stated assumptions are genuine hypotheses of the theorem rather than unacknowledged limitations.","major_comments":[],"minor_comments":[{"comment":"In the fourth-order correlator of Eq. (4.11), the first term should read δ_ij δ_kl n_i(T) n_k(T) rather than δ_ij δ_kl n_i(T) n_j(T), and the final phase factor should be exp[i(ω_i+ω_k−ω_j−ω_l)t] rather than exp[i(ω_i+ω_k−2ω_j)t]. As written the expression is not consistent with the pairing rule (4.9); the corrected form is needed if Eq. (4.11) is used for variance calculations.","section":"Section 4.1, Eq. (4.11)"},{"comment":"In Eq. (2.14), the second equality for J3,ij is off by a central term: from the first expression one obtains (1/2)a†_i a_j + (1/4)δ_ij, not (1/2)(a†_i a_j + δ_ij). The commutators in Eq. (2.15) are unchanged because the central contributions cancel, but the displayed identity should be corrected. Also, in the last line of Eq. (2.15) the symbol J± should be J+ in both terms.","section":"Section 2.3, Eqs. (2.14)-(2.15)"},{"comment":"When passing from Eq. (3.1) to Eq. (3.4), the term ℏ Σ_{i,j} h_ij(t) a†_i a_j in the interaction picture should carry a phase e^{i(ω_i−ω_j)t} unless h(t) has been redefined to absorb that phase. The authors should state this redefinition explicitly, since the time dependence of h(t) is part of the input to the linear ODEs (3.19).","section":"Section 3.1, Eqs. (3.1)-(3.4)"},{"comment":"After deriving the linear equations (3.19), the proof would be fully explicit if the authors noted that any solution (μ,ν) satisfying the constraints (3.13) can be written in the form (3.12) for suitable ζ and φ, via the polar decomposition of a symplectic matrix. This standard fact is the step that closes the existence of the factorization U'_I = D S P; adding one sentence would make the argument self-contained.","section":"Section 3.2, after Eq. (3.19)"},{"comment":"There are several typographical issues: 'Analitical' in the caption of Figure 2, 'Galuber' in Section 3.1, 'reaction operators' in Section 3.2, and the left-hand side of the first line of Eq. (3.28) should be β(t) rather than β_i(t) if the sum over i is retained in that expression.","section":"Various"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a good fit for hep-ph; the axion and graviton applications are motivating, and while the graviton part is presented as motivation rather than a full derivation, the generalized Glauber theorem itself is the deliverable. The manuscript is careful and cites the relevant quantum-optics literature. The equation typos in Section 4.1 and Section 2.3 should be fixed before publication, but they do not undermine the central theorem. I recommend minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news in this paper is not the theorem—the generalized Glauber factorization has been in the quantum-optics literature since Schumaker and Ma/Rhodes, and the authors say so themselves. The value is in the presentation and one genuine technical improvement: the parameter evolution is reduced to the linear equations (3.19), which are substantially simpler than the nonlinear equations in the earlier references. That makes the formalism actually usable for particle-physics calculations, and the paper gives the derivations in full. Good.\n\nThe proof in Sec. 3.2 holds up. The factorization ansatz is checked by matching Heisenberg equations for the transformed operators, and uniqueness of the linear initial-value problem gives the equality up to a phase. The symplectic condition on M(t) is consistent with the constraint equations. I went through the key steps and found no circularity or hidden fitting. The thermal-state extension in Sec. 4 is also solid and practically useful: the thermal coherent state with the shifted Gaussian probability density gives a clean way to compute normal-ordered expectation values, and the point about variance mixing thermal and quantum contributions is worth making for axion-experiment simulations. The comparison with the naive Dyson computation and the demonstration that time ordering matters is pedagogically valuable, though it is not new physics.\n\nSoft spots, in proportion: the graviton production from black-hole mergers is motivation, not a derivation. The paper does not provide the black-hole-merger Hamiltonian or show how the quadratic couplings arise; it cites recent work and says the formalism applies. Flagging that clearly would help readers. There is also a typo in Eq. (4.11)—the first δijδkl term should multiply ni nk, and the phase should be ωi + ωk − ωj − ωl. It is cosmetic and does not affect the variance discussion. The paper is long and partly a review, so readers familiar with quantum optics will skim large parts of Sec. 2. That is fine, but it explains the modest novelty score.\n\nThe citation pattern is honest: the prior quantum-optics theorems are credited, the axion and gravitational-wave experiment literature is current, and the recent papers on nonclassical axion dark matter are cited. Nothing looks self-serving.\n\nWho should read this: particle physicists working on axion haloscopes, microwave-cavity quantum optics, or graviton detection who want the quantum-optics machinery in a language they can use directly. It is also useful for people writing Monte Carlo simulation codes for quantum-enhanced axion searches.\n\nYes, this deserves a serious referee. The central derivation is sound, the linear parameter equations are a real simplification, and the thermal-state treatment is a concrete contribution. I would send it to review and expect minor revisions.","headline":"A careful, honest restatement of a known theorem with genuinely simpler linear parameter equations and a useful thermal-state extension; graviton application is asserted, not derived.","tokens_in":26386,"tokens_out":1495,"would_cite":true,"duration_ms":23369,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V80","81R30","22E70"],"pacs":["14.80.Va","04.30.-w","42.50.Dv","03.65.-w"],"model":"deepseek-v4-flash","headline":"For any Hamiltonian quadratic in creation and annihilation operators acting on classical sources, time evolution factors exactly into displacement, squeezing, and rotation operators.","keywords":["generalized Glauber theorem","displacement operator","squeeze operator","rotation operator","axion haloscope","thermal coherent state","graviton production","Dyson series"],"falsifier":"Compute the exact unitary time evolution for a finite set of modes by exponentiating the full quadratic Hamiltonian numerically, for a case with nonzero f, g, and h, and compare expectation values such as the mode occupancy and the pair amplitude against the factorized prediction solved through the linear equations; any mismatch beyond numerical precision would falsify the theorem. Alternatively, in an axion haloscope, measure the photon second-order correlation function across the signal line and check whether it matches the displaced-thermal-coherent prediction including the mixed thermal-quantum terms; a statistically significant deviation would indicate that the classical-source quadratic Hamiltonian is not the right description.","tokens_in":25454,"feed_emoji":"⚛️","tokens_out":6308,"duration_ms":69592,"temperature":0.7,"pith_summary":"This paper establishes a generalized Glauber theorem: when a Hamiltonian contains terms up to quadratic order in creation and annihilation operators, with coefficients that are ordinary complex functions of time (classical sources), the full time-evolution operator factorizes exactly into a displacement operator, a squeezing operator, and a rotation operator, up to a global phase. The paper derives simple linear differential equations for the parameters of these three operators, so the quantum dynamics reduces to solving c-number equations rather than a perturbative Dyson expansion. This matters for experiments because axion-photon conversion in microwave haloscopes emerges as the purely linear (displacement) limit, while squeezed graviton production from black-hole mergers follows from the quadratic terms. The same machinery works for initial thermal states, giving displaced thermal coherent states that describe real microwave cavities and support Monte Carlo simulation of quantum-enhanced axion searches.","feed_headline":"Quadratic interactions factor into displacement, squeezing, rotation","feed_subtitle":"Axion haloscope signals and black-hole merger gravitons are two limits of one exact quantum formula.","key_machinery":"The machinery is the closed Lie algebra generated by $\\hat{a}_i$, $\\hat{a}_i^\\dagger$, and the quadratic operators $\\hat{J}_+$, $\\hat{J}_-$, $\\hat{J}_3$, together with the generalized Baker-Campbell-Hausdorff theorem. From this algebra, the unitary is written as $\\hat{U}'_I(t) = \\hat{D}\\hat{S}\\hat{P}$, and the parameter dynamics are fixed by demanding that the Heisenberg-picture operators $\\hat{b}'(t) = \\hat{U}'^\\dagger \\hat{a} \\hat{U}'$ satisfy the same linear equation as the exact $\\hat{b}(t)$. The resulting system $(d/dt)[\\mu^*; -\\nu^*] = M [\\mu^*; -\\nu^*]$ and $(d/dt)(\\alpha; \\alpha^*) = M(\\alpha; \\alpha^*) + (f; f^*)$, with $M = [[-ih, -g], [-g^*, ih^*]]$, is linear and solved with initial conditions $\\mu(t_0)=I$, $\\nu(t_0)=0$, $\\alpha(t_0)=0$.","core_discovery":"On the paper's own terms, the central claim is that for the interaction-picture Hamiltonian $$H_I(t) = i\\hbar \\sum_i [f_i(t)\\hat{a}_i^\\dagger - f_i^*(t)\\hat{a}_i] + \\frac{i}{2}\\hbar \\sum_{i,j}[g^*_{ij}(t)\\hat{a}_i \\hat{a}_j - g_{ij}(t)\\hat{a}_i^\\dagger \\hat{a}_j^\\dagger] + \\hbar \\sum_{i,j} h_{ij}(t)\\hat{a}_i^\\dagger \\hat{a}_j$$ with $f,g,h$ c-number functions, the interaction-picture evolution operator equals $\\hat{D}(\\{\\alpha(t)\\})\\hat{S}(\\{\\zeta(t)\\})\\hat{P}(\\{\\phi(t)\\})$ up to an irrelevant phase. The proof proceeds by writing the Heisenberg equations for the transformed operators and matching them to the same linear system $(d/dt)(\\hat{b},\\hat{b}^\\dagger)^T = M(t)(\\hat{b},\\hat{b}^\\dagger)^T + (f,f^*)^T$, giving linear equations for the parameters. The paper recovers the original Glauber coherent state as the special case with only $f(t)$ nonzero, recovers squeezed vacuum for real diagonal $g(t)$, and shows that a naive 'integrate the exponent first' application of the Dyson series fails when time ordering matters, with deviations shown numerically.","pith_inferences":["If the theorem is right, the same factorization should apply to any bosonic system with classical drivings and quadratic couplings, e.g. phonons in optomechanics or magnons in cavity spintronics, wherever the source stays classical.","A testable consequence is that the second-order correlation function $g^{(2)}$ of microwave photons in an axion haloscope should interpolate between thermal statistics and coherent statistics according to the displaced-thermal formula; measuring $g^{(2)}$ at signal-on and signal-off would probe the cross terms directly.","The classical probability density found for thermal coherent states suggests that existing axion-search Monte Carlo codes could be upgraded by sampling complex amplitudes from that distribution instead of heuristic noise power estimates; the paper does not run such a simulation itself.","For black-hole merger gravitons, the formalism implies that the quantum state's squeezing parameters are computable from the classical metric perturbation alone; whether such squeezed graviton states are detectable remains open, but their statistics are now defined."],"forward_implications":["Axion-photon conversion in haloscopes is exactly the f-only limit, so the emitted microwave field is a (thermal) coherent state whose displacement is set by the axion coupling.","Graviton production from a classical black-hole merger includes squeezing and mode mixing, so the produced gravitons are in a displaced squeezed state rather than a coherent state.","For thermal initial states, the expectation values factor into thermal and coherent parts, but the variance contains mixed thermal-times-quantum cross terms, so noise estimates cannot be obtained by adding the two contributions independently.","The linear differential equations for the parameters allow the full quantum state to be computed without truncating the Dyson series, so the formalism is a practical tool for the parameter regimes of realistic experiments.","A classical probability density over coherent-state amplitudes is provided, giving a direct sampling rule for Monte Carlo simulations of axion detectors."],"supporting_citations":[{"why":"Establishes the original Glauber theorem that linear coupling to a classical current produces coherent states, the special case this paper generalizes.","marker":"[38]"},{"why":"Shows that quadratic Hamiltonians lead to Gaussian states and introduces the displacement-squeeze-rotation decomposition that the paper re-derives with linear parameter equations.","marker":"[43]"},{"why":"Provides the multimode squeeze-operator formalism and squeezed states that the generalized theorem builds on.","marker":"[44]"},{"why":"Pursues a similar squeezed-state description for gravitons, which this paper extends to a general and comprehensive treatment.","marker":"[45]"},{"why":"Introduces the inverse Gertsenshtein effect by which gravitational waves convert to photons, the process behind the graviton application.","marker":"[37]"},{"why":"Documents the HAYSTAC demonstration of squeezing in an axion haloscope, the experimental context that motivates the thermal-coherent-state treatment.","marker":"[18]"},{"why":"Gives a recent coherent-state description of gravitational waves from binary black holes, the application this paper places under the generalized Glauber theorem.","marker":"[42]"}],"fun_headline_variants":["Generalized Glauber theorem unifies axion and graviton detection","Exact factorization unifies axion haloscope and graviton emissions","Quadratic quantum dynamics reduced to c-number equations","Generalized Glauber theorem: one formula for axion and graviton"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the interaction Hamiltonian is exactly quadratic in creation and annihilation operators with no higher-order terms and with all coefficients being ordinary complex numbers (classical sources), so that a single closed algebra governs the evolution; if the source has quantum fluctuations or the Hamiltonian has cubic or higher terms, the exact factorization into displacement, squeezing, and rotation no longer holds.","fun_headline_variants_meta":{"raw":{"variants":["Generalized Glauber theorem unifies axion and graviton detection","Exact factorization unifies axion haloscope and graviton emissions","Quadratic quantum dynamics reduced to c-number equations","Generalized Glauber theorem: one formula for axion and graviton"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000629,"raw_usage":{"total_tokens":2953,"prompt_tokens":1038,"completion_tokens":1915,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":1841}},"tokens_in":654,"tokens_out":1915,"duration_ms":16193,"temperature":1.0,"reasoning_tokens":1841,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T05:41:52.355780+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact unitary time evolution for a finite set of modes by exponentiating the full quadratic Hamiltonian numerically, for a case with nonzero f, g, and h, and compare expectation values such as the mode occupancy and the pair amplitude against the factorized prediction solved through the linear equations; any mismatch beyond numerical precision would falsify the theorem. Alternatively, in an axion haloscope, measure the photon second-order correlation function across the signal line and check whether it matches the displaced-thermal-coherent prediction including the mixed thermal-quantum terms; a statistically significant deviation would indicate that the classical-source quadratic Hamiltonian is not the right description.","supporting_citations":[{"cited_title":"Schumaker,Quantum mechanical pure states with gaussian wave functions,Phys","cited_arxiv_id":null,"evidence_quote":"Shows that quadratic Hamiltonians lead to Gaussian states and introduces the displacement-squeeze-rotation decomposition that the paper re-derives with linear parameter equations."},{"cited_title":"Ma and W","cited_arxiv_id":null,"evidence_quote":"Provides the multimode squeeze-operator formalism and squeezed states that the generalized theorem builds on."},{"cited_title":"Gertsenshtein,Wave resonance of light and gravitational waves,Zhurnal Eksperimental’noi i Teoreticheskoi Fiziki41(1961) 113","cited_arxiv_id":null,"evidence_quote":"Introduces the inverse Gertsenshtein effect by which gravitational waves convert to photons, the process behind the graviton application."},{"cited_title":"Kanno, J","cited_arxiv_id":null,"evidence_quote":"Gives a recent coherent-state description of gravitational waves from binary black holes, the application this paper places under the generalized Glauber theorem."}],"review_version":1}