---
title: Zero-Frequency Ermakov Equation
url: https://www.emergentmind.com/topics/zero-frequency-ermakov-equation
type: topic
---

# Zero-Frequency Ermakov Equation

Searching arXiv for recent and core papers on Ermakov equations and zero-/free-frequency specializations.
First, I’ll confirm the main cited paper and then gather related Ermakov/Pinney references that explicitly or implicitly contain the zero-frequency/free-particle limit.
The zero-frequency Ermakov equation is the frequency-free specialization of the Ermakov–Pinney family, obtained by suppressing the linear restoring term that ordinarily accompanies a harmonic oscillator. In the notation standard across the cited literature, the parent equation is written as
\[
\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3},
\]
or in equivalent sign conventions,
\[
z''=Mz+\frac{a}{z^3}.
\]
The zero-frequency case is the specialization in which the coefficient of the linear term vanishes, so that one obtains
\[
\ddot{\rho}=\frac{c}{\rho^3},
\qquad\text{or equivalently}\qquad
z''=\frac{a}{z^3}.
\]
Although the exact phrase “zero-frequency Ermakov equation” is not uniformly used across the literature, the case is repeatedly present as the free-particle limit of oscillator-based formulations, the \(M=0\) member of constant-\(M\) Ermakov–Pinney systems, the \(k(t)=0\) specialization of invariant-based constructions, or the vanishing-\(B(z)\) specialization of Schwarzian and Bäcklund frameworks [1002.2545], [1012.5374], [1805.08194], [1609.00248], [1510.08992], [2201.02267].

## 1. Definition and notational variants

The zero-frequency Ermakov equation is most directly understood as the case in which the oscillator term disappears from the Ermakov–Pinney equation. One standard form used in the literature is
\[
\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3}.
\]
Setting \(\omega(t)=0\) yields
\[
\ddot{\rho}=\frac{c}{\rho^3}.
\]
Kim and Kim use the equivalent notation
\[
\ddot{x}(t)+\omega^2(t)x(t)=\frac{L^2}{x^3(t)},
\]
so that the zero-frequency specialization becomes
\[
\ddot{x}(t)=\frac{L^2}{x^3(t)}.
\]
They also give the general Pinney representation
\[
x(t)=\big(Au^2(t)+2B\,u(t)v(t)+Cv^2(t)\big)^{1/2},
\]
with \(u,v\) solving the corresponding linear equation and the constants constrained by
\[
AC-B^2=\frac{L^2}{(\mathrm{Wr}[u,v])^2}
\]
[1609.00248].

A second common notation is
\[
\frac{d^2 z}{dx^2}=Mz+\frac{a}{z^3},
\]
used in a constant-\(M\) analytical study. In that convention, the zero-frequency case corresponds by inference to \(M=0\), giving
\[
z''=\frac{a}{z^3}
\]
[1012.5374].

A third notation appears in KvN mechanics and Lewis–Ermakov invariant constructions:
\[
\ddot{\rho}+k(t)\rho=\frac{1}{\rho^3},
\]
so the zero-frequency case is
\[
\ddot{\rho}=\frac{1}{\rho^3}
\]
[1805.08194]. The same zero-frequency specialization follows from the amplitude equation
\[
\ddot{\xi}+\omega^2(t)\xi=\frac{1}{4\xi^3}
\]
by setting \(\omega(t)=0\), giving
\[
\ddot{\xi}=\frac{1}{4\xi^3}
\]
[1609.00248].

The sign convention is not completely uniform across the literature. For example, Schwarzian-based work writes the Ermakov equation as
\[
u''=B(z)u+\frac{I}{u^3},
\]
so the frequency-free case is obtained by setting \(B(z)=0\):
\[
u''=\frac{I}{u^3}
\]
[2201.02267]. Likewise, work on nonlinear transformations of Liénard systems produces an Ermakov–Pinney instance
\[
\ddot{x}+w^2x+\frac{1}{x^3}=0,
\]
whose zero-frequency specialization is
\[
\ddot{x}+\frac{1}{x^3}=0
\]
[1905.00610].

## 2. Historical and structural position within Ermakov–Pinney theory

