Localization of an impurity particle on a boson Mott insulator background

Localization of an impurity particle on a boson Mott insulator background

Sevilay Sevinçli    R. O. Umucalılar    M. Ö. Oktel oktel@fen.bilkent.edu.tr Department of Physics, Bilkent University, 06800 Ankara, Turkey
July 12, 2019

We investigate the behavior of a single particle hopping on a three dimensional cubic optical lattice in the presence of a Mott insulator of bosons in the same lattice. We calculate the critical interaction strength between the impurity and background bosons, beyond which there is bound state (polaron) formation. We give exact results in the limit of a perfect Mott insulator, where polaron formation is equivalent to impurity localization. We calculate the effects of lattice anisotropy, higher impurity bands, and fluctuations of the Mott insulator on the localization threshold. We argue that our results can be checked experimentally by RF spectroscopy of impurity particles.

03.75.Lm, 71.10.Fd

I Introduction

Recent experiments using cold atoms in optical lattices have displayed remarkable versatility, creating highly controllable, isolated, low temperature environments where ideas about many particle quantum mechanics can be tested greiner (); bloch1 (); ketterle (); esslinger (); phillips (); weiss (). With improvement over the control of system parameters, as well as advancement of novel measurement techniques such as noise correlations altman (), it seems plausible that a great variety of models will be realized in optical lattices. The precision of optical lattice experiments presents another challenge for theory. Measurements of effects that are beyond the simplest descriptions such as mean field theory, are becoming possible bloch2 ().

One of the new classes of quantum models realized by cold gas experiments is the mixture of different species of atoms jin (). Boson-boson mixtures as well as boson-fermion mixtures have been created. It is now also possible to selectively turn on optical lattice potentials for any of the species forming the mixture inguscio ().

The possibility of experimental realization of atomic gas mixtures stimulated a lot of theoretical interest molmer (); albus (); lewenstein (); sacha (); pethick (); kuklov (); gavish (). From the theory point of view, mixtures are appealing as they can display some qualitatively new physics stemming from simple ideas in many body theory. There is possibility of pairing due to mediated interactions pethick (), formation of composite particles similar to molecules lewenstein (), large counterflows of different species or even countersuperfluidity kuklov (). There are also ideas to simulate random potentials using one species as the disorder potential for the others gavish (). Two recent experiments about boson-fermion mixtures have shown that the optical lattice experiments are advanced enough so that all these theoretical ideas may now be tested in the laboratory gunter (); ospelkaus ().

The parameter space for mixtures is very large, with many possible phases and a complicated phase diagram lewenstein (); sacha (). The phase diagram of mixtures changes remarkably depending on the strength and type of optical lattices, number of components, nature of interactions, and the identity of the particles. While more general investigations of this parameter space, which point out novel phases, are very useful, in this paper, we concentrate on a simple limit and provide some exact results. Experimentally, it would be easy to check such results and gain a better understanding of the various effects that may be present in the system.

We investigate the limit where one bosonic species interacts with a single particle of another species. While the identity of the external particle does not matter in what follows, we refer to it as the fermion for convenience. All the results of the paper are valid for a bosonic impurity as well. We assume that both species share the same lattice potential and that the lattice potential is deep enough so that only one band of the lattice is populated. The Hubbard-type Hamiltonian jaksch () for this system can be written as


where and are the bosonic and impurity creation and annihilation operators respectively; is the on-site number operator for the bosons and represents a sum over nearest neighbors. The strength of the tunnelling terms are characterized by hopping matrix elements for bosons and for the fermion. and are the on-site interaction strengths between bosons and between a boson and a fermion respectively.

Depending on the ratio between the the tunnelling strength and the interaction for bosons, the system may either form a superfluid (SF) or a Mott insulator (MI) state fisher (). While the localization problem can be studied in both regimes, we can provide exact results only for the Mott insulator case. In this paper, we concentrate on the problem in the Mott regime as the problem is qualitatively different from that of a superfluid interacting with an impurity. Thus, we consider , give exact results only in the limit and calculate corrections due to finite boson hopping. For a Mott insulator with bosons per site, the chemical potential is constrained to


