# Photonic Wannier-Stark Ladder from Coupled Electromagnetic Cavities

Authors: Shahzad Anwar; Sucheng Li; Weixin Lu; Bo Hou
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

We have investigated the photonic Wannier-Stark ladder in the system of coupled electromagnetic cavities, which consists of a stack of metallic plates structured with subwavelength apertures and where the tilted potential effect is mimicked by imposing the gradient variation of refractive index. Making an analogy to its quantum counterpart and assuming the translational property of its solutions, we have shown the photonic ladder has the eigenenergies, that is, frequencies, in a geometrical series. Within the approximation of small gradient, the ladder states manifest the equidistant frequency spacing in the spectrum. By both analytical derivation and numerical simulation, we have illustrated the geometrically progressed energies of the photonic Wannier-Stark ladder.


---

## Body

## 1. Introduction

A lot of similar physical properties have been revealed between electrons in a solid and photons in a periodic structure [1]. Besides band gap, another fascinating example is the concepts of Wannier-Stark ladder (WSL) and Bloch oscillation (BO) [2]. When a solid is biased by a static electric field, the electron eigenstates do not belong to the extended Bloch states because the linear potential imposed by the static electric field breaks the translational symmetry of the periodical crystal potential but becomes the WSL which is a family of localized wavefunctions with equidistant energies in the energy spectrum. The Bloch oscillation, in which the electrons driven by the DC bias oscillate periodically in semiclassical picture, is just the time domain counterpart of the quantum-mechanical WSL. Due to the potential applications toward high speed electronics, both WSL and BO have been theoretically investigated and experimentally observed in the semiconductor superlattices [3]. In contrast, WSL and BO in natural crystals are difficult to observe because the coherence of electron wave packets is very fragile to various decoherence processes, for example, disorder, in the crystals.Recently, the two concepts advance to electromagnetic (EM) waves in an analogy way, driven by the physical similarity of wave systems and the phenomenological robustness of classical waves [4]. For EM waves, the gradient refractive index or thickness of the dielectric layers introduced to a photonic structure plays the same role as the DC field in electronic cases and gives rise to both the photonic WSL in the frequency domain and the photonic BO in the time domain [5–7]. Besides the conventional photonic crystals [8–11], the microcavity systems from metamaterial with negative index have been shown to display the photonic BO [12–16]. In addition, the BO can be demonstrated equivalently inside the transversely coupled waveguide arrays undergoing the gradient modulation along the transverse direction [17, 18] and manifests the rich physics of wave propagation, such as quasi and fractional versions [19, 20]. More recently, the WSL has also been observed in acoustic/elastic systems [21–24].In this work, we have investigated the photonic WSL in the system of coupled EM cavities, where the biased potential effect is mimicked by imposing the gradient variation of refractive index. A numerical model is constructed, using a stack of metallic plates patterned with subwavelength apertures. By calculating the transmission, we have observed that the transmission peaks of the ladder states have the equidistant frequency spacing in the spectrum, that is, their frequencies (eigenenergies) showing an arithmetical series, in the case of small gradient. Specifically, the frequencies of the peaks evolve into a geometrical series in the case of large gradient, which points out the generalized photonic WSL with geometrically progressed energies.


## 2. Materials and Methods

