Interband resonant transitions in
two-dimensional hexagonal lattices:
Rabi oscillations, Zener tunnelling, and
tunnelling of phase dislocations

doi:10.1364/OE.16.014076

« Show journal navigation

Interband resonant transitions in two-dimensional hexagonal lattices: Rabi oscillations, Zener tunnelling, and tunnelling of phase dislocations

Valery S. Shchesnovich, Anton S. Desyatnikov, and Yuri S. Kivshar »View Author Affiliations

Optics Express, Vol. 16, Issue 18, pp. 14076-14094 (2008)
http://dx.doi.org/10.1364/OE.16.014076

View Full Text Article

Acrobat PDF (1390 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (48)
Cited By
Figures (8)
Multimedia (6)
Metrics

Abstract

We study, analytically and numerically, the dynamics of interband transitions in two-dimensional hexagonal periodic photonic lattices. We develop an analytical approach employing the Bragg resonances of different types and derive the effective multi-level models of the Landau-Zener-Majorana type. For two-dimensional periodic potentials without a tilt, we demonstrate the possibility of the Rabi oscillations between the resonant Fourier amplitudes. In a biased lattice, i.e., for a two-dimensional periodic potential with an additional linear tilt, we identify three basic types of the interband transitions or Zener tunnelling. First, this is a quasi-one-dimensional tunnelling that involves only two Bloch bands and occurs when the Bloch index crosses the Bragg planes away from one of the high-symmetry points. In contrast, at the high-symmetry points (i.e., at the M and Γ points), the Zener tunnelling is essentially two-dimensional, and it involves either three or six Bloch bands being described by the corresponding multi-level Landau-Zener-Majorana systems. We verify our analytical results by numerical simulations and observe an excellent agreement. Finally, we show that phase dislocations, or optical vortices, can tunnel between the spectral bands preserving their topological charge. Our theory describes the propagation of light beams in fabricated or optically-induced two-dimensional photonic lattices, but it can also be applied to the physics of cold atoms and Bose-Einstein condensates tunnelling in tilted two-dimensional optical potentials and other types of resonant wave propagation in periodic media.

1. Introduction

Electrons in crystalline solids or semiconductor superlattices, cold atoms in optical lattices, light beams in photonic crystals or waveguide arrays have the energies confined to the Bloch bands separated by band gaps. The system response to a weak linear potential (i.e., a weak constant tilt) is in the form of Bloch oscillations [1

F. Bloch, “Ü ber die Quantenmechanik der Elektronen in Kristallgittern,” Z. Phys. 52, 555–600 (1928).

] and the interband transitions known as Zener tunnelling [2

C. Zener, “Non-Adiabatic Crossing of Energy Levels,” Proc. R. Soc. London A 137, 696–702 (1932).

] (see also [3

G. H. Wannier, “Wave Functions and Effective Hamiltonian for Bloch Electrons in an Electric Field,” Phys. Rev. 117, 432–439 (1960). [CrossRef]

]). The case of an avoided crossing of two Bloch bands is described by a two-level model due to Zener [2

C. Zener, “Non-Adiabatic Crossing of Energy Levels,” Proc. R. Soc. London A 137, 696–702 (1932).

], Landau [4

L. D. Landau, “On the theory of transfer of energy at collisions II,” Phys. Z. Sowjetunion 2, 46 (1932).

] and Majorana [5

E. Majorana, “Orientated atoms in a variable magnetic field,” Nuovo Cimento 9, 43–50 (1932). [CrossRef]

]. Examples of Zener tunnelling include electrical breakdown in Zener diodes [6

L. Esaki, “Long journey into tunneling,” Rev. Mod. Phys. 46, 237–244 (1974). [CrossRef]

], electrical conduction in nanotubes [7

B. Bourlon, D. C. Glattli, B. Placais, J. M. Berroir, C. Miko, L. Forro, and A. Bachtold, “Geometrical Dependence of High-Bias Current in Multiwalled Carbon Nanotubes,” Phys. Rev. Lett. 92, 026804 (2004). [CrossRef] [PubMed]

] and super lattices [8

A. Sibille, J. F. Palmier, and F. Laruelle, “Zener Interminiband Resonant Breakdown in Superlattices,” Phys. Rev. Lett. 80, 4506–4509 (1998). [CrossRef]

], tunnelling of the Bose-Einstein condensates (BECs) in optical lattices [9

M. Jona-Lasinio, O. Morsch, M. Cristiani, N. Malossi, J. H. Müller, E. Courtade, M. Anderlini, and E. Arimondo, “Asymmetric Landau-Zener Tunneling in a Periodic Potential,” Phys. Rev. Lett. 91, 230406 (2003). [CrossRef] [PubMed]

, 10

M. Jona-Lasinio, O. Morsch, M. Cristiani, E. Arimondo, and C. Menotti, “Nonlinear Effects for BoseEinstein Condensates in Optical Lattices,” Laser Phys. 15, 1180–1188 (2005).

], and an optical analog of tunnelling in arrays of optical waveguides and photonic crystals [11

G. Malpuech, A. Kavokin, G. Panzarini, and A. Di Carlo, “Theory of photon Bloch oscillations in photonic crystals,” Phys. Rev. B 63, 035108 (2001).

, 12

V. Agarwal, J. A. del Río., G. Malpuech, M. Zamfirescu, A. Kavokin, D. Coquillat, M. Scalbert, M. Vladimirova, and B. Gil, “Photon Bloch Oscillations in Porous Silicon Optical Superlattices,” Phys. Rev. Lett . 92, 097401 (2004). [CrossRef] [PubMed]

, 13

H. Trompeter, T. Pertsch, F. Lederer, D. Michaelis, U. Streppel, A. Bräuer, and U. Peschel, “Visual Observation of Zener Tunneling,” Phys. Rev. Lett. 96, 023901 (2006). [CrossRef] [PubMed]

The above studies, however, consider only the one-dimensional case. In the recent experiment [14

H. Trompeter, W. Krolikowski, D. N. Neshev, A. S. Desyatnikov, A. A. Sukhorukov, Yu. S. Kivshar, T. Pertsch, U. Peschel, and F. Lederer, “Bloch Oscillations and Zener Tunneling in Two-Dimensional Photonic Lattices,” Phys. Rev. Lett. 96, 053903 (2006). [CrossRef] [PubMed]

], the interband transitions have been observed for the first time in a twodimensional periodic structure of square symmetry (an optically-induced photonic lattice). Multi-dimensional optical lattices are also routinely used for trapping of ultracold atoms [15

M. B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, “Bloch Oscillations of Atoms in an Optical Potential,” Phys. Rev. Lett. 76, 4508–4511 (1996). [CrossRef] [PubMed]

, 16

E. Peik, M. B. Dahan, I. Bouchoule, Y. Castin, and C. Salomon, “Bloch oscillations of atoms, adiabatic rapid passage, and monokinetic atomic beams,” Phys. Rev. A 55, 2989–3001 (1997).

] and BEC [17

B. P. Anderson and M. A. Kasevich, “Macroscopic Quantum Interference from Atomic Tunnel Arrays,” Science 282, 1686 (1998). [CrossRef] [PubMed]

, 18

F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, “Josephson Junction Arrays with Bose-Einstein Condensates,” Science 293, 843 (2001). [CrossRef] [PubMed]

] of degenerate quantum gases (see [19

L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford: Clarendon Press, 2003).

] and references therein), where more sophisticated lattice geometries for trapping of ultracold atoms have been shown experimentally, e.g., the Kagomé lattice [20

L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, H. Fehrmann, and M. Lewenstein, “Atomic Quantum Gases in Kagomé Lattices,” Phys. Rev. Lett. 93, 030601 (2004). [CrossRef] [PubMed]

], the dice lattice [21

A. A. Burkov and E. Demler, “Vortex-Peierls States in Optical Lattices,” Phys. Rev. Lett. 96, 180406 (2006). [CrossRef] [PubMed]

], and quasi-periodic lattices [22

L. Sanchez-Palencia and L. Santos, “Bose-Einstein condensates in optical quasicrystal lattices,” Phys. Rev. A 72, 053607 (2005).

]. Simple theories are especially important for understanding the wave dynamics in the periodic structures. Recently, the theory of Zener tunnelling in the square lattices has been developed in [23

V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, “Zener tunneling in two-dimensional photonic lattices,” Phys. Rev. E. 74, 056602 (2006).

], moreover, the three-fold Bragg resonance in the hexagonal lattices has been studied as well [24

A. S. Desyatnikov, Yu. S. Kivshar, V. S. Shchesnovich, S. B. Cavalcanti, and J. M. Hickmann, “Resonant Zener tunneling in two-dimensional periodic photonic lattices,” Opt. Lett. 32, 325–327 (2007). [CrossRef] [PubMed]

The purpose of this paper is to give a complete theory of the Zener tunnelling in the hexagonal periodic lattices and test it by extensive numerical simulations. Though the perturbation theory for band-gap structures we use is a subject of textbooks on the solid state physics (see, for instance, Ref. [25

N. W. Ashcroft and N. D. Mermin, Solid state physics (NewYork: Holt, Rinehart and Winston, 1976).

]), any detailed study of Zener tunnelling in band gap structures in more than one spatial dimension has never been attempted. In particular, only the case of avoided crossing of two Bloch bands, long known since the works of Zener, Landau and Majorana, is discussed in literature. We choose as the illustrative example the hexagonal lattice for the two reasons. First, the hexagonal lattice, as distinct from the square lattice, cannot be separable, i.e. it cannot be represented as a sum of two one-dimensional potentials. Second, the hexagonal periodic lattice is advantageous in the two dimensional case and most of the two-dimensional photonic crystals have the hexagonal symmetry.

Due to the experimental advances in fabrication of periodic structures for both light and matter waves the interest in wave dynamics in two-dimensional lattices is growing. For instance, numerical studies of Bloch oscillations and Zener tunnelling [26

A. R. Kolovsky and H. J. Korsch, “Bloch oscillations of cold atoms in two-dimensional optical lattices,” Phys. Rev. A 67, 063601 (2003).

, 27

D. Witthaut, F. Keck, H. J. Korsch, and S. Mossmann, “Bloch oscillations in two-dimensional lattices,”New J. Phys. 6, 41 (2004). [CrossRef]

] were recently performed. It was also shown that the nonlinearity of the governing equations, in the matter waves case, is responsible for many new features. Some of them can be studied already in the one-dimensional lattices: asymmetry of the tunnelling probabilities, found in Refs. [28

B. Wu and Q. Niu, “Nonlinear Landau-Zener tunneling,” Phys. Rev. A 61, 023402 (2000).

, 29

O. Zobay and B. M. Garraway, “Time-dependent tunneling of Bose-Einstein condensates,” Phys. Rev. A 61, 033603 (2000).

] and shown experimentally [9

, 10

M. Jona-Lasinio, O. Morsch, M. Cristiani, E. Arimondo, and C. Menotti, “Nonlinear Effects for BoseEinstein Condensates in Optical Lattices,” Laser Phys. 15, 1180–1188 (2005).

], and modulational instability of Bloch waves [30

V. A. Brazhnyi, V. V. Konotop, and V. Kuzmiak, “Nature of the Intrinsic Relation between Bloch-Band Tunneling and Modulational Instability,” Phys. Rev. Lett. 96, 150402 (2006). [CrossRef] [PubMed]

] (see also Ref. [31

V. V. Konotop, P. G. Kevrekidis, and M. Salerno, “Landau-Zener tunneling of Bose-Einstein condensates in an optical lattice,” Phys. Rev. A 72, 023611 (2005).

]), also resulting in asymmetric nonlinear tunnelling. Moreover, a recent study of nonlinear tunnelling in the square two-dimensional lattices [32

V. A. Brazhnyi, V. V. Konotop, V. Kuzmiak, and V. S. Shchesnovich, “Nonlinear tunneling in two-dimensional lattices,” Phys. Rev. A 76, 023608 (2007).

] shows the existence of intraband tunnelling (see also Ref. [33

V. S. Shchesnovich and V. V. Konotop, “Nonlinear tunneling of Bose-Einstein condensates in an optical lattice: Signatures of quantum collapse and revival,” Phys. Rev. A 75, 063628 (2007).

, 34

V. S. Shchesnovich and V. V. Konotop, “Nonlinear intraband tunneling of a Bose-Einstein condensate in a cubic three-dimensional lattice,” Phys. Rev. A 77, 013614 (2008).

] for the quantum case), absent in the linear multi-dimensional case and in the one-dimensional nonlinear case.

Our attention is concentrated on the linear wave dynamics in hexagonal lattices, the nonlinear case is relegated to a future study. This imposes an upper bound on the nonlinearity, which is discussed below (see section 2). We make a first significant step forward in understanding of the phenomenon of multi-dimensional Zener tunnelling by reducing the Schrödinger equation to a system of ordinary differential equations. It is shown, for instance, that the two-dimensional tunnelling in the hexagonal lattices is described by either the three- or six-level Landau-Zener- Majorana (LZM) model with various interesting reductions in the special cases. Moreover, we derive analytical expressions for some of the transition probabilities by using the recent results on the multi-level LZM models [35

S. Brundobler and V. Elser, “S-matrix for generalized Landau-Zener problem,” J. Phys. A: Math. Gen. 26, 1211–1227 (1993). [CrossRef]

, 36

Yu. N. Demkov and V. N. Ostrovsky, “Multipath interference in a multistate Landau-Zener-type model,” Phys. Rev. A. 61, 032705 (2000). [CrossRef]

, 37

A. V. Shytov, “Landau-Zener transitions in a multilevel system: An exact result,” Phys. Rev. A 70, 052708 (2004).

The paper is organized as follows. In section 2 we introduce the model equation describing tunnelling in tilted two-dimensional lattices and discuss the limits of its applicability to BEC tunnelling and photonic interband transitions. Then, we discuss the shallow-lattice approximation and show that it allows to reduce the most general hexagonal periodic lattice to a simpler form. In sections 3 and 4 we discuss, respectively, the Rabi oscillations and Zener tunnelling. We derive simple LZM-type models and obtain the analytical formulae for the tunnelling probabilities. In section 5 we consider the interband tunnelling of vortices. Throughout the paper, the analytical models are compared with the numerical simulations. Section 6 provides some concluding remarks.

2. Shallow-lattice approximation

We consider the model describing the propagation of paraxial optical beams in a planar periodic photonic structure, e.g., a two dimensional lattice optically induced in photorefractive crystal [14

, 38

Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (San Diego: Academic, 2003).

i \partial_{t} Ψ + \frac{1}{2} \nabla_{⊥}^{2} Ψ + ϰ Δ n Ψ = 0,

(1)

where 𝝒 is the normalized nonlinear refractive index, ∇²=∂ ² _x+∂ ² _y and

Δ n = \frac{I_{g} (x) + I_{m} (x)}{1 + I_{g} (x) + I_{m} (x)}

, where x=(x,y), I _g (x) describes the optical lattice, and I _m (x) is a lattice tilt. We define the lattice potential as V(x)≡-𝝒 I _g (x) and assume the tilt to be linear, 𝝒 I _m (x)=-αx. Equation (1) in the shallow lattice approximation, |V(x)|≪1 (the main approximation used below), and for a weak linear tilt, |α‖d|/2≪|V(d/2)-V(0)|, reduces to the linear Schrödinger equation

i \partial_{t} Ψ = - \frac{1}{2} \nabla^{2} Ψ + [V (x) + α x] Ψ,

(2)

with an inessential constant term (see [23

V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, “Zener tunneling in two-dimensional photonic lattices,” Phys. Rev. E. 74, 056602 (2006).

] for more details). The dimensionless evolution variable t is the propagation distance in the case of an optical beam in a periodic photonic structure and the normalized time in the case of BEC. In the model, V(x) is the periodic lattice potential, V(x+d)=V(x), where d is one of the lattice periods, see Fig. 1(b),α is the acceleration of the lattice in the case of BEC and the steepness of the refractive index tilt in the case of the optical beam propagation. Below, we denote by Q a reciprocal lattice vector.

In the case of BEC, the Gross-Pitaevskii equation for the order parameter of BEC in a twodimensional optical lattice (for further details consult [19

L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford: Clarendon Press, 2003).

]) can be also reduced to Eq. (2) if we allow for the following conditions. The parabolic trap along the transverse z-direction is assumed to be much stronger than the nonlinearity of BEC, so that the condensate performs the ground-state quantum motion along this direction, while the trap in the (x,y)-plane is assumed to be weak and is neglected, i.e., it contains a large number of the lattice periods. The nonlinearity of BEC can be neglected if the Bloch oscillation period t _B is much less than the characteristic time t _nonl of the modulational instability development (see, for instance, [30

, 39

V. V. Konotop and M. Salerno, “Modulational instability in Bose-Einstein condensates in optical lattices,” Phys. Rev. A 65, 021602 (2002).

]) t _B ≪t _nonl~1/γ, where γ is proportional to the s-wave scattering length multiplied by the number of BEC atoms per lattice cell [19

L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford: Clarendon Press, 2003).

]. Since t _B ~|Q|/|α _‖| with Q=b _j (see Fig. 1(a)) and α _‖ the component of the acceleration parallel to Q, we obtain the condition in the form γ≪|αα‖d| (since |b _j |~1/|d|). Under the above conditions Eq. (2) describes the condensate tunneling in a titled two-dimensional optical lattice.

The shallow lattice approximation allows one to develop a fully analytical approach to the Zener interband transitions and reduce the governing equation to the multi-level LZM models. This approximation was already successfully used in [23

V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, “Zener tunneling in two-dimensional photonic lattices,” Phys. Rev. E. 74, 056602 (2006).

, 24

] and, moreover, showed an excellent qualitative agreement with the experimental results [14

] for an arbitrary lattice. On the other hand, for a deep lattice, which is the other limiting case, the continuous Eq. (2) is replaced by a set of coupled discrete equations for the amplitudes of Wannier states, which is the so-called tight-binding approximation [26

A. R. Kolovsky and H. J. Korsch, “Bloch oscillations of cold atoms in two-dimensional optical lattices,” Phys. Rev. A 67, 063601 (2003).

, 40

F. Keck and H. J. Korsch, “Infinite-variable Bessel functions in two-dimensional Wannier-Stark systems,” J. Phys. A: Math. Gen. 35, L105–L116 (2002). [CrossRef]

, 41

G. A. Alfimov, P. G. Kevrekidis, V. V. Konotop, and M. Salerno, “Wannier functions analysis of the nonlinear Schrdinger equation with a periodic potential,” Phys. Rev. E 66, 046608 (2002).

]. Clearly, the two limits are very different and assume different localization properties of the wave Ψ(x).

The interband transitions in a shallow lattice occur at the Bragg resonances, thus we assume that Ψ(x) has a narrow Fourier spectrum,

Ψ (x, 0) = \int_{D} d k C (k) e^{i kx},

(3)

where we require that the size of D be much smaller than the size of the Brillouin zone. In other terms, the initial condition is broad enough and covers at least several lattice periods. This means using broad optical beams in photonic crystals (as it was done in the recent experiment [14

]) and preparing the ultracold atoms or BEC in the state with a small momentum spread (as in the experiments of [9

, 10

M. Jona-Lasinio, O. Morsch, M. Cristiani, E. Arimondo, and C. Menotti, “Nonlinear Effects for BoseEinstein Condensates in Optical Lattices,” Laser Phys. 15, 1180–1188 (2005).

]). In this case, C(k) can be modelled by a sum of the Dirac delta-functions and the Zener tunnelling is described by an LZM model.

It is well known that at least three monochromatic plane waves

E_{l} = E_{0} e^{i k_{l} x - i ω t}

are necessary to produce a two-dimensional optical lattice, V=|∑E _l |². To have a hexagonal lattice the wave vectors should satisfy the corresponding symmetry, i.e., they should transform under the rotation by π/3 as follows: k ₁→k ₂, k ₂→k ₃, k ₃→-k ₁. The general expression for a hexagonal lattice (not necessarily created by an interference) can be given as the following Fourier series, infinite in general,

V = I_{0} {ε_{1} Σ_{l = 1}^{3} \cos (b_{l} x) + ε_{2} Σ_{l = 1}^{3} \cos (2 b_{l} x) +

ε_{3} [\cos ([b_{1} + b_{2}] x) + \cos ([b_{2} + b_{3}] x) + \cos ([b_{3} - b_{1}] x)] + \dots},

(4)

where the base vectors are b ₁=k ₁-k ₂, b ₂=k ₁-k ₃, and b ₃=k ₂-k ₃. Note that b ₃=b ₂-b ₁. On the r.h.s. of Eq. (4) we denote by the symbol “…” the omitted higher-order terms in the base vectors b _k, i.e., the cosines of n b _l x,n≥3, and the cosines of various sums of b ₁ x, b ₂ x, and b ₃ x, involving more than two terms (all cosines of the same kind have the same amplitude to satisfy the hexagonal rotation symmetry). The hexagonal lattice formed by using three plane waves, for instance, is given just by the first sum on the r.h.s. of Eq. (4): V=V ₀∑³ _l=₁cos(b l x).

Fig. 1. (a) The first Brillouin zone of the hexagonal lattice with b=|b _l |=2, where the vectors b _l correspond to translations between the points M, M′ and M″. (b) The hexagonal cos-lattice Eq. (5) and (c) corresponding “triangular” sin-lattice. The vectors d ₁ and d ₂ give two fundamental periods of the lattice. (d) The Bloch band structure (the first 9 bands) of the hexagonal cos-lattice in (b) with I ₀=0.1.

For ε ₁≠0 the fundamental periods of the lattice (4) can be selected as

d_{l} = \frac{2 π}{b} e_{1}

, l=1, 2, where b=|b _l | and the vectors e ₁ and e ₂ give the reciprocal basis for

n_{1} = \frac{b_{1}}{b}

and

n_{2} = \frac{b_{2}}{b}

: e _l n _m =δ _l , _m . Note also that

∣ e_{l} ∣ = \frac{2}{\sqrt{3}}

and

e_{1} e_{2} = - \frac{2}{3}

. We set b=2 (which can be arranged by an appropriate scaling of the spatial co-ordinate). In this case the vectors e ₁ and e ₂ connect the Γ-point with the M-points on the border of the Brillouin zone, see Fig. 1(a). We have b ₁=(2,0), b ₂=(1,√3), and b ₃=(-1,√3). Two important special cases of the general hexagonal lattice are considered below. In the first case,

V = 8 I_{0} \cos^{2} (\frac{b_{1} x}{2}) \cos^{2} (\frac{b_{2} x}{2}) \cos^{2} (\frac{b_{3} x}{2})

(5)

see Fig. 1(b), and a “triangular” lattice, which can be obtained from the lattice (5) by replacing all cosines by sines, see Fig. 1(c). The lattice (5) results from the setting ε ₁=3/2, ε ₂=1/4, and ε ₃=1/2 (all other terms are set to zero), while the “triangular” lattice corresponds to ε ₁=-1/2, ε ₂=-1/4, and ε ₃=1/2. Because the latter has the amplitude ε ₁ three times smaller, the efficiency of interband coupling is much lower in certain cases for the triangular lattices. Therefore, while our analytical theory is applicable to both lattices, we give numerical simulations for the hexagonal cosine lattice (5) only. Its Bloch band structure is given in Fig. 1(d).

3. Rabi oscillations

In the absence of the external bias α=0 the interband transition can be observed in the form of periodic Rabi oscillations. In the limit of shallow lattice, instead of expanding the solution over the Bloch waves, it is more convenient to use the Fourier version of the Schrödinger equation (2) directly. Setting Θ(x,t)=∫dk C(k,t)e ⁱ ^kx , expanding the lattice potential into the Fourier series,

V (x) = \sum_{Q} {\hat{V}}_{Q} e^{i Q x},

(6)

and representing the Fourier index as k=q-Q with q lying in the first Brillouin zone, one obtains the equation for the Fourier coefficients

i \partial_{t} C (q - Q) = E_{0} (q - Q) C (q - Q) + \sum_{Q'} {\hat{V}}_{Q' - Q} C (q - Q'),

(7)

where

E_{0} (q) = \frac{1}{2} q^{2}

. The usual resonant perturbation theory can be used to determine the Bloch bands [25

N. W. Ashcroft and N. D. Mermin, Solid state physics (NewYork: Holt, Rinehart and Winston, 1976).

]. In essence, the band gaps appear at the resonant Bragg planes defined as points q q where E ₀(q-Q)=E ₀(q). Most of the Fourier amplitudes are nonresonant and can be neglected in the first-order approximation. Then Eq. (7) predicts oscillations between the resonant peaks in Fourier space defined by Bragg condition E ₀(q-Q)=E ₀(q-Q′) [42

V. S. Shchesnovich and S. Chávez-Cerda, “Bragg-resonance-induced Rabi oscillations in photonic lattices,” Opt. Lett. 32, 1920–1922 (2007). [CrossRef] [PubMed]

]. To observe these oscillations the input wave (or beam) must be placed on the boundary of the first Brillouin zone. In the simplest case, i.e., away from the high-symmetry points, there are just two equivalent points on the boundary (Q/2 and -Q/2 below) and the oscillations are equivalent to those in the two-level system:

i \partial_{t} C (δ q + \frac{Q}{2}) = E_{0} (δ q + \frac{Q}{2}) C (δ q + \frac{Q}{2}) + {\hat{V}}_{Q} C (δ q - \frac{Q}{2}),

(8)

i \partial_{t} C (δ q - \frac{Q}{2}) = E_{0} (δ q - \frac{Q}{2}) C (δ q - \frac{Q}{2}) + {\hat{V}}_{Q}^{*} C (δ q + \frac{Q}{2}),

(9)

where we have used V̂ _-Q =V ^̂* _Q and introduced the running index δ _q for the points of the two resonant Fourier peaks. Since the Bloch waves are approximated as linear combinations of plane waves, the system (8)–(9) corresponds to the inter band oscillations between two Bloch bands. By analogy, these oscillations can be called Rabi oscillations. Rabi oscillations between several Bloch bands can be realized by placing the wave on one of the high-symmetry points on the boundary of the first Brillouin zone.

Fig. 2. Rabi oscillations between two X-points, see also movie Rabi.avi. Intensity of the initially Gaussian beam after propagation of t=6 in real (a) and Fourier (b) spaces. The dynamics of the powers P ₁,₂ of two oscillating beams is shown in (c), see text for the details. [Media 1]

Figure 2 and the movie Rabi.avi demonstrate the results of the simulations of Eq. (2) with a Gaussian beam as the initial condition, namely we use exp(-x ²/2w ²+i q ₀ x) with w=20 and q ₀=b ₁/₂, i.e., in Fourier space the beam is initially at the right X-point, see Fig. 1(a). We use spilt-step beam propagation method and monitor the dynamics also in the Fourier domain, see Fig. 2(b). The energy is periodically transferred between two X-points (±b ₁/2) and, for a quantitative comparison with the predictions of the LZM system (8–9), we integrate the intensities of two interacting beams in the Fourier domain to obtain their normalized powers P ₁,₂. The results presented in Fig. 2(c) allow to estimate, roughly, the period of these oscillations to be ~28.8. At the same time, the system (8–9) has a solution (δ q=0) in terms of harmonic functions with the period T=π/V̂Q=π/(ε ₁ I ₀)≈41.9, for the used parameters of the lattice I ₀=0.1 and ε ₁=3/2. Most probably, the disagreement is due to the transitional dynamics in Fig. 2 and averaging over longer propagation time provides better comparison with analytical predictions, similar as in the square lattice case [42

V. S. Shchesnovich and S. Chávez-Cerda, “Bragg-resonance-induced Rabi oscillations in photonic lattices,” Opt. Lett. 32, 1920–1922 (2007). [CrossRef] [PubMed]

In general, Rabi oscillations also show up in the real space as the oscillations of the average position of the beam. Besides performing oscillations, however, the beam propagates in the lattice due to the non-zero width of its Fourier image. Indeed, for δ q=0 there is the frequency mismatch between the corresponding Fourier amplitudes: E ₀(Q/2+δ q)-E ₀(-Q/2+δ _q)=Q δ q resulting in a higher average Fourier power of the input wave over the Bragg reflected one (see also [42

V. S. Shchesnovich and S. Chávez-Cerda, “Bragg-resonance-induced Rabi oscillations in photonic lattices,” Opt. Lett. 32, 1920–1922 (2007). [CrossRef] [PubMed]

]). In our simulations in Fig. 2 we choose relatively broad initial beam so that the oscillations and drift of its position are minimized.

4. Zener tunnelling

Since we consider the states (waves packets or beams)Ψ(x) with a narrowFourier spectrum, we can use the expansion over the Bloch waves with a definite t-dependent Bloch index q=q(t), similar to Houston’s approach [43

W. V. Houston, “Acceleration of Electrons in a Crystal Lattice,” Phys. Rev. 57, 184–186 (1940). [CrossRef]

] in the theory of accelerating electrons. In particular, we represent the solution to Eq. (2) in the following convenient form

Ψ (x, t) = \int_{D} d k C (k, t) e^{i (k - α t) x} .

(10)

This representation, derived from the plane wave solution in a linear potential Ψ(x,t)=exp{i(k ₀-α t)x-i∫ ^t dτ(k ₀-ατ)²/2}, is an alternative to having a partial derivative in the governing equation in Fourier space, which would account for the linear potential according to the rule x→i ∂ _k . Indeed, we get from Eq. (10)

i \partial_{t} C (q - Q) = E_{0} (q - Q - α t) C (q - Q) + \sum_{Q'} {\hat{V}}_{Q' - Q} C (q - Q^{'}) .

(11)

By switching to the “interaction picture” of the perturbation theory (where the lattice is considered a perturbation) by the transformation C(q-Q,t)=e ^-i ∫ ^t dτE ₀(q-Q-ατ)C̃(q-Q,t) one sees that any two coefficients C(q-Q) and C( q - Q ′) are effectively coupled on the time interval where the property E ₀( q - Q -α t)≈E ₀(q-Q′- α t) is satisfied, otherwise the coupling coefficients are oscillating about zero. Thus the interband transitions take place on the Bragg resonance planes. It is convenient to explicitly account for the Bragg resonance by defining the resonant point q _res by setting q - Q - α t= q _res- Q - α (t-t ₀) in the energy E ₀. For simplicity, below we set t ₀=0 (in this case one obtains the governing LZM-type models in the standard form). There are three types of Zener interband transitions in the hexagonal lattices: (i) the quasi one-dimensional tunnelling (section 4.1), (ii) tunnelling between three Bloch bands at the M-point (section 4.2), and (iii) between the six Bloch bands at the Γ-point (section 4.3).

4.1. Quasi-one-dimensional Zener transitions

The quasi one-dimensional Bragg resonance takes place when the Bloch index crosses the Bragg planes outside the small neighborhoods of the high-symmetry points Γ and M (of radius R~I ₀). This is a transition between two Bloch bands at an avoided crossing, which takes place along one of the borders of the irreducible Brillouin zone, i.e. the triangle ΓMX in Fig. 1, to which the Bragg plane is equivalent (after performing rotations by multiples of π/3 and translations by the reciprocal lattice vectors). Let us briefly recall the derivation, which is similar as in the case of the square lattice [23

V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, “Zener tunneling in two-dimensional photonic lattices,” Phys. Rev. E. 74, 056602 (2006).

] (see also [28

B. Wu and Q. Niu, “Nonlinear Landau-Zener tunneling,” Phys. Rev. A 61, 023402 (2000).

, 29

O. Zobay and B. M. Garraway, “Time-dependent tunneling of Bose-Einstein condensates,” Phys. Rev. A 61, 033603 (2000).

, 44

V. S. Shchesnovich and S. B. Cavalcanti, “Finite-dimensional model for the condensate tunnelling in an accelerating optical lattice,” J. Phys. B: At. Mol. Opt. Phys. 39, 1997–2011 (1997). [CrossRef]

] for the one-dimensional case). By keeping only the resonant terms in the potential and in the Bloch wave, we arrive at the following expressions:

V_{j} = \frac{ε_{j} I_{0}}{2} (e^{i Q_{j} x} + e^{- i Q_{j} x}), Ψ_{j} = C_{1} (t) {e^{i q}}^{(t) x} + C_{2} (t) e^{i [q (t) - Q_{j}] x} .

(12)

The Bragg resonance point is q _j =Q _j /2. Substituting the expression for the Bloch wave into Eq. (2) and requiring cancellation of the terms linear in x, which gives q=-α(the dot denotes derivative with respect to t), we get a system of coupled equations for the incident C ₁ and Bragg reflected C ₂ amplitudes of the Bloch wave. By setting q =Q _j /2- α t and (C ₁,C ₂)=e ^{i
ϕ(t)}(c ₁,c ₂), with ϕ=(q ²(t)+[q(t)- Q _j ]²)/4, the latter system is cast in the form of Zener [2

C. Zener, “Non-Adiabatic Crossing of Energy Levels,” Proc. R. Soc. London A 137, 696–702 (1932).

], Landau [4

L. D. Landau, “On the theory of transfer of energy at collisions II,” Phys. Z. Sowjetunion 2, 46 (1932).

], and Majorana [5

E. Majorana, “Orientated atoms in a variable magnetic field,” Nuovo Cimento 9, 43–50 (1932). [CrossRef]

i {\dot{c}}_{1} = - \frac{(Q_{j} α) t}{2} c_{1} + \frac{ε_{j} I_{0}}{2} c_{2}, i {\dot{c}}_{2} = - \frac{(Q_{j} α) t}{2} c_{2} + \frac{ε_{j} I_{0}}{2} c_{1} .

(13)

System (13) is Hamiltonian with the adiabatic energy levels corresponding to the two Bloch bands E at an avoided crossing:

E_{1, 2} = \mp \sqrt{q^{2} + {(\frac{ε_{j} I_{0}}{2})}^{2}}

, where q≡(Q _j α )t/2 is the running band parameter. All quasi one-dimensional cases are treated in the same way: there is a single resonant reciprocal lattice vector Q which defines a Bragg plane and the probability of tunnelling is given by the LZM formula:

P = \exp {- \frac{2 π {∣ {\hat{V}}_{Q} ∣}^{2}}{∣ Q ∣ ∣ α_{⊥} ∣}},

(14)

with α _⊥ being the component of the bias perpendicular to the respective Bragg plane. Thus, Eq. (14) gives the probability of transition between the two Bloch bands E ₁,₂(q), defined for the j-th band as P _j =|c _j(∞)|² for the initial condition |c _j (-∞)|=1. Due to the symmetry t→-t and c ₁→c ^* ₂, c ₂→c ^* ₁ the two probabilities coincide. Note also that, in the case of a sufficiently small bias | α |≪(ε ₁ I ₀)², Zener tunnelling is negligible when crossing the border of the first Brillouin zone (away from the high-symmetry point M). This means that the wave (beam) is totally Bragg-reflected at the boundary of the first Brillouin zone, i.e., it performs Bloch oscillations. There are three distinct cases corresponding to the three inequivalent Bragg planes, each one being equivalent to one side of the irreducible Brillouin zone. Recall that in the derivation of the LZM model (13) we have subtracted the mean value of the Bloch energy Ē=(Q _j )²/8 (by using the phase transformation) which specifies the two bands involved in the tunnelling transition.

I. Zener transitions on the XM-border. In this case the resonant reciprocal lattice vector is Q=b ₁ (equivalently b ₂ or b ₃). The resonant Fourier components of the lattice are those with the coefficient ε ₁ in formula (4) and the LZM formula (14) becomes

P = \exp {- \frac{π ε_{1}^{2} I_{0}^{2}}{4 ∣ α_{1} ∣}}

. This is the transition between the first and the second Bloch bands (see Fig. 1(d)), in fact, it occurs when the Fourier index crosses the border of the first Brillouin zone away from the M point.

II. Zener transitions on the ΓM-border. The resonant reciprocal lattice vector is Q=2b ₁ (or its equivalent). The resonant Fourier components of the lattice are those with the coefficient ε ₂. The shift by - b ₁ and rotation by -π/3 transforms the resonance point to a point lying on the ΓM-border. The probability of tunnelling reads

P = \exp {- \frac{π ε_{2}^{2} I_{0}^{2}}{4 ∣ \sqrt{3} α_{2} - α_{1} ∣}}

. This is the transition between the second and the third Bloch bands along the ΓM-border, in contrast to the above transition, it occurs outside the first Brillouin zone in the extended zone picture.

III. Zener transitions on the ΓX-border. The resonant reciprocal lattice vector is equivalent to Q= b ₁+ b ₂ and the resonant Fourier components are those with the coefficient ε ₃. The rotation by π/3 followed by the shift by -b ₃ transforms the resonance point to a point lying on the ΓX-border. We obtain the probability as follows

P = \exp {- \frac{π ε_{3}^{2} I_{0}^{2}}{4 \sqrt{3} ∣ α_{2} ∣}}

. This is the transition between the third and the fourth Bloch bands along the ΓX-border and it also occurs outside the first Brillouin zone in the extended zone picture.

Fig. 3. Quasi-one-dimensional tunnelling through the X-point with the tilt directed along b ₁, see also movie 1D.avi. (a,b) Intensity in the real (top) and Fourier (bottom) spaces are shown for t=0 in (a) and t=10 in (b), corresponding dynamics of beam powers is shown in (d, solid lines). Solutions to the LZM system Eq. (13) are shown in (c) and (d, dashed lines). See text for the details and parameter values. [Media 2]

To illustrate the dynamics of the one-dimensional Zener tunnelling we choose the simples case of the tunnelling through the X-point (case I. above), similar to Rabi oscillations in Fig. 2. We show in Fig. 3 and the movie 1D.avi results of the simulations with the parameters |b _j |=2, I ₀=0.1,| α |=0.05, and ε ₁=3/2. The LZM result here is P=exp(-9π/80)≃0.7 and it is shown with horizontal dashed line in Fig. 2(c), together with the solution to Eq. (13) with initial condition c ₁(-∞)=1 and c ₂(-∞)=0. The numerical simulations with a finite beam in Fig. 3(d, solid lines) reproduce the qualitative features of the asymptotic LZM result in Fig. 3(c). To have a quantitative comparison with the numerics one has to match the initial conditions. Instead of the asymptotic initial condition at minus infinity, a better account is given by the LZM system Eq. (13) with the initial condition at t ₀=(q _X-q ₀)/| α |, i.e. c ₁(t ₀)=1 and c ₂(t ₀)=0. By doing so, we were able to find an excellent agreement between direct numerical results (solid lines in Fig. 2(d)) and the solution to LZM system (dashed lines), including the value of the tunnelling efficiency ~0.77.

4.2. Zener transitions between three Bloch bands

Bragg resonance at the M-point is described by the three-fold resonance. Indeed, from Fig. 1 one can see that an M-point of the hexagonal lattice has exactly two equivalent points, M′ and M″, such that the translations between all these points are given by the reciprocal lattice vectors b _j . To derive the corresponding LZM model we note that the resonant terms in the lattice potential and in the Bloch wave can be cast as follows:

V_{res} = \frac{ε_{1} I_{0}}{2} [e^{i (q_{M'} - q_{M}) x} + e^{i (q_{M ″} - q_{M}) x} + e^{i (q_{M ″} - q_{M'}) x} + c . c .],

(15)

Ψ = C_{1} (t) e^{i q (t) x} + C_{2} (t) e^{i [q (t) + q_{M'} - q_{M}] x} + C_{3} (t) e^{i [q (t) + q_{M ″} - q_{M}] x},

(16)

where q _M, q _M′, and q _M″ are the vectors connecting the Γ-point with the corresponding M-point. The Fourier amplitudes of the Bloch wave are related to the M-points as follows: C ₁→M, C ₂→M′, and C ₃→M″. Setting C _j =e ^iϕ ( ^t )c _j , with ϕ=(q ² _M+ α ² t ²)/2, and proceeding as in the derivation of the system (13) we obtain the three-level system

i {\dot{c}}_{1} = - (v_{1} t) c_{1} + \frac{ε_{1} I_{0}}{2} (c_{2} + c_{3}),

(17)

{i \dot{c}}_{2} = - (v_{2} t) c_{2} + \frac{ε_{1} I_{0}}{2} (c_{1} + c_{3}),

(18)

{i \dot{c}}_{3} = - (v_{3} t) c_{3} + \frac{ε_{1} I_{0}}{2} (c_{1} + c_{2}),

(19)

where ν ₁= q _M α, ν ₂=q _M′ α , and ν ₃= q _M″ α . The coefficients ν _j satisfy the obvious constraint ν₁+ν₂+ν₃=0. The invariance of the system (17)–(19) with respect to the π/3-rotation is evident from the corresponding transformation ν₁→ν₂→ν₃→ν₁. The coefficients νj can be cast as follows: ν ₁=ℓ| α |cosθ, ν ₂=-ℓ| α |(cosθ-√3sinθ)/2, and _ν 3=-ℓ|α|(cosθ+√3sinθ)/2, whereℓ=| q _M|=b/√3 with b≡| b _l | (in our case b=2) is the length of the hexagon side in Fig. 1 and θ is the polar angle of the bias direction counting from the ΓM-line (the choice of a particular ΓM-line just sets the ordering of the Fourier coefficients c _j ). Since the transformation θ→-θ, i.e., c ₂→c ₃ and c ₃→c ₂, leaves the system (17)–(19) invariant, the polar angle can be restricted to 0≤θ≤π/3.

The first three Bloch bands along the direction of the tilt in the vicinity of the M-point (see also Fig. 1(d)) are given by the set of the adiabatic energy levels of the system (17)–(19). Though the analytical result for the band structure in the vicinity of the M-point is quite complicated, at the M-point itself the values are given by (compare with Fig. 1(d))

E_{res} (M) = \frac{b^{2}}{6} + [- \frac{ε_{1} I_{0}}{2}; - \frac{ε_{1} I_{0}}{2}; ε_{1} I_{0}],

(20)

where we have used that | q _M|=b/√3. Hence, in the case of the hexagonal lattice (4), two of the three Bloch bands collide at the M-point (in the case of the lattice (5) the first two bands collide, while in the case of the “triangular” lattice the second and third bands collide at the M-point). It is not known whether the n-level LZM system admits an analytical solution for all transition probabilities in the general case (see [35

S. Brundobler and V. Elser, “S-matrix for generalized Landau-Zener problem,” J. Phys. A: Math. Gen. 26, 1211–1227 (1993). [CrossRef]

, 36

Yu. N. Demkov and V. N. Ostrovsky, “Multipath interference in a multistate Landau-Zener-type model,” Phys. Rev. A. 61, 032705 (2000). [CrossRef]

, 37

A. V. Shytov, “Landau-Zener transitions in a multilevel system: An exact result,” Phys. Rev. A 70, 052708 (2004).

] for further discussion). However, the probabilities of two transitions are known in the analytical form. These are the transitions m→m, where m is the diagonal index of the minimal or maximal value of the coefficient at t in the system (which must be unique). The probabilities are given by the general formula [35

S. Brundobler and V. Elser, “S-matrix for generalized Landau-Zener problem,” J. Phys. A: Math. Gen. 26, 1211–1227 (1993). [CrossRef]

, 37

A. V. Shytov, “Landau-Zener transitions in a multilevel system: An exact result,” Phys. Rev. A 70, 052708 (2004).

] (P _m → _m ≡|c _m (∞)|² for |c _m(-∞)|=1)

P_{m \to m} = e x p {- 2 π \sum_{j \neq m} \frac{{∣ Δ_{m j} ∣}^{2}}{∣ v_{m} - v_{j} ∣}},

(21)

where ν _j is the coefficient at t in the j-th equation in the system and Δ _mj is the coupling coefficient between the j-th and m-th equations. In the generic case, i.e., when 0<θ<π/3, we have ν1>ν ₂>ν ₃ in Eqs. (17)–(19) and the transition probabilities P ₁→₁ and P ₃→₃ are

P_{1 \to 1} = e x p {- \frac{π ε_{1}^{2} I_{0}^{2}}{2} [\frac{1}{v_{1} - v_{2}} + \frac{1}{v_{1} - v_{3}}]},

(22)

P_{3 \to 3} = e x p {- \frac{π ε_{1}^{2} I_{0}^{2}}{2} [\frac{1}{v_{1} - v_{3}} + \frac{1}{v_{2} - v_{3}}]} .

(23)

Below we consider in more details two special cases, θ=0 and θ=π/3, when Eqs. (17)–(19) can be reduced to the LZM system for two amplitudes.

Fig. 4. (a) The structure of the Bloch bands along the ΓM-direction in a neighborhood of the M-point. The distance between the quasi-degenerate Bloch bands (E ₂-E ₁ for -t≫1 and E ₃-E ₁ for t≫1) is equal to 3ε ₁ I ₀/4. (b–d) Symmetric tunnelling through the M-point [24

] with the tilt directed along b ₁+b ₂, see text and the movie . [Media 3]

I. Special case θ=0: symmetric three-fold Zener transitions.

In this case (the tilt is along the ΓM-line) we have ν≡-ν ₂=-ν ₃=ν ₁/2 with ν=l| α |/2 and it is convenient to introduce new amplitudes related to the original Bloch amplitudes by orthogonal transformation b ₁=(c ₂-c ₃)/√2, b ₂=(c ₂+c ₃)/√2, and b ₃=c ₁. The amplitude b ₁ decouples from the system, while b ₂ and b ₃ remain coupled:

i {\dot{b}}_{1} = (v t - \frac{ε_{1} I_{0}}{2}) b_{1},

(24)

i {\dot{b}}_{2} = (v t + \frac{ε_{1} I_{0}}{2}) b_{2} + \frac{ε_{1} I_{0}}{\sqrt{2}} b_{3},

(25)

i {\dot{b}}_{3} = - 2 v t b_{3} + \frac{ε_{1} I_{0}}{\sqrt{2}} b_{2} .

(26)

Equation (24) is readily solved, b ₁(t)=b ₁(0)exp(-i νt ²/2+i ε ₁ I ₀ t/2), and the system (25)–(26) has an additional integral of motion, |b ₂|²+|b ₃|²≡C=const. We employ the initial conditions c ₁(-∞)=1 and c ₂,₃(-∞)=0, so that C=1. It follows from the initial condition that b ₁(t)=0 and, hence, c ₂=c ₃=b ₂/√2 for all t. The conservation of norm can be expressed now as |c ₂,₃|²=(1-|c ₁|²)/2 and we refer to this case as the symmetric tunnelling. It was analyzed numerically in [24

] and here we would like to stress several additional details.

Along the ΓM-line the Bloch bands have the structure about the M-point shown in Fig. 4(a), where E ₁=(ν ₁ t-ε ₁ I ₀)/2,

E_{2, 3} = \frac{(- ν_{1} t + ε_{1} I_{0} \mp 3 \sqrt{D (t)})}{4}

, with D(t)=ν ² ₁ t ²+ε ² ₁ I ² ₀+2ε ₁ ν ₁ t/3 (D(t)>0 for ε ₁ I ₀<3, i.e., for a shallow lattice). The tunnelling probability between the Bloch bands corresponding to the system (25)–(26) is still given by the probability P ₁→₁ of Eq. (22), which in this case becomes the LZM result

P = \exp {- \frac{2 π ε_{1}^{2} I_{0}^{2}}{3 ν_{1}}}

. An example of tunnelling through the M-point is given in Fig. 4(b,c), see also [24

] and the movie 3D0.avi. It is interesting to note that the three-beam interference in Fig. 4(c) reproduces the wave used to optically induce the lattice itself [45

A. S. Desyatnikov, N. Sagemerten, R. Fischer, B. Terhalle, D. Träger, D. N. Neshev, A. Dreischuh, C. Denz, W. Krolikowski, and Yu. S. Kivshar “Two-dimensional self-trapped nonlinear photonic lattices,” Opt. Express 14, 2851–2863 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-7-2851. [CrossRef] [PubMed]

], together with its regular phase structure shown in Fig. 4(d). The honeycomb pattern of phase dislocations, or optical vortices [46

M. S. Soskin and M. V. Vasnetsov, “Singular Optics,” Progress in Optics 42, 219 (2001). [CrossRef]

], is the signature of the corresponding Bloch wave [47

T. J. Alexander, A. S. Desyatnikov, and Yu. S. Kivshar, “Multivortex solitons in triangular photonic lattices,” Opt. Lett. 32, 1293–1295 (2007). [CrossRef] [PubMed]

] at the M-point of hexagonal lattice.

We conclude that Zener tunnelling can be employed for the excitation of specific Bloch waves and thus characterization of photonic structures. In the nonlinear regime optical vortices can form localized states, or vortex solitons [48

A. S. Desyatnikov, Yu. S. Kivshar, and L. Torner, “Optical Vortices and Vortex Solitons,” Progress in Optics 47, 291–391 (2005). [CrossRef]

], and novel multi-vortex states bifurcate from the vortex lattice [47

T. J. Alexander, A. S. Desyatnikov, and Yu. S. Kivshar, “Multivortex solitons in triangular photonic lattices,” Opt. Lett. 32, 1293–1295 (2007). [CrossRef] [PubMed]

]. Similar effects should be expected to occur with BEC and can be of interest for the excitation of vortex lattices.

Fig. 5. Tunnelling of the Rabi oscillations. Solutions to the LZM system in (a) and (d) are compared to numerical simulations in (b), (c), and (e); see also movie 3D3.avi. [Media 4]

II. Special case θ=π/3: Zener tunnelling of Rabi oscillations. We have in this case ν ₁=ν ₂=-ν ₃/2=-ν (note that now ν=-l| _α |/2). Introducing the amplitudes b ₁=(c ₁-c ₂)/√2, b ₂=(c ₁+c ₂)/√2, and b ₃=c ₃ we obtain formally the same system (24)–(26) (in fact, complex conjugate to it with inverted time t→-t, if one takes into account change of sign of ν). The initial conditions now read b ₁(0)=b ₂(0)=1/√2 and b ₃(0)=0, so that |b ₂|²+|b ₃|²=1/2. It follows from the solution to Eq. (24) that c ₁(t)=c ₂(t)+exp(-iνt ²/2+i ε ₁ I ₀ t/2).

Although now the reduction of essentially asymmetric three-fold resonance to the two-level LZM system is not so obvious as in previous case, its meaning becomes clear from the solution presented in Fig. 5(a). The Fourier components c ₁ and c ₂ form a Rabi oscillating state, described by the two-level system (17)–(18) (equivalent up to a phase transformation to the system (8)–(9) with Q=b ₁/2 and V̂ _Q =ε ₁ I ₀/2) before and after the three-fold resonance crossing at t~0, where they couple to the third Fourier amplitude c ₃ (away from the resonance the coupling is non-resonant and can be neglected, i.e., c ₃ rapidly varies with respect to c ₁ and c ₂). In the equivalent Bloch band population interpretation the Rabi oscillations are between the first and the second Bloch bands in Fig. 1(d) (along the XM-line) before the resonance crossing and between the second and third bands (now along the ΓM-line) after the crossing. Before the resonance, the first two band population amplitudes are, in fact, given by b ₁ and b ₂ (and the third by b ₃=c ₃). After the resonance crossing b ₁ and b ₂ now give the second and third band populations, while b ₃ is the population of the first band. As the result of the resonance, a part of the Rabi energy (given by |b ₃|²) is left in the first band and the amplitude of oscillations is reduced (recall the conservation law |b ₂|²=1/2-|b ₃|²). Thus the tunnelling of the oscillating state, described by the LZM system (25)–(26), can be thought of as Zener tunnelling of Rabi oscillations.

Fig. 6. (a) The six resonant Γ-points (indicated by Γ _j with _j =1, …, 6) forming an extended Brillouin zone (solid hexagon). The reciprocal lattice vectors relating the Γ₁-point to the other resonant Γ-points are indicated. (b–d) solutions to the LZM systems: (b) corresponds to Eq. (32), (c) to system 1 of (35) and (d) to system 2 of (35). The dashed lines show the analytical results; c ₁(0)=1 in all cases.

Figure 5(b, c) shows the results of numerical simulations of the Gaussian beam of a finite size, while the tunneled beam powers are plotted in Fig. 5(e), see also movie 3D3.avi. Although the short time which our main moving beam spends near M-point is not sufficient to observe Rabi oscillations, we could successfully match these results with corresponding solution to the LZM system in Fig. 5(e). Note that after the tunnelling starts to develop, the LZM system shows the onset of oscillations between c ₁ and c ₂. In contrast, the powers P ₁ and P ₂ in (e) remain close because these two beams remain no longer in the vicinity of resonant M-points. Nevertheless, the efficiency of tunnelling to the level P3~0.277 is in perfect agreement with LZM model.

4.3. Zener transitions between six Bloch bands

In general, the Bragg resonance at the Γ-point is six-fold due to the fact that each of the reciprocal lattice vectors {± b ₁,± b ₂,± b ₃} if taken from the center of the Brillouin zone would point at one of the six nearest equivalent Γ-points, see Fig. 6(a), and the Bragg vectors giving the translations between these Γ-points correspond to non-zero Fourier amplitudes of the lattice (4) if all ε _j ≠0.

Let us assume that the resonance is at t=0 and the Bloch index is q(0)=q _Γ1=ΓΓ₁, i.e., q=q _Γ ₁-α t. Then the corresponding resonant terms of the Bloch wave can be cast as follows

Ψ = C_{1} (t) e^{i q (t) x} + \sum_{l = 2}^{6} C_{l} (t) e^{i {[q (t) + q_{Γ_{l}} - q_{Γ_{1}}]}^{x}},

(27)

where q _Γ _l is the reciprocal lattice vector connecting the center of the Brillouin zone with the Γ _l -point. Setting C _j =e ^iϕ ^(t) c _j , with ϕ=( q ² _Γ ₁ +α ² t ²)/2, and following the same procedure as above we get the system of equations for the Fourier amplitudes c _j of the Bloch wave. Defining C=(c ₁, …,c ₆) ^T we have

i \dot{C} = (Λ t + \frac{I_{0}}{2} H) C, H = (\begin{matrix} A & B \\ B & A \end{matrix}),

(28)

with Λ=diag(-λ ₁,-λ ₂,-λ ₃,λ ₁,λ2,λ ₃), where λ _j = b _j α, _j =1,2, 3, and

A = (\begin{matrix} 0 & ε_{1} & ε_{3} \\ ε_{1} & 0 & ε_{1} \\ ε_{3} & ε_{1} & 0 \end{matrix}), B = (\begin{matrix} ε_{2} & ε_{3} & ε_{1} \\ ε_{3} & ε_{2} & ε_{3} \\ ε_{1} & ε_{3} & ε_{2} \end{matrix}) .

The matrix H is a circulant, i.e., it is a Toeplitz matrix, H _i , _j =h(i-_j), with the property h(l±6)=h(l). This, together with the transformation λ ₁→λ ₂→λ ₃→-λ ₁ under the rotation by π/3, guarantees the invariance of the system (28) under the rotations multiple of π/3. Due to b ₃=b ₂-b ₁, the coefficients λ _j satisfy the identity λ ₃=λ ₂-λ ₁. They can be given as follows: λ ₁=b| α |cosφ, λ ₂=b| α |(cosφ+√3sinφ)/2, and λ ₃=b| α |(-cosφ+√3sinφ)/2, with φ being the polar angle of the bias direction counting from the ΓΓ ₁-line. The angle φ can be restricted to 0≤φ<π/3 by the rotation symmetry.

The Bloch band structure about the Γ-point corresponding to the system (28), i.e., the six resonant bands, cannot be obtained in the analytical form. However, at the Γ-point itself the six resonant Bloch bands assume the following values

E_{res} (Γ) = \frac{b^{2}}{2} + \frac{I_{0}}{2} [Δ_{1}; Δ_{1}; Δ_{2}; Δ_{2}; Δ_{3}; Δ_{4}],

(29)

with Δ₁=ε ₂-ε ₁-ε ₃, Δ₂=ε ₁-ε ₂-ε ₃, Δ₃=2(ε ₃-ε ₁)-ε ₂, and Δ₄=2(ε ₃+ε ₁)+ε ₂. For the lattice (5) we get (we use b=2, compare with Fig. 1(d))

E = 2 + [- \frac{9 I_{0}}{8}; - \frac{7 I_{0}}{8}; - \frac{7 I_{0}}{8}; \frac{3 I_{0}}{8}; \frac{3 I_{0}}{8} \frac{17 I_{0}}{8}]

while for the “triangular” lattice we have

E = 2 + [- \frac{3 I_{0}}{8}; - \frac{3 I_{0}}{8}; - \frac{I_{0}}{8}; - \frac{I_{0}}{8}; - \frac{I_{0}}{8}; \frac{9 I_{0}}{8}] .

In the case when there is a single minimal or maximal value of ±λ _j two tunnelling probabilities corresponding to the system (28) can be given in the analytic form by using formula (21). Noticing that λ ₁>λ ₃ and λ₂≥-λ₃ (the equality sign for φ=0) we have the two possible orderings of the coefficients:

λ1>λ2≥-λ3>λ3≥-λ2>-λ1 for 0≤φ<π/6,

λ2≥λ1>λ3≥-λ3>-λ1≥-λ2 for π/6≤φ<π/3.

Thus, from the general result Eq. (21) we obtain:

P_{1 \to 1} = P_{4 \to 4} = e x p {- \frac{π I_{0}^{2}}{2} [\frac{ε_{1}^{2}}{λ_{1} - λ_{2}} + \frac{ε_{1}^{2}}{λ_{1} - λ_{3}} + \frac{ε_{3}^{2}}{λ_{1} - λ_{3}} + \frac{ε_{3}^{2}}{λ_{1} - λ_{2}} + \frac{ε_{2}^{2}}{2 λ_{1}}]},

(30)

for 0<φ<π/6 and

P_{2 \to 2} = P_{5 \to 5} = e x p {- \frac{π I_{0}^{2}}{2} [\frac{ε_{1}^{2}}{λ_{2} - λ_{1}} + \frac{ε_{1}^{2}}{λ_{2} - λ_{3}} + \frac{ε_{3}^{2}}{λ_{1} - λ_{2}} + \frac{ε_{3}^{2}}{λ_{2} - λ_{3}} + \frac{ε_{2}^{2}}{2 λ_{2}}]} .

(31)

for π/6<φ<π/3. Under the rotation by π/3, due to the transformationλ ₁→λ ₂→λ ₃→-λ ₁, the probabilities in Eq. (30) go into those of Eq. (31), as it should be.

The most interesting cases of the six-fold resonance are for the special (i.e., symmetric in the lattice) bias directions, when the six-level system (28) decouples into two subsystems. There just two such cases. In the case of φ=0 (the tilt is along the ΓΓ₁-line, i.e., in the direction of b ₁) the system decouples into a two-level system and a four-level system, whereas in the case φ=π/6 (the tilt is in the direction of b ₁+b ₂) the general system decouples into two three-level systems. Since the special cases could be relevant for future experiments, we consider them in detail.

Special case φ=0. In this case we haveλ₁=b|α|, λ ₂=λ₁/2 andλ₃=-λ₁/2. Introducing new amplitudes b _l by an orthogonal transformation, b ₁=c ₁, b ₂=(c ₂+c ₆)/√2, b ₃=(c ₃+c ₅)/√2, b ₄=c ₄, b ₅=(c ₂-c ₆)/√2 and b ₆=(c ₃-c ₅)/√2, we get for B ₀=(b ₁,b ₂,b ₃,b ₄) ^T a four-level system

i \dot{B_{0}} = (Λ_{0} t + I_{0} H_{0}) B_{0}, H_{0} \equiv (\begin{matrix} 0 & \frac{ε_{1}}{\sqrt{2}} & \frac{ε_{3}}{\sqrt{2}} & \frac{ε_{2}}{2} \\ \frac{ε_{1}}{\sqrt{2}} & \frac{ε_{3}}{2} & \frac{ε_{1} + ε_{2}}{2} & \frac{ε_{3}}{\sqrt{2}} \\ \frac{ε_{3}}{\sqrt{2}} & \frac{ε_{1} + ε_{2}}{2} & \frac{ε_{3}}{2} & \frac{ε_{1}}{\sqrt{2}} \\ \frac{ε_{2}}{2} & \frac{ε_{3}}{\sqrt{2}} & \frac{ε_{1}}{\sqrt{2}} & 0 \end{matrix})

(32)

with Λ₀=diag(-λ ₁,-λ ₁/2,λ ₁/2,λ ₁). The rest two amplitudes, b ₅ and b ₆, are decoupled from the system (32), they evolve according to the LZM system:

i \dot{b_{5}} = (- \frac{λ_{1} t}{2} - \frac{ε_{3} I_{0}}{2}) b_{5} + \frac{I_{0}}{2} (ε_{1} - ε_{2}) b_{6},

(33)

i \dot{b_{6}} = (\frac{λ_{1} t}{2} - \frac{ε_{3} I_{0}}{2}) b_{6} + \frac{I_{0}}{2} (ε_{1} - ε_{2}) b_{5} .

(34)

Fig. 7. Symmetric tunnelling through the Γ-point with the tilt directed along b ₁. As before and in Fig. 6 we use ε ₁=3/2, ε ₂=1/4, ε ₃=1/2, I ₀=0.1, and | α |=0.05. [Media 5]

We see that Zener transitions take place between either four or two Bloch bands, depending on the initial populations, i.e., the initial values of b _j . In the case of tunnelling between the two bands the probability is given by the LZM result

P = \exp {- \frac{π I_{0}^{2}}{2 λ_{1}} {(ε_{1} - ε_{2})}^{2}}

. More interesting is the case of tunnelling involving four bands with initial condition c ₁(0)=b ₁(0)=1, so that, solving Eqs. (33)–(34), we obtain b ₅,₆=0 and c ₂=c ₆, c ₃=c ₅. The probability of the transition P ₁→₁=P ₄→₄ is known, and, in fact, follows from Eq. (30) which becomes

P = \exp {- \frac{2 π I_{0}^{2}}{λ_{1}} [ε_{1}^{2} + \frac{ε_{2}^{2}}{8} + \frac{ε_{3}^{2}}{3}]}

. This case is illustrated in Fig. 6(b) and the results of the direct numerical simulations are presented in Fig. 7 and movie 6D.avi. It is clearly seen that tunnelling is indeed symmetric with respect to reflection in x-axis (i.e., c ₂=c ₆, c ₃=c ₅; x-axis corresponds to the b ₁-direction).

Special case φ=π/6. We have

λ_{1} = \frac{\sqrt{3}}{2} b ∣ α ∣, λ_{2} = λ_{1}

and λ ₃=0. Therefore, in this case system (28) decouples into two three-level systems for the amplitudes given by the following orthogonal transformation b ₁=(c ₁+c ₂)/√2, b ₂=(c ₃+c ₆)/√2, b ₃=(c ₄ +c ₅)/√2, b ₄= (c ₁-c ₂)/√2, b ₅=(c ₃-c ₆)/√2, b ₆=(c ₄-c ₅)/√2. We get:

i \dot{B_{j}} = (Λ_{j} z + \frac{I_{0}}{2} H_{j}) B_{j}, j = 1, 2,

(35)

with B ₁=(b ₁,b ₂,b ₃) ^T , B ₂=(b ₄,b ₅,b ₆) ^T , Λ ₁=Λ ₂=diag(-λ ₁,0,λ₁), and

H_{1} = (\begin{matrix} ε_{1} & ε_{1} + ε_{3} & ε_{2} + ε_{3} \\ ε_{1} + ε_{3} & ε_{2} & ε_{1} + ε_{3} \\ ε_{2} + ε_{3} & ε_{1} + ε_{3} & ε_{1} \end{matrix}), H_{2} = (\begin{matrix} - ε_{1} & ε_{3} - ε_{1} & ε_{2} - ε_{3} \\ ε_{3} - ε_{1} & - ε_{2} & ε_{1} - ε_{3} \\ ε_{2} - ε_{3} & ε_{1} - ε_{3} & - ε_{1} \end{matrix}) .

In both cases, the probability of the transition P ^(j)≡P ^(j) ₁→₁=P ^(j) ₃→₃ is known, it reads

P^{(j)} = \exp {- \frac{π I_{0}^{2}}{2 λ_{1}} [{(ε_{1} \pm ε_{3})}^{2} + \frac{1}{2} {(ε_{2} \pm ε_{3})}^{2}]},

(36)

with the plus sign for j=1 and the minus for j=2. For the initial condition c ₁(-∞)=1 the tunnelling involves all six Bloch bands since in this case we have b ₁(-∞)=b ₄(-∞)=1, see Fig. 6(c,d).

5. Vortex tunnelling

The theory developed in section 4 can also be extended to the interband transitions of optical beams carrying phase dislocations, e.g. optical vortices [46

M. S. Soskin and M. V. Vasnetsov, “Singular Optics,” Progress in Optics 42, 219 (2001). [CrossRef]

]. Consider, for instance, a Gaussian beam carrying a single vortex of the charge n,

Ψ_{n} = \frac{{[(x - x_{0}) + i (y - y_{0})]}^{n}}{\sqrt{2 π} σ} \exp {i k_{0} (x - x_{0}) - \frac{{(x - x_{0})}^{2}}{2 σ^{2}}} .

(37)

The main point is based on the fact that a vortex in the real space corresponds to a vortex of the same vorticity (or topological charge) in the Fourier space according to the rule

F {Ψ_{n}} (k) = {[i \partial_{k} - \partial_{λ} - (x_{n} + {iy}_{n})]}^{n} F {Ψ_{0}} (k),

(38)

e.g. the Gaussian beam Eq. (37) in the Fourier space reads

F {Ψ_{n}} = \sqrt{2 π} σ {(- i σ^{2})}^{n} {[(k - k_{0}) + i (λ - λ_{0})]}^{n} \exp {- i x_{0} (k - k_{0}) - \frac{σ^{2}}{2} {(k - k_{0})}^{2}},

(39)

where k=k n _x +λ n _y . In the absence of the lattice, but with the linear potential alone, V _lin= α x, the wave (37) expands (diffracts), σ²(t)=σ²+it, acquiring a chirp and an overall phase according to the law x ₀(t)=x ₀+∫t ₀ k ₀(t)dt, k ₀(t)=k ₀-α t. Note that a vortex initially placed close to the beam center x=x ₀ remains close to x ₀(t) for all times.

First, we demonstrate that the Zener tunnelling of a beam carrying a phase dislocation in a periodic lattice results in several output beams with the phase dislocations of the same charge. As an example, we consider the case of a single n-vortex in a Gaussian beam, i.e., the n-order Laguerre-Gauss beam. The Bragg resonance will result in several identical Gaussian peaks in the Fourier space with the indices shifted by the resonant reciprocal lattice vectors. For instance, near the resonance at t=0, for the output Fourier amplitude with the index q-Q we have

C(q-Q,t)=-iV̂QC(q,0)t +O(t2),

where C(q,0) is the Fourier image of the input wave (in this case, q=k0, it has the Fourier index inside the first Brillouin zone). Therefore, from Eq. (11) and Eq. (39) we conclude that every output peak in the Fourier space tagged by q-Q will have an n-vortex placed on it at its center index k=k ₀(t)-Q. It is important to note that the position of the output beamC(q-Q), x ₀(t)=x ₀+(k ₀-Q)t-α t ²/2, is carrier-dependent through Q. Returning to the real space, we obtain an n-vortex sitting on each of the output waves (whose velocities are different). Thus the phase dislocations resulting from the tunnelling are identical and preserved in the evolution between the resonances. The vortices can be replaced by anti-vortices without any affect on the distribution of power in the Fourier space.

Fig. 8. Three-fold symmetric Zener tunnelling of a vortex beam, see the movie V-1.avi. (a) Intensity and (b) phase of the initial beam with the topological charge m=-1 in the vicinity of the M-point, see (c). Intensities in the Fourier domain (c,d) are identical for both, vortex, m=+1, and antivortex, m=-1, while they differ significantly in the real space, compare (e,g) and corresponding phases (f,h). [Media 6]

We have performed numerical simulations using the Laguerre-Gaussian vortex beams for the case of symmetric three-fold resonance, similar to Fig. 4 and [24

]; the results are presented in Fig. 8. In the Fourier space, the intensity distribution does not depend on the vortex topological charge, see Fig. 8(c,d), and the vortices have the same charge as the initial beam. We do not observe any difference in the tunnelling efficiency, simply repeating the results with the Gaussian beams in Fig. 4. In the real space, however, due to interference, the two solutions corresponding to the tunnelling of vortices and that of anti-vortices are substantially different, see Fig. 8(e,g), since the phases of the carriers are the same while the vortex phases are conjugated, when a vortex is replaced by an anti-vortex. At the same time, the phase of the tunneled vortex beam is simply of the opposite charge, see the characteristic fork-type phase dislocations in Fig. 8(f,h) similar to the one in Fig. 8(b). Note also that the dynamics of anti-vortices is equivalent to that of vortices in the complex conjugate Schrödinger equation, which requires the inversion of time t→-t.

6. Conclusions

We have analyzed the interband transitions, such as Zener tunnelling and Rabi oscillations, in hexagonal photonic lattices and derived several different multi-level LZM systems that capture the essence of these phenomena. We have demonstrated that the direction of the tilt together with the Fourier coefficients of the lattice potential determine how many of the lowest-order Bloch bands of the photonic bandgap spectrum are involved in the interband transitions. We have identified three general regimes of the Zener tunnelling in the hexagonal photonic lattice: (i) quasi one-dimensional Zener tunnelling (or, equivalently, simple Bragg resonance involving only two Bloch bands) which occurs when the Bloch index crosses the Bragg planes far from one of the high-symmetry points; (ii) three-fold Bragg resonance at the high-symmetry M-point with the Zener transitions between the three Bloch bands; and (iii) the six-fold Bragg resonance at the high-symmetry Γ-point with the Zener tunnelling involving, in general, six Bloch bands. For some special directions of the tilt, the tunnelling trough the Γ-point is described by either two three-level systems or by a combination of one four-level and one two-level systems; the six Bloch bands decouple into two groups of either two pairs of three bands or four bands and two bands, respectively, with the Zener tunnelling taking place only between the bands of the same group. In the special symmetric cases, we have found reductions of the general LZM models and found several interesting features of the tunneling dynamics, for instance, we have found a new effect, which can be called Zener tunneling of Rabi oscillations. In addition, we have shown that tunnelling of phase dislocations (or optical vortices) results in the output waves carrying the same phase dislocations. In the real space, there is an asymmetry in the interference pictures resulting from the tunnelling of vortices and anti-vortices, since the dynamics of vortices is equivalent to that of anti-vortices in the complex-conjugated equation with the inverted time.

Our theoretical results have immediate applications to a number of different phenomena such as the light beam propagation in optically-induced photonic lattices, tunnelling of Bose-Einstein condensates in optical lattices, and other types of wave propagation in tilted hexagonal periodic potentials. Our results may also provide an additional field for the search for the integrableLZM systems (see, e.g., [35

S. Brundobler and V. Elser, “S-matrix for generalized Landau-Zener problem,” J. Phys. A: Math. Gen. 26, 1211–1227 (1993). [CrossRef]

, 36

Yu. N. Demkov and V. N. Ostrovsky, “Multipath interference in a multistate Landau-Zener-type model,” Phys. Rev. A. 61, 032705 (2000). [CrossRef]

, 37

A. V. Shytov, “Landau-Zener transitions in a multilevel system: An exact result,” Phys. Rev. A 70, 052708 (2004).

] and references therein). Indeed, one would expect that the multimode LZM systems obtained from the highly-symmetric hexagonal lattices are integrable, in the sense that they allow an analytical expression for the S-matrix (as defined in [35

S. Brundobler and V. Elser, “S-matrix for generalized Landau-Zener problem,” J. Phys. A: Math. Gen. 26, 1211–1227 (1993). [CrossRef]

]), so that the transition probabilities between all levels could be found in an analytic form.

We should notice, however, that in the derivation of the LZM models, we have neglected the effects of nonlinearity since our aim was to elucidate the role played by the lattice symmetry on Zener interband transitions. In the case of optical beams propagating in the photonic crystals, the effect of nonlinearity can be negligible indeed, while nonlinearity of BEC has been found to cause the breakdown of adiabaticity in the Zener tunnelling [28

B. Wu and Q. Niu, “Nonlinear Landau-Zener tunneling,” Phys. Rev. A 61, 023402 (2000).

] and change the symmetry between the upper-to-lower vs. lower-to-upper interband transitions [30

, 28

B. Wu and Q. Niu, “Nonlinear Landau-Zener tunneling,” Phys. Rev. A 61, 023402 (2000).

, 29

O. Zobay and B. M. Garraway, “Time-dependent tunneling of Bose-Einstein condensates,” Phys. Rev. A 61, 033603 (2000).

, 31

V. V. Konotop, P. G. Kevrekidis, and M. Salerno, “Landau-Zener tunneling of Bose-Einstein condensates in an optical lattice,” Phys. Rev. A 72, 023611 (2005).

]. Nonlinear tunnelling in the hexagonal photonic lattices will be considered in one of our future publications (see also Ref. [32

V. A. Brazhnyi, V. V. Konotop, V. Kuzmiak, and V. S. Shchesnovich, “Nonlinear tunneling in two-dimensional lattices,” Phys. Rev. A 76, 023608 (2007).

] where the nonlinear tunnelling has been considered for the case of a square lattice).

Acknowledgments

The work has been supported by the CAPES and FAPESP of Brazil (VSS) and the Australian Research Council in Australia (ASD and YSK). VSS also thanks the Nonlinear Physics Center of the Australian National University for support and hospitality.

References and links

1.	F. Bloch, “Ü ber die Quantenmechanik der Elektronen in Kristallgittern,” Z. Phys. 52, 555–600 (1928).
2.	C. Zener, “Non-Adiabatic Crossing of Energy Levels,” Proc. R. Soc. London A 137, 696–702 (1932).
3.	G. H. Wannier, “Wave Functions and Effective Hamiltonian for Bloch Electrons in an Electric Field,” Phys. Rev. 117, 432–439 (1960). [CrossRef]
4.	L. D. Landau, “On the theory of transfer of energy at collisions II,” Phys. Z. Sowjetunion 2, 46 (1932).
5.	E. Majorana, “Orientated atoms in a variable magnetic field,” Nuovo Cimento 9, 43–50 (1932). [CrossRef]
6.	L. Esaki, “Long journey into tunneling,” Rev. Mod. Phys. 46, 237–244 (1974). [CrossRef]
7.	B. Bourlon, D. C. Glattli, B. Placais, J. M. Berroir, C. Miko, L. Forro, and A. Bachtold, “Geometrical Dependence of High-Bias Current in Multiwalled Carbon Nanotubes,” Phys. Rev. Lett. 92, 026804 (2004). [CrossRef] [PubMed]
8.	A. Sibille, J. F. Palmier, and F. Laruelle, “Zener Interminiband Resonant Breakdown in Superlattices,” Phys. Rev. Lett. 80, 4506–4509 (1998). [CrossRef]
9.	M. Jona-Lasinio, O. Morsch, M. Cristiani, N. Malossi, J. H. Müller, E. Courtade, M. Anderlini, and E. Arimondo, “Asymmetric Landau-Zener Tunneling in a Periodic Potential,” Phys. Rev. Lett. 91, 230406 (2003). [CrossRef] [PubMed]
10.	M. Jona-Lasinio, O. Morsch, M. Cristiani, E. Arimondo, and C. Menotti, “Nonlinear Effects for BoseEinstein Condensates in Optical Lattices,” Laser Phys. 15, 1180–1188 (2005).
11.	G. Malpuech, A. Kavokin, G. Panzarini, and A. Di Carlo, “Theory of photon Bloch oscillations in photonic crystals,” Phys. Rev. B 63, 035108 (2001).
12.	V. Agarwal, J. A. del Río., G. Malpuech, M. Zamfirescu, A. Kavokin, D. Coquillat, M. Scalbert, M. Vladimirova, and B. Gil, “Photon Bloch Oscillations in Porous Silicon Optical Superlattices,” Phys. Rev. Lett . 92, 097401 (2004). [CrossRef] [PubMed]
13.	H. Trompeter, T. Pertsch, F. Lederer, D. Michaelis, U. Streppel, A. Bräuer, and U. Peschel, “Visual Observation of Zener Tunneling,” Phys. Rev. Lett. 96, 023901 (2006). [CrossRef] [PubMed]
14.	H. Trompeter, W. Krolikowski, D. N. Neshev, A. S. Desyatnikov, A. A. Sukhorukov, Yu. S. Kivshar, T. Pertsch, U. Peschel, and F. Lederer, “Bloch Oscillations and Zener Tunneling in Two-Dimensional Photonic Lattices,” Phys. Rev. Lett. 96, 053903 (2006). [CrossRef] [PubMed]
15.	M. B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, “Bloch Oscillations of Atoms in an Optical Potential,” Phys. Rev. Lett. 76, 4508–4511 (1996). [CrossRef] [PubMed]
16.	E. Peik, M. B. Dahan, I. Bouchoule, Y. Castin, and C. Salomon, “Bloch oscillations of atoms, adiabatic rapid passage, and monokinetic atomic beams,” Phys. Rev. A 55, 2989–3001 (1997).
17.	B. P. Anderson and M. A. Kasevich, “Macroscopic Quantum Interference from Atomic Tunnel Arrays,” Science 282, 1686 (1998). [CrossRef] [PubMed]
18.	F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, “Josephson Junction Arrays with Bose-Einstein Condensates,” Science 293, 843 (2001). [CrossRef] [PubMed]
19.	L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford: Clarendon Press, 2003).
20.	L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, H. Fehrmann, and M. Lewenstein, “Atomic Quantum Gases in Kagomé Lattices,” Phys. Rev. Lett. 93, 030601 (2004). [CrossRef] [PubMed]
21.	A. A. Burkov and E. Demler, “Vortex-Peierls States in Optical Lattices,” Phys. Rev. Lett. 96, 180406 (2006). [CrossRef] [PubMed]
22.	L. Sanchez-Palencia and L. Santos, “Bose-Einstein condensates in optical quasicrystal lattices,” Phys. Rev. A 72, 053607 (2005).
23.	V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, “Zener tunneling in two-dimensional photonic lattices,” Phys. Rev. E. 74, 056602 (2006).
24.	A. S. Desyatnikov, Yu. S. Kivshar, V. S. Shchesnovich, S. B. Cavalcanti, and J. M. Hickmann, “Resonant Zener tunneling in two-dimensional periodic photonic lattices,” Opt. Lett. 32, 325–327 (2007). [CrossRef] [PubMed]
25.	N. W. Ashcroft and N. D. Mermin, Solid state physics (NewYork: Holt, Rinehart and Winston, 1976).
26.	A. R. Kolovsky and H. J. Korsch, “Bloch oscillations of cold atoms in two-dimensional optical lattices,” Phys. Rev. A 67, 063601 (2003).
27.	D. Witthaut, F. Keck, H. J. Korsch, and S. Mossmann, “Bloch oscillations in two-dimensional lattices,”New J. Phys. 6, 41 (2004). [CrossRef]
28.	B. Wu and Q. Niu, “Nonlinear Landau-Zener tunneling,” Phys. Rev. A 61, 023402 (2000).
29.	O. Zobay and B. M. Garraway, “Time-dependent tunneling of Bose-Einstein condensates,” Phys. Rev. A 61, 033603 (2000).
30.	V. A. Brazhnyi, V. V. Konotop, and V. Kuzmiak, “Nature of the Intrinsic Relation between Bloch-Band Tunneling and Modulational Instability,” Phys. Rev. Lett. 96, 150402 (2006). [CrossRef] [PubMed]
31.	V. V. Konotop, P. G. Kevrekidis, and M. Salerno, “Landau-Zener tunneling of Bose-Einstein condensates in an optical lattice,” Phys. Rev. A 72, 023611 (2005).
32.	V. A. Brazhnyi, V. V. Konotop, V. Kuzmiak, and V. S. Shchesnovich, “Nonlinear tunneling in two-dimensional lattices,” Phys. Rev. A 76, 023608 (2007).
33.	V. S. Shchesnovich and V. V. Konotop, “Nonlinear tunneling of Bose-Einstein condensates in an optical lattice: Signatures of quantum collapse and revival,” Phys. Rev. A 75, 063628 (2007).
34.	V. S. Shchesnovich and V. V. Konotop, “Nonlinear intraband tunneling of a Bose-Einstein condensate in a cubic three-dimensional lattice,” Phys. Rev. A 77, 013614 (2008).
35.	S. Brundobler and V. Elser, “S-matrix for generalized Landau-Zener problem,” J. Phys. A: Math. Gen. 26, 1211–1227 (1993). [CrossRef]
36.	Yu. N. Demkov and V. N. Ostrovsky, “Multipath interference in a multistate Landau-Zener-type model,” Phys. Rev. A. 61, 032705 (2000). [CrossRef]
37.	A. V. Shytov, “Landau-Zener transitions in a multilevel system: An exact result,” Phys. Rev. A 70, 052708 (2004).
38.	Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (San Diego: Academic, 2003).
39.	V. V. Konotop and M. Salerno, “Modulational instability in Bose-Einstein condensates in optical lattices,” Phys. Rev. A 65, 021602 (2002).
40.	F. Keck and H. J. Korsch, “Infinite-variable Bessel functions in two-dimensional Wannier-Stark systems,” J. Phys. A: Math. Gen. 35, L105–L116 (2002). [CrossRef]
41.	G. A. Alfimov, P. G. Kevrekidis, V. V. Konotop, and M. Salerno, “Wannier functions analysis of the nonlinear Schrdinger equation with a periodic potential,” Phys. Rev. E 66, 046608 (2002).
42.	V. S. Shchesnovich and S. Chávez-Cerda, “Bragg-resonance-induced Rabi oscillations in photonic lattices,” Opt. Lett. 32, 1920–1922 (2007). [CrossRef] [PubMed]
43.	W. V. Houston, “Acceleration of Electrons in a Crystal Lattice,” Phys. Rev. 57, 184–186 (1940). [CrossRef]
44.	V. S. Shchesnovich and S. B. Cavalcanti, “Finite-dimensional model for the condensate tunnelling in an accelerating optical lattice,” J. Phys. B: At. Mol. Opt. Phys. 39, 1997–2011 (1997). [CrossRef]
45.	A. S. Desyatnikov, N. Sagemerten, R. Fischer, B. Terhalle, D. Träger, D. N. Neshev, A. Dreischuh, C. Denz, W. Krolikowski, and Yu. S. Kivshar “Two-dimensional self-trapped nonlinear photonic lattices,” Opt. Express 14, 2851–2863 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-7-2851. [CrossRef] [PubMed]
46.	M. S. Soskin and M. V. Vasnetsov, “Singular Optics,” Progress in Optics 42, 219 (2001). [CrossRef]
47.	T. J. Alexander, A. S. Desyatnikov, and Yu. S. Kivshar, “Multivortex solitons in triangular photonic lattices,” Opt. Lett. 32, 1293–1295 (2007). [CrossRef] [PubMed]
48.	A. S. Desyatnikov, Yu. S. Kivshar, and L. Torner, “Optical Vortices and Vortex Solitons,” Progress in Optics 47, 291–391 (2005). [CrossRef]

OCIS Codes
(050.0050) Diffraction and gratings : Diffraction and gratings
(050.4865) Diffraction and gratings : Optical vortices

ToC Category:
Diffraction and Gratings

History
Original Manuscript: May 19, 2008
Revised Manuscript: August 19, 2008
Manuscript Accepted: August 20, 2008
Published: August 26, 2008

Citation
Valery S. Shchesnovich, Anton S. Desyatnikov, and Yuri S. Kivshar, "Interband resonant transitions in two-dimensional hexagonal lattices: Rabi oscillations, Zener tunnelling, and tunnelling of phase dislocations," Opt. Express 16, 14076-14094 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-18-14076

Sort: Author | Year | Journal | Reset

References

F. Bloch, "Uber die Quantenmechanik der Elektronen in Kristallgittern," Z. Phys. 52, 555-600 (1928).
C. Zener, "Non-Adiabatic Crossing of Energy Levels," Proc. R. Soc. London A 137, 696-702 (1932).
G. H. Wannier, "Wave Functions and Effective Hamiltonian for Bloch Electrons in an Electric Field," Phys. Rev. 117, 432-439 (1960). [CrossRef]
L. D. Landau, "On the theory of transfer of energy at collisions II," Phys. Z. Sowjetunion 2, 46 (1932).
E. Majorana, "Orientated atoms in a variable magnetic field," Nuovo Cimento 9, 43-50 (1932). [CrossRef]
L. Esaki, "Long journey into tunneling," Rev. Mod. Phys. 46, 237-244 (1974). [CrossRef]
B. Bourlon, D. C. Glattli, B. Placais, J. M. Berroir, C. Miko, L. Forro, and A. Bachtold, "Geometrical Dependence of High-Bias Current in Multiwalled Carbon Nanotubes," Phys. Rev. Lett. 92, 026804 (2004). [CrossRef] [PubMed]
A. Sibille, J. F. Palmier, and F. Laruelle, "Zener Interminiband Resonant Breakdown in Superlattices," Phys. Rev. Lett. 80, 4506-4509 (1998). [CrossRef]
M. Jona-Lasinio, O. Morsch, M. Cristiani, N. Malossi, J. H. M¨uller, E. Courtade, M. Anderlini, and E. Arimondo, "Asymmetric Landau-Zener Tunneling in a Periodic Potential," Phys. Rev. Lett. 91, 230406 (2003). [CrossRef] [PubMed]
M. Jona-Lasinio, O. Morsch, M. Cristiani, E. Arimondo, and C. Menotti, "Nonlinear Effects for BoseEinstein Condensates in Optical Lattices," Laser Phys. 15, 1180-1188 (2005).
G. Malpuech, A. Kavokin, G. Panzarini, and A. Di Carlo, "Theory of photon Bloch oscillations in photonic crystals," Phys. Rev. B 63, 035108 (2001).
V. Agarwal, J. A. del Rıo, G. Malpuech, M. Zamfirescu, A. Kavokin, D. Coquillat, D. Scalbert, M. Vladimirova, and B. Gil, "Photon Bloch Oscillations in Porous Silicon Optical Superlattices," Phys. Rev. Lett. 92, 097401 (2004). [CrossRef] [PubMed]
H. Trompeter, T. Pertsch, F. Lederer, D. Michaelis, U. Streppel, A. Brauer, and U. Peschel, "Visual Observation of Zener Tunneling," Phys. Rev. Lett. 96, 023901 (2006). [CrossRef] [PubMed]
H. Trompeter, W. Krolikowski, D. N. Neshev, A. S. Desyatnikov, A. A. Sukhorukov, Yu. S. Kivshar, T. Pertsch, U. Peschel, and F. Lederer, "Bloch Oscillations and Zener Tunneling in Two-Dimensional Photonic Lattices," Phys. Rev. Lett. 96, 053903 (2006). [CrossRef] [PubMed]
M. B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, "Bloch Oscillations of Atoms in an Optical Potential," Phys. Rev. Lett. 76, 4508-4511 (1996). [CrossRef] [PubMed]
E. Peik, M. B. Dahan, I. Bouchoule, Y. Castin, and C. Salomon, "Bloch oscillations of atoms, adiabatic rapid passage, and monokinetic atomic beams," Phys. Rev. A 55, 2989-3001 (1997).
B. P. Anderson and M. A. Kasevich, "Macroscopic Quantum Interference from Atomic Tunnel Arrays," Science 282, 1686 (1998). [CrossRef] [PubMed]
F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, "Josephson Junction Arrays with Bose-Einstein Condensates," Science 293, 843 (2001). [CrossRef] [PubMed]
L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford: Clarendon Press, 2003).
L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, H. Fehrmann, and M. Lewenstein, "Atomic Quantum Gases in Kagome Lattices," Phys. Rev. Lett. 93, 030601 (2004). [CrossRef] [PubMed]
A. A. Burkov and E. Demler, "Vortex-Peierls States in Optical Lattices," Phys. Rev. Lett. 96, 180406 (2006). [CrossRef] [PubMed]
L. Sanchez-Palencia and L. Santos, "Bose-Einstein condensates in optical quasicrystal lattices," Phys. Rev. A 72, 053607 (2005).
V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, "Zener tunneling in two-dimensional photonic lattices," Phys. Rev. E. 74, 056602 (2006).
A. S. Desyatnikov, Yu. S. Kivshar, V. S. Shchesnovich, S. B. Cavalcanti, and J. M. Hickmann, "Resonant Zener tunneling in two-dimensional periodic photonic lattices," Opt. Lett. 32, 325-327 (2007). [CrossRef] [PubMed]
N. W. Ashcroft and N. D. Mermin, Solid state physics (NewYork: Holt, Rinehart and Winston, 1976).
A. R. Kolovsky and H. J. Korsch, "Bloch oscillations of cold atoms in two-dimensional optical lattices," Phys. Rev. A 67, 063601 (2003).
D. Witthaut, F. Keck, H. J. Korsch, and S. Mossmann, "Bloch oscillations in two-dimensional lattices," New J. Phys. 6, 41 (2004). [CrossRef]
B. Wu and Q. Niu, "Nonlinear Landau-Zener tunneling," Phys. Rev. A 61, 023402 (2000).
O. Zobay and B. M. Garraway, "Time-dependent tunneling of Bose-Einstein condensates," Phys. Rev. A 61, 033603 (2000).
V. A. Brazhnyi, V. V. Konotop, and V. Kuzmiak, "Nature of the Intrinsic Relation between Bloch-Band Tunneling and Modulational Instability," Phys. Rev. Lett. 96, 150402 (2006). [CrossRef] [PubMed]
V. V. Konotop, P. G. Kevrekidis, and M. Salerno, "Landau-Zener tunneling of Bose-Einstein condensates in an optical lattice," Phys. Rev. A 72, 023611 (2005).
V. A. Brazhnyi, V. V. Konotop, V. Kuzmiak, and V. S. Shchesnovich, "Nonlinear tunneling in two-dimensional lattices," Phys. Rev. A 76, 023608 (2007).
V. S. Shchesnovich and V. V. Konotop, "Nonlinear tunneling of Bose-Einstein condensates in an optical lattice: Signatures of quantum collapse and revival," Phys. Rev. A 75, 063628 (2007).
V. S. Shchesnovich and V. V. Konotop, "Nonlinear intraband tunneling of a Bose-Einstein condensate in a cubic three-dimensional lattice," Phys. Rev. A 77, 013614 (2008).
S. Brundobler and V. Elser, "S-matrix for generalized Landau-Zener problem," J. Phys. A: Math. Gen. 26, 1211- 1227 (1993). [CrossRef]
Yu. N. Demkov and V. N. Ostrovsky, "Multipath interference in a multistate Landau-Zener-type model," Phys. Rev. A. 61, 032705 (2000). [CrossRef]
A. V. Shytov, "Landau-Zener transitions in a multilevel system: An exact result," Phys. Rev. A 70, 052708 (2004).
Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (San Diego: Academic, 2003).
V. V. Konotop and M. Salerno, "Modulational instability in Bose-Einstein condensates in optical lattices," Phys. Rev. A 65, 021602 (2002).
F. Keck and H. J. Korsch, "Infinite-variable Bessel functions in two-dimensional Wannier-Stark systems," J. Phys. A: Math. Gen. 35, L105-L116 (2002). [CrossRef]
G. A. Alfimov, P. G. Kevrekidis, V. V. Konotop, and M. Salerno, "Wannier functions analysis of the nonlinear Schrdinger equation with a periodic potential," Phys. Rev. E 66, 046608 (2002).
V. S. Shchesnovich and S. Chavez-Cerda, "Bragg-resonance-induced Rabi oscillations in photonic lattices," Opt. Lett. 32, 1920-1922 (2007). [CrossRef] [PubMed]
W. V. Houston, "Acceleration of Electrons in a Crystal Lattice," Phys. Rev. 57, 184-186 (1940). [CrossRef]
V. S. Shchesnovich and S. B. Cavalcanti, "Finite-dimensional model for the condensate tunnelling in an accelerating optical lattice," J. Phys. B: At. Mol. Opt. Phys. 39, 1997-2011 (1997). [CrossRef]
A. S. Desyatnikov, N. Sagemerten, R. Fischer, B. Terhalle, D. Tr¨ager, D. N. Neshev, A. Dreischuh, C. Denz, W. Krolikowski, and Yu. S. Kivshar, "Two-dimensional self-trapped nonlinear photonic lattices," Opt. Express 14, 2851-2863 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-7-2851. [CrossRef] [PubMed]
M. S. Soskin and M. V. Vasnetsov, "Singular Optics," Progress in Optics 42, 219 (2001). [CrossRef]
T. J. Alexander, A. S. Desyatnikov, and Yu. S. Kivshar, "Multivortex solitons in triangular photonic lattices," Opt. Lett. 32, 1293-1295 (2007). [CrossRef] [PubMed]
A. S. Desyatnikov, Yu. S. Kivshar, and L. Torner, "Optical Vortices and Vortex Solitons," Progress in Optics 47, 291-391 (2005). [CrossRef]

J. Phys. A: Math. Gen.

S. Brundobler and V. Elser, "S-matrix for generalized Landau-Zener problem," J. Phys. A: Math. Gen. 26, 1211- 1227 (1993). [CrossRef]
F. Keck and H. J. Korsch, "Infinite-variable Bessel functions in two-dimensional Wannier-Stark systems," J. Phys. A: Math. Gen. 35, L105-L116 (2002). [CrossRef]

J. Phys. B: At. Mol. Opt. Phys.

V. S. Shchesnovich and S. B. Cavalcanti, "Finite-dimensional model for the condensate tunnelling in an accelerating optical lattice," J. Phys. B: At. Mol. Opt. Phys. 39, 1997-2011 (1997). [CrossRef]

Laser Phys.

M. Jona-Lasinio, O. Morsch, M. Cristiani, E. Arimondo, and C. Menotti, "Nonlinear Effects for BoseEinstein Condensates in Optical Lattices," Laser Phys. 15, 1180-1188 (2005).

New J. Phys.

D. Witthaut, F. Keck, H. J. Korsch, and S. Mossmann, "Bloch oscillations in two-dimensional lattices," New J. Phys. 6, 41 (2004). [CrossRef]

Nuovo Cimento

E. Majorana, "Orientated atoms in a variable magnetic field," Nuovo Cimento 9, 43-50 (1932). [CrossRef]

Opt. Express

A. S. Desyatnikov, N. Sagemerten, R. Fischer, B. Terhalle, D. Tr¨ager, D. N. Neshev, A. Dreischuh, C. Denz, W. Krolikowski, and Yu. S. Kivshar, "Two-dimensional self-trapped nonlinear photonic lattices," Opt. Express 14, 2851-2863 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-7-2851. [CrossRef] [PubMed]

Opt. Lett.

Phys. Rev.

G. H. Wannier, "Wave Functions and Effective Hamiltonian for Bloch Electrons in an Electric Field," Phys. Rev. 117, 432-439 (1960). [CrossRef]
W. V. Houston, "Acceleration of Electrons in a Crystal Lattice," Phys. Rev. 57, 184-186 (1940). [CrossRef]

Phys. Rev. A

A. V. Shytov, "Landau-Zener transitions in a multilevel system: An exact result," Phys. Rev. A 70, 052708 (2004).
A. R. Kolovsky and H. J. Korsch, "Bloch oscillations of cold atoms in two-dimensional optical lattices," Phys. Rev. A 67, 063601 (2003).
V. V. Konotop and M. Salerno, "Modulational instability in Bose-Einstein condensates in optical lattices," Phys. Rev. A 65, 021602 (2002).
E. Peik, M. B. Dahan, I. Bouchoule, Y. Castin, and C. Salomon, "Bloch oscillations of atoms, adiabatic rapid passage, and monokinetic atomic beams," Phys. Rev. A 55, 2989-3001 (1997).
V. V. Konotop, P. G. Kevrekidis, and M. Salerno, "Landau-Zener tunneling of Bose-Einstein condensates in an optical lattice," Phys. Rev. A 72, 023611 (2005).
V. A. Brazhnyi, V. V. Konotop, V. Kuzmiak, and V. S. Shchesnovich, "Nonlinear tunneling in two-dimensional lattices," Phys. Rev. A 76, 023608 (2007).
V. S. Shchesnovich and V. V. Konotop, "Nonlinear tunneling of Bose-Einstein condensates in an optical lattice: Signatures of quantum collapse and revival," Phys. Rev. A 75, 063628 (2007).
V. S. Shchesnovich and V. V. Konotop, "Nonlinear intraband tunneling of a Bose-Einstein condensate in a cubic three-dimensional lattice," Phys. Rev. A 77, 013614 (2008).
B. Wu and Q. Niu, "Nonlinear Landau-Zener tunneling," Phys. Rev. A 61, 023402 (2000).
O. Zobay and B. M. Garraway, "Time-dependent tunneling of Bose-Einstein condensates," Phys. Rev. A 61, 033603 (2000).
L. Sanchez-Palencia and L. Santos, "Bose-Einstein condensates in optical quasicrystal lattices," Phys. Rev. A 72, 053607 (2005).

Phys. Rev. A.

Yu. N. Demkov and V. N. Ostrovsky, "Multipath interference in a multistate Landau-Zener-type model," Phys. Rev. A. 61, 032705 (2000). [CrossRef]

Phys. Rev. B

G. Malpuech, A. Kavokin, G. Panzarini, and A. Di Carlo, "Theory of photon Bloch oscillations in photonic crystals," Phys. Rev. B 63, 035108 (2001).

Phys. Rev. E

G. A. Alfimov, P. G. Kevrekidis, V. V. Konotop, and M. Salerno, "Wannier functions analysis of the nonlinear Schrdinger equation with a periodic potential," Phys. Rev. E 66, 046608 (2002).

Phys. Rev. E.

V. S. Shchesnovich, S. B. Cavalcanti, J. M. Hickmann, and Yu. S. Kivshar, "Zener tunneling in two-dimensional photonic lattices," Phys. Rev. E. 74, 056602 (2006).

Phys. Rev. Lett.

L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, H. Fehrmann, and M. Lewenstein, "Atomic Quantum Gases in Kagome Lattices," Phys. Rev. Lett. 93, 030601 (2004). [CrossRef] [PubMed]
A. A. Burkov and E. Demler, "Vortex-Peierls States in Optical Lattices," Phys. Rev. Lett. 96, 180406 (2006). [CrossRef] [PubMed]
V. A. Brazhnyi, V. V. Konotop, and V. Kuzmiak, "Nature of the Intrinsic Relation between Bloch-Band Tunneling and Modulational Instability," Phys. Rev. Lett. 96, 150402 (2006). [CrossRef] [PubMed]
V. Agarwal, J. A. del Rıo, G. Malpuech, M. Zamfirescu, A. Kavokin, D. Coquillat, D. Scalbert, M. Vladimirova, and B. Gil, "Photon Bloch Oscillations in Porous Silicon Optical Superlattices," Phys. Rev. Lett. 92, 097401 (2004). [CrossRef] [PubMed]
H. Trompeter, T. Pertsch, F. Lederer, D. Michaelis, U. Streppel, A. Brauer, and U. Peschel, "Visual Observation of Zener Tunneling," Phys. Rev. Lett. 96, 023901 (2006). [CrossRef] [PubMed]
H. Trompeter, W. Krolikowski, D. N. Neshev, A. S. Desyatnikov, A. A. Sukhorukov, Yu. S. Kivshar, T. Pertsch, U. Peschel, and F. Lederer, "Bloch Oscillations and Zener Tunneling in Two-Dimensional Photonic Lattices," Phys. Rev. Lett. 96, 053903 (2006). [CrossRef] [PubMed]
M. B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, "Bloch Oscillations of Atoms in an Optical Potential," Phys. Rev. Lett. 76, 4508-4511 (1996). [CrossRef] [PubMed]
B. Bourlon, D. C. Glattli, B. Placais, J. M. Berroir, C. Miko, L. Forro, and A. Bachtold, "Geometrical Dependence of High-Bias Current in Multiwalled Carbon Nanotubes," Phys. Rev. Lett. 92, 026804 (2004). [CrossRef] [PubMed]
A. Sibille, J. F. Palmier, and F. Laruelle, "Zener Interminiband Resonant Breakdown in Superlattices," Phys. Rev. Lett. 80, 4506-4509 (1998). [CrossRef]
M. Jona-Lasinio, O. Morsch, M. Cristiani, N. Malossi, J. H. M¨uller, E. Courtade, M. Anderlini, and E. Arimondo, "Asymmetric Landau-Zener Tunneling in a Periodic Potential," Phys. Rev. Lett. 91, 230406 (2003). [CrossRef] [PubMed]

Phys. Z. Sowjetunion

L. D. Landau, "On the theory of transfer of energy at collisions II," Phys. Z. Sowjetunion 2, 46 (1932).

Proc. R. Soc. London A

C. Zener, "Non-Adiabatic Crossing of Energy Levels," Proc. R. Soc. London A 137, 696-702 (1932).

Progress in Optics

A. S. Desyatnikov, Yu. S. Kivshar, and L. Torner, "Optical Vortices and Vortex Solitons," Progress in Optics 47, 291-391 (2005). [CrossRef]
M. S. Soskin and M. V. Vasnetsov, "Singular Optics," Progress in Optics 42, 219 (2001). [CrossRef]

Rev. Mod. Phys.

L. Esaki, "Long journey into tunneling," Rev. Mod. Phys. 46, 237-244 (1974). [CrossRef]

Science

B. P. Anderson and M. A. Kasevich, "Macroscopic Quantum Interference from Atomic Tunnel Arrays," Science 282, 1686 (1998). [CrossRef] [PubMed]
F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, "Josephson Junction Arrays with Bose-Einstein Condensates," Science 293, 843 (2001). [CrossRef] [PubMed]

Z. Phys.

F. Bloch, "Uber die Quantenmechanik der Elektronen in Kristallgittern," Z. Phys. 52, 555-600 (1928).

Other

L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford: Clarendon Press, 2003).
N. W. Ashcroft and N. D. Mermin, Solid state physics (NewYork: Holt, Rinehart and Winston, 1976).
Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (San Diego: Academic, 2003).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1.	Fig. 2.	Fig. 3.


Fig. 4.	Fig. 5.	Fig. 6.


Fig. 7.	Fig. 8.

Multimedia

Multimedia Files	Recommended Software
» Media 1: AVI (2621 KB)	QuickTime
» Media 2: AVI (841 KB)	QuickTime
» Media 3: AVI (1716 KB)	QuickTime
» Media 4: AVI (1722 KB)	QuickTime
» Media 5: AVI (2124 KB)	QuickTime
» Media 6: AVI (3015 KB)	QuickTime

« Previous Article | Next Article »

OSA is a member of CrossRef.

Select a Conference
Session or Paper #
Year

OpticsInfoBase

Optics Express

Interband resonant transitions in two-dimensional hexagonal lattices: Rabi oscillations, Zener tunnelling, and tunnelling of phase dislocations

Article Tools

Article Navigation

Abstract

1. Introduction

2. Shallow-lattice approximation

3. Rabi oscillations

4. Zener tunnelling

4.1. Quasi-one-dimensional Zener transitions

4.2. Zener transitions between three Bloch bands

4.3. Zener transitions between six Bloch bands

5. Vortex tunnelling

6. Conclusions

Acknowledgments

References and links

References

J. Phys. A: Math. Gen.

J. Phys. B: At. Mol. Opt. Phys.

Laser Phys.

New J. Phys.

Nuovo Cimento

Opt. Express

Opt. Lett.

Phys. Rev.

Phys. Rev. A

Phys. Rev. A.

Phys. Rev. B

Phys. Rev. E

Phys. Rev. E.

Phys. Rev. Lett.

Phys. Z. Sowjetunion

Proc. R. Soc. London A

Progress in Optics

Rev. Mod. Phys.

Science

Z. Phys.

Other

Cited By

Figures

Multimedia