With these considerations for bosons, the fermion can show two qualitatively different behaviors. It can either behave like a free fermion with its wavefunction stretching throughout the system, or it may create a defect in the Mott insulator and form a bound state with this defect. We calculate the critical interaction strength that separate these two regimes. In Section II, we study the case of a perfect Mott insulator and obtain the phase diagram for localization of the impurity; in Section III, we generalize this exact result to an anisotropic lattice where hopping strength in one direction is different from the other two. We consider the effect of higher impurity bands in Section IV. In Section V, we calculate the change in the fermion hopping strength due to fluctuations in the Mott insulator background. We summarize our results and discuss experimental methods to measure this bound state formation in Section VI.

Ii Localization in a perfect Mott insulator

In this section, we calculate the critical interaction strength for bound state (polaron) formation in the limit that the boson Mott insulator is perfect, i.e. the hopping strength for bosons is zero (). We relax this assumption in Section V, and analyze the effects of fluctuations in the Mott background. When boson hopping is neglected the Mott insulator background becomes almost inert for the fermion, presenting a spatially independent mean field energy shift. In this case the fermion will move with the dispersion relation


where is the crystal momentum of the fermion in the direction and is the lattice constant. We cannot expect this delocalized behavior of the fermion to continue if the interactions between the background bosons and the fermion become very strong. Let us assume that the fermion and the bosons attract each other, i.e. . In this case, it would be energetically favorable to put more bosons at a lattice site and bind the fermion to these bosons. This will only happen at a critical interaction strength, beyond which the energy gained by boson-fermion attraction is larger than the sum of the kinetic energy cost of localizing the fermion and the interaction energy cost of introducing more bosons.

This critical interaction strength can then be found by investigating the single particle Hamiltonian


where we take the site at which the defect is formed as the origin and assume that the defect represents a localized attractive potential to the fermion. We can then ask for which value of there will be a bound state in the spectrum. This is the discrete version of the problem of existence of a bound state for a localized potential well landau (). Just like the continuum version, in one and two dimensional lattices there is a bound state for an infinitesimally small attractive potential. In three dimensions, however, there is a certain critical value below which there is no bound state. In the case of a finite attractive potential, the dispersion relation takes the following implicit form


where and are the scaled quantities and is the Fourier transform of the particle’s wave function. We obtain the relation between the binding energy and the attractive potential by taking the Fourier transform of Eq. (5)


where and we take


with being the binding energy. When , i.e. at the localization threshold, the above integral can be evaluated exactly delves () and the critical value at which the localization takes place is


where and is the complete elliptic integral of the first kind. For nonzero , the integral was evaluated by Joyce joyce (). Using this result, the potential is obtained as


where , , and . In the limit of large binding energy, , we obtain a linear relation , which is expected as the particle is strongly localized at a single lattice site.

The exact evaluation of the integral above allows us to calculate not only the critical boundary but also the binding energy of the bound state (Fig. 1).

Figure 1: Binding energy of the impurity as a function of . The critical interaction strength where the localization begins can also be obtained from the figure, i.e. for .

We now translate our one particle results to the many particle case. Let us first assume that the fermion-boson interaction is attractive, . In this case the simplest defect would be to introduce one more boson, thus the attractive potential seen by the fermion will be . However, to introduce one more boson to a Mott insulator with particles per lattice site would cost energy, . Thus, the phase boundary between the free fermion state and the bound state of the fermion and one bosonic defect (polaron) is given by


This is by no means the only defect that can be created in the Mott insulator. If the boson-fermion attraction is strong enough, it becomes energetically favorable to attract more bosons and form a bound state of two bosons and one fermion. The phase boundary for such a defect can be decided by comparing the energy of this state with the energy of the bound state of one boson and one fermion. Thus, the equation for phase boundary is


One can similarly find the boundaries for bound states with higher number of bosons.

Another kind of defect would present itself for repulsive interactions, . For sufficiently strong repulsive interactions it would be preferable to create a hole in the Mott insulator state and bind the fermion to this hole. The corresponding phase boundary is given by


