## Calculation, normalization and perturbation of quasinormal modes in coupled cavity-waveguide systems

Philip Trøst Kristensen, 1 2 , ∗ Jakob Rosenkrantz de Lasson, 1 and Niels Gregersen 1

1 DTU Fotonik, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark 2 Institut f¨r Physik, Humboldt Universit¨t zu Berlin, D-12489 Berlin, Germany u a ∗ Corresponding author: philip.kristensen@physik.hu-berlin.de

## Compiled August 5, 2014

We show how one can use a non-local boundary condition, which is compatible with standard frequency domain methods, for numerical calculation of quasinormal modes in optical cavities coupled to waveguides. In addition, we extend the definition of the quasinormal mode norm by use of the theory of divergent series to provide a framework for modeling of optical phenomena in such coupled cavity-waveguide systems. As an example, we apply the framework to study perturbative changes in the resonance frequency and Q value of a photonic crystal cavity coupled to a defect waveguide.

Quasinormal modes (QNMs) [1-3] provide a natural framework for modeling light propagation and lightmatter interaction in resonant electromagnetic material systems, such as optical cavities [4-10] and plasmonic nanoparticles [11-15]. The QNMs of localized electromagnetic resonators in an otherwise homogeneous medium may be calculated as eigenmodes of the sourcefree Maxwell equations augmented with the SilverM˜ A ¼ ller radiation condition [3]. The radiation condition admits only fields that propagate away from the resonator at large distances. This, in turn, leads to a discrete spectrum of complex resonance frequencies ˜ ω µ = ω µ -i γ µ from which the Q value may be calculated directly as Q = ω / γ µ 2 µ [16]. In numerical calculations, the radiation condition is often modeled using perfectly matched layers (PMLs) [17] in either frequency- or timedomain calculations, but alternatives are also available, for example in the form of integral equations [9,12].

In this Letter we present an alternative to the FMM based on a non-local boundary condition, which can be used with standard frequency domain methods to calculate QNMs in coupled cavity-waveguide systems. Figure 1 shows the QNM of a side-coupled cavity in a PC with lattice constant a made from high-index cylinders ( glyph[epsilon1] rods = 8 9) in air. This QNM - which is identical to . the QNM that was calculated with the FMM in Ref. [22] -wascalculated using the finite element method (FEM) with the special boundary condition to properly model the waveguide radiation.

In many technologically relevant applications, light coupling to and from the cavity is controlled by waveguides, such as fibers or line defects in photonic crystal (PC) [18] circuits. In such coupled cavity-waveguide systems, the coupling to the waveguides represents by far the largest decay channel for light in the cavity; indeed, any decay through other channels than the waveguide often represents unwanted and detrimental losses and is typically minimized through careful engineering. Calculation of QNMs in these coupled systems is non-trivial because of the need for a suitable radiation condition. In particular, the Silver-M ˜ A ller radiation condition applies only ¼ to problems with a homogeneous background material, and if the waveguides have discrete translational symmetry - as in PCs, for example - one cannot use PMLs to avoid reflections from the calculation domain boundary. One suitable calculation method for such problems is the Fourier Modal Method (FMM) [19-21], which can easily incorporate the radiative dissipation of a QNM leaking from the resonator through the waveguide [22].

Fig. 1. Top: Absolute value of the QNM | ˜ f c ( r ) | in a cavity side-coupled to an infinite waveguide in a PC with lattice constant a . The QNM has a complex resonance frequency of ˜ ω a/ π c 2 c = 0 3969 . -0 0014i, corresponding . to Q = 146. Bottom: Real part (red full), imaginary part (blue dashed) and absolute value (black dashed-dotted) of the QNM along the line y = 0 in the center of the waveguide. In both panels, the QNM is scaled to unity at r = r c in the center of the cavity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0609v1-figure-be1c63a0-ec67-438b-88d2-0ab47db48760.jpg)