The zero-frequency equation is not an isolated curiosity but a degenerate member of the general Ermakov–Pinney structure. The classical solution formula attributed to Pinney expresses a nonlinear solution in terms of two independent solutions of the associated linear equation. In the notation
\[
\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3},
\qquad
\ddot{u}+\omega^2(t)u=0,
\qquad
\ddot{v}+\omega^2(t)v=0,
\]
the nonlinear solution takes the form
\[
\rho(t)=\big(Au^2+2Buv+Cv^2\big)^{1/2},
\]
with the usual determinant constraint on \(A,B,C\) [1609.00248].

The zero-frequency specialization makes the associated linear equation trivial:
\[
\ddot{u}=0,\qquad \ddot{v}=0.
\]
A convenient independent pair is
\[
u(t)=1,\qquad v(t)=t,
\]
with unit Wronskian. Substituting these into Pinney’s formula gives the exact zero-frequency family
\[
\rho(t)=\sqrt{A+2Bt+Ct^2},
\qquad AC-B^2=c
\]
[1609.00248]. In the \(z''=Mz+a/z^3\) notation, the same family appears as
\[
z(x)=\pm \sqrt{C(x-x_0)^2+\frac{a}{C}},
\qquad (C\neq 0),
\]
or equivalently
\[
z(x)=\sqrt{Dx^2+2Bx+A},
\qquad AD-B^2=a
\]
[1012.5374].

This exact solvability is one reason the zero-frequency equation remains central. It is the nonlinear companion of the free linear equation \(u''=0\), and it preserves the characteristic square-root-of-a-quadratic structure of the Pinney solution. The equation is therefore a particularly transparent instance of nonlinear superposition, symmetry reduction, and invariant-based reconstruction [1012.5374], [1510.08992].

## 3. Free-particle and wavepacket-width interpretation

One of the clearest physical realizations of the zero-frequency Ermakov equation occurs in quantum mechanics through Gaussian wavepacket dynamics. In a Madelung treatment of the time-dependent Schrödinger equation for a harmonic oscillator,
\[
i\hbar \partial_t \psi = \left(-\frac{\hbar^2}{2m}\partial_x^2+\frac{m\omega_0^2 x^2}{2}\right)\psi,
\]
the wave function is written as
\[
\psi=\sqrt{\rho}\exp(iS/\hbar),
\]
leading to the continuity equation and a force-balance equation with Bohm potential. For a Gaussian density ansatz with rms width \(\sigma(t)\), one obtains
\[
m\ddot{\sigma}+m\omega_0^2 \sigma = \frac{\hbar^2}{4m\sigma^3}.
\]
This equation governs the width of the packet, not the particle coordinate [1002.2545].

In this framework, the zero-frequency limit means \(\omega_0=0\), i.e. the free-particle case. The width equation becomes
\[
m\ddot{\sigma}=\frac{\hbar^2}{4m\sigma^3},
\qquad
\ddot{\sigma}=\frac{\hbar^2}{4m^2\sigma^3},
\]
with exact solution
\[
\sigma^2(t)=\sigma_0^2+\left(\frac{\hbar t}{2m\sigma_0}\right)^2
\]
[1002.2545]. In this context, the zero-frequency Ermakov equation is the equation of ballistic spreading of a free Gaussian quantum packet.

A closely related paraxial-optics construction develops a propagation-dependent invariant formalism for the free-space paraxial wave equation. There the auxiliary scaling function \(\rho(\tau)\) satisfies the general Ermakov equation
\[
\ddot{\rho}+\omega^2(\tau)\rho=\frac{b_0^2}{\rho^3},
\]
and free-space propagation corresponds to the zero-frequency case
\[
\ddot{\rho}=\frac{b_0^2}{\rho^3}.
\]
With initial conditions
\[
\rho(0)=1,\qquad \dot{\rho}(0)=0,
\]
the solution is
\[
\rho(\tau)=\sqrt{1+b_0^2\tau^2}
\]
[2510.00277]. This suggests a direct optical analogue of free-particle width evolution, with the beam-width scale replacing the quantum rms width.

## 4. Invariant formulations and amplitude equations

