---
title: Quantum Duffing Oscillator Dynamics
url: https://www.emergentmind.com/topics/quantum-duffing-oscillator
type: topic
---

# Quantum Duffing Oscillator Dynamics

Searching arXiv for recent and foundational papers on the quantum Duffing oscillator to ground the article in published work.
arxiv_search(query="quantum Duffing oscillator Josephson driven dissipative Floquet metastability", max_results=10, sort_by="relevance")
The quantum Duffing oscillator is a driven nonlinear quantum oscillator whose dynamics are governed by a quartic anharmonicity or, in rotating-frame form, by an effective Kerr term. In the literature summarized here, it appears both as a direct quantization of a periodically driven quartic oscillator and as a reduced model for superconducting nonlinear resonators and Josephson circuits. Its defining problems are the structure of quasienergy spectra under periodic drive, the role of dissipation in producing unique steady states or long-lived metastable manifolds, and the extent to which classical Duffing notions such as bistability, hysteresis, and dynamical switching survive in the quantum regime [1003.1007, 2206.06338, 1812.02839].

## 1. Canonical definitions and model classes

A common laboratory-frame definition is the driven quartic oscillator
\[
\hat{H}_{\rm DO}(t)=\frac{\hat{P}_y^2}{2M}+\frac{M\Omega^2}{2}\hat{y}^2+\frac{\alpha}{4}\hat{y}^4+\hat{y}F\cos(\omega_{ex} t),
\]
with mass \(M\), bare frequency \(\Omega\), quartic nonlinearity \(\alpha\), and monochromatic drive of amplitude \(F\) and frequency \(\omega_{ex}\). In the hard-nonlinearity case \(\alpha>0\), the static potential is monostable, while the driven response acquires the characteristic nonlinear resonance structure of the Duffing problem [1003.1007].

In superconducting implementations, the same driven nonlinear dynamics is often written directly in a rotating frame as a driven-dissipative Kerr oscillator,
\[
H/\hbar = \Delta a^{\dagger}a + U a^{\dagger}a^{\dagger}aa + \xi \left(a+a^{\dagger}\right),
\]
supplemented by single-photon loss \((\gamma/2)\mathcal D[a]\). Here \(a\) and \(a^\dagger\) act on the resonator mode, \(\Delta\) is the detuning, \(U\) the Kerr nonlinearity, \(\xi\) the drive strength, and \(\gamma\) the total energy dissipation rate. This formulation is the operative quantum Duffing model in superconducting nonlinear resonators [2206.06338].

A second, widely used dimensionless form emphasizes a double-well geometry and the semiclassical control parameter:
\[
\ddot{x} + 2 \Gamma \dot{x} +\beta^2x^3 -x = \frac{g}{\beta}\cos(\Omega t).
\]
Its quantum Hamiltonian contains the quartic Duffing term, the inverted quadratic term generating the double well, the periodic drive, and an explicit damping-related symmetrized \(QP\) contribution chosen so that expectation values recover the damped classical limit as \(\beta\to0\). This scaling makes \(\beta\) the effective size or quantumness parameter: small \(\beta\) is the classical limit, while \(\beta\sim1\) is more deeply quantum [2011.00118].

These model classes are not interchangeable in all regimes. The quartic or Kerr form is the quantum Duffing oscillator proper, but some physical systems only reduce to it after controlled approximations. That distinction is central in Josephson implementations.

## 2. Rotating-frame and Floquet structure

Because the Duffing oscillator is periodically driven, Floquet theory provides a natural description. The Schrödinger solutions are written as
\[
|\psi_j(t)\rangle=e^{-i\epsilon_j t/\hbar}|\phi_j(t)\rangle,
\]
with periodic Floquet modes \(|\phi_j(t+T_{\omega_{ex}})\rangle=|\phi_j(t)\rangle\) and quasienergies \(\epsilon_j\). In Sambe space, the quasienergies are defined only modulo \(\hbar\omega_{ex}\), so the relevant physics is organized by avoided crossings and multiphoton resonances rather than by the static energy spectrum alone [1003.1007].

Two complementary perturbative Floquet schemes have been developed for the deep quantum regime. One treats the driven linear oscillator exactly and the quartic nonlinearity perturbatively; it is exact in drive amplitude and valid away from resonance. The other uses Van Vleck quasi-degenerate perturbation theory in Floquet space, perturbative in both drive and nonlinearity, and is designed for the vicinity of an \(N\)-photon resonance. Their overlap provides a controlled description across the driving-frequency axis [1003.1007].