The leaky nature of the QNMs leads to an exponential divergence at large distances, which means that they cannot be normalized in the same way as the so-called normal modes calculated using Dirichlet (or periodic) boundary conditions. Instead, the QNMs may be normalized by an alternative prescription that explicitly compensates for the divergence. This normalization was first derived for spherically symmetric material systems in Refs. [4-8] and applied to general optical cavities in Ref. [9], which also introduced a generalized effective mode volume for leaky optical cavities. Using the Lorentz reciprocity theorem, Sauvan et al. [11] later derived an alternative formulation of the norm in an elegant and very transparent way. When the Silver-M ˜ A ller ¼ radiation condition applies, the two formulations of the norm can be shown to be identical [14]. As noted above, this is not the case for cavities coupled to PC waveguides. Whereas the norm in Refs. [4-8] rely explicitly on the Silver-M ˜ A ller ¼ condition, the formulation in Ref. [11] does not. Therefore, one could expect this formulation to be general enough to also handle QNMs leaking through waveguides but, as we show below, a direct application of the norm in Ref. [11] leads to a divergent integrand. As a remedy, we show in this Letter how one can use the theory of divergent series [23] to assign a finite value to the integral, and we use the normalization in an example problem of perturbation theory for the PC cavitywaveguide system in Fig. 1.

We consider systems of cavities coupled to waveguides which are defined in general by a periodic relative permittivity distribution glyph[epsilon1] r ( r ) for which glyph[epsilon1] r ( r + R ) = glyph[epsilon1] r ( r ), where R is a lattice vector in the direction of the waveguide. In addition, we limit the analysis to non-magnetic, isotropic and dispersionless materials, and to cavities coupled to waveguides with a single band of defect waveguide modes in the frequency range of interest. This is the case, for example, for the material system in Fig. 1, cf. chapter 5 of Ref. [18]. We focus on electric field QNMs defined as solutions to the wave equation

$$\nabla \times \nabla \times \tilde { \mathbf f } _ { \mu } ( \mathbf r ) - \left ( \frac { \tilde { \omega } _ { \mu } } { c } \right ) ^ { 2 } \epsilon _ { \mathbf r } ( \mathbf r ) \tilde { \mathbf f } _ { \mu } ( \mathbf r ) = 0, \quad ( 1 ) \quad \text{ arou} \\ \text{corr} \quad \text{$\tilde{.}$} \quad \text{$\tilde{.}$}$$

where c is the speed of light, subject to a suitable radiation condition describing the light-propagation through the waveguide away from the cavity. In the periodic waveguides, Bloch-Floquet theory [18] ensures that the solutions to the wave equation may be written as