Invariant-based approaches treat the zero-frequency equation as the auxiliary nonlinear sector underlying exact propagation, rather than as a standalone nonlinear ODE. In KvN mechanics, the time-dependent-frequency harmonic oscillator is controlled by an Ermakov–Lewis invariant built from an auxiliary function \(\rho(t)\) satisfying
\[
\ddot{\rho}+k(t)\rho=\frac{1}{\rho^3}.
\]
The zero-frequency specialization is simply
\[
\ddot{\rho}=\frac{1}{\rho^3}
\]
[1805.08194]. Even though that paper does not analyze \(k(t)=0\) explicitly, the formal construction survives: the invariant, the time-dependent transformations, and the reduced Liouvillian remain meaningful once \(\rho\) solves the zero-frequency equation [1805.08194].

The same pattern appears in Lewis–Ermakov operator formulations of the time-dependent harmonic oscillator. In one amplitude-phase treatment, the amplitude \(\rho\) satisfies
\[
\ddot \rho+\Omega^2(t)\rho=\frac{1}{\rho^3}.
\]
The corresponding zero-frequency equation is
\[
\ddot\rho=\frac{1}{\rho^3}.
\]
That work stresses that while the invariant can remain finite in the \(\omega\to 0\) limit, the number/phase representation becomes singular because formulas such as
\[
\hat n(t)=\frac{1}{G\omega(t)}\,\hat a^\dagger\hat a
\]
blow up as \(\omega(t)\to 0\) [1309.1498]. This clarifies a common misconception: the difficulty at zero frequency is not necessarily the failure of the exact Ermakov equation itself, but the failure of certain adiabatic or number-based interpretations.

A closely related exact amplitude equation appears in Kim and Kim’s invariant construction:
\[
\ddot{\xi}+\omega^2(t)\xi=\frac{1}{4\xi^3}.
\]
The zero-frequency form
\[
\ddot{\xi}=\frac{1}{4\xi^3}
\]
has the exact solution
\[
\xi(t)=\sqrt{A+2Bt+Ct^2},
\qquad AC-B^2=\frac14
\]
[1609.00248]. There, the amplitude determines invariant operators and exact free-particle wavefunctions through a time-dependent width and quadratic phase.

## 5. Dissipative, thermal, and generalized zero-frequency equations

The zero-frequency equation admits several nonconservative generalizations. In the quantum Gaussian-width formulation, adding linear friction \(-bV\) to the Madelung force balance yields
\[
m\ddot{\sigma}+b\dot{\sigma}+m\omega_0^2 \sigma = \frac{\hbar^2}{4m\sigma^3}.
\]
At zero frequency this becomes
\[
m\ddot{\sigma}+b\dot{\sigma}=\frac{\hbar^2}{4m\sigma^3}.
\]
In the strong-friction regime, neglecting \(m\ddot{\sigma}\) gives
\[
b\dot{\sigma}=\frac{\hbar^2}{4m\sigma^3},
\]
which integrates to
\[
\sigma^4(t)=\sigma^4(0)+\frac{\hbar^2}{mb}t,
\]
and asymptotically yields the sub-diffusive law
\[
\sigma^2 = \hbar\sqrt{\frac{t}{mb}}
\]
[1002.2545].

The same paper introduces temperature-dependent generalizations. One finite-temperature model has
\[
m\ddot{\sigma}+b\dot{\sigma}+m\omega_0^2 \sigma
=
2\,\partial_\beta\!\left(\frac{1}{\sigma}\right)_{b}
+\frac{\hbar^2}{4m\sigma^3},
\]
so the zero-frequency form is
\[
m\ddot{\sigma}+b\dot{\sigma}
=
2\,\partial_\beta\!\left(\frac{1}{\sigma}\right)_{b}
+\frac{\hbar^2}{4m\sigma^3}.
\]
In the high-temperature limit, the thermal equations reduce to
\[
m\ddot{\sigma}+b\dot{\sigma}+m\omega_0^2 \sigma
=
\frac{k_BT}{\sigma}+\frac{\hbar^2}{4m\sigma^3},
\]
and thus at zero frequency
\[
m\ddot{\sigma}+b\dot{\sigma}
=
\frac{k_BT}{\sigma}+\frac{\hbar^2}{4m\sigma^3}
\]
[1002.2545].

