Madelung–Bohm trajectories in the paraxial and Schrödinger wave equations

In this work, we show that the Madelung–Bohm trajectories determined by any solution, Ψ( ξ , ζ ), of the paraxial wave equation in an arbitrary two-dimensional optical medium, where ζ is the propagation coordinate and ξ the transversal one, are given by the contours of the function G ( ξ , ζ )≡∫|Ψ( ξ , ζ )| 2 d ξ . Furthermore, we remark that an analogous result follows directly for the one-dimensional Schrödinger equation for a particle evolving in an arbitrary classical potential, where ξ gives the position of the particle at time ζ . We apply this result to several examples in free and quadratic optical media. They include the Airy beam, Hermite–Gaussian beams, and some of their superpositions. Additionally, this work provides an alternative derivation of Hermite–Gaussian (HG) beams in free space by employing the quantum potential equations and the framework of quantum trajectories. By analyzing known trajectories, we demonstrate that the fundamental parameters characterizing these beams—the Gouy phase, beam radius, and wavefront curvature—emerge naturally within the quantum potential formalism.