$$f _ { k } ( r ) = \exp ( i { \mathbf k } \cdot { \mathbf r } ) { \mathbf u } _ { k } ( { \mathbf r } ), \quad \ \ ( 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

in which k is the wave vector and u k ( r + R ) = u k ( r ). At sufficient distances from the cavity, the leaky cavity mode may be written in the form of a single Bloch mode as in Eq. (2) but with a complex wave vector ˜ k µ corresponding to the resonance frequency ˜ ω µ and pointing in the direction of the waveguide away from the cavity. We assume the calculation domain boundary ∂V to be a plane perpendicular to the waveguide direction (a line in two dimensions) with outward pointing unit vector n . In this case, Eq. (2) can be used to define the waveguide radiation condition as a nonlocal boundary condition of the form

$$\tilde { f } _ { \tilde { k } } ( { \mathbf r } ) \Big | _ { { \mathbf r } \in \partial V } = \exp ( i \tilde { \mathbf k } \cdot { \mathbf n } a ) \tilde { f } _ { \tilde { k } } ( { \mathbf r } - { \mathbf n } a ). \quad \ ( 3 )$$

This non-local boundary condition is similar to the one introduced in Ref. [24], but the application is different in that there is no incoming field; the QNMs are defined as the solutions to the wave equation with no sources. Using Eq. (3), one can calculate QNMs of coupled cavity-waveguide systems using standard frequency domain methods. For most problems of practical interest, however, the dispersion of the waveguide modes has no closed form expression, and this complicates the use of Eq. (3) in QNM calculations for which the complex resonance frequency is the eigenvalue of interest. In practice, therefore, we calculate the complex wave vector ˜ k µ using a Taylor expansion approximation of the dispersion along the real frequency axis, which may then be readily extended to the complex frequency plane by analytic continuation.

As with other frequency domain calculations of QNMs, the radiation condition results in a non-linear problem because it depends on the complex resonance frequency. Therefore, one must calculate the QNMs by an iterative procedure in which a fixed frequency ˜ ω guess is used to set up the equation system which is then subsequently solved to find the frequency closest to ˜ ω guess . The iteration continues until the difference ∆˜ is less ω than some prescribed tolerance ∆˜ ω max . All calculations in this Letter were performed using FEM (Comsol Multiphysics 4.3a), and the iterations were performed with an absolute tolerance of ∆˜ ω max a/ c = 10 -5 . As for discretization, the number of elements was kept sufficiently high so as not to influence the results to the quoted number of digits. The analytic continuation of the dispersion curve into the complex frequency plane was based on a fourth order polynomial fit to the real dispersion data around the point ˜ ω R a/ π 2 c = 0 395. Comparing to the . corresponding 5th order expansion reveals a difference ∆˜ ka/ π &lt; 2 10 -6 over an interval | ˜ ω -˜ ω R | a/ π 2 c &lt; 0 01. . This interval is more than an order of magnitude larger than the imaginary part of the calculated resonance frequency, wherefore we expect the error in the complex dispersion to be negligible for the calculations presented here. Finally, an important source of error in the calculations is linked to the size of the calculation domain and can be easily described in the language of the FMM [22]. The FMM approach is based on an expansion using the full set of Bloch modes in the system - with either propagating, growing or decaying characteristics. In the waveguides, the field is described by several decaying Bloch modes as well as a single growing Bloch mode [22], and Eq. (3) thus requires all purely decaying Bloch mode components of the QNM to be negligibly small at the calculation domain boundary. Figure 2 illus-

trates the convergence of the QNM resonance frequency as the calculation domain size L x is increased. Specifically, it shows the error δω ˜ L = ˜( | ω L x ) -˜( ω L max ) | when comparing to a fixed domain size L max = 18 a as well as the error δω ˜ L = ˜( | ω L x ) -˜( ω L x + a ) | when comparing the results obtained using successively larger domain sizes. In both cases, we used a fixed calculation domain height of L y = 18 . Because of the exponential decay a of the QNM along the y -direction, this was found to be large enough to not influence the results to the quoted number of digits. The error in the resonance frequency decreases exponentially with the domain size, consistent with the behavior of the non-propagating Bloch mode components of the QNM, which decay exponentially into the waveguide.

Fig. 2. Change in resonance frequency with increasing calculation domain width L x as defined in the inset. The red solid curve shows the difference to a reference calculation with fixed domain size L max = 18 a , and the blue dashed curve shows the difference between results using successively larger domain sizes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0609v1-figure-2d8c2b5f-5fb9-4385-b873-90acd7851e71.jpg)

The procedure outlined above was applied to calculate the QNM in Fig. 1, which we denote by µ = c and which has only an out-of-plane component, ˜ ( f c r ) = ˜ ( f c r e ) z . The QNM has a complex resonance frequency of ˜ ω a/ π c 2 c = 0 3969 . -0 0014i, . corresponding to Q = 146. From analytic continuation of the dispersion we find that the magnitude of the complex wave vector in the direction of the waveguide is ˜ k a/ c 2 π = 0 2837 . -0 0026i. .

Having calculated the QNM, we next turn to the question of normalization. Following Ref. [11], we start by writing the norm of the QNM ˜ f µ as

