Dissipative Linear Parametric Oscillator (DLPO)

• DLPO is a linear oscillator with time-modulated stiffness and damping, serving as a critical theory in modeling discrete time crystal instabilities.
• Research leverages various formulations—damped Mathieu, Langevin, and Caldirola–Kanai—to analyze stability, threshold behavior, and resonance structures.
• DLPO frameworks provide practical insights into quantum state dynamics as well as applications in optical parametric systems and nonlinear oscillator linearizations.

Searching arXiv for recent and foundational papers on dissipative linear parametric oscillators and closely related formulations. A dissipative linear parametric oscillator (DLPO) is a linear oscillator in which the restoring coefficient is modulated in time and dissipation is present in the dynamical equation. In an explicit formulation used as an effective critical theory for dissipative discrete time crystals, the DLPO is

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$

where $\gamma$ is the dissipation rate, $\Omega$ is the natural frequency, $\omega_d$ is the external drive frequency, and $A$ is the drive amplitude; the modulation acts on the oscillator “spring constant” or restoring coefficient rather than as an additive force (Jr. et al., 25 Jul 2025). Closely related formulations across the literature include damped Mathieu equations, Kramers/Fokker–Planck generators for thermal-bath-coupled oscillators, Caldirola–Kanai Hamiltonians, non-Markovian influence-functional models, and exact second-order equations with time-dependent damping and frequency (Chaki et al., 2021, Zhu et al., 2016, Pachón et al., 2012, Rosu et al., 2012).

1. Definition and canonical forms

In the narrow sense, the DLPO is the damped parametrically driven linear oscillator

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$

used as a generic effective description for systems near a discrete-time-crystal instability. In that setting, the control parameter is the modulation amplitude $A$, the critical point is a threshold value $A_c$, and the resonant condition is $\omega_d=2\Omega$ (Jr. et al., 25 Jul 2025).

A broader second-order representation, used repeatedly in the DLPO literature, is

$\ddot x+2\gamma(t)\dot x+\Omega^2(t)x=0.$

In the Darboux-generated family of nonsingular parametric oscillators, this appears with explicit time-dependent damping and frequency,

$\gamma$0

or equivalently

$\gamma$1

with $\gamma$2 and $\gamma$3 (Rosu et al., 2012). In the Caldirola–Kanai-based quantum treatment, the same structure is written as

$\gamma$4

which is identified as a damped Mathieu-type equation (Chaki et al., 2021). In the thermal-bath Langevin setting, the oscillator takes the form

$\gamma$5

with white-noise forcing fixed by fluctuation–dissipation (Zhu et al., 2016).

The literature also uses first-order complex-amplitude formulations. For the linearized signal fluctuations of a degenerate optical parametric oscillator, the effective equation is

$\gamma$6

which is exactly a dissipative linear parametric oscillator equation in bosonic-envelope form (Navarrete-Benlloch et al., 2014). A plausible implication is that the DLPO is best understood not as a single notation, but as a model class whose members are equivalent under second-order, first-order-envelope, or phase-space formulations when the dynamics remain linear.

Form	Representative equation	Context
Resonant second-order DLPO	$\gamma$7	DTC onset and KZM (Jr. et al., 25 Jul 2025)
Damped Mathieu form	$\gamma$8	CK quantum/classical correspondence (Chaki et al., 2021)
General time-dependent damped oscillator	$\gamma$9	Darboux-related nonsingular oscillators (Rosu et al., 2012)

2. Instability, threshold, and resonance structure

For the resonantly driven DLPO, the disordered state $\Omega$0 and the resonant state are separated by the critical drive amplitude

$\Omega$1

Near threshold, for $\Omega$2, the response takes the envelope form

$\Omega$3

and the relaxation time diverges as

$\Omega$4

This divergence is the critical-slowing-down ingredient used to establish the adiabatic–impulse approximation and the Kibble–Zurek scaling with $\Omega$5 (Jr. et al., 25 Jul 2025).

In the damped Mathieu formulation, stability is analyzed through Floquet theory and the Monodromy matrix. Writing the problem as a first-order periodic system $\Omega$6 with period $\Omega$7, the motion is stable if $\Omega$8, where $\Omega$9 are the Floquet exponents, equivalently if the Floquet multipliers $\omega_d$0 satisfy $\omega_d$1. Numerical Monodromy analysis shows that the dynamical instability generated by parametric driving is reduced by dissipation: stable bands widen and unstable regions shrink as $\omega_d$2 increases (Chaki et al., 2021).

