Abstract
Faraday rotation is arguably the most widely studied nonreciprocal phenomenon, observable in the form of polarization rotation as waves travel in magneto-optical materials. It is at the basis of the realization of various forms of isolators and circulators. However, magneto-optical materials are bulky and difficult to integrate, and inherently associated with losses. Here we leverage nonlinear phenomena in conventional optical crystals to mimic Faraday rotation in a resonant cavity excited with two pumps with opposite circular polarizations and slightly detuned frequencies. The effect can be observed in basic resonant structures with no need for material resonances or tailored dispersion, arbitrary frequency of the pumps, and low loss, yielding an interesting path towards the realization of compact magnet-free nonreciprocal optical elements.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. INTRODUCTION
The realization of magnetless nonreciprocal components has been a major theme in recent photonics research. Nonreciprocity is useful for source protection and necessary to improve optical chip packing density. Magnetic materials, while popular in the microwave regime, are known to be lossy and difficult to integrate with current photonic technologies, though there has been recent progress with silicon-compatible materials, such as graphene and Ce:YIG [1,2]. Spatiotemporal modulation has been explored as an effective way to create nonreciprocal microwave components, with interesting pathways towards the implementation of magnet-free nonreciprocity [3–6]. In optics, beyond electro-optical modulation which is fundamentally limited in terms of speed, this effect can be realized with a traveling pump wave in nonlinear media [7]. For example, a Kerr or optomechanical ring resonator can be used to break degeneracy between counter-propagating modes, therefore breaking time-reversal symmetry [8–10]. Nonlinear resonators [11–15] can provide isolation without external modulation, but have limited dynamic range, and are fundamentally limited by time-reversal symmetry and dynamic reciprocity [16,17]. Given the limitations of current technologies, research in improved nonreciprocal responses that do not rely on magneto-optics is essential for progress in nanophotonic and quantum technology.
The Faraday effect is a widely known non-reciprocal phenomenon, wherein a magnetized sample rotates the polarization angle of a propagating field. The rotation handedness depends on the applied bias and not on the direction of propagation, hence leading to a nonreciprocal response. Due to microscopic reversibility, a circularly polarized optical pump can also induce magnetization in magneto-optical materials, a phenomenon known as inverse Faraday effect [18,19]. Such a response can in principle be used to rotate the polarization of a probe signal with a single optical pump, but it still requires magnetic media, which are generally not ideal for integration and practical implementation. A circularly polarized pump can also induce Faraday rotation in non-magnetic media, but it requires either a coherent interaction between a pump at the same frequency as the probe [20] or breaking Kleinmann symmetry, which is only possible in lossy, resonant media, posing an important tradeoff between the amount of rotation and the transmitted power through the media [20–23]. Recently, this concept has been explicitly explored by coupling to phonons [24,25]. While the effect may be strong, it requires careful tuning of pump and probe frequencies and the mechanical resonance, and it is narrowband in nature, being limited by the mechanical resonance frequency.
Here, we introduce an alternative approach to magnetless Faraday rotation in media with instantaneous third-order nonlinearity, by biasing the system with two optical pumps with opposite circular polarizations and detuned frequencies. Our approach does not require material loss, nor phase-matching between pump and probe signals, making it insensitive to the actual frequency of the pump signal and material dispersion, and widely applicable to waves propagating in opposite directions. We discuss a possible implementation in a compact Fabry–Perot resonator supporting degenerate modes, but alternative configurations are also available, for instance, in microresonators supporting degenerate modes and coupled to a waveguide. In addition to its theoretical significance as an alternative way to induce Faraday rotation via optical pumping, the proposed approach may be used in the design of magnetless isolators and other nonreciprocal devices for nanophotonic and quantum applications.
2. THEORY
We start our analysis with a brief discussion of the different mechanisms that can lead to Faraday rotation via optical pumping. First, consider the fundamental conservation laws in a scattering event: energy, angular momentum, and linear momentum. Linear momentum is dependent on the dispersive material properties, but photon frequency and spin are well-defined and sufficient to capture the behavior we describe below. Polarization is an important consideration for nonlinear processes and results in selection rules for a given lattice geometry [23]. We further assume that no energy or momentum is transferred to the lattice, as in a lossless isotropic medium. We decompose the linear signal into right-handed (RH, ↻) and left-handed (LH, ↺) components, defined with respect to a fixed axis instead of the wave propagation direction, and assume a RH pump at [20,26]. Consider the first transition in Fig. 1(a): the relevant component is , where the L component arises because conjugating a circularly polarized wave reverses its handedness. The second transition in Fig. 1(b) is described by . Away from resonances, when Kleinmann symmetry holds, we can perform the symmetry permutations

