The nonlinear Schrödinger equation in cylindrical geometries

The nonlinear Schrödinger equation was originally derived in nonlinear optics as a model for beam propagation, which naturally requires its application in cylindrical coordinates. However, that derivation was performed in Cartesian coordinates for linearly polarized fields with the Laplacian $\Delta_{\perp} = \partial_{x}^{2} + \partial_{y}^{2}$ transverse to the beam z-direction, and then, tacitly assuming covariance, extended to axisymmetric cylindrical setting. As we show, first with a simple example and next with a systematic derivation in cylindrical coordinates for axisymmetric and hence radially polarized fields, $\Delta_{\perp} = \partial_{r}^{2} + \frac{1}{r} \partial_{r}$ must be amended with a potential $V(r) = \frac{1}{r^{2}}$ , which leads to a Gross-Pitaevskii equation instead. Hence, results for beam dynamics and collapse must be revisited in this setting.