Entanglement Dynamics of Separable Squeezed States in Finite Memory Structured Reservoir
Pith reviewed 2026-05-19 15:02 UTC · model grok-4.3
The pith
Separable squeezed vacuum states generate and control entanglement in structured reservoirs through non-Markovian mechanisms unavailable in Markovian baths.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors show that three mechanisms absent from Markovian dynamics govern the entanglement: a detuning condition freezes entanglement trajectories across reservoir correlation times; birth death and revival of entanglement occur from orthogonal squeezed inputs; and integer-locked beating with square-wave oscillations arises from periodic detuning. The dynamics are tracked with Gaussian covariance matrices evolved under approximate non-Markovian quantum state diffusion and finite-temperature pseudomode embeddings. All three mechanisms remain effective at finite temperature with deviations no larger than 5 percent in cryogenic regimes and 20 percent at moderate thermal occupations.
What carries the argument
Gaussian covariance matrix evolution under the approximate non-Markovian quantum state diffusion (QSD) method with finite-temperature pseudomode embeddings, which incorporates reservoir memory effects for two bosonic modes.
If this is right
- Entanglement can be generated from initially separable squeezed inputs solely through reservoir structure and detuning choices.
- Periodic detuning modulation produces controllable square-wave entanglement oscillations locked to integer multiples.
- The identified mechanisms survive with small errors at finite temperatures typical of cryogenic cavity and optomechanical experiments.
- Structured spectral densities function as tunable resources for creating and modulating continuous-variable entanglement.
Where Pith is reading between the lines
- The detuning-freezing condition could be used to stabilize entanglement against decoherence in platforms where spectral densities are engineered.
- Similar revival and beating patterns may appear in other open quantum systems with structured baths such as superconducting circuits or trapped ions.
- Testing the mechanisms with time-dependent detuning waveforms beyond periodic cases could reveal additional control over entanglement sudden death and birth times.
Load-bearing premise
The reported dynamics rest on an approximate non-Markovian quantum state diffusion method together with pseudomode embeddings whose accuracy for the entanglement measures has not been checked against exact master-equation solutions.
What would settle it
Exact numerical simulation of the two-mode master equation for a chosen Lorentzian spectral density and specific detuning values would show whether the predicted times of entanglement freezing and revival match the approximate QSD results within the stated temperature deviation bounds.
Figures
read the original abstract
Entanglement in continuous-variable Gaussian systems is a key resource, and common reservoirs can both suppress and generate correlations. Existing work focused on pre-entangled states or Markovian baths, leaving open whether separable squeezed inputs entangle in structured environments or under modulation. We study two bosonic modes coupled to a common reservoir, each initialized in a separable squeezed vacuum. Dynamics are analyzed utilizing Gaussian covariance methods, evolved under approximate Non-Markovian quantum state diffusion (QSD), finite-temperature pseudomode embeddings, and Bures-based non-Markovian diagnostics. We identify three mechanisms absent in Markovian dynamics: (1) A detuning condition that freezes entanglement trajectories across reservoir correlation times; (2) birth, death, and revival of entanglement from orthogonal inputs; and (3) integer-locked beating with square-wave oscillations produced by periodic detuning. All mechanisms persist at finite temperature, with deviations bounded within 5% in cryogenic regimes and 20% at moderate occupations. These deviation bounds align with cryogenic cavity, phononic, and optomechanical platforms, where structured spectral densities and detuning modulation are already accessible. Structured reservoirs are shown to emerge as tunable entanglement resources for continuous-variable quantum technologies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the entanglement dynamics of two bosonic modes initialized in separable squeezed vacuum states and coupled to a common structured non-Markovian reservoir. Using Gaussian covariance evolution under an approximate Non-Markovian quantum state diffusion (QSD) method combined with finite-temperature pseudomode embeddings and Bures-based diagnostics, the authors identify three mechanisms absent from Markovian dynamics: (1) a detuning condition that freezes entanglement trajectories across reservoir correlation times, (2) birth, death, and revival of entanglement from orthogonal inputs, and (3) integer-locked beating with square-wave oscillations under periodic detuning. These effects are reported to persist at finite temperature with deviations bounded by 5% in cryogenic regimes and 20% at moderate occupations, positioning structured reservoirs as tunable resources for continuous-variable quantum technologies.
Significance. If the reported mechanisms prove robust under exact dynamics, the work would extend entanglement control in Gaussian systems beyond pre-entangled states and Markovian baths, offering concrete protocols (detuning freeze, periodic modulation) relevant to optomechanical and phononic platforms. The use of efficient Gaussian covariance methods for continuous-variable systems is a methodological strength that enables exploration of finite-memory effects.
major comments (1)
- [Methods (QSD implementation) and Results (entanglement trajectories)] The central claims rest on the accuracy of the approximate Non-Markovian QSD scheme and pseudomode embeddings for the entanglement measures (Bures distance and logarithmic negativity). No quantitative benchmarking against exact methods such as HEOM or direct integration of the extended master equation is provided for the plotted trajectories or the reported oscillation amplitudes. This validation is load-bearing for asserting that the three mechanisms are genuine non-Markovian features rather than truncation or approximation artifacts (see the methods description of the QSD unraveling and the results sections presenting the detuning-freeze and beating dynamics).
minor comments (2)
- [Introduction or Methods] Clarify the precise definition and implementation of the 'Bures-based non-Markovian diagnostics' when first introduced, including any relation to the covariance matrix evolution.
- [Results (finite-temperature section)] The abstract states deviation bounds of 5% and 20%; the main text should include explicit error estimates or sensitivity plots showing how these bounds were obtained across the parameter space.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed report. The main concern is the absence of quantitative benchmarking of the approximate QSD scheme against exact solvers. We address this point directly below and indicate the revisions we will make to strengthen the manuscript.
read point-by-point responses
-
Referee: The central claims rest on the accuracy of the approximate Non-Markovian QSD scheme and pseudomode embeddings for the entanglement measures (Bures distance and logarithmic negativity). No quantitative benchmarking against exact methods such as HEOM or direct integration of the extended master equation is provided for the plotted trajectories or the reported oscillation amplitudes. This validation is load-bearing for asserting that the three mechanisms are genuine non-Markovian features rather than truncation or approximation artifacts (see the methods description of the QSD unraveling and the results sections presenting the detuning-freeze and beating dynamics).
Authors: We agree that explicit benchmarking would increase in the quantitative amplitudes. The approximate QSD unraveling with pseudomode embedding was selected for its ability to evolve the full Gaussian covariance matrix efficiently over long times and at finite temperature, where HEOM or direct integration of a high-dimensional extended master equation becomes prohibitive. The three reported mechanisms (detuning freeze, birth-death-revival, and integer-locked beating) originate from the structure of the spectral density and the detuning condition; they appear already at the level of the memory kernel and are therefore expected to survive under exact dynamics. Nevertheless, to meet the referee’s standard we will add a dedicated validation subsection in the Methods. This will include (i) short-time comparisons of the covariance evolution against the exact pseudomode master equation for representative parameter sets, (ii) convergence tests with respect to the number of pseudomodes, and (iii) error estimates based on the neglected higher-order terms in the QSD expansion. These additions will be placed before the results sections so that readers can assess the reliability of the plotted trajectories. revision: yes
Circularity Check
No circularity detected; mechanisms emerge from numerical evolution of Gaussian covariances under QSD and pseudomode methods
full rationale
The paper evolves separable squeezed vacuum states via Gaussian covariance matrices under an approximate Non-Markovian QSD scheme combined with finite-temperature pseudomode embeddings. The reported mechanisms (detuning-induced freezing, birth-death-revival from orthogonal inputs, and integer-locked square-wave beating) are extracted as dynamical outcomes of this evolution rather than being presupposed by the inputs, fitted parameters, or self-referential definitions. No load-bearing self-citations, uniqueness theorems imported from prior author work, or ansatzes smuggled via citation appear in the provided text. The analysis is self-contained as a forward simulation whose results are independent of the target quantities by construction.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Zero-Temperature entanglement trajectories We first examine the Markovian limit in which two un- coupled bosonic modes couple symmetrically to a memo- ryless reservoir, with entanglement quantified by the log- arithmic negativity defined in Eq. 8. Figure 1(a) shows trajectories ofE N(t) for modeAfixed ats a = +1 and modeBprepared withs b =−0.5 (blue),s b ...
-
[2]
Thermal suppression of steady state entanglement Including finite thermal occupation rescales the partial-transpose symplectic spectrum by √1 + 2¯n, thereby tightening the separability condition and sup- pressing steady correlations. Figure 2 shows that aligned inputs retain nonzero entanglement up to ¯n≈1.5, while weakly orthogonal inputs become separabl...
-
[3]
Figure 3(a) compares the entan- glement trajectories forγ/κ= 0.5 (blue) andγ/κ= 5 (red)
Zero-temperature freezing law We next consider an Ornstein–Uhlenbeck reservoir at resonance, where memory effects are controlled by the correlation time 1/γ. Figure 3(a) compares the entan- glement trajectories forγ/κ= 0.5 (blue) andγ/κ= 5 (red). Despite a tenfold change in memory scale, the curves nearly overlap, oscillating betweenE N ≈0.2 and EN ≈1.7 a...
-
[4]
Figure 4 shows the relative de- viation ofE N(κt) as a function of bath occupation ¯n
Finite-temperature robustness To assess stability against thermal noise, we embed the Ornstein–Uhlenbeck kernel into a pseudomode rep- resentation that extends the ˆO0 closure to finite tem- perature [28, 44, 47]. Figure 4 shows the relative de- viation ofE N(κt) as a function of bath occupation ¯n. Deviations remain within 5% for cryogenic ¯n≤0.05 and wi...
-
[5]
Figure 5(a) shows that the entanglement response depends sensitively on detuning
Revival and long-lived correlations in structured baths We now initialize the system with orthogonal squeez- ing,s a = +1 ands b =−1, and examine the dynamics in an Ornstein–Uhlenbeck reservoir. Figure 5(a) shows that the entanglement response depends sensitively on detuning. For (δ AB, δAE)/κ= (−0.88,0.1) the loga- rithmic negativity peaks near 0.65 befo...
-
[6]
5(b) shows the same orthogonal input evolved in a memoryless reservoir
Markovian baseline comparison For comparison, Fig. 5(b) shows the same orthogonal input evolved in a memoryless reservoir. In this case, the entanglement dynamics reduce to brief birth–death events. Atδ AB/κ=−1.62 a small peak ofE N ≈0.17 appears nearκt≈2 before vanishing within four time units, while atδ AB/κ=−0.88 the response is limited to a minor plat...
-
[7]
11 provides a complementary probe of environmental memory
Non-Markovian witness and entanglement mismatch The integrated witnessNdefined in Eq. 11 provides a complementary probe of environmental memory. Fig- ure 5(c) shows its dependence on mode detuning. For δAE/κ= 10 the witness is identically zero, even though Fig. 5(a) exhibits substantial and long-lived entangle- ment. Atδ AE/κ= 1 the witness becomes finite...
-
[8]
Beating dynamics at integer modulation amplitudes We now introduce a sinusoidal modulation of the mode detuning,δ AB(t) =δ 0(sin(t)+1), applied to the orthogo- nal input state. Figure 6 demonstrates that robust beat- ing envelopes appear only when the modulation ampli- tude takes integer values. Forδ 0/κ=−2, the logarith- mic negativity exhibits persisten...
-
[9]
Integer locking of time-averaged entanglement The long-time average of the logarithmic negativity, defined as⟨E N ⟩T = 1 T R T 0 EN(t)dt, provides a global measure of modulation-induced stabilization. Figure 7(a) shows that for large environment detuning (δAE/κ= 10) this average develops sharp peaks at integerδ 0/κvalues between−1 and−4, consistent with t...
-
[10]
Backaction response under modulation The modulation-induced memory is further assessed using the integrated witnessNdefined in Eq. 11. Fig- ure 7(b) shows that forδ AE/κ= 10 the witness van- ishes identically, even in regimes where Fig. 7(a) reveals strong integer-locked entanglement. At moderate detun- ing,δ AE/κ= 1,Nexhibits quasi-periodic oscillations,...
-
[11]
Figure 8 reports the relative deviation ofE N(t) from the zero-temperature baseline atδ 0/κ=−2
Finite-temperature persistence of beating Finally, we test the stability of periodically driven beating against thermal excitations. Figure 8 reports the relative deviation ofE N(t) from the zero-temperature baseline atδ 0/κ=−2. In the cryogenic regime (¯n≤ 0.05), deviations remain below 5%, increasing to at most FIG. 8: Relative deviation %E N =|∆E N /Eb...
-
[12]
II A, which acts symmetrically on both modes
Bright-Dark Mode Decomposition The system couples to the reservoir through the col- lective operator of Sec. II A, which acts symmetrically on both modes. We make this structure explicit by introduc- ing bright and dark modes through a unitary rotation of the ladder operators [3, 49, 50] ˆd± = 1√ 2(ˆa±ˆb).(A1) In this basis, the free Hamiltonian separates...
-
[13]
Simon Separability Criteria Entanglement in Gaussian states is certified by the Peres–Horodecki partial transpose condition, which has an explicit covariance-matrix form [1, 13, 14]. For the squeezed-vacuum inputs studied here, the pairing of a vacuum bright mode and a displaced squeezed thermal dark mode allows one to express the relevant steady state dy...
-
[14]
Separability holds when the smallest symplectic eigen- value satisfies ˜ν− ≥ 1 2 [13]
Symplectic Eigenvalue Formulas Entanglement in Gaussian systems is diagnosed by the spectrum of the partially transposed covariance matrix. Separability holds when the smallest symplectic eigen- value satisfies ˜ν− ≥ 1 2 [13]. This condition underlies logarithmic negativity, a computable monotone defined directly from ˜ν− [1, 41]. For the squeezed-coheren...
-
[15]
17] simplify and admit both a transient so- lution and a universal steady-state limit
Mode Resonant Steady-State At resonance (δ AB = 0), the coefficient equations of motion [Eq. 17] simplify and admit both a transient so- lution and a universal steady-state limit. The steady- state formF 1,2(∞) andχ(γ, δ AE) are given in Methods Eqs. 17, 18. The transient solution was used to plot ana- lytic trajectories for resonant freezing, avoiding re...
-
[16]
Methods already list F1,2(∞) Eq
Freezing Law and Stability At exact resonance (δ AB = 0), the convolution coef- ficients admit analytic inspection. Methods already list F1,2(∞) Eq. 17,χ(γ, δ AE) Eq. 18, and the freezing in- equalities Eqs. 19, 20. Here we sketch how these results follow from the transient solution and show why they are stable. These derivations underlie the trajectory o...
-
[17]
Detuning-Based Phase Accumulation A direct demonstration of phase accumulation fol- lows from the detuning term in the coefficient equations Eq. 17. SubstitutingF 2(t)→e iδAB t ˜F2(t) isolates the oscillatory part and yields an evolution equation with a time-dependent source and a global phase. d dt ˜F2(t) = γ√κ 2 e−iδAB t+ (B7) (√κ(F1+eiδAB t ˜F2)−(γ+iδ ...
-
[18]
Entanglement in continuous variable systems: Recent advances and current perspectives
G. Adesso and F. Illuminati, J. Phys. A: Math. Theor. 40, 7821 (2007), arXiv:quant-ph/0701221
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[19]
C. Weedbrook, S. Pirandola, R. Garcia-Patron, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, Rev. Mod. Phys.84, 621 (2012), arXiv:1110.3234 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[20]
Gaussian states in continuous variable quantum information
A. Ferraro, S. Olivares, and M. G. A. Paris, Gaus- sian states in continuous variable quantum information (2005), arXiv:quant-ph/0503237
work page internal anchor Pith review Pith/arXiv arXiv 2005
- [21]
-
[22]
F. Benatti, R. Floreanini, and M. Piani, Phys. Rev. Lett. 91, 070402 (2003)
work page 2003
-
[23]
Entangled states and collective nonclassical effects in two-atom systems
Z. Ficek and R. Tanas, Physics Reports372, 369 (2002), arXiv:quant-ph/0302082
work page internal anchor Pith review Pith/arXiv arXiv 2002
-
[24]
J. S. Prauzner-Bechcicki, Two-mode squeezed vacuum state coupled to the common thermal reservoir (2004), arXiv:quant-ph/0211114
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[25]
A. Wolf, G. D. Chiara, E. Kajari, E. Lutz, and G. Morigi, EPL95, 60008 (2011), arXiv:1102.1838 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2011
- [26]
- [27]
-
[28]
H.-P. Breuer, E.-M. Laine, J. Piilo, and B. Vacchini, Rev. Mod. Phys.88, 021002 (2016)
work page 2016
- [29]
-
[30]
R. Simon, Phys. Rev. Lett.84, 2726 (2000), publisher: American Physical Society
work page 2000
-
[31]
L.-M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, Phys. Rev. Lett.84, 2722 (2000), publisher: American Physical Society
work page 2000
-
[32]
T. Linowski, C. Gneiting, and L. Rudnicki, Phys. Rev. A 102, 042405 (2020)
work page 2020
- [33]
- [34]
-
[35]
Y. Zhao, F. Zheng, J. Liu, and Y. Yao, Physica A: Sta- tistical Mechanics and its Applications423, 80 (2015)
work page 2015
- [36]
- [37]
- [38]
-
[39]
A. A. Valido, D. Alonso, and S. Kohler, Phys. Rev. A88, 042303 (2013), arXiv:1307.4571 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [40]
-
[41]
T. Yu, L. Di´ osi, N. Gisin, and W. T. Strunz, Phys. Rev. A60, 91 (1999), arXiv:quant-ph/9902043
work page internal anchor Pith review Pith/arXiv arXiv 1999
- [42]
-
[43]
S. M. Barnett and S. Stenholm, Phys. Rev. A64, 033808 (2001)
work page 2001
-
[44]
S. Einsiedler, A. Ketterer, and H.-P. Breuer, Phys. Rev. A102, 022228 (2020)
work page 2020
- [45]
-
[46]
C.-H. Bai, D.-Y. Wang, S. Zhang, S. Liu, and H.- F. Wang, Annalen der Physik531, 1800271 (2019), arXiv:1811.06171 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[47]
J. C. Gonzalez-Henao, E. Pugliese, S. Euzzor, S. Ab- dalah, R. Meucci, and J. A. Roversi, Sci Rep5, 13152 (2015)
work page 2015
-
[48]
J. Yang, T.-X. Lu, M. Peng, J. Liu, Y.-F. Jiao, and H. Jing, Opt. Express32, 785 (2024)
work page 2024
-
[49]
M. J. Woolley and A. A. Clerk, Phys. Rev. A89, 063805 (2014)
work page 2014
-
[50]
C.-G. Liao, R.-X. Chen, H. Xie, and X.-M. Lin, Phys. Rev. A97, 042314 (2018)
work page 2018
-
[51]
H.-P. Breuer, E.-M. Laine, and J. Piilo, Phys. Rev. Lett. 103, 210401 (2009)
work page 2009
- [52]
- [53]
-
[54]
H. Krauter, C. A. Muschik, K. Jensen, W. Wasilewski, J. M. Petersen, J. I. Cirac, and E. S. Polzik, Phys. Rev. Lett.107, 080503 (2011)
work page 2011
-
[55]
C. F. Ockeloen-Korppi, E. Damsk¨ agg, J.-M. Pirkkalainen, M. Asjad, A. A. Clerk, F. Massel, M. J. Woolley, and M. A. Sillanp¨ a¨ a, Nature556, 478 (2018)
work page 2018
-
[56]
T. Park, H. Stokowski, V. Ansari, S. Gyger, K. K. S. Multani, O. T. Celik, A. Y. Hwang, D. J. Dean, F. Mayor, T. P. McKenna, M. M. Fejer, and A. Safavi-Naeini, Sci. Adv.10, eadl1814 (2024). 15
work page 2024
-
[57]
Y. Zhao and G. Chen, Physica A: Statistical Mechanics and its Applications317, 13 (2003)
work page 2003
-
[58]
G. Vidal and R. F. Werner, Phys. Rev. A65, 032314 (2002), publisher: American Physical Society
work page 2002
-
[59]
S. Maniscalco, S. Olivares, and M. G. A. Paris, Phys. Rev. A75, 062119 (2007)
work page 2007
- [60]
-
[61]
G. Pleasance, B. M. Garraway, and F. Petruccione, Phys. Rev. Research2, 043058 (2020)
work page 2020
-
[62]
H.-P. Breuer and F. Petruccione,The theory of open quantum systems, 1st ed. (Clarendon Press, Oxford, 2009)
work page 2009
-
[63]
Lindblad, Communications in Mathematical Physics 48, 119 (1976)
G. Lindblad, Communications in Mathematical Physics 48, 119 (1976)
work page 1976
-
[64]
S. Luo, N. Lambert, P. Liang, and M. Cirio, PRX Quan- tum4, 030316 (2023)
work page 2023
-
[65]
B. M. Garraway, Phys. Rev. A55, 2290 (1997)
work page 1997
- [66]
-
[67]
C.-H. Chou, T. Yu, and B. L. Hu, Phys. Rev. E77, 011112 (2008)
work page 2008
- [68]
-
[69]
C. Enss and H. S.,Low-Temperature Physics(Springer- Verlag, Berlin/Heidelberg, 2005)
work page 2005
-
[70]
S. Krinner, S. Storz, P. Kurpiers, P. Magnard, J. Heinsoo, R. Keller, J. L¨ utolf, C. Eichler, and A. Wallraff, EPJ Quantum Technology6, 2 (2019)
work page 2019
- [71]
-
[72]
A. Serafini, M. Lostaglio, S. Longden, U. Shackerley- Bennett, C.-Y. Hsieh, and G. Adesso, Physical Review Letters124, 010602 (2020), arXiv:1909.06123 [quant-ph]
- [73]
-
[74]
J. P. Paz and A. J. Roncaglia, Phys. Rev. A79, 032102 (2009), arXiv:0809.1676 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[75]
J. P. Paz and A. J. Roncaglia, Phys. Rev. Lett.100, 220401 (2008)
work page 2008
-
[76]
S. Direkci, K. Winkler, C. Gut, M. Aspelmeyer, and Y. Chen, Universality of stationary entanglement in an optomechanical system driven by non-Markovian noise and squeezed light (2025), arXiv:2502.02979 [quant-ph]
-
[77]
Teklu, Physics Letters A432, 128022 (2022), arXiv:2202.06313 [quant-ph]
B. Teklu, Physics Letters A432, 128022 (2022), arXiv:2202.06313 [quant-ph]
- [78]
-
[79]
K. Mui, A. Couvertier, and T. Yu, APL Quantum2, 016126 (2025)
work page 2025
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.