Ermakov Invariants in Stationary Quantum Mechanics: A Bohm-Madelung and Hamilton-Jacobi Perspective

1. Introduction

The stationary state description of Quantum Mechanics developed through de Broglie Bohm’s method is introduced as a simple application in Peter Holland’s, Quantum theory of Motion, Chapter 4, Sec 4.1 [1]. The author describes the general properties of stationary state and its correspondence to WKB’s version and concludes with an approximate solution in one dimension, where in reference to Ermakov’s differential equation. There had been very little extension of the solution in the rest of the monograph, however, the solution is a first integrable invariant and has interesting applications to description of Bohmian’s Hamilton–Jacobi (HJ) equations. The present work generalizes the role of EP invariants in stationary quantum mechanics by showing that they arise naturally within the BM formulation for diagonal and separable Hamiltonians, while preserving the standard probabilistic structure of quantum mechanics.

The EP equation and its associated invariant originate in the nineteenth-century work of Ermakov [2] and were rediscovered and systematized by Lewis and Riesenfeld in the context of time-dependent harmonic oscillators [3]. Subsequent mathematical developments clarified the nonlinear superposition principle and invariant structure of Ermakov systems, notably in the work of Reid and collaborators [4,5], including the concise analysis of Reid and Ray on nonlinear oscillator invariants [6].

Connections between Ermakov-type invariants and quantum mechanics have appeared sporadically, often within trajectory-based or hydrodynamic formulations. In particular, Nassar developed an invariant formulation of quantum systems within the Bohmian and quantum potential framework, emphasizing auxiliary equations closely related to the EP system [7,8]. Despite these contributions, explicit treatments of Ermakov invariants in stationary BM equations and spatially separable Hamiltonians remain limited.

The present work addresses this gap by demonstrating that Ermakov invariants arise naturally in stationary quantum mechanics whenever the Hamiltonian is diagonal and separable, with the EP equation emerging as the amplitude-space counterpart of HJ separability.

• The Ermakov–Pinney (EP) equation

  $${\rho }^{″}+{\Omega }^2\left(x\right)\rho ={\displaystyle \frac{k}{{\rho }^3}}$$

  (1)

  plays a central role in the theory of time-dependent oscillators [9,10] and admits a conserved Ermakov–Lewis invariant (or "invariant") [11]. In contrast, stationary quantum mechanics is usually presented as an eigenvalue problem, with little emphasis on hidden invariants beyond the energy.

Ermakov–Lewis invariant structures in quantum mechanics have been discussed in earlier works, particularly by Reinisch [12] and Schuch [13,14] in Bohm–Madelung and Riccati formulations. The present work extends the formalism to separable coordinate systems and demonstrates its application suitable for a wide class of boundary value problems.

In this work we demonstrate that this distinction is largely historical rather than structural. When stationary quantum systems are expressed in the BM representation, the same mathematical structure underlying Ermakov invariants are identical to that seen with time, while the spatial coordinate playing the role of evolution parameter. This occurs whenever the Hamiltonian is diagonal and separable.

• Scope of the present work.

The analysis presented here is restricted to stationary quantum systems with diagonal and separable Hamiltonians, for which the BM equations admit a coordinate-wise decomposition and a natural EP structure at the level of the amplitude dynamics. Issues related to angular-variable subtleties, non-separable Hamiltonians, and broken invariance due to non-diagonal (or gauge-coupled) kinetic terms are intentionally not pursued in this manuscript.

Given the broad scope of the present work, spanning coordinate-separable structures, and nonlinear representations of quantum mechanics, it is possible that some related contributions have not been cited. Any such omissions are unintentional, and the authors welcome further references to relevant literature.

2. Stationary BM Formulation and EP Structure

Consider the stationary Schrödinger equation

$$\widehat{H}\psi =E\psi ,$$

(2)

with the wave function $\psi$ in polar decomposition

$$\psi (q_1,\dots ,q_n)=R(q_1,\dots ,q_n)e^{iS(q_1,\dots ,q_n)/\hslash }.$$

(3)

We assume a diagonal Hamiltonian of the form

$$H=\sum _{i=1}^n{\displaystyle \frac{p_i^2}{2m}}+V(q_1,\dots ,q_n),$$

(4)

and separability,

$$R=\prod _i{R}_i\left(q_i\right),S=\sum _i{S}_i\left(q_i\right).$$

(5)

Substitution into the Schrödinger equation yields a set of continuity equations and separated energy equations for each degree of freedom.

