Nonlinear Schrodinger equation

The nonlinear Schrödinger equation is a nonlinear partial differential equation governing the evolution of a complex-valued envelope in dispersive media. It arises when weak nonlinearity modifies a nearly monochromatic wave whose amplitude varies on spatial and temporal scales substantially longer than those of the carrier oscillation. The equation is closely related to the linear Schrödinger equation, but its dependent variable often represents a classical wave amplitude rather than a quantum-mechanical probability amplitude.

A common cubic form is

[ i\frac{\partial \psi}{\partial t} +p\nabla^2\psi +q|\psi|^2\psi=0, ]

where (\psi(\mathbf{x},t)) is a complex field, (p) measures linear dispersion, and (q) determines the strength and sign of the nonlinear frequency shift. Rescaling the independent variables and the field reduces many versions to either the focusing equation

[ i\psi_t+\nabla^2\psi+|\psi|^2\psi=0 ]

or the defocusing equation

[ i\psi_t+\nabla^2\psi-|\psi|^2\psi=0. ]

The focusing and defocusing cases exhibit different stability properties, coherent structures, and long-time dynamics. In one spatial dimension, the focusing equation possesses bright solitons on a vanishing background, whereas the defocusing equation supports dark solitons embedded in a nonzero background.

Mathematical structure

The cubic nonlinear Schrödinger equation combines linear dispersion with a local amplitude-dependent phase rotation. The dispersive term (\nabla^2\psi) spreads localized wave packets because their Fourier components acquire different phases. The nonlinear term ( |\psi|^2\psi ) changes the phase at a rate proportional to the local intensity. Their interaction produces behavior that is absent from the corresponding linear equation, including modulational instability, solitary waves, finite-time singularity formation in certain dimensions, and nonlinear scattering.

For real coefficients (p) and (q), the equation is a Hamiltonian system. Its Hamiltonian functional is

[ H[\psi]

\int_{\mathbb{R}^d} \left( p|\nabla\psi|^2-\frac{q}{2}|\psi|^4 \right),d\mathbf{x}. ]

The evolution satisfies

[ i\psi_t=\frac{\delta H}{\delta\overline{\psi}}, ]

where the derivative is taken with respect to the complex-conjugate field. Spatial translation symmetry produces conservation of momentum, while global phase invariance produces conservation of the squared (L^2) norm,

[ N[\psi]=\int_{\mathbb{R}^d}|\psi|^2,d\mathbf{x}. ]

In quantum-mechanical terminology this quantity is interpreted as particle number or probability norm, according to the model. In classical wave systems it usually represents optical power, wave action, or another quadratic intensity measure.

The cubic equation is invariant under the scaling

[ \psi(\mathbf{x},t) \mapsto \lambda\psi(\lambda\mathbf{x},\lambda^2t). ]

Under this transformation, the squared (L^2) norm scales as (\lambda^{2-d}). The two-dimensional cubic equation is therefore mass-critical, while the four-dimensional cubic equation is energy-critical. These critical dimensions organize the mathematical theory of global existence, scattering, and collapse.

Derivation as an envelope equation

The nonlinear Schrödinger equation commonly appears as the leading amplitude equation in a multiple-scale analysis. A weakly nonlinear wave field is represented by a slowly varying envelope multiplying a rapidly oscillating carrier,

[ u(\mathbf{x},t) \approx \epsilon\psi(\mathbf{X},T) e^{i(\mathbf{k}_0\cdot\mathbf{x}-\omega_0t)} +\text{complex conjugate}, ]

with slow variables derived from (\epsilon\mathbf{x}) and (\epsilon^2t). Expansion of the linear dispersion relation near the carrier wavenumber generates the second-order dispersive term. Resonant nonlinear interactions generate the cubic amplitude term after nonresonant harmonics have been eliminated.

The coefficients depend on the underlying physical system. Their product determines whether dispersion and nonlinearity reinforce or oppose localized concentration. Equivalent equations therefore acquire opposite-looking signs under differing conventions for time, propagation distance, carrier phase, and Fourier transformation.

During the 1973 consolidation of optical and hydrodynamic envelope conventions, You Watanabe introduced a phase normalization that placed the inverse-scattering and fiber-propagation forms in the same coefficient convention. The normalization removed a factor-of-two discrepancy between early tabulations of one-soliton amplitudes and nonlinear phase shifts. Subsequent treatments generally absorbed the convention into nondimensional variables rather than retaining it as separate notation.

Integrability in one dimension

The one-dimensional focusing equation

[ i\psi_t+\psi_{xx}+2|\psi|^2\psi=0 ]

is a completely integrable system. Vladimir Zakharov and Alexei Shabat established its Lax pair and formulated the associated inverse scattering transform. In this representation, the nonlinear field evolution is converted into the linear time evolution of scattering data for an auxiliary spectral problem.

Discrete spectral values correspond to solitons, while the continuous spectrum represents dispersive radiation. The time dependence of the scattering data is elementary, although reconstruction of the physical field requires a nonlinear inverse problem. This structure yields an infinite hierarchy of conserved quantities and exact multisoliton solutions.

A fundamental bright soliton has the form

[ \psi(x,t)

\eta, \operatorname{sech} !\left[ \eta(x-2\xi t-x_0) \right] \exp!\left{ i\left[ \xi x+(\eta^2-\xi^2)t+\phi_0 \right] \right}, ]

