The numerical models of lasers and the most of nonlinear optical systems stem from Maxwell–Bloch equations (MBE). This full set of Partial Differential Equations includes Maxwell equations for electromagnetic field and semiclassical equations of the two-level (or multilevel) atoms. For this reason the simplified theoretical approaches were developed for numerical simulation of laser beams formation and their propagation since the early years of laser era. The Slowly varying envelope approximation of MBE follows from the standard nonlinear wave equation with nonlinear polarization P NL {\displaystyle \mathbf {P} ^{\text{NL}}} as a source:
∇ 2 E ( r → , t ) − n 2 c 2 ∂ 2 ∂ t 2 E ( r → , t ) = 1 ε 0 c 2 ∂ 2 ∂ t 2 P NL , {\displaystyle \nabla ^{2}{\cal {E({\vec {r}},t)}}-{\frac {n^{2}}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}{\cal {E({\vec {r}},t)}}={\frac {1}{\varepsilon _{0}c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}\mathbf {P} ^{\text{NL}},} where : E ( r → , t ) ∝ E ( r → ⊥ , z , t ) e i ( n k 0 ) ( z − c t ) + c.c. {\displaystyle {\cal {E}}({\vec {r}},t)\propto E({\vec {r}}_{\perp },z,t)e^{i(nk_{0})(z-ct)}+{\text{c.c.}}} resulting in the standard "parabolic" wave equation:
∂ E ∂ z + ω 0 k 0 c 2 ∂ E ∂ t − 1 2 k 0 i ∇ ⊥ 2 E = 1 ε 0 k 0 c 2 ∂ 2 ∂ t 2 P NL {\displaystyle {\frac {\partial E}{\partial z}}+{\frac {\;\omega _{0}\ }{k_{0}c^{2}}}{\frac {\partial E}{\partial t}}-{\tfrac {1}{2k_{0}}}\ i\ \nabla _{\perp }^{2}E={\frac {1}{\varepsilon _{0}k_{0}c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}\mathbf {P} ^{\text{NL}}~} , under conditions :
… excerpt ends here. Continue reading the full article.