2.1. Continuity Equation and Ermakov Invariant Connection

The separated continuity equations take the form

$${\displaystyle \frac{d}{d{q}_i}}\left(R_i^2{p}_i\right)=0,$$

(6)

which integrate immediately to

$$p_i={\displaystyle \frac{C_i}{R_i^2}},$$

(7)

where $C_i$ are constants subject only to global conservation constraints (for example, ${\sum }_i{C}_i=0$ for stationary states).

2.2. The EP Equation and the Ermakov-Lewis Invariant

The separated energy equations are

$${\displaystyle \frac{p_i^2}{2m}}+V_i\left(q_i\right)+Q_i=E_i,\sum _i{E}_i=E,$$

(8)

with quantum potential

$$Q_i=-{\displaystyle \frac{{\hslash }^2}{2m}}{\displaystyle \frac{R_i^{″}}{R_i}}.$$

(9)

Substituting Eq. (7) and scaling the coefficients yields

$$R_i^{″}+{\Omega }_i^2\left(q_i\right)R_i={\displaystyle \frac{k_i}{R_i^3}},$$

(10)

which is precisely the EP equation for each coordinate $q_i$.

Thus, any diagonal and separable stationary quantum system admits EP dynamics in amplitude space.

The nonlinear term proportional to $1/R^3$ arises directly from the stationary Bohm continuity constraint and distinguishes the EP equation from its linear counterpart. In addition to the differential equation itself, the EP system admits a first integral, commonly referred to as the Ermakov–Lewis invariant.

Specifically, if $y\left(x\right)$ is any solution of the associated linear equation

$$y^{″}\left(x\right)+{\Omega }^2\left(x\right)y\left(x\right)=0,$$

(11)

and $R\left(x\right)$ satisfies the EP equation,

$$R^{″}\left(x\right)+{\Omega }^2\left(x\right)R\left(x\right)={\displaystyle \frac{k}{R^3\left(x\right)}},$$

(12)

then the quantity

$$\begin{array}{|c|}\hline I={\displaystyle \frac{1}{2}}\left[{\left(R{y}^{\prime }-R^{\prime }y\right)}^2+{\displaystyle \frac{k{y}^2}{R^2}}\right]\\ \hline\end{array}$$

(13)

is independent of the coordinate x. Equation (13) therefore provides a first invariant of the stationary amplitude dynamics.

Equation (7) is the key link between stationary quantum mechanics and Ermakov theory.

3. Canonical One-Dimensional Examples

The following examples illustrate how the EP structure manifests for the three canonical one-dimensional stationary quantum systems: the free particle, the harmonic oscillator, and the Coulomb potential.

3.1. Free Particle: Bohm Continuity ⇒ Ermakov Amplitude and Invariant

For the free particle $V\left(x\right)=0$ and $E={\hslash }^2{k}_0^2/\left(2m\right)$, hence

$${\Omega }^2\left(x\right)=k_0^2.$$

(14)

The Bohmian amplitude therefore satisfies the constant-frequency EP equation

$$R^{″}+k_0^{2}R={\displaystyle \frac{k}{R^3}},k={\displaystyle \frac{C^2}{{\hslash }^2}}.$$

(15)

A particularly transparent solution is the constant-amplitude solution $R\left(x\right)=R_0$, for which (15) gives

$$k_0^2{R}_0={\displaystyle \frac{k}{R_0^3}}\Rightarrow R_0^4={\displaystyle \frac{k}{k_0^2}}={\displaystyle \frac{C^2}{{\hslash }^2{k}_0^2}}\Rightarrow R_0={\left({\displaystyle \frac{\left|C\right|}{\hslash k_0}}\right)}^{1/2}.$$

(16)

The Bohmian momentum $p=S^{\prime }$ then becomes

$$p\left(x\right)={\displaystyle \frac{C}{R_0^2}}=\mathrm{sgn}\left(C\right)\hslash k_0,$$

(17)

recovering the constant momentum of the plane wave. The corresponding probability density is ${\left|\psi \right|}^2=R_0^2=\mathrm{const}$, consistent with standard quantum mechanics for scattering states.

To connect explicitly with the linear partner equation, note that the associated linear equation $y^{″}+k_0^{2}y=0$ admits the independent solutions

$$y_1\left(x\right)=cos\left(k_{0}x\right),y_2\left(x\right)=sin\left(k_{0}x\right),W=y_1{y}_2^{\prime }-y_1^{\prime }{y}_2=k_0.$$

(18)

Hence the general Ermakov amplitude can be written as

