# Klein Tunneling of Light in Fiber Bragg Gratings

Authors: Stefano Longhi
Journal: Physics Research International
 (2010)

---

## Abstract

A photonic analogue of Klein tunneling (KT), that is, of the exotic property of relativistic electrons to
pass a large repulsive and sharp potential step, is proposed for pulse propagation in a nonuniform
fiber Bragg grating with an embedded chirped region. KT can be simply observed as the opening of
a transmission window inside the grating stop band, provided that the impressed chirp is realized
over a length of the order of the analogue of the Compton wavelength.


---

## Body

## 1. Introduction

A remarkable prediction of the Dirac equation is that a below-barrier electron can pass a large repulsive andsharp potential step without the exponential damping expected for a nonrelativistic particle. Such a transparency effect, originally predicted by Klein [1] and referred to as Klein tunneling (KT), arises from the existence of negative-energy solutions of the Dirac equation and requires a potential step height ΔV of the order of twice the rest energy mc2 of the electron [2]. Relativistic tunneling across asmooth potential step, which describes the more physical situation of a constant electric field E in a finite region of space of length l, was subsequently studied by Sauter [3]. Sauter showed that to observe barrier transparency, the potential increase ΔV≃eEl should occur over a distance l of the order or smaller than the Compton wavelength λC=ℏ/(mc), the transmission probability rapidly decaying toward zero for a smoother potential increase [2–4]. The required field corresponds to the critical field for e+e- pair production in vacuum, and its value is extremely strong making the observation of relativistic KT for electrons very challenging. Therefore, growing efforts have been devoted to find experimentally accessible systems to investigate analogs of relativistic KT [5]. Recently, great interest has suscitated the proposal [6] and first experimental evidences [7, 8] of KT for non-relativistic electrons in graphene, which behave like massless Dirac fermions. On the other hand, optics has offered on many occasions a test bed to investigate the dynamical aspects embodied in a wide variety of coherent quantum phenomena (see, for instance, [9] and references therein). In optics, several proposals of KT analogs have been suggested as well, including light propagation in deformed honeycomb photonic lattices [10] whose band structure is similar to the one of graphene [11, 12], light refraction at the interface between positive-index and negative-index media [13], spatial light propagation in binary waveguide arrays [14], and stationary light pulses in an atomic ensemble with electromagnetically induced transparency [15]. The experimental implementations of such schemes, however, might be a nontrivial matter, and an experimental observation of KT for photons is still lacking. On the other hand, multilayer and Bragg dielectric structures, such as fiber Bragg gratings (FBGs), are rather simple photonic devices with flexible design that have been successfully demonstrated to provide an accessible laboratory tool to investigate photonic analogues of non-relativistic tunneling phenomena [16–18]. Here it is shown that an optical analogue of KT can be achieved in a nonuniform FBG composed by two periodic sections linked by a chirped section which mimics an external potential step in the Dirac equation. Such a FBG-based system might be considered the simplest system proposed so far in order to observe Klein tunneling in any optical system.


## 2. Quantum-Optical Analogy

