Abstract
We study radiationmatter interaction in a system of ultracold atoms trapped in an optical lattice in a Mott insulator phase. We develop a fully general quantum model, and we perform calculations for a onedimensional geometry at normal incidence. Both two and threelevel atomic configurations are studied. The polariton dispersion and the reflectivity spectra are characterized in the different regimes, for both semiinfinite and finitesize geometries. We apply this model to propose a photon energy lifter experiment: a device which is able to shift the carrier frequency of a slowly travelling wavepacket without affecting the pulse shape nor its coherence.
Light propagation in atomic Mott Insulators
Francesco Bariani
BECCNRINFM and Dipartimento di Fisica, Università di Trento,
Via Sommarive 14, I38050 Povo, Italy
bariani@science.unitn.it
Iacopo Carusotto
BECCNRINFM and Dipartimento di Fisica, Università di Trento,
Via Sommarive 14, I38050 Povo, Italy
1 Introduction
It is a general fact of electrodynamics in continuous media that matter acts as a sort of potential for light: the dielectric constant in Maxwell equations plays in fact a role mathematically equivalent to the external potential in the Schrödinger equation [1]. In particular, remarkable features appear for light propagation through a system with a periodic modulation of the dielectric properties on the scale of the radiation wavelength [2, 3]: in analogy with electron propagation in crystalline solids, these structures are then called photonic crystals (PCs). According to the FloquetBloch theorem, the discrete translational symmetry garantees in fact the conservation of the Bloch wave vector in the first Brillouin zone, the spatial periodicity of the eigenfunctions in the elementary cell, and the separation of the energy spectrum in bands and forbidden gaps [4, 5]. Depending on the frequency dependence of the dieletric response of the material, PCs are generally classified into two broad categories: passive, if the underlying media are dispersionless, and resonant [6, 7].
An interesting example of resonant PC can be realized using ultracold alkali atoms trapped in the periodic potential of an optical lattice [8, 9]. In particular, a milestone in this field has consisted of the achievement of a Mott Insulator (MI) phase [10, 11], in which an integer and constant number of atoms are trapped at each site of the periodic optical potential formed by the interference pattern of several laser beams. If temperature is low enough, atoms are frozen in the lowest vibrational level of each well and the system periodicity is extremely regular: almost no impurities are in fact present (i.e. missing or extra atoms), and there are no phonons. Thanks to the simple level structure of alkali atoms, one can selectively address specific transitions, so to realize e.g. two or threelevel models. From this point of view, the MI can be seen as an extremely resonant PC, where the Bragg scattering processes due to the periodic arrangement of atoms have a strong interplay with the atomic optical resonances. The cleanness of the system, and the weakness of spurious nonradiative effects guarantees that coherence can be preserved for very long times during radiationmatter interaction. This is crucial to manipulate the photon propagation without losing its coherence.
In this paper, we develop a general and fully quantum theory to describe the radiationmatter interaction in these systems. Although the model is perfectly general, we specialize our analysis to the 1D case with normally incident light. Both a twolevel and a threelevel configurations are considered. The band dispersion is characterized as a function of the relative position of the atomic resonance and the Bragg frequency corresponding to the lattice periodicity. Two main regimes are identified depending on whether the two frequency scales are close or well separated. The band dispersion is then used in a calculation of the reflection spectra for both semiinfinite and finite systems. These results are the starting point to propose a photon energy lifter experiment [12], where the carrier frequency of a slowly travelling wavepacket can be continuously tuned while fully preserving its pulse shape and its coherence.
The paper is organized as follows. In Sec.2, we review the Hopfield approach to radiationmatter interaction and its application to the case of a lattice of twolevel atoms. In Sec.3, the Hopfield formalism is applied to get predictions for the dispersion of the elementary excitations of the system, the socalled polaritons. Reflection spectra for semiinfinite and finite 1D systems are presented. In Sec.4 we extend the model to the case of threelevel atoms in a level configuration: we implement the dressed atom approach into the Hopfield Hamiltonian and we obtain predictions for the polariton dispersion and reflection spectra. In Sec.5, the application of threelevel atomic systems as photon energy lifter is discussed for realistic parameters of Rb systems. The main experimental issues are discussed in Sec.6. Conclusions are drawn in Sec.7.
2 Hopfield approach
We consider the interaction between light and a collection of twolevel atoms trapped in the periodic potential of an optical lattice [8, 9]. The atoms are assumed to be in a perfect Mott Insulator state where one and only one atom is present at each lattice site [10, 11]. With a suitable choice of the lattice laser frequency and polarization, the minima of the optical potential felt by the different atomic states can be made to coincide and to form a cubic lattice of spacing . Provided trapping at each site is tight enough, the atomic centerofmass motion results frozen in the motional ground state around the potential minimum (LambDicke regime) and no atomic tunneling is allowed between adjacent sites. The energies of the ground and excited states of the atom have therefore to include the light shift due to the optical potential as well as the zeropoint kinetic energy: the resulting energy separation is indicated by .
Singleatom excitations on the atom are respectively destroyed and created by the operators and . In a solidstate terminology, these localized excitations can be seen as an extreme kind of Frenkel excitons where the electronic excitation is confined to a single atom or molecule of the crystal. This, in contrast to weakly bound Wannier excitons whose electron and hole can be spatially separated by a distance much larger than the lattice spacing [4].
In addition to the coupling to the transverse e.m. field that we shall discuss at length in what follows, excitons are coupled to the longitudinal e.m. field. This Coulomb interaction between the dipoles of distinct atoms allows for the transfer of excitation from one site to another according to Hamiltonian terms of the form [14] . To take the most advantage of the translational symmetry of the lattice, one can construct creation operators for delocalized excitations with a welldefined Bloch wave vector :
(1) 
An integration box of size is assumed with periodic boundary conditions. is the total number of atoms in the lattice (which is assumed to fill the whole integration box). The exciton operators satisfy the following approximate Bosonic commutation rules
(2) 
where is the number of excitations present in the system. Excitons therefore behave as bosons, at least in the “linear optics” limit in which the probability of a double excitation of the same atom is negligible. Throughout the whole paper, we shall stick to this limit.
Thanks to the translational invariance of the system, the Hamiltonian describing the internal atomc dynamics is diagonal in the operators:
(3) 
The dependence of the exciton band is a consequence of the Coulomb dipoledipole interaction and describes the exciton propagation [21]. As this dependence is quite weak in the present system, we will for simplicity neglect it in the following and take a constant value throughout the whole first Brillouin zone (fBz) of the reciprocal lattice. In physical terms, the difference between and corresponds to the ClausiusMossotti local field correction [14].
The quantized transverse electromagnetic field is represented by the e.m. vector potential operator [13]
(4) 
where and are the photon annihilation and creation operators for the different modes, labelled by their polarization state and their wave vector. This latter sum is split into the sum over Bloch wave vectors spanning over the fBz of the reciprocal lattice, and vectors belonging to the reciprocal lattice. In the cubic lattice geometry under consideration here, the reciprocal lattice is itself cubic, with a lattice constant equal to [4]. The freefield Hamiltonian has the usual form
(5) 
with the vacuum frequency of the photon . is here the speed of light and the dielectric constant of vacuum. Throughout the whole paper SI units are in use.
In addition to the terms describing the internal atomic dynamics and the noninteracting field , the total Hamiltonian has to include terms that couple the matter to the transverse e.m. field. This is most simply done by means of the standard minimal coupling replacement [13]. As we are considering optical fields with wavelengths much bigger than the atomic size, a dipolar approximation can be performed within . Atoms are represented as electric dipoles of dipole moment . For the sake of simplicity, we restrict to a single polarization state parallel to the dipole moment of the atoms. In this way, we can drop the polarization index and write in the final, compact form:
(6) 
We can easily recognize the different interacting terms: the former is the usual dipole excitonphoton coupling, with strength proportional to
(7) 
and the latter is a “diamagnetic” photonphoton coupling induced by the presence of atoms, with a coupling strength proportional to
(8) 
For typical values of the system parameters, the coupling coefficient is generally small, i.e. , which implies .
Thanks to the discrete translational symmetry of the lattice, Bloch momentum is conserved by all terms in (6). This means that only states with the same Bloch wave vector are mixed by the lightmatter interaction. As this set is discrete, radiative decay (i.e. spontaneous emission) is forbidden and energy coherently oscillates between the photonic and the excitonic degrees of freedom [14]. This observation is crucial in simplifying the physical analysis of the system.
The normal modes of the system, the socalled polaritons are therefore superposition of the original photonic and excitonic modes; they are classified by their Bloch momentum , and by a band index . Thanks to the quadratic structure of the Hamiltonian, the polaritonic operators can be obtained by means of a HopfieldBogoliubov transformation in the general form [14]:
(9) 
Provided the correct normalization is chosen, the polaritonic operators satisfy bosonic commutation rules
(10) 
The coefficients are determined by solving the eigenvalue problem associated with the commutator
(11) 
where is the polariton energy. The solution of this system is equivalent to the diagonalization of the following Bogoliubov matrix
(12) 
with the Hermitian matrix
(13) 
and the symmetric matrix
(14) 
The first row and column correspond to the matter excitation. For notational simplicity, only two photonic modes of wave vectors and have been explicitely shown here, but the matrices are intended to contain rows and columns for each reciprocal lattice vector .
While the diagonal blocks of are hermitian, the nondiagonal ones break the hermiticity of the matrix in the usual sense: the matrix is in fact hermitian, in the sense that , where defines the Bogoliubov metric. Physically, this property is related to the fact that the blocks correspond to the antiresonant terms in the Hamiltonian (6), which do not conserve the number of excitations [15, 16, 17, 18, 19]. Thanks to the small value of the lightmatter interaction coefficients and in the atomic systems under consideration here, most of the physics under investigation here can be obtained by neglecting the antiresonant terms and truncating the matrix to a small number of photonic modes [20, 21].
3 1D lattice of twolevel atoms
3.1 Band structure
To get a simple physical understanding of the system, it is useful to concentrate our discussion on the simplest case of a 1D geometry: most effects related to resonant lightmatter interaction are in fact indipendent from the dimensionality of the system under consideration [6, 20, 22]. Interesting discussions of the optical properties of different kinds of 1D resonant PCs can be found in [23, 24, 25, 26, 7, 27].
Two frequency scales are to be considered: the atomic resonance frequency , and the Bragg frequency of the lattice which carries information on the periodicity of the lattice and derives from the famous Bragg condition for diffraction spectroscopy in a reflection geometry [4]. The width of the frequency region in which states are effectively mixed and their splittings are determined by the interacting Hamiltonian (6). At leading order, this is proportional to . Starting from this consideration, two different regimes can be distinguished according to the ratio between the difference and the characteristic mixing .
The first regime, let’s call it purely excitonic regime corresponds to the case when the resonance frequency and the Bragg frequency are well separated . An example of polaritonic dispersion for this regime is shown in Fig.1 for . For the sake of clarity, a somehow exaggerated value of has been used in the figure.
On one hand, the exciton and a single photonic mode intersect in the interior of the fBz; their mixing results in an anticrossing of the polariton modes with a Rabi splitting equal to . Far from this region, the polaritonic modes tend to almost purely photon and exciton modes. On the other hand, pair of photon modes intersect at the edge of the fBz at frequencies multiple of ; their mixing is due to Bragg scattering processes on the atomic lattice.
Gaps open at the edges of the fBz: the lower one (around ) is the usual polaritonic gap of resonant dielectrics [28], while the upper one (just above ) is due to Bragg scattering processes. Differently from the bulk case, the former gap extends on both sides of because of the limited size of the fBz. The latter one is instead located strictly above . Its lower edge is exactly at and corresponds to a purely photonic state, unaffected by the presence of the atoms which are located at the electric field nodes. As usual, the polaritonic density of states vanishes inside the gaps, and radiative propagation at these frequencies is forbidden.
From a formal point of view, the lower gap mainly originates from the excitonphoton coupling, while the upper one contains a contribution from the direct photonphoton coupling as well. As , all these terms are of the same order. This picture is preserved in 2D and 3D, as all other photonic modes that may participate are far away in energy.
The other regime, let’s call it mixed excitonBragg regime, is characterized by the condition . In this case three modes are simultaneously mixed, namely two photon branches (the incoming one at and the first Bragg diffracted at ) and the excitonic state. Differently from the previous case, the splittings (again of the order of ) are now located close to the edges of the fBz. As one can see in Fig. 1, this results in much wider forbidden gaps of the order of the splitting. In addition to this, the “central” band between and results squeezed and shows a very flat dispersion over most of the fBz.
In the mixing region around the band edge, the upper and lower polaritons are mixture of exciton and photons with almost equal weights, while the central polariton band at is mostly photonic. This effect is easily explained in terms of the electric field showing nodes at the atomic locations, as in the previous regime.
3.2 Reflectivity spectra
The band dispersion introduced in the previous section is a complete description of the photon propagation inside the lattice. Most spectroscopic experiments, however, involve light beams which are incident onto finite systems and therefore require a description of the interfaces between regions of different optical properties, namely the external vacuum and the atomic lattice. This allows to calculate crucial properties of the system, such as its reflectivity [7, 26, 25, 24, 23].
This can be done by imposing suitable boundary conditions on the electromagnetic fields: Maxwell theory imposes in fact the continuity of both the electric field and its spatial derivative. In the following, we shall restrict to the most relevant case of a coherent, monochromatic excitation at . Outside the lattice, we have purely photonic, coherent plane waves with wave vectors such that . Inside the lattice, the field propagates as coherent polaritonic states of modes such that .
The electric field is the expectation value of the electric field operator
(15) 
Outside the lattice, the electric field associated with a photonic coherent state of wave vector is given by
(16) 
Inside the lattice, one is to consider a polaritonic coherent state with Bloch wave vector in the band satisfying
(17) 
time evolution of such a state under the system Hamiltonian corresponds to
(18) 
The electric field of such a state is then:
(19) 
The amplitudes and are then determined by matching the fields (16) and (19) at the interfaces for the given geometry under consideration.
The plots in Fig. 2 and Fig. 3 show the reflectivity spectra in respectively the purely excitonic and the mixed excitonBragg regimes. Two geometries will be considered: a semiinfinite [panels (a)], and a finite one [panels (b)]. Note that these predictions exactly coincide with the ones of semiclassical calculations where matter is described in terms of a dielectric polarizability [6, 24, 22, 7, 27, 25, 26].
3.2.1 Semiinfinite geometry
In the first case, Fig. 2 and Fig. 3, there is a single interface, dividing the space in two semiinfinite regions: vacuum and lattice.
We first consider the input problem with an incoming and a reflected plane wave in the vacuum and a single transmitted polariton Bloch mode in the lattice: the Bloch momentum in the fBz and the band index of this latter are chosen in order to satisfy energy conservation . The reflection amplitude is then expressed in terms of the electric field of the polariton and its spatial derivative as
(20) 
The point lies on the edge of the elementary cell, i.e. midway between neighbouring atoms ^{1}^{1}1Actual numerical calculations are performed by truncating the full matrix (12) to a finite number of modes; in order to obtain a smooth convergence, a gaussian cutoff has also been added on higher photonic modes. This cutoff physically mimicks finite size atoms. Because of the singularity in the field at the atomic position , evaluating the fields at rather than at ensures faster convergence. .
The mismatch between the photonic components of the polariton state and the incoming wave determines the reflectivity shown in the figures: this is significant around the gaps where the incoming wave is strongly mixed with either an exciton [around the lower gap in Fig.2], a photon [around the upper gap in Fig.2], or both [around both gaps of Fig.3]. In the mixed excitonBragg regime, note that the reflectivity remains quite large in between the two gaps: the flatter the middlepolariton branch, the higher the corresponding reflectivity. In the present semiinfinite geometry, reflectivity is complete inside the gaps where the wave vector becomes imaginary and the field inside the lattice consists of an evanescent wave.
The output problem corresponds to two counterpropagating Bloch modes with the same energy inside the lattice and a single transmitted plane wave in the external vacuum. In this situation, the reflectivity is given by with
(21) 
where the and signs refer to the propagation versa of the polaritons.
The system being invariant under time reversal and spatial parity, it is easy to prove that the coefficients and entering in the formula (19) for polariton field share the same phase and can be chosen to be all real. The electric fields of counterpropagating polaritons are therefore complex conjugates of each other . Plugging this into (21) and comparing the result with (20), it is immediate to see that .
3.2.2 Finite slab
Reflectivity spectra are shown in Figs. 2 and 3 for a finite system containing a quite large number of elementary cells. Note that the present approach is not able to calculate the reflectivity inside the gaps where the field consists of an evanescent wave. The main difference with respect to the semiinfinite case is the presence of fast oscillations on top of the reflectivity spectrum around the main gaps. This can be explained as follows.
Two interfaces at respectively and now separate three regions of space: the vacuum with the incident and reflected photons, the finitesize lattice with counterpropagating polaritons, and again vacuum with now only a transmitted photon. The field in the last cell () is determined by the output problem to be
(22) 
As both are Bloch states, the field in the first cell (taking ) has the simple form
(23) 
By solving the continuity conditions at the front interface at , we get
(24) 
Because of the phase factors in (23), fast oscillations occur in the reflectivity (24) due to the FabryPerotlike interference of Bloch waves which undergo multiple reflections at the lattice boundaries. The period of these oscillations is fixed by the group velocity and the total length of the system ,
(25) 
the slower , the closer the peaks.
To compare the envelope of this oscillations with the spectrum in the semiinfinite geometry, we can consider a simplified model where the lattice is replaced by a bulk medium of refractive index . In this case, the reflectivity for a single interface separating vacuum and medium is
(26) 
For a slab of thickness , the reflectivity is [29]
(27) 
FabryPerot oscillations are apparent, with a maximum reflectivity at the peaks equal to
(28) 
In the limit , the ratio : this is due to the presence of two counterpropagating Bloch modes in the slab as compared to the single propagating mode in the semiinfinite case. This factor provides a good approximation in the lattice case as well, as one can easily see in the lowreflectivity tails of the spectra of Figs.2 and 3.
The peak at Bragg frequency, which is more related to the photonic components, is generally smaller than the resonant one. This difference is dramatically enhanced in the purely excitonic regime as you can see in Fig.2.
4 1D lattice of threelevel atoms
Systems of threelevel atoms are of great relevance for both quantum and nonlinear optics. A most remarkable phenomenon in this respect is the so called Electromagnetically Induced Transparency (EIT) [30] which, among other properties, can lead to ultraslow light propagation [31, 32, 33] in spectral regions where absorption and reflectivity are also very low. This fact makes threelevel systems extremely interesting systems to study and manipulate light in new regimes.
We consider the threelevel configuration shown in Fig.4: in addition to the previously considered twolevel scheme, there is a metastable state which is coupled to the excited state by a classical (laser) field. The transition is optically forbidden.
The coupling beam dresses the metastable state, so to give a new mixed excitation formed by the atom promoted to the metastable state and an extra photon correspondingly emitted into the dressing beam: this new state has an energy , where and are the frequencies of respectively the bare metastable state (with respect to the ground state ) and the dressing laser light [13]. Its detuning from the excited state is then . The Hamiltonian describing the dressing is
(29) 
where is the Rabi frequency of the dressing beam and the operators , destroy and create the dressed metastable excitation for the atom .
As above, the introduction of the radiationmatter interaction is performed by means of the minimal coupling replacement [19]:
(30) 
The operators remain instead unchanged as the metastable state is not directly coupled to the quantized field. After the minimal coupling replacement, the dressing Hamiltonian has the form:
(31) 
Note the appearance of terms directly coupling the dressed metastable state to the light field. These terms are however of the order of (with ), and therefore much smaller than all the other coupling strengths in the system.
By following the same procedure as above, new excitonic delocalized operators can be constructed for the dressed state, and then included in the Bogoliubov matrix: the resulting polaritonic bands and the corresponding reflectivity spectra are shown in Figs. 5(a,b) in the purely excitonic regime and for a resonant dressing . Note that only the region around the atomic resonance is shown in the figure.
For the sake of simplicity, we concentrate our attention onto a purely excitonic case where the lattice spacing is much smaller than the resonance optical wavelength. The signature of the threelevel scheme is the presence of a flat EIT miniband corresponding to the “central” polariton. The most interesting region is around the Raman resonance between the quantized radiation and the dressed state. This spectral region can be simply explored in terms of the reduced Bogoliubov matrix
(32) 
where only the incident photon (1st row) mode, the excited state (2nd row), and the dressed metastable state (3rd row) are considered. is the detuning of the photon mode from atomic resonance. As we are considering a small region , we have replaced with in the lightmatter coupling terms. We have also neglected the and coupling terms as they are much smaller than the other terms.
A noticeable property of (32) exactly at Raman resonance is that for any value of different parameters there exists an eigenvalue corresponding to the central polariton. Its eigenvector has a vanishing exciton component and a photonic weight
(33) 
so that the group velocity is
(34) 
a significant reduction of the polariton group velocity is then observed as soon as a weak dressing amplitude is used (see Fig.5(c)). Furthermore, reflection at the lattice interfaces is vanishing in the central region of the central polariton band, which allows for easy injection of polaritons in the system (as shown in Fig.5(b)).
A simple model to quantify the width of the region where reflection is low can be developed as follows. Defining an effective index for the lattice as , and applying classical reflection result (28), the reflectivity envelope for a semiinfinite system turns out to be approximately given by:
(35) 
This approximation is accurate at the center of the dip where the reflectivity is small. Comparing the two main results (34) and (35) of the present section, a tradeoff is apparent: the slower the light, the narrower the reflectivity dip [32].
5 Application: photon energy lifter
Obtaining a coherent and widely tunable frequency conversion of an optical signal is a central task in optical telecommunications [34]. Several techniques have been developed during the years to perform this operation, but most of them suffer from significant limitations in their application range, or are disturbed by spurious effects.
A recent proposal is based on the socalled dynamic photonic structures (for a review see e.g. [35]), whose optical properties are varied in real time while the optical wavepacket is propagating inside them. The basic idea of the photon lifter consists in the adiabatic shift of the photonic band on which the photon is located. This was orginally proposed for solidstate photonic structures [12], but it is interesting to explore the potential of cold atom systems to this purpose: the very long coherence times in a Mott insulator state and the easy tunability by external electric or magnetic fields makes them very promising for this kind of applications.
As a specific example, we shall consider in what follows a system of Rb atoms, trapped in a cubic optical lattice of spacing nm. The optical properties are varied by means of an external magnetic field (taken as parallel to the axis) which shifts the atomic levels via the Zeeman effect [36].
In the following, we concentrate on the transition at a frequency [37].
As we are interested in substantial shifts, we can concentrate our attention in the high field regime () where the atomic nucleus is decoupled from the electronic degrees of freedom, and the energy shift mostly comes from the electronic total angular momentum only: , where is the Bohr magneton, is the Landé factor of the considered level and is the component of the total angular momentum of the electron.
We use the sublevels of the electronic ground state as respectively ground and metastable states, and the sublevel of the electronic excited state as excited state. The corresponding Landé factors are and .
The nucleus is not affected by the optical process and maintains the same polarization it had in the initial state: in the absolute atomic ground state, the nuclear spin is e.g. polarized antiparallel to the electron spin of the state.
A polarization is used for the dressing light beam that couples the and the states, and a one is used to probe the polariton dispersion on the transition. Using tabulated values for the electric dipole moment of the transition, the excitonphoton coupling (7) for the system under consideration is of the order of .
To maximize the available time to perform the lifter operation, it is useful to have a very slow group velocity, which in turns requires a small dressing amplitude. In the following, we shall choose . This value corresponds to times the radiative linewidth of the line of Rb atoms.
The dressing frequency is chosen in a way to have at the initial value of the magnetic field: the corresponding polariton dispersion is the one shown in Fig.5(a). The light pulse is injected into the system in proximity of the resonant point , where the interface reflectivity goes to zero, and injection is most effective (see the circle in Fig.5(a)): the width of this dip results from (35) to be of the order of and the group velocity (34) is , i.e. .
The magnetic field variation is performed while the light pulse to be shifted is completely contained in the lattice and is propagating through an effectively bulk system. As the magnetic field is varied in a spatially homogeneous way, the Bloch wave vector is conserved during the process. If the field variation is slow enough as compared to the frequency difference of neighbouring bands, the polaritons will adiabatically follow the band and their frequency at the end of the process will be accordingly shifted (see the circle in Fig.5(d)).
As an example, we propose to tune the magnetic field from up to T: this results in the metastable and excited states being shifted by respectively and with respect to the ground state. For light initially injected in proximity of , the shift of the photon frequency results approximately equal to , which amounts to the quite sizeable value . As the lifter operation is based on an adiabatic shift of the polariton dispersion, it completely preserves the pulse shape and the coherence properties of the incident wavepacket, both at classical and at quantum level.
6 Experimental issues
To verify the actual feasibility of such a promising experiment, it is important to mention the main practical difficulties that may arise in an actual experiment, and discuss how these can be overcome.