$$R^2\left(x\right)=A{cos}^2\left(k_{0}x\right)+B{sin}^2\left(k_{0}x\right)+2Dsin\left(k_{0}x\right)cos\left(k_{0}x\right),AB-D^2={\displaystyle \frac{k}{k_0^2}}.$$

(19)

The constant-amplitude solution corresponds to the symmetric choice $A=B=R_0^2$ and $D=0$, which satisfies $AB-D^2=R_0^4=k/k_0^2$.

Finally, for any linear solution y and the Ermakov solution R, the invariant (13) is constant in x.

3.2. Harmonic Oscillator: Weber Basis and Bohmian Ermakov Amplitude

For the one-dimensional harmonic oscillator with $V\left(x\right)={\displaystyle \frac{1}{2}}m{\omega }^2{x}^2,$ the stationary BM reduction yields the EP equation for the amplitude $R\left(x\right)$,

$$R^{″}\left(x\right)+\left({\displaystyle \frac{2mE}{{\hslash }^2}}-{\displaystyle \frac{m^2{\omega }^2}{{\hslash }^2}}{x}^2\right)R\left(x\right)={\displaystyle \frac{C^2}{{\hslash }^2}}{\displaystyle \frac{1}{R^3\left(x\right)}}.$$

(20)

The associated linear partner equation is obtained by setting $C=0$:

$$y^{″}\left(x\right)+\left({\displaystyle \frac{2mE}{{\hslash }^2}}-{\displaystyle \frac{m^2{\omega }^2}{{\hslash }^2}}{x}^2\right)y\left(x\right)=0.$$

(21)

Introducing the dimensionless coordinate

$$\xi =\sqrt{{\displaystyle \frac{m\omega }{\hslash }}}x,$$

(22)

Eq. (21) reduces to the Weber equation

$${\displaystyle \frac{d^{2}y}{d{\xi }^2}}+\left(\nu +{\displaystyle \frac{1}{2}}-{\displaystyle \frac{{\xi }^2}{4}}\right)y=0,\nu :={\displaystyle \frac{E}{\hslash \omega }}-{\displaystyle \frac{1}{2}}.$$

(23)

A natural independent basis is provided by parabolic cylinder functions,

$$y_1\left(\xi \right)=D_{\nu }\left(\xi \right),y_2\left(\xi \right)=D_{\nu }(-\xi ).$$

(24)

The general Bohmian amplitude solving the nonlinear EP equation (20) can then be written in Pinney form using the same basis:

$$R^2\left(\xi \right)=A{y}_1^2\left(\xi \right)+B{y}_2^2\left(\xi \right)+2D{y}_1\left(\xi \right)y_2\left(\xi \right),AB-D^2={\displaystyle \frac{k}{W^2}},$$

(25)

where $W=y_1{y}_2^{\prime }-y_1^{\prime }{y}_2$ is the (constant) Wronskian and $k=C^2/{\hslash }^2$ as in Eq. (12). The invariant defined in Eq. (13) is therefore constant along x for any choice of y solving the Weber equation and R given by (25).

• Why Weber Basis?

The construction of the Ermakov invariant requires a choice of a fundamental pair of independent solutions of the associated linear equation. Naturally, for the harmonic oscillator, this linear partner equation arises in the corresponding separable coordinate system. A convenient representation of its independent solutions is provided by the parabolic cylinder functions $D_{\nu }\left(\xi \right)$ and $D_{\nu }(-\xi )$. In contrast to the conventional treatment, where the harmonic oscillator is formulated directly in terms of Hermite functions after imposing normalizability, the Ermakov–Pinney framework retains access to the full space of independent solutions prior to the application of boundary conditions. Accordingly, the amplitude admits the invariant quadratic form (25).

Here the Wronskian W ensures linear independence, while the invariant guarantees consistency of the amplitude construction. Boundary and normalization conditions then select particular invariant combinations, rather than eliminating solutions at the outset. For discrete values $\nu =n$, the Weber functions reduce to Hermite functions, $D_n\left(\xi \right)\propto e^{-{\xi }^2/4}{H}_n(\xi /\sqrt{2})$, recovering the standard harmonic-oscillator spectrum.

3.3. Coulomb Potential

