Bistability, multistability and nonreciprocal light propagation in Thue-Morse multilayered structures
- 1 -
Bistability, multistability and nonreciprocal light propagation in Thue-Morse multilayered structures
Victor Grigoriev and Fabio Biancalana
Nonlinear Photonic Nanostructures Research Group, Max Planck Institute for the Science of Light, Guenther-Scharowsky-Str. 1 / Bau 26 (R440), Erlangen 91058, Germany E-mail: [email protected]
Abstract . The nonlinear properties of quasiperiodic photonic crystals based on the Thue-Morse sequence are investigated. The intrinsic spatial asymmetry of these one-dimensional structures for odd generation numbers results in bistability thresholds which are sensitive to the propagation direction. Along with resonances of perfect transmission, this feature allows to achieve strongly non-reciprocal propagation and to create an all-optical diode. The salient qualitative features of such optical diode action is readily explained through a simple coupled resonator model. The efficiency of a passive scheme, which does not necessitate of an additional short pump signal, is compared to an active scheme, where such a signal is required.
PACS : 42.65.Pc, 78.67.Pt, 42.79.Ta
1. Introduction
The experimental discovery of quasicrystals [1] has given rise to a sudden breakthrough in the area of solid state physics. It had been assumed for a long time that a periodic arrangement of unit cells was a main prerequisite for the specific properties of crystalline matter. The structure of quasicrystals does exhibit long-range order or correlations between distant parts, but at the same time there is no underlying periodicity, in the sense that a shifted copy of the crystal never matches exactly the original one. In fact, quasicrystals represent an intermediate stage between random media and periodic crystals, effectively combining both localization properties as a result of short-range disorder and the presence of band gaps due to long-range correlations [2]. The concept of quasiperiodicity was readily transferred to photonic crystals and proved to be of great value for many practical purposes [3-6]. Photonic quasicrystals are deterministically generated dielectric structures with a non-periodic modulation of the refractive index. In the one-dimensional (1D) case, they can be formed by stacking together dielectric layers of several different types according to the substitutional sequence under investigation (Cantor, Fibonacci, Rudin-Shapiro, Thue-Morse, etc.) [7]. The Fibonacci sequence is of particular importance, since it leads to the existence of two incommensurable periods in the spatial spectrum of the structure. Such behavior is typical of sequences with a so-called pure point spectrum, which makes the Fibonacci sequence truly quasiperiodic, as a consequence of the appearance of Bragg-like peaks in the spatial spectrum [8]. This - 2 - property has been demonstrated to be very valuable for nonlinear optical applications, such as third harmonic generation [9]. In fact, the latter is a two-stage process, and for each stage different phase matching conditions should be met. The Fibonacci sequence allows to fulfill both of them simultaneously and on the same crystal [10]. In contrast to the Fibonacci sequence, the Thue-Morse (ThM) sequence possesses a singular continuous spectrum, which is neither continuous nor singular [11-12]. Thus, strictly speaking, the ThM sequence is not quasiperiodic, but rather deterministically aperiodic. To an even different class belongs the Rudin-Shapiro sequence, which shows evidence of a continuous spectrum, analogous to the one exhibited by random sequences [13]. With a slight abuse of terminology, we shall always refer to the above sequences as quasiperiodic, and we shall focus our attention specifically to the ThM sequence, for reasons explained below. The bandgaps of quasiperiodic multilayered structures, also called ‘pseudo band gaps’ [14] or ‘fractal bandgaps’ in the literature [15], often contain resonant states, which can be considered as a manifestation of numerous defects distributed along the structure. Their quality factors are not very large when compared with those of symmetrically placed defects inside conventional Bragg structures, but the mode profiles associated to these resonances are extended in space and, as a consequence, very suitable for the enhancement of nonlinear effects throughout the whole body of the crystal. In order to investigate the specific properties of bistability and multistability in quasicrystals, in the following we shall restrict our attention exclusively to 1D structures. As a representative case, we focus our attention to the ThM aperiodic multilayer, and in making this choice we were guided by the following two reasons. Firstly, the spatial asymmetry of the ThM sequence (of odd generation numbers) interplays with the nonlinearity and makes the transmission sensitive to the propagation direction [16]. Secondly, almost all resonances in ThM quasicrystals are resonances of complete (i.e. 100%) transmission, with some of them being situated well inside the pseudo band gap regions. Taken together, the above two properties provide favorable conditions for the design of an all-optical diode, i.e. a device that shows a strong contrast in the transmission between forward and backward incidence. A proper understanding of bistability phenomena in quasiperiodic crystals is very important for the understanding of the optical diode action described above. The paper is organized as follows. In section 2, we briefly review the well-known linear properties of ThM quasicrystals and classify their resonances by using the method of trace maps. It is shown that the field profiles at resonance frequencies follow the pattern of Thue-Morse sequence and the classification suggested can be used as a measure of this self-similarity. In section 3, we show that the spatial asymmetry inherent to ThM quasicrystals of odd generation numbers can interplay with Kerr-nonlinearity in such a way that transmission become sensitive to the propagation direction. It is very important that the dynamics of this interplay can have a sudden jump due to bistability when the intensity of the incident light changes, and the thresholds for such jumps are different for the incidence from the left and from the right of the multilayer. In section 4, we apply a phenomenological model based on a coupled resonator model to explain this behavior analytically for both bistable and multistable cases. The general formulas describing nonlinear transmission spectra are derived and their relationship to the level of self-similarity in the field profiles is emphasized. In section 5, we propose a design of nonlinear optical diode based on the differences in bistability thresholds between forward and backward propagation. The efficiency of two schemes is compared: passive, when only one pulse is used, and active, when an additional short pump signal is applied to facilitate switching between stable branches of hysteresis. The last section summarizes the results and presents the conclusions. The appendix provides some details about the FDTD method used in the paper for the dynamical simulations of multilayered structures with instantaneous Kerr-nonlinearity.
2. Thue-Morse multilayers
There are two main approaches to generate 1D quasiperiodic sequences. The first makes use of a projection from a higher-dimensional space, while the second employs the so-called substitutional sequences [17]. Being strictly quasiperiodic, the Fibonacci sequence can be obtained by using both the - 3 - above methods, and the presence of two incommensurable periods in their spatial spectra is implicitly related to the projection of a two-dimensional grid onto a one-dimensional line [18]. However, these two approaches are not equivalent to each other, and the use of substitutional sequences tends to be a more general procedure [19]. In particular, this is applicable to ThM sequences and makes it distinct from strictly quasiperiodic sequences. Sometimes, it is called a deterministic aperiodic sequence, in order to emphasize that it has more disorder than quasiperiodic ones and stands closer to random sequences [20]. On the other hand, the development of resonances when changing generation number in ThM quasicrystals resembles the analogous development in periodic structures to a certain extent and brings these two kinds of crystal closer to each other [21]. Two dielectric layers of different materials are required to compose ThM quasicrystals, which shall be indicated with the letters 'A' and 'B' . They should be arranged in the same way as in literal ThM sequence, which is governed by the following substitutional (‘inflation’) rules: 'A' 'AB' → , 'B' 'BA' → . Thus, starting from a single layer 'A' S = , which is defined to be the ThM quasicrystal of 0th generation or ThM , one obtains 'AB' S = , 'ABBA' S = , 'ABBABAAB' S = and so forth, with each step giving a sequence of generation number increased by one (figure 1a). It can be shown that an additional recurrence relation follows from this definition, which holds for ThM sequences as single blocks: n n n S S S + = (cid:4) , where n is the generation number. In this notation n S (cid:4) indicates a sequence ‘conjugated’ to n S , where all letters are interchanged to their opposite as in the rule 'A' 'B' ↔ . It is interesting to note that there is only a slight difference between the ThM and Bragg (or periodic) sequences as far as inflation rules are concerned. More specifically, Bragg sequences possess the inflation rule 'A' 'AB' → , 'B' 'AB' → . Therefore, the growth rate of these sequences is the same, and the total number of layers is subject to a rapid (exponential) increase as n . The occurrence of Figure 1. ( a ) Application of inflation rules 'A' 'AB' → , 'B' 'BA' → for ThM sequence. Numbers near the notches designate the boundaries of sequences with corresponding generation numbers. ( b ) Field profile at one of the resonance frequencies ( / 0.756 Ω Ω = ) of a ThM quasicrystal showing self-similar pattern of ThM sequence. The frequency Ω corresponds to quarter wavelength condition taken at 0.7 μ m. The linear refractive indices of the materials used are 1.55 A n = (polydiacetylene 9-BCMU) and 2.3 B n = (rutile). Background shows alternation of layers inside the structure. - 4 - layers 'A' and 'B' is also the same for both ThM and Bragg structures for any generation number. Nevertheless, such a small modification of the inflation rules leads to completely distinct transmission and localization properties. The primary method to compute transmission spectra and field profiles for multilayered structures is the transfer matrix method [22]. If normal incidence is assumed, the electric E and magnetic H fields at the opposite sides of each single layer can be related in the following way: i i mm m E EH H − ⎛ ⎞ ⎛ ⎞=⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠ M , (1) cos( ) sin( )sin( ) cos( ) m m m m mm m m m m m n d n n dn n d n d − ⎡ ⎤Ω Ω= ⎢ ⎥− Ω Ω⎢ ⎥⎣ ⎦ M , (2) where m d is the thickness of the layer with refractive index m n , and / c Ω =ω is the wavenumber in vacuum. In the representation of equation (1) the imaginary unit is deliberately kept with the magnetic field, in order to show explicitly that transfer matrices of nonabsorbing layered structures can be characterized by real matrix elements. In general, the refractive indices of the materials can depend on frequency, but for the sake of simplicity we will assume that all materials are dispersionless. By multiplying together all transfer matrices of subsequent layers, it is possible to construct the total transfer matrix of the structure ...
L L − = M M M M M , (3) and to find the transmission spectrum in terms of its matrix elements
211 22 21 12
4( ) | i( ) |
T M M M M
Ω = + + − . (4) Another – more efficient – approach to obtain transfer matrices of ThM quasicrystals is to make use of the following two-step recurrence relations { 1} { } { } n n n + = M M M (cid:4) , (5) { 1} { } { } n n n + = M M M (cid:4) (cid:4) , (6) supported by the initial conditions {0} A = M M and {0} B = M M (cid:4) , which are the transfer matrices for layers 'A ' and 'B' , respectively. The subscripts in curly brackets indicate generation numbers of corresponding ThM quasicrystals. Formulas (5)-(6) take their origin in the recurrence relations of ThM sequence explained above, and stress explicitly the interplay between ‘ordinary’ n S and ‘conjugated’ n S (cid:4) counterparts of the sequence. The self-similarity of ThM quasicrystals is naturally embedded in this approach, and this simplifies considerably the computation of transmission spectra for large generation numbers. However, even this modified representation of the transfer matrix method still contains redundant information about the structure, if one is interested only in the analysis of the properties of the resonance frequencies. In this case, it is sufficient to operate with traces of transfer matrices. The trace is defined as
11 22
Tr( ) x M M = = + M . (7) For ThM quasicrystals of fixed generation numbers, these traces are unchanged under conjugation { } { 1} { 1} { 1} { 1} { } Tr( ) Tr( ) Tr( ) Tr( ) n n n n n n n x − − − − = = = = M M M M M M (cid:4) (cid:4) (cid:4) , (8) and for adjacent generation numbers, an additional relation or ‘trace map’ holds [23]:
Tr( ) Tr( ) ( 2) 2 n n n n n n n n n x x x + + = = = − +
M M M M M M (cid:4) (cid:4) (cid:4) . (9) - 5 -
The formulas (8)-(9) are valid for 1 n ≥ , and in derivation of (9) the Cayley-Hamilton theorem was applied to prove that any unimodular matrix satisfies the relation U x = − U U I , (10) where U x is the trace of a unimodular matrix U , and I is the identity matrix. To use the equation for the trace map (7), two initial conditions are required. They can be found directly from the multiplication of transfer matrices having the form given by equation (2): A B B A x n n n n = − + α β α β , (11) A B B A x n n n n = − + α β α β , (12) where
A A n d = Ω α and B B n d = Ω β . Several important conclusions can be extracted from the trace map of ThM quasicrystals. The analysis simplifies greatly in absence of absorption, so that the imaginary parts of the refractive indices are negligible. In this case, all traces will be real, and there will be frequencies in the spectrum in correspondence to which n x = . Independent of the value of n x + , this makes n x + = , which coincides with the trace of identity matrix. This simply means that frequencies for which 0 n x = correspond to frequencies of perfect transmission. To prove this, an auxiliary matrix should be constructed, the trace of which gives the difference of diagonal elements of the original transfer matrix:
11 22
Tr( ) z z M M = = − σ M , (13) where z σ is the third (diagonal) Pauli matrix. Provided that n x = , it can be readily shown that the diagonal elements of { 2} n + M are equal: Tr( ) Tr[ ] 0 n z n n n n z z + = = = σ M M M M σ (cid:4) (cid:4) . (14) Therefore, the resonances of complete transmission in ThM quasicrystals can be characterized by identity transfer matrices, and they are preserved from one generation to another. Only additional resonances can appear in new generations, and the overall transmission spectrum of these quasicrystals shows a fractal nature. As generation number increases, the spectrum reveals self-similarity and characteristic trifurcation of resonances [23]. In what follows, we will consider that layers 'A ' and 'B' have the same optical thickness, A A B B n d n d = , and we will introduce the reference wavenumber Ω , defined as / 2 A A n d
Ω = π , which corresponds to the quarter wavelength condition. The ratio / Ω Ω will be called ‘normalized frequency’ and will make the comparison of different structures easier. For example, a normalized frequency / 1 Ω Ω = corresponds to the middle of bandgap region for Bragg crystals, but it is the centre of a wide ‘pseudo-pass band’ for ThM quasicrystals, a region made of densely located shallow resonances, not always of perfect transmission (figure 2). This band exists only for the specific case of equal optical thicknesses of the layers, and it starts to appear from the second generation, for which formulas (8)-(9) are not valid. The above condition makes them applicable for 0 n ≥ . Moreover, as long as the condition is fulfilled, the transmission spectra of all 1D structures will be translationally invariant along frequency axis, ( / 2) ( / ) T T
Ω Ω + = Ω Ω and symmetric with respect to / 1 Ω Ω = . - 6 - Figure 2. ( a ) Linear transmission spectra as a function of normalized frequency ( / Ω Ω ) for ThM quasicrystal with the same set of parameters as in figure 1b. Overlapped diagrams show trifurcation of resonances as generation number increases from 3 to 7 (see vertical axes on the right). ( b ) Same as ( a ), but where band gap regions and their relative strength are emphasized in logarithmic scale.
3. Bistability and multistability in Thue-Morse multilayers
The bistability and multistability properties of Bragg gratings have been studied in depth by a number of authors [24-26]. The knowledge of the field profiles is crucial for the proper understanding of bistability and multistability in ThM quasicrystals. In correspondence to the resonance frequencies these field profiles show a self-similar pattern which strongly resembles that of a ThM sequence of smaller generation number (figure / 0.756 Ω Ω = , and, according to the trifurcation diagram shown in figure 2, this resonance is inherited from the 4th generation. Therefore, the field profile mentioned above has the level of similarity equal to 7 4 3 − = , and it is possible to distinguish = independent blocks. Moreover, these blocks are arranged according to the ThM sequence of 3rd generation. We have found that for any ThM structure with layers of equal optical thickness, the level of self-similarity is maximal for normalized frequencies belonging to the following series: / {0, 1 , 1, 1 , 2, 3 , 3, 3 , ...} σ σ σ σ Ω Ω = − + − + , (15) where 2 21 arctan 0.512( / ) ( / )
A B B A n n n n σ π ⎛ ⎞= − ≈⎜ ⎟⎜ ⎟+⎝ ⎠ . (16) - 7 -
However, the localization strength is relatively small at these frequencies, and thus the nonlinear response cannot be enhanced significantly, because the electric field profile is almost flat along the structure. For this purpose, resonances located near the edges of pseudo band gaps are more suitable. In general, the localization strength is inversely proportional to the level of self-similarity, and when the length of the structure increases, the resonant states gradually change from localized to extended ones [28]. The case of Kerr (or cubic) nonlinearity was considered in this work, so that a nonlinear refractive index was additionally taken into account for each type of layer. This leads to an intensity-dependent self-phase modulation, which is able to shift resonance frequencies [25]. The direction of this shift is determined by the sign of the Kerr coefficient. The most interesting cases evidently occur when the sign of the nonlinearity is such that resonances in the spectrum of transmission bend towards the bandgap regions (figures 3-4). In the following, we always choose positive nonlinear coefficients for the two materials, A n − = ⋅ cm / MW and B n − = ⋅ cm / MW (see also figures 3-4) [29-30]. One of the most serious advantages provided by ThM sequence is that for odd generation numbers the corresponding photonic structures are intrinsically asymmetric, and nonlinearity is capable of making transmission sensitive to the propagation direction [16]. This feature is completely absent in nonlinear Bragg structures, where hysteresis curves are the same regardless of the direction of incidence. Although a similar nonreciprocal behavior can be achieved in the framework of linear optics, the latter typically requires making use of magneto-optical media with externally applied static magnetic fields [31] or chiral media such as cholesteric liquid crystals [32]. Moreover, the polarization Figure 3. ( a ) Electric field profiles at resonance frequency / 0.740 Ω Ω = for linear ThM structure with the same set of parameters as in figure 1b. ( b ) Hysteresis of transmission at / 0.738 Ω Ω = . The nonlinear Kerr coefficients used are A n − = ⋅ cm / MW and B n − = ⋅ cm / MW. The bistability thresholds are different for the forward (blue) and backward (red) incidence. ( c ) Nonlinear transmission spectra for fixed input intensity 10 MW / cm . - 8 - Figure 4. ( a ) Electric field profiles at resonance frequency / 0.810 Ω Ω = for linear ThM structure with the same set of parameters as in figure 1b. ( b ) Hysteresis of transmission at / 0.808 Ω Ω = . The nonlinear Kerr coefficients used are A n − = ⋅ cm / MW and B n − = ⋅ cm / MW. The switching thresholds are different for the forward (blue) and backward (red) incidence. ( c ) Nonlinear transmission spectra for fixed input intensity 10 MW / cm . state is very important for the correct operation of these linear devices, while the design suggested here is free from such complications. The level of self-similarity in the field profiles is also related to the type of hysteresis that can be observed near a particular resonance. In presence of several independent localization centers inside the structure, the interplay between them gives rise to multistability. It is related mainly to the fact that ThM structures of higher generation numbers can be decomposed into those of lower ones. Therefore, to single out a particular hysteresis type in its pure form, it is necessary to check that the corresponding resonance does not occur in previous generations. Two characteristic examples are shown in figures 3-4. The first case (figure 3) corresponds to a resonance with no self-similarity in the field profile and clearly demonstrates bistability of transmission, whereas in the second case (figure 4) two independent localization centers can be distinguished. The hystereses are shown in figure 3b and 4b, the nonlinear transmission spectra are given in figure 3c and 4c. A detailed explanation of this behavior will be given in the next section. Notice that similarly to Bragg structures, the maxima of electric field can be located mostly in the layers of one type, so that in principle it would be sufficient to use only one nonlinear material.
4. Coupled resonator model
In order to provide a qualitative analytical understanding for the nonlinear behavior of ThM quasicrystals, we apply the coupled resonator model, which has been recognized as a useful tool for describing the dynamical properties of single-mode cavities in 2D photonic crystal waveguides [33- - 9 -
20 0 d | |i 1d u v
A A A u vt P + + ⎡ ⎤⎛ ⎞= − − + + +⎢ ⎥⎜ ⎟⎜ ⎟⎢ ⎥⎝ ⎠⎣ ⎦ ω γ κ κ , (17) u uv v r tu u At rv v − +− + ⎛ ⎞ ⎛ ⎞⎛ ⎞ ⎛ ⎞= +⎜ ⎟ ⎜ ⎟⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠⎝ ⎠ ⎝ ⎠κκ . (18) In essence, they represent an extension of the scattering matrix method to the time domain. The first equation (17) describes the temporal evolution of a resonant mode in the vicinity of a resonant frequency ω . The amplitude of this mode A is normalized in such a way that | | A gives the energy inside the cavity. Energy is carried to the cavity from the two ports or channels, which are denoted as U and V , located on the left and the right sides, respectively (figure 5a). The amplitudes of the incident signals are denoted as u + and v + . They are normalized in such a way that | | u + and | | v + give the power flow. The interaction between these signals and the cavity is described by the two coupling coefficients u κ and v κ , respectively. Since the cavity not only accumulates the incident energy, but may also partly absorb it or return it back to the channels, there is a damping constant γ , which characterizes the efficiency of these processes. It can be related to the quality factor of the resonator /(2 ) Q ω γ = . As to the influence of nonlinearity, it results in the shift of the resonant frequency. In equation (17), we introduce a characteristic energy P that allows to take this shift into account [37], which can be either positive or negative depending on the common sign of the Kerr coefficients in the multilayered structure. The second equation (18) describes the output from the cavity. The right hand side of this equation consists of two parts. The first part characterises the background reflection from the cavity as if the amplitude of the mode inside it were zero. In fact, it has exactly the form of a scattering matrix. The background transmission coefficient t is equal for both directions even for nonsymmetric structures, while the reflection coefficients in general have different phases and are denoted as u r and v r . As to the second part of equation (18), it is proportional to the amplitude of the mode, and thus provides the corrections caused by the interaction of incident signals with the cavity. Assuming that the system operates at a fixed frequency, it is possible to derive from equations (17), (18) an explicit expression for the nonlinear transmission spectrum along the forward (backward) propagation direction: ,, out, , , 2out 0in u vu v u v u v u v IT I I I ≡ = + + ηδ , (19) where | | u I u + = ( | | v I v + = ) is the input intensity, | | u I v − = ( | | v I u − = ) is the output intensity, ( ) / = −δ ω ω γ is the normalized detuning from the resonance frequency, and | | / u v = η κκ γ corresponds to the maximum allowed transmission. The characteristic intensities ,0 u v I , and this is a crucial observation for the present work, are in general different for the incidence from the left and from the right , 20 0 ,0 ( / ) | | u v u v I P = γ ω κ . (20) - 10 - Figure 5. ( a ) The scheme describing the quantities used in the equations (17)-(18) for the coupled resonator model. ( b ) The typical form of nonlinear transmission spectrum following from formula (19). In the case of single Dirac-delta-like nonlinear layer in the system, formula (19) coincides with a previously investigated exact analytical solution [38]. However, taking into account a continuously distributed nonlinearity gives rise to a new feature: if the structure is spatially asymmetric, the characteristic intensity can depend on the direction of incidence, as it is evident from equation (20). Since the formula (19) is equivalent to a third-order polynomial in the output intensity ,out u v I , by definition it fails to describe resonances with multistability, but it gives a good approximation for resonances with a strictly bistable response (compare for instance figure 5b with figure 3b-c). In the more complex case, when the field profiles consist of several independent parts, one simply needs to increase the number of equations in the coupled resonator model. There should be exactly one set of equations (17)-(18) for each separable resonance. For example, the case shown in figure 4 can be described by the following system of equations d | |i 1 ( )d a aa A A A u st P + + ⎡ ⎤⎛ ⎞= − − + + +⎢ ⎥⎜ ⎟⎜ ⎟⎢ ⎥⎝ ⎠⎣ ⎦ ω γ κ , (21) d | |i 1 ( )d b bb B B B v st P + − ⎡ ⎤⎛ ⎞= − − + + +⎢ ⎥⎜ ⎟⎜ ⎟⎢ ⎥⎝ ⎠⎣ ⎦ ω γ κ , (22) exp(i ) a a u u As s − +− + ⎛ ⎞ ⎛ ⎞= +⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠ϕ κ , (23) exp(i ) b b v v Bs s − ++ − ⎛ ⎞ ⎛ ⎞= +⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠ϕ κ , (24) - 11 -
Figure 6. ( a ) The scheme describing quantities used in the equations (21)-(24) of the coupled resonator model. ( b , c ) The typical form of nonlinear transmission spectra following from the formula (25), when the dimensionless parameters are set to ω = , 2.89 a γ = , 1.44 b γ = , 120 a P = , 80 b P = . The only difference between figures ( b ) and ( c ) is in the phase mismatch a b ϕ ϕ π + = and a b ϕ ϕ π + = , respectively. The amplitudes of incident signals are set to 2 u + = for the propagation in forward direction (blue lines) and 2 v + = for the backward direction (red lines). The linear transmission spectra are given by gray lines. where A and B give the amplitudes of the two separable resonances, which interact through a virtual port S (figure 6a). These resonances have the same resonance frequency ω , but in general their damping constants , a b γ and characteristic energies , a b P will be different. Since both separable resonances shown in figure 4 are symmetric, the coupling to either of them can be described by single coefficient , , a b a b κ γ= . Moreover, because the common resonant frequency is located in the bandgap (or pseudo band gap) region, the background transmission can be neglected, and the scattering matrix entering equation (18) can be replaced by a reflection coefficient, which is equal for both directions of incidence and has the form , exp(i ) a b ϕ . The system (21)-(24) can be readily solved in the frequency domain. We found the following formula, which describes the nonlinear transmission spectra with multistable behaviour:
1| (1 i )(1 i ) exp[ i( )] | a b a b a b T = + − − − +ξξ ξ ξ ϕ ϕ . (25) Expression (25) contains two parameters, a ξ and b ξ , which can be calculated given the output signal v − (or u − ) by using the following ladder: / b B v − = κ , (26) | | b b b b BP = + ωξ δ γ , (27) - 12 - (1 i i exp[i( )])( / ) b b a b a A v − = − + +ξ ξ ϕ ϕ κ , (28) | | a a a a AP = + ωξ δ γ . (29) The typical nonlinear transmission spectra calculated by using formula (25) are shown in figure 6b-c. The main parameter which determines the shape of nonlinear transmission spectra is the phase mismatch a b + ϕ ϕ . It describes the efficiency of interaction between separable resonances and is responsible for the splitting of initially identical resonant frequencies even in linear case. The only difference between figures 6b and 6c is in the phase mismatch, and the two well-separated resonances shown in the first case become undistinguishable for small intensities in the second case (see gray lines). As the intensity of incident signal increases, this degeneracy disappears giving rise to multistability (see red and blue lines for the forward and backward incidence, respectively). Qualitatively, the nonlinear transmission spectra shown in figure 6c corresponds to that in figure 4c. It is worth noting that if 0 a b + = ϕ ϕ , the formula (25) reduces to
11 ( ) a b T = + + ξ ξ , (30) which means that the two resonances do not interact with each other and the resulting nonlinear transmission spectra reveals only a bistable behaviour.
5. Nonlinear optical diode action in Thue-Morse multilayers
The difference in bistability thresholds can be small, but it creates favorable conditions for unidirectional propagation, so that ThM quasicrystals can be used as all-optical diodes [16]. Similar to electronic circuits, these devices are indispensable, when there is a need to suppress the flow of light in one direction or to avoid problems caused by unwanted reflections. Various types of optical diodes have been proposed and realized, which can work both in linear [31-32] and nonlinear regime [39-42]. The primary figure of merit, which determines the efficiency of this device, is the contrast ratio f b / C T T = between transmission along the forward ( f T ) and the backward ( b T ) directions. C can be very large in ThM structures due to the fact that this device can operate on two different hysteresis branches depending on propagation direction (figure 3b). It is possible to use the intensities of both up and down transitions to achieve strongly nonreciprocal transmission. In the former case, the scheme works in passive mode, but the maximum value of transmission is limited [43]. In the latter case, the transmission can be almost perfect, but the scheme requires an additional short pump signal in order to switch to the higher stable branch of hysteresis [44]. The FDTD method was applied to demonstrate the switching dynamics in time domain (figure 7). To maintain the second order accuracy of the Yee-scheme in case of multilayered structures, nonuniform spatial mesh was used with electric field nodes aligned to the boundaries between layers (see Appendix). Since the structure can accumulate and release energy, the sum of reflected (given by red dotted line in figure 7) and transmitted intensities (blue dashed line) should not necessarily give the incident one (gray solid line). This constraint becomes valid only after all the transient dynamics is over. - 13 - Figure 7.
Numerical simulations of nonlinear ThM structure acting as an optical diode for the forward ( a ) and backward ( b ) incidence at the frequency 0.738 Ω . Left parts of the figures refer to the passive scheme and right ones to the active scheme when an additional short pump signal is applied to facilitate switching. The intensities of incident (solid), reflected (dotted), and transmitted (dashed) signals are given as a function of time. The input intensity was set first to the value, which is sufficiently large for the up switching in the forward direction, but at the same time is smaller than corresponding value for the backward direction. The maximum transmission obtained was f T = with the contrast ratio C = , which is quite close to the theoretical limit of max C = following from (31) in the assumption of passive scheme [43]. In the second stage, the input intensity was set just above the down switching threshold for the forward direction, which ensures transmission of almost f T = in this direction, while in the opposite direction transmission will be strongly suppressed due to the presence of the pseudo band gap. The contrast ratio achieved was 20 C = , but it is not limited in principle. The major drawback of this scheme is that it is necessary to force the system to switch to the upper branch of hysteresis whether by applying an auxiliary pump signal, or by temporarily increasing the input intensity.
6. Conclusions and future work
The nonlinear properties of quasicrystals based on ThM sequence were investigated. It was shown that the interplay between the intrinsic spatial asymmetry of these structures for odd generation numbers and Kerr nonlinearity makes the switching thresholds induced by bistability and multistability sensitive to the propagation direction. The role of self-similarity was emphasized to explain the shape of hysteresis curves observed near a particular resonance, and conditions necessary to achieve highly non-reciprocal propagation were formulated. Coupled resonator model was used as a phenomenological model. FDTD simulations have been used to confirm the results obtained by the nonlinear transfer matrix method, and to show how ThM multilayers can exhibit an effective optical diode action, both in passive and active operation.
Acknowledgments
This work was supported by the German Max Planck Society for the Advancement of Science (MPG). - 14 -
Appendix
Maxwell's equations in differential form ( , )( , ) 1 yz H x tE x tx c t ∂∂ =∂ ∂ , (A.1) ( , ) ( , )1 y z
H x t D x tx c t ∂ ∂=∂ ∂ , (A.2) can be approximated by finite differences
1/ 2 1/ 2 ( 1/ 2) ( 1/ 2) ( 1) ( ) n n n ny y z z c tH i H i E i E ix + −
Δ ⎡ ⎤+ = + + + −⎣ ⎦Δ , (A.3) ( ) ( ) ( 1/ 2) ( 1/ 2) n n n nz z y y c tD i D i H i H ix + + +
Δ ⎡ ⎤= + + − −⎣ ⎦Δ , (A.4) where integer indices n and i serve as temporal and spatial coordinates for the Yee-scheme [45]. Courant stability condition requires c t x Δ ≤ Δ , and it will be assumed that / 2 c t x
Δ = Δ . The spatial grid should be fine enough to capture the smallest wavelength which can be excited in the system. It is usually chosen as min /10 x Δ = λ for uniform grids. However, nonuniform grids can be more advantageous for simulating multilayered structures. To tailor two different uniform grids, it is convenient to use the integral representation of Maxwell's equations 1( , ) ( , ) ( , ) RL xz L z R yx E x t E x t H x t dxc t ∂− = − ∂ ∫ , (A.5)
1( , ) ( , ) ( , ) RL xy L y R zx H x t H x t D x t dxc t ∂− + = ∂ ∫ . (A.6) After integration over the boundary between two layers, the following finite difference equations can be obtained (the symbol used are defined in figure A1)
1/ 2 1/ 2 ( ) ( ) [ ( ) ( )] n n n ny b y b z b z aa b c tH H E E + −
Δ= + − −+β β ρ ρρ ρ , (A.7) ( ) ( ) ( ) ( ) n n n nz b z a y b y aa a b b c tE E H H + + +
Δ ⎡ ⎤= − + − −⎣ ⎦+ρ ρ β βε β ε β . (A.8) It can be proven that these equations are valid up to the second order, if the nodes of electric field are placed on the boundaries between layers, and the steps of the spatial grids inside layers of different types satisfy ( ) ( ) ( ) a a b b x x x
Δ = Δ = Δε ε . This condition can be exactly fulfilled if layers composing the structure are of equal optical thickness. In general case, the second-order error of the scheme with nonuniform mesh is estimated to be ( / ) b b b x m Δ ε δ , (A.9) where m is the number of nodes inside an unmatched layer, and δ is a missing or superfluous thickness that caused the mismatch. When the scheme with uniform mesh is used, the error is proportional to the difference in permittivities and cannot be turned into zero ( )( / 2) b a x − Δε ε . (A.10) - 15 - Figure A1.
The scheme describing quantities used in the tailoring of two uniform grids at the boundary between dielectric media. The constituent relation with instantaneous Kerr-nonlinearity is considered
3L K z z z
D E E = +ε ε . (A.11) According to (A.6) and (A.8), the linear L ε and nonlinear K ε contributions to permittivity should be averaged at the boundaries between layers L,eff L, L, ( ) ( ) / a a b b i x x x = Δ + Δ Δε ε ε , (A.12)
K,eff K, K, ( ) ( ) / a a b b i x x x = Δ + Δ Δε ε ε . (A.13) The equation (A.11) can be written in finite differences as [46] [ 2 ( ) ] /[ 3 ( ) ] n n n n
E D E E − − = + +ε ε ε , (A.14) which is equivalent to the solution of the cubic equation (A.11) by iteration method and is always stable for positive L ε and K ε . In contrast to equation
2L K (3/ 4) | | z z z z
D E E E = +ε ε , (A.15) the usage of constituent relation (A.11) allows taking into account the generation of third harmonic, but attention should be paid to the correct definition of K ε . To separate incident and reflected waves from the structure, the so-called Total-Field Scattered-Field technique is used on the left (and of the right) of the structure
1/ 2 1/ 2F F F ,inc F F ( 1/ 2) ( 1/ 2) ( ) ( ) ( 1) n n n n ny y z z z c tH i H i E i E i E ix + −
Δ ⎡ ⎤= ± − −⎣ ⎦Δ ∓ ∓ ∓ , (A.16) ( ) ( ) ( 1/ 2) ( 1/ 2) ( 1/ 2) n n n n nz z y y y c tD i D i H i H i H ix + + + +
Δ ⎡ ⎤= + + − −⎣ ⎦Δ ∓ ∓ . (A.17) Sources of electromagnetic field on the left (and on the right) are introduced as ,inc S ,inc ( ) ( ) nz z
E i E n c t = Δ , (A.18)
1/ 2,inc S ,inc ( 1/ 2) (( 3/ 2) ) ny z
H i E n c t + = + Δ ∓ ∓ . (A.19) - 16 - The Absorbing Boundary Conditions in case when / 2 c t x
Δ = Δ take the following form on the left (and on the right) boundary of the simulation region
1/ 2 3/ 2B B ( 1/ 2) ( 1/ 2) n ny y
H i H i + − = ± ∓ . (A.20) Smooth envelops for switching on and off monochromatic signals are defined by splines of 6th order, which ensure the continuity of fields up to the second derivative [47] ( / ) for 0 ,( ) 1 for ,[( ) / ] for , u t tU t t Tu T t T t T < <⎧⎪= ≤ ≤ −⎨⎪ − − < <⎩ τ ττ ττ τ (A.21) where T is the total time when the signal is nonzero, τ is the duration of the switching, and ( ) 10 15 6 u x x x x = − + . The intensities of reflected and transmitted waves A can be extracted from the simulation through the following integration (1/ ) cos( ) exp(i ) ( / 2) exp(i ) t P nmt P A m t n t dt A + − = ∫ ω ϕ ω ϕδ , (A.22) where P = π ω is the period of carrying wave. References [1]
Shechtman D, Blech I, Gratias D and Cahn J W 1984 Metallic phase with long-range orientational order and no translational symmetry
Phys. Rev. Lett. Janot C 1994
Quasicrystals (Oxford: Clarendon Press) [3]
Kohmoto M, Sutherland B and Iguchi K 1987 Localization of optics: quasiperiodic media
Phys. Rev. Lett. Chan Y S, Chan C T and Liu Z Y 1998 Photonic band gaps in two dimensional photonic quasicrystals
Phys. Rev. Lett. Zoorob M E, Charlton M D B, Parker G J, Baumberg J J and Netti M C 2000 Complete photonic bandgaps in 12-fold symmetric quasicrystals
Nature
Romero-Vivas J, Chigrin D, Lavrinenko A and Sotomayor Torres C 2005 Resonant add-drop filter based on a photonic quasicrystal
Opt. Express Albuquerque E L and Cottam M G 2004
Polaritons in Periodic and Quasiperiodic Structures (Amsterdam: Elsevier) [8]
Luck J M, Godreche C, Janner A and Janssen T 1993 The nature of the atomic surfaces of quasiperiodic self-similar structures
J. Phys. A Boyd R W 2008
Nonlinear Optics (Boston: Academic Press) [10]
Zhu S, Zhu Y Y and Ming N B 1997 Quasi-phase-matched third-harmonic generation in a quasi-periodic optical superlattice
Science
Cheng Z, Savit R and Merlin R 1988 Structure and electronic properties of Thue-Morse lattices
Phys. Rev. B Liviotti E 1996 A study of the structure factor of Thue-Morse and period-doubling chains by wavelet analysis
J. Phys. Cond. Matt. Luck J M 1989 Cantor spectra and scaling of gap widths in deterministic aperiodic systems
Phys. Rev. B Boriskina S V, Gopinath A and Dal Negro L 2008 Optical gap formation and localization properties of optical modes in deterministic aperiodic photonic structures
Opt. Express Jiang X Y, Zhang Y G, Feng S L, Huang K C, Yi Y H and Joannopoulos J D 2005 Photonic band gaps and localization in the Thue-Morse structures
Appl. Phys. Lett. Biancalana F 2008 All-optical diode action with quasiperiodic photonic crystals
J. Appl. Phys. - 17 - [17]
Werchner M et al.
Opt. Express Steurer W and Deloudi S 2009
Crystallography of Quasicrystals: Concepts, Methods and Structures (Berlin: Springer) [19]
Kolar M 1993 New class of one-dimensional quasicrystals
Phys. Rev. B Negro L D, Yi J H, Nguyen V, Yi Y, Michel J and Kimerling L C 2005 Spectrally enhanced light emission from aperiodic photonic structures
Appl. Phys. Lett. Riklund R, Severin M and Youyan L 1987 The Thue-Morse aperiodic crystal, a link between the Fibonacci quasicrystal and the periodic crystal
Int. J. Mod. Phys. B Yeh P 2005
Optical Waves in Layered Media (New York: Wiley) [23]
Liu N H 1997 Propagation of light waves in Thue-Morse dielectric multilayers
Phys. Rev. B Lidorikis E, Busch K, Li Q M, Chan C T and Soukoulis C M 1998 Wave propagation in linear and nonlinear structures
Physica D
He J and Cada M 1991 Optical bistability in semiconductor periodic structures
IEEE J. Quantum Electron. Chen W and Mills D L 1987 Optical response of nonlinear multilayer structures: bilayers and superlattices
Phys. Rev. B Ryu C S, Oh G Y and Lee M H 1992 Extended and critical wave functions in a Thue-Morse chain
Phys. Rev. B Ostlund S and Pandit R 1984 Renormalization-group analysis of the discrete quasiperiodic Schroedinger equation
Phys. Rev. B Molyneux S, Kar A K, Wherrett B S, Axon T L and Bloor D 1993 Near-resonant refractive nonlinearity in polydiacetylene 9-BCMU thin films
Opt. Lett. Adair R, Chase L L and Payne S A 1989 Nonlinear refractive index of optical crystals
Phys. Rev. B Amemiya T, Shimizu H, Yokoyama M, Hai P N, Tanaka M and Nakano Y 2007 1.54- μ m TM-mode waveguide optical isolator based on the nonreciprocal-loss phenomenon: device design to reduce insertion loss Appl. Opt. Hwang J, Song M H, Park B, Nishimura S, Toyooka T, Wu J W, Takanishi Y, Ishikawa K and Takezoe H 2005 Electro-tunable optical diode based on photonic bandgap liquid-crystal heterojunctions
Nature Mater. Manolatou C, Khan M J, Fan S H, Villeneuve P R, Haus H A and Joannopoulos J D 1999 Coupling of modes analysis of resonant channel add-drop filters
IEEE J. Quantum Electron. Maes B, Bienstman P and Baets R 2005 Switching in coupled nonlinear photonic-crystal resonators
J. Opt. Soc. Am. B Mingaleev S F, Miroshnichenko A E, Kivshar Y S and Busch K 2006 All-optical switching, bistability, and slow-light transmission in photonic crystal waveguide-resonator structures
Phys. Rev. E Fan S, Suh W and Joannopoulos J D 2003 Temporal coupled-mode theory for the Fano resonance in optical resonators
J. Opt. Soc. Am. A Soljacic M, Ibanescu M, Johnson S G, Fink Y and Joannopoulos J D 2002 Optimal bistable switching in nonlinear photonic crystals
Phys. Rev. E Lidorikis E, Busch K, Li Q M, Chan C T and Soukoulis C M 1997 Optical nonlinear response of a single nonlinear dielectric layer sandwiched between two linear dielectric structures
Phys. Rev. B Tocci M D, Bloemer M J, Scalora M, Dowling J P and Bowden C M 1995 Thin-film nonlinear-optical diode
Appl. Phys. Lett. Gallo K, Assanto G, Parameswaran K R and Fejer M M 2001 All-optical diode in a periodically poled lithium niobate waveguide
Appl. Phys. Lett. - 18 - [41] Konorov S O, Sidorov-Biryukov D A, Bugar I, Bloemer M J, Beloglazov V I, Skibina N B, Chorvat D, Scalora M and Zheltikov A M 2004 Experimental demonstration of a photonic-crystal-fiber optical diode
Appl. Phys. B Philip R, Anija M, Yelleswarapu C S and Rao D 2007 Passive all-optical diode using asymmetric nonlinear absorption
Appl. Phys. Lett. Lin X S and Lan S 2005 Unidirectional transmission in asymmetrically confined photonic crystal defects with Kerr nonlinearity
Chin. Phys. Lett. Lin X S, Yan J H, Wu L J and Lan S 2008 High transmission contrast for single resonator based all-optical diodes with pump-assisting
Opt. Express Taflove A and Hagness S 2005
Computational Electrodynamics: The Finite-Difference Time-Domain Method (Boston: Artech House) [46]
Sullivan D M 1995 Nonlinear FDTD formulations using Z-transforms
IEEE Trans. Microw. Theory Tech. Ziolkowski R W and Heyman E 2001 Wave propagation in media having negative permittivity and permeability
Phys. Rev. E64