$$\langle \langle \tilde { \mathfrak { f } } _ { \mu } | \tilde { \mathfrak { f } } _ { \mu } \rangle \rangle & = \frac { 1 } { 2 } \int _ { V } \epsilon _ { \mathrm r } ( \mathrm r ) \tilde { \mathfrak { f } } _ { \mu } ( \mathrm r ) \cdot \tilde { \mathfrak { f } } _ { \mu } ( \mathrm r ) & & \text{volum} \\ & \quad + \left ( \frac { \tilde { \omega } _ { \mu } } { c } \right ) ^ { 2 } ( \nabla \times \tilde { \mathfrak { f } } _ { \mu } ( \mathrm r ) ) \cdot ( \nabla \times \tilde { \mathfrak { f } } _ { \mu } ( \mathrm r ) ) \mathrm d \mathrm r, \ \ ( 4 ) & & \text{the c$} \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where the integral is formally over the entire volume of space (area in two dimensions), and we have used the source-free Maxwell curl equation to express the integrand in terms of ˜ f µ only. When the QNM leaks into homogeneous surroundings, the integrand in Eq. (4) tends to zero at large distances from the cavity sufficiently fast to render the integral convergent. When the QNM leaks through a periodic waveguide this is not necessarily so; for the QNM in Fig. 1 we find that in the waveguide, the integrand oscillates around zero and increases (exponentially) in magnitude as a function of the distance from the cavity. The QNM in Fig. 1 is clearly symmetric wherefore we limit the below calculations to x &gt; 0, but we note that the procedure generalizes immediately to non-symmetric material systems. Moreover, we split the x -integration into two parts at x = x 0 which is chosen sufficiently large that for x &gt; x 0 the QNM in the waveguide is well described by Eq. (2). For x &lt; x 0 , the integrand is well-behaved, and the integral can easily be evaluated. For the problematic part over the semi-infinite region x &gt; x 0 , we rewrite the integral as

$$I _ { x > x _ { 0 } } = I _ { a } ( x _ { 0 } ) \sum _ { m = 0 } ^ { \infty } e ^ { 2 i \tilde { k } _ { \mu } m a },$$

where