A radiation-reaction-like extension appends a jerk term,
\[
m\ddot{\sigma}-r\dddot{\sigma}+m\omega_0^2\sigma
=
\frac{\hbar^2}{4m\sigma^3},
\]
with zero-frequency specialization
\[
m\ddot{\sigma}-r\dddot{\sigma}
=
\frac{\hbar^2}{4m\sigma^3}
\]
[1002.2545].

Other generalizations change the meaning of the “frequency term” rather than merely suppressing it. In the Schwarzian derivative framework, the proper Ermakov equation is
\[
u''=B(z)u+\frac{I}{u^3}.
\]
There, the zero-frequency specialization is \(B(z)=0\), yielding
\[
u''=\frac{I}{u^3}
\]
[2201.02267]. In work on generalized Liénard equations, an Ermakov-Pinney-type equation appears as
\[
\ddot{x}
+\left(\frac{\ddot{\alpha}}{\alpha}-2\frac{\dot{\alpha}^2}{\alpha^2}\right)x
=
\frac{\dot{\alpha}^2}{\alpha^2x^3}-w^2x,
\]
so setting \(w=0\) gives the zero-frequency version
\[
\ddot{x}
+\left(\frac{\ddot{\alpha}}{\alpha}-2\frac{\dot{\alpha}^2}{\alpha^2}\right)x
=
\frac{\dot{\alpha}^2}{\alpha^2x^3}
\]
[1905.00610].

## 6. Symmetry, Bäcklund, and geometric perspectives

The zero-frequency Ermakov equation retains substantial algebraic structure. In the autonomous Ermakov–Pinney equation
\[
\ddot{x}+\omega^2x=\frac{h^2}{x^3},
\]
setting \(\omega=0\) gives
\[
\ddot{x}=\frac{h^2}{x^3}.
\]
A Lie-symmetry analysis shows that the full three-dimensional \(sl(2,\mathbb R)\) symmetry is preserved in this limit. The generators take the projective form
\[
\partial_t,\qquad
t\partial_t+\frac12x\partial_x,\qquad
t^2\partial_t+t x\partial_x
\]
[1510.08992]. This directly contradicts a possible misconception that zero frequency is intrinsically symmetry-breaking. What breaks the full \(sl(2,\mathbb R)\) structure in that analysis is time dependence in the nonlinear coefficient \(G(t)\), not setting the frequency term to zero [1510.08992].

Bäcklund and Schwarzian constructions sharpen this picture. The Schwarzian-based treatment writes the proper Ermakov equation as
\[
u''=B(z)u+\frac{I}{u^3},
\]
with associated Schwarzian equation
\[
\{\Omega,z\}=-2B(z).
\]
At zero frequency, \(B(z)=0\), hence
\[
\{\Omega,z\}=0,
\]
so \(\Omega\) is Möbius. The general solution then reduces to
\[
u(z)^2 = 2r\,\frac{(az+b)(cz+d)}{ad-bc},
\qquad ad-bc\neq 0,
\]
equivalently the square root of a quadratic polynomial [2201.02267]. This shows that the zero-frequency equation is the Möbius-invariant sector of the Schwarzian construction.

A related Bäcklund analysis of the class
\[
y\,y''=F(z,y^2)
\]
includes the Ermakov–Pinney equation in the form
\[
y''=Q(z)y-\frac{\alpha}{y^3}.
\]
Setting \(Q(z)=0\) yields the zero-frequency equation
\[
y''=-\frac{\alpha}{y^3}.
\]
The auto-Bäcklund condition simplifies to vanishing Schwarzian,
\[
\{f,z\}=0,
\]
so the admissible transformations are exactly Möbius maps
\[
f(z)=\frac{az+b}{cz+d},
\qquad ad-bc\neq 0
\]
[1711.04304]. This suggests that zero frequency is not merely a degenerate limit but a particularly transparent projective sector of the broader Ermakov family.

