Beam propagation method: Difference between revisions

Revision as of 15:19, 26 June 2012

Latest revision as of 00:02, 12 September 2023

(31 intermediate revisions by 25 users not shown)

Line 1:

The '''beam propagation method''' ('''BPM''') is an approximation technique for simulating the propagation of [[light]] in [[slowly varying opticalenvelope approximation|slowly varying]] [[optical waveguide]]s. It is essentially the same as the so-called [[Parabolicparabolic Equationequation]] (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 one-way models. These one-way 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.<ref>{{citation|title=Integrated Photonics|author=Clifford R. Pollock, Michal. Lipson|year= 2003|publisher=Springer |url=httphttps://books.google.com/books?id=DNJEoypcI6oC&pg=RA1-PA210&lpg=RA1-PA210&dq=%22Beam+propagation+method%22&pg=RA1-PA210|isbn=978-1-4020-7635-0}}</ref>

'''Beam Propagation Method''' (BPM) refers to a computational technique in [[Electromagnetics]], used

to solve the [[Helmholtz equation]] under conditions of a [[time-harmonic wave]]. BPM works under [[slowly varying envelope approximation]], for linear and nonlinear equations.

The original BPM and PE were derived from the [[slowly varying envelope approximation]] and they are the so-called paraxial one-way models. Since then, a number of improved one-way models are introduced. They come from a one-way model involving a square root operator. They are obtained by applying rational approximations to the square root operator. After a one-way 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 so-called one-way propagator (the exponential of the square root operator) directly. The rational approximations are not trivial. Standard diagonal Padé approximants have trouble with the so-called 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 energy-conserving one-way model or the single-scatter one-way model.

The '''beam propagation method''' (BPM) is an approximation technique for simulating the propagation of [[light]] in slowly varying optical [[waveguide]]s. It is essentially the same as the so-called [[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 one-way models. These one-way 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.<ref>{{citation|title=Integrated Photonics|author=Clifford R. Pollock, Michal. Lipson|year= 2003|publisher=Springer |url=http://books.google.com/?id=DNJEoypcI6oC&pg=RA1-PA210&lpg=RA1-PA210&dq=%22Beam+propagation+method%22|isbn=978-1-4020-7635-0}}</ref>

The original BPM and PE were derived from the slowly varying envelope approximation and they are the so-called paraxial one-way models. Since then, a number of improved one-way models are introduced. They come from a one-way model involving a square root operator. They are obtained by applying rational approximations to the square root operator. After a one-way 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 so-called one-way propagator (the exponential of the square root operator) directly. The rational approximations are not trivial. Standard diagonal Padé approximants have trouble with the so-called 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 energy-conserving one-way model or the single-scatter one-way model.

==Principles==

Line 14 ⟶ 11:

</math>

with the field written as,

:<math>E(x,y,z,t)=\psi(x,y,z)\exp(-j\omega t)</math>.<!-- Where is the variable z in the right hand side of equation? It should be in the complex amplitude: Phi(x,y,z). -->

Now the spatial dependence of this field is written according to any one [[Transverse mode|TE or TM]] polarizations

:<math>\psi(x,y) = A(x,y)\exp(+jk_o\nu y)

</math>,<!-- Confusion between a transverse direction called y and the propagation direction called y too. The propagation direction should be called z and replace y in the exponential term. -->

</math>,

with the envelope

:<math>A(x,y)

Line 24 ⟶ 22:

:<math>

\frac{\partial^2( A(x,y) )}{\partial y^2} = 0

</math><!-- The amplitude should vary slowly along the propagation, hence the derivative should be done with respect to z. -->

</math>

Now the solution when replaced into the Helmholtz equation follows,

:<math>

\left[\frac{\partial^2 }{\partial x^2} + k_0^2(n^2 - \nu^2) \right]A(x,y) = \pm 2 jk_0 \nu \frac{\partial A_k(x,zy)}{\partial zy}

</math><!-- Same problem. The left operator should contain a double derivative with respect to y and the derivative in the right hand side of the equation should be realized with respect to z. -->

</math>

With the aim to calculate the field at all points of space for all times, we only need to compute the function

<math>A(x,y)</math> for all space, and then we are able to reconstruct <math>\psi(x,y,z)</math>. Since the solution

is for the time-harmonic Helmholtz equation, we only need to calculate it over one time period. We can

visualize the fields along the propagation direction, or the cross section waveguide modes.

==Numerical Methodsmethods==

There are bothBoth ''spatial ___domain'' methods, and ''frequency (spectral) ___domain'' methods are available for the numerical solution of the discretized master equation. Upon discretization into a grid, (using various [[central difference|centralized difference]], [[Crank–Nicolson method]], FFT-BPM etc.) and field values rearranged in a causal fashion, the field evolution is computed through iteratioiteration, along the propagation direction. The spatial ___domain method computes the field at the next step (in the propagation direction) by solving a linear equation, whereas the spectral ___domain methods use the powerful forward/inverse [[Fast Fourier transform|DFT]] algorithms. Spectral ___domain methods have the advantage of stability even in the presence of nonlinearity (from refractive index or medium properties), while spatial ___domain methods can possibly become numerically unstable.

___domain method computes the field at the next step (in the propagation direction) by solving a linear equation, whereas the spectral ___domain

methods use the powerful forward/inverse [[Fast Fourier Transform|DFT]] algorithms. Spectral ___domain methods have the advantage of stability

even in the presence of nonlinearity (from refractive index or medium properties), while spatial ___domain methods can possibly become numerically unstable.

==Applications==

BPM is a quick and easy method of solving for fields in integrated optical devices. It is typically used only in solving for intensity and modes within shaped (bent, tapered, terminated) waveguide structures, as opposed to scattering problems. These structures typically consist of [[isotropic]] optical materials, but the BPM has also been extended to be applicable to simulate the propagation of light in general [[anisotropic]] materials such as [[liquid crystals]]. This allows one to [httphttps://wwwdoi.elisorg/10.ugent1364/OE.be/ELISgroups/lcd/research/bpm17.php010895 analyze] e.g. the polarization rotation of light in anisotropic materials, the tunability of a directional coupler based on liquid crystals or the light diffraction in LCD pixels.

BPM is a quick and easy method of solving for fields in integrated optical devices. It is typically

used only in solving for intensity and modes within shaped (bent, tapered, terminated) waveguide

structures, as opposed to scattering problems. These structures typically consist of [[isotropic]] optical materials, but the BPM has also been extended to be applicable to simulate the propagation of light in general [[anisotropic]] materials such as [[liquid crystals]]. This allows to [http://www.elis.ugent.be/ELISgroups/lcd/research/bpm.php analyze] e.g. the polarization rotation of light in anisotropic materials, the tunability of a directional coupler based on liquid crystals or the light diffraction in LCD pixels.

==Limitations of BPM==

The Beam Propagation Method relies on the [[slowly varying envelope approximation]], and is inaccurate for the modelling of discretely or fastly varying structures. ItBasic isimplementations are also inaccurate for the modelling of structures in which light propagespropagates in large range of angles and for devices with high refractive-index contrast, commonly found for instance in [[silicon photonics]]. Advanced implementations, however, mitigate some of these limitations allowing BPM to be used to accurately model many of these cases, including many silicon photonics structures.

The BPM method can be used to model bi-directional propagation, but the reflections need to be implemented iteratively which can lead to convergence issues.

More rigorous alternatives include [[Eigenmode Expansion]] and [[FDTD]].

==BPM software==

* RSoft's [http://www.rsoftdesign.com/products.php?sub=Component+Design&itm=BeamPROP BeamPROP]: vector; commercial, free trial possible

* Optiwave's [http://www.optiwave.com/products/bpm_overview.html OptiBPM]: commercial, free trial possible

* [http://www.elis.ugent.be/ELISgroups/lcd/research/bpm.php FEAB] (Finite Element Anisotropic Beam propagation method): academic, free version available

==See also==

Line 65 ⟶ 51:

*[[Finite element method]]

*[[Maxwell's equations]]

*[[MoL|Method of Lineslines]]

*[[Light]]

*[[Photon]]

Line 72 ⟶ 58:

<references/>

[[Category:Computational scienceelectromagnetics]]

[[Category:Electrodynamics]]

[[Category:Physical optics]]