A Bogoliubov-ratio framework for quantum-information diagnostics of time-dependent two-mode Boson System

AI summary

Overview

Research area: Quantum information and quantum field theory — specifically two-mode squeezing, Bogoliubov transformations, and the reduced-state (purity/entropy) diagnostics of time-dependent bosonic pair-production systems, with applications to cosmological perturbations and nonlinear/parametric optics.

Technical level: Advanced. The paper assumes familiarity with Bogoliubov coefficients, SU(1,1) squeezing algebra, Gaussian covariance-matrix methods, and von Neumann/Rényi entropies, though its central object — a single complex ratio — is conceptually simple.

Scope (one sentence): The paper introduces the Bogoliubov ratio λ_k = β_k/α_k as a single closed dynamical variable that governs both the evolution and the one-mode reduced-state purity and entropies of time-dependent two-mode bosonic pair states, and extends it exactly to non-Gaussian SU(1,1) Fock sectors.

What This Paper Is About

When a two-mode system undergoes time-dependent quadratic evolution, excitations are produced in correlated pairs (momenta k and −k), producing a two-mode squeezed state. Describing how much information is lost when an observer can access only one mode of the pair conventionally requires tracking two real squeezing variables — an amplitude r_k and a phase φ_k — which obey coupled nonlinear equations and can oscillate rapidly, and then rebuilding the reduced density matrix separately for every model.

This paper asks whether a single complex number, the ratio of the two Bogoliubov coefficients λ_k = β_k/α_k, can carry all of that model-dependent dynamics and simultaneously determine the reduced-state spectrum. The goal is a compact "two-stage" calculational scheme in which model physics enters once, through λ_k, and universal quantum-information quantities follow automatically.

Key Contributions

1. A closed dynamical equation for the ratio. For a general Hermitian quadratic two-mode Hamiltonian with effective frequency ω_k(η) and pair-production coupling g_k(η), the authors show that λ_k obeys a single complex Riccati-type equation, dλ_k/dη = g_k(1 − λ_k²) − 2iω_k λ_k, replacing the coupled real equations for squeezing amplitude and phase.

2. A universal map from dynamics to information measures. After tracing out the −k mode, the reduced spectrum is the geometric distribution p_n = (1 − q_k) q_k^n with q_k = |λ_k|². Purity, linear entropy, Rényi-2 entropy, and von Neumann entropy are all written directly in terms of |λ_k|², giving the chain {ω_k, g_k} → λ_k → q_k → {μ_k, S_L,k, S_2,k, S_vN,k}. The authors state explicitly that this is an exact re-parametrization of the Gaussian pair sector, mathematically equivalent to the density-matrix and covariance-matrix formalisms, not a new set of information measures.

3. A radial/phase decomposition with a clear physical reading. Writing λ_k = √q_k e^{iθ_k} separates pair production (which drives radial growth of q_k through cos θ_k) from frequency rotation and detuning (which act through the phase θ_k and can suppress coherent squeezing accumulation).

4. An exact non-Gaussian SU(1,1) extension. For a lowest-weight Fock seed with conserved number difference d, the same ratio equation governs the evolution, while the reduced spectrum becomes a negative-binomial distribution determined by the pair (q_k, d). The Gaussian results are recovered at d = 0.

Main Findings

• One variable suffices for the Gaussian sector: λ_k combines the information normally carried by the two real squeezing variables r_k and φ_k, and its modulus squared alone fixes the reduced-state spectrum and all four information measures.

• Explicit entropy formulas in terms of |λ_k|²: purity μ_k = (1 − |λ_k|²)/(1 + |λ_k|²); linear entropy S_L,k = 2|λ_k|²/(1 + |λ_k|²); Rényi-2 entropy S_2,k = ln[(1 + |λ_k|²)/(1 − |λ_k|²)]; and von Neumann entropy S_vN,k = −ln[1 − |λ_k|²] − (|λ_k|²/(1 − |λ_k|²)) ln|λ_k|². Vacuum (λ_k = 0) gives unit purity and vanishing entropies.

• The measures are not independent: S_L,k = 1 − μ_k and S_2,k = −ln μ_k, and S_vN,k is likewise a function of the same q_k.

• Consistency with standard Gaussian tools: the mean occupation is n̄_k = q_k/(1 − q_k), the one-mode symplectic eigenvalue is ν_k = (1 + q_k)/[2(1 − q_k)], and μ_k = 1/(2ν_k), S_2,k = ln(2ν_k), while S_vN,k takes the standard form (n̄_k + 1)ln(n̄_k + 1) − n̄_k ln n̄_k — showing the Bogoliubov-ratio, reduced-density-matrix, and covariance-matrix descriptions contain the same information.

