{"id":"3f8ea908-c6f0-4c3e-8f29-8c7d1edc7bc1","arxiv_id":"2504.19164","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Adiabatic number operators and vacua for coupled bosonic systems are constructed order by order from a linear recursion for complex structures, generalizing WKB and Lewis-Riesenfeld invariant methods beyond a single mode.","lead":"The paper builds adiabatic vacuum states for quantum systems with many coupled oscillators and time-dependent couplings, using a geometric tool called a linear complex structure. It generalizes the standard WKB trick from one oscillator to many, with applications to cosmology, many-body physics, and quantum thermodynamics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Definition 1 leaves the remainder R(bar n) unspecified; different choices satisfying J0^2 = -1 yield different finite-order adiabatic vacua (Fig. 5), so the central claim is underdetermined without a canonical remainder prescription.","rationale":"The paper's core algebraic recursion (38)-(43) is coherent: Proposition 1's uniqueness argument for the coefficients Jn and zeta_n is sound, and the d=1 reduction reproduces the WKB adiabatic vacuum (Proposition 2) with the expected O(epsilon^{bar n+1}) discrepancy (Examples 5-6). The external benchmarks — Airy asymptotic match, Bunch-Davies power spectrum, Lewis-Riesenfeld comparison — support the physical relevance of the underlying idea. However, the central construction as stated is underdetermined. Definition 1 introduces R(bar n) with the single condition J0^2 = -1, and the paper provides no existence proof, no order bound, and no canonical choice. The paper's own Example 3 and Figure 5 demonstrate that different R(2) choices yield different J0 and different particle number at order bar n+1. This is load-bearing because the central claim asserts that the order-bar n adiabatic number operator and vacuum are expressible in terms of the order-bar n adiabatic complex structure; if that structure is ambiguous, the state and all derived predictions are ambiguous. This is not an external disagreement with consensus; it is an internal underdetermination of the algorithm. The concrete test above would establish whether the ambiguity is benign (confined to O(epsilon^{bar n+1}), in which case a canonical prescription such as a projection onto the Kaehler-structure manifold would resolve it without changing the advertised accuracy) or malignant (affecting the leading-order prediction). Given that the needed fix is a specification rather than a rejection of the recursion, the appropriate verdict remains CONDITIONAL rather than REJECT.","tokens_in":28718,"tokens_out":17830,"duration_ms":174291,"concrete_test":"Evaluate order bar n=2 adiabatic vacua for a one-parameter family of slow protocols (e.g., omega(t)=omega0(1+epsilon t) with epsilon -> 0) using two remainder prescriptions: (i) the diagonal ansatz of Eqs. (79)-(80), and (ii) the WKB complex structure of Eq. (114) with W0, dotW0 from Eqs. (104)-(105). Compute the particle number (177) at a fixed t1 and the norm difference ||J0^(i)-J0^(ii)||_g. If the difference and the particle-number discrepancy scale as epsilon^3 (i.e., epsilon^{bar n+1}), the remainder ambiguity is confined to the advertised truncation error and a canonical prescription could fix Definition 1; if they scale as epsilon^2 or slower, the order-2 adiabatic vacuum is not well-defined by the recursion. For a d>1 cross-check, repeat with the magnetic-field Hamiltonian of Sec. V.D and verify J0^T Omega J0 = Omega and G>0 for each prescription.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The construction's central output, the adiabatic complex structure J0 of order bar n, is defined in Definition 1 (Eq. 51) as the sum of the uniquely determined Jn(t0) plus a remainder R(bar n) that enforces J0^2 = -1. The recursion (38)-(43) of Proposition 1 determines the Jn uniquely, but it does not determine R(bar n): the equation (S+R)^2 = -1 for S = sum_{n=0}^{bar n} Jn(t0) has multiple real solutions in any neighborhood of the instantaneous structure J0. The paper itself acknowledges this ('The solution is not unique', Eqs. 79-80) and Figure 5 shows two choices of R(2) — the quadratic ansatz (79)-(80) and the WKB prescription — producing different particle-number predictions at the same order. Since the central claim is that the adiabatic number operator and vacuum of order bar n are expressed in terms of the adiabatic complex structure of the same order, non-uniqueness of that complex structure makes the order-bar n state, and any physical prediction derived from it (e.g., Eq. 60, Eq. 86), ambiguous. The paper neither proves that all admissible remainders differ only at order epsilon^{bar n+1} in a norm sense nor supplies a canonical prescription (e.g., symplectic polar decomposition of S, or projection onto the Kaehler-structure manifold). This is not a convergence issue: even in the formal/asymptotic regime, Definition 1 is incomplete. The Discussion's claim of 'a well-defined framework with a unique solution' overstates Proposition 1, which proves uniqueness of the Jn, not of J0.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a construction of adiabatic vacua for finite-dimensional systems with time-dependent quadratic Hamiltonians, using linear complex structures on phase space. Starting from the invariant equation for a Gaussian number operator, the authors expand the complex structure J(t) and the displacement z(t) in formal power series in a slow-time parameter λ and derive a first-order linear recursion, Eqs. (38)-(43), for the coefficients J_n and ζ_n. Truncating at order n and adding a remainder R(n) that enforces J0^2 = -1 defines the adiabatic initial condition (J0, z0), from which the adiabatic number operator and Gaussian vacuum are defined. The paper compares the construction to the WKB method for a single oscillator, proving equivalence up to order n+1 (Proposition 2), and tests the formalism on several examples: a driven oscillator, an Airy frequency profile, a tanh profile with exact Bogoliubov coefficients, a charged particle in a time-dependent magnetic field, and cosmological scalar perturbations. In de Sitter space the series is claimed to self-truncate at order 2 and to reproduce the Bunch-Davies vacuum.","tokens_in":29134,"tokens_out":21993,"duration_ms":198518,"significance":"If the construction is made fully precise, the recursion is a valuable contribution: it generalizes WKB and Lewis-Riesenfeld invariants to d coupled bosonic modes through linear first-order equations, and it is accompanied by explicit checks against exact solutions (Airy, de Sitter, tanh Bogoliubov). The proof of Proposition 1 is careful, and the comparison with WKB in Proposition 2 is explicit and informative. The main limitation is that Definition 1 leaves the remainder R(n) unspecified, so the finite-order adiabatic complex structure, and hence the finite-order adiabatic vacuum, is not uniquely determined by the recursion as it stands; this does not invalidate the recursion itself but does affect the paper's central claim of a well-defined unique finite-order state.","major_comments":[{"comment":"The recursion (38)-(43) determines the coefficients J_n and ζ_n uniquely, but Definition 1 adds a remainder R(n) that is not determined by the recursion. The paper acknowledges the non-uniqueness ('The solution is not unique', Eq. 79-80) and Figure 5 shows two admissible choices of R(2) giving different particle-number predictions at the same order n=2. Since J0 is used in Definition 2 and Definition 3 to define the adiabatic number operator and the adiabatic vacuum, the finite-order state is underdetermined. The manuscript should either give a canonical prescription for R(n) (for example, a symplectic or polar projection of the truncated sum onto the manifold of complex structures) or prove that all admissible remainders produce the same physical predictions at order n.","section":"Sec. II.D, Definition 1 (Eq. 51); Example 3 (Eqs. 79-80); Fig. 5"},{"comment":"The assertion R(n) = O(epsilon^{n+1}) is made without proof. While the recursion enforces J^2 = -1 order by order, the existence of an exact complex structure within O(epsilon^{n+1}) of the truncated sum S = sum_{m=0}^{n} J_m(t0) is a nontrivial statement, especially for d > 1. The explicit check in Example 3 covers only a quadratic ansatz in d = 1. This gap is load-bearing because the claimed smallness of the remainder is what justifies ignoring the ambiguity in the truncation. Please add a general existence argument, or state precisely the smoothness and positivity conditions under which such a remainder exists.","section":"Sec. II.D, Definition 1 (Eq. 51); Sec. VI"},{"comment":"The Discussion states that the construction yields 'a well-defined framework with a unique solution.' This statement overstates Proposition 1, which proves uniqueness only for the formal coefficients J_n and ζ_n, not for the adiabatic initial condition J0 of Definition 1, because of the non-unique remainder. The manuscript should either qualify this claim or supply the missing canonical remainder prescription so that the uniqueness claim becomes accurate.","section":"Sec. VI (Discussion)"}],"minor_comments":[{"comment":"'reminder' should be 'remainder'.","section":"Example 3 (after Eq. 77)"},{"comment":"The typos 'appoach' and 'spatimes' should be corrected to 'approach' and 'spacetimes'.","section":"Abstract and Section I"},{"comment":"The caption should state explicitly that the black curve uses the WKB-based remainder and the orange curve uses the ansatz of Eq. (79).","section":"Figure 5 caption"},{"comment":"The notation W0 is used both for the order-zero WKB frequency and for the infinite-order sum in Eq. (164); please use distinct notation, for example W_infty.","section":"Section V.B (Eq. 164)"},{"comment":"The proof would be easier to follow if it showed explicitly that the solution of Eq. (48) indeed satisfies both Eqs. (41) and (42); currently the reversibility of the derivation is implicit.","section":"Proof of Proposition 1 (Eq. 48)"},{"comment":"The claim that the adiabatic series self-truncates at order n=2 in de Sitter space is stated but not demonstrated; a short argument, for example showing J_n = 0 for n >= 3, would strengthen the example.","section":"Example 7 (Eq. 150)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the recursion in Proposition 1 is a solid and useful result. The main obstacle to acceptance is the unspecified remainder in Definition 1, which directly affects the uniqueness claim in the abstract and Discussion. This is fixable within the scope of the manuscript. The reliance on [17] for Gaussian-state machinery is legitimate reuse of earlier work by two of the authors."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this paper generalizes WKB/Lewis-Riesenfeld adiabatic vacua to d coupled bosonic modes via a linear recursion for complex structures, and that generalization is real. The main soft spot is Definition 1: the remainder R(n-bar) that enforces J^2 = -1 is not specified, and the paper itself admits different choices give different finite-order vacua (Figure 5). That gap is fixable, but the Discussion's claim of a 'unique solution' overstates what Proposition 1 actually proves.\n\nWhat's new: Eqs. (38)-(43) give a first-order linear recursion for the coefficients J_n and zeta_n. Previous WKB and Lewis-Riesenfeld methods handle one mode; the multimode case with time-dependent couplings had no systematic adiabatic construction of this kind. The authors check against exact solutions (Airy, de Sitter, tanh profile) and against WKB, and the checks line up. That external benchmarking is real evidence the recursion is correct.\n\nWhere it's soft: the remainder R(n-bar) in Definition 1 is non-unique, and the paper acknowledges this in Example 3 and Figure 5. Two reasonable choices give different particle-number predictions at the same order. So the advertised claim—adiabatic number operator and vacuum of order n are determined by the adiabatic complex structure of order n—is underdetermined until you specify a canonical prescription (symplectic polar decomposition, projection, or something). The paper asserts R = O(epsilon^{n+1}) without a proof, and the Bunch-Davies identification in Example 7 is stated rather than demonstrated. These are real gaps, but they are localized; they don't sink the recursion itself. Proposition 1's uniqueness is only for the J_n, not for the truncated J_0.\n\nWho benefits: people working on adiabatic regularization in cosmology, many-body systems with time-dependent Hamiltonians, and shortcuts to adiabaticity. It deserves a serious referee. I'd ask the referee to focus on Definition 1 and whether a canonical remainder choice can be proved to preserve the claimed order. If the authors fix that, this is a solid contribution.\n\nRecommendation: accept for peer review.","headline":"Genuine multimode generalization of adiabatic vacua with a real but fixable ambiguity in the remainder.","tokens_in":29601,"tokens_out":1780,"would_cite":true,"duration_ms":17339,"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":"The paper shows that finite-order adiabatic vacua of any quadratic time-dependent Hamiltonian can be built from a linear recursion on phase-space complex structures, generalizing WKB and Lewis-Riesenfeld invariants to d coupled bosonic…","keywords":["adiabatic vacuum","linear complex structure","quadratic time-dependent Hamiltonian","WKB approximation","Lewis-Riesenfeld invariant","Gaussian states","particle production","adiabatic subtraction"],"falsifier":"Compute the adiabatic vacuum of order 2 for the oscillator with time-dependent frequency using each of the two remainder choices described in the paper and evaluate the expectation value of an order-3 number operator; if the answers differ at order $ε^{3}$, then order-(n+1) predictions are not determined by the recursion alone. A more direct test is to search for a smooth frequency profile where the quadratic remainder equation has no solution with the required smallness, which would break Definition 1.","tokens_in":28522,"feed_emoji":"🌀","tokens_out":5913,"duration_ms":55184,"temperature":0.7,"pith_summary":"The paper claims that for any quadratic time-dependent Hamiltonian with d coupled bosonic degrees of freedom, the adiabatic vacuum and adiabatic number operator of a given order can be constructed directly from a linear complex structure, a phase-space matrix that squares to minus one and encodes a Gaussian state. The construction is a recursion that produces, order by order, a complex structure J0 and a displacement vector z0 from time derivatives of the Hamiltonian alone. Because the equations are linear and first order in time, it avoids the nonlinear WKB equation needed for a single oscillator. The authors prove the recursion has a unique solution at each order, show it reproduces WKB and Lewis-Riesenfeld results when those apply, and exhibit examples in cosmology, driven oscillators, and a charged particle in a time-dependent magnetic field.","feed_headline":"Complex structures turn adiabatic vacua into a linear recursion","feed_subtitle":"One recursion handles coupled bosonic modes, generalizing WKB and exact invariants to many degrees of freedom.","key_machinery":"The central object is the linear complex structure J, a 2d×2d real matrix with $J^{2}$=-1 that, together with the symplectic form Ω, defines a positive metric G=-JΩ and a Gaussian vacuum via the number operator N_{J,z}=12(ξ-z)·$Ω^{{-1}}$·(J-i1)·(ξ-z). The mechanism is the slow-time reparametrization t→(t-t0)/λ+t0: analyticity in λ of the invariant number operator turns the dynamical equations λ Jdot=[K,J] and λ zdot=Kz+F into the order-by-order recursion (38)-(43). Because the equations are linear and first order in time, the recursion only needs derivatives of K and F and a single matrix anti-commutator inversion at each step; no nonlinear equation of Ermakov type must be solved.","core_discovery":"With a time-dependent quadratic Hamiltonian written as H(t)=12ξ·h(t)·ξ+f(t)·ξ+c(t), any Gaussian state is fixed by a pair (J,z) with $J^{2}$=-1. The central claim is that imposing invariance of the associated number operator under the reparametrized dynamics, order by order in a slow-time parameter λ, reduces to linear algebraic equations (38)-(43) for the Taylor coefficients J_n and ζ_n. Proposition 1 states these equations have a unique solution at every order whenever K=Ωh and F=Ωf are sufficiently differentiable; the zeroth-order terms are J0=|$K^{{-1}}$|K and ζ0=-$K^{{-1}}$F, and the higher orders are built by inverting an anti-commutator. The adiabatic initial conditions at a reference time are then the truncated sums (51)-(52), with a remainder R(n) that enforces $J0^{2}$=-1 exactly. The resulting number operator and its ground state are the paper's finite-order adiabatic number operator and adiabatic vacuum.","pith_inferences":["The freedom in the remainder R(n) deserves to be treated as a completion rule: Example 3 and Figure 5 show two different remainder choices give different particle-number curves, so any prediction at order n+1 depends on a choice the recursion alone does not fix.","If combined with optimal-truncation ideas, the complex-structure formulation could turn single-mode results on superadiabatic particle number and Stokes phenomena into a multi-mode algorithm; the paper mentions this as future work but does not carry it out.","Because the construction is phrased entirely in phase space, it should transfer to fermionic quadratic Hamiltonians through the same Kähler-structure language; the paper notes this possibility.","For quantum fields, the formalism may provide a route to define adiabatic vacua directly from the background's real-time history rather than through mode-by-mode WKB, potentially simplifying renormalization in inhomogeneous spacetimes."],"forward_implications":["For any Hamiltonian in the class, the adiabatic vacuum of order n can be computed by an explicit linear recursion, so the method applies to coupled multi-mode systems where WKB on individual Fourier modes fails.","When restricted to a single oscillator with time-dependent frequency, the resulting state agrees with the WKB adiabatic vacuum up to order n+1 (Proposition 2), tying the new method to the established literature.","The same construction defines adiabatic subtraction for the energy: the renormalized Hamiltonian (156) has exactly zero expectation value in the adiabatic vacuum of the same order, matching the order n≥4 needed in curved-spacetime renormalization.","For special frequency profiles the expansion self-truncates and produces exact adiabatic states of infinite order; examples include de Sitter's Bunch-Davies vacuum and an oscillator with ω(t)=ω0/(1-t/2τ)^2.","In general the adiabatic vacuum of infinite order does not evolve into itself under exact unitary evolution; the tanh-frequency example shows residual particle production even at infinite order."],"supporting_citations":[{"why":"Supplies the Kähler-structure formalism for Gaussian states, including the definition of the number operator N_{J,z} and vacuum correlations used throughout.","marker":"[17]"},{"why":"Introduces the exact Lewis-Riesenfeld invariant for a time-dependent harmonic oscillator and a charged particle in a time-dependent electromagnetic field, which the paper generalizes.","marker":"[18]"},{"why":"Establishes the WKB adiabatic initial conditions as preferred instantaneous vacuum in cosmology, the benchmark the new construction extends.","marker":"[9]"},{"why":"Original treatment of particle creation in expanding universes and the notion of adiabatic vacuum for cosmological modes.","marker":"[10]"},{"why":"Provides the standard definition of adiabatic particle number and vacuum for a single oscillator, used for comparison.","marker":"[12]"},{"why":"Supplies the adiabatic theorem background that motivates defining initial states which follow a slowly-changing Hamiltonian.","marker":"[13]"},{"why":"Williamson's theorem on symplectic normal forms is used to show that J0=|K^{-1}|K is a complex structure with paired imaginary eigenvalues.","marker":"[33]"},{"why":"Defines the Bunch-Davies vacuum in de Sitter space, which the paper identifies as the infinite-order adiabatic vacuum in that example.","marker":"[4]"},{"why":"Discusses superadiabatic particle number and optimal truncation order, which the paper cites for the asymptotic nature of its expansion.","marker":"[29]"}],"fun_headline_variants":["Linear complex structures unlock adiabatic vacua for many modes","Adiabatic vacua via linear complex structures: generalized","Generalizing WKB and invariants: adiabatic vacua from complex structures","Linear recursion for adiabatic vacua in coupled bosonic modes","One method for many coupled modes: adiabatic vacua via complex structures"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction assumes that a remainder term R(n) always exists that restores the exact condition $J0^{2}$=-1 while being of order n+1; the paper shows one way to choose it in a single example but does not prove existence or give a canonical choice for all Hamiltonians.","fun_headline_variants_meta":{"raw":{"variants":["Linear complex structures unlock adiabatic vacua for many modes","Adiabatic vacua via linear complex structures: generalized","Generalizing WKB and invariants: adiabatic vacua from complex structures","Linear recursion for adiabatic vacua in coupled bosonic modes","One method for many coupled modes: adiabatic vacua via complex structures"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000849,"raw_usage":{"total_tokens":3677,"prompt_tokens":913,"completion_tokens":2764,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":529,"completion_tokens_details":{"reasoning_tokens":2675}},"tokens_in":529,"tokens_out":2764,"duration_ms":16163,"temperature":1.0,"reasoning_tokens":2675,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T06:00:18.585852+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the adiabatic vacuum of order 2 for the oscillator with time-dependent frequency using each of the two remainder choices described in the paper and evaluate the expectation value of an order-3 number operator; if the answers differ at order $ε^{3}$, then order-(n+1) predictions are not determined by the recursion alone. A more direct test is to search for a smooth frequency profile where the quadratic remainder equation has no solution with the required smallness, which would break Definition 1.","supporting_citations":[{"cited_title":"White, Asymptotic Analysis of Differential Equations (Imperial College Press, 2010)","cited_arxiv_id":null,"evidence_quote":"Supplies the Kähler-structure formalism for Gaussian states, including the definition of the number operator N_{J,z} and vacuum correlations used throughout."},{"cited_title":"Cosmological particle production and the precision of the WKB approximation","cited_arxiv_id":"gr-qc/0510001","evidence_quote":"Introduces the exact Lewis-Riesenfeld invariant for a time-dependent harmonic oscillator and a charged particle in a time-dependent electromagnetic field, which the paper generalizes."},{"cited_title":"Parker and S","cited_arxiv_id":null,"evidence_quote":"Establishes the WKB adiabatic initial conditions as preferred instantaneous vacuum in cosmology, the benchmark the new construction extends."},{"cited_title":"Birrell, The application of adiabatic regularization to calculations of cosmological interest, Proc","cited_arxiv_id":null,"evidence_quote":"Original treatment of particle creation in expanding universes and the notion of adiabatic vacuum for cosmological modes."},{"cited_title":"Parker, Quantized fields and particle creation in ex- panding universes","cited_arxiv_id":null,"evidence_quote":"Supplies the adiabatic theorem background that motivates defining initial states which follow a slowly-changing Hamiltonian."},{"cited_title":"de Gosson, Symplectic Geometry and Quantum Me- chanics (Springer, 2006)","cited_arxiv_id":null,"evidence_quote":"Williamson's theorem on symplectic normal forms is used to show that J0=|K^{-1}|K is a complex structure with paired imaginary eigenvalues."},{"cited_title":"Fulling, Aspects of Quantum Field Theory in Curved Space-time (Cambridge University Press, 1989)","cited_arxiv_id":null,"evidence_quote":"Defines the Bunch-Davies vacuum in de Sitter space, which the paper identifies as the infinite-order adiabatic vacuum in that example."}],"review_version":1}