🌀 nonlinear sciences

Dimensional reduction for optical beams with thermal nonlocal nonlinearity

This paper extends dimensional reduction techniques to thermal nonlocal nonlinear media by coupling the optical field to temperature-induced refractive index changes and employing a Laguerre-Gauss variational ansatz to derive an efficient analytic model that captures mode-dependent nonlinear interactions.

Original authors: Francesco Lorenzi, Luca Salasnich

Published 2026-06-23

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Francesco Lorenzi, Luca Salasnich

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 you are shining a flashlight through a thick, slightly warm piece of glass. Usually, light just travels in a straight line or spreads out a little. But in certain materials, the light actually changes the glass itself as it passes through. The light gets absorbed, heats up the glass, and that heat changes how the glass bends light. This creates a feedback loop: the light shapes the glass, and the shaped glass guides the light.

This paper is about figuring out a simpler way to predict exactly how that light beam will behave, without needing a supercomputer to solve millions of tiny equations at once.

Here is the breakdown of their work using everyday analogies:

The Problem: Too Much Math for a Simple Beam

When light travels through these "thermal" materials, it doesn't just interact with the spot it's currently hitting. Because heat spreads out (diffuses) like steam rising from a cup of coffee, the light at one point is affected by the heat generated by the light everywhere else in the beam.

Scientists call this "nonlocal nonlinearity." It's like if you pushed a heavy mattress; the whole mattress moves, not just the spot under your hand. To simulate this accurately, you usually have to run a massive, 3D computer simulation that calculates the heat and light interaction at every single point in space. This is slow and computationally expensive.

The Solution: A "Shadow" Model

The authors wanted to create a "dimensional reduction." Think of this like taking a complex, 3D sculpture and casting a shadow of it on a wall. The shadow is much simpler to draw and analyze, but if you choose the right angle, it still tells you everything important about the sculpture's shape and movement.

They took the complex 3D equations governing the light and heat and "squashed" them down into a simpler 1D model. Instead of tracking the light's shape across a wide area, they track just a few key numbers:

How bright the beam is (Amplitude).
How wide the beam is (Width).
How the beam is curving (Phase/Chirp).

The Secret Ingredient: The "Ring" Shape

To make this shadow accurate, they had to guess what the beam looks like. They used a specific mathematical shape called a "Laguerre-Gauss" mode.

The Analogy: Imagine a donut or a ring of light rather than a solid circle. These rings can have different sizes and "twists" (like a corkscrew).
The Twist: They realized that just guessing the size of the ring wasn't enough. They also had to include a "radial phase" term. Think of this as the beam having a slight "twist" or "slope" in its wave pattern as it moves outward. Without this twist, their simplified model would fail to predict how the beam spreads or focuses correctly.

What They Found

By using this "ring" guess and the "twist" factor, they derived a new set of equations. These equations act like a rulebook for the beam's behavior:

The Balance: The beam is constantly fighting a tug-of-war. On one side, the light wants to spread out (diffraction). On the other, the heat it creates pulls it together or pushes it apart (thermal nonlinearity).
The Result: Their simplified model successfully predicts how the beam's width changes as it travels. It captures how the beam reacts to different "twists" (vortex orders) and how "far-reaching" the heat effect is.

Did It Work?

They tested their new "shadow" model against the heavy-duty 3D computer simulations (the "real" physics).

The Good News: For most situations, especially when the beam is moving, their simple model matched the complex simulation almost perfectly. It correctly predicted how the beam would breathe (expand and contract) and how it would lose energy.
The Limitation: When the heat effect is extremely long-range (like a very large, slow-moving wave), the actual beam shape gets a bit weird and distorted. In those specific cases, the simple "ring" guess isn't a perfect match for the messy reality, but it still gets the general behavior right.

The Takeaway

This paper doesn't invent a new laser or a new material. Instead, it provides a shortcut. It gives scientists a much faster, easier-to-use mathematical tool to understand how light beams behave in heat-sensitive materials. It's like having a reliable map instead of having to walk every single step of the journey to know where you are going. This is particularly useful for designing optical fibers (the cables that carry internet data) where managing these light beams is crucial.

Technical Summary: Dimensional Reduction for Optical Beams with Thermal Nonlocal Nonlinearity

Problem Statement
Optical nonlinearities arising from the thermorefractive effect provide a long-range material response governed by heat diffusion and absorption. In graded-index (GRIN) media, this nonlocality fundamentally alters modal interactions, such as stabilizing vortex beams and ring solitons that would otherwise undergo azimuthal instability in local cubic media. However, accurately modeling these dynamics starting from the full spatial nonlinear Schrödinger equation (NLSE) coupled to a heat diffusion equation is computationally demanding. While reduced-dimensional models exist for local Kerr nonlinearities (e.g., the Nonpolynomial Schrödinger Equation, or NPSE), they fail to account for the nonlocal nature of thermal responses. The authors address the need for a tractable, reduced description of beam dynamics in thermal nonlocal media that captures the dependence of nonlinear interactions on mode order and the degree of nonlocality.

