Parametric Dynamical Casimir Effect

• Parametric-DCE is the generation of real quanta from vacuum fluctuations through nonadiabatic modulation of system parameters like mirror position and refractive index.
• The phenomenon is modeled by a parametric-oscillator equation and Bogoliubov transformations, capturing resonance-enhanced pair production and vacuum squeezing.
• Various platforms—from SQUID metamaterials to optomechanical cavities—demonstrate DCE, paving the way for advances in quantum computing, sensing, and metrology.

Parametric dynamical Casimir effect (Parametric-DCE) is the generation of real quanta from quantum vacuum fluctuations by a time-dependent modulation of a system parameter or boundary condition, so that a cavity or mode behaves as a parametric amplifier rather than as a static resonator. In contemporary formulations the modulated quantity can be a mirror position, cavity length, effective refractive index, Josephson boundary inductance, qubit transition frequency, mechanical spring constant, or an interaction parameter such as the $s$-wave collision frequency of a Bose–Einstein condensate. The common structure is a nonadiabatic time dependence of $\omega_k(t)$, Bogoliubov mixing between annihilation and creation operators, and resonance-enhanced pair production near $\Omega \simeq 2\omega_k$ or related multimode sum rules (Hsiang et al., 2024, Motazedifard et al., 2017).

1. Conceptual and mathematical structure

A minimal description of Parametric-DCE is the parametric-oscillator equation

$\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$

which appears both in cavity DCE and in cosmological particle creation. The corresponding in/out operators are related by a Bogoliubov transformation,

$$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$

with created-particle number $N_k=|\beta_k|^2$. In this sense, Parametric-DCE is not defined by any specific platform but by the mode-mixing mechanism that converts zero-point fluctuations into real excitations (Xie et al., 2023, Hsiang et al., 2024).

For a single cavity mode, several equivalent Hamiltonian forms recur. A direct time-dependent-frequency formulation is

$$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$

while the standard effective parametric form is

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$

Both encode the same physical content: the appearance of anomalous $a^2$ and $a^{\dagger 2}$ terms that squeeze the vacuum and create pairs (Dodonov, 2012, Macrì et al., 2017).

A persistent misconception is that DCE requires literal relativistic mirror motion. The literature represented here treats rapid modulation of effective boundaries or medium parameters as fully equivalent at the Hamiltonian level. Flux-modulated SQUID metamaterials, dispersion-oscillating fibres, driven optomechanical boundaries, and time-dependent matter parameters all realize Parametric-DCE through the same parametric-amplification mechanism. Conversely, time dependence alone is not sufficient: in the exact $\omega_k(t)$0-dimensional conformal ring model, the massless conformal scalar exhibits no particle creation, $\omega_k(t)$1, even though the geometry is dynamical; the backreaction there is anomaly-driven rather than particle-creation-driven (Lähteenmäki et al., 2011, Vezzoli et al., 2018, Xie et al., 2023).

2. Representative physical platforms

Parametric-DCE has been implemented or proposed across microwave, optical, and hybrid matter platforms. In a Josephson metamaterial embedded in a microwave cavity, flux modulation of a SQUID-based transmission line changes the effective wave velocity by a few percent and produces energy-correlated photons with a bimodal spectrum around a $\omega_k(t)$2 GHz cavity resonance, in agreement with DCE theory (Lähteenmäki et al., 2011). In a dispersion-oscillating photonic crystal fibre, periodic spatial modulation of the group-velocity dispersion becomes a purely temporal modulation in the co-moving frame of a short pump pulse; the observed correlated sidebands at $\omega_k(t)$3 nm and $\omega_k(t)$4 nm, a best $\omega_k(t)$5, and heralded $\omega_k(t)$6 were interpreted as an optical analogue of DCE (Vezzoli et al., 2018).

A different route uses electro-mechanical boundary modulation. In a piezoelectric FBAR directly coupled to a superconducting cavity, the capacitance

$\omega_k(t)$7

implements a time-dependent electromagnetic boundary, leading to a parametric interaction

