Zero-Frequency Ermakov Equation

• Zero-frequency Ermakov equation is the inverse-cubic nonlinear specialization of the Ermakov–Pinney family achieved by eliminating the linear oscillator term.
• It provides exact solutions that underlie free-particle Gaussian wavepacket spreading in quantum mechanics and free-space dynamics in paraxial optics.
• Its significance is underscored by preserved sl(2,ℝ) symmetry and versatile formulations including Schwarzian, Bäcklund, and relativistic extensions.

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 (Tsekov, 2010, Ershkov et al., 2010, Ramos-Prieto et al., 2018, Kim et al., 2016, Morris et al., 2015, Carillo et al., 2022).

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

$z''=Mz+\frac{a}{z^3}.$0

so that the zero-frequency specialization becomes

$z''=Mz+\frac{a}{z^3}.$1

They also give the general Pinney representation

$z''=Mz+\frac{a}{z^3}.$2

with $z''=Mz+\frac{a}{z^3}.$3 solving the corresponding linear equation and the constants constrained by

$z''=Mz+\frac{a}{z^3}.$4

(Kim et al., 2016).

A second common notation is

$z''=Mz+\frac{a}{z^3}.$5

used in a constant-$z''=Mz+\frac{a}{z^3}.$6 analytical study. In that convention, the zero-frequency case corresponds by inference to $z''=Mz+\frac{a}{z^3}.$7, giving

$z''=Mz+\frac{a}{z^3}.$8

(Ershkov et al., 2010).

A third notation appears in KvN mechanics and Lewis–Ermakov invariant constructions: $z''=Mz+\frac{a}{z^3}.$9 so the zero-frequency case is

$$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$0

(Ramos-Prieto et al., 2018). The same zero-frequency specialization follows from the amplitude equation

$$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$1

by setting $$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$2, giving

$$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$3

(Kim et al., 2016).

The sign convention is not completely uniform across the literature. For example, Schwarzian-based work writes the Ermakov equation as

$$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$4

so the frequency-free case is obtained by setting $$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$5: $$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$6 (Carillo et al., 2022). Likewise, work on nonlinear transformations of Liénard systems produces an Ermakov–Pinney instance

$$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$7

whose zero-frequency specialization is

$$\ddot{\rho}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$8

(Guha et al., 2019).

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}=\frac{c}{\rho^3}, \qquad\text{or equivalently}\qquad z''=\frac{a}{z^3}.$$9

the nonlinear solution takes the form

$M=0$0

with the usual determinant constraint on $M=0$1 (Kim et al., 2016).

The zero-frequency specialization makes the associated linear equation trivial: $M=0$2 A convenient independent pair is

$M=0$3

with unit Wronskian. Substituting these into Pinney’s formula gives the exact zero-frequency family

$M=0$4

(Kim et al., 2016). In the $M=0$5 notation, the same family appears as

$M=0$6

or equivalently

$M=0$7

(Ershkov et al., 2010).

This exact solvability is one reason the zero-frequency equation remains central. It is the nonlinear companion of the free linear equation $M=0$8, 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 (Ershkov et al., 2010, Morris et al., 2015).

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,

$M=0$9

the wave function is written as

$M$0

leading to the continuity equation and a force-balance equation with Bohm potential. For a Gaussian density ansatz with rms width $M$1, one obtains

$M$2

This equation governs the width of the packet, not the particle coordinate (Tsekov, 2010).

In this framework, the zero-frequency limit means $M$3, i.e. the free-particle case. The width equation becomes

$M$4

with exact solution

$M$5

(Tsekov, 2010). 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 $M$6 satisfies the general Ermakov equation

$M$7

and free-space propagation corresponds to the zero-frequency case

$M$8

With initial conditions

$M$9

the solution is

$k(t)=0$0

(Huerta-Sandoval et al., 30 Sep 2025). 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 $k(t)=0$1 satisfying

$k(t)=0$2

The zero-frequency specialization is simply

$k(t)=0$3

(Ramos-Prieto et al., 2018). Even though that paper does not analyze $k(t)=0$4 explicitly, the formal construction survives: the invariant, the time-dependent transformations, and the reduced Liouvillian remain meaningful once $k(t)=0$5 solves the zero-frequency equation (Ramos-Prieto et al., 2018).

The same pattern appears in Lewis–Ermakov operator formulations of the time-dependent harmonic oscillator. In one amplitude-phase treatment, the amplitude $k(t)=0$6 satisfies

$k(t)=0$7