For the one-dimensional Coulomb potential $V\left(x\right)=-\alpha /x$ on the half-line, assume bound states $E<0$ and define the decay scale $\beta =\sqrt{-2mE/{\hslash }^2}.$ After the standard rescaling, the stationary Schrödinger equation reduces to Whittaker’s differential equation, which provides the natural representation of solutions for Sturm–Liouville problems with a Coulomb-type singularity [15]. We therefore adopt Whittaker functions as the fundamental solution basis. With the Whittaker variable $z=2\beta x$, a fundamental solution pair is given by $M_{\kappa ,1/2}\left(z\right)$ and $W_{\kappa ,1/2}\left(z\right)$, where the Whittaker index is $\kappa =m\alpha /\left({\hslash }^2\beta \right).$ This leads to the general Ermakov amplitude

$${\rho }^2\left(2\beta x\right)=A{M}_{\kappa ,1/2}^2\left(2\beta x\right)+B{W}_{\kappa ,1/2}^2\left(2\beta x\right)+2C{M}_{\kappa ,1/2}\left(2\beta x\right)W_{\kappa ,1/2}\left(2\beta x\right),$$

(26)

with $AB-C^2=k/{\mathcal{W}}^2$, where $\mathcal{W}$ is the Wronskian of the above solution pair. Imposing the bound-state condition $\kappa =n+1$ yields $E_n=-m{\alpha }^2/\left(2{\hslash }^2{(n+1)}^2\right)$ and reduces $M_{\kappa ,1/2}$ to the associated Laguerre form,

$$M_{n+1,1/2}\left(2\beta x\right)\propto 2\beta x{e}^{-\beta x}{L}_n^{\left(1\right)}\left(2\beta x\right).$$

Thus, the textbook Coulomb bound states are recovered as a constrained Ermakov limit.

4. Invariant Structure, Implications and Extensions

4.1. Invariant Structure and Scope

Beyond establishing the existence of an Ermakov–Pinney equation for separable stationary problems, it is instructive to examine the nature and scope of the associated invariant in the examples considered above. The Ermakov–Lewis invariant can be written in the form

$$I={\displaystyle \frac{1}{2}}\left[{(R{y}^{\prime }-R^{\prime }y)}^2+k{\left({\displaystyle \frac{y}{R}}\right)}^2\right],$$

(27)

where y is any solution of the associated linear equation and R solves the Ermakov–Pinney equation.

In the linear limit $k=0$, the invariant reduces to half the squared Wronskian of a fundamental solution pair of the linear equation. This corresponds to a spatial analogue of an adiabatic invariant, which is implicitly present in standard stationary quantum mechanics but is not usually identified as such. From this perspective, the Ermakov–Pinney formulation makes explicit a conserved structure that is otherwise absorbed into basis selection.

For nonzero k, the invariant retains a bounded value even when individual solutions of the linear equation exhibit divergent behavior. This occurs because the invariant depends on a Wronskian-like combination rather than on individual solutions themselves. Consequently, the Ermakov–Pinney framework naturally accommodates general linear combinations of independent solutions and is well suited to boundary-value formulations, without requiring the a priori exclusion of solutions based on asymptotic considerations.

Within this perspective, the Bohm–Madelung equations describe the underlying stationary dynamics, while the Ermakov–Lewis invariant supplies an additional integrability structure that constrains and organizes those dynamics. The invariant acts in amplitude space as a dual to Hamilton–Jacobi separability, providing bounds and consistency conditions without altering the physical content of the theory. Thus, while the Ermakov–Bohm framework reorganizes the description of amplitude dynamics, it preserves the probabilistic predictions of stationary quantum mechanics.

4.2. Implications, Limitations and Extensions

This explicit coordinate-level analysis confirms that the appearance of EP dynamics in stationary quantum mechanics is a structural consequence of separability and not an artifact of a particular coordinate choice. It also clarifies why such invariants fail to exist in systems with non-diagonal kinetic terms or gauge-induced momentum couplings, such as motion in magnetic fields.

The EP framework discussed above captures invariant amplitude structures associated with separable stationary problems. Additional structure may emerge from a variational treatment of the BM equations, resulting in closed relations that regularize the wave function and impose global constraints. These aspects are explored further in [16].

5. Conclusions

We have shown that Ermakov–Pinney equations and their associated invariants arise naturally in stationary quantum mechanics whenever the Hamiltonian is diagonal and separable, within the BM formulation. In this setting, the spatial coordinate plays the role of an evolution parameter, and the Ermakov invariant encodes a hidden geometric structure of the amplitude dynamics.

The framework unifies the free particle, harmonic oscillator, and Coulomb potential within a common invariant-based description and preserves the standard probabilistic interpretation of quantum mechanics. While additional structures emerge for non-separable systems and angular variables, these are beyond the scope of the present work and are addressed in an extended Hamilton–Jacobi formulation.