• Division of labor between coupling and frequency: the pair-production coupling g_k supplies the radial drive, whereas ω_k rotates the pair phase and influences q_k only indirectly, through the phase alignment entering cos θ_k. Frequency rotation and detuning can therefore suppress coherent squeezing accumulation.

• Exact state-level Bogoliubov representation of the deformed state: the normalized Meixner-polynomial two-mode state is shown to equal √(1 − |λ_k|²) Σ_n λ_k^n |n_k, n_−k⟩ with λ_k = −√|1 − u₁²| e^{2iφ_k} tanh r_k / (1 + u₂ tanh r_k), and the undeformed limit u₁ = u₂ = 0 reproduces the canonical coefficients α_k^(0) = cosh r_k, β_k^(0) = −e^{2iφ_k} sinh r_k.

• Krylov/Meixner dictionary: the generalized Lanczos coefficients b_n² = |1 − u₁²| n(n − 1 + χ) and c_n = iu₂(2n + χ) map onto the Hamiltonian parameters, with χ = d + 1; in the vacuum pair sector d = 0 one has b_n = |g_k| n, −ic_n = ω_k(2n + 1), hence u₂ = ω_k and |1 − u₁²| = |g_k|² (or u₂ = ω_k/Ω₀, |1 − u₁²| = |g_k/Ω₀|² under normalization by a reference frequency Ω₀).

• Cosmological realization: for a modified dispersion function f(η) and scale factor a(η), u₂ = (k/2)[f²(η) + 1] = ω_k(η) and |1 − u₁²| = (k²/4)[f²(η) − 1]² + [a′(η)/a(η)]² = |g_k(η)|², so both deformation parameters generally evolve with the background rather than being fixed constants.

• Optical consistency check: in the resonant constant-coupling limit the ratio equation gives λ_opt = tanh τ, reproducing the standard two-mode squeezing solution. For a finite Gaussian pump with time-dependent detuning, correlation production is localized near the pulse center, the chirp suppresses subsequent accumulation, and the reduced-state quantities show a delayed response approaching finite late-time values.

• Domain of validity is stated explicitly: reduced-state diagnostics are evaluated only within the normalized Bogoliubov domain |λ_k| < 1, i.e. |α_k|² − |β_k|² = 1.

• What is not reported: the provided text contains no numerical values, benchmark comparisons, tabulated results, or quantitative plots for the cosmological and optical cases; those are deferred to Sec. VI and the appendices, which are not included here.

Methodology in Plain English

The authors start from the most general quadratic two-mode Hamiltonian one can write for a pair of opposite-momentum modes: a term that rotates the phase of each mode at a time-dependent frequency ω_k, plus a term that creates or destroys pairs at a time-dependent strength g_k. Because the operator content is fixed, entirely different physical systems — an expanding universe, a pumped nonlinear crystal — differ only in which two functions ω_k and g_k they supply.

They then take the standard expansion of a two-mode squeezed state in paired number states |n, n⟩, whose coefficients form a geometric series, and repackage the whole state into the standard Bogoliubov form with coefficients α_k and β_k. Crucially, the deformed (Meixner-polynomial) state they inherited from earlier work can be brought into exactly this form, with the deformation parameters u₁ and u₂ absorbed into λ_k = β_k/α_k.

Next, they insert this ansatz into the Schrödinger equation and match coefficients of each paired number state. Instead of two coupled equations for the squeezing strength and phase, out drops one complex Riccati equation for λ_k.

Tracing out one mode of a geometric state gives a geometric probability distribution, so the reduced density matrix is immediately diagonal with entries (1 − q_k)q_k^n. Every entropy then becomes a closed-form function of q_k = |λ_k|² — no matrix reconstruction or diagonalization per model is needed.

Finally, they check the machinery against two known settings: standard cosmological perturbation theory (where the background supplies ω_k and g_k) and an optical parametric amplifier (where the resonant constant-coupling limit must give the textbook tanh τ). They then generalize the same group-theoretic structure from the vacuum seed (d = 0) to excited "lowest-weight" seeds with fixed number imbalance d, which yields a negative-binomial rather than geometric spectrum.

Why This Matters