Methodology
The authors extend the NPSE framework to incorporate thermally mediated nonlocal nonlinearities. The approach involves the following steps:

Governing Equations: The system is modeled by coupling the paraxial spatial NLSE for the optical field envelope $\Phi$ with a steady-state heat equation for the temperature-induced refractive index change. Under the assumption of weak absorption and specific boundary conditions, the heat equation is reduced to a 2D screened Poisson equation. The resulting optical propagation equation includes a nonlocal operator defined by a convolution with a modified Bessel function kernel ( $K_0$ ).
Variational Ansatz: To perform dimensional reduction, the authors employ a variational approach based on a Lagrangian formulation of the lossless system. They assume the transverse field profile is a single Laguerre-Gauss (LG) mode with arbitrary azimuthal order $m$ (where $S=|m|$ ) and zero radial index ( $p=0$ ), though the formalism allows for general $p$ . Crucially, the ansatz includes a radial phase chirp term, $b(z)r^2/2$ , which the authors identify as necessary to properly describe diffraction in nonlocal media.
Derivation of Effective Model: By integrating the Lagrangian density over the transverse plane, the authors derive an effective one-dimensional action. This process yields explicit analytic expressions for the nonlocal contribution, which depends on the mode indices and the variational transverse width $\sigma(z)$ . The nonlocal term is evaluated using the Hankel transform of the intensity profile, leading to a closed-form expression for the fundamental Gaussian mode ( $p=0, m=0$ ) and a general integral form for higher orders.
Equations of Motion: The Euler-Lagrange equations derived from this effective action result in a coupled system:
- An evolution equation for the axial amplitude $A(z)$ , which includes a $\sigma$ -dependent nonlinear term.
- An algebraic equation determining the transverse width $\sigma(z)$ based on the balance between local diffraction, the GRIN potential, and the nonlocal thermal pressure.
- An equation for the radial phase chirp $b(z)$ .

Key Results

Analytic Expressions: The paper derives explicit expressions for the nonlocal integral $I_{pS}(\sigma)$ , which quantifies the thermal interaction strength. For the fundamental Gaussian mode, this is expressed in terms of the exponential integral function $Ei$.
Role of Radial Chirp: The study demonstrates that the standard NPSE ansatz (without a radial phase term) is insufficient for nonlocal media. The inclusion of the chirp parameter $b(z)$ is essential to capture the correct beam dynamics, particularly the breathing oscillations of the beam width.
Stationary States: The authors analyze the effective potential $V_{eff}(\sigma)$ governing the beam width. They show that in the defocusing regime, increasing the nonlocality length $L_{eff}$ shifts the potential minimum toward larger widths, resulting in broader stationary beams. In the focusing regime, nonlocality creates a deeper potential well at small widths.
Numerical Validation: The reduced model is validated against full 2D split-step Fourier simulations of the paraxial diffraction equation.
- In regimes of short-range nonlocality ( $L_{eff} = 0.5$ ), the variational ansatz shows excellent agreement with exact stationary radial intensity profiles for various vortex orders.
- In regimes of strong nonlocality ( $L_{eff} = 5.0$ ), the exact beam profiles deviate from the assumed LG functional form. Consequently, the variational model shows discrepancies in the detailed shape of the stationary intensity profile.
- Despite these static profile discrepancies, the reduced model accurately captures the dynamical evolution of the beam, including the decay of power due to absorption and the damped breathing oscillations of the transverse width, even for large nonlocality lengths ( $L_{eff} = 40$ ).

Significance and Claims
The paper claims to provide a generalized framework for thermally mediated interactions in GRIN fibers. Its primary contribution is the derivation of a reduced 1D model that explicitly incorporates the dependence of nonlinear interactions on mode order and transverse width. The authors emphasize that while the model cannot perfectly reproduce the exact transverse profile in highly nonlocal regimes (where the beam shape distorts away from the LG ansatz), it remains a powerful tool for understanding the overall beam dynamics, such as stability and propagation evolution. The work extends the applicability of variational dimensional reduction techniques to systems with screened-Poisson type nonlocalities, which are relevant not only to thermal media but also to liquid crystals, electrostrictive effects, and plasma interactions. The authors note that the current model is restricted to continuous waves and conservative systems (with absorption added phenomenologically), and future generalizations could address multimode fields, saturable responses, and time-dependent thermal relaxation for pulsed systems.