Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

All-optical light storage in bound states in the continuum and release by demand

Open Access Open Access

Abstract

In the framework of the temporal coupled mode theory we consider bound states embedded in the continuum (BSC) of photonic crystal waveguide as a capacity for light storage. A symmetry protected BSC occurs in two off-channel microresonators positioned symmetrically relative to the waveguide. We demonstrate that the symmetry protected BSC captures a fraction of a light pulse due to the Kerr effect as the pulse passes by the microresonators. However the amount of captured light is found to be strongly sensitive to the parameters of the gaussian light pulse such as basic frequency, duration and intensity. In contrast to the above case the BSC resulted from a full destructive interference of two eigenmodes of a single microresonator accumulates a fixed amount of light dependent on the material parameters of the microresonator but independent of the light pulse. The BSCs in the Fabry-Perot resonator show similar effects. We also show that the accumulated light can be released by a secondary pulse. These phenomena pave a way for all-optical storage and release of light.

© 2015 Optical Society of America

1. Introduction

One of the most technologically attractive approaches for light storage is the use of slow wave designs of optical periodic structures to slow down the light pulse [1–3 ]. Various methods and designs such as direct-coupled resonators [2, 4, 5] and photonic crystal (PhC) waveguides [6] have been proposed to slow down the group velocity of light. Slow light waveguides composed of left-handed materials have been also proposed [7, 8]. In these waveguides, the time average power in the left-handed material layer flows in an opposite direction of that in the normal dielectric layer. For the case when light flows are almost equal, slow light and even trapping can be achieved [9]. The frequency of the slow light regime varies with the thickness of the layer resulting in so-called trapped rainbow [10–12 ].

In this paper we propose to employ bound states in the continuum (BSC) as a novel capacity for light storage. Recently the BSCs have attracted much interest in photonics owing to a possibility to trap light in various photonic systems [13–30 ]. The BSC is a localized solution of the Maxwell equations with a discrete eigenfrequency coexisting with the extended modes of continuous spectrum of the PhC waveguide or radiation continuum. Our main goal is to demonstrate that the BSCs are able to capture light due the Kerr effect in the microresonators. Moreover we demonstrate that the application of a secondary light pulset releases trapped light making the BSC potentially interesting for all-optical light storage and release.

The physics of the light storage by the BSCs is the following. Assume for a moment that the amplitude of the injected wave is so small that we can neglect the nonlinearity. Because the BSC is completely decoupled from the continuum it can not be probed by the incoming wave. However with the increase of the injected power the Kerr effect of the microresonators becomes important killing two birds with one stone. First, there is no necessity for tuning material parameters because the Kerr effect in the microresonators results in occurrence of self-induced BSCs [31–33 ]. Second, the nonlinearity couples the BSC with the continuum so that the injected wave excites the BSC transforming it into a quasi-BSC. Once the pulse has passed by the microeresonator the quasi-BSC is again decoupled from the waveguide and becomes a true BSC. As a result some amount of light is trapped in the true BSC opening an opportunity for light storage. Finally, the application of a secondary pulse again transforms the true BSC into a quasi-BSC with finite life-time and releases the light with substantial or even 100% efficiency. Below we inspect these ideas on the basis of three simplest 2D PhC designs which realize the self-induced BSCs.

2. Symmetry protected BSC

We begin with the most obvious system in which light can be trapped by the symmetry protected BSC [16, 20, 23, 34]. For clear visualization of such a BSC we present in Fig. 1 2D photonic crystal (PhC) layout with electric field profile of the antisymmetric mode which is the BSC [16]. However interested reader can easily design different layouts based for instance on nano resonators with extremely high Q-factors in air-hole 2D PhC system [35–37 ]. In order to present the mechanism of light storage in the BSC in the most comprehensive form we assume that among all eigenmodes of the microresonator only the lowest monopole TM mode has the eigenfrequency ω 0 embedded into the first TM propagation band of the directional PhC waveguide. Then the eigenmodes of the total system of two resonators can be classified as symmetric and antisymmetric modes

As,a=(A1±A2)/2,ωs,a=ω0±u,
where variable Aj, j = 1,2 presents the amplitude of the monopole mode in the j-th microresonator and u is due to overlapping of these modes. For simplicity we disregard the dispersion properties of the waveguide and consider that light with the amplitude Ein is injected into the waveguide as shown in equivalent coupled mode theory layout in Fig. 1(b). We assume the injected light is symmetrical relative to the cental line of the waveguide. Then the temporal coupled mode theory (CMT) equations have the following form [38, 39]
iA˙s=(ωs+2iγ)AsiγEineiωt,iA˙a=ωaAa,
where the term γ is responsible for the coupling of the off-channel resonators with the waveguide. Thus, the temporal CMT equations clearly demonstrate that the antisymmetric mode Aa is a symmetry protected BSC.

 figure: Fig. 1

Fig. 1 (a). The single row of the rods is removed from the PhC to form a directional photonic waveguide which supports a band of guided even TM mode. Two defect rods shown in gray with the Kerr effect are inserted symmetrically relative to the waveguide to form identical off-channel optical microresonators coupled with the waveguide. (b). Equivalent CMT presentation for PhC layout. (c). Electric field profile in the symmetry protected BSC.