Similar to the attractive interactions, it is possible to form bound states of the fermion with more holes. One can continue to deplete the Mott state until all the bosons are removed from the defect site. After this point it would be preferable to deplete bosons from the neighboring sites. We have, however, not included such states in our phase diagram. In Figs. 2, 3, and 4, we present three phase diagrams for different chemical potentials. Fig. 2 indicates that when , phase diagram is symmetric around . This is expected as this value of corresponds to lobe centers of the Bose-Hubbard phase diagram where there is particle-hole symmetry. One can also notice from this diagram that when is close to zero, even for small values, the fermion can be bound to a large number of bosons. If we take as in Fig. 3, the symmetry around is broken and for the repulsive interactions it is harder to attract holes. Fig. 4 represents the opposite case () where stronger interactions are required to attract particles.

Figure 2: Phase diagram for . Numbers in each region show how many extra particles () or holes () are attracted to the localization site. The region marked as contains all the phases with seven or more extra bosons (holes). Phase diagram for this value of is symmetric around . For small boson-boson repulsion , even for small values, large number of bosons are attracted. While this phase diagram is independent of , the number of holes that are attracted is limited by .
Figure 3: Phase diagram for . Symmetry in Fig. 2 is broken and particle attraction is easier than the hole attraction, since the chemical potential is increased with respect to the symmetry point. To attract a hole one needs higher boson-fermion interaction for the same .
Figure 4: Phase diagram for . As opposed to Fig. 3, to attract a particle one needs higher boson-fermion interaction.

We believe that the phase diagram can be checked experimentally. While it would be possible to modify by an interspecies Feshbach resonance, an easier route would be to change which is controlled by the strength of the optical lattice. The localized impurity states can be distinguished from free fermion states as their mean field shifts would be different; in principle, RF spectroscopy hadzibabic (); grimm () would directly detect the difference in the mean field shift. Although the calculation was carried out for a single impurity, we expect these results to be quantitatively correct for a small density of fermionic impurities over a bosonic Mott insulator background. Essentially, if the inverse of Fermi momentum is much larger than the lattice spacing, then the fermions would hardly effect each other’s behavior.

In this section, we obtained the phase diagram for the interaction of one impurity particle with a perfect Mott insulator. In the next three sections, we investigate how this ideal situation is affected by lattice anisotropy, higher impurity bands, and fluctuations of the Mott insulator.

Iii Effects of Lattice Anisotropy

In the optical lattice experiments, it is possible to change the strength of the laser beams forming the lattice, hence realize a model system where the lattice is not isotropic. For a quantitative comparison with experiment, it is necessary to take this effect into account. We assume that the hopping strength for the fermion is different in one direction compared to the other two directions and calculate the effect of such anisotropy on the phase diagram of the previous section.

The localization threshold for the anisotropic case can also be calculated analytically. Thus, in the following discussion we need not assume that the anisotropy of the lattice is small. As a simple limit, we obtain the two dimensional lattice localization problem as the layers are decoupled. For that limit, even the smallest attraction between the fermion and bosons causes bound state formation. As the coupling between the two dimensional layers is increased the critical interaction needed for bound state formation monotonically increases.

Because of the anisotropy, the single particle Hamiltonian in Eq. (4) is modified as


where we take the hopping term in the direction to be . We obtain the relation between and as


where and . For this integral can be evaluated exactly delves (). The critical value for the localization is found to be




As increases, i.e. the anisotropy of the lattice increases, the critical value for the potential increases and the localization becomes more difficult (Fig. 5). As , the system becomes two dimensional and there is no threshold for localization, as expected. For large , , which gives in the one dimensional limit, . Moreover, it is possible to evaluate the integral for nonzero delves (), yielding


where , and

Figure 5: The critical value of interaction as a function of lattice anisotropy characterized by (Eq. (15)). If the hopping parameter increases, anisotropy of the lattice increases and the localization becomes more difficult. gives for the isotropic case.

Using this exact result, the phase diagram can be obtained for arbitrary . In Fig. 6, we display the phase diagram for . Comparing Fig. 6 with Fig. 2 (isotropic case) we see that the phase boundaries are closer to the axis.