Fig. 1. (a), (b) Allowed four-wave interactions for a linear signal () with a right-handed circularly polarized (↻) pump () in an isotropic medium. (c), (d) Allowed interactions with oppositely polarized pumps with frequency spacing . (e), (f) Model of the resonator model in (e) transmission mode and (f) reflection mode.
Equation (1) is the general permutation relation for any nonlinear process, while Eq. (2) follows from Kleinmann symmetry [21]. The material is isotropic, , which follows after performing a 180° rotation, hence the transitions in Figs. 1(a) and 1(b) are identical. This is unsurprising when one thinks in terms of - and - polarization components, due to symmetry. It is known that the refractive-index increase from a strong pump in the presence of instantaneous third-order nonlinearity is independent of the relative phase of pump and probe. Thus, the phase difference between the - and - components of the pump should not matter, and there can be no first-order all-optical Faraday effect in such a medium.
Two other possible nonlinear processes are the conversion of a LH probe into a RH photon at , and the conversion of a RH probe to a LH photon at . These new photons will oscillate at their own frequencies and can then be converted back to photons at , in principle interfering with the signal at . In order for this effect to be significant, the up/down-converted signals need to be phase-matched with the signal at or, in the case of a resonator, they need to couple to one of the resonances. Such an approach requires carefully calibrated resonant frequency spacings and widths, which may be difficult to realize in practice. One idea to overcome this problem may be to use a small , so that the converted signals fall within the same resonance as the original signal, but generally there are not powerful sources at low enough frequency to induce such an effect.
We propose instead a different approach to induce strong Faraday rotation based on nonlinear processes and optical pumping based on two detuned pumps. The idea is to induce a nonlinear polarization current in the resonator that rotates in time, analogous to non-reciprocity induced by synthetic angular momentum bias using spatio-temporal modulation [27]. Rather than creating the artificial rotation with a low frequency electro-optical bias needing its own integration into the structure as proposed in [25], we show that beating between two optical waves can create a similar angular-momentum bias in the nonlinear medium, inducing strong polarization conversion and nonreciprocity based on all-optical modulation. Consider the interaction with two circularly polarized pumps with arbitrary frequency, opposite handedness, and with a small frequency detuning . The second pump induces the additional parametric interactions illustrated in Figs. 1(c) and 1(d), and with them photons at . When is comparable with the resonance linewidth, the converted signals couple to the same resonance mode as the original signal. As shown in the following, this approach offers a number of benefits: small footprints and the absence of phase-matching requirements compared to traveling wave structures, as well as complete freedom in the choice of the absolute frequencies of the pumps, since only their difference matters.
In order to analyze the scenario described above, we use coupled-mode theory. The two pump waves are defined as
where are the circular unit vectors, and are the spatial beam profiles. We define handedness here and in the following solely based on the relative phase of the x- and y-components of the fields, and not their direction of propagation. The physical fields are given by [21]. The fields in the th cavity mode are with normalization , so that gives the energy in mode with the given polarization. The RH and LH components are initially degenerate because of isotropy and symmetry.Using the standard assumption of slowly varying , the temporal evolution of the system is governed by
In order to evaluate the nonlinear polarization components, we assume an isotropic medium with where [21]. Utilizing , , and , and keeping only the terms near and first-order in , we find with and . The top line in Eq. (7) corresponds to a resonance shift due to the intensity-dependent refractive index, and is not interesting for our purpose. We will neglect it, pointing out that in general will shift slightly as the pump power increases. The second and third lines, on the other hand, describe the coupling of the left-handed cavity mode to a right-handed component and vice versa, proportional to the product of the two pump waves.After taking the second derivative of (7) (again neglecting ), plugging into Eq. (5), we find the closed-cavity coupled-mode equations
where is the coupling constant, and the resonance frequency is denoted as instead of for clarity, since we will use a single resonance for the rest of the paper. Note that is the only factor in the integrand that is not guaranteed to be positive and real, and its spatial behavior is important to determine the total coupling strength. If two pump waves are co-propagating, they should have a small wave-vector mismatch and thus should have a slowly varying phase over the interaction volume, ensuring that the contributions to the integral are phase-matched and contribute constructively. On the other hand, if the pumps travel in opposite directions, the situation may get more complicated, and special care should be paid to make sure that is large.Equation (8) manifests the photon interactions of Figs. 1(c) and 1(d), and illuminates the underlying phenomena behind our approach, which we will now study in the time-domain. Since the system is linear in the small-signal approximation, we can decompose the cavity components as , where are the primary cavity fields, excited from the incident fields and therefore sharing their handedness, while are the secondary fields excited only through coupling terms in Eq. (8) and thus of opposite handedness. The secondary fields are therefore zero before the moment at which the optical pump is turned on ( in Fig. 2). For RH and LH excitations, the initial conditions are and , respectively. This approach highlights the effect of the pumps on right- and left-handed signals (initial fields) independently. For simplicity, we assume , and so we can neglect the first term in the bracket on the right-hand side of Eq. (8). Hence, the impulse response of the envelopes for RH and LH excitations is given by
where . We plot the results for the envelope functions in a normalized time domain in Fig. 2 for two regimes.
Fig. 2. (a), (b) Time-domain envelope amplitudes for , ; (c), (d) time-domain phases for the fundamentals , for cavities with , , and (a), (c) or (b), (d) . The time is normalized to the low-detuning envelope oscillation period .
Figures 2(a) and 2(b) show the amplitude of the envelope functions when the initial field is entirely right-handed. It can be seen that, as the system evolves in time, there is a periodic exchange of power between the primary RH and the secondary LH signals with frequency . When is small compared to , there is a large transfer of power into the sideband , whereas this transfer is smaller as increases. This power exchange between the primary and secondary signals is accompanied by rapid changes in the phase of the primary signals, as seen in Figs. 2(c) and 2(d), showing the arguments of for the same as in Figs. 2(a) and 2(b). Importantly, the phase accumulation is opposite for RH and LH signals, leading to a nonreciprocal response, driven by the pumps breaking time-reversal symmetry. This property can be better understood by realizing that, for RH excitation, the RH photons in the cavity are coupled by the pump to lower frequency LH photons, and vice versa. The transition rates are essentially identical for the two processes. However, during their duration in the sidebands, the up-converted photons acquire more phase than the photons at the original frequency, while the down-converted photons acquire less phase. Both sidebands then couple back to the original frequency, bringing the accumulated phase with them and resulting in opposite phases for RH and LH signals, as in Figs. 2(c) and 2(d). The difference in phase, and therefore the strength of the nonreciprocal response, is both a function of the frequency difference and the proportion of photons that participate in the up- and down- conversion processes.
We can then use this intuition to explain the response of the system under a linearly polarized excitation. A linear input will have equal RH and LH components. When the cavity is driven close to , the RH and LH resonances will be excited roughly equally, and the field at the fundamental frequency will rotate in time with an angle given by . After the fields spend a time in the cavity equal to the modal lifetime determined by the -factor, the mode will leak out. We can phenomenologically add the coupling into and out of the cavity using coupled-mode theory [28], assuming negligible material loss. For a two-sided cavity, with modal input power from one side of the cavity, and leakage rate , we have
3. RESULTS
For a linear input at ( is the detuning of the probe frequency from the resonance frequency of the resonator, accounting for the Kerr frequency shift), we can solve for steady-state Fourier amplitude that couples to , and then separately for and . For a system with no direct pathway and no material loss, the transmission coefficient , so
and reflection . Figures 3(a) and 3(b) show versus the input frequency and the pump power. It can be seen that the presence of the pumps breaks the degeneracy of the RH and LH modes in the cavity, which now resonate at different frequencies. Once sufficient splitting of the RH and LH transmission resonances are achieved, quarter-wave plates and polarization-dependent beam splitters can be added to build a 4-port circulator, as described in [24]. Similarly, a two-port isolator can be built with material loss or absorption of the fields on the second side of the cavity. Since these setups have been previously proposed, we simply note potential use here and focus the rest of the paper on polarization effects.
Fig. 3. Transmission magnitudes for a symmetric cavity with and silicon-like nonlinearity (see main text for details). is the power in one of the (equally strong) two-pump beams, in . This definition is used in all subsequent plots. The pump beams are detuned as . The pumps break the degeneracy between LH and RH resonances, which can create polarization rotation ().
Exactly at the original resonance (), both RH and LH components have identical transmission and reflection coefficient magnitudes [Figs. 3(a) and 3(b)] but in general different phases. A linearly polarized input wave at this center frequency will thus remain linearly polarized, but with a different polarization angle. This rotation is calculated in Figs. 3(c) and 3(d), which show the transmission coefficients for an incident -polarized wave, highlighting the emergence of a cross-polarized component of the transmitted field. The corresponding polarization rotation is independent of the direction of propagation, hence nonreciprocal. In the figures, in order to define , for the cavity modes we have used the standing wave fields in a Fabry–Perot resonator (the longitudinal mode order does not matter) with refractive index and , close to values for silicon, and we have assumed the pumps to be nonresonant, collinear traveling waves. For such a system, , with being the optical intensity of the pump beams. More fundamentally, these powers represent the intensities of the pump circulating in the cavity. One benefit of the described approach, also stressed above, is that the absolute frequency of the pumps does not matter. Hence, it should be quite feasible to enable a different, wider-linewidth resonance for the pumps to enable strong field enhancement in the cavity for the pump fields, while maintaining a detuning larger than the resonance linewidth of the probe.
One drawback of this design is the fundamental limit in polarization conversion efficiency stemming from symmetries. An ideal Faraday rotator has scattering coefficients and . On the other hand, in the structure described above, for an -polarized incident field, light will leak equally to both sides of the above symmetric cavity, so that . To achieve better rotation efficiency, the longitudinal symmetry of the cavity must be broken. If we consider a one-sided cavity, in which transmission is eliminated, as in Fig. 1(f), the coupled-mode equations are only minorly modified, and the new reflection coefficients are
In this case, there is no limitation on the rotation efficiency for the reflected fields. In Fig. 4, we show the behavior of the reflected wave for and , as a function of pump power and frequency. We define the polarization angle as the major axis of the polarization ellipse, and we utilize the 1st and 4th quadrants [29]. The ellipticity angle magnitude is defined as the inverse tangent of the ratio of minor and major axes, with the sign determined by the RH(+) or LH(−) nature. For smaller detunings, the required pump power levels for a given rotation are smaller, but the reflection at the fundamental frequency is also small, since most of the power is scattered not into y-polarized light, but converted to other frequencies. With increased detuning, more power is needed for any rotation angle, but the reflection efficiency at the signal frequency improves. At the trivial singularity , all the modes under consideration are at the same frequency and the pump is merely a single linearly polarized signal. In this case we correctly recover that the nonreciprocal polarization conversion goes to zero, yet there may be a reciprocal component due to nonlinearity-induced birefringence [21], which we have ignored from Eq. (6) in our previous analysis. At extremely large pump powers, the rotation saturates at 0°, independently of the detuning.
Fig. 4. (a), (d) Rotation angle, (b), (e) ellipticity, and (c), (f) fraction of power reflected at the input frequency for (a)–(c) and (d)–(f) with . The black lines indicate the contours for , around the typically interesting value of 45°.
To better map the efficiency as a function of pump detuning, we show in Fig. 5 the scattering coefficients for and , with equal intensities in each pump beam, versus detuning and pump power. For a given pump detuning, the rotation angle monotonically increases with increasing pump power, while the reflected wave remains linearly polarized. This allows efficient, nonreciprocal polarization shift at any angle with sufficient detuning and pump power. We stress that in all these examples we have assumed a non-resonant pump; utilizing resonance enhancement for the pump can increase the coupling coefficient with the same input power.