Download Full Size | PDF

If the waveguide can support only the symmetrical propagating mode then injected light can not probe the BSC. However that is true only for the linear case. If the microresonators are nonlinear owing to the Kerr effect there is a nonlinear shift of the eigenfrequency ω0 → ω0 + λ|Aj|2 where |Aj|2 is the intensity of light in the j-th microresonator. Substituting Eq. (1) into the CMT equations (2) we obtained the CMT equations modified with account of the Kerr effect as follows

iA˙1=(ω0+λ|A1|2)A1+iγ(A1+A2)+uA2iγEin(t)eiωt,iA˙2=(ω0+λ|A2|2)A2+iγ(A1+A2)+uA1iγEin(t)eiωt,
where λ is the nonlinear coefficient [35, 40]
λ=ω0n0c2n228πa2Em4(x,y)dxdy,Nm=ε(x,y)Em2(x,y)dxdy=a2cn2
with integration over the cross-section of the defect rod. Em(x,y) is the monopole eigenmode normalized via the constant Nm. The profile of the monopole mode is clearly seen in the BSC mode formed by two defect rods shown in Fig. 1. We present estimations based on the PhC structure in Fig. 1(a) which consists of GaAs rods with n 0 = 4.3,a = 0.5μm, radius 0.18μm. The single row of the rods is removed from the PhC to form a directional photonic waveguide which supports a band of guided even TM mode spanning from 0.302 to the upper band edge 0.444 [41]. Two nonlinear defect rods with the same radius 0.18μm, n 0 = 2 and n 2 = 2 · 10−12 W/cm 2 shown in gray are inserted symmetrically relative to the waveguide to form identical off-channel optical microresonators coupled with the waveguide. Estimations give the following values of the relevant parameters for the equivalent CMT layout: ω 0 = 0.352, γ ≈ 0.01, u ≈ 0.001, λ ≈ 10−4.

We rewrite Eq. (3) by use of As,a(t) = as a(t)eiωt as follows

ia˙s=(ν+λ[|as|2+2|aa|2])as+λaa2as*+2iγasiγEin(t),ia˙a=(ν+λ[|aa|2+2|as|2])aa+λas2aa*.
where ν = ω 0ω is the detuning parameter. One can see from the above equations that although the antisymmetric mode aa is decoupled from the injected light it can be excited via the nonlinear terms.

In order to trap light we apply a gaussian pulse of light

Ein(t)=E0exp(t2/2σ2).

Figure 2(a) shows that the injected gaussian pulse excites the symmetric mode as which follows the amplitude of the injected light pulse. When the symmetric mode is excited sufficiently it triggers excitation of the antisymmetric mode through the nonlinear terms in Eq. (5). After the pulse has passed by the resonator the symmetric mode leaks into the waveguide and the antisymmetric mode becomes isolated from the waveguide, i.e., becomes a true BSC with trapped amplitude |aa(∞)| as shown in Fig. 2(a). Figures 2(b) and 2(c) show that there are windows in the amplitude and frequency of gaussian pulse in which light storage takes place. The duration of gaussian pulse σ is to be of order 1/ν or larger to result in high intensity of the captured light as shown in Fig. 2(d). Therefore the robust symmetric BSC is capable for light trapping, however the amount of storage energy proportional to |aa(∞)|2 [35] is highly sensitive to the parameters of the injected gaussian pulse as seen from Figs. 2(b)–2(d). The application of the secondary pulse again opens the BSC according to Eq. (5) and the symmetry protected BSC leaks into the waveguide as shown in Fig. 2(e). The parameters of the secondary gaussian pulse σ = 9885, E 0 = 0.7741, ν = 0.0001 give optimal effect of release. One can also see high dynamical response of the BSC mode when the secondary pulse is passed by [42].

 figure: Fig. 2

Fig. 2 (a). Amplitudes of eigenmodes of resonators under application of a gaussian pulse in the PhC waveguide shown in Fig. 1(a). Green dash line shows the profile of the injected pulse with parameters: E 0 = 0.575, ν = 10−4, σ = 10000, dash-dot line shows the symmetric mode amplitude |as| and red solid line shows the amplitude of the antisymmetric |aa| BSC mode. The amplitude of trapped light aa(∞) in the symmetry protected BSC vs: (b). Amplitude E 0, (c) The detuning parameter ν, and (d) The duration of impulse σ for E 0 = 0.575 of the gaussian pulse Eq. (6). (e). The release of light from BSC after application of a secondary gaussian pulse with σ = 10797, ν = 0.0001, E 0 = 0.7741. The parameters of the resonators are the following: with the parameters ω 0 = 0.352, γ = 0.01, u = 0.001, λ = 0.0001.

Download Full Size | PDF

3. Single microcavity