For research: The framework provides a common coordinate in which cosmological particle production and laboratory squeezing can be compared directly, since both enter only through the two functions ω_k and g_k. It also reduces the computational cost of quantum-information diagnostics: instead of evolving two real variables through rapidly oscillating phases and then reconstructing a reduced density matrix for each model, one solves a single complex equation and reads off purity and entropies algebraically. The non-Gaussian extension shows the same ratio dynamics survives beyond the Gaussian vacuum sector, where the scalar ratio alone is not a complete state descriptor.

Real-world applications (as motivated in the paper):

• Cosmology: quantifying the quantum-to-classical transition of primordial perturbations, their entropy, decoherence, and observable quantum correlations during inflation, including modifications from nontrivial propagation speed, K-essence dynamics, and modified dispersion relations.
• Nonlinear and quantum optics: pulsed squeezing experiments, multimode optical analyses of parametric amplifiers, and chirped or quasi-phase-matched parametric amplification protocols.
• Superconducting circuits and dynamical Casimir physics: parametrically generated photon pairs in circuit-QED devices and the dynamical Casimir effect, both of which share the same quadratic pairing algebra.
• Continuous-variable quantum information: Gaussian-state descriptions of squeezing and mode correlations used in CV protocols.

Industry relevance: squeezed-light sources are deployed in precision metrology and interferometric sensing (e.g., gravitational-wave detectors), and parametric amplifiers are workhorses of quantum optics. A compact diagnostic that predicts one-mode purity and entropy directly from device parameters (pump profile, chirp, detuning) is useful for designing and benchmarking such sources, and for estimating how much quantum information is inaccessible when only one half of a correlated pair is measured.

Future Directions

• Completing the quantitative case studies: the cosmological and chirped-OPA sections are promised to contain numerical solutions of the Riccati equation; reporting explicit trajectories of λ_k, purity, and entropies for realistic backgrounds and pulse profiles is the natural next step.

• Non-Gaussian sectors in practice: the d ≥ 1 negative-binomial reductions are established formally, but their behaviour under specific driving protocols (and the limits of the scalar ratio as a state descriptor when extra representation labels matter) remains to be worked out.

• Beyond quadratic pairing: the authors note the framework may extend to other time-dependent quadratic bosonic systems, provided the evolved state retains a normalized paired Bogoliubov representation — a condition worth testing against states with deformation, open-system effects, or non-Gaussianity from other sources.

• Connecting to other formalisms: the paper points to Krylov-space techniques and two-mode information geometry as complementary descriptions; making the ratio-based diagnostics interoperable with those tools, and with covariance-matrix invariants, is an open avenue.

Target Audience

Advanced graduate students and researchers in quantum field theory in curved spacetime, quantum optics, and continuous-variable quantum information, particularly those already working with Bogoliubov coefficients, two-mode squeezing, or Gaussian-state entanglement measures. The paper is also useful for cosmologists who compute inflationary purity and entropy and want a lighter-weight calculational route, and for quantum-optics experimentalists seeking a compact parameter-level diagnostic for pulsed squeezing. Beginners will find the state expansion and entropy formulas accessible, but the Riccati dynamics and SU(1,1) representation theory require prior background.

Authors’ abstract

We formulate a compact dynamical representation of quantum-information diagnostics for time-dependent two-mode bosonic systems in terms of the Bogoliubov ratio $λ_k(η)=β_k(η)/α_k(η)$. For vacuum-evolved pure Gaussian pair states generated by a Hermitian quadratic Hamiltonian, $λ_k$ obeys a closed complex Riccati equation. After one partner mode is traced out, the reduced-state spectrum is determined by $q_k=|λ_k|^2$, from which the purity, linear entropy, Rényi-2 entropy, and von Neumann entropy follow directly. These Gaussian results are mathematically equivalent to those obtained from the standard density-matrix or covariance-matrix formalism; the advantage of the ratio representation is a compact separation between model-dependent dynamics and universal state diagnostics. Writing $λ_k=\sqrt{q_k}e^{iθ_k}$ further exposes how pair production drives radial growth, whereas frequency rotation and detuning act through the phase and can suppress coherent squeezing accumulation. We also establish an exact extension to non-Gaussian $SU(1,1)$ sectors. For a lowest-weight Fock seed with conserved number difference $d$, the same ratio equation governs the evolution, while the reduced spectrum becomes a negative-binomial distribution determined by $(q_k,d)$. The Gaussian formulas are recovered at $d=0$. Cosmological perturbations and a chirped-pulse optical parametric amplifier illustrate the common dynamical mechanism and its information-theoretic consequences.