Polarization-Induced Beam Bending: Mathematical Model, Discretization, and Algorithm
This work provides a numerical tool for studying polarization-induced beam bending over long propagation distances, which is relevant for applications in nonlinear optics and beam manipulation.
The authors developed a conservative numerical scheme for a reduced hydrodynamic model of paraxial vector beam propagation, enabling long-distance simulations over tens of meters with millimeter resolution. The method reproduces the analytically predicted quadratic beam bending at short distances and reveals systematic deviations at longer ranges due to nonlinear phase accumulation and dispersive effects.
We study a reduced hydrodynamic formulation of paraxial vector beam propagation in which the beam intensity, optical phase, and spatially-dependent polarization are coupled through a nonlinear dispersive system. While prior analytical work derived a solution for the beam path valid for short propagation distances, a fully resolved numerical treatment of the model over long ranges has not previously been available. Here we present a conservative numerical scheme for the coupled system, combining a finite-volume discretization of the intensity equation with monotone Hamilton--Jacobi (H-J) solvers for the phase dynamics and upwind transport of polarization. The method preserves the nonnegativity of the intensity and remains stable under long-distance propagation. We perform large-scale simulations over propagation distances of tens of meters, while resolving millimeter-scale transverse structure. The numerical results reproduce the analytically predicted and experimentally observed quadratic beam bending at short distances and reveal systematic deviations beyond the asymptotic regime. These deviations arise from nonlinear phase accumulation and dispersive effects captured by the full model but are neglected in the short-distance approximation.