Miniaturization of photonic elements and efficient transport of light are two paramount issues required for the design of integrated photonic circuits. To address the above issues instead of two microresonators we use a single microresonator with two eigenmodes embedded into the propagation band of the waveguide. Then the localization occurs due to a full destructive interference of corresponding resonant modes [43]. Specifically, we exploit the design shown in Fig. 3(a) [25]. Four linear defect rods of the same radius but with dielectric constant ε different from the dielectric constant of host rods shown by green open circles are placed at the vertexes of a square. On the both sides of the resonator four additional rods are inserted in the waveguide in order to suppress the coupling of the resonator with the waveguide. Totally, these nine defect rods form a microresonator embedded into the PhC waveguide whose eigenfrequencies versus the dielectric constant ε are plotted in Fig. 3(b) together with the eigenmodes [41]. Note, both eigenmodes are symmetric relative to inversion of the y-axis shown in Fig. 3(a), and therefore are coupled with the symmetrical propagating mode in the waveguide channel. The eigenfrequencies of the monopole and quadrupole-diag modes cross while effect of other eigenmodes, two dipole and quadrupole-xy modes are not important [25] and are not shown in Fig. 3(b).

 figure: Fig. 3

Fig. 3 (a). The same PhC structure as shown in Fig 1(a) however with single in-channel microresonator. Four linear defect rods of the same radius with dielectric constant ε are shown by green circles. (b). The eigenfrequencies of two eigenmodes (monopole and quadrupole-diag) shown as dependent on the dielectric constant of the central defect rod. (c). Light trapping by the BSC for the first pulse with amplitude E 0 = 0.2 and duration 2 × 104. The release of light by the secondary pulse with amplitude E 0 = 0.006 and duration 2 × 104. Inset shows details of excitation dynamics of eigenmode amplitudes after application of the first pulse. The parameters substituted into the CMT equations (9) are taken from Ref. [25]: ω 1 = 0.36, ω 2 = 0.365, γ 1 = 3 · 10−5, γ 2 = 1.3 · 10−4, u = 1.77 · 10−4, λ 11 = 10−4.

Download Full Size | PDF

Assume, we can vary the dielectric constant ε of the defect rods. Then as shown in Fig. 3(b) the eigenfrequencies are crossing. It is important to note that both eigenmodes have the same parity as the parity of electromagnetic mode propagating over the waveguide. Therefore each mode would decay into the waveguide becoming a resonant mode with a finite resonant width. The reader can find specific values of the resonant widths in Ref. [25]. When the eigenfrequencies cross the resonances undergo avoiding crossing resulting in that one of the resonant widths turns to zero forming a BSC [25, 43]. Experimentally this scenario for BSC occurrence meets with difficulties of tuning the dielectric constant of defect rods that makes the BSC non robust [23, 25]. It is remarkable that the Kerr effect in the microresonator removes this difficulty. The robust BSC appears in a self-adaptive way due to the nonlinear shift of the dielectric constant [25, 31–33 ]. Therefore the design in Fig. 3 offers a capacity for light storage.

The temporary CMT equations for the relevant modes of the resonator take the following form with account of instantaneous Kerr effect

iA˙1=(ω1+λ11|A1|2+iγ1)A1+(u+λ12(A1A2*+A1*A2)+iγ1γ2)A2iγ1Ein(t),iA˙2=(ω2+λ22|A1|2+iγ2)A2+(u+λ12(A1A2*+A1*A2)+iγ1γ2)A1iγ2Ein(t),
where the subscripts 1, 2 refer to the amplitudes of monopole and quadrupole-diag eigenmodes, and nonlinear coefficients equal [35, 40]
λmn=ω0n0c2n228πa2Em2(x,y)En2(x,y)dxdy
with normalization constants given by Eq. (4). The eigenmodes Em(x,y) are shown in Fig. 3(b). Integration includes host rods within a region of localization of the eigenmodes Em(x,y).

For the stationary process Ein (t) = E 0 eiωt, Aj(t) = Ajeiωt the CMT equations (7) takes the following form

(ω1+λ11|A1|2ω+iγ1u+λ12(A1A2*+A1*A2)+iγ1γ2u+λ12(A1A2*+A1*A2)+iγ1γ2ω2+λ22|A2|2ω+iγ2)(A1A2)=i(γ1γ2)E0.

The BSC occurs when the inverse of the matrix in Eq. (9) does not exist, i.e., when, the determinant of the matrix and respectively, one of the complex eigenvalues turns to zero. For the linear case formation of the BSC requires tuning of material parameters and occurs at [25]

ω2ω1=u(γ2γ1)γ1γ2,ωBSC=ω1+uγ2γ1,

For the nonlinear case Eq. (10) can be satisfied in a self-induced way [25,31] without necessity to tune material parameters: ωjωj +λjj|Aj|2, uu+λ12(A1A2*+A1*A2).

Two factors substantially weaken the nonlinear contribution of the quadrupole-diag mode. First, as seen from Fig. 3(a) two nodal lines of the quadrupole mode go through the central defect rod which enormously decreases the integrals in Eq. (8). Second, γ 1γ 2 that suppresses excitation of the quadrupole-diag mode amplitude A 2 compared to the monopole mode amplitude A 1. Therefore for simplicity we take λ 12 = 0, λ 22 = 0 and obtain from Eq. (10) that the BSC is achieved at fixed values of the mode intensities [25]

λ11|A1c|2=ω2ω1+u(γ2γ1)γ1γ2,|A2c|2=γ1γ2|A1c|2
when the intensities of monopole and quadrupole-diag modes become equal, respectively,
I2c=|A2c|2=γ1γ2I1c,I1,cI2c,ωBSCω2uγ1γ2.

