Abstract
Mie theory is a powerful method to model electromagnetic scattering from a multilayered sphere. Usually, the incident beam is expanded to its vector spherical harmonic representation defined by beam shape coefficients, and the multilayer sphere scattering is obtained by the T-matrix method. However, obtaining the beam shape coefficients for arbitrarily shaped incident beams has limitations on source locations and requires different methods when the incident beam is defined inside or outside the computational domain or at the scatterer surface. We propose a 3D angular spectrum method for defining beam shape coefficients from arbitrary source field distributions. This method enables the placement of the sources freely within the computational domain without singularities, allowing flexibility in beam design. We demonstrate incident field synthesis and spherical scattering by comparing morphology-dependent resonances to known values, achieving excellent matching and high accuracy. The proposed method has significant benefits for optical systems and inverse beam design. It allows for the analysis of electromagnetic forward/backward propagation between optical elements and spherical targets using a single method. It is also valuable for optical force beam design and analysis.
Published by Optica Publishing Group under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
1. Introduction
Electromagnetic scattering from a homogeneous sphere illuminated by an arbitrary incident beam can be computed with conventional methods, including full-wave simulations [1], geometrical optics [2], or physical optics [3,4]. These techniques are well-studied and accurate, given the model fidelity and a suitable wavelength range. However, they cannot assess the internal and scattered electric fields from multilayered spherical objects without considerable computational effort. Especially when the sphere’s radius is of the order of wavelengths, complete classical electromagnetic wave theory is needed [5].
Mie theory and the generalized Lorentz-Mie theory are accurate methods to evaluate the internal and scattered fields from the multilayered spheres [6–8], where the incident beam is presented in vector spherical harmonics (VSH) expansion defined by beam shape coefficients (BSCs) [9]. However, obtaining the BSCs for arbitrarily shaped incident beams can be difficult; often needing a combination of several complex methods [6,10]. BSCs can be computed from a known function or field distribution with certain constraints, such as polarization, source shape, and location limitations [11,12]. In these cases, BSCs for the incident beam are obtained by the Bromwich method [13] or multipole expansion [14–16]. These methods define the sources inside a closed volume, and the VSH expanded fields can be computed only outside of this area, limiting the source location.
BSCs can also be computed from known electric and magnetic field distributions on the surface of the spherical scatterer using closed surface orthogonality (CSO). CSO also has limitations as the fields must be defined across the entire spherical surface [17]. For example, CSO cannot be applied to fields defined only on the spherical subregion, thus excluding it from many inverse beam synthesis tasks, where the beam only illuminates a small area of the sphere [18,19]. On the other hand, BSCs can be computed from the standard 2D ASM, which does not have these location restrictions due to its eigenfunction property [20,21]. However, 2D ASM is limited to the planar surface distributions [22–24]; see the comparison of the methods in Table 1.
In this article, we adapt the curved boundary integral method (CBIM) presented in [27] to a source-free method allowing beams to propagate through the boundary surface. We simply generalize the 2D ASM for curved surfaces following the idea presented in the proof of Theorem 2.1 [27] and call this method as 3D angular spectrum method (3D ASM). The geometrical framework is the same as in [27].
We aim to create a VSH-expanded electromagnetic beam from a known electric field distribution on an arbitrarily shaped and positioned surface. Then the incident beam can be modified in terms of propagating field distribution, polarization, and local phase variation [19]. Also, the source location can be positioned inside the simulation area without restrictions. This goal can be achieved by using 3D ASM approach to accommodate arbitrary fields defined on arbitrarily shaped surfaces and construct modified BSCs for an incident beam VSH presentation.
3D ASM present the incident field in the angular spectrum domain as the sum of differently oriented plane waves. When this angular summation of plane waves on a source plane is presented in a direction cosine coordinate system, these angles can be used directly to compute the BSCs for the VSH expansion. Additionally, the internal and scattered fields from a multilayered dielectric sphere can be mapped with the extended boundary condition method (EBCM), which is referred as the T-matrix method in this article [21]. The advantages of the presented method are non-singularity at the source points, allowing the synthesis of a continuous beam through the source surface from any known electric field distribution. Finally, the incident beam can be expanded into a VSH presentation to compute scattered fields from multilayered spheres.
2. Theory
This chapter presents a method to compute the forward/backscatter of multilayered spheres under arbitrary beam illumination. The chapter is divided into three sections introducing the 3D ASM, global coordinate mapping system, and modified BSCs coefficients for the VSH expansion.
Section 2.1 derives the 3D ASM from the curved boundary integral method, which is polarized along local xz-plane and is suitable for obtaining BSCs.
Section 2.2 presents a system to map local coordinates into a global origin-centered coordinate system with each source point’s orientation and position information, which is later used to construct the BSCs.
Section 2.3 presents the VSH expansion from an arbitrary surface with modified BSCs. These BSCs are modified to include the mapped source’s orientation and position information while preserving the needed spherical symmetry for the VSH expansion.
2.1 3D angular spectrum method
Let the global Cartesian base-vectors be ($\mathbf {e}_x,\mathbf {e}_y,\mathbf {e}_z$) in space $\mathbb {R}^3$, where position vector is $\mathbf {r}=x\mathbf {e}_x+y\mathbf {e}_y+z\mathbf {e}_z$ and global coordinates are marked as $\mathbf {r}_\text {glob}=(x,y,z)$. Consider a regular, compact surface $\Omega$ in $\mathbb {R}^3$, which has a continuously differentiable parametrization with parameters $p$ and $q$ as
Keeping point $\mathbf {o}$ as a local origin, unit vectors $\mathbf {e}_1$, $\mathbf {e}_2$ and $\mathbf {e}_3$ define an orthogonal base for the local coordinate system, where $\mathbf {e}_1$ and $\mathbf {e}_2$ are in the tangent plane, and $\mathbf {e}_3$ is along the surface normal, see Fig. 1. Polarization is along $(\mathbf {e}_1,\mathbf {e}_3)-$plane.
Let vector $\mathbf {r}$ be presented in this base as $\mathbf {r}=\bar {x}\mathbf {e}_1+\bar {y}\mathbf {e}_2+\bar {z}\mathbf {e}_3$. Its coordinate triple is marked briefly as $\mathbf {r}_{\text {loc}}=(\bar {x},\bar {y},\bar {z})$. When the vector $\mathbf {r}$ is expressed with both $(\mathbf {e}_x,\mathbf {e}_y,\mathbf {e}_z$) and $(\mathbf {e}_1,\mathbf {e}_2,\mathbf {e}_3$) bases, as it is known, the dependence of the coordinates on each other is determined by an orthogonal transformation matrix
Let us introduce compact notation $\mathbf {\bar {r}}=\mathbf {r}-\mathbf {o}$. The dependence of the corresponding coordinates can be written using the transformation matrix Eq. (2) as
Let $E_0$ be a continuous electric field on the surface $\Omega$. In the article [27], we have presented integral formulas for an incident electric field created by $E_0$. Because our aim in this article is to utilize the Mie theory, we only use frequencies in a disc $\bar {B}(\mathbf {o};k)=\big \{(k_{\bar {x}},k_{\bar {y}}) \ | \ k_{\bar {x}}^2+k_{\bar {y}}^2\leq k^2 \big \}$, where $k$ is the wavenumber. Thus, we ignore some fast-decaying ones. Accordingly, our incident electric field takes a form
Moreover, let’s denote $\mathbf {s}=\frac {1}{k}\mathcal {k}$. Hence the exponent in Eqs. (6)–(7) can be written as $i(k_{\bar {x}}\bar {x}+k_{\bar {y}}\bar {y}+k_{\bar {z}}\bar {z})=i\mathbf {k}\cdot \mathbf {r}=ik(\mathbf {s}\cdot \mathbf {r})$. As in [27], the radiated electric field $\mathbf {E}(\mathbf {r};\mathbf {o})=E_1(\mathbf {r};\mathbf {o})\mathbf {e}_1+E_3(\mathbf {r};\mathbf {o})\mathbf {e}_3$ can be written as
2.2 Global coordinate mapping
Equations (6)–(9) presents a electric field in the base ($\mathbf {e}_1,\mathbf {e}_2,\mathbf {e}_3$). When expanding to VSH presentation, this source must be presented using position vector $\mathbf {r}$ to preserve the spherical symmetry of VSH. Thus, a global mapping of local coordinates is needed.
We obtain according to Eq.(4) that $(\bar {x},\bar {y},\bar {z})=\mathbf {\bar {r}}_{loc}=\Theta ^T(\mathbf {r}-\mathbf {o})_{glob}$, $k_{\bar {z}}=[k^2-k_{\bar {x}}^2-k_{\bar {y}}^2]^{1/2}$. Let’s write briefly $\mathbf {t}=\mathbf {s}_{loc}$ a coordinate triple of vector $\mathbf {s}$. Consequently, we can write
2.3 VSH expansion with modified beam shape coefficients
Modified BSCs coefficients for each source point are derived for presenting any parameterized field distribution in a global coordinate system, including the locations and propagation direction of the sources. These coefficients map complicated coordinate geometries into one presentation. This allows synthesizing a total electromagnetic beam as a superposition of VSH expanded fields.
Integrated wavenumbers $k_{\bar {x}}$ and $k_{\bar {y}}$ in Eq. (8) have direction cosine presentation analogous to the BSC angles [21]. Thus Eq. (8) is expressed in the direction cosines $\xi$ and $\zeta$, where $(k\cos {\xi }=k_{\bar {x}}$ and $k\cos {\zeta }=k_{\bar {y}}$ as
The polarization $\mathbf {e}_0$ is given as
Then, Eq. (12) is evaluated by numerical integration using the trapezoidal rule method, which is sufficiently accurate for periodic function Eq. (12) [28], with uniform-width $l$ [21]
The incident electric field can be written at an arbitrary point $\mathbf {r}$ by VSH expansion where local plane wave angles $(i,j)$ are converted to even and odd VSH modes $(m,n)$ as [21]
The BSCs for each $ij-$ plane wave $a^t_{{\text {emn}}_{ij}}, a^t_{{\text {omn}}_{ij}}, b^t_{{\text {emn}}_{ij}}$ and $b^t_{{\text {omn}}_{ij}}$ are defined in the Appendix. Now, the incident magnetic field is obtained by rearranging the modified BSCs and by multiplying by constant $-i/\eta$ as
The total incident electric and magnetic field, accounting for polarization as in Eq. (13), is written as the integration of vector spherical harmonics expansions Eqs. (17) and (19) as
Similar scattered field presentations mapped with the T-matrix method are presented in the Appendix.
In summary, we have expanded 3D ASM incident field synthesis to the VSH expansion for obtaining Mie scattering from spherical objects. We have mapped source points’ global location and orientation information into the modified BSCs and obtained the total fields as a superposition of VSH-expanded source points. This approach satisfies the Helmholtz equation and allows us to adjust the incident field’s polarization and wavefront from the parametrization.
3. Results
The results are divided into two sections; first, we demonstrate the practicality and accuracy of the presented theory by comparing the Gaussian beam’s morphology-dependent resonances from a sphere with the traditional 2D ASM and the presented 3D ASM, and show a high degree of agreement. In the second section, we synthesize an electromagnetic beam from an ellipsoidal surface, verify the incident field with physical optics simulations, and compute scattered fields from a 100-dielectric layer sphere.
3.1 Morphology-dependent resonance test
Mie scattering region is known for frequency-dependent backscatter intensity. This phenomenon is due to the scattering resonances, which trap energy temporally inside the sphere. Backscatter intensity is typically reported as a function of size parameter $k\alpha$, where $k=k_0n$ is the refractive index dependent wavenumber inside the sphere and $\alpha$ is the sphere radius. These backscatter resonances are called morphology-dependent resonances (MDRs) [29,30]. MDRs are highly responsive to simulation errors and serve as a good test of the accuracy of the derived 3D ASM combined with VSH expansion [31].
To test the 3D ASM simulation accuracy, an ensemble of source points representing the Gaussian beam in [21] was defined via the 2D ASM of [32], and differences in predicted MDRs were assessed. First, a $\text {TEM}_{00}$ beam waist was defined such that the waist was coplanar with $z = 0$ plane and the beam centroid intersected a homogenous sphere center of curvature (plane 1 in Fig. 2) located at $(x, y, z) = (0, 0, 0)$. Then, the beam was backpropagated $5\alpha$ to plane 2 (Fig. 2) via the 2D ASM utilized in [32]. The field at plane 2 was discretized into a square grid of $500\times 500$ points. Source points generated via the 3D ASM were placed at these points, and the field magnitude, phase, polarization, and Poynting vector were defined by the 2D ASM computed field. The fields were then propagated from plane 2 to the sphere via 3D ASM, and backscatter intensity was analyzed at a point on the optical axis $(x, y, z) = (0, 0, -500\alpha )$ as the size parameter range was tuned from 32 to 36. Backscatter at this point was also computed via 2D ASM from plane 1. Finally, the backscatter from a single source point on the axis at plane 2 $(0, 0, -5\alpha )$ was also assessed at $(0, 0, -500\alpha )$ to demonstrate the difference between a single source point and the source point ensemble.
The coated sphere simulations were done in the same way, except that the beam axis was placed tangential to the outer coating surface. The waist (plane 1) was still coplanar with the $z = 0$ plane, but the beam centroid was displaced vertically to $y = \alpha$. This configuration couples to whispering gallery-like modes, and discrepancies in the predicted MDRs are even more sensitive to simulation error than the above-described homogeneous sphere configuration. See simulation flow with source planes 1 and 2 in Fig. 3.
Simulations are executed by illuminating homogeneous and coated spheres by Gaussian beam in size parameter range $k\alpha \in [32,36]$ with beam waist radius $\omega _0=1.5\lambda$. Homogeneous spheres have a refractive index of $n_r=1.36$, and coated spheres have a $n_r=1.36$ core with a $n_c=1.5$ shell with radius of $0.7\ \alpha$. These values were selected from article [32] for comparing results.
The Cartesian electric and magnetic field components of the incident field at plane 2 are presented in Fig. 4. These components are synthesized from plane 1.
Gaussian beam MDRs synthesized from both source plane scenarios are computed at -500 $\alpha$ distance from the origin [21] and compared in Fig. 5.
The resonances are labeled as $TE_{n,l}$ or $TM_{n,l}$, where $n$ indicates the mode number, and $l$ the number of radial peaks in the angle-averaged internal energy density distribution [21]. Electromagnetic backscatter intensities show an excellent match between MDRs in both beam focus scenarios. The resonance shape, location, and peak spacing in Fig. 5 match with values from [32]. In Fig. 5(a), the peaks correspond to modes $TE_{39,1},TE_{40,1}$, and $TE_{41,1}$ respectively, and the spacing between the peaks satisfies the $\Delta k\alpha =\pi /2m_r\approx 1.15$ condition. Also, in Fig. 5(b), the MDRs peaks from the coated sphere correspond to the modes $TM_{39,1},TM_{40,1},TM_{41,1}$ and $TM_{42,1}$ respectively.
The MDRs of an individual source point do not match the Gaussian beam resonances. However, a total beam created as a superposition of the source points from the surface $\Omega$ synthesizes the original beam with highly matching MDRs. This leads to the conclusion that the presented 3D ASM combined with the VSH expansion synthesizes electromagnetic scattering from multilayered spheres with high accuracy.
3.2 Elliptical source distribution
An elliptical surface was used as an example of a source distribution to synthesize a focused linearly polarized incident field using $\lambda /6$ sampling, see Fig. 6. The incident field is compared to the Physical optics simulations at the yz-plane with an excellent agreement to less than a -39 dB difference in amplitude and less than the 0.1-degree difference in phase. Also, the scattered fields are computed from a $\alpha =$ 7.8 mm radius sphere in xz-plane with 100 dielectric layers with linearly varying permittivity $\epsilon _r\in [1-i0.001,\ 3-i0.01]$ from the surface to the center.
Base vectors for elliptical surfaces are derived from the Eqs. (1)–(4) with the parametrization
Figure 7 illustrate the amplitude and phase of the total field $\mathbf {E}_{tot}=\mathbf {E}_{inc}+\mathbf {E}_{scat}$ on xz-plane where $x,z\in [-5a,5a]$. Simulations are performed with $500\times 500$ evaluation points with size parameter of $k\alpha =13$ and VSH modes of $N=70$, see Eq. (26) in the Appendix.
4. Conclusions
We propose a method to compute scattering from multilayered spheres illuminated by incident fields that can be modified freely by adjusting the source distribution’s shape, amplitude, phase, and location. The incident field is presented in a vector spherical harmonic (VSH) expansion defined by the beam shape coefficients (BSCs), and the scattered field is obtained by mapping the incident field BSCs with the T-matrix method.
The goal is achieved by computing BSCs with the 3D angular spectrum method (3D ASM), which we derive by adjusting the curved boundary integral method. The proposed method approximates incident field synthesis from an arbitrary surface, considering the wavelength compared to the surface details. Then, the scattered fields are rigorously computed from approximated incident fields with the Mie theory.
Location and orientation information of 3D ASM sources is transformed in a global coordinate system and mapped into the BSCs to preserve spherical symmetry for VSH expansion. As a result, the total incident field is obtained as a superposition of VSH-expanded source points, in which orientation, locations, and complex amplitude can be adjusted. Each source point affects the total source distribution with an essential point in this approach. If the original $E_0$ is a non-physical distribution, it will be filtered to a physical one, which always satisfies Maxwell’s equations. On the other hand, the final form of the more complex source distributions should be verified by simulations.
Simulations verify the proposed method’s incident field synthesis and scattering accuracy. First, we verify the scattering accuracy by comparing the Morphology-dependent resonances (MDRs) from spheres illuminated by the Gaussian beam created from the beam waist by nominal 2D ASM and by more complex electric field distribution with the proposed 3D ASM. The MDRs resonances of the Gaussian beam synthesized by the 3D ASM were in excellent agreement with the reference values, validating the high accuracy on the scattered fields. Hereafter, we demonstrate the method’s practicality by synthesizing an incident field from an elliptical surface and computing scattered fields from a 100-layer sphere.
The novelty of this approach lies in a straightforward simulation algorithm, where the incident field can be defined at any parametrized surface. Furthermore, the parameterized surface can be located inside the computational domain, on the surface, or inside the spherical scatterer without restrictions and source singularities. Unrestricted placement of the source distribution has clear advantages, especially in the inverse beam design, where the desired beam is defined on the sub-region on the spherical scatterer. Additionally, the synthesized beam can be reradiated from any evaluation surface to another, enabling beam simulation between optical elements with the ability to consistently compute the scattered fields from multilayered spheres.
Appendix
The VSH beam shape coefficients for each $ij-$ plane wave of the incident field is defined as [33]
The required number of the VSH modes ($N_{req}$) for the series to convergence is defined as a function of size parameter $x_l=kr$ and $N_{req}=max(N_{stop},|n_lx_l|,|n_lx_{l-1}|)+15$ when $l=1,2,..,L$ and [22]
Funding
Agencia Estatal de Investigación (PID2019-107885GB282 C31/AEI/10.13039, PRE2018-084326); Assessment of the Graft 281 Rejection Using Millimeter Waves; Academy of Finland (327640).
Acknowledgment
The authors would like to thank Dr. Lasse Lamberg for his mathematical support for the 3D ASM theory.
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. E. L. Roetman, R. P. Kochhar, and G. L. Hower, “An integral representation for the fields in electromagnetic scattering problems,” Electromagn. 12(1), 1–15 (1992). [CrossRef]
2. M. Gharghi, H. Bai, G. Stevens, et al., “Modeling and simulation of spherical solar cells,” Conference Record of the Thirty-first IEEE Photovoltaic Specialists Conference, (2005). 1177–1180.
3. N. Keerativoranan, P. Hanpinitsak, K. Saito, et al., “Circular mesh based physical optics for scattered field prediction,” Antennas Wirel. Propag. Lett. 20(3), 391–395 (2021). [CrossRef]
4. J. R. Sanford, “Scattering by spherically stratified microwave lens antennas,” IEEE Trans. Antennas Propag. 42(5), 690–698 (1994). [CrossRef]
5. A. A. R. Neves and C. L. Cesar, “Analytical calculation of optical forces on spherical particles in optical tweezers: tutorial,” J. Opt. Soc. Am. B 36(6), 1525–1537 (2019). [CrossRef]
6. G. Gouesbet, J. Lock, and G. Gréhan, “Generalized lorenz-mie theories and description of electromagnetic arbitrary shaped beams: Localized approximations and localized beam models, a review,” J. Quant. Spectrosc. Radiat. Transfer 112(1), 1–27 (2011). [CrossRef]
7. G. Gouesbet, “A scientific story of generalized lorenz-mie theories with epistemological remarks,” J. Quant. Spectrosc. Radiat. Transfer 126, 7–15 (2013). [CrossRef]
8. Z. S. Wu, L. X. Guo, K. F. Ren, et al., “Improved algorithm for electromagnetic scattering of plane waves and shaped beams by multilayered spheres,” Appl. Opt. 36(21), 5188–5198 (1997). [CrossRef]
9. C. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley Science Paperback Series, 1998).
10. G. Gouesbet, G. Gréhan, and B. Maheu, “Localized interpretation to compute all the coefficients gnm in the generalized lorenz–mie theory,” J. Opt. Soc. Am. A 7(6), 998–1007 (1990). [CrossRef]
11. F. Onofri, G. Gréhan, and G. Gouesbet, “Electromagnetic scattering from a multilayered sphere located in an arbitrary beam,” Appl. Opt. 34(30), 7113–7124 (1995). [CrossRef]
12. L. Mees, K. F. Ren, G. Gréhan, et al., “Scattering of a gaussian beam by an infinite cylinder with arbitrary location and arbitrary orientation: numerical results,” Appl. Opt. 38(9), 1867–1876 (1999). [CrossRef]
13. B. Maheu, G. Gouesbet, and G. Grehan, “A concise presentation of the generalized lorenz-mie theory for arbitrary location of the scatterer in an arbitrary incident profile,” J. Opt. 19(2), 59–67 (1988). [CrossRef]
14. J. A. Lock, “Partial-wave expansions of angular spectra of plane waves,” J. Opt. Soc. Am. A 23(11), 2803–2809 (2006). [CrossRef]
15. A. V. Zvyagin and K. Goto, “Mie scattering of evanescent waves by a dielectric sphere: comparison of multipole expansion and group-theory methods,” J. Opt. Soc. Am. A 15(12), 3003–3008 (1998). [CrossRef]
16. J. M. Taylor and G. D. Love, “Multipole expansion of bessel and gaussian beams for mie scattering calculations,” J. Opt. Soc. Am. A 26(2), 278–282 (2009). [CrossRef]
17. G. Gouesbet, L. A. Ambrosio, and J. A. Lock, “On an infinite number of quadratures to evaluate beam shape coefficients in generalized lorenz-mie theory and the extended boundary condition method for structured em beams,” J. Quant. Spectrosc. Radiat. Transfer 242, 106779 (2020). [CrossRef]
18. J. Lamberg, F. Zarrinkhat, A. Tamminen, et al., “Investigation of radial top-hat electric field distributions for corneal reflectometry using modified fourier optics method,” in 46th International Conference on Infrared, Millimeter and Terahertz Waves (IRMMW-THz), (2021), pp. 1–2.
19. J. Lamberg, F. Zarrinkhat, A. Tamminen, et al., “Wavefront-modified vector beams for thz cornea spectroscopy,” Opt. Express (2023).
20. J. W. Goodman, Introduction to Fourier optics (Roberts, Englewood, CO, 2005).
21. E. E. M. Khaled, S. C. Hill, and P. W. Barber, “Scattered and internal intensity of a sphere illuminated with a gaussian beam,” IEEE Trans. Antennas Propag. 41(3), 295–303 (1993). [CrossRef]
22. F. Zarrinkhat, J. Lamberg, A. Tamminen, et al., “Vector spherical harmonic analysis and experimental validation of spherical shells illuminated with broadband, millimeter wave gaussian beams: applications to corneal sensing,” Biomed. Opt. Express 13(7), 3699–3722 (2022). [CrossRef]
23. F. Zarrinkhat, M. Baggio, J. Lamberg, et al., “Calibration alignment sensitivity in corneal terahertz imaging,” Sensors 22(9), 3237 (2022). [CrossRef]
24. F. Zarrinkhat, J. Lamberg, A. Tamminen, et al., “Fourier analysis of submillimeter-wave scattering from the human cornea,” in 15th European Conference on Antennas and Propagation (EuCAP), (2021), pp. 1–5.
25. C. Y. Hwang, S. Oh, I. K. Jeong, et al., “Stepwise angular spectrum method for curved surface diffraction,” Opt. Express 22(10), 12659–12667 (2014). [CrossRef]
26. L. Ebers, M. Hammer, and J. Förstner, “Light diffraction in slab waveguide lenses simulated with the stepwise angular spectrum method,” Opt. Express 28(24), 36361–36379 (2020). [CrossRef]
27. J. Lamberg, F. Zarrinkhat, A. Tamminen, et al., “Curved boundary integral method for electromagnetic fields,” Opt. Express (2023).
28. L. N. Trefethen and J. A. C. Weideman, “The exponentially convergent trapezoidal rule,” SIAM Rev. 56(3), 385–458 (2014). [CrossRef]
29. B. R. Johnson, “Theory of morphology-dependent resonances: shape resonances and width formulas,” J. Opt. Soc. Am. A 10(2), 343–352 (1993). [CrossRef]
30. D. Q. Chowdhury, S. C. Hill, and P. W. Barber, “Morphology-dependent resonances in radially inhomogeneous spheres,” J. Opt. Soc. Am. A 8(11), 1702–1705 (1991). [CrossRef]
31. J. A. Lock, “Excitation efficiency of a morphology-dependent resonance by a focused gaussian beam,” J. Opt. Soc. Am. A 15(12), 2986–2994 (1998). [CrossRef]
32. E. E. M. Khaled, S. C. Hill, and P. W. Barber, “Light scattering by a coated sphere illuminated with a gaussian beam,” Appl. Opt. 33(15), 3308–3314 (1994). [CrossRef]
33. P. W. Barber and S. C. Hill, Light Scattering By Particles: Computational Methods (Advanced Series in Applied Physics, 1990).
34. U. P. O. Peña, “Scattering of electromagnetic radiation by a multilayered sphere,” IEEE Antennas and Propagation Magazine 57, 69–116 (2015).