Lewis–Ermakov Invariant in Oscillator Dynamics

• The Lewis–Ermakov invariant is a conserved quantity derived from the Ermakov equation that governs time-dependent harmonic oscillator dynamics in both classical and quantum frameworks.
• It is constructed using a positive auxiliary function ρ(t) that satisfies a nonlinear differential equation, ensuring exact conservation despite explicit time dependence.
• Its applications extend to quantum ladder operators, optics, cosmology, and invariant-based control, highlighting its broad methodological and practical significance in physics.

The Lewis–Ermakov invariant, also called the Ermakov–Lewis invariant or Lewis invariant, is a conserved quantity associated with nonautonomous oscillator dynamics, most prominently the time-dependent harmonic oscillator (TDHO). In its standard unit-mass form, for

$H(t)=\frac{p^2}{2}+\frac{1}{2}\omega^2(t)x^2,$

one introduces a positive auxiliary function $\rho(t)$ satisfying the Ermakov equation

$$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$

and obtains the invariant

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$

Its defining feature is exact conservation despite explicit time dependence of the Hamiltonian, and this property underlies a broad body of work in classical mechanics, quantum mechanics, optics, cosmology, Koopman–von Neumann mechanics, and modern invariant-based control (Sen et al., 2020).

1. Historical setting and standard definition

The modern theory combines two strands. Ermakov’s nineteenth-century nonlinear auxiliary equation provided the amplitude variable $\rho$, while Lewis and Lewis–Riesenfeld established the invariant method for the quantum TDHO. In the standard one-dimensional case, $x(t)$ solves

$\ddot{x}(t)+\omega^2(t)x(t)=0,$

and $\rho(t)$ solves the associated Ermakov equation; the Lewis–Ermakov invariant is then the quadratic form written above. In the quantum theory the same structure becomes a Hermitian operator $\hat I_{\mathrm{EL}}$ satisfying

$$\frac{d\hat I_{\mathrm{EL}}}{dt}-i[\hat I_{\mathrm{EL}},\hat H(t)]=0,$$

so its eigenvalues are time-independent even though $\rho(t)$0 is not (Sen et al., 2020).

A commonly used generalization includes a time-dependent mass. For

$\rho(t)$1

the invariant is written as

$\rho(t)$2

with $\rho(t)$3 now determined by

$\rho(t)$4

This form is the natural extension for variable-mass quadratic systems and appears directly in optical analog models of time-dependent oscillators (Lara et al., 2014).

The standard invariant is exact rather than adiabatic. Several sources in the literature emphasize that $\rho(t)$5 is closely related to the adiabatic quantity $\rho(t)$6 for slowly varying frequency, but unlike the adiabatic invariant it remains conserved for arbitrary time dependence of $\rho(t)$7 once $\rho(t)$8 satisfies the appropriate Ermakov equation (Ramos-Prieto et al., 2018).

2. Classical constructions and equivalent formulations

A canonical construction uses two real linearly independent solutions $\rho(t)$9 of the linear TDHO equation and their Wronskian

$$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$0

Writing

$$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$1

one obtains

$$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$2

This orthogonal-functions representation identifies $$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$3 as an amplitude and $$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$4 as a phase, with $$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$5. In this language the classical invariant can be expressed not only in the usual $$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$6 form but also as an amplitude-phase quantity,

$$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$7

which makes explicit the interplay between envelope and phase rotation (Guasti et al., 2013).

Another influential formulation interprets the invariant as ordinary conserved energy after a time-dependent change of variables. With

$$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$8

and $$\ddot{\rho}(t)+\omega^2(t)\rho(t)=\frac{1}{\rho^3(t)},$$9 chosen to satisfy the appropriate Ermakov-type equation, the original TDHO is mapped to a constant-frequency oscillator in the new variables. The conserved energy of that transformed system becomes, when re-expressed in the original variables, the Lewis–Ermakov invariant. This viewpoint was used to “demystify” the constancy of the invariant and to connect it to a broader family of Ray–Reid-type invariants generated by more general transformed potentials (Padmanabhan, 2017).

The same classical framework also extends beyond the standard Pinney nonlinearity. For Reid’s $$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$0-th order generalized Ermakov systems, the nonlinear equation is equivalent to an integrable Emden–Fowler equation, and a closed formula for the higher-order invariant can be written. In the standard $$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$1 case the invariant is

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$2

while for $$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$3 the invariant acquires the characteristic higher-power term

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$4

This preserves the basic “Wronskian-like term plus nonlinear potential term” architecture of the Lewis–Ermakov construction (Mancas et al., 2014).

3. Quantum operator theory, ladder structures, and exact solution methods

In the quantum TDHO, the invariant is not merely conserved in expectation value; it organizes the exact solution space. A standard operator realization is

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$5