Thus the amplitudes Am of the BSC eigenmode are fixed and depend only on the design of the optical microresonator and material parameters. Therefore, the amount of trapped light by the BSC is also fixed. The amount does not depend on the basic frequency and profile of a light pulse. Therefore for successful light storage the amount of pulse power should be enough to populate the BSC with the intensity given by Eqs. (11) and (12). The process of trapping numerically computed and shown in Fig. 3(c) agrees with this consideration. One can see that at first both mode amplitudes Am,m = 1,2 follow the pulse. After the time of order 104 the amplitudes stick to the BSC amplitudes given by Eq. (11) shown by dash lines in Fig. 3(c). The light accumulated in the resonator remains trapped until the next gaussian pulse is applied.

4. Fabry-Perot BSC

The next design is close to that considered in section 2, however two defect nonlinear rods are positioned at some distance L between each other as shown in Fig. 4(a). Each defect rods presents off-channel nonlinear microresonator which gives rise to a full reflection at the frequency dependent on the intensity of light [44]. Therefore these two off-channel resonators form a Fabry-Perot resonator. Then if the distance between the microresonators is tuned to fulfill the integer half-wavelength condition [32] the system supports a BSC. For simplicity we consider that only defect rods shown by gray in Fig. 4(a) are subjected to the Kerr effect and neglect the nonlinear phase shift in the waveguide. The CMT equations have the following form [13, 45]

iA˙1=(ω0+λ|A1|2+iγ)A1+iγA2eiθiγEin(t)eiωt,iA˙2=(ω0+λ|A2|2+iγ)A2+iγA1eiθiγEin(t)eiωt+iθ,
where θ is the phase shift incurred as the wave travels from the first defect to the second one. For a monochromatic wave θ = k(ω)L. However for a light pulse which is an expansion over monochromatic waves these equations (13) are not valid. Therefore we have to modify Eq. (13) taking into consideration the delay time τ for travelling from one microcavity to another:
iA˙1=(ω0+λ|A1|2+iγ)A1+iγA2(tτ)iγEin(t)eiωt,iA˙2=(ω0+λ|A2|2+iγ)A2+iγA1(tτ)iγEin(tτ)eiω(tθ).

 figure: Fig. 4

Fig. 4 (a). Layout of 2D PhC which forms the Fabry-Perot resonator. Two defect rods shown by gray with the Kerr effect are inserted near by PhC waveguide. (b). Equivalent CMT presentation for PhC layout with equivalent parameters: γ = 0.01, λ = 0.0001, ω 0 = 0.352, u = 0.

Download Full Size | PDF

For stationary process with aj(t) = Ajeiωt, θ = τω. Equation (14) can be obtained from Eq. (13) after the substitution ω with the differential operator it. Then we have Aj(tτ)=exp[iτt]Aj(t)=exp(iθ)Aj(t). The BSC occurs if the determinant of the matrix in the stationary CMT equations

(ω0+λ|A1|2ω+iγiγeiτωiγeiτωω0+λ|A2|2ω+iγ)
equals zero [32]. That condition gives us
τω=πn,|A1c|2=|A2c|2=1λ(πnτω0),
where n is an integer. For the parameters listed in Fig. 4 and ω = ω 0 +ν, τ = π/ω we obtain the amplitudes of the intensities of the BSC take discrete value defined by the ratio νλ for n = 1. They are are shown in Fig. 5(a) by black dash lines. The numerical solution of Eq. (14) shows as injected gaussian pulse is captured at this BSC as shown in Fig. 5(a) provided that the amplitude of the injected pulse exceeds the threshold similar to the previous design with single microresonator. Application of the secondary pulse releases captured light. In Fig. 5(b) we present the most optimal process which occurs at the parameters of the gaussian pulse listed in figure caption. As different from the case of single microresonator in section 3, the BSC intensities are discretized by integer n as given by Eq. (15). Therefore application of larger pulse intensities can give rise to trapping of light in the next BSCs with n > 1.

 figure: Fig. 5

Fig. 5 Time evolution of mode amplitudes |A 1| (blue solid line) and |A 2| (red dash line) after the injection of a first gaussian pulse shown by green dash-dot line with σ = 10000, E 0 = 0.3522, ν = 10−5 (a) and a secondary pulse with σ = 4705, E 0 = 0.545, ν = 10−5. Black dash lines show the BSC amplitude given by Eq. (15).

Download Full Size | PDF

5. Discussion and conclusions

In regard to stability the approaches based on stopping light in linear PhC designs briefly discussed in the Introduction are preferable compared to the light storage in BSCs in nonlinear PhC structures. Moreover the linear system can operate at any light intensity while the light storage in the self-induced BSCs is possible only if the amplitude of the injected pulse exceeds a certain threshold sufficient to populate discrete value of the BSC intensity. The process of light slowing demands tuning of operational frequency of light to the specific points of the frequency spectrum of waveguide where the group velocity of light tends to zero. While the only requirement in the present approach is that the BSC frequencies reside in the propagation band of PhC waveguide. The next fundamental advantage of the BSC designs is a ultra-high compactness. The trapping of light in a single defect rod due to a full destructive interference of the monopole and quadrupole-diag modes [43] occurs in the volume of order of a 2 h where a is the period of 2D PhC and h is the length of the rods. While the former linear designs for light storage require much more volume. The second advantage of exploiting the self-induced BSC is that the light can be easily released by a secondary pulse rendering the process of storage and release by demand all-optical. Instead of the Kerr effect resulting in self-induced BSCs we can utilize a local heating or local deformation of the PhC to reach the condition for the BSC in controllable manner. However these approaches deprive the advantage of all-optical manipulation.