The corresponding zero-frequency equation is

$k(t)=0$8

That work stresses that while the invariant can remain finite in the $k(t)=0$9 limit, the number/phase representation becomes singular because formulas such as

$B(z)$0

blow up as $B(z)$1 (Guasti et al., 2013). 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: $B(z)$2 The zero-frequency form

$B(z)$3

has the exact solution

$B(z)$4

(Kim et al., 2016). 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 $B(z)$5 to the Madelung force balance yields

$B(z)$6

At zero frequency this becomes

$B(z)$7

In the strong-friction regime, neglecting $B(z)$8 gives

$B(z)$9

which integrates to

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

and asymptotically yields the sub-diffusive law

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

(Tsekov, 2010).

The same paper introduces temperature-dependent generalizations. One finite-temperature model has

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

so the zero-frequency form is

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

In the high-temperature limit, the thermal equations reduce to

$\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3}.$4

and thus at zero frequency

$\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3}.$5

(Tsekov, 2010).

A radiation-reaction-like extension appends a jerk term,

$\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3}.$6

with zero-frequency specialization

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

(Tsekov, 2010).

Other generalizations change the meaning of the “frequency term” rather than merely suppressing it. In the Schwarzian derivative framework, the proper Ermakov equation is

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

There, the zero-frequency specialization is $\ddot{\rho}+\omega^2(t)\rho=\frac{c}{\rho^3}.$9, yielding

$\omega(t)=0$0

(Carillo et al., 2022). In work on generalized Liénard equations, an Ermakov-Pinney-type equation appears as

$\omega(t)=0$1

so setting $\omega(t)=0$2 gives the zero-frequency version

$\omega(t)=0$3

(Guha et al., 2019).

6. Symmetry, Bäcklund, and geometric perspectives

The zero-frequency Ermakov equation retains substantial algebraic structure. In the autonomous Ermakov–Pinney equation

$\omega(t)=0$4

setting $\omega(t)=0$5 gives

$\omega(t)=0$6

A Lie-symmetry analysis shows that the full three-dimensional $\omega(t)=0$7 symmetry is preserved in this limit. The generators take the projective form

$\omega(t)=0$8

(Morris et al., 2015). This directly contradicts a possible misconception that zero frequency is intrinsically symmetry-breaking. What breaks the full $\omega(t)=0$9 structure in that analysis is time dependence in the nonlinear coefficient $\ddot{\rho}=\frac{c}{\rho^3}.$0, not setting the frequency term to zero (Morris et al., 2015).

Bäcklund and Schwarzian constructions sharpen this picture. The Schwarzian-based treatment writes the proper Ermakov equation as

$\ddot{\rho}=\frac{c}{\rho^3}.$1

with associated Schwarzian equation

$\ddot{\rho}=\frac{c}{\rho^3}.$2

At zero frequency, $\ddot{\rho}=\frac{c}{\rho^3}.$3, hence

$\ddot{\rho}=\frac{c}{\rho^3}.$4

so $\ddot{\rho}=\frac{c}{\rho^3}.$5 is Möbius. The general solution then reduces to

$\ddot{\rho}=\frac{c}{\rho^3}.$6

equivalently the square root of a quadratic polynomial (Carillo et al., 2022). 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

$\ddot{\rho}=\frac{c}{\rho^3}.$7

includes the Ermakov–Pinney equation in the form

$\ddot{\rho}=\frac{c}{\rho^3}.$8

Setting $\ddot{\rho}=\frac{c}{\rho^3}.$9 yields the zero-frequency equation

$z''=Mz+\frac{a}{z^3}.$00

The auto-Bäcklund condition simplifies to vanishing Schwarzian,

$z''=Mz+\frac{a}{z^3}.$01

so the admissible transformations are exactly Möbius maps

$z''=Mz+\frac{a}{z^3}.$02

(Carillo et al., 2017). 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

$z''=Mz+\frac{a}{z^3}.$03

with the spatial coordinate acting as the evolution parameter. The zero-frequency analogue is the special case $z''=Mz+\frac{a}{z^3}.$04, giving

$z''=Mz+\frac{a}{z^3}.$05

That work emphasizes that $z''=Mz+\frac{a}{z^3}.$06 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 (Kumar, 31 Jan 2026).

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 $z''=Mz+\frac{a}{z^3}.$07 crossing zero, the complex mode function $z''=Mz+\frac{a}{z^3}.$08 satisfies the linear equation

$z''=Mz+\frac{a}{z^3}.$09