Figure 6: Phase diagram in the presence of lattice anisotropy (for and ). To be compared with Fig. 2 (, ). One can see that anisotropy with causes the localization threshold to move to higher values of .

Iv Effects of higher impurity bands

An important point one always has to keep in mind about the optical lattice experiments is that the effective Hubbard models, such as Eq. (4), are obtained by projecting the system into the lowest band of the lattice jaksch (). This procedure is expected to describe the low energy physics as long as the band gaps are larger than the temperature and interaction scales in the problem. In the equivalent language of Wannier functions, this condition corresponds to requiring the Wannier function of each lattice site to be undisturbed by interactions.

Figure 7: Schematic representation of the effect of higher impurity bands to the hopping parameter. If the localized impurity attracts extra particles (holes) to the localization site, the local wave function of the impurity particle changes. Then the hopping parameter for this site is different from that for the other sites.

In the context of the current problem, we discussed the critical hopping strength that is needed to localize the impurity particle to a small region, which is of the order of one lattice site. The precise determination of the Hubbard model parameters such as depends on the microscopic model one starts from. For the Hubbard model to work correctly, the Wannier functions for the impurity must be unchanged even if the impurity particle is localized to one lattice site. As a localized impurity attracts (or repels) extra particles (holes) to its localization site, one may expect the on-site wave function of the localized particle to be different from the Wannier functions at other lattice sites. This is essentially considering the coupling of the localized particle to higher impurity bands, and should be a small effect controlled by the parameter , where is the width of the first band gap of the impurity bands. Thus, the effect we are considering in this section would be important only if the impurity particle is highly mobile in the lattice, while the interaction between the background particles and the impurity is strong enough to localize the particle (Fig. 7).

In such a case, the system can still be modelled by a Hubbard model where the hopping strength between the localization site and its neighbors () is different from the hopping strength between any other neighboring sites in the lattice (). These hopping strengths can once again be calculated by looking at the overlaps of the localized wavefunctions between neighboring lattice sites jaksch ().

In this case, we take the single particle Hamiltonian as


Calculations similar to those performed in the previous sections yield


where . Evaluating this integral exactly joyce2 (), we obtain the relation


where , , and . When , the critical value for the potential is found to be


As can be seen in Fig. 8, the coupling to higher impurity bands can substantially change the critical value for localization.

Figure 8: The critical value of interaction as a function of (Eq. (22)). As the ratio of the hopping parameters increases, localization occurs for smaller values of the interaction.

If , we expect a narrowing of the local wave function (as in Fig. 7), then and consequently localization is harder . Similarly if we expect easier localization. In general one would then expect each different polaron state to have a different value. Still to gain a basic understanding of this effect we obtain the phase diagram (Fig. 9) using constant value of .

Figure 9: Phase diagram obtained when the effect of higher impurity bands is taken into account (, ). This effect is modelled by the parameter , which is the ratio of the hopping strength between the localization site and its neighbors to the one between any other neighboring sites. Compare this figure with Fig. 2 (, ). One can see that localization is easier if .

V Effects of Fluctuations of the Mott insulator

The ideal Mott insulator state is achieved only when the boson hopping term in the Bose-Hubbard Hamiltonian is zero. When there is a small but nonzero hopping probability for the bosons, the ground state of the MI contains virtual particle and hole excitations. When gets larger, these excitations gain more amplitude and will finally destroy the order in the MI, and cause a transition to the SF state.

In the previous sections, we considered the MI to be a perfect insulator by setting boson hopping to zero. In this section, we consider the effects of small but nonzero hopping on our calculations. We always work in the limit that the boson hopping is the smallest energy scale of the system and require that the system be away from the MI-SF transition boundaries. Under these conditions, we can use perturbation theory to investigate the effects of particle and hole fluctuations on the localization of the impurity particle.