In the present paper we considered the phenomenon of light storage and release by the bound states in the continuum in the framework of temporal coupled mode theory with the parameters borrowed from numerical calculations for real 2D PhC structures. Direct study of these processes in PhC designs requires extremely time consuming computations even for the case when the nonlinearity is accounted for only the defect rods. The temporal CMT allows us easily asses merits and disadvantages of the PhC designs considered in sections 2–4. The first design consists of a waveguide with two off-channel nonlinear microresonators weakly coupled with the waveguide. The design holds the antisymmetric mode shown in Fig. 1(c) which is decoupled from the symmetric propagating mode and thereby is determined as a symmetry protected bound state. Application of a light pulse damages this BSC as the pulse passes by the microresonators because of nonlinear interaction of the BSC with the symmetrical mode as seen from Eq. (5). As the result some amount of light is trapped in the BSC after the pulse has been passed. However the next pulse is again coupled with the BSC giving rise to release of captured light. Figures 2(b)–2(d) show that the light storage is highly sensitive to the parameters of the gaussian pulse. That is result of the BSC intensity or equivalently the BSC frequency depends on to which extend the BSC is populated.

In the next section we considered different PhC design with a single in-channel microresonator in which the BSC is a result of full destructive interference of two resonant modes. In the linear case the occurrence of BSC would require tuning of the material parameters. It is remarkable that due to the Kerr effect BSCs occur by self-induced way as it was shown in Refs. [25,31,33] if the intensity of accumulated light achieves critical values given by Eq. (12). Finally we considered in the framework of the temporal CMT the Fabry-Perot PhC structure shown in Fig. 4. The temporal CMT equations were modified to account for the delay times to yield Eq. (14). This design also supports self-induced BSCs as was shown in Ref. [32] however with multiple discretized BSC intensities as given by Eq. (15). That gives more freedom in light storage compared to the previous design.

There are processes of irreversible light emission due to finite thickness of 2D PhC as well as due to imperfectness of the structure which take place for any approach. That will restrict the storage in time. Also the nonlinear resonator has a risk of bifurcation into modes coupled with different waveguide channels or radiation modes for example because of modulation instability of enough long nonlinear defect rod. These processes can be readily involved into the CMT approach by expansion of number of modes of the microresonator and number of channels coupled with the modes. But it is clear that a criteria to exclude these unwanted modes is the frequency distance between useful modes and unwanted modes should substantially exceed the perturbation due to the Kerr effect. In the present paper we consider this criteria is fulfilled.

Let us estimate the values of the power of injected light which are necessary to achieve the all-optical storage and release by demand. The nonlinearity constant is estimated in order of magnitude as 10−4. Then the injected amplitude in terms of CMT equations E 0 = 1 corresponds to 625W/a where a = 0.5μm is the lattice unit. Therefore in order to store light in the symmetry protected BSC (section 2) by gaussian pulse its power should exceed E02625W/a. For the length of rods compared to a we obtain that power injected into the waveguide channel should be around 180 W. However this power can be reduced further when the tuning parameter ν is decreased as it is seen from Fig. 2(b). Similar estimates can be applied for other PhC designs in Sections 3 and 4. However these designs demand that the injected amplitude exceed the threshold values. For the single microresonator in section 3 the injected power should exceed 25W. For the Fabry-Perot design we obtain the value around 100 W however this power also can be reduced with decereasing of the tuning parameter ν. It is very important to notice the Kerr effect is small λ ∼ 10−4. Therefore to compensate for small λ we have to tune ν to be the same order in magnitude as λ. The same refers to the duration of the gaussian pulse σ which has to be of order 1/ν.

Acknowledgments

This work was supported by Russian Scientific Foundation through Grant 14-12-00266. We acknowledge discussions with D.N. Maksimov.

References and links

1. T. F. Krauss, “Slow light in photonic crystal waveguides,” J. Phys. D 40, 2666–2670 (2007). [CrossRef]  

2. T. Baba, “Slow light in photonics crystals,” Nat. Photonics , 2, 465–473 (2008). [CrossRef]  

3. M. Notomi, “Manipulating light with strongly modulated photonic crystals,” Rep. Progr. Phys. 73, 096501 (2010). [CrossRef]  

4. A. Melloni, F. Morichetti, and M. Martinelli, “Optical slow wave structures,” Opt. Photon. News 14, 44–48 (2003). [CrossRef]  

5. F. Morichetti, C. Ferrari, A. Canciamilla, and A. Melloni, “The first decade of coupled resonator optical waveguides: bringing slow light to applications,” Laser Photonics Rev. 6, 74–96 (2012). [CrossRef]  

6. H. Gersen, T. J. Karle, R. J. P. Engelen, W. Bogaerts, J. P. Korterik, N. F. van Hulst, T. F. Krauss, and L. Kuipers, “Real-space observation of ultraslow light in Photonic Crystal Waveguides,” Phys. Rev. Lett. 94, 073903 (2005). [CrossRef]   [PubMed]  