and one may define invariant and time-dependent ladder operators by

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$6

with

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$7

The two ladder sets are related by a unitary phase transformation

$$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$8

so the explicit time dependence is a phase rotation generated by the invariant itself (Guasti et al., 2013).

This operator structure supports a number-phase representation. A generalized Turski phase operator $$I_{\mathrm{EL}}=\frac{1}{2}\left[\left(\rho p-\dot{\rho}x\right)^2+\left(\frac{x}{\rho}\right)^2\right].$$9 satisfies

$\rho$0

and the coordinate operator can be written in amplitude-phase form using Dirac-type relations

$\rho$1

Within this formulation the invariant constrains the co-evolution of excitation number and instantaneous frequency. In particular, for each mode the paper on Manley–Rowe relations states

$\rho$2

which becomes the basis for a quantum derivation of power-balance relations in nonlinear optics (Guasti et al., 2013).

A distinct quantum reinterpretation identifies the invariant with second-order quantum moments. For Gaussian wave packets in a time-dependent harmonic well, the width $\rho$3 obeys an Ermakov equation,

$\rho$4

and the quadratic uncertainty Casimir

$\rho$5

is conserved. For Gaussian states this Casimir reproduces the Lewis–Ermakov invariant up to fixed factors and a constant term. The same work further interprets the wave-packet width as an intrinsic clock through

$\rho$6

which maps the mean trajectory to a constant-frequency oscillator in the $\rho$7 variable (Livine, 2022).

4. Interpretations: angular momentum, symmetry, and hidden geometry

One classical interpretation, due to Eliezer and Gray and reused in later work, embeds the problem into an auxiliary two-dimensional isotropic oscillator. Writing $\rho$8, $\rho$9, one has conserved angular momentum

$x(t)$0

and the Ermakov–Milne–Pinney equation becomes the radial equation

$x(t)$1

The Lewis–Ermakov invariant then reduces to

$x(t)$2

so, in this picture, the invariant is simply angular momentum squared over two (Haas, 2021).

That same angular-momentum interpretation survives in a relativistic extension. For the relativistic Ermakov–Milne–Pinney system, the invariant becomes

$x(t)$3

and again evaluates to $x(t)$4. The novelty lies not in the conserved quantity’s physical meaning but in the modified radial equation, where relativistic factors introduce additional nonlinearities in both the restoring term and the inverse-cubic term (Haas, 2021).

A complementary symmetry-based interpretation uses time reparameterizations and the Virasoro action on the TDHO. Under the synchronizing change of variables $x(t)$5 and $x(t)$6, with $x(t)$7 solving the Ermakov equation, the dynamics becomes that of a constant-frequency oscillator in $x(t)$8. The invariant is then the oscillator energy in the synchronized frame, while the underlying quadratic moments generate an $x(t)$9 or $\ddot{x}(t)+\omega^2(t)x(t)=0,$0-type algebra. This places the Lewis–Ermakov invariant among the conserved charges associated with the hidden conformal structure of the problem (Livine, 2022).

A further geometric interpretation appears in stationary Bohm–Madelung quantum mechanics. After Liouville normalization of separated Sturm–Liouville sectors, each sector obeys

$\ddot{x}(t)+\omega^2(t)x(t)=0,$1

while the Bohm amplitude satisfies

$\ddot{x}(t)+\omega^2(t)x(t)=0,$2

The sectorwise invariant

$\ddot{x}(t)+\omega^2(t)x(t)=0,$3

is then independent of the spatial coordinate $\ddot{x}(t)+\omega^2(t)x(t)=0,$4. In this formulation the quantum potential is encoded as a curvature contribution of the self-adjoint Sturm–Liouville operator rather than as an additional dynamical term (Kumar, 31 Jan 2026).

5. Reformulations and generalizations beyond the standard oscillator

A notable nonstandard setting is Koopman–von Neumann (KvN) mechanics, the Hilbert-space formulation of classical mechanics. For the TDHO the KvN Hamiltonian can be written

$\ddot{x}(t)+\omega^2(t)x(t)=0,$5

and, after a canonical transformation in the extended KvN phase space, it splits as

$\ddot{x}(t)+\omega^2(t)x(t)=0,$6

with each $\ddot{x}(t)+\omega^2(t)x(t)=0,$7 an ordinary quantum TDHO Hamiltonian. This permits a direct import of two standard Ermakov–Lewis invariants and yields the KvN invariant

$\ddot{x}(t)+\omega^2(t)x(t)=0,$8

An earlier KvN construction derived the invariant from a quadratic ansatz and a coupled ODE system; the later simplification avoided those coupled equations but was explicitly described as less general because it relied on the Hamiltonian splitting into two standard TDHOs (Ramos-Prieto et al., 2018, Sen et al., 2020).