Fig. 5. (a) Polarization rotation angle and (b) reflection at the fundamental harmonic for probes at the resonance center as a function of pump detuning and power, with and silicon-like nonlinearity. Non-unitary reflection coefficients are due to conversion to frequency sidebands.
To confirm the validity of our approach, we implemented a one-and-a-half-dimension (one spatial dimension and two polarizations) finite-difference time-domain (FDTD) solver utilizing auxiliary differential equations [30]. A dielectric resonator was chosen to mimic lossless silicon, and the resonator was backed by a perfect-electric conductor (PEC), while the front mirror was realized as a thin slab with Drude material dispersion, so that it is reflective at the probe frequency but transparent to the pump waves. While this particular example was chosen for numerical simulation simplicity, similar effects may be reasonably achieved experimentally with frequency-selective surfaces or Bragg mirrors. The probe frequency was chosen to be at the center of a Fabry–Perot resonance peak, with -factor . Figure 6 shows numerical calculations of this structure, which match well with the coupled-mode analysis. The probe frequency was chosen to track the resonance center with increasing pump power due to the Kerr effect, a frequency shift of order of a linewidth for the highest pump powers used. The spectrum (not shown) confirms that sidebands exist at only , as required by conservation laws, minimizing insertion loss due to frequency conversion. Overall, this approach offers tunability over the entire range of rotation angles while maintaining low insertion loss and ellipticity, important properties of Faraday rotators.

