Abstract
The resonant-state expansion, a recently developed powerful method in electrodynamics, is generalized here for open optical systems containing magnetic, chiral, or bi-anisotropic materials. It is shown that the key matrix eigenvalue equation of the method remains the same, but the matrix elements of the perturbation now contain variations of the permittivity, permeability, and bi-anisotropy tensors. A general normalization of resonant states in terms of the electric and magnetic fields is presented.
Published by The Optical Society 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.
The resonant state expansion (RSE) is a novel powerful theoretical method that has been recently developed in electrodynamics [1]. The RSE is a rigorous perturbation theory that is not limited to small perturbations and warrants an efficient calculation of all resonant states (RSs) of an open optical system in an arbitrarily selected spectral range. This calculation is based on knowing the RSs of another, so-called basis system, which is usually (but not necessarily) simpler than the system of interest, ideally having an exact analytic solution. The RSE was verified and tested on optical systems of different shape and dimensionality [2–5], demonstrating its superior computational efficiency [3,5] compared to available numerical methods, such as finite-difference time-domain [6,7], finite element [8], and the aperiodic Fourier modal method [9,10].
Being originally introduced in nuclear physics almost a century ago [11,12], RSs in electrodynamics present eigensolutions of Maxwell’s equations satisfying outgoing boundary conditions, which correspond to electromagnetic excitations decaying in time, with the electromagnetic energy leaking out of the system. This leakage, however, causes an exponential growth of the RS wave function with distance, so that the standard normalization used, for example, for bound states in quantum mechanics or for waveguide modes in optics, diverges. While the correct normalization for scalar fields was known [13], expressions for the normalization of the electromagnetic vector fields of the RSs, intensively used in the literature [14,15], are only approximate, as has been recently clarified [16,17]. The correct normalization of RSs in finite optical systems, which is a cornerstone of the RSE, was presented in the very first paper on the method [1] and was later generalized to arbitrary systems with frequency dispersion of the permittivity [18]. Recently, it has been used to formulate an exact theory of the Purcell effect [16]. Furthermore, the exact normalization was extended to photonic crystal structures [19,20] and applied to resonantly enhanced refractive index sensing using the RSE with only one and two RSs in the basis.
The RSE has also been generalized to optical systems with frequency dispersion of the permittivity [18] without affecting the computational complexity, which is a very important step towards describing realistic materials and specifically plasmonic effects. This was achieved by treating the dispersion as an analytical function with a finite number of simple poles in the lower half-plane of the complex frequency, known in the literature as the generalized Drude–Lorentz model [21].
So far, the RSE was applied to non-magnetic optical systems (), which are fully described by the wave equation containing solely the electric field, the permittivity tensor, and an electric current. Naturally, the RS normalization and the RSE itself dealing with perturbations of the permittivity were formulated in terms of the electric field only. However, the most general materials with local responses are bi-anistropic and have non-zero magnetic susceptibility and coupling tensors between the electric and magnetic fields, including the chiral optical activity and circular dichroism [22]. Describing such systems, which include but are not limited to metamaterials [23], chiral plasmonics [24,25], and chiral sensors, is of growing interest. Therefore, it is crucial to have a general formulation of the RSE and the RS normalization, in which the electric and magnetic fields contribute as equal partners, and the local linear response of an optical system is taken in the most general form. This is done in the this Letter below by using a novel approach for deriving the RSE that is based on Maxwell’s equations.
Maxwell’s equations and Green’s dyadic: An arbitrary linear optical system is described by Maxwell’s equations in a medium:
where is the wave number in vacuum, and is the frequency of the electromagnetic field. Generally, for systems with a spatially local linear response, one can write with frequency dependent tensors of permittivity and permeability , and bi-anisotropy tensors and . In the following, we concentrate on systems satisfying the reciprocity relation, leading additionally to , , and . Equations (1) and (2) can be written in the following compact symmetric way: with matrix operator , where and . The electric and magnetic fields, as well as the currents, are now represented by six-dimensional vectors: respectively, where , and the magnetic current is introduced artificially for symmetry purposes.We now introduce a generalized dyadic Green’s function (GF) with outgoing boundary conditions in the regions outside the optical system, satisfying the equation
in which is the identity matrix. The GF has simple poles [2,3] at , which are the wave numbers of the RSs of the system. The RSs, in turn, are the eigensolutions of the homogeneous Maxwell’s equations: satisfying outgoing boundary conditions, where the index is used to label the RSs. Owing to the reciprocity principle and the Mittag-Leffler (ML) theorem, the GF is represented as a series [2]: determining the normalization of RSs that is considered below. Note that Eq. (8) is valid within the system, or rather within a minimal convex volume including it.Closure relation and sum rules: Substituting the ML expansion Eq. (8) into Eq. (6) for the GF and using Eq. (7), we obtain
In the absence of dispersion, Eq. (9) immediately results in the following closure relation: In the case of a frequency dispersion described by a generalized Drude–Lorentz model [18,21], the matrix becomes having complex poles at with generalized conductivities playing the role of the residues in the ML expansion of the permittivity, permeability, and bi-anisotropy tensors. Examples of such an expansion for plasmonic and chiral materials are given in Refs. [18] and [24], respectively. Substituting Eq. (11) into Eq. (9) and using the algebraic identity yields Since the Lorentzian functions are linearly independent, Eq. (13) splits into sum rules and a closure relation similar to the non-dispersive one, Eq. (10). Summing Eq. (14) over all and adding it to Eq. (15), we can reformulate the closure relation asNormalization of resonant states: As already mentioned, the form of the GF Eq. (8) determines the normalization of the RS wave functions . To derive this normalization, we introduce an analytic continuation of RS field in the complex -plane around the point. satisfies Eq. (3), which can be solved with the help of the GF. Using Eq. (8), we obtain
The requirement that in the limit results in the following dependence of the current: , where is normalized such that While is an arbitrary function of and , vanishing outside the system, it can be chosen to be independent of . Equation (18) then provides the normalization of the RSs. Indeed, multiplying Eq. (3) with , Eq. (7) with , and taking the difference between the two, yields where the and dependencies are dropped for brevity of notations. The third and the fourth terms in the left-hand side of Eq. (19) can be written as Integrating Eq. (19) over an arbitrary volume containing the system, using the divergence theorem, and taking the limit , we obtain a general formula for the RS normalization: where is the boundary of , and the prime means the derivative with respect to taken at . The differentiation of the matrix is straightforward, whereas the derivatives of the analytic continuation of the fields outside the system can be expressed as [1,16] in which is either or . The normalization Eq. (20) then takes an explicit form in terms of the electric and magnetic fields of a given RS: This general normalization is fully consistent with the analytic normalizations we have previously used in terms of the electric field [1,3,16,18,19], for systems described by the permittivity, as we demonstrate below. We note, however, that writing the GF as in Eq. (8), the electric field of the normalized RS is a factor of smaller than the one used in our previous works. Furthermore, as we also show below, the general normalization Eq. (22) is suited for both static and non-static RSs, which is consistent with two different expressions used previously for these cases [3].Let us show that for non-magnetic materials, described by only the permittivity, the general normalization Eq. (22) reduces to the one previously used in terms of the electric field only [1,3,16,18,19]. In this case, and , where is a identity matrix, and Eq. (22) becomes
where we have taken the surface term again in the form of the field derivatives, as in Eq. (20). Using the Poynting theorem for the RS wave function, we can transform the second volume integral in Eq. (23) into For the surface integral in Eq. (23), we obtain using vector identities, as well as and the fact that outside the system. Collecting all terms, we obtain the normalization condition for RSs with : where means the spatial derivative along the surface normal, and , according to Eq. (21).For static electric modes with , the condition leads to the volume term in Eq. (23) with the magnetic field vanishing. Since the electric field of a static mode far away from the system [3] and the surface of integration can be chosen as any closed surface including the system, one can get rid of the surface integral, ending up with the volume integral of the electric field over the entire space:
Both results Eq. (26) and Eq. (27) are identical to the normalization of RSs in non-magnetic materials obtained in [1,3,16,18], with the already noted factor of 2 introduced in this Letter.Resonant-state expansion: Let us now consider a perturbed system described by a general frequency dependent perturbation of the permittivity, permeability, and bi-anisotropy tensors. The Maxwell equation for a perturbed RS then takes the form
where is the perturbed eigenvalue. Note that the unperturbed system and the perturbation are chosen in such a way that the perturbation is included in the minimal convex volume containing the unperturbed system. Solving Eq. (28) with the help of the GF, we obtain Let us first assume a non-dispersive perturbation . Substituting the ML expansion Eq. (8) into Eq. (29) and expanding the perturbed field inside the system into the unperturbed RSs as we obtain where the matrix elements of the perturbation are given by Equating coefficients at the basis functions , Eq. (31) reduces to which is the standard non-dispersive RSE equation [1,3].Taking into account the dispersion of the perturbation in a generalized Drude–Lorentz form,
Eq. (29) becomes where we have added in the second line zeros in the form of the sum rules defined by Eq. (14): Using again the ML expansion Eq. (8) of the GF and the algebraic identity Eq. (12), we arrive, after equating coefficients at the basis functions , at the linear eigenvalue equation of the dispersive RSE: with the matrix elements of the dispersive perturbation given by Note that Eq. (37) has the same form as that developed in Ref. [18] and, in the case of no frequency dispersion, it reduces back to Eq. (33). However, the matrix elements Eq. (38) now have the most general form, which can be written explicitly as The matrix elements Eq. (39) are expressed in terms of the electric and magnetic fields of basis RSs and and generally dispersive changes of the tensors of the permittivity , permeability , and bi-anisotropy couplings and between the electric and magnetic fields. Solving the matrix eigenvalue problem Eq. (37) of the RSE determines the wave numbers of the perturbed RSs and the coefficients of the expansion of the perturbed wave functions into the known RSs of a basis system. Presently, this is the most efficient and intuitive computational approach for finding the RSs of open optical systems, as demonstrated in numerous publications [1–5,18–20]. This approach is now generalized to bi-anisotropic systems.In conclusion, we have generalized the resonant-state expansion for open optical systems containing arbitrary reciprocal bi-anisotropic materials or metamaterials, including those having magnetic and chiral optical activity, as well as circular dichroism. We have presented the theory in the most general, compact, and symmetrized way, with the electric and magnetic field vectors contributing on equal footing. We have addressed both cases of non-dispersive systems and systems having frequency dispersion described by a generalized Drude-Lorentz model. We have derived a general compact formula for the normalization of resonant states, expressed in terms of the electric and magnetic fields, and shown its equivalence to the one used previously for systems fully described by their permittivity and expressed in terms of the electric field only. The presented theory has the widest spectrum of applications, ranging from the modeling and optimization of chirality sensors to the accurate description of the optics of magnetic and metamaterial systems.
Funding
Engineering and Physical Sciences Research Council (EPSRC) (EP/M020479/1); Sêr Cymru National Research Network in Advanced Engineering and Materials; Russian Foundation for Basic Research (RFBR) (16-29-03283); Deutsche Forschungsgemeinschaft (DFG) (DFG SPP 1839); VW Foundation; Ministerium für Wissenschaft, Forschung und Kunst Baden-Württemberg (MWK).
Acknowledgment
E. A. Muljarov acknowledges discussions with W. Langbein.
REFERENCES
1. E. A. Muljarov, W. Langbein, and R. Zimmermann, Europhys. Lett. 92, 50010 (2010). [CrossRef]
2. M. B. Doost, W. Langbein, and E. A. Muljarov, Phys. Rev. A 87, 043827 (2013). [CrossRef]
3. M. B. Doost, W. Langbein, and E. A. Muljarov, Phys. Rev. A 90, 013834 (2014). [CrossRef]
4. L. J. Armitage, M. B. Doost, W. Langbein, and E. A. Muljarov, Phys. Rev. A 89, 053832 (2014). [CrossRef]
5. S. V. Lobanov, G. Zoriniants, W. Langbein, and E. A. Muljarov, Phys. Rev. A 95, 053848 (2017). [CrossRef]
6. K. S. Yee, IEEE Trans. Antennas Propag. 14, 302 (1966). [CrossRef]
7. K. S. Kunz and R. J. Luebbers, The Finite Difference Time Domain Method for Electromagnetics (CRC Press, 1993).
8. G. Dhatt, G. Touzot, and E. Lefrançois, Finite Element Method (ISTE Ltd, 2012).
9. P. Lalanne and E. Silberstein, Opt. Lett. 25, 1092 (2000). [CrossRef]
10. M. Pisarenco, J. Maubach, I. Setija, and R. Mattheij, J. Opt. Soc. Am. A 27, 2423 (2010). [CrossRef]
11. G. Gamow, Zeitschrift für Physik 51, 204 (1928). [CrossRef]
12. A. J. F. Siegert, Phys. Rev. 56, 750 (1939). [CrossRef]
13. E. E. Shnol, Theor. Math. Phys. 8, 729 (1971). [CrossRef]
14. P. T. Leung and K. M. Pang, J. Opt. Soc. Am. B 13, 805 (1996). [CrossRef]
15. P. T. Kristensen, C. V. Vlack, and S. Hughes, Opt. Lett. 37, 1649 (2012). [CrossRef]
16. E. A. Muljarov and W. Langbein, Phys. Rev. B 94, 235438 (2016). [CrossRef]
17. E. A. Muljarov and W. Langbein, Phys. Rev. A 96, 017801 (2017). [CrossRef]
18. E. A. Muljarov and W. Langbein, Phys. Rev. B 93, 075417 (2016). [CrossRef]
19. T. Weiss, M. Mesch, M. Schäferling, H. Giessen, W. Langbein, and E. A. Muljarov, Phys. Rev. Lett. 116, 237401 (2016). [CrossRef]
20. T. Weiss, M. Schäferling, H. Giessen, N. A. Gippius, S. G. Tikhodeev, W. Langbein, and E. A. Muljarov, Phys. Rev. B 96, 045129 (2017). [CrossRef]
21. H. S. Sehmi, W. Langbein, and E. A. Muljarov, Phys. Rev. B 95, 115444 (2017). [CrossRef]
22. I. Lindell, A. Sihvola, S. Tretyakov, and A. Viitanen, Electromagnetic Waves in Chiral and Bi-Isotropic Media (Artech House, 1994).
23. C. E. Kriegler, M. S. Rill, S. Linden, and M. Wegener, IEEE J. Sel. Top. Quantum Electron. 16, 367 (2010). [CrossRef]
24. A. O. Govorov and Z. Fan, Chem. Phys. Chem. 13, 2551 (2012). [CrossRef]
25. M. L. Nesterov, X. Yin, M. Schäferling, H. Giessen, and T. Weiss, ACS Photon. 3, 578 (2016). [CrossRef]