The results presented here emphasize that Ermakov invariants are not restricted to time-dependent problems but are intrinsic to the structure of stationary quantum mechanics under suitable separability conditions.

Appendix A.1. Relation to Hamilton–Jacobi Separability

In classical mechanics, Hamilton–Jacobi separability guarantees integrability of the phase function S. In Bohmian mechanics, this structure is mirrored:

• the phase $S_i$ satisfies Hamilton–Jacobi–type equations,
• the amplitude $R_i$ satisfies EP equations.

Appendix A.2. General Form of the Quantum Potential

The Bohm quantum potential is defined as

$$Q=-{\displaystyle \frac{{\hslash }^2}{2m}}{\displaystyle \frac{{\nabla }^{2}R}{R}},$$

(A1)

where R is the amplitude in the polar decomposition $\psi =R{e}^{iS/\hslash }$.

Assume separability of the amplitude,

$$R(q_1,\dots ,q_n)=\prod _{i=1}^n{R}_i\left(q_i\right),$$

(A2)

in an orthogonal coordinate system $\left\{q_i\right\}$ with diagonal metric coefficients $h_i\left(q_i\right)$.

Substituting the product form of R yields

$${\displaystyle \frac{{\nabla }^{2}R}{R}}={\displaystyle \frac{1}{{\prod }_{j=1}^n{h}_j}}\sum _{i=1}^n{\displaystyle \frac{\partial }{\partial q_i}}\left({\displaystyle \frac{{\prod }_{j=1}^n{h}_j}{h_i^2}}{\displaystyle \frac{1}{R}}{\displaystyle \frac{\partial R}{\partial q_i}}\right)={\displaystyle \frac{1}{{\prod }_{j=1}^n{h}_j}}\sum _{i=1}^n{\displaystyle \frac{\partial }{\partial q_i}}\left({\displaystyle \frac{{\prod }_{j=1}^n{h}_j}{h_i^2}}{\displaystyle \frac{R_i^{\prime }}{R_i}}\right).$$

(A3)

demonstrating that the quantum potential separates additively,

$$Q=\sum _{i=1}^n{Q}_i\left(q_i\right).$$

(A4)

Equivalently, expanding the derivative gives a form that isolates the usual second-derivative term and a “measure” contribution from the coordinate system:

$${\displaystyle \frac{{\nabla }^{2}R}{R}}=\sum _{i=1}^n{\displaystyle \frac{1}{h_i^2}}{\displaystyle \frac{R_i^{\prime \prime }}{R_i}}+\sum _{i=1}^n{\displaystyle \frac{1}{h_i^2}}\left[{\displaystyle \frac{\partial }{\partial q_i}}ln\left({\displaystyle \frac{{\prod }_{j=1}^n{h}_j}{h_i^2}}\right)\right]{\displaystyle \frac{R_i^{\prime }}{R_i}}.$$

(A5)

This additive structure is essential for the emergence of independent EP equations for each degree of freedom.

• Additive structure and emergence of EP equations.

For diagonal and separable Hamiltonians in orthogonal coordinates, the Bohm quantum potential admits an additive decomposition,

$$Q=-{\displaystyle \frac{{\hslash }^2}{2m}}{\displaystyle \frac{{\nabla }^{2}R}{R}}=\sum _i{Q}_i\left(q_i\right),$$

(A6)

provided the amplitude factorizes as $R={\prod }_i{R}_i\left(q_i\right)$ and the coordinate operators are in Sturm–Liouville form. This additive structure is essential, as it allows each coordinate contribution to enter the quantum Hamilton–Jacobi equation independently.

At the same time, separation of the stationary continuity equation implies that each Bohmian momentum component satisfies

$$p_i\left(q_i\right)={\displaystyle \frac{C_i}{R_i^2\left(q_i\right)}},$$

(A7)

where $C_i$ are separation constants subject only to a global constraint. Substitution into the quantum Hamilton–Jacobi equation then yields, for each degree of freedom, a nonlinear second-order equation of the EP type,

$$R_i^{″}\left(q_i\right)+{\Omega }_i^2\left(q_i\right)R_i\left(q_i\right)={\displaystyle \frac{1}{R_i^3\left(q_i\right)}}{\displaystyle \frac{C_i^2}{{\hslash }^2}}.$$

(A8)

Thus, the EP structure arises directly from the combination of additive separability and the Bohm continuity constraint.