We start with the quantum WSL problem. Consider a particle of massm and charge q moving in the simplified crystal potential of spatial periodicity d (Kronig-Penney model [25]) tilted by a static electric field ξ, as depicted in Figure 1(a). For the quantum system, the Schrodinger equation has the following expression:(1)-ħ22md2ψxdx2+Vxψx=Eψx,where ψ(x) is the wavefunction, E is the eigenenergy, and the position-dependent potential has the property V(x)=V(x+d)+qξd. When an electron moves in the tilted periodic potential, it will localize somewhere due to the translational symmetry breaking. Assume such a state ψ(x) localized near well 0 and with energy E0, and then the above equation reads(2)-ħ22md2ψxdx2+Vxψx=E0ψx.Change the variable x to x-jd (j being an integer), and we obtain(3)-ħ22md2ψx-jddx2+Vx-jdψx-jd=E0ψx-jd;that is,(4)-ħ22md2ψx-jddx2+Vx-jd-qξjdψx-jd=E0-qξjdψx-jd,or(5)-ħ22md2ψx-jddx2+Vxψx-jd=Ejψx-jd,withEj=E0-qξjd.Therefore, a new solution ψ(x-jd) is localized near the well j and has the energy Ej. It is seen that a family of such states can be obtained with the consistent wavefunctions ψ, except that their locations translate with respect to each other. The family of states is just the WSL where their eigenenergies compose an arithmetical progression, Ej=E0-qξjd, with common difference qξd.Figure 1
(a) The Kronig-Penney potential (coupled quantum wells) titled by a static electric field, whered is the spatial periodicity between the neighboring wells and the red line depicts schematically the quantum WSL state localized in the different well. (b) The one-dimensional coupled EM cavities, that is, “photonic wells,” with increasing the refractive indexes of the dielectric layers, α, β, γ, and δ, along the x-direction.
(a)
(b)The photonic WSL has been proved in one-dimensional coupled EM cavities [5, 6], where a ramp-up of the refractive index along the stacking direction is used to mimic the effect of the electric field ξ in the quantum case, as illustrated in Figure 1(b). Here, we reformulate the EM problem in the same style as solving the quantum system. For explicitness and simplification, the individual cavity is assumed including four dielectric layers, and the solution also holds for the cavity with more or less layers. The Helmholtz wave equation for the photonic system is written as(6)d2dx2Fx+n2xω2c2Fx=0,where the refractive index profile (7)nx=niα,x∈id,id+Lα,niβ,x∈id+Lα,id+Lα+Lβ,niγ,x∈id+Lα+Lβ,id+Lα+Lβ+Lγ,niδ,x∈id+Lα+Lβ+Lγ,i+1d,with the integer i specifying the cavity cell and the Greek alphabet specifying the layer, and the cavity length d=Lα+Lβ+Lγ+Lδ with Lα, Lβ, Lγ, and Lδ being the thicknesses of the four layers, α, β, γ, and δ, respectively. Following n(x), the photonic state, F(x), is supposed to be sectioned. Similar to the quantum case, we hypothesize a solution, F(x), which is centralized somewhere, to say near cell 0, and has the photon energy ħω0 and the field function:(8)Fx=Fiαx,x∈id,id+Lα,Fiβx,x∈id+Lα,id+Lα+Lβ,Fiγx,x∈id+Lα+Lβ,id+Lα+Lβ+Lγ,Fiδx,x∈id+Lα+Lβ+Lγ,i+1d,where the cell number i may range within the integers [-M,M], depending on the localizing degree. In fact, the positive integer M characterizes the localization scale of the solution (M→0 representing a state localized inside only one cell and M→∞ meaning the weak localization). Without losing generality, the solution in the region x∈[id+Lα,id+Lα+Lβ) satisfies (9)d2dx2Fiβx+niβ2ω02c2Fiβx=0.Changing the variable, we obtain(10)d2dx2Fiβx-jd+niβ2ω02c2Fiβx-jd=0,or (11)d2dx2Fiβx-jd+nj+iβ2ωj2c2Fiβx-jd=0,where nj+iβωj=niβω0 and x∈[(j+i)d+Lα,(j+i)d+Lα+Lβ). With applying to the other layers and cells, we can see that the new field function, F(x-jd), its segmented expressions being Fiα,β,γ,δ(x-jd), actually represents another state which centralizes near the cell j and has the energy ħωj, where ωj=niω0/nj+i. For the strongly localized photonic state (M→0), for example, the one which allows the tight binding description of the coupling effect between the EM cavities, the energy relation niω0=nj+iωj can be simplified as n0ω0=njωj. This relation can also be justified directly from (6), if the photonic WSL takes a series of solutions which are spatially translational to each other, like the quantum case.It can be proved that a geometric series,ni=n0qi and ωj=ω0/qj (q being the common ratio), meets exactly the relation niω0=nj+iωj or n0ω0=njωj. Particularly, in the case of q=1+η with the approximation η≪1, ni=n0(1+η)i≈n0+in0η and ωj=ω0(1+η)-j≈ω0-jω0η, where ω0η is identified as Bloch oscillation frequency ωB which manifests the equidistant frequency spacing in spectra.


## 3. Results and Discussion