7. I. V. Shadrivov, A. A. Sukhorukov, and Y. S. Kivshar, “Guided modes in negative-refractive-index waveguides,” Phys. Rev. E 67, 057602 (2003). [CrossRef]  

8. J. He, Y. Jin, Z. Hong, and S. He, “Slow light in a dielectric waveguide with negative-refractive-index photonic crystal cladding,” Opt. Express , 16, 11077–11082 (2008). [CrossRef]   [PubMed]  

9. S. He, Y. He, and Y. Jin, “Revealing the truth about trapped rainbow storage of light in metamaterials,” Sci. Rep. 583, 1–9 (2012).

10. K. L. Tsakmakidis, A. D. Boardman, and O. Hess, “Trapped rainbow storage of light in metamaterials,” Nature 450, 397–401 (2007). [CrossRef]   [PubMed]  

11. V. N. Smolyaninova, I. I. Smolyaninov, A. V. Kildishev, and V. M. Shalaev, “Experimental observation of the trapped rainbow,” Appl. Phys. Lett. 96, 211121 (2010). [CrossRef]  

12. R. Yang, W. Zhu, and J. Li, “Realization of trapped rainbow in 1D slab waveguide with surface dispersion engineering,” Opt. Express 23(5), 6326–6335 (2015). [CrossRef]   [PubMed]  

13. S. Fan, P. R. Villeneuve, J. D. Joannopoulos, M. J. Khan, C. Manolatou, and H. A. Haus, “Theoretical analysis of channel drop tunneling processes,” Phys. Rev. B 59(24), 15882–15892 (1999). [CrossRef]  

14. S. P. Shipman and S. Venakides, “Resonant transmission near nonrobust periodic slab modes,” Phys. Rev. E 71, 026611 (2005). [CrossRef]  

15. D. C. Marinica, A. G. Borisov, and S. V. Shabanov, “Bound States in the Continuum in Photonics,” Phys. Rev. Lett. 100, 183902 (2008). [CrossRef]   [PubMed]  

16. E. N. Bulgakov and A. F. Sadreev, “Bound states in the continuum in photonic waveguides inspired by defects,” Phys. Rev. B 78, 075105 (2008). [CrossRef]  

17. S. Longhi, “Transfer of light waves in optical waveguides via a continuum,” Phys. Rev. A 78, 013815 (2008). [CrossRef]  

18. S. Longhi, “Optical analogue of coherent population trapping via a continuum in optical waveguide arrays,” J. Mod. Opt. 56, 729–737 (2009). [CrossRef]  

19. F. Dreisow, A. Szameit, M. Heinrich, R. Keil, S. Nolte, A. Tünnermann, and S. Longhi, “Adiabatic transfer of light via a continuum in optical waveguides,” Opt. Lett. 34, 2405–2407 (2009). [CrossRef]   [PubMed]  

20. Y. Plotnik, O. Peleg, F. Dreisow, M. Heinrich, S. Nolte, A. Szameit, and M. Segev, “Experimental observation of optical bound states in the continuum,” Phys. Rev. Lett. 107, 183901 (2011). [CrossRef]   [PubMed]  

21. G. Corrielli, G. DellaValle, A. Crespi, R. Osellame, and S. Longhi, “Observation of surface states with algebraic localization,” Phys. Rev. Lett. 111, 220403 (2013). [CrossRef]   [PubMed]  

22. S. Weimann, Y. Xu, R. Keil, A. E. Miroshnichenko, A. Tunnermann, S. Nolte, A. A. Sukhorukov, A. Szameit, and Yu. S. Kivshar, “Compact surface Fano states embedded in the continuum of waveguide arrays,” Phys. Rev. Lett. 111, 240403 (2013). [CrossRef]  

23. C. W. Hsu, B. Zhen, S.-L. Chua, S. G. Johnson, J. D. Joannopoulos, and M. Soljačić, “Bloch surface eigenstates within the radiation continuum,” Light: Sci. Appl. 2, 1–5 (2013).

24. C. W. Hsu, B. Zhen, J. Lee, S.-L. Chua, S. G. Johnson, J. D. Joannopoulos, and M. Soljačić, “Observation of trapped light within the radiation continuum,” Nature 499, 188–191 (2013). [CrossRef]   [PubMed]  

25. E. N. Bulgakov and A. F. Sadreev, “Robust bound state in the continuum in a nonlinear microcavity embedded in a photonic crystal waveguide,” Opt. Lett. 39(17), 5212–5215 (2014). [CrossRef]   [PubMed]  

26. Y. Yang, C. Peng, Y. Liang, Z. Li, and S. Noda, “Analytical perspective for bound states in the continuum in photonic crystal slabs,” Phys. Rev. Lett. 113, 037401 (2014). [CrossRef]   [PubMed]  

27. B. Zhen, C. W. Hsu, L. Lu, A. D. Stone, and M. Soljačić, “Topological nature of bound states in the radiation continuum,” Phys. Rev. Lett. 113, 257401 (2014). [CrossRef]  

28. E. N. Bulgakov and A. F. Sadreev, “Bloch bound states in the radiation continuum in a periodic array of dielectric rods,” Phys. Rev. A 90, 053801 (2014). [CrossRef]  

