
The beam propagation method (BPM) is an approximation technique for simulating the propagation of light in slowly varying optical waveguides. It is essentially the same as the socalled parabolic equation (PE) method in underwater acoustics. Both BPM and the PE were first introduced in the 1970s. When a wave propagates along a waveguide for a large distance (larger compared with the wavelength), rigorous numerical simulation is difficult. The BPM relies on approximate differential equations which are also called the oneway models. These oneway models involve only a first order derivative in the variable z (for the waveguide axis) and they can be solved as "initial" value problem. The "initial" value problem does not involve time, rather it is for the spatial variable z. The original BPM and PE were derived from the slowly varying envelope approximation and they are the socalled paraxial oneway models. Since then, a number of improved oneway models are introduced. They come from a oneway model involving a square root operator. They are obtained by applying rational approximations to the square root operator. After a oneway model is obtained, one still has to solve it by discretizing the variable z. However, it is possible to merge the two steps (rational approximation to the square root operator and discretization of z) into one step. Namely, one can find rational approximations to the socalled oneway propagator (the exponential of the square root operator) directly. The rational approximations are not trivial. Standard diagonal Padé approximants have trouble with the socalled evanescent modes. These evanescent modes should decay rapidly in z, but the diagonal Padé approximants will incorrectly propagate them as propagating modes along the waveguide. Modified rational approximants that can suppress the evanescent modes are now available. The accuracy of the BPM can be further improved, if you use the energyconserving oneway model or the singlescatter oneway model. ==Principles== BPM is generally formulated as a solution to Helmholtz equation in a timeharmonic case, 〔Okamoto K. 2000 Fundamentals of Optical Waveguides (San Diego, CA: Academic)〕 〔EE290F: BPM course slides, Devang Parekh, University of Berkeley, CA〕 :$$ (\nabla^2 + k_0^2n^2)\psi = 0 with the field written as, :$E(x,y,z,t)=\backslash psi(x,y)\backslash exp(j\backslash omega\; t)$. Now the spatial dependence of this field is written according to any one TE or TM polarizations :$\backslash psi(x,y)\; =\; A(x,y)\backslash exp(+jk\_o\backslash nu\; y)$ , with the envelope :$A(x,y)$ following a slowly varying approximation, :$$ \frac = 0 Now the solution when replaced into the Helmholtz equation follows, :$$ \left(ウィキペディア（Wikipedia）』 ■ウィキペディアで「Beam propagation method」の詳細全文を読む スポンサード リンク