The high-frequency Magnus treatment of the thermal-bath-coupled oscillator yields a different but related picture. At second order, the effective generator produces a positive restoring-force renormalization,

$\omega_d$3

which corresponds to

$\omega_d$4

Higher orders further renormalize the oscillator frequency and create a weakly renormalized effective temperature,

$\omega_d$5

At fourth order, however, a cutoff-dependent term appears,

$\omega_d$6

which the paper interprets as a precursor to instability and as an indication of the breakdown of the hierarchy $\omega_d$7 (Zhu et al., 2016).

A stochastic analogue arises when the stiffness is modulated by dichotomic noise,

$\omega_d$8

There, stability is characterized by the Lyapunov exponent

$\omega_d$9

with $A$0 stable and $A$1 unstable. In the weak-damping, weak-noise regime,

$A$2

leading to the threshold curve

$A$3

This is a kinetic-statistical instability criterion rather than a Floquet one (Schirdewahn, 2020).

3. Dissipation models and bath structure

The simplest dissipative description uses local friction and white noise. In the Langevin model for a thermal bath,

$A$4

with

$A$5

the corresponding Kramers generator is

$A$6

supplemented by the parametric drive operator $A$7 (Zhu et al., 2016).

A Hamiltonian representation of linear dissipation is provided by the Caldirola–Kanai model,

$A$8

which reproduces the classical equation

$A$9

In this formulation, damping is encoded directly into the Hamiltonian parameters rather than through an explicit bath or Lindblad operator (Chaki et al., 2021).

A more general effective-Hamiltonian treatment begins from the time-dependent quadratic Hamiltonian

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$0

for which the classical center satisfies

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$1

with

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$2

This gives a direct coefficient map from Hamiltonian data to DLPO form (Zhang et al., 2016).

The most general dissipative treatment among the cited works is non-Markovian. For a charged quantum oscillator with arbitrary time-dependent frequency $$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$3, coupled to a thermal bath and blackbody radiation, the stationary paths satisfy

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$4

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$5

where $$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$6. The thermal bath may be chosen Ohmic–Drude,

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$7

while blackbody radiation uses

$$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$8

Here the DLPO becomes a linear parametric oscillator with memory friction and colored quantum noise rather than a local $$\ddot{\theta} + \gamma \dot{\theta}+ \Omega^{2}\left[ 1 + A \sin(\omega_{d}t) \right]\theta = 0,$$9 term (Pachón et al., 2012).

4. Quantum dynamics, invariants, and state dependence

A recurrent result is that quantum DLPO dynamics can be reduced to auxiliary classical equations. In the Liouville–von Neumann treatment of the Caldirola–Kanai oscillator, the invariant ladder operator coefficient $A$0 satisfies

$A$1

the same damped Mathieu equation as the classical coordinate. The ground-state wavefunction is

$A$2

and the Wigner function is a Gaussian whose widths and tilt are controlled by $A$3. Stable classical dynamics corresponds to bounded or reviving localization, while unstable dynamics corresponds to persistent stretching and delocalization (Chaki et al., 2021).

The Lie-transformation approach reaches a related separation from a general quadratic Hamiltonian. Writing

$A$4

with a Heisenberg–Weyl part and an $A$5 part, the parameters $A$6 satisfy

$A$7

and reproduce the external classical motion of the wave-packet center, while $A$8 govern squeezing, dilation, chirp, and spreading. The exact evolved state is

$A$9

This separates external classical dynamics from internal quantum deformation (Zhang et al., 2016).

A more specialized but conceptually important quantum result concerns parametric resonance with time-dependent damping and time-dependent frequency. There the auxiliary functions $A_c$0 and $A_c$1 satisfy

$A_c$2

$A_c$3

Two asymptotic regimes are distinguished. If $A_c$4 remains bounded while $A_c$5 undergoes parametric resonance, then

$A_c$6

so the initial quantum state is forgotten. If $A_c$7 and $A_c$8 undergoes unbounded growth, then

$A_c$9

