🌀 nonlinear sciences

Using the Ehrenfest theorem for determining the self-focusing and self-trapping of nonlinear beams

This paper generalizes the Ehrenfest theorem by modeling the nonlinear potential as a power- and width-dependent quantum harmonic oscillator to provide a simple, intuitive framework for analyzing the stability, breathing, and self-focusing properties of nonlinear waves described by the nonlinear Schrödinger equation, with results validated against numerical simulations.

Original authors: Chandroth P. Jisha, Stefan Nolte, Alessandro Alberucci

Published 2026-07-23

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Chandroth P. Jisha, Stefan Nolte, Alessandro Alberucci

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine light not just as a beam that travels in a straight line, but as a living, breathing thing that can change its own shape. In the world of physics, this happens when light gets really intense. Normally, if you shine a flashlight through a lens, the beam spreads out like butter on toast; this is called diffraction. But in certain materials, intense light can do the opposite: it can bend the material so that the material acts like a lens, pulling the light back together. This is called self-focusing. If this pull is just right, the light can get trapped in a tight, stable bundle that travels forever without spreading out or collapsing. Scientists call these bundles "solitons." Understanding how to predict exactly when light will spread, when it will collapse, and when it will get trapped is a huge deal for making better lasers, faster internet, and even for understanding how particles behave in the quantum world.

For decades, physicists have used complex math to predict these behaviors, often treating the light beam like a cloud of particles or using heavy calculus that requires supercomputers. But a new paper by researchers Chandroth P. Jisha, Stefan Nolte, and Alessandro Alberucci offers a clever shortcut. They decided to treat the width of a light beam like a tiny particle bouncing inside a special kind of box. By borrowing a famous rule from quantum mechanics called the "Ehrenfest theorem"—which basically says that the average motion of a quantum particle follows the same rules as a classical ball rolling on a hill—they created a simple, intuitive model. Instead of solving a massive, scary equation for every single point in the beam, they showed that you can describe the entire beam's behavior by watching just one number: its width.

Here is the magic trick they used: they imagined the light beam is trapped inside a "quantum harmonic oscillator." Think of this as a springy box that gets stiffer or looser depending on how much power the light has and how wide the beam is. If the light is too strong, the spring pulls so hard the beam collapses into a tiny dot (a disaster for lasers). If the light is too weak, the spring is too loose, and the beam spreads out. But if the power is just right, the beam finds a sweet spot where it bounces back and forth in a stable rhythm, or stays perfectly still. The authors proved that this simple "spring" model works for many different types of materials, including those where the light interacts with its neighbors (nonlocal materials) and those with complex, multi-layered effects.

The paper's main finding is that this "spring" model can predict the exact behavior of light beams with surprising accuracy, matching up perfectly with complex computer simulations. They showed that for simple materials (like standard glass), there is a specific "critical power" where the beam will inevitably collapse. However, in more complex materials (like those with a mix of focusing and defocusing effects), the beam can find a stable home, but only if the power is within a specific "Goldilocks zone"—not too low, not too high. The researchers also discovered that the starting conditions of the beam (how wide it is when it enters the material) don't just change the starting point; they actually change the shape of the "spring" itself. This means that by tweaking how you launch the beam, you can fundamentally alter how it behaves inside the material.

The authors didn't just guess this; they tested it rigorously. They compared their simple spring equations against heavy-duty computer simulations (called Beam Propagation Method) for various scenarios, including pure focusing, mixtures of focusing and defocusing, and materials where the effect spreads out over a distance. In every case, their simple model matched the complex simulations almost perfectly. They even found that for certain mixtures of materials, there is a maximum size a stable light beam can be, no matter how much power you pump into it.

So, what does this mean for the real world? It means scientists now have a much simpler, more intuitive way to design optical systems. Instead of needing a supercomputer to figure out if a laser beam will survive a trip through a special crystal, they can use this "spring" math to get a very good answer quickly. This could help in designing better lasers for cutting materials, improving fiber-optic communications, and even in understanding how light behaves in exotic states of matter. The paper suggests that this approach is versatile enough to handle even more complex situations in the future, like ultra-short pulses of light or materials that lose energy, opening the door to new ways of controlling the most powerful beams of light we can create.

Technical Summary: Ehrenfest Theorem for Nonlinear Beam Dynamics

Problem Statement
The propagation of nonlinear optical waves, particularly those governed by the Nonlinear Schrödinger Equation (NLSE), presents significant challenges in modeling self-focusing and self-trapping phenomena. While established methods such as the method of moments, variational approaches, and ray-optics models exist, they often require complex coupled systems of equations or specific ansatzes (e.g., Gaussian profiles) that may not fully capture the dynamics of beam width evolution, especially in regimes involving catastrophic collapse or complex nonlinear responses (local, nonlocal, higher-order). There is a need for a unified, intuitive framework that can describe the evolution of the beam width $w(z)$ across diverse nonlinear media while maintaining analytical tractability.