29. M. Song, H. Yu, C. Wang, N. Yao, M. Pu, J. Luo, Z. Zhang, and X. Luo, “Sharp Fano resonance induced by a single layer of nanorods with perturbed periodicity,” Opt. Express 23(3), 2895–2903 (2015). [CrossRef]   [PubMed]  

30. M. Zhang and X. Zhang, “Ultrasensitive optical absorption in graphene based on bound states in the continuum,” Sci. Rep. 5, 8266 (2015). [CrossRef]   [PubMed]  

31. E. N. Bulgakov and A. F. Sadreev, “Resonance induced by a bound state in the continuum in a two-level nonlinear Fano-Anderson model,” Phys. Rev. B 80, 115308 (2009); [CrossRef]  

32. E. N. Bulgakov and A. F. Sadreev, “Bound states in photonic Fabry-Perot resonator with nonlinear off-channel defects,” Phys. Rev. B 81, 115128 (2010). [CrossRef]  

33. M. I. Molina, A. E. Miroshnichenko, and Yu. S. Kivshar, “Surface bound states in the continuum,” Phys. Rev. Lett. 108, 070401 (2012). [CrossRef]   [PubMed]  

34. E. N. Bulgakov, K. N. Pichugin, and A. F. Sadreev, “Symmetry breaking for transmission in a photonic waveguide coupled with two off-channel nonlinear defects,” Phys. Rev. B 83, 045109 (2011). [CrossRef]  

35. J. Joannopoulos, R. D. Meade, and J. Winn, Photonic Crystals: Molding the Flow of Light (Princeton University Press, 1995).

36. D. Englund, I. Fushman, and J. Vuckovic, “General recipe for designing photonic crystal cavities,” Opt. Express 13, 5961–5975 (2005). [CrossRef]   [PubMed]  

37. M. Notomi, E. Kuramochi, and H. Taniyama, “Ultrahigh-Q nanocavity with 1D photonic gap,” Opt. Express 16, 11096–11102 (2009).

38. Sh. 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]  

39. W. Suh, Z. Wang, and S. Fan, “Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities,” IEEE J. Quantum Electron. 40(10), 1511–1518 (2004). [CrossRef]  

40. E. N. Bulgakov and A. F. Sadreev, “All-optical diode based on dipole modes of Kerr microcavity in asymmetric L-shaped photonic crystal waveguide,” Opt. Lett. 39(7), 1787–1790 (2014). [CrossRef]   [PubMed]  

41. K. Busch, S. F. Mingaleev, A. Garcia-Martin, M. Schillinger, and D. Hermann, “The Wannier function approach to photonic crystal circuits,” J. Phys.: Cond. Mat. 15, R1233–R1256 (2003).

42. K. N. Pichugin and A. F. Sadreev, “Frequency comb generation for wave transmission through the nonlinear dimer,” J. Phys.: Conf. Ser. 574012037 (2015).

43. H. Friedrich and D. Wintgen, “Interfering resonances and bound states in the continuum,” Phys. Rev. A 32, 3231–3242 (1985). [CrossRef]   [PubMed]  

44. A. E. Miroshnichenko, S. F. Mingaleev, S. Flach, and Yu. S. Kivshar, “Nonlinear Fano resonance and bistable wave transmission,” Phys. Rev. E 71, 036626 (2005). [CrossRef]  

45. Z. Wang and S. Fan, “Compact all-pass filters in photonic crystals as the building block for high-capacity optical delay lines,” Phys. Rev. E 68, 066616 (2003). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1
Fig. 1 (a). The single row of the rods is removed from the PhC to form a directional photonic waveguide which supports a band of guided even TM mode. Two defect rods shown in gray with the Kerr effect are inserted symmetrically relative to the waveguide to form identical off-channel optical microresonators coupled with the waveguide. (b). Equivalent CMT presentation for PhC layout. (c). Electric field profile in the symmetry protected BSC.
Fig. 2
Fig. 2 (a). Amplitudes of eigenmodes of resonators under application of a gaussian pulse in the PhC waveguide shown in Fig. 1(a). Green dash line shows the profile of the injected pulse with parameters: E 0 = 0.575, ν = 10−4, σ = 10000, dash-dot line shows the symmetric mode amplitude |as | and red solid line shows the amplitude of the antisymmetric |aa | BSC mode. The amplitude of trapped light aa (∞) in the symmetry protected BSC vs: (b). Amplitude E 0, (c) The detuning parameter ν, and (d) The duration of impulse σ for E 0 = 0.575 of the gaussian pulse Eq. (6). (e). The release of light from BSC after application of a secondary gaussian pulse with σ = 10797, ν = 0.0001, E 0 = 0.7741. The parameters of the resonators are the following: with the parameters ω 0 = 0.352, γ = 0.01, u = 0.001, λ = 0.0001.
Fig. 3
Fig. 3 (a). The same PhC structure as shown in Fig 1(a) however with single in-channel microresonator. Four linear defect rods of the same radius with dielectric constant ε are shown by green circles. (b). The eigenfrequencies of two eigenmodes (monopole and quadrupole-diag) shown as dependent on the dielectric constant of the central defect rod. (c). Light trapping by the BSC for the first pulse with amplitude E 0 = 0.2 and duration 2 × 104. The release of light by the secondary pulse with amplitude E 0 = 0.006 and duration 2 × 104. Inset shows details of excitation dynamics of eigenmode amplitudes after application of the first pulse. The parameters substituted into the CMT equations (9) are taken from Ref. [25]: ω 1 = 0.36, ω 2 = 0.365, γ 1 = 3 · 10−5, γ 2 = 1.3 · 10−4, u = 1.77 · 10−4, λ 11 = 10−4.
Fig. 4
Fig. 4 (a). Layout of 2D PhC which forms the Fabry-Perot resonator. Two defect rods shown by gray with the Kerr effect are inserted near by PhC waveguide. (b). Equivalent CMT presentation for PhC layout with equivalent parameters: γ = 0.01, λ = 0.0001, ω 0 = 0.352, u = 0.
Fig. 5
Fig. 5 Time evolution of mode amplitudes |A 1| (blue solid line) and |A 2| (red dash line) after the injection of a first gaussian pulse shown by green dash-dot line with σ = 10000, E 0 = 0.3522, ν = 10−5 (a) and a secondary pulse with σ = 4705, E 0 = 0.545, ν = 10−5. Black dash lines show the BSC amplitude given by Eq. (15).