The starting point of our analysis if provided by a standard model of light propagation in an FBG with a longitudinal refractive indexn(z′)=n0+Δnm(z′)cos[2πz′/Λ+2ϕ(z′)], where n0 is the effective mode index in absence of the grating, Δn≪n0 is the peak index change of the grating, Λ is the nominal period of the grating defining the reference frequency ωB=πc/(Λn0) of Bragg scattering, c is the speed of light in vacuum, and m(z'), 2ϕ(z') describe the slow variation, as compared to the scale of Λ, of normalized amplitude and phase, respectively, of the index modulation. Note that the local spatial frequency of the grating is k(z')=2π/Λ+2(dϕ/dz'), so that the local chirp rate is C=dk/dz'=2(d2ϕ/dz'2). The periodic index modulation leads to Bragg scattering between two counterpropagating waves at frequencies close to ωB. By letting E(z',t)=φ1(z',t)exp[-iωBt+ikBz'+iϕ(z')]+φ2(z',t)exp[-iωBt-ikBz'-iϕ(z')]+c.c. for the electric field in the fiber, where kB=π/Λ, the envelopes φ1 and φ2 of counterpropagating waves satisfy the coupled-mode equations [19](1)i[∂z'+(1vg)∂t]φ1=(dϕdz′)φ1-κ(z')φ2,i[-∂z'+(1vg)∂t]φ2=(dϕdz′)φ2-κ(z')φ1,
where κ(z')≡[kBm(z')Δn]/(2n0) and vg~c/n0 is the group velocity at the Bragg frequency. The analogy between pulse propagation in the FBG and the Dirac equation in presence of an electrostatic field is at best captured by introducing the dimensionless variables z=z'/Z and τ=t/T, with characteristic spatial and time scales Z=2n0/(kBΔn) and T=Z/vg, and the new envelopes ψ1,2(z')=[φ1(z')∓φ2(z')]/2. In this way, (1) can be cast in the Dirac form(2)i∂τψ=-iσ1∂zψ+mσ3ψ+V(z)ψ
for the spinor wave function ψ=(ψ1,ψ2)T, where V(z)=(dϕ/dz) and σ1,3 are the Pauli matrices, defined by(3)σ1=(0110),σ3=(100-1).In its present form, (2) is formally analogous to the one-dimensional Dirac equation with ℏ=c=1 in presence of an external electrostatic potential V(z), m playing the role of a dimensionless (and generally space-dependent) rest mass (see, for instance, [2, 4]). As is wellknown, a nonvanishing mass m is responsible for the existence of a forbidden energy region, which separates the positive- and negative-energy branches of the massive Dirac equation. The optical analogue of the forbidden energy region is precisely the photonic stop band of the periodic grating. As the refractive index modulation of the grating, that is, the mass term m in the Dirac equation (2), is decreased, the stop band region shrinks and the limit of a massless Dirac equation (similar to the one describing the dynamics of electrons in graphene near a Dirac point) is attained. The additional external potential V in (2), related to the chirp of the grating according to V(z)=(dϕ/dz), changes the local position of the forbidden energy region. Therefore, pulse propagation in an FBG with a suitably designed chirp profile can be used to mimic the relativistic tunneling of a wave packet in a potential step V(z). It should be noted that, as compared to other photonic analogues of KT recently proposed in [10, 14] and based on spatial light propagation in periodic photonic structures, the phenomenon of KT occurring in FBGs and discussed in the following section involves the temporal (rather than the spatial) light dynamics and can be therefore simply investigated in the frequency domain by spectrally resolved transmission measurements.


## 3. Klein Tunneling