In the undriven nonlinear oscillator, the quartic term removes harmonic level equidistance:
\[
E_j=\hbar\Omega\left(j+\frac12\right)+\frac{3}{8}\alpha y_0^4\left(j(j+1)+\frac12\right),
\qquad
y_0=\sqrt{\frac{\hbar}{M\Omega}}.
\]
Consequently, the resonance condition becomes level dependent. For a one-photon resonance,
\[
\omega_{ex}=\Omega_j=\Omega+\frac{3\alpha y_0^4(j+1)}{4\hbar}.
\]
Near such resonances, the stationary dissipative response can exhibit antiresonant dips rather than resonant peaks; at weak damping the antiresonance is asymmetric, while stronger damping restores an ordinary broadened resonance centered at the shifted frequency [1003.1007].

This Floquet viewpoint changes the interpretation of the driven quantum Duffing oscillator. In the deep quantum regime, the central objects are quasienergy doublets, multiphoton avoided crossings, and bath-controlled stationary Floquet populations, not the classical folded response curve.

## 3. Duffing as a controlled limit of Josephson dynamics

A current-driven Josephson junction provides a physically important case in which Duffing behavior is emergent rather than fundamental. The circuit degree of freedom is the superconducting phase difference \(\varphi\), conjugate to the charge \(q\), with laboratory-frame Hamiltonian
\[
H_{\rm JJ} = \frac{q^2}{2C} - E_J \cos\varphi - F \cos(\omega_d t)\,\varphi .
\]
Introducing \(\varphi=\kappa(a+a^\dagger)\), moving to a frame rotating at \(\omega_d\), and applying the rotating-wave approximation yields an exact-in-amplitude rotating-frame quasienergy operator whose nonlinearity is governed by a Bessel function,
\[
H_{\rm JJ}^{\rm RF} = \hbar \delta\omega\, \hat n -\frac{F}{2}\kappa (a+a^\dagger) - E_J^*\left[\, :J_0\!\left(2\kappa \sqrt{\hat n}\right): + \kappa^2 \hat n \right].
\]
The Duffing Hamiltonian appears only after a further small-amplitude expansion of this Bessel structure [1812.02839].

The controlled Duffing regime requires four simultaneous conditions: small rotating-frame amplitude \(\kappa\sqrt{\langle \hat n\rangle}\ll1\), weak enough drive such that \(F/E_J\ll \tilde\gamma\) when \(\langle \hat n\rangle\sim(\kappa F/\hbar\gamma)^2\), weak quantum fluctuations \(\kappa\ll1\), and near-resonant driving with \( |\delta\omega| \ll \omega_d \) and \(F/E_J<1\). When these conditions hold, the junction is quantitatively described by a driven Kerr oscillator. When they fail, especially for \(\kappa\sim1\) and \(\epsilon\gtrsim2\), the full cosine nonlinearity becomes essential [1812.02839].

The qualitative failure of the Duffing approximation is sharp. In the standard driven softening regime, the Duffing oscillator supports at most bistability, whereas the full Josephson cosine can generate multiple stable dynamical states because its nonlinearity is periodic rather than merely quartic. In the cited study, parameters \(\kappa=1\), \(\epsilon=2.05\), \(f=0.8\), and \(\tilde\gamma=0.05\) produce multiple peaks in the steady-state Fock distribution around \(n=0,4,26,66\) in the full Josephson model, while the corresponding Duffing approximation gives only a single peak near \(n=0\). The same framework shows that quantum fluctuations renormalize the resonance frequency and Josephson energy, and that fluctuations around higher-amplitude fixed points depend nontrivially on Bessel functions \(J_0\) and \(J_2\), which are beyond-Duffing effects [1812.02839].

A persistent misconception is therefore that a strongly driven Josephson junction is generically a Duffing oscillator. The cited analysis shows instead that it is a more general quantum nonlinear Josephson oscillator whose low-amplitude sector is Duffing-like.

## 4. Dissipation, metastability, and dissipative phase transitions

For finite-size driven-dissipative quantum Duffing oscillators, the asymptotic state is unique under general conditions. In Liouvillian language, if \(\delta_n=-{\rm Re}(\lambda_n)\) and the eigenvalues are ordered so that \(\delta_n<\delta_{n+1}\), then \(\delta_0=0\) and \(\delta_1>0\). The smallest nonzero decay rate \(\delta_1\) is the Liouvillian gap and sets the final relaxation time \(1/\delta_1\). When \(\delta_1\ll\delta_2\), the dynamics first relaxes rapidly onto a metastable manifold and only later converges to the unique steady state [2206.06338].