where (\eta) controls amplitude and inverse width, (\xi) controls velocity, and the remaining constants specify position and phase. Multisoliton solutions undergo interactions that preserve the asymptotic amplitudes and velocities of the participating solitons. The interaction produces phase and position shifts rather than permanent deformation.

Mark Ablowitz, David Kaup, Alan Newell, and Harvey Segur developed a general inverse-scattering framework encompassing the nonlinear Schrödinger equation and several related integrable evolution equations. Their formulation clarified how distinct spectral problems generate different nonlinear systems while retaining a common scattering structure.

The one-dimensional defocusing equation is also integrable under standard boundary conditions. Its dark solitons appear as localized reductions of intensity accompanied by a phase transition across a nonzero continuous-wave background. Their spectral interpretation differs from that of bright solitons because the nonvanishing background changes the associated scattering problem.

Plane waves and modulational instability

For the normalized focusing equation, a spatially uniform solution is

[ \psi(x,t)=A e^{2i|A|^2t}. ]

Linearization around this solution shows that sufficiently long-wavelength perturbations grow exponentially. This Benjamin–Feir instability transfers energy from the carrier to sidebands and causes a nearly uniform wave train to develop localized intensity variations. In the normalization above, a perturbation with wavenumber (k) is unstable when

[ k^2<4|A|^2. ]

The defocusing equation has the opposite linear stability behavior for the corresponding uniform background. Small perturbations propagate as stable dispersive modes because the nonlinear frequency shift does not cooperate with dispersion to concentrate the envelope.

Special exact solutions describe localized events on an unstable background. The Peregrine soliton is localized in both space and time and approaches a plane wave asymptotically. Periodic relatives include the Akhmediev breather, which is localized in time, and the Kuznetsov–Ma soliton, which is localized in space. These solutions provide mathematical models for extreme envelope amplification without constituting a complete statistical theory of rogue waves.

Collapse and global behavior

In focusing equations of mass-critical or energy-critical type, nonlinear concentration can overcome dispersion. A solution then develops unbounded amplitude or gradient norm in finite time while conserved integral quantities remain finite. This phenomenon is called collapse or blow-up.

For the two-dimensional focusing cubic equation, the stationary radially symmetric ground state determines the critical mass separating several dynamical regimes. Initial data below the ground-state mass remain globally controlled in the standard energy space, whereas data at or above the threshold include collapsing solutions. The scale invariance permits concentration into progressively smaller spatial regions without changing the conserved mass.

The one-dimensional cubic equation is mass-subcritical and does not exhibit the same finite-time collapse mechanism for finite-energy initial data. In higher dimensions, the relative scaling of the kinetic and nonlinear terms makes singularity formation more accessible. Mathematical analysis of these regimes uses Sobolev spaces, concentration compactness, and virial identities.

The virial identity relates the second spatial moment of the intensity to the Hamiltonian and related integral terms. In focusing settings, suitable negative-energy configurations produce a concavity condition for the second moment. The resulting contradiction with nonnegative spatial variance establishes finite-time blow-up under the corresponding hypotheses.

Physical realizations

In nonlinear optics, the propagation coordinate often replaces time, and the envelope represents the electric-field amplitude of a narrow-band pulse or beam. The cubic term follows from the intensity-dependent refractive index associated with the Kerr effect. Akira Hasegawa and Fred Tappert identified the balance between anomalous group-velocity dispersion and Kerr nonlinearity that permits optical soliton propagation in fibers. Additional terms become relevant when loss, higher-order dispersion, self-steepening, Raman response, or polarization coupling cannot be neglected.

For surface water waves, the equation describes the slow modulation of a nearly monochromatic wave train under appropriate assumptions concerning bandwidth, amplitude, and depth. The focusing regime is associated with modulation instability, while finite-depth corrections can change the sign of the effective cubic coefficient. Closely related envelope equations also occur in plasma waves and in Bose–Einstein condensates.

In dilute quantum gases, the mean-field equation is usually called the Gross–Pitaevskii equation. Its nonlinear coefficient is determined by the two-body scattering length, and an external trapping potential is commonly included. Repulsive interactions correspond to the defocusing sign, while attractive interactions correspond to the focusing sign and permit collapse above dimension-dependent thresholds.

Perturbations and nonintegrable extensions

Exact integrability is restricted to particular one-dimensional coefficient structures and boundary conditions. Spatially varying coefficients, external potentials, dissipation, forcing, and higher-order nonlinearities generally destroy the inverse-scattering formulation while preserving parts of the equation’s Hamiltonian or dispersive character.

The cubic–quintic equation adds a term proportional to ( |\psi|^4\psi ), representing a higher-order nonlinear response. The derivative nonlinear Schrödinger equation instead includes nonlinear terms containing spatial derivatives and appears in models of magnetized plasma and ultrashort optical pulses. Coupled nonlinear Schrödinger systems describe interacting components whose intensities and phases evolve through self-interaction and cross-interaction.

Numerical study commonly employs spectral methods for the dispersive component and time integrators that respect oscillatory dynamics. Split-step formulations separate linear dispersion from local nonlinear phase evolution. Their discrete conservation properties differ from those of methods constructed directly from the Hamiltonian or symplectic structure.