$$\underset { 1 2 } { \longrightarrow } \Big | \quad & I _ { a } ( x _ { 0 } ) = \frac { 1 } { 2 } \int _ { \perp } \int _ { x _ { 0 } } ^ { x _ { 0 } + a } \epsilon _ { r } ( { \mathbf r } ) \tilde { \mathbf f } _ { \mu } ( { \mathbf r } ) \cdot \tilde { \mathbf f } _ { \mu } ( { \mathbf r } ) \\ & \quad + \left ( \frac { \tilde { \omega } _ { \mu } } { c } \right ) ^ { 2 } ( \nabla \times \tilde { \mathbf f } _ { \mu } ( { \mathbf r } ) ) \cdot \left ( \nabla \times \tilde { \mathbf f } _ { \mu } ( { \mathbf r } ) \right ) \mathrm d x \mathrm d { \mathbf r } _ { \perp }, \ ( 6 ) \\ \quad \cdot \quad \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{$$

in which ⊥ denotes the dimension(s) perpendicular to the x -axis (area in three dimensions and line in two dimensions). The sum in Eq. (5) is a geometric series of the form ∑ m b m , which is formally divergent for | b | ≥ 1. The cavity mode in Fig. 1 is an example of such a divergent case, for which we find that exp(2i ˜ k a c ) = -0 9421 . -0 4243i. From Eq. (5), this can . now be seen as the origin of the oscillating and divergent behavior of the integrand. Although the series is formally divergent, it is possible to assign a finite value to the right hand side of Eq. (5) using Borel summation (for Re { } b &lt; 1) or Lindel¨f or Mittag-Leffler summation o (for b ∈ C \ [1 , ∞ [) [23] , in which cases the sum evaluates to the same expression as in the case of | b | &lt; 1, and one thus finds

$$I _ { x > x _ { 0 } } = \frac { I _ { a } } { 1 - e ^ { 2 i \tilde { k } _ { \mu } a } }. \quad \quad ( 7 ) \\ \intertext { m a t i n d } \max \min \min \epsilon \, T \, = \, 1 \circ \max \max \max \, \text{background}$$

Using a calculation domain size of L x = 18 a we choose x 0 = 8 a and scale the QNM to unity in the cavity center r c , in which case the generalized effective mode volume [9] evaluates to v Q ≡〈〈 ˜ f c | ˜ f c 〉〉 / ˜ f 2 c ( r c ) a 2 = 1 441 . -0 055i (the size of the contribution from . | x &gt; x | 0 in this case is 2 | I x&gt;x 0 | / ˜ f 2 c ( r c ) a 2 = 0 005). We note, that when . the cavity mode leaks into homogeneous surroundings as well as into a guided waveguide mode, as in PC membranes in three dimensions, for example, the integrand in Eq. (4) will tend to zero in regions of homogeneous

material far from the cavity, whereas it will oscillate and increase in magnitude for positions inside the waveguide. Therefore, we expect the above prescription to generalize immediately to such systems as well.

Last, as an example application, we use the norm in Eqs. (4) and (5) to study the influence of small material changes. Using perturbation theory, the first order change in the resonance frequency due to a small change in the permittivity distribution is given as

$$\Delta \tilde { \omega } _ { c } ^ { ( 1 ) } = - \frac { \tilde { \omega } _ { c } } { 2 \langle \langle \tilde { f } _ { c } | \tilde { f } _ { c } \rangle \rangle } \int _ { V } \Delta \epsilon _ { r } ( r ) \tilde { f } _ { c } ( r ) \cdot \tilde { f } _ { c } ( r ) d r, \quad ( 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \sim \dots \quad \Lambda \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \end{cases},$$

where ∆ ( ) is the position dependent change in permitglyph[epsilon1] r tivity. Eq. (8) is identical to the result for QNMs leaking to homogeneous media [4,8]; it may be derived, for example, using the Hellmann-Feynman theorem [25] in an operator formulation of Maxwell's equations similar to Ref. [26] but with a non-Hermitian approach using the QNM norm in Eq. (4). We consider perturbations to the material system in which an additional low-permittivity rod is inserted in the center of the cavity. As shown in Fig. 3, the perturbation result of Eq. (8) agrees well with full numerical FEM calculations even for relatively large shifts of several linewidths.

Fig. 3. Perturbative change in the resonance frequency ∆˜ ω c = ∆ ω c -i∆ γ c as a function of permittivity change ∆ glyph[epsilon1] r of a cylinder in the center of the cavity as indicated in the inset. Red and blue circles (left axis) show changes in resonance frequency and decay rate, respectively, using full numerical calculations, and dashed black lines show the first order perturbation result of Eq. (8). Gray shading (right axis) shows the error in the perturbation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0609v1-figure-e8605481-b131-44be-958a-f458672838c2.jpg)

In summary, we have described an approach for numerical calculation of QNMs in cavities coupled to waveguides. The approach relies on a non-local boundary condition to correctly model the radiation condition for light leaking through the waveguide. In addition, we have shown how to normalize these QNMs by use of the theory of divergent series to assign a finite value to the divergent integral appearing in the norm. Finally, we have used the normalization procedure for perturbation theory calculations to illustrate how it leads to the correct first order predictions of the resonance shift.

This work was supported by the Carlsberg Foundation and the Danish Council for Independent Research (FTP 10-093651)

## References

- 1. E. S. C. Ching, P. T. Leung, A. Maassen van den Brink, W. M. Suen, S. S. Tong, and K. Young. Rev. Mod. Phys. , 70, 1545 (1998).
- 2. E. S. C. Ching, P. T. Leung, and K. Young. Optical processes in microcavities-the role of quasi-normal modes, R. K. Chang and A. J. Campillo Eds. World Scientic (1996).
- 3. P. T. Kristensen and S. Hughes. ACS Photonics , 1, 2 (2014).
- 4. H. M. Lai, P. T. Leung, K. Young, P. W. Barber, and S. C. Hill. Phys. Rev. A , 41, 5187 (1990).
- 5. P. T. Leung, S. Y. Liu, and K. Young. Phys. Rev. A , 49, 3057 (1994).
- 6. P. T. Leung and K. M. Pang. J. Opt. Soc. Am. B , 13, 805 (1996).
- 7. K. M. Lee, P. T. Leung, and K. M. Pang. J. Opt. Soc. Am. B , 16, 1409 (1999).
- 8. K. M. Lee, P. T. Leung, and K. M. Pang. J. Opt. Soc. Am. B , 16, 1418 (1999).
- 9. P. T. Kristensen, C. Van Vlack, and S. Hughes. Opt. Lett. , 37, 1649 (2012).
- 10. B. Maes, J. Petr´ˇek, S. Burger, P. Kwiecien, J. Luksch ac and I. Richter. Opt. Ex. , 21, 6794 (2013).
- 11. C. Sauvan, J. P. Hugonin, I. S. Maksymov, and P. Lalanne. Phys. Rev. Lett. , 110, 237401 (2013).
- 12. J. R. de Lasson, J. Mørk, and P. T. Kristensen. J. Opt. Soc. Am. B , 30, 1996 (2013).
- 13. Q. Bai, M. Perrin, C. Sauvan, J-P Hugonin, and P. Lalanne. Opt. Ex. , 21, 27371 (2013).
- 14. R.-C. Ge, P. T. Kristensen, J. F. Young, and S. Hughes. arXiv:1312.2939 [quant-ph] (2013).
- 15. R.-C. Ge and S. Hughes. Opt. Lett. , 39, 4235 (2014).
- 16. P. Lalanne, C. Sauvan, and J.P. Hugonin. Laser &amp; Phot. Rev. , 2, 514 (2008).
- 17. J.-P. Berenger. J. Comp. Phys. , 114, 185 (1994).
- 18. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade. Photonic Crystals - Molding the Flow of Light, second edition . Princeton University Press (2008).
- 19. M. G. Moharam and T. K. Gaylord. J. Opt. Soc. Am. , 71, 811 (1981).
- 20. E. Noponen and J. Turunen. J. Opt. Soc. Am. A , 11, 2494 (1994).
- 21. G. Lecamp, J. P. Hugonin, and P. Lalanne. Opt. Ex. , 15, 11042 (2007).
- 22. J. R. de Lasson, P. T. Kristensen, J. Mørk, and Niels Gregersen. arXiv:1405.2595 [physics.optics] (2014).
- 23. G. H. Hardy. Divergent Series . Oxford at the Clarendon Press (1949).
- 24. Z. Hu and Y. Y. Lu. J. Opt. Soc. Am. B , 29, 1356 (2012).
- 25. R. E. Allen and M. Menon. Phys. Rev. B , 33, 5611 (1986).
- 26. B. A. Daniel, D. N. Maywar, and G. P. Agrawal. J. Opt. Soc. Am. B , 28, 2207 (2011).
Philip Trøst Kristensen Jakob Rosenkrantz de Lasson Niels Gregersen
Calculation, normalization and perturbation of quasinormal modes in coupled cavity-waveguide systems
We show how one can use a non-local boundary condition, which is compatible
with standard frequency domain methods, for numerical calculation of
quasinormal modes in optical cavities coupled to waveguides. In addition, we
extend the definition of the quasinormal mode norm by use of the theory of
divergent series to provide a framework for modeling of optical phenomena in
such coupled cavity-waveguide systems. As an example, we apply the framework to
study perturbative changes in the resonance frequency and Q value of a photonic
crystal cavity coupled to a defect waveguide.