To realize the analogue of KT, let us first assume that the optical pulse propagates in a region of the grating wherem(z) is uniform and equal to one, and let us assume a chirp profile that mimics a step potential with an increase from V=0 to V=ΔV which occurs over a length l (see Figure 1). Since for the Dirac equation (2) written in dimensionless units the Compton length is λC=1 and the rest energy is mc2=1; according to Sauter's analysis KT is expected to be observable for l smaller than ~1 and for a potential height ΔV larger than 2 [2–4]. The process of KT and tunneling inhibition for a smooth potential step can be simply explained by a graphical analysis of the space-energy diagrams (z,Ω) of the one-dimensional Dirac equation [2], which are shown in Figures 1(a) and 1(b) for a sharp and for a smooth potential step, respectively. For the sake of clearness, in the figures the potential V(z) has been chosen to yield a nonvanishing and constant chirp rate over a length l; different forms for the potential step, such as the profile V(z)=(ΔV/2)[1+tanh(z/l)] considered in the seminal work by Sauter [3], can be assumed as well without changing the main results.Energy diagrams of the Dirac equation (2) for (a) a sharp and (b) a smooth potential step V(z) of height ΔV. The shaded regions are the forbidden energies that separate the electron and positron states, the dotted horizontal line is the energy Ω0 of the incoming wave packet, and the dashed curve is the shape of the potential step V(z). (c) Schematic of the grating structure that realizes the optical analogue of relativistic tunneling across a potential step. The grating comprises five sections, denoted by roman numbers (from I to V) that are defined by the amplitude m(z) and phase gradient V(z)=(dϕ/dz) grating profiles.
(a)(b)(c)The space-energy diagrams of Figures1(a) and 1(b) schematically show the behavior of the energy spectrum of (2) versus z, which is composed by two branches—the electron and positron energy branches of the Dirac equation—separated by a gap of width 2m(z) and centered along the curve Ω=V(z). The gap regions are visualized in the diagrams by the shaded areas. A wave packet (optical pulse) in the electron branch with an initial mean energy Ω0 (1<Ω0<ΔV-1) coming from z→-∞ tunnels into the z>0 region after crossing a forbidden energy region, indicated by the bold segment AB in Figure 1(b). The forbidden gap region AB vanishes for a sharp potential step (l=0, see Figure 1(a)). According to Sauter's analysis [2, 3], the tunneling probability is appreciable provided that l is smaller than ~1. In the FBG context, the energy diagrams of Figure 1 are equivalent to the band-reflection diagrams introduced by Poladian for a graphical analysis of nonuniform gratings [20], where the energy Ω represents the frequency detuning of the incoming wave from the Bragg frequency ωB. The Sauter's condition l<~1 for KT can be derived following the analysis of [20] by computation of the transmittance of the effective grating associated to the evanescent region AB shown in Figure 1(b) (see Section 5.1 of [20]). In the previous discussion, we assumed m(z)=1, however for a grating with finite spatial extent one has m(z)→0 as z→±∞. To inject and to eject the optical pulse into the m(z)=1 grating region around z=0, an input and an output apodization sections can be introduced, which adiabatically convert the input and output wave packets from the m(z)=0 regions into the m(z)=1 grating region (see Figure 1). Therefore, the general structure of the FBG that realizes a photonic analogue of relativistic tunneling across a potential step consists of five sections, as shown in Figure 1(c): two boundary apodization sections (regions I and V), and two uniform sections (regions II and IV) separated by a central chirped section of length ~l (region III). In Figures 2(a)–2(c) typical examples of pulse tunneling across the potential step V(z)=(ΔV/2)[1+tanh(z/l)] are presented, showing KT for a sharp potential step (Figure 2(a)) and inhibition of tunneling as the step gets smooth (Figures 2(b) and 2(c)). The figures depict the temporal evolution of |ψ1|2+|ψ2|2=|φ1|2+|φ2|2—which is proportional to the field intensity averaged in time over a few optical cycles and in space over a few wavelengths—as obtained by numerical analysis of (1) for a grating length of z=160 with a quarter-cosine apodization profile (see Figure 2(d)), ΔV=6, and for a few values of l. A forward-propagating Gaussian pulse φ1 of mean energy Ω0=2 coming from z→-∞ and of duration (FWHM in intensity) τp=5 has been assumed as an initial condition. For typical parameter values n0=1.45, Δn=3.3×10-4, and λB≡2n0Λ=1560 nm, which apply to FBGs used in optical communications, the spatial and temporal scales in Figure 2 are Z≃1.5 mm and T≃7.3 ps, respectively. Hence, in physical units the grating length is L≃24 cm whereas the optical analogue of the Compton length is λC=Z≃1.5 mm. Such nonuniform FBG structures should be realizable with current FBG technology based on UV continuous laser writing [21]. It should be finally noticed that, as in an experiment the tracing of pulse evolution in the grating (Figure 2) can be a nontrivial task, the signatures of KT can be simply obtained from standard spectral transmission measurements of the grating. In fact, for a given value of ΔV>2 and according to the band diagram of Figure 1(a), in the KT regime a transmission window at 1<Ω<ΔV-1 (the segment AB in Figure 1(a)), embedded into the two gaps ΔV-1<Ω<ΔV+1 and -1<Ω<1 (the segments BC and AD in Figure 1(a)) should be observed in the transmission spectrum, the suppression of KT for a smooth potential corresponding to the lowering of such a transmission window. This is clearly shown in Figure 3, where the spectral transmittance of the FBGs corresponding to the simulations of Figures 2(a), 2(b), and 2(c) are depicted. Note that, as the length l of the chirped region is increased (from Figures 3(a) to 3(c)), the transmission window embedded in the two adjacent gaps disappears, which is the signature of KT inhibition.(a)–(c) Pulse propagation in a FBG with an chirp profileV(z)=(ΔV/2)[1+tanh(z/l)] for ΔV=6 and for (a) l=0.1, (b) l=1.5, and (c) l=5. (d) Profiles of V=(dϕ/dz) (upper plot) and of grating amplitude m(z) (apodization profile, lower plot). The potential V is shown for l=5, corresponding to the simulation of Figure 2(c).
(a)(b)(c)(d)Numerically computed spectral transmittance of FBGs used in numerical simulations of Figure2 for n0=1.45, Δn=3.3×10-4, and λB=1560 nm, corresponding to a grating length L≃24 cm, and for increasing values of l: (a) l=0.1, (b) l=1.5, and (c) l=5.
(a)(b)(c)


## 4. Conclusions

In conclusion, a photonic analogue of Klein tunneling based on pulse propagation in nonuniform fiber Bragg gratings has been proposed. As compared to other photonic analogues of KT recently proposed in [10, 14] and based on spatial light propagation in periodic photonic structures, the phenomenon of KT in FBGs suggested in this work can be simply observed in the frequency domain as the opening of a transmission window inside the grating stop band, provided that the impressed chirp is realized over a length of the order of the analogue of the Compton wavelength. Such a FBG-based system might be thus considered to be the simplest optical analogue proposed so far to observe KT.



---

*Source: 645106-2010-08-31.xml*