A stationary Bohm–Madelung reformulation of separable quantum mechanics provides a different geometric perspective. There the sector amplitudes satisfy
\[
\rho_i''(q_i)+\Omega_i^2(q_i)\rho_i(q_i)=\frac{k_i}{\rho_i^3(q_i)},
\]
with the spatial coordinate acting as the evolution parameter. The zero-frequency analogue is the special case \(\Omega_i^2(q_i)=0\), giving
\[
\rho_i''(q_i)=\frac{k_i}{\rho_i^3(q_i)}.
\]
That work emphasizes that \(\Omega_i^2\) includes both physical and geometric contributions after Liouville normalization, so “zero frequency” means cancellation of the full normal-form coefficient, not merely absence of an external oscillator potential [2602.00507].

## 7. Adiabatic zero crossings and limits of the usual invariant picture

The most delicate modern use of the zero-frequency limit concerns systems whose frequency passes through zero. In the study of slowly varying oscillators with \(\omega(t)\) crossing zero, the complex mode function \(\varepsilon(t)\) satisfies the linear equation
\[
\ddot{\varepsilon}+\omega^2(t)\varepsilon=0,
\]
and its amplitude \(\rho(t)=|\varepsilon(t)|\) satisfies the Ermakov equation
\[
\ddot\rho+\omega^2(t)\rho=\rho^{-3}.
\]
That work does not propose a new nonlinear zero-frequency equation; instead, it analyzes what happens when the standard adiabatic law
\[
\frac{\mathcal E(t)}{\omega(t)}=\text{const}
\]
fails at \(\omega=0\) [2303.08299].

The exact energy formula in terms of the Ermakov amplitude is
\[
\langle {\cal E}\rangle_t =
\frac{\langle {\cal E}\rangle_{-\tau}}{2\omega_0}
\left\{
\omega^2(t)\rho^2(t)+[\dot\rho(t)]^2+[\rho(t)]^{-2}
\right\}.
\]
At \(\omega=0\), the energy remains finite because it is controlled by \(\dot\rho^2+\rho^{-2}\), not by the singular ratio \(\mathcal E/\omega\) [2303.08299]. After a slow crossing, the proportionality \(\mathcal E\propto \omega\) is restored but with a renormalized coefficient determined by mode mixing. For profiles
\[
\omega^2(t)\sim |t|^n,
\qquad
\nu=\frac{1}{n+2},
\]
the post-crossing coefficient is
\[
\beta=\frac{1+\cos^2(\nu\pi)}{\sin^2(\nu\pi)}.
\]
In the case \(\omega(t)\sim |t|\) near zero, corresponding to \(n=2\), one finds \(\beta=3\), so the mean energy triplicates after a single adiabatic crossing [2303.08299].

This suggests that the zero-frequency Ermakov regime is less a failure of exact invariant theory than a failure of naive adiabatic interpretation. The nonlinear amplitude equation remains valid; what changes is the asymptotic branch structure of its solutions.

## 8. Related free and relativistic extensions

The zero-frequency equation also appears in relativistic and field-theoretic generalizations. A relativistic Ermakov–Milne–Pinney system derived from planar relativistic oscillator motion has nonrelativistic form
\[
\ddot{\rho}+\kappa^2(t)\rho=\frac{J^2}{\rho^3},
\]
so the ordinary zero-frequency limit is
\[
\ddot{\rho}=\frac{J^2}{\rho^3}.
\]
The relativistic version reads
\[
\ddot{\rho}+\frac{\kappa^2}{\gamma}\left(1-\frac{\dot\rho^2}{c^2}\right)\rho
=
\frac{J^2}{\gamma^2\rho^3},
\]
with
\[
\gamma=
\left(
\frac{1+J^2/(c^2\rho^2)}{1-\dot\rho^2/c^2}
\right)^{1/2},
\]
so the zero-frequency relativistic specialization is
\[
\ddot{\rho}=\frac{J^2}{\gamma^2\rho^3}
\]
[2102.09613]. There the ordinary inverse-cubic centrifugal barrier is dressed by Lorentz-factor nonlinearities. The associated relativistic Ermakov–Lewis invariant remains meaningful because it does not depend explicitly on \(\kappa\) [2102.09613].