This spectral picture resolves the classical language of bistability and hysteresis. In a superconducting Duffing oscillator realized by a coplanar-waveguide resonator interrupted by a weakly asymmetric SQUID, the two classically identified “steady states” were measured to be long-lived metastable states rather than true stationary states. At \(U/2\pi=-132\ \mathrm{kHz}\), delay-dependent measurements reproduced familiar hysteresis for \(\tau=3.25\ \mu\mathrm{s}\), yet increasing the waiting time made the branches collapse onto a single steady state. Around a critical drive \(\xi^*/2\pi=1.51\ \mathrm{MHz}\), the loop area in the complex reflection coefficient decayed as \(A\propto \exp(-\eta\tau)\) with two rates, \(\eta_1 = 0.74\ \mu\mathrm{s}^{-1}\) and \(\eta_2 = 0.04\ \mu\mathrm{s}^{-1}\), reflecting fast approach to the metastable manifold and slow relaxation to the true steady state. The Liouvillian gap dropped from approximately \(\delta_1\approx 3.79\ \mu\mathrm{s}^{-1}\) away from the transition to \(\delta_1^{\rm min}\approx 0.02\ \mu\mathrm{s}^{-1}\) near criticality, which is the finite-size signature of critical slowing down [2206.06338].

Under the scaling
\[
U \rightarrow \frac{U_0}{N}, \qquad \xi \rightarrow \sqrt{N}\,\xi_0,
\]
with \(U_0\equiv -\gamma\) at \(N=1\) and \(\Delta=3\gamma\), the same device emulates the mean-field description of an \(N\)-site driven-dissipative Bose-Hubbard lattice. The critical rescaled drive \(\xi_0^*/2\pi = 0.58\ \mathrm{MHz}\) is essentially independent of \(N\), while the transition step sharpens as \(N\) increases. Quantum state tomography identified the two sides of the transition as a coherent-state phase and a squeezed-state phase, with a two-component mixed Wigner function across \(0.52\ \mathrm{MHz} \le \xi_0/2\pi \le 0.64\ \mathrm{MHz}\) [2206.06338].