The invariant also survives nontrivially under deformations of the underlying oscillator. For singular parametric oscillators generated by one-parameter Darboux deformations of the classical harmonic oscillator, the associated Lewis–Ermakov invariants were shown to be independent of the Darboux deformation parameter $\ddot{x}(t)+\omega^2(t)x(t)=0,$9 and singularity-free, even though the deformed frequencies and linear modes possess finite-time singularities. In that setting the invariant retains the undeformed harmonic-oscillator form

$\rho(t)$0

for both the odd and even Darboux families (Rosu et al., 2024).

An algebraic extension reaches time-dependent two-level systems. Exploiting the shared complexification $\rho(t)$1 of the $\rho(t)$2 algebra of two-level systems and the $\rho(t)$3 algebra of quadratic oscillators, a Lewis–Ermakov-type invariant was constructed for

$\rho(t)$4

The auxiliary scaling function $\rho(t)$5 satisfies a complex Ermakov–Pinney equation, and the invariant $\rho(t)$6 generates a closed-form evolution operator for arbitrary time-dependent coupling $\rho(t)$7 and detuning $\rho(t)$8. The same framework was used for Landau–Zener transitions, non-Hermitian dissipative processes, adiabatic rapid passage, and inverse-engineered shortcuts to adiabaticity (Huerta-Sandoval et al., 21 Jun 2026).

6. Applications, methodological significance, and limitations

In optics the invariant appears in several distinct guises. In nonlinear optics, writing the quantum invariant in terms of number and phase operators yields modewise relations of the form $\rho(t)$9, which were used to derive Manley–Rowe-type power-balance relations from a quantum invariant rather than from semiclassical transport equations (Guasti et al., 2013). In photonic lattices, a class of waveguide arrays was mapped to a quantum harmonic oscillator with propagation-dependent mass and frequency, so the lattice inherits Lewis–Ermakov symmetry through the invariant

$\hat I_{\mathrm{EL}}$0

with $\hat I_{\mathrm{EL}}$1 satisfying the generalized Ermakov equation (Lara et al., 2014).

A spatial version of the invariant was developed for quantum light propagation in separable longitudinally inhomogeneous media. After generalized canonical transformations to a comoving frame, the generator of propagation becomes a physical Lewis–Ermakov invariant in space,

$\hat I_{\mathrm{EL}}$2

and the net propagation effect is a quantum Gouy phase rather than the unphysical squeezing produced by naive quantization in the original variables (Barral et al., 2015). A related paraxial-wave treatment showed that finite-energy free-space beams behave as if subject to an effective quadratic GRIN-like confinement, with a zero-frequency Ermakov equation

$\hat I_{\mathrm{EL}}$3

and invariant

$\hat I_{\mathrm{EL}}$4

That work also stressed a limitation: for any $\hat I_{\mathrm{EL}}$5, $\hat I_{\mathrm{EL}}$6 and the free-space Hamiltonian do not commute, and exact commutativity would require an infinitely wide, nonsquare-integrable field (Huerta-Sandoval et al., 30 Sep 2025).

Cosmological and gravitational applications arise because many mode equations and scale-factor equations reduce to EMP form. Within an Eisenhart–Duval lift, the auxiliary function $\hat I_{\mathrm{EL}}$7 obeys

$\hat I_{\mathrm{EL}}$8

and the usual invariant

$\hat I_{\mathrm{EL}}$9

emerges as the conserved energy of a conformally related static oscillator metric. In the cosmological extension of that construction, the scale factor itself satisfies an EMP equation,

$$\frac{d\hat I_{\mathrm{EL}}}{dt}-i[\hat I_{\mathrm{EL}},\hat H(t)]=0,$$0

linking the Ermakov structure to Friedmann-type dynamics and Dmitriev–Zel'dovich equations (Cariglia et al., 2018).

Several methodological caveats recur across the literature. In the stochastic setting, numerical experiments with multiplicative noise suggested that the classical invariant is not affected when the same Wiener process is added to both members of the Ermakov system, but the same study explicitly noted the absence of a rigorous stochastic proof and found strong deviations when the noise was applied asymmetrically (Cervantes-Lopez et al., 2014). In virialised stellar dynamics, the Ermakov–Lewis–Leach invariant was expanded asymptotically in small parameters to generate an infinite series of conservation laws; the resulting constants were stated to depend on the virialised radius and virialised mass independently of the potential used, but the construction is explicitly asymptotic and tied to short-time regimes (Lecian et al., 2023). These examples show that the invariant is both robust and structurally adaptable, but also that the form of the auxiliary equation, the choice of variables, and the admissible class of transformations determine how far any given construction extends.