A different sort of degeneration appears in exact wave-function constructions for generalized quadratic Hamiltonians. There the autonomous target problem is selected by a discrete parameter \(c_0\): \(c_0=1\) corresponds to harmonic oscillator reduction and yields a genuine Ermakov-type equation for the scaling function \(\mu(t)\), whereas \(c_0=0\) corresponds to free-particle reduction and the inverse-cubic term disappears. The “zero-frequency” case in that framework is therefore not another nonlinear EP equation but the degenerate branch in which the Ermakov structure collapses to a linear characteristic equation
\[
\mu''-\tau(t)\mu'+4\sigma(t)\mu=0
\]
[1102.5119]. This suggests that “zero frequency” can mean either retention of the inverse-cubic equation with vanishing linear term, or degeneration of an Ermakov-type reduction to its underlying linear sector, depending on how the auxiliary variables are defined.

## 9. Misconceptions and interpretive cautions

A recurring source of confusion is the identification of “zero frequency” with “zero nonlinearity.” The literature distinguishes these clearly. Setting the linear coefficient to zero in
\[
\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3}
\]
gives a still nonlinear equation,
\[
\ddot{\rho}=\frac{c}{\rho^3}.
\]
By contrast, setting the flux or nonlinear constant to zero gives back the linear equation. In the stationary Bohm–Madelung setting, \(C_i=0\) eliminates the inverse-cubic term, whereas \(\Omega_i^2=0\) eliminates the effective frequency term [2602.00507]. These are distinct limits.

A second misconception is that the zero-frequency equation is dynamically trivial because the corresponding linear equation is free. The nonlinear sector is not trivial: its exact solution still contains nontrivial width evolution, invariant structure, and projective symmetry. In free-particle quantum spreading,
\[
m\ddot{\sigma}=\frac{\hbar^2}{4m\sigma^3},
\]
the dynamics encode ballistic broadening of the packet width [1002.2545].

A third misconception is that zero frequency necessarily destroys the invariant formalism. Several papers indicate the opposite. The KvN and paraxial-optics constructions retain the invariant structure after setting the frequency coefficient to zero, provided the auxiliary function solves the zero-frequency Ermakov equation [1805.08194], [2510.00277]. What may fail is not the invariant itself but number-operator, adiabatic, or oscillator-mode interpretations [1309.1498], [2303.08299].

## 10. Summary

The zero-frequency Ermakov equation is the inverse-cubic nonlinear equation obtained by removing the linear oscillator term from the Ermakov–Pinney family:
\[
\ddot{\rho}=\frac{c}{\rho^3}.
\]
In alternate notations it appears as
\[
z''=\frac{a}{z^3},\qquad
u''=\frac{I}{u^3},\qquad
m\ddot{\sigma}=\frac{\hbar^2}{4m\sigma^3}.
\]
Its exact general solution is the square root of a quadratic polynomial,
\[
\rho(t)=\sqrt{A+2Bt+Ct^2},
\qquad AC-B^2=c,
\]
or equivalent reparametrizations thereof [1609.00248], [1012.5374], [1510.08992].

Physically, the equation describes free Gaussian wavepacket spreading in quantum mechanics [1002.2545], free-space scale evolution in paraxial optics [2510.00277], and the zero-frequency or free-motion limit of oscillator-based invariant constructions [1805.08194], [1309.1498]. Structurally, it occupies a privileged position in Pinney superposition theory, \(sl(2,\mathbb R)\) symmetry analysis, Schwarzian and Bäcklund constructions, and relativistic extensions [1510.08992], [2201.02267], [1711.04304], [2102.09613].

The modern literature also shows that the exact equation remains well defined when adiabatic frequency profiles pass through zero; what fails in that regime is the naive interpretation of \(\mathcal E/\omega\) as an invariant, not the Ermakov amplitude equation itself [2303.08299]. For that reason, the zero-frequency Ermakov equation serves both as the simplest nontrivial member of the Ermakov–Pinney hierarchy and as a diagnostic limit in which the relation between invariant structure, adiabaticity, and free evolution becomes especially transparent.

Source: https://www.emergentmind.com/topics/zero-frequency-ermakov-equation