Fig. 6. FDTD results for (a) polarization and ellipticity angle and (b) fundamental amplitude as a function of pump power. The cavity is filled with a silicon-like nonlinear material (, ), considering a detuning , and a coupling constant , due to the reflection of the pump beams. is in .
4. CONCLUSION
In this paper, we have explored the use of instantaneous nonlinear effects to mimic Faraday rotation in resonators, with a number of benefits over other approaches based on inverse Faraday phenomena or nonresonant pumping. In steady-state, these structures are linear and non-reciprocal to small signals, unlike static nonlinear isolators. The beating of two pumps in the nonlinear resonators establishes an effective synthetic angular momentum bias in the structure, which lifts degeneracy and replaces a magnetic bias to break reciprocity. The interactions we describe can be used to preferentially transmit RH or LH light and create a circulator, or as a polarization rotator. The proposed platform is easy to adapt to a wide range of settings, as the cavity structure can be simple and there are no requirements for the dispersion of the materials or the relative frequencies of pump and probe beams. This may be particularly important for utilization with silicon, which has a strong nonlinearity but suffers from two-photon absorption for pumps at telecom wavelengths. Furthermore, unlike nonreciprocal devices based on traveling wave modulation in which the nonreciprocal effect essentially exists in only one direction due to phase-matching, in our system the rotation may be applied in both directions in the transmitting resonator. Because it does not rely on directional-dependent phase-matching, the proposed structure can fit in small footprints that lend themselves to modern nanophotonic systems. We believe that these findings open important pathways for integrated nonreciprocal devices in classical and quantum computing and processing systems.
Funding
Air Force Office of Scientific Research; National Science Foundation; Simons Foundation.
REFERENCES
1. D. Huang, P. Pintus, C. Zhang, P. Morton, Y. Shoji, T. Mizumoto, and J. E. Bowers, “Dynamically reconfigurable integrated optical circulators,” Optica 4, 23–30 (2017). [CrossRef]
2. T. Mizumoto, R. Baets, and J. E. Bowers, “Optical nonreciprocal devices for silicon photonics using wafer-bonded magneto-optical garnet materials,” MRS Bull. 43(6), 419–424 (2018). [CrossRef]
3. N. A. Estep, D. L. Sounas, J. Soric, and A. Alù, “Magnetic-free non-reciprocity and isolation based on parametrically modulated coupled-resonator loops,” Nat. Phys. 10, 923–927 (2014). [CrossRef]
4. T. Dinc, M. Tymchenko, A. Nagulu, D. Sounas, A. Alu, and H. Krishnaswamy, “Synchronized conductivity modulation to realize broadband lossless magnetic-free non-reciprocity,” Nat. Commun. 8, 795 (2017). [CrossRef]
5. A. Kord, D. L. Sounas, and A. Alu, “Magnet-less circulators based on spatiotemporal modulation of bandstop filters in a delta topology,” IEEE Trans. Microw. Theory Tech. 66, 911–926 (2018). [CrossRef]
6. D. L. Sounas and A. Alù, “Non-reciprocal photonics based on time modulation,” Nat. Photonics 11, 774–783 (2017). [CrossRef]
7. Z. Yu and S. Fan, “Complete optical isolation created by indirect interband photonic transitions,” Nat. Photonics 3, 91–94 (2009). [CrossRef]
8. L. Del Bino, J. M. Silver, M. T. M. Woodley, S. L. Stebbings, X. Zhao, and P. Del’Haye, “Microresonator isolators and circulators based on the intrinsic nonreciprocity of the Kerr effect,” Optica 5, 279–282 (2018). [CrossRef]
9. D. B. Sohn, S. Kim, and G. Bahl, “Time-reversal symmetry breaking with acoustic pumping of nanophotonic circuits,” Nat. Photonics 12, 91–97 (2018). [CrossRef]
10. F. Ruesink, J. P. Mathew, M.-A. Miri, A. Alù, and E. Verhagen, “Optical circulation in a multimode optomechanical resonator,” Nat. Commun. 9, 1798 (2018). [CrossRef]
11. S. Lepri and G. Casati, “Asymmetric wave propagation in nonlinear systems,” Phys. Rev. Lett. 106, 164101 (2011). [CrossRef]
12. L. Fan, J. Wang, L. T. Varghese, H. Shen, B. Niu, Y. Xuan, A. M. Weiner, and M. Qi, “An all-silicon passive optical diode,” Science 335, 447–450 (2012). [CrossRef]
13. I. V. Shadrivov, V. A. Fedotov, D. A. Powell, Y. S. Kivshar, and N. I. Zheludev, “Electromagnetic wave analogue of an electronic diode,” New J. Phys. 13, 33025 (2011). [CrossRef]
14. H. Zhou, K.-F. Zhou, W. Hu, Q. Guo, S. Lan, X.-S. Lin, and A. V. Gopal, “All-optical diodes based on photonic crystal molecules consisting of nonlinear defect pairs,” J. Appl. Phys. 99, 123111 (2006). [CrossRef]
15. D. L. Sounas, J. Soric, and A. Alù, “Broadband passive isolators based on coupled nonlinear resonances,” Nat. Electron. 1, 113–119 (2018). [CrossRef]
16. D. L. Sounas and A. Alù, “Fundamental bounds on the operation of Fano nonlinear isolators,” Phys. Rev. B 97, 115431 (2018). [CrossRef]
17. Y. Shi, Z. Yu, and S. Fan, “Limitations of nonlinear optical isolators due to dynamic reciprocity,” Nat. Photonics 9, 388–392 (2015). [CrossRef]
18. J. P. van der Ziel, P. S. Pershan, and L. D. Malmstrom, “Optically-induced magnetization resulting from the inverse Faraday effect,” Phys. Rev. Lett. 15, 190–193 (1965). [CrossRef]
19. P. S. Pershan, J. P. van der Ziel, and L. D. Malmstrom, “Theoretical discussion of the inverse Faraday effect, Raman scattering, and related phenomena,” Phys. Rev. 143, 574–583 (1966). [CrossRef]
20. Y. P. Svirko and N. I. Zheludev, “Coherent and incoherent pump-probe specular inverse Faraday effect in media with instantaneous nonlinearity,” J. Opt. Soc. Am. B 11, 1388–1393 (1994). [CrossRef]
21. R. W. Boyd, Nonlinear Optics (Academic, 2008).
22. S. V. Popov, Y. P. Svirko, and N. I. Zheludev, “Coherent and incoherent specular inverse Faraday effect: χ(3) measurements in opaque materials,” Opt. Lett. 19, 13–15 (1994). [CrossRef]
23. C. L. Tang and H. Rabin, “Selection rules for circularly polarized waves in nonlinear optics,” Phys. Rev. B 3, 4025–4034 (1971). [CrossRef]
24. R. Duggan, J. del Pino, E. Verhagen, and A. Alù, “Optomechanically induced birefringence and optomechanically induced Faraday effect,” Phys. Rev. Lett. 123, 023602 (2019). [CrossRef]
25. M. Lawrence and J. A. Dionne, “Nanoscale nonreciprocity via photon-spin-polarized stimulated Raman scattering,” Nat. Commun. 10, 3297 (2019). [CrossRef]
26. M. S. Wismer, M. I. Stockman, and V. S. Yakovlev, “Ultrafast optical Faraday effect in transparent solids,” Phys. Rev. B 96, 224301 (2017). [CrossRef]
27. D. L. Sounas, C. Caloz, and A. Alù, “Giant non-reciprocity at the subwavelength scale using angular momentum-biased metamaterials,” Nat. Commun. 4, 2407 (2013). [CrossRef]
28. S. Fan, W. Suh, and J. D. Joannopoulos, “Temporal coupled-mode theory for the Fano resonance in optical resonators,” J. Opt. Soc. Am. A 20, 569–572 (2003). [CrossRef]
29. E. Collett, Field Guide to Polarization (SPIE, 2005).
30. A. Taflove and S. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, 2005).