There are two main consequences of turning on the boson hopping. First, the compound object formed by the impurity and the extra bosons (polaron) will become mobile. The width of the polaron band will be proportional to , where is the number of extra bosons (or holes) forming the polaron. As is the smallest energy scale in the problem, this polaron mobility effect will be small, especially for polarons containing more than one boson or hole. In other words as , the fermion is much more mobile than the boson (or the combined polaron). The mobility of the polaron will then have a small effect on the localization problem considered above.

The second effect is the change of effective mass (or effective hopping ) of the impurity particle due to scattering from the fluctuations of the MI. The first contribution to such processes will be proportional to for the fermion to scatter from a particle-hole pair, and for the pair to be excited. Because of the presence of which is large in the limit that we are interested in, this second effect will be dominant over polaron mobility. We should also mention that it is only lowest orders of the perturbation theory that these two effects can be considered separate processes.

To the lowest order in the boson wavefunction can be written as: . The state of the system without boson-fermion interaction is then , where is the wavefunction of a fermion moving with lattice momentum .

By treating the boson-fermion interaction as a perturbation


we can write the second order energy shift as


where the only non-vanishing term is


Again considering the lowest non-vanishing order contribution we get:


Hence, we see that the first effect of the fluctuations of the MI is to renormalize the hopping strength of the impurity particle to


Essentially a decreasing means the effective mass of the particle is larger, and it becomes easier to localize the particle. Thus, the effects of boson hopping on all the situations discussed in previous sections can be obtained by scaling the phase diagrams by the renormalized value of .

Vi Conclusions

We study the localization of a single impurity particle in a lattice containing a Mott insulator of bosons. This is a simple limit of the complex physics of the Bose-Hubbard model for mixtures of different atomic species. The impurity particle has two distinct types of behavior; it may either move freely throughout the lattice or may choose to localize at a certain lattice site by attracting extra bosons or holes. In the limit of a perfect Mott insulator, we calculate the boundary between these two phases as well as the number of extra bosons (holes) forming the bound state exactly. Our result for the phase diagram is given in Fig. 2.

We believe that this phase diagram can be checked experimentally. In recent experiments boson-fermion mixtures were created in the parameter regimes that we consider in this paper gunter (); ospelkaus (). If the density of fermions is small enough, one can disregard the many-particle nature of the fermions and consider the localization problem for one of them. To determine whether a fermion is localized, RF spectroscopy would be an ideal tool. The mean-field shift of a localized fermion would be larger as it sees more background bosons on the average. Equivalently, one can say that the binding energy of the fermion calculated in Section II would be directly reflected in the RF spectra for fermions. Different polaron states, containing multiple numbers of bosons or holes can be distinguished similarly.

The calculations presented in this paper were carried out with the assumption of a homogenous infinite system. In cold gas experiments there is always a confining trap, which in general complicates the correspondence between the infinite system predictions and experimental results. However, when a lattice boson system is driven deep into the MI regime, the density profile of the system consists of MI plateaus separated by thin SF regions. Within each MI region, the fermion would see a flat interaction potential, and when the fermion is localized, the width of the wavefunction of the fermion becomes of the order of few lattice sites. Thus, we expect our calculations to accurately reflect the transition boundaries, especially for transitions between different polaron states.

After the ideal case, we consider three effects which may play a role in the experiments. The first case we consider is the possibility of tunnelling anisotropy. As the strength of the lattice is determined by the laser intensity (and beam waist), experiments occasionally have such anisotropy. We generalize the exact results to this case and find that it becomes easier to localize the impurity when the system becomes more two dimensional, as expected. Next, we consider the effect of higher bands of the impurity particle. We argue that this effect can be taken into account by modifying the hopping strengths between the localization site and its neighbors, and obtain the phase diagram. Finally, we consider the effect of small but finite tunnelling strength for the bosons forming the MI, and calculate the effective hopping strength for the impurity particle in the presence of fluctuations.

We believe that our exact results about impurity localization on a Mott insulator background provide a starting point for the investigation of the complex phase diagram of mixtures in optical lattices.