Equations (16)

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

A s , a = ( A 1 ± A 2 ) / 2 , ω s , a = ω 0 ± u ,
i A ˙ s = ( ω s + 2 i γ ) A s i γ E i n e i ω t , i A ˙ a = ω a A a ,
i A ˙ 1 = ( ω 0 + λ | A 1 | 2 ) A 1 + i γ ( A 1 + A 2 ) + u A 2 i γ E i n ( t ) e i ω t , i A ˙ 2 = ( ω 0 + λ | A 2 | 2 ) A 2 + i γ ( A 1 + A 2 ) + u A 1 i γ E i n ( t ) e i ω t ,
λ = ω 0 n 0 c 2 n 2 2 8 π a 2 E m 4 ( x , y ) d x d y , N m = ε ( x , y ) E m 2 ( x , y ) d x d y = a 2 c n 2
i a ˙ s = ( ν + λ [ | a s | 2 + 2 | a a | 2 ] ) a s + λ a a 2 a s * + 2 i γ a s i γ E i n ( t ) , i a ˙ a = ( ν + λ [ | a a | 2 + 2 | a s | 2 ] ) a a + λ a s 2 a a * .
E i n ( t ) = E 0 exp ( t 2 / 2 σ 2 ) .
i A ˙ 1 = ( ω 1 + λ 11 | A 1 | 2 + i γ 1 ) A 1 + ( u + λ 12 ( A 1 A 2 * + A 1 * A 2 ) + i γ 1 γ 2 ) A 2 i γ 1 E i n ( t ) , i A ˙ 2 = ( ω 2 + λ 22 | A 1 | 2 + i γ 2 ) A 2 + ( u + λ 12 ( A 1 A 2 * + A 1 * A 2 ) + i γ 1 γ 2 ) A 1 i γ 2 E i n ( t ) ,
λ m n = ω 0 n 0 c 2 n 2 2 8 π a 2 E m 2 ( x , y ) E n 2 ( x , y ) d x d y
( ω 1 + λ 11 | A 1 | 2 ω + i γ 1 u + λ 12 ( A 1 A 2 * + A 1 * A 2 ) + i γ 1 γ 2 u + λ 12 ( A 1 A 2 * + A 1 * A 2 ) + i γ 1 γ 2 ω 2 + λ 22 | A 2 | 2 ω + i γ 2 ) ( A 1 A 2 ) = i ( γ 1 γ 2 ) E 0 .
ω 2 ω 1 = u ( γ 2 γ 1 ) γ 1 γ 2 , ω B S C = ω 1 + u γ 2 γ 1 ,
λ 11 | A 1 c | 2 = ω 2 ω 1 + u ( γ 2 γ 1 ) γ 1 γ 2 , | A 2 c | 2 = γ 1 γ 2 | A 1 c | 2
I 2 c = | A 2 c | 2 = γ 1 γ 2 I 1 c , I 1 , c I 2 c , ω B S C ω 2 u γ 1 γ 2 .
i A ˙ 1 = ( ω 0 + λ | A 1 | 2 + i γ ) A 1 + i γ A 2 e i θ i γ E i n ( t ) e i ω t , i A ˙ 2 = ( ω 0 + λ | A 2 | 2 + i γ ) A 2 + i γ A 1 e i θ i γ E i n ( t ) e i ω t + i θ ,
i A ˙ 1 = ( ω 0 + λ | A 1 | 2 + i γ ) A 1 + i γ A 2 ( t τ ) i γ E i n ( t ) e i ω t , i A ˙ 2 = ( ω 0 + λ | A 2 | 2 + i γ ) A 2 + i γ A 1 ( t τ ) i γ E i n ( t τ ) e i ω ( t θ ) .
( ω 0 + λ | A 1 | 2 ω + i γ i γ e i τ ω i γ e i τ ω ω 0 + λ | A 2 | 2 ω + i γ )
τ ω = π n , | A 1 c | 2 = | A 2 c | 2 = 1 λ ( π n τ ω 0 ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.