The differential equation d x d t = x − x 2 {\displaystyle {\frac {dx}{dt}}=x-x^{2}} has two stationary (time-independent) solutions: x = 0 and x = 1. The linearization at x = 0 has the form d x d t = x {\displaystyle {\frac {dx}{dt}}=x} . The linearized operator is A0 = 1. The only eigenvalue is λ = 1 {\displaystyle \lambda =1} . The solutions to this equation grow exponentially; the stationary point x = 0 is linearly unstable.
To derive the linearization at x = 1, one writes d r d t = ( 1 + r ) − ( 1 + r ) 2 = − r − r 2 {\displaystyle {\frac {dr}{dt}}=(1+r)-(1+r)^{2}=-r-r^{2}} , where r = x − 1. The linearized equation is then d r d t = − r {\displaystyle {\frac {dr}{dt}}=-r} ; the linearized operator is A1 = −1, the only eigenvalue is λ = − 1 {\displaystyle \lambda =-1} , hence this stationary point is linearly stable.
The nonlinear Schrödinger equation i ∂ u ∂ t = − ∂ 2 u ∂ x 2 − | u | 2 k u , {\displaystyle i{\frac {\partial u}{\partial t}}=-{\frac {\partial ^{2}u}{\partial x^{2}}}-|u|^{2k}u,} where u(x,t) ∈ C and k > 0, has solitary wave solutions of the form ϕ ( x ) e − i ω t {\displaystyle \phi (x)e^{-i\omega t}} .4 To derive the linearization at a solitary wave, one considers the solution in the form u ( x , t ) = ( ϕ ( x ) + r ( x , t ) ) e − i ω t {\displaystyle u(x,t)=(\phi (x)+r(x,t))e^{-i\omega t}} . The linearized equation on r ( x , t ) {\displaystyle r(x,t)} is given by ∂ ∂ t [ Re r Im r ] = A [ Re r Im r ] , {\displaystyle {\frac {\partial }{\partial t}}{\begin{bmatrix}{\text{Re}}\,r\\{\text{Im}}\,r\end{bmatrix}}=A{\begin{bmatrix}{\text{Re}}\,r\\{\text{Im}}\,r\end{bmatrix}},} where A = [ 0 L 0 − L 1 0 ] , {\displaystyle A={\begin{bmatrix}0&L_{0}\\-L_{1}&0\end{bmatrix}},} with L 0 = − ∂ ∂ x 2 − k ϕ 2 − ω {\displaystyle L_{0}=-{\frac {\partial }{\partial x^{2}}}-k\phi ^{2}-\omega } and L 1 = − ∂ ∂ x 2 − ( 2 k + 1 ) ϕ 2 − ω {\displaystyle L_{1}=-{\frac {\partial }{\partial x^{2}}}-(2k+1)\phi ^{2}-\omega } the differential operators. According to Vakhitov–Kolokolov stability criterion,5 when k > 2, the spectrum of A has positive point eigenvalues, so that the linearized equation is linearly (exponentially) unstable; for 0 < k ≤ 2, the spectrum of A is purely imaginary, so that the corresponding solitary waves are linearly stable.
It should be mentioned that linear stability does not automatically imply stability; in particular, when k = 2, the solitary waves are unstable. On the other hand, for 0 < k < 2, the solitary waves are not only linearly stable but also orbitally stable.6
V.I. Arnold, Ordinary Differential Equations. MIT Press, Cambridge, MA (1973) ↩
P. Glendinning, Stability, instability and chaos: an introduction to the theory of nonlinear differential equations. Cambridge university press, 1994. ↩
V.V. Nemytskii, V.V. Stepanov, "Qualitative theory of differential equations", Princeton Univ. Press (1960) ↩
H. Berestycki and P.-L. Lions (1983). "Nonlinear scalar field equations. I. Existence of a ground state". Arch. Rational Mech. Anal. 82 (4): 313–345. Bibcode:1983ArRMA..82..313B. doi:10.1007/BF00250555. S2CID 123081616. /wiki/Bibcode_(identifier) ↩
N.G. Vakhitov and A.A. Kolokolov (1973). "Stationary solutions of the wave equation in the medium with nonlinearity saturation". Radiophys. Quantum Electron. 16 (7): 783–789. Bibcode:1973R&QE...16..783V. doi:10.1007/BF01031343. S2CID 123386885. /wiki/Bibcode_(identifier) ↩
Manoussos Grillakis, Jalal Shatah, and Walter Strauss (1987). "Stability theory of solitary waves in the presence of symmetry. I". J. Funct. Anal. 74: 160–197. doi:10.1016/0022-1236(87)90044-9.{{cite journal}}: CS1 maint: multiple names: authors list (link) https://doi.org/10.1016%2F0022-1236%2887%2990044-9 ↩