and the relative fluctuations approach a finite $\omega_d=2\Omega$0-dependent constant, so memory of the initial state persists over arbitrary long time scales (Ferrari, 2019).

5. Relation to nonlinear parametric oscillators

Several influential models are not DLPOs in the strict sense but reduce to them after linearization. In the degenerate optical parametric oscillator, the full nonlinear open-cavity system contains pump and signal modes, but below threshold the signal fluctuation block decouples and obeys

$\omega_d=2\Omega$1

or, in dimensional variables,

$\omega_d=2\Omega$2

The threshold condition is $\omega_d=2\Omega$3, corresponding to $\omega_d=2\Omega$4 at the classical level. The paper’s main contribution is a self-consistent linearization, equivalent to a Gaussian-state closure, that regularizes the effective linear description near threshold (Navarrete-Benlloch et al., 2014).

A related envelope model arises in coupled parametric-oscillator Ising machines. The starting slow-flow equation is

$\omega_d=2\Omega$5

For DLPO interpretation, the natural linear limit is obtained by dropping the cubic saturation,

$\omega_d=2\Omega$6

The paper states explicitly that this linearized equation is not called a DLPO model there; this is an inference rather than a direct claim. What survives in that limit are dissipation or gain, parametric drive, and both coupling channels, while the amplitude-limiting mechanism disappears (Khan et al., 28 Oct 2025).

The same caution applies to the Kerr parametric oscillator. Its laboratory-frame equation,

$\omega_d=2\Omega$7

reduces in the strict linear limit $\omega_d=2\Omega$8 to

$\omega_d=2\Omega$9

Under the rotating-wave approximation, the linearized slow quadratures satisfy

$\ddot x+2\gamma(t)\dot x+\Omega^2(t)x=0.$0

with

$\ddot x+2\gamma(t)\dot x+\Omega^2(t)x=0.$1

The paper emphasizes that the stable phase states and Kramers-turnover physics require the nonlinear Kerr term and therefore do not survive in a strict DLPO (Boneß et al., 21 Apr 2026).

6. Applications, scope, and common distinctions

The DLPO functions as a universal linear critical theory in several settings. In dissipative discrete time crystals, a huge class of systems that can form DTCs can be mapped onto a DLPO, including the low-temperature, small-amplitude Sine-Gordon model after linearization and mode decomposition,

$\ddot x+2\gamma(t)\dot x+\Omega^2(t)x=0.$2

and the open Dicke lattice model after expansion about the normal phase, linearization of fluctuations, decomposition into momentum modes, reduction to polaritonic normal modes, and use of the truncated Wigner approximation (Jr. et al., 25 Jul 2025).

A different application is exact model generation. Darboux-related nonsingular oscillators produce closed-form linear parametric equations with explicit time-dependent damping and stiffness, but the damping coefficient changes sign periodically: $\ddot x+2\gamma(t)\dot x+\Omega^2(t)x=0.$3 The paper therefore emphasizes that these are better described as periodically dissipative / periodically amplifying linear parametric oscillators, or as sign-changing gain–loss oscillators, rather than as standard uniformly dissipative oscillators (Rosu et al., 2012).

The stochastic-frequency problem establishes that a DLPO need not be periodically driven to become unstable. With dichotomic stiffness modulation, the relevant criterion is the Lyapunov exponent rather than a Floquet multiplier, and finite-correlation multiplicative noise can overcome viscous damping even when the mean modulation vanishes (Schirdewahn, 2020). Conversely, the damped Mathieu and Magnus analyses show that ordinary dissipation suppresses coherent parametric instability by shrinking unstable tongues or widening stable regions (Chaki et al., 2021, Zhu et al., 2016).

A recurring distinction in the literature is between strict DLPOs and nonlinear parent systems. Degenerate optical parametric oscillators, Stuart–Landau parametric oscillators, and Kerr parametrons can yield effective linear parametric-oscillator equations near threshold or after linearization, but their stable above-threshold phase states, amplitude saturation, and activated switching phenomena are not properties of the strict linear model itself (Navarrete-Benlloch et al., 2014, Khan et al., 28 Oct 2025, Boneß et al., 21 Apr 2026). This suggests a useful taxonomy: the DLPO is the linear instability kernel, while finite-amplitude bistability and phase selection generally require additional nonlinear saturation or external limiting mechanisms.