Methodology
The authors propose a generalized application of the Ehrenfest theorem from quantum mechanics to optical beam propagation. The core methodology involves:

Quantum Harmonic Oscillator (QHO) Analogy: The nonlinear potential $V_{NL}$ is modeled as a parabolic potential (a QHO) where the oscillator strength $a$ depends on the beam power $P$ and the beam width $w$ .
Derivation of the Width Equation: By applying the Ehrenfest theorem to the second-order moments of position ( $\hat{r}_T^2$ ) and momentum ( $\hat{p}^2$ ), the authors derive a single fourth-order ordinary differential equation (ODE) for the squared beam width $w^2$ .
Generalized Ansatz: The strength of the QHO is defined via the ansatz $a(z, P) = \gamma P^N / w^{2M}$ , where $N$ and $M$ are integers determined by the specific nonlinearity type (e.g., Kerr, cubic-quintic, nonlocal).
Reduction to Second-Order ODE: Through integration and the application of appropriate boundary conditions derived from the input beam parameters (waist $w_0$ , curvature radius $R_0$ , and power $P$ ), the fourth-order equation is reduced to a second-order ODE (Eq. 17). This equation governs the convexity of $w^2$ and is analogous to the equation of motion for a particle in an effective potential $U(w^2)$ .
Boundary Conditions: The model rigorously incorporates launch conditions, including the interplay between linear diffraction and the nonlinear phase shift induced at the interface, to determine integration constants.

Key Contributions

Unified Framework: The paper establishes a generalized formalism that reduces the complex dynamics of nonlinear wave propagation to a 1D mechanical problem (particle in a potential), applicable to local, nonlocal, and higher-order nonlinearities.
Analytical Solutions for Kerr Media: For pure Kerr media, the model yields closed-form solutions for the nonlinear beam waist and focal position, accurately predicting the critical power $P_{cr}$ and the onset of catastrophic collapse. It clarifies the distinction between the critical power for collapse and a secondary threshold related to the concavity of the beam width curve.
Cubic-Quintic Media Analysis: The approach successfully models the interplay between self-focusing (Kerr) and self-defocusing (higher-order) nonlinearities. It derives the existence conditions for stable self-trapped waves (solitons) and breathing solitons, identifying specific power intervals ( $P_{inf} < P < P_{sup}$ ) where stable confinement occurs.
Nonlocal Media Modeling: The method is extended to nonlocal media (e.g., nematic liquid crystals) by approximating the nonlocal response as a series expansion within the parabolic potential framework, capturing the transition from highly nonlocal to local limits.

Results

Kerr Nonlinearity: The theoretical predictions for Gaussian beams in Kerr media show excellent agreement with numerical simulations (Beam Propagation Method). The model correctly predicts that catastrophic collapse occurs at $P_{cr}$ , regardless of focusing conditions, and accurately describes the evolution of the beam waist $w_{0NL}$ and the position of collapse $z_c$ .
Cubic-Quintic Nonlinearity: The model predicts the existence of shape-preserving solitons and breathing oscillations. It identifies that stable self-trapping requires a balance where the defocusing term counteracts the Kerr effect, leading to a local minimum in the effective potential. Numerical comparisons confirm that the theoretical soliton width and existence intervals are accurate, with deviations less than 10% in flat regions, though the model predicts a finite existence interval while numerical simulations find solutions for higher powers (often with flat-top profiles).
Breathing Dynamics: The effective potential $U(w^2)$ explains the anharmonic breathing of self-trapped waves. The amplitude and period of breathing are shown to depend on the input power and the initial beam waist relative to the potential minimum.
Nonlocal Effects: The model demonstrates that nonlocality modifies the effective potential, stabilizing solitons and altering the power-width relationship. In the highly nonlocal limit, the results converge to previously established logarithmic potential forms.

Significance and Claims
The paper claims that this Ehrenfest-based approach offers a relatively simple and intuitive manner to explain the main properties of nonlinear confined waves, such as stability and breathing, without the computational overhead of full numerical simulations or the complexity of coupled variational equations.

Key claims regarding the significance of the work include:

Versatility: The method is applicable to a wide range of nonlinear responses (local, nonlocal, higher-order) within a single mathematical framework.
Physical Insight: By mapping the problem to a particle in a potential, the approach provides clear physical intuition regarding the balance between diffraction and self-phase modulation.
Predictive Power: The model provides closed-form expressions for critical parameters (critical power, soliton width, breathing characteristics) that agree well with numerical benchmarks.
Future Directions: The authors suggest that while the current work focuses on lossless continuous waves, the framework can be generalized to include nonlinear losses, ultrashort pulses (incorporating temporal dispersion and plasma effects), and non-cylindrically symmetric beams.

The paper does not propose new experimental setups but positions the theoretical model as a tool for interpreting existing phenomena and potentially guiding the measurement of nonlinear coefficients (e.g., the Kerr coefficient) via beam dynamics analysis prior to collapse.

Technical Summary: Ehrenfest Theorem for Nonlinear Beam Dynamics

More like this