$\omega_k(t)$8

For the reported baseline parameters, the predicted output at $\omega_k(t)$9 is of order $\Omega \simeq 2\omega_k$0 photons/s per unit bandwidth, with a much smaller purely mechanical contribution in the same device (Sanz et al., 2017).

In reversed-dissipation-regime circuit optomechanics, coherent modulation of the drive frequency allows adiabatic elimination of the mechanical mode and yields an effective intracavity Hamiltonian

$\Omega \simeq 2\omega_k$1

with $\Omega \simeq 2\omega_k$2 and $\Omega \simeq 2\omega_k$3. Here the DCE term is intrinsically dressed by an optomechanical Kerr nonlinearity, which saturates the radiated Casimir-photon number and can induce oscillatory dynamics, sub-Poissonian windows, negative Wigner function, and quadrature squeezing (Solki et al., 1 Aug 2025).

Platform	Modulated quantity	Characteristic feature
Josephson metamaterial cavity	SQUID flux / effective wave velocity	Bimodal spectrum and energy-correlated photons
Dispersion-oscillating fibre	Co-moving temporal GVD modulation	Correlated sidebands and heralded antibunching
FBAR–superconducting cavity	Piezoelectric capacitance modulation	Measurable microwave photon generation rate
RDR circuit optomechanics	Drive-frequency modulation after mechanical elimination	Kerr-modified Casimir photons and nonclassical microwave output

These realizations suggest that Parametric-DCE is best understood as a family of squeezing-driven nonequilibrium instabilities rather than as a single moving-boundary protocol.

3. Dual-phononic Parametric-DCE in a hybrid BEC–optomechanical cavity

A particularly elaborate realization is the coupled BEC–optomechanical Fabry–Pérot cavity with one movable end mirror and an interacting cigar-shaped Bose–Einstein condensate. The relevant degrees of freedom are a cavity mode $\Omega \simeq 2\omega_k$4, a mechanical mode $\Omega \simeq 2\omega_k$5, and a Bogoliubov mode $\Omega \simeq 2\omega_k$6, with red-detuned operating point

$\Omega \simeq 2\omega_k$7

weak linearized couplings $\Omega \simeq 2\omega_k$8, and good-cavity condition $\Omega \simeq 2\omega_k$9. Two coherent modulations are applied: a mechanical spring modulation

$\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$0

and an $\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$1-wave collision-frequency modulation

$\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$2

Under the rotating-wave approximation these produce degenerate parametric drives of $\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$3 and $\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$4, directly amplifying mechanical and atomic vacuum fluctuations and generating mechanical-type and Bogoliubov-type Casimir phonons (Motazedifard et al., 2017).

The crucial feature is that the cavity field is not driven parametrically in any direct sense. Instead, the two phononic pumping channels induce an effective cavity gain in frequency space,

$\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$5

which multiplies $\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$6 in the cavity Langevin equation. Casimir photons are therefore generated indirectly, through phonon-mediated amplification of intracavity vacuum fluctuations. The same formalism also yields cross-induced phonon channels: atomic modulation contributes to the effective gain of the mechanical mode, and mechanical modulation contributes to the effective gain of the Bogoliubov mode.

The operative distinction is between a coherent and a dissipative cavity-gain regime. When the optomechanical couplings, or equivalently the cooperativities,

$\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$7

are strongly mismatched, $\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$8 or $\ddot f_k(t)+\omega_k^2(t)f_k(t)=0,$9, the effective cavity gain becomes dominantly real,

$$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$0

and the cavity behaves like a coherent optical parametric amplifier. This is the regime favoring strong Casimir-photon production. By contrast, for $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$1 the effective gain is dissipative-dominated and photon generation is suppressed, although phonon generation persists. Stability is controlled by the dressed damping rates $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$2 and the Routh–Hurwitz bounds

$$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$3

Within the stable amplification window, the generated particle numbers are externally controllable by the modulation amplitudes $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$4, $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$5, and the cooperativity hierarchy.