S.S. is supported by TUBITAK Grant No. 106T052. R.O.U. is supported by TUBITAK. M.Ö.O. is supported by TUBA-GEBIP and TUBITAK Kariyer Grant No. 104T165. M.Ö.O. wishes to thank Hui Zhai and B. Tanatar for useful discussions and acknowledge the hospitality of Institut Henri Poincaré-Centre Émile Borel, where part of the work was completed.


  • (1) M. Greiner, O. Mandel, T. Esslinger, T.W. Hansch, and I. Bloch, Nature 415, 39 (2002)
  • (2) B. Paredes, A. Widera, V. Murg, O. Mandel, S. Folling, I. Cirac, G.V. Shlyapnikov, T.W. Hasch, and I. Bloch, Nature 429, 277 (2004)
  • (3) J.K. Chin, D.E. Miller, Y. Liu, C. Stan, W. Setiawan, C. Sanner, K. Xu, and W. Ketterle, Nature 443, 961 (2006)
  • (4) T. Stoferle, H. Moritz, K. Gunter, M. Kohl, and T. Esslinger, Phys. Rev. Lett. 96, 030401 (2006)
  • (5) M. Anderlini, P.J. Lee, B.L. Brown, J. Sebby-Strabley, W.D. Phillips, and J.V. Porto, Nature 448, 452 (2007)
  • (6) K.D. Nelson, X. Li, and D.S. Weiss, Nature Phys. 3, 556 (2007)
  • (7) E. Altman, E. Demler, and M.D. Lukin, Phys. Rev. A 70, 013603 (2004)
  • (8) I. Bloch, J. Dalibard, and W. Zwerger, arXiv:07043011
  • (9) J. Goldwin, S.B. Papp, B. DeMarco, and D.S. Jin, Phys. Rev. A 65, 021402(R) (2002)
  • (10) H. Ott, E. de Mirandes, F. Ferlaino, G. Roati, G. Modugno, and M. Inguscio, Phys. Rev. Lett. 92, 160601 (2004)
  • (11) K. Mølmer, Phys. Rev. Lett. 80, 1804 (1998)
  • (12) H. Heiselberg, C.J. Pethick, H. Smith, and L. Viverit , Phys. Rev. Lett. 85, 2418 (2000)
  • (13) A. Albus, F. Illuminati, and J. Eisert, Phys. Rev. A 68, 023606 (2003)
  • (14) A.B. Kuklov and B.V. Svistunov, Phys. Rev. Lett. 90, 100401 (2003)
  • (15) M. Lewenstein, L. Santos, M.A Baranov, and H. Fehrmann, Phys. Rev. Lett. 92, 050401 (2004)
  • (16) U. Gavish and Y. Castin, Phys. Rev. Lett. 95, 020401 (2005)
  • (17) K. Sacha and E. Timmermans, Phys. Rev. A 73, 063604 (2006)
  • (18) K. Gunter, T. Stoferle, H. Moritz, M.Kohl, and T. Esslinger, Phys. Rev. Lett. 96, 180402 (2006)
  • (19) S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, K. Sengstock, and K. Bongs, Phys. Rev. Lett. 96, 180403 (2006)
  • (20) D. Jaksch, C. Bruder, J.I. Cirac, C.W. Gardiner, and P. Zoller, Phys. Rev. Lett. 81, 3108 (1998)
  • (21) M.P.A. Fisher, P.B. Weichman, G. Grinstein, and D.S. Fisher , Phys. Rev. B 40, 546 (1989)
  • (22) L.D. Landau and L.M. Lifshitz, Quantum Mechanics vol 3 , (Pergamon Press) p.162
  • (23) R.T. Delves and G.S. Joyce, J. Phys. A: Math. Gen. 34, L59 (2001), and references therein
  • (24) G.S. Joyce, J. Phys. A: Math. Gen. 5, L65 (1972)
  • (25) M.W. Zwierlein, Z. Hadzibabic, S. Gupta, and W. Ketterle, Phys. Rev. Lett. 91, 250404 (2003)
  • (26) C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, J.H. Denschlag, and R. Grimm, Science 305, 1128 (2004)
  • (27) G.S. Joyce, J. Phys. A: Math. Gen. 35, 9811 (2002)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description