We have verified that the transmittivity of the lattice interfaces is close to for both the injection and the extraction process (Figs.5(b,e)). The pulse is injected into the lattice at a frequency corresponding to the EIT reflectivity dip around Raman resonance. The extraction takes place in close proximity of the Raman resonance where reflectivity is again very low. This, in spite of the fact we are very close to a gap: thanks to the now significant detuning , the metastable state is in fact weakly coupled to light, and the corresponding crossing point is displaced slightly away from the light line.

In order to have a reasonably long time to vary the magnetic field, we have verified that the group velocity of the polariton states involved in the lifter operation is slow. Light initially propagates on the EIT slow light branch, which is deformed during the lifter operation. At the end, the wavepacket is found on the very flat region below the gap where the group velocity is low.

The wavepacket has to be shorter than the lattice length, still its frequency spectrum has to fit in the reflectivity dip at both injection and extraction. A lattice of cells is able to accomodate pulses with at most . From panel (a), this corresponds to a lower bound on the frequency width of the incoming wavepacket . One can easily see in panel (b) that this frequency spread still fits within the injection window where reflectivity is low. The same can be verified on panels (de) for the extraction process.

In order for the pulse shape not to be affected, dispersion of the group velocity should be small for the wavevector window under examination. Initially, this is not a problem, as we are working close to the center of the EIT branch where the group velocity has a weak dispersion. The situation can be more critical on extraction, because of the strong squeezing of the polariton band in the region just below the gap. The importance of this effect can be reduced by choosing pulses initially tuned just above the Raman resonance.
One major constraint that still exists on the experimental parameters concerns the speed at which the magnetic field has to be actually varied. As this has to be done while the wavepacket is inside the lattice, a very slow group velocity and a long lattice are required. Using values for stateoftheart MIs. namely , and , one obtains that one disposes of a time of approximately to perform the magnetic field variation. This means that a variation of requires a very large rate of .
As this can pose serious difficulties in an actual experiment, it is worth briefly exploring alternative strategies. An interesting possibility is to further reduce the dressing amplitude . As the polariton group velocity is proportional to the square of the dressing amplitude, the value used in a recent slow light experiment [33] already leads to a group velocity of the order of which corresponds to an available time of . In the highfield regime considered here, a photon frequency shift of then requires a magnetic field variation of in , a rate routinely used in cold atom experiments.
It is important to note that the reduction of the dressing amplitude implies a squeezing of the reflectivity dip at injection and an enhancement of the dispersion at extraction. These, apparently serious problems are overcome thanks to the fact that a reduction in the group velocity implies a spatial shortening of the pulse in the lattice, and therefore a reduced frequency spread for a given length.
7 Conclusions
In conclusion, we have developed a fully quantum description of radiationmatter interaction in a gas of ultracold atoms trapped in the Mott insulator phase of an optical lattice. The coherent interaction between photons and atomic excitations gives rise to new, mixed polaritonic excitations.
In the case of twolevel atoms, two different regimes are identified. In the purely excitonic regime, where the atomic resonance is far from the Bragg frequency of the lattice, two gaps appear in the energy spectrum. In the mixed excitonBragg regime, the interplay of the atomic resonance and the lattice periodicity enhances the gap amplitudes and gives rise to a flat band between them. The consequences of the polariton dispersion on the reflection properties of finite lattices have been investigated.
The theory is then extended to a system of threelevel atoms in a configuration. The slowlight band which joins weak reflection to ultraslow group velocity turns out to be the most promising in view of applications as a photon lifter, i.e. a device able to shift the carrier frequency of an optical pulse without affecting its shape nor spoiling its coherence properties. Advantages and disadvantages of using atomic gases as photon lifters are discussed.
Future work will address the application of slow as well as stopped polaritons as nondestructive probes of the microscopic properties of ultracold gases, e.g. the behaviour of matter wave coherence across the superfluid to Mott insulator transition.
8 Acknowledgements
We thank A. Recati, G. C. La Rocca, V. Guarrera and S. De Liberato for stimulating discussions. FB acknowledges hospitality at the Institut Henri PoincareCentre Emile Borel and financial support from ESFQUDEDIS network through Short Visit Grant, Ref. Num. 1802. IC acknowledges hospitality at the Institut Henri PoincareCentre Emile Borel and financial support from CNRS. FB dedicates his first paper in Physics to his family for the continuous and strong support received.
References
 [1] P. S. Russell, T. Birks, and F. Dominic LloydLucas, “Photonic Bloch waves and photonic band gap”, in Confined electrons and photons (Plenum Press, 1995), pp. 585–633.
 [2] E. Yablonovitch, “Inhibited spontaneous emission in solidstate physics and electronics”, Phys. Rev. Lett. 58 2059 (1987).
 [3] S. John, “Strong localization of photons in certain disordered dielectric superlattices”, Phys. Rev. Lett. 58 2486 (1987).
 [4] N. Ashcroft and N. Mermin, Solid state physics (Brooks/Cole, 1976).
 [5] P. S. Russell, “Interference of integrated FloquetBloch waves”, Phys. Rev. A 33 3232 (1986).
 [6] D. Coevorden, R. Sprik, A. Tip, and A. Lagendijk, “Photonic band structure of atomic lattices”, Phys. Rev. Lett. 77 2412 (1996).
 [7] M. Erementchouk, L. Deych, and A. Lisyansky, “Optical properties of onedimensional photonic crystals based on multiplequantumwell structures”, Phys. Rev. B 71 235335 (2005).
 [8] C. Adams and E. Riis, “Laser Cooling and Trapping of Neutral Atoms”, Prog. Quant. Elet. 21 1 (1997).
 [9] O. Morsch and M. Oberthaler, “Dynamics of BoseEinstein condensates in optical lattices”, Rev. Mod. Phys. 78 179 (2006).
 [10] D. Jaksch, C. Bruder, J. Cirac, C. Gardiner, and P. Zoller, “Cold bosonic atoms in optical lattices”, Phys. Rev. Lett. 81 3108 (1998).
 [11] M. Greiner, O. Mandel, T. Esslinger, T. Hansch, and I. Bloch, “Quantum phase transition from a supefluid to a Mott insulator in a gas of ultracold atoms”, Nature 415 39 (2002).
 [12] Z. Gaburro, M. Ghulinyan, F. Riboli, L. Pavesi, A. Recati, and I. Carusotto, “Photon energy lifter”, Opt. Express 14 7270 (2006).
 [13] C. CohenTannoudji, J. DupontRoc, and G. Grynberg, Processus d’interaction entre photon et atomes (CNRS Editions, 1996).
 [14] J. Hopfield, “Theory of contribution of excitons to the complex dielectric constant of crystals”, Phys. Rev. 112 1555 (1958).
 [15] M. Artoni and J. Birman, “Quantumoptical properties of polariton waves”, Phys. Rev. B 44 3736 (1991).
 [16] C. Ciuti, G. Bastard, and I. Carusotto, “Quantum vacuum properties of the intersubband cavity polariton field”, Phys. Rev. B 72 115303 (2005).
 [17] C. Ciuti and I. Carusotto, “Inputoutput theory of cavities in the ultrastrong coupling regime: The case of timeindependent cavity parameters”, Phys. Rev. A 74 033811 (2006).
 [18] S. De Liberato, C. Ciuti, and I. Carusotto, “Quantum Vacuum Radiation Spectra form a Semiconductor Microcavity with a TimeModulated Vacuum Rabi Frequency”, Phys. Rev. Lett. 98 103602 (2007).
 [19] I. Carusotto, M. Antezza, F. Bariani, S. De Liberato, and C. Ciuti, “Atomic insulators as novel optical media: from slow light to the dynamical Casimir effect”, (to be published).
 [20] Y. Chong, D. Pritchard, and M. Soljačić, “Quantum theory of a resonant photonic crystal”, Phys. Rev. B 75 235124 (2007).
 [21] H. Zoubi and H. Ritsch, “Excitons and cavity polaritons for ultracold atoms in an optical lattice”, Phys. Rev. A 76 013817 (2007).
 [22] K. Kempa, R. Ruppin, and J. Pendry, “Electromagnetic response of a pointdipole crystal”, Phys. Rev. B 72 205103 (2005).
 [23] T. Ikawa and K. Cho, “Fate of the superradiant mode in a resonant Bragg reflector”, Phys. Rev. B 66 085338 (2002).
 [24] I. Deutsch, R. Spreeuw, S. Rolston, and W. Phillips, “Photonic band gaps in optical lattices”, Phys. Rev. A 52 1394 (1995).
 [25] M. Artoni, G. La Rocca, and F. Bassani, “Resonantly absorbing onedimensional photonic crystals”, Phys. Rev. E 72 046604 (2005).
 [26] M. Artoni and G. La Rocca, “Optically Tunable Photonic Stop Bands in Homogeneous Absorbing Media”, Phys. Rev. Lett. 96 073905 (2006).
 [27] M. Erementchouk, L. Deych, and A. Lisyansky, “Spectral properties of exciton polaritons in onedimensional resonant photonic crystals”, Phys. Rev. B 73 115321 (2006).
 [28] J. Jackson, Classical electrodynamics (Wiley, New York, 1998).
 [29] J. Wait, Electromagnetic waves in stratified media (Pergamon, 1970).
 [30] E. Arimondo, “Coherent Population Trapping in Laser Spectroscopy”, Prog. Optics 35 257 (1996).
 [31] S. Harris, “Electromagnetically Induced Transparency”, Phys. Today 50 36 (1997).
 [32] A. Matsko, O. Kocharovskaya, Y. Rostovtsev, G. Welch, A. Zibrov, and M. Scully, “Slow, Ultraslow, Stored and Frozen Light”, Adv. At. Mol. Opt. Phys. 46 191 (2001).
 [33] L. Hau, S. Harris, Z. Dutton, and C. Behroozi, “Light speed reduction to 17 metres per second in an ultracold atomic gas”, Nature 397 594 (1999).
 [34] S. Yoo, ‘‘Wavelength Conversion Technologies for WDM Network Applications”, J. Lightwave Technol. 14 955 (1996).
 [35] M. Yanik and S. Fan, “Dynamic Photonic Structures: Stopping, Storage and Time Reversal of Light”, Studies in Applied Mathematics 115 233 (2005).
 [36] J. Brehm and W. Mullin, Introduction to the structure of matter (Wiley, 1989).
 [37] See, e.g., D. Steck, Alkali D line data, available at the web page: http://steck.us/alkalidata/.