The reported scaling trends are correspondingly specific. With mechanical-only modulation, mechanical vacuum amplification begins at $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$6, whereas substantial photon generation requires $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$7 for $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$8 and high $$a_k^{\mathrm{out}}=\alpha_k a_k^{\mathrm{in}}+\beta_k a_k^{\mathrm{in}\dagger},$$9. With atomic-only modulation, photons remain negligible when $N_k=|\beta_k|^2$0 but grow strongly for $N_k=|\beta_k|^2$1. With both modulations present, both phonon species and photons are enhanced further, and photon generation is maximized in the coherent regime of strongly unequal cooperativities.

4. Strong-coupling and subharmonic regimes

The weak-coupling parametric picture does not exhaust Parametric-DCE. In a fully quantum optomechanical treatment without linearization or perturbative reduction, the Hamiltonian

$N_k=|\beta_k|^2$2

generates a ladder of vacuum Casimir–Rabi splittings. Resonances occur at

$N_k=|\beta_k|^2$3

so photon pairs can be created even for $N_k=|\beta_k|^2$4. The dynamics then becomes a coherent exchange between dressed states, for example $N_k=|\beta_k|^2$5, rather than the unbounded growth of an ideal lossless DPA. This nonperturbative regime supports vacuum Casimir–Rabi oscillations, detectable photon fluxes in the weak-coupling limit, and mirror–field entanglement quantified by the negativity (Macrì et al., 2017).

A complementary extension arises from higher-order modulation harmonics in moving-boundary DCE. Expanding the cavity-length modulation beyond first order produces resonance conditions of the form

$N_k=|\beta_k|^2$6

so that the mechanical frequency can be smaller than the lowest cavity-mode frequency by a factor $N_k=|\beta_k|^2$7. The corresponding Floquet exponent is

$N_k=|\beta_k|^2$8

with instability threshold

$N_k=|\beta_k|^2$9

This establishes a subharmonic Parametric-DCE regime: lower pump frequencies become accessible, but the effective gain is suppressed factorially with $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$0 (Ordaz-Mendoza et al., 2016).

In ultrastrong light–matter coupling, fast switching of the interaction itself becomes the parametric drive. For the time-dependent quantum Rabi Hamiltonian, sudden on/off switching of $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$1 activates counter-rotating channels and generates photons from the vacuum. In the state-transfer protocol analyzed in this setting, the fidelity becomes non-monotonic in $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$2, with maxima at $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$3 and minima at $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$4, while the mean photon number is anticorrelated with fidelity. Parametric-DCE therefore appears ոչ as a resource only, but also as a fundamental limit on fast Rabi-oscillation-based protocols in the ultrastrong regime (Benenti et al., 2019).

5. Dissipation, backreaction, and self-consistent dynamics

Once the modulated boundary or medium is treated as dynamical rather than externally prescribed, Parametric-DCE becomes a backreaction problem. In the self-consistent ring model with metric $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$5, adiabatic regularization yields the renormalized energy density

$$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$6

and the effective action for the circumference $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$7,

$$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$8

In this $$\hat H_c=\omega_t \hat n+i\chi_t(\hat a^{\dagger 2}-\hat a^2), \qquad \chi_t=(4\omega_t)^{-1}\frac{d\omega_t}{dt},$$9-dimensional conformal case there is no particle creation; the backreaction comes entirely from the trace anomaly and accelerates collapse relative to static Casimir evolution. In the $$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$0-dimensional rectangular box, by contrast, genuine particle creation produces a radiation-reaction force that opposes motion, implementing the quantum Lenz law (Xie et al., 2023).

This distinction is central to broader foundational analyses. The DCE/CPC analogy places Parametric-DCE within quantum field theory in curved spacetime, where the same mode equation, Bogoliubov structure, and adiabatic regularization recur. In the open-system description, integrating out the quantum field generates nonlocal dissipation and colored noise for the source coordinate,

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$1

with fluctuation–dissipation relations linking the noise kernel to the retarded response. The mirror-oscillator-field microphysical model goes further by replacing ideal boundary conditions with an internal oscillator $$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$2 coupled to the field,

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$3