and its amplitude $z''=Mz+\frac{a}{z^3}.$10 satisfies the Ermakov equation

$z''=Mz+\frac{a}{z^3}.$11

That work does not propose a new nonlinear zero-frequency equation; instead, it analyzes what happens when the standard adiabatic law

$z''=Mz+\frac{a}{z^3}.$12

fails at $z''=Mz+\frac{a}{z^3}.$13 (Dodonov et al., 2023).

The exact energy formula in terms of the Ermakov amplitude is

$z''=Mz+\frac{a}{z^3}.$14

At $z''=Mz+\frac{a}{z^3}.$15, the energy remains finite because it is controlled by $z''=Mz+\frac{a}{z^3}.$16, not by the singular ratio $z''=Mz+\frac{a}{z^3}.$17 (Dodonov et al., 2023). After a slow crossing, the proportionality $z''=Mz+\frac{a}{z^3}.$18 is restored but with a renormalized coefficient determined by mode mixing. For profiles

$z''=Mz+\frac{a}{z^3}.$19

the post-crossing coefficient is

$z''=Mz+\frac{a}{z^3}.$20

In the case $z''=Mz+\frac{a}{z^3}.$21 near zero, corresponding to $z''=Mz+\frac{a}{z^3}.$22, one finds $z''=Mz+\frac{a}{z^3}.$23, so the mean energy triplicates after a single adiabatic crossing (Dodonov et al., 2023).

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

$z''=Mz+\frac{a}{z^3}.$24

so the ordinary zero-frequency limit is

$z''=Mz+\frac{a}{z^3}.$25

The relativistic version reads

$z''=Mz+\frac{a}{z^3}.$26

with

$z''=Mz+\frac{a}{z^3}.$27

so the zero-frequency relativistic specialization is

$z''=Mz+\frac{a}{z^3}.$28

(Haas, 2021). 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 $z''=Mz+\frac{a}{z^3}.$29 (Haas, 2021).

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 $z''=Mz+\frac{a}{z^3}.$30: $z''=Mz+\frac{a}{z^3}.$31 corresponds to harmonic oscillator reduction and yields a genuine Ermakov-type equation for the scaling function $z''=Mz+\frac{a}{z^3}.$32, whereas $z''=Mz+\frac{a}{z^3}.$33 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

$z''=Mz+\frac{a}{z^3}.$34

(Lanfear et al., 2011). 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

$z''=Mz+\frac{a}{z^3}.$35

gives a still nonlinear equation,

$z''=Mz+\frac{a}{z^3}.$36

By contrast, setting the flux or nonlinear constant to zero gives back the linear equation. In the stationary Bohm–Madelung setting, $z''=Mz+\frac{a}{z^3}.$37 eliminates the inverse-cubic term, whereas $z''=Mz+\frac{a}{z^3}.$38 eliminates the effective frequency term (Kumar, 31 Jan 2026). 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,

$z''=Mz+\frac{a}{z^3}.$39

the dynamics encode ballistic broadening of the packet width (Tsekov, 2010).

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 (Ramos-Prieto et al., 2018, Huerta-Sandoval et al., 30 Sep 2025). What may fail is not the invariant itself but number-operator, adiabatic, or oscillator-mode interpretations (Guasti et al., 2013, Dodonov et al., 2023).

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: $z''=Mz+\frac{a}{z^3}.$40 In alternate notations it appears as

$z''=Mz+\frac{a}{z^3}.$41

Its exact general solution is the square root of a quadratic polynomial,

$z''=Mz+\frac{a}{z^3}.$42

or equivalent reparametrizations thereof (Kim et al., 2016, Ershkov et al., 2010, Morris et al., 2015).

Physically, the equation describes free Gaussian wavepacket spreading in quantum mechanics (Tsekov, 2010), free-space scale evolution in paraxial optics (Huerta-Sandoval et al., 30 Sep 2025), and the zero-frequency or free-motion limit of oscillator-based invariant constructions (Ramos-Prieto et al., 2018, Guasti et al., 2013). Structurally, it occupies a privileged position in Pinney superposition theory, $z''=Mz+\frac{a}{z^3}.$43 symmetry analysis, Schwarzian and Bäcklund constructions, and relativistic extensions (Morris et al., 2015, Carillo et al., 2022, Carillo et al., 2017, Haas, 2021).

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 $z''=Mz+\frac{a}{z^3}.$44 as an invariant, not the Ermakov amplitude equation itself (Dodonov et al., 2023). 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.