In order to illustrate the photonic WSL states in a specific system, we employ a metallic plate patterned with narrow H-fractal slits as EM cavity. The resonant state supported by the exotic fractal pattern has been found to be responsible for transmission enhancement of EM wave through such metallic plates [26]. When a number of the metallic plates with H-fractal slits are cascaded, their resonant modes will couple together; the transmission bands would be developed in a way analog to the formation of miniband in the superlattice and furthermore comply with the photonic tight binding description [27]. Note that our fractal slits should not be considered as a sole candidate, while all other shaped apertures supporting the local resonance are applicable, too.The coupled cavity system is 8 stacked metallic plates on which a periodic array of fractal slits was cut with lattice constant smaller than the relevant incident wavelength to avoid the grating effect, illustrated schematically in Figure2(a). The unit cell of the array includes a five-level H-fractal slit structure where the longest slit is 10 mm and the width of each slit is 0.8 mm [26], and the plates have a thickness of 3 mm and a spacing of 9 mm. Numerically, we employ a commercial finite-difference time-domain (FDTD) software package to investigate the microwave spectra [28]. The microwave polarization was such that the electric field is perpendicular to the shortest slits of the fractal pattern. In the simulation, we assume periodic boundary condition in the transverse direction and perfect conductor approximation for metal.Figure 2
(a) The 8 stacked metallic plates on which a periodic array of fractal slits was generated. In the schematic picture, the plate has 2 × 2 unit cells, each of which is a five-level H-fractal slit. (b) The 8 plates are embedded, respectively, inside 8 dielectrics, which have the refractive index,nj=1+dnj, where j takes the denoted numbers. The index gradient is along the propagation direction, as pointed by the arrow. The dash lines depict a single EM cavity cell composed of a structured plate and a homogeneous dielectric. (c) The simulated transmissions at normal incidence of microwave, where three cases, no, small, and large gradient of refractive index, are considered.
(a)
(b)
(c)To impose a tilted potential effect, we embed the plates in the dielectrics, which have the geometrically increasing refractive index,nj=1+dnj, where the integer j denotes the plate number from 0 to 7 and the varying quantity dn tunes the magnitude of the common ratio of the geometrical series. Each structured plate is embedded inside the dielectric in a symmetrical way, as illustrated in Figure 2(b).Whendn=0, which is referred to as the periodic case because there is no refractive index gradient, the simulated result shows that a passband is formed between 2.5  and 3.1 GHz in the same physics as the miniband in the semiconductor superlattice, seeing the upper curve in Figure 2(c). The peaks originate from the coupling of the resonant modes in the plates at different phase delay, analogous to the development of energy band (or miniband) in crystal (or superlattice) within the tight binding description [25, 27]. It is seen that the peaks are more crowded at the two band edges than in the middle region, which is the typical density-of-state characteristics of the miniband.After a small variation,dn=0.01, is introduced, the magnitudes and positions of the peaks change drastically, seeing the middle curve in Figure 2(c). The underlying resonance states, responsible for transmission, become localized somewhere in the stacking direction. For example, the peak at the lower band edge has degraded into a spectral shoulder (~2.45 GHz) with low transmission (~10%). This is because the fractal slits which support the peak are being located at the far end plate of the stacking system, where the refractive index is highest, and the other plates with the higher resonance frequencies behave like the potential barrier. Consequently, the incident wave has to tunnel through them to excite the eigenmode of the farthest one, similar to the Zener tunneling [23, 24]. Likewise, the situation also occurs at the upper band edge, the tiny peak at 3.0 GHz. Note that the peaks’ spacing in frequency begin to appear even.While a large variation,dn=0.1, is considered, an obvious downshift of the miniband region is seen in the spectrum, as plotted in Figure 2(c), originating from the overall increasing in the refractive index. At the same time, the transmission magnitude is suppressed significantly, as the titled potential is much steeper and the states have become more localized, so that only five peaks are identified, seeing Figure 3 for their frequencies. The decrease of the number of spectra peaks is a consequence of larger dn or steeper titled potential, because a few of WSL states (e.g., the farthest one) cannot be efficiently excited any more via the tunneling. Likewise, the spectra spacing of the peaks become wider with increasing dn, due to the larger eigenenergy differences between neighboring WSL states under steeper potential.Figure 3
(a) The peak frequencies in the three cases, periodic, smalldn, and large dn. The vertical axis is in linear scale. (b) The same as (a) but with the vertical axis being logarithmic scale. The horizontal axis in (a) and (b) represents the peak number. The numbers, 3, 4, 5, 6, and 7, of the five peaks in the case of large dn are arbitrarily assigned, being not relevant to the line shape recognition. To see the linear line shape, a straight line connecting two end points is drawn, respectively, in (a) and (b).
(a)
(b)We have labeled the frequency values of all identified peaks from Figure2(c) and plotted them in Figure 3. With the linearly scaled frequency axis, as shown in Figure 3(a), the frequency dots in the periodic case display a curved trend, consistent with the tight binding description of the miniband [27]. In contrast, the line shapes of the latter two cases can be approximately regarded as being linear. However, if the axis is switched to be logarithmic scale, seeing Figure 3(b), the frequency dots of the gradient cases are more close to a linear curve, especially for the large dn case, where the dots of peaks 4, 5, and 6 fall completely on the straight line connecting two end points (the dots of peaks 3 and 7). The linear line shape in logarithmic scale reveals that these frequency values are actually in a geometric sequence, which is consistent with the analytical derivation.From the experimental perspective, the critical requirement for observing the proposed WSL is to fabricate the dielectrics with designed refractive indexnj=1+dnj. The feasible technique route of well-controlling refractive index or dielectric constant is to employ composite material whose permittivity can be fine-tuned by the component concentration according to effective medium theory [29–31]. The composite suitable to our model may be the plastic foam where the porosity (air concentration) can be engineered to realize the desired dielectrics. We will conduct the experimental investigation along with this technique consideration in the future.


## 4. Conclusions

In conclusion, by making an analogy to the quantum counterpart and assuming the translational properties of the solutionF(x), we have shown the photonic WSL in the system of coupled EM cavities, where the tilted potential effect is mimicked by imposing the gradient variation of refractive index along the propagation direction, may have the eigenenergies (frequencies) in geometrical progression. Within the approximation of small gradient, the WSL states have the equidistant frequency spacing in the spectrum. The numerical simulations of the model system, the stack of structured plates, have illustrated the generalized photonic WSL state with geometrically progressed energies.



---

*Source: 360359-2015-10-04.xml*