## Abstract

The fast and accurate propagation of general optical fields in free space is still a challenging task. Most of the standard algorithms are either fast or accurate. In the paper we introduce a new algorithm for the fast propagation of non-paraxial vectorial optical fields without further physical approximations. The method is based on decomposing highly divergent (non-paraxial) fields into subfields with small divergence. These subfields can then be propagated by a new semi-analytical spectrum of plane waves (SPW) operator using fast Fourier transformations. In the target plane, all propagated subfields are added coherently. Compared to the standard SPW operator, the numerical effort is reduced drastically due to the analytical handling of linear phase terms arising after the decomposition of the fields. Numerical results are presented for two examples demonstrating the efficiency and the accuracy of the new method.

© 2012 OSA

Full Article |

PDF Article
### Equations (25)

Equations on this page are rendered with MathJax. Learn more.

(1)
$${A}_{\ell}\left(\mathit{\kappa},0\right)=\mathcal{F}\left[{V}_{\ell}\left(\mathit{\rho},0\right)\right]=\frac{1}{2\pi}{\iint}_{-\infty}^{\infty}{V}_{\ell}\left(\mathit{\rho},0\right){\text{e}}^{-\text{i}\mathit{\kappa}\cdot \mathit{\rho}}\text{d}\mathit{\rho}.$$
(2)
$${V}_{\ell}\left(\mathit{\rho},z\right)={\mathcal{F}}^{-1}\left[{A}_{\ell}\left(\mathit{\kappa},0\right){\text{e}}^{\text{i}{k}_{z}(\mathit{\kappa})z}\right]$$
(3)
$$\mathit{k}=k{\widehat{\mathit{s}}}_{0}=k{\left({s}_{0x},{s}_{0y},{s}_{0z}\right)}^{\text{T}}=k{\left(\text{sin}{\theta}_{0}\text{cos}{\varphi}_{0},\text{sin}{\theta}_{0}\text{sin}{\varphi}_{0},\text{cos}{\theta}_{0}\right)}^{\text{T}}$$
(4)
$${\mathit{\kappa}}_{0}={\left({k}_{0x},{k}_{0y}\right)}^{\text{T}}=k{\left({s}_{0x},{s}_{0y}\right)}^{\text{T}}$$
(5)
$${V}_{\ell}\left(\mathit{\rho},0\right)={V}_{\ell}^{\prime}\left(\mathit{\rho},0\right){\text{e}}^{\text{i}{\mathit{\kappa}}_{0,\ell}\cdot \mathit{\rho}}$$
(6)
$${A}_{\ell}\left(\mathit{\kappa},0\right)={A}_{\ell}^{\prime}\left(\mathit{\kappa}-{\mathit{\kappa}}_{0,\ell},0\right)$$
(7)
$${V}_{\ell}^{\prime}\left(\mathit{\rho},0\right)={\mathcal{F}}^{-1}\left[{A}_{\ell}^{\prime}\left(\mathit{\kappa},0\right)\right].$$
(8)
$${k}_{z,\ell}(\mathit{\kappa})={\gamma}_{0,\ell}+{\mathit{\gamma}}_{1,\ell}\cdot \mathit{\kappa}+{\gamma}_{\ell}\left(\mathit{\kappa}-{\mathit{\kappa}}_{0,\ell}\right)$$
(9)
$${\gamma}_{0,\ell}={k}_{0z,\ell}+\frac{1}{{k}_{0z,\ell}}\left({k}_{0x,\ell}^{2}+{k}_{0y,\ell}^{2}\right)-\frac{{k}_{0x,\ell}^{2}{k}_{0y,\ell}^{2}}{{k}_{0z,\ell}^{3}},$$
(10)
$${\mathit{\gamma}}_{1,\ell}={\left[\frac{{k}_{0x,\ell}}{{k}_{0z,\ell}}\left(\frac{{k}_{0y,\ell}^{2}}{{k}_{0z,\ell}^{2}}-1\right),\frac{{k}_{0y,\ell}}{{k}_{0z,\ell}}\left(\frac{{k}_{0x,\ell}^{2}}{{k}_{0z,\ell}^{2}}-1\right)\right]}^{\text{T}}$$
(11)
$${k}_{0z,\ell}=\sqrt{{k}^{2}-{k}_{0x,\ell}^{2}-{k}_{0y,\ell}^{2}}.$$
(12)
$${\gamma}_{\ell}\left({\mathit{\kappa}}^{\prime}\right)={k}_{z,\ell}\left({\mathit{\kappa}}^{\prime}+{\mathit{\kappa}}_{0,\ell}\right)-{\gamma}_{0,\ell}-{\mathit{\gamma}}_{1,\ell}\cdot \left({\mathit{\kappa}}^{\prime}+{\mathit{\kappa}}_{0,\ell}\right)$$
(13)
$${k}_{z,\ell}\left({\mathit{\kappa}}^{\prime}+{\mathit{\kappa}}_{0,\ell}\right)=\sqrt{{k}^{2}-{\left({k}_{x}^{\prime}+{k}_{0x,\ell}\right)}^{2}-{\left({k}_{y}^{\prime}+{k}_{0y,\ell}\right)}^{2}}.$$
(14)
$${V}_{\ell}\left(\mathit{\rho},z\right)={V}_{\ell}^{\prime}\left(\mathit{\rho}+z{\mathit{\gamma}}_{1,\ell},z\right)\text{exp}\left[\text{i}\left({\gamma}_{0,\ell}z+{\mathit{\gamma}}_{1,\ell}\cdot {\mathit{\kappa}}_{0,\ell}z+\mathit{\rho}\cdot {\mathit{\kappa}}_{0,\ell}\right)\right]$$
(15)
$${V}_{\ell}^{\prime}\left(\mathit{\rho},z\right)={\mathcal{F}}^{-1}\left[{A}_{\ell}^{\prime}\left(\mathit{\kappa},0\right){\text{e}}^{\text{i}{\gamma}_{\ell}(\mathit{\kappa})z}\right]$$
(16)
$$d:=\frac{{\sum}_{x,y}{\left|{V}_{\ell ,\text{operator}}\left(x,y,z\right)-{V}_{\ell ,\text{reference}}\left(x,y,z\right)\right|}^{2}}{{\sum}_{x,y}{\left|{V}_{\ell ,\text{reference}}\left(x,y,z\right)\right|}^{2}}$$
(17)
$$\eta :=\frac{\text{Complexity}\hspace{0.17em}\text{of}\hspace{0.17em}\text{SPW}\hspace{0.17em}\text{algorithm}}{\text{Complexity}\hspace{0.17em}\text{of}\hspace{0.17em}\text{semi-analytical}\hspace{0.17em}\text{SPW}\hspace{0.17em}\text{algorithm}}$$
(18)
$${A}_{\ell}\left(\mathit{\kappa},0\right)=\sum _{i,j=1}^{N}{\mathrm{\Phi}}_{i,j}(\mathit{\kappa}){A}_{\ell}\left(\mathit{\kappa},0\right)$$
(19)
$$\begin{array}{ll}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathrm{\Phi}\left({k}_{x}\right)=0.5\left[\text{sin}\left(\frac{\pi}{a}\left({k}_{x}-\frac{a}{2}\right)\right)+1\right]\hfill & \text{for}\hspace{0.17em}0<{k}_{x}\le a,\hfill \\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathrm{\Phi}\left({k}_{x}\right)=1\hfill & \text{for}\hspace{0.17em}a<{k}_{x}<b-a,\hfill \\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathrm{\Phi}\left({k}_{x}\right)=1-0.5\left[\text{sin}\left(\frac{\pi}{a}\left({k}_{x}-\frac{2b-a}{2}\right)\right)+1\right]\hfill & \text{for}\hspace{0.17em}b-a\le {k}_{x}<b,\hfill \\ \text{and}\hspace{0.17em}\mathrm{\Phi}\left({k}_{x}\right)=0\hfill & \text{anywhere}\hspace{0.17em}\text{else},\hfill \end{array}$$
(20)
$${\mathrm{\Phi}}_{i,j}(\mathit{\kappa})=\mathrm{\Phi}\left({k}_{x}-{k}_{i,x}\right)\mathrm{\Phi}\left({k}_{y}-{k}_{j,y}\right)$$
(21)
$${V}_{\ell}\left(\mathit{\rho},0\right)={\widehat{V}}_{\ell}\left(\mathit{\rho},0\right){\text{e}}^{\text{i}{\phi}_{\ell}\left(\mathit{\rho},0\right)}$$
(22)
$${\phi \left(\mathit{\rho},0\right)|}_{\left(\mathit{\rho}\approx {\mathit{\rho}}_{0,\ell}\right)}={\epsilon}_{\ell}+{\mathit{\kappa}}_{0,\ell}\cdot \mathit{\rho}+{\mathrm{\Delta}}_{\ell}\left(\mathit{\rho},0\right)$$
(23)
$${\epsilon}_{\ell}={\phi}_{\ell}\left({\mathit{\rho}}_{0,\ell},0\right)+{\frac{{\partial}^{2}{\phi}_{\ell}}{\partial x\partial y}|}_{{\mathit{\rho}}_{0,\ell}}\cdot \left({x}_{0,\ell}{y}_{0,\ell}\right)-{\frac{\partial {\phi}_{\ell}}{\partial x}|}_{{\mathit{\rho}}_{0,\ell}}\cdot {x}_{0,\ell}-{\frac{\partial {\phi}_{\ell}}{\partial y}|}_{{\mathit{\rho}}_{0,\ell}}\cdot {y}_{0,\ell}$$
(24)
$${\mathit{\kappa}}_{0,\ell}={\left({k}_{0x,\ell},{k}_{0y,\ell}\right)}^{\text{T}}={\left({\frac{\partial {\phi}_{\ell}}{\partial x}|}_{{\mathit{\rho}}_{0,\ell}}-{\frac{{\partial}^{2}{\phi}_{\ell}}{\partial x\partial y}|}_{{\mathit{\rho}}_{0,\ell}}\cdot {y}_{0,\ell},{\frac{\partial {\phi}_{\ell}}{\partial y}|}_{{\mathit{\rho}}_{0,\ell}}-{\frac{{\partial}^{2}{\phi}_{\ell}}{\partial x\partial y}|}_{{\mathit{\rho}}_{0,\ell}}\cdot {x}_{0,\ell}\right)}^{\text{T}}.$$
(25)
$${V}_{\ell}\left(\mathit{\rho},0\right)={\widehat{V}}_{\ell}\left(\mathit{\rho},0\right)\text{exp}\left[\text{i}({\epsilon}_{\ell}+{\mathrm{\Delta}}_{\ell}(\mathit{\rho}-{\mathit{\rho}}_{0,\ell},0)\right]\text{exp}\left[\text{i}\left(\mathit{\rho}\cdot {\mathit{\kappa}}_{0,\ell}\right)\right]$$