A separate nonequilibrium result concerns global stability in the bistable Duffing regime. For the rotating-frame Lindblad description of the Duffing oscillator, the long-time stationary state can be reduced to an effective two-state model with transition rates \(\Gamma_\uparrow\) and \(\Gamma_\downarrow\) between low- and high-amplitude branches. The crossover point is where \(\Gamma_\uparrow=\Gamma_\downarrow\), but that point is non-unique: the limits \(\hbar\to0\), \(T\to0\), and \(\tau\to\infty\) do not commute. For \(\alpha'=6\), the cited study finds \(\Omega'\approx3.292\) on the classical thermal route and \(\Omega'\approx3.222\) on the quantum zero-temperature route [1208.1701]. This indicates a specifically nonequilibrium notion of tunneling and global stability, distinct from the static double-well criterion of equal potential minima.

## 5. Quantum dynamics near attractors

Near a stable classical state of the driven Duffing oscillator, attractor-centered descriptions are more accurate than bare linearization. Starting from
\[
H=(\delta\omega+\chi)a^\dagger a+\chi (a^\dagger a)^2+\epsilon(a^\dagger+a),
\]
one can displace the oscillator by \(D[\alpha]\) to the attractor center and then squeeze by \(S(\xi)\) to eliminate local anomalous quadratic mixing. The resulting renormalized master equation contains an effective thermal occupation
\[
\bar N=\bar n |v|^2+(1+\bar n)|u|^2,
\]
where \(u\) and \(v\) are the squeezing coefficients. Even at \(\bar n=0\), one has \(\bar N=|u|^2\), so the local oscillator is thermally populated in the squeezed frame by quantum fluctuations associated with the attractor [2502.16803].

This framework shows precisely what standard linearization misses. Around the high-amplitude stable state, the orbit amplitude scales as \(1/\sqrt{\lambda}\), so cubic and quartic fluctuation terms remain parametrically important even when the effective Planck constant \(\lambda\) is small. A refined double perturbation theory for
\[
H=H_0+\gamma H_1+\gamma^2 H_2,\qquad \gamma=\sqrt{\lambda},
\]
yields renormalized local level spacings, corrected eigenvectors, and level-dependent orbital displacements. The local quasienergy spectrum is not exactly harmonic near the attractor bottom: higher-order corrections make \(\Delta E_n=|E_{n+1}-E_n|\) level dependent and decrease the spacing toward the saddle point. Likewise, the average orbital center \(\langle N|a|N\rangle\) becomes level dependent, whereas harmonic linearization would predict \(\langle N|a|N\rangle=\alpha\) for all local levels [2502.16803].

The same treatment gives a balance equation for local populations,
\[
\frac{dp_{n'}}{dt} = \kappa \sum_{m'} \left( W_{n',m'}p_{m'}-W_{m',n'}p_{n'} \right),
\]
valid in the weak-damping regime \(\kappa \ll E_n-E_{n+1}\). In harmonic approximation the adjacent-level ratio is Bose-like,
\[
\frac{p_{n'+1}}{p_{n'}}=\frac{\bar N}{1+\bar N},
\]
but higher-order corrections make the effective occupation number level dependent. Strong damping slightly renormalizes the effective level spacing extracted from the emission spectrum, and pure dephasing \(\eta^{ph}\mathcal D[a^\dagger a]\rho\) raises the effective local temperature through
\[
\tilde N = |u|^2+\bar n(|u|^2+|v|^2)+\frac{\eta^{ph}}{\kappa}|\alpha^*u+\alpha v^*|^2.
\]
Thus dephasing is not only coherence loss; after attractor-dependent displacement and squeezing it also acts as local heating [2502.16803].

## 6. Semiclassical correspondence, escape, and related realizations

In semiclassical open-system studies, quantum-classical differences in the Duffing oscillator are not controlled simply by how chaotic the classical motion is. Using Quantum State Diffusion, a Gaussian-moment semiclassical reduction, and noise-added classical surrogates, one study found that classically regular orbits can show the largest deviations from classical behavior, and that the growth of quantum-classical differences is not determined by the classical Lyapunov exponent. In an intermediate dissipation range, a parametrically invariant meta-attractor emerges at
\[
\beta_{\text{conv}} = 0.0341,
\]
where attractors across different \(\Gamma\) values become nearly indistinguishable in temporal and spatial complexity. Below that scale, noise-added classical models deviate strongly from semiclassical dynamics [2011.00118].

A complementary semiclassical theory treats switching in the weak-dissipation regime as quasienergy diffusion over a dynamical barrier in the rotating frame rather than as coordinate-space escape. For the driven Duffing Hamiltonian
\[
H_S(t)=\frac{p^2}{2M}+\frac{1}{2}M\omega_0^2 q^2-\frac{1}{4}\Gamma q^4+Fq\cos(\omega_d t),
\]
the rotating-frame quasienergy landscape has two stable extrema and an unstable saddle for \(0<\alpha<4/27\). In this regime, finite reflection at the phase-space barrier suppresses switching, but bath-induced quantum fluctuations generated by the rotating-frame \(Q\)- and \(P\)-couplings enhance it, and the cited analysis finds that the environmental enhancement overcompensates the reflection-induced reduction. Near the bifurcation threshold, the barrier softens so strongly that the local dynamics becomes effectively overdamped, and the quantum energy scale associated with friction can exceed the thermal scale [1009.2401].

Experimental and structural generalizations show how the quantum Duffing problem connects to other nonlinear platforms. A single \(^{88}\mathrm{Sr}^{+}\) ion in a linear RF Paul trap realizes a driven nonlinear mechanical oscillator whose steady-state motion is described by
\[
\ddot{x} + 2\mu\dot{x}+\gamma\dot{x}^{3} + \omega_{0}^{2}x + \alpha x^{3} = k \cos(\omega t).
\]
The experiment demonstrated hardening Duffing response, hysteresis, jump phenomena, phase discontinuities, and tunable linear and nonlinear damping through Doppler cooling, but it did not yet realize the deep quantum Duffing regime of nonclassical quasienergy states or quantum tunneling between attractors [1003.1577]. In a different direction, two weakly coupled classical Duffing oscillators in 1:1 resonance reduce after averaging to boson Josephson junction equations for an energy partition index \(P\) and relative phase \(\Delta\), with effective Hamiltonian
\[
H_{\mathrm{eff}} =-\varepsilon\Omega\left(\frac{\kappa}{2}P^2+\cos\Delta\sqrt{1-P^2}\right).
\]
Within that reduced model, local modes correspond to macroscopic quantum self-trapping states, giving a semiclassical bridge between nonlinear vibration localization and Josephson dynamics [1205.4842].

Taken together, these results define the quantum Duffing oscillator as a family of closely related driven-dissipative quantum nonlinear problems rather than a single universal approximation. In low-amplitude weakly anharmonic regimes it is well described by a Kerr Hamiltonian; near resonances it is governed by Floquet quasienergy structure; in driven superconducting circuits it exhibits metastability and finite-size precursors of dissipative phase transitions; near attractors it supports squeezed, effectively thermal local modes with renormalized spectra; and in more general Josephson settings its Duffing description fails once the full periodic nonlinearity of the cosine becomes dynamically relevant.

Source: https://www.emergentmind.com/topics/quantum-duffing-oscillator