thereby incorporating finite reflectivity, dispersion, and absorption into the boundary dynamics (Hsiang et al., 2024).

Dissipation also reshapes the instability landscape in engineered photonic environments. For an optomechanical cavity embedded in a photonic crystal, projection onto a reduced Floquet sector yields an energy-dependent effective Liouvillian

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$4

with self-energy

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$5

Here parametric instability from the off-diagonal squeezing term competes with resonance instability from the band-induced self-energy. The resulting stationary mode is selected by the condition that the imaginary part of the complex Floquet eigenvalue vanishes, and a narrow-band photonic crystal with a bandgap can reduce the required pump frequency substantially (Tanaka et al., 2019).

6. Correlations, information-theoretic uses, and frontier directions

Because Parametric-DCE is a squeezing process, entanglement production is not incidental but structural. At parametric resonance in a hyperrectangular cavity, the entanglement entropy exhibits a strong dimensional dependence: in $$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$6 dimensions,

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$7

whereas in higher dimensions

$$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$8

with $$H_{\mathrm{eff}}(t)\approx \hbar \omega_c a^\dagger a+\frac{\hbar \epsilon}{2}\left[a^2 e^{-i2\Omega t}+a^{\dagger 2}e^{i2\Omega t}\right].$$9 identifiable with the Lyapunov exponent of the resonant instability. The logarithmic $a^2$0 prefactor is specific to the bosonic field-theory setting with strongly interacting mode ladders, while the linear law reflects a genuine Floquet instability of a decoupled resonant block (Romualdo et al., 2019).

The nonclassical output of Parametric-DCE has been proposed as a computational and metrological resource. In two multimode coplanar-waveguide resonators coupled by a modulated SQUID boundary, blue-sideband drives implement two-mode squeezers and red-sideband drives implement beam splitters, so the DCE platform realizes a Gaussian boson sampler and naturally performs scattershot boson sampling; the output probabilities in the collision-free regime scale with hafnians of submatrices of the Gaussian adjacency matrix (Peropadre et al., 2016). In a single-qubit circuit-QED realization, slow chirping of the qubit-frequency modulation across $a^2$1 and $a^2$2 generates nonclassical cavity states from vacuum whose phase and displacement quantum Fisher information exceed the classical bounds at the same average photon number, even in the presence of dissipation (Costa et al., 13 Aug 2025).

Microwave platforms also make explicit that ideal linear squeezing is not the only possible output. In reversed-dissipation-regime circuit optomechanics, the Kerr-dressed Hamiltonian

$a^2$3

supports two qualitatively distinct regimes: exponential amplification for $a^2$4, and coherent oscillations for $a^2$5. The Kerr term both saturates growth and creates non-Gaussian windows with sub-Poissonian statistics, negative Wigner function, and quadrature squeezing. The proposed source was explicitly identified as having potential applications in quantum information processing, quantum computing, and microwave quantum sensing (Solki et al., 1 Aug 2025).

At the frontier of external perturbations, a three-dimensional ideal cavity with one oscillating mirror and a gravitational wave background exhibits GW-induced sideband resonances in addition to the standard mechanical DCE peak. The mode frequencies acquire the anisotropic modulation

$a^2$6

and on resonance the created-particle number grows as

$a^2$7

The sidebands are spectrally distinct from the mechanical peak, which suggests a clean conceptual signature, but the predicted rates are suppressed by realistic GW strains $a^2$8–$a^2$9, so observation in standard cavities is not currently feasible (Oliveira et al., 28 Jan 2026).

Parametric-DCE thus spans a continuum from linear vacuum squeezing to Kerr-limited non-Gaussian generation, from prescribed modulation to self-consistent backreaction, and from laboratory cavity engineering to analogue-gravity and gravitational-wave thought experiments. Across these settings, its defining content remains the same: time-dependent mode structure, anomalous pair terms, and the controlled conversion of vacuum fluctuations into correlated quanta.