PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ncommsLink to Publisher's site
 
Nat Commun. 2016; 7: 11279.
Published online 2016 April 14. doi:  10.1038/ncomms11279
PMCID: PMC4834636

Doublon dynamics and polar molecule production in an optical lattice

Abstract

Polar molecules in an optical lattice provide a versatile platform to study quantum many-body dynamics. Here we use such a system to prepare a density distribution where lattice sites are either empty or occupied by a doublon composed of an interacting Bose-Fermi pair. By letting this out-of-equilibrium system evolve from a well-defined, but disordered, initial condition, we observe clear effects on pairing that arise from inter-species interactions, a higher partial-wave Feshbach resonance and excited Bloch-band population. These observations facilitate a detailed understanding of molecule formation in the lattice. Moreover, the interplay of tunnelling and interaction of fermions and bosons provides a controllable platform to study Bose-Fermi Hubbard dynamics. Additionally, we can probe the distribution of the atomic gases in the lattice by measuring the inelastic loss of doublons. These techniques realize tools that are generically applicable to studying the complex dynamics of atomic mixtures in optical lattices.

Polar molecules with long-ranged dipolar interactions are ideally suited to the exploration of strongly correlated quantum matter and intriguing phenomena such as quantum magnetism, exotic superfluidity and topological phases1,2,3,4,5,6,7,8. The recent observation of the dipole-mediated spin-exchange interaction in an optical lattice9 and the demonstration of the many-body nature of the spin-exchange dynamics10 mark important steps for the use of polar molecules to study strongly correlated matter. While this initial work was done with a molecular filling fraction of only ~5% in a three-dimensional (3D) lattice9, more recent work has demonstrated a quantum synthesis scheme for molecule production in the lattice that relies on careful preparation of the initial atomic gases11. This led to a reduction in the final entropy of polar molecules by a factor of ~4, and, correspondingly, a much higher filling fraction of ~25% that opens up the possibility for studying non-equilibrium, many-body spin dynamics in a quantum gas of polar molecules where every molecule is connected to others.

The quantum synthesis approach reported in ref. 11 starts by preparing atomic insulator states that depend on atomic interactions, quantum statistics and low temperature12. However, realizing the full potential of this approach requires not only control over the atomic distributions, but also a detailed understanding of the molecule creation process.

Here we investigate this important step by leveraging our capability of molecule production in an optical lattice to create a clean system of doublons13. This technique allows us to additionally study 3D Bose-Fermi Hubbard dynamics. After creating ground-state molecules, we efficiently remove all unpaired atoms from the lattice and convert the molecules back to free atoms (in their lowest hyperfine states of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m1.jpg for 40K and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m2.jpg for 87Rb, where An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m3.jpg denotes the hyperfine state and its projection onto the magnetic field). This realizes a lattice where the sites are either empty or occupied by individual doublons that comprise a pair of bosonic and fermionic atoms. This well-defined initial state allows us to directly address limitations in the molecule creation process by probing the efficiency with which these doublons are converted back to molecules under various experimental conditions that affect atomic tunnelling rates, higher Bloch-band populations and the adiabaticity of a magnetic-field sweep through a higher partial-wave Feshbach resonance. Furthermore, this well-initialized, non-equilibrium state of a disordered doublon distribution provides an ideal platform to explore the many-body dynamics of a lattice-confined Bose-Fermi mixture in a regime that is beyond the current simulation capabilities.

Results

Preparing the doublons

The experiment proceeds in steps as depicted schematically in Fig. 1. To prepare the doublons, we create a sample of molecules in their ro-vibrational ground state in the lattice as described in ref. 11 and then remove unpaired atoms with resonant light, so that all lattice sites are either empty or contain a single molecule. We then transfer the ground-state molecules back to a weakly bound Feshbach molecule state, followed by a magnetic-field (B) sweep to above the resonance to create a clean system of doublons. The solid black line in the upper panel of Fig. 1 shows schematically B relative to the s-wave Feshbach resonance (dashed line) that is used to manipulate the atomic inter-species interactions and to create molecules. After this preparation, the doublons are left to evolve in the lattice for a variable time An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m4.jpg. Our measurement then consists of sweeping B to below the resonance to associate atoms into Feshbach molecules and determining the fraction of K atoms that form molecules. Specifically, we measure the molecule number using the following protocol. We first apply radio frequency (rf) to spin-flip the unpaired K atoms to another hyperfine state, which renders the unpaired K atoms invisible for subsequent molecular detection. We then sweep B back above the resonance to dissociate the molecules, and measure the number of resulting K atoms by spin-selective resonant absorption imaging. The conversion efficiency is determined by dividing this molecule number by the total number of K atoms measured when we do not apply the rf.

Figure 1
A schematic of the experiment.

d-wave Feshbach resonance

We begin by investigating a narrow d-wave Feshbach resonance14,15,16,17 that is located less than 0.1 mT above the 0.3-mT-wide s-wave resonance that is used for making molecules (Fig. 2a). With a pair of atoms confined in the same lattice site, the on-site density is orders of magnitude higher than that in ordinary optical traps, and thus this narrow Feshbach resonance can adversely affect the magneto-association process, where B is swept down from above the s-wave resonance to create molecules. Crossing the d-wave resonance too slowly will produce d-wave molecules, which will not be coupled to the ground state by the subsequent STIRAP laser pulses, as the process is weak and off resonance. If B is swept sufficiently fast to be diabatic for this narrow resonance (but still slow enough to be adiabatic for the broad s-wave resonance), crossing the d-wave resonance has no impact; however, the high effective densities at each site in an optical lattice can make it challenging to sweep fast enough. Although we study here specific resonances for the K-Rb system, the possibility of having to cross other Feshbach resonances and the issue of sweep speeds are general to magneto-association of atoms in an optical lattice.

Figure 2
The d-wave Feshbach resonance.

In the experiment illustrated in Fig. 1, we investigate the d-wave resonance by varying the rate, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m5.jpg, and the final value, Bhold, of the sweep that creates doublons. We then measure the subsequent molecule conversion efficiency after An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m6.jpg using a fast 1.68-mT ms−1 magneto-association sweep. Figure 2b illustrates relevant states, above and below the resonance, for two atoms in a lattice site18,19: at low fields these are the s-wave molecule (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m83.jpg), d-wave molecule (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m87.jpg) and unbound atoms (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m86.jpg), and at high fields these are unbound atoms in the ground band of the lattice (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m84.jpg) and atoms with a band excitation in their relative motion (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m85.jpg). For simplicity, we illustrate states for a harmonic potential whose trap frequency An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m7.jpg is the same for both atoms, with eigenstates of relative motion denoted by An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m8.jpg. The dashed arrows show the diabatic (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m83.jpgAn external file that holds a picture, illustration, etc.
Object name is ncomms11279-m84.jpg) and adiabatic (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m83.jpgAn external file that holds a picture, illustration, etc.
Object name is ncomms11279-m85.jpg) trajectories for the dissociation of s-wave Feshbach molecules when crossing the d-wave resonance, while the solid arrows show the diabatic trajectories (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m84.jpgAn external file that holds a picture, illustration, etc.
Object name is ncomms11279-m83.jpg and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m85.jpgAn external file that holds a picture, illustration, etc.
Object name is ncomms11279-m86.jpg) for the subsequent, fast magneto-association sweep.

Figure 2c shows the measured molecule conversion efficiency as a function of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m9.jpg when sweeping across the d-wave resonance from 54.56 to 56.24 mT (there are no other resonances in this field range), while Fig. 2d shows the effect of the final field B for a relatively slow, 0.018-mT ms−1, sweep. The data are taken for lattice depths of Vlatt=35ER (circles) and 30ER (diamonds), where ER=ħ2k2/(2mRb) is the recoil energy for Rb, mRb is the Rb atom mass, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m10.jpg and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m11.jpg. For our highest sweep rates, or when Bhold is below the d-wave resonance, the measured molecule conversion efficiency is near unity. This high conversion of doublons20,21,22 is crucial for the quantum synthesis approach to producing molecules with a high filling fraction in the lattice. The near-unity conversion also provides an excellent starting point for diagnosing potential limitations to molecule production, and the data in Fig. 2c,d clearly show the negative effect that the d-wave resonance can have on magneto-association in the lattice.

The lines in Fig. 2c,d show fits used to extract the width (Δd) and position of the d-wave resonance. We use a Landau-Zener formalism23 where the probability to cross the d-wave resonance diabatically, and therefore create s-wave Feshbach molecules in the subsequent magneto-association step, is An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m12.jpg, where A depends on the on-site densities and the Feshbach resonance parameters. By approximating the sites in the deep optical lattice as harmonic oscillator potentials, we extract Δd using An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m13.jpg, where An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m14.jpg is the harmonic trap frequency for relative motion of the two atoms (see Methods) and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m15.jpg is the harmonic oscillator length with the doublon reduced mass An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m16.jpg (ref. 24) (we note that the right-hand side of equation (26) in this reference is missing a factor of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m17.jpg). From an exponential fit (line in Fig. 2c), A=0.110(7) mT ms−1, and using a background scattering length of abg=−187(5)a0 (ref. 25), where a0 is the Bohr radius, we extract a width of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m18.jpg. By fitting an error function (line) to the data in Fig. 2d, we determine the location of the resonance to be 54.747(1) mT, which is consistent with previous experiments where atom loss was observed14,16.

The precise determination of the width of the d-wave resonance allows us to gauge its significance in molecule creation. Our typical sweep rate of 0.34 mT ms−1 for magneto-association, which has remained the same since the first creation of KRb molecules in an optical lattice22, gives ~70% probability of being diabatic when crossing the d-wave resonance. This suggests that we create a substantial fraction of d-wave molecules that are dark to our detection (in Fig. 2b). These d-wave molecules may have played a role in limiting the lattice filling fraction for polar molecules achieved in ref. 11.

Short-time tunnelling dynamics

Tunnelling dynamics of doublons in the lattice26 can also affect molecule production. In the quantum synthesis approach, achieving a high lattice filling for molecules requires not only the preparation of a large fraction of lattice sites that have doublons, but also that these doublons are not lost due to tunnelling and/or collisions prior to conversion to molecules. In our system, K feels a lattice depth that, in units of recoil energy, is 2.6 times weaker than for Rb due to differences in atomic mass and polarizability. Consequently, K tunnels faster than Rb. While a sufficiently deep lattice can prevent tunnelling of both K and Rb, practically this may not be possible in all cases, especially for polar molecule production using two atomic species that have large differences in mass and polarizability.

Figure 3 illustrates doublon dynamics due to the interplay between tunnelling and interactions, which we control by varying the lattice depth, interspecies scattering length aK-Rb and band population. The fraction of doublons that remain after An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m19.jpg is essentially equal to the measured molecule conversion efficiency described above. We note that for aK-Rb>−850a0, the B sweep crosses the d-wave Feshbach resonance with a An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m20.jpg that varies from 0.5 to 1.9 mT ms−1. Using our measured width of the d-wave resonance, the data presented in Fig. 3 have been multiplied by a factor that increases the doublon fraction to account for the finite An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m21.jpg when crossing the d-wave resonance. Figure 3a shows the effect of the lattice depth for An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m22.jpg at three different values of Bhold, corresponding to different values of aK-Rb. This timescale is relevant for both molecule production and K tunnelling dynamics. We observe that the remaining doublon fraction is highly sensitive to the lattice depth for weak interspecies interactions, for example, aK-Rb=−220a0, with a lower doublon fraction for shallower lattices that exhibit higher tunnelling rates. For stronger interactions, the dependence on lattice depth becomes less significant and almost disappears in the strongly interacting regime, for example, aK-Rb=−1,900a0. Similar behaviour is observed if we fix the lattice depth but vary the interspecies interactions, as shown in Fig. 3b.

Figure 3
Interaction and tunnelling dynamics of doublons in the lattice.

Modelling

The data in Fig. 3 clearly show evidence of decay of doublons due to tunnelling that is affected by both the lattice depth and interspecies interactions. We can model these doublon dynamics with the following Hamiltonian:

An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m23.jpg

where η=0 and 1 denote, respectively, the ground and the first excited lattice bands. The first and second terms are the kinetic energy of the K and Rb atoms, respectively. Here, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m24.jpg is the bosonic annihilation (creation) operator for a Rb atom at lattice site i in the lowest band, and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m25.jpg is the fermionic annihilation (creation) operator for a K atom at lattice site i and band An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m26.jpg. We use An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m27.jpg to indicate nearest-neighbour hopping between sites i and j with matrix element An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m28.jpg with An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m29.jpg or Rb. The third term describes the inter-species on-site interactions with matrix element An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m30.jpg. The last term is the on-site intra-species interaction between ground-band Rb atoms with strength An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m31.jpg, with An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m32.jpg as the occupation of site i.

The tunnelling rates and interaction energies are calculated for a particular Vlatt and aK-Rb (ref. 27). For example, for Vlatt=10ER, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m33.jpg, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m34.jpg. The solid curves in Fig. 3a,b show the calculations based on the Hamiltonian given in equation (1), where we have neglected Rb tunnelling by setting An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m35.jpg. We start with a single doublon, evolve the K for a hold time An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m36.jpg, and then extract the doublon fraction from the probability that the K atom remains on the same site as the Rb atom. In this treatment, we ignore the role of the magnetic-field sweeps. Calculations for a single doublon (solid lines), where the initial decay scales as An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m37.jpg, agree well with the data, except at doublon fractions below ~30%, where the disagreement arises from the finite probability in the experiment that a K atom finds a different Rb partner. Simulating a Gaussian distribution of doublons with 10% peak filling accounts for this effect (dashed lines) (see Methods). The good agreement of these calculations with the data shows that tunnelling of K, which is suppressed for deeper lattices, is the dominant mechanism for the reduction of the doublon fraction at short (~1 ms) times. The on-site interaction with Rb suppresses the K tunnelling when the interaction energy becomes larger than the width of the K Bloch band28.

When studying doublon dynamics measured for two different initial atom conditions, we find indirect evidence for excited-band molecules. Here, we compare results for our usual molecule preparation using atomic insulators to a case where we start with a hotter initial atom gas mixture at a temperature above that for the Rb Bose-Einstein condensation transition. Using a band-mapping technique, we measure the initial population of K in the ground and first excited band, as shown in Fig. 3di and ii (see Methods). We find that 11(2)% of the K atoms occupy the first excited band for the colder initial atom gas (these conditions are similar to those in ref. 11 and are used in all the measurements described in this work, except for the green squares in Fig. 3c). When starting with the hotter atom gas, we measure a significantly higher K excited-band population of 31(6)%. When looking at doublon dynamics for these two cases (Fig. 3c), we observe a lower doublon fraction for the hotter initial gas for Vlatt≤25ER. These data are taken for 1.68 mT ms−1 sweeps, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m38.jpg and aK-Rb=−220a0.

The lower doublon fraction can be explained by excited-band K atoms, which have a high tunnelling rate (An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m39.jpg and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m40.jpg are 89.3 and 1110 Hz, respectively, for Vlatt=25ER). The presence of excited-band K atoms suggests that the B sweeps for magneto-association (and dissociation) couple excited-band K atoms (plus a ground-band Rb atom) to excited-band Feshbach molecules. Moreover, the data suggest that the conversion efficiency for the excited-band Feshbach molecules is still high for Vlatt≤25ER since the observed difference in the initial excited-band K atoms is similar to the observed difference (roughly 20%) in the doublon fraction (Fig. 3c). Since, in our preparation scheme, the doublons are directly formed from the dissociation of ro-vibrational ground-state molecules, these results further indicate that a polar molecule sample prepared from a finite-temperature atom gas can contain a small fraction of molecules in an excited motional state in the lattice. We also observed a Rb excited-band population of 31(5)% after loading the thermal gas in the lattice; however, even for the excited band, the off-resonant Rb tunnelling is slow compared to the 1-ms time scale of the measurements presented in Fig. 3.

The green dashed curve in Fig. 3c shows the theoretical results for a K excited-band fraction of 24%. For comparison, the red solid curve, which is the same as the red curve in Fig. 3a, includes no excited-band population. The estimated excited band fraction ignores the effects of harmonic confinement on tunnelling, which are more significant for the hotter initial atom gas, where the resulting molecular cloud is also larger. For the hotter initial atom gas, the green dashed curve overlaps the data at the shallower lattice depths, but deviates from the measured doublon fraction at larger lattice depths (the excited band fraction of the initial K gas is independent of lattice depth). This may be expected since in the limit of a very deep lattice and a fully adiabatic magneto-association sweep, one expects that only the heavier atom (Rb) in excited bands (plus a ground-band K atom) will couple to excited-band Feshbach molecules. Future studies of the magneto-association process in a lattice for systems such as K-Rb where centre-of-mass and relative motion are coupled29 would be interesting and relevant to polar molecule preparation.

Long-time tunnelling dynamics

In Fig. 4, we present data taken for An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m41.jpg up to 40 ms, in order to look for the effects of Rb tunnelling. Measurements of the remaining doublon fraction are shown for two lattice depths (10ER and 15ER) and two values of aK-Rb (−910a0 and −1900a0). In Fig. 4, the doublon fraction has been normalized by the measured value for An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m42.jpg in order to remove the effect of the shorter-time dynamics that are presented in Fig. 3a,b. Similar to the shorter-time dynamics, at the longer hold times we observe a reduction in the doublon fraction that is suppressed for a deeper lattice and for strong inter-species interactions. Modeling these dynamics is theoretically challenging, and the lines in Fig. 4 are exponential fits that are intended only as guides to the eye. Compared to doublons composed of identical bosons13 or fermions in two-spin states30, the heteronuclear system has the additional complexities of two particle masses, two tunnelling rates and two relevant interaction energies. For example, for large aK-Rb, the interspecies interactions will strongly suppress Rb tunnelling from a doublon to a neighbouring empty site. Similarly, tunnelling of a doublon to an empty lattice site is a slow second-order process at the rate An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m43.jpg due to the energy gap of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m44.jpg (Fig. 4 inset i). However, Rb tunnelling between two neighbouring doublons, which creates a triplon (Rb-Rb-K) on one site and a lone K atom on the other site (Fig. 4 inset ii), may occur on a faster time scale due to a much smaller energy gap of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m45.jpg, which is smaller than the K tunnelling bandwidth. While the theoretical description is complicated, we observe that the time scale of the doublon decay roughly matches An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m46.jpg.

Figure 4
Longer time dynamics.

Measuring atomic distributions with doublon detection

The studies discussed thus far demonstrate that the Feshbach molecule conversion that we use to detect doublons could potentially underestimate the doublon fraction. For example, the conversion efficiency of doublons containing excited-band atoms is complicated to calculate and is likely to be less than 1. In addition, the efficiency of converting doublons to Feshbach molecules depends on the magnetic-field sweep rate, and, as shown in Fig. 2c, a very slow sweep does not always yield a unity conversion efficiency. Finally, Feshbach molecules can suffer losses from inelastic collisions with other Feshbach molecules or unpaired K atoms31, which could reduce the measured number. Given these factors and the importance of measuring the doublon fraction as a powerful diagnostic for optimizing molecule production from ultracold atoms in a lattice, we have implemented a second, complementary approach for measuring the doublon fraction using inelastic collisional loss in the initial atomic mixture, without the molecular purification step. In our system, inelastic collisions are initiated by transferring the Rb atoms from the An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m47.jpg to the An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m48.jpg state. Collisions of the An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m49.jpg Rb atoms with K can result in spin relaxation back to the Rb F=1 manifold. At B=55 mT, the An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m50.jpg state is higher in energy by h × 8.1 GHz; this inelastic collision releases a large amount of energy compared to the trap depth and therefore results in atom loss from the trap. At a collision energy corresponding to 1 μK, the calculated inelastic collision rate using the coupled channels model of ref. 15 is An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m51.jpg, and using the on-site densities in a Vlatt=25ER lattice, the resulting doublon lifetime is ~2 ms.

In Fig. 5a,b, we show example data for the number of Rb atoms as a function of time after a 2.1-ms rf sweep that transfers Rb atoms to the An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m52.jpg state. We observe a fast loss on the time scale of a few ms, followed by slower loss. We attribute the fast loss to inelastic collisions of Rb atoms in lattice sites shared with K, and the slow loss to tunnelling of atoms followed by inelastic collisions. The dashed lines in Fig. 5a,b show a fit to the sum of two exponential decays with different time constants. We can extract the fraction of Rb that is lost on the short timescale from the fits. We compared this technique with Feshbach molecule formation, and found that the two measurements generally agree.

Figure 5
Measuring the initial atomic distributions with spin-changing collisions.

As a further demonstration of the inelastic collision technique, we use this to probe the initial atomic distribution in the lattice before molecule formation, providing quantitative information on the Rb Mott insulator. Figure 5c shows the fraction of Rb atoms that are lost quickly from a Vlatt=25ER lattice after the Rb atoms are transferred to the An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m53.jpg state. For these data, we vary the initial number of Rb atoms that form a Mott insulator in the optical lattice prior to molecule creation. In Fig. 5c, the blue diamonds correspond to the data shown in Fig. 5a,b. For the data shown in circles, the fraction lost is determined by comparing the Rb number measured before to that measured 8 ms after the rf transfer. The solid curve shows a calculation of the expected loss for a Mott insulator with a temperature An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m54.jpg and a total radial harmonic confinement of 33 Hz for Rb. At low Rb number, where we expect only one Rb atom per site in the Mott insulator, the fraction lost is just the K filling fraction, assuming no double occupancy for K. For higher Rb number, double (and eventually triple and higher) occupancy in Mott shells causes a reduction in the fractional loss under the assumption of one Rb and one K lost per inelastic collision. The shaded area indicates a 10% uncertainty in the harmonic trapping frequency and 30% uncertainty in T. From the fit, we extract a K filling fraction of 0.77(2), which is in excellent agreement with the measured peak K filling reported in ref. 11. We note the previously measured fraction of Rb converted to Feshbach molecules at low Rb atom number was significantly less, at about 0.5(1) (ref. 11). This disagreement may now be attributed to the large number of K atoms present in the lattice after molecule formation, which can induce losses through inelastic collisions, and to the effect of the d-wave resonance when making molecules, as discussed above.

Discussion

Our investigation of heteronuclear doublons and their conversion to molecules by magneto-association reveals the important roles played by the lattice depth for both atomic species, the inter-species interactions, the population in excited motional states of the lattice and the magnetic-field sweep rate. The doublon dynamics uncovered in this study provides insights into the universal mechanism of their decay and atomic mixture dynamics in a 3D optical lattice, and allows preparation of optimal conditions for producing polar molecules. The highly non-equilibrium state of doublons that we use for these studies also provides an intriguing system for exploring the Hubbard dynamics of a Bose-Fermi mixture, where the behaviour of the many-body system can depend on two different tunnelling rates and two different interaction strengths28,32. This system sets the stage for performing benchmarking experiments in 1D for theory, and investigating the thermalization of an isolated many-body quantum system, including novel phases such as quasi-crystallization and many-body localization in higher dimensions33,34,35,36.

Methods

Optical trapping potentials

The preparation of the atomic gas in a 3D lattice, with a wavelength of 1,064 nm, as well as the creation of ground-state polar molecules, follows the procedures described in ref. 11. The lattice is superimposed on a crossed-beam optical dipole trap that is cylindrically symmetric. The dipole trap alone has an axial trap frequency of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m55.jpg in the vertical direction and a radial trap frequency of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m56.jpg for Rb. The measured optical trap frequencies for K are An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m57.jpg and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m58.jpg.

Width of the d-wave resonance

The scattering lengths reported in Fig. 3 have been calculated using An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m59.jpg with abg=−187(5)a0, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m60.jpg and Δs=0.304 mT (ref. 25). Including the d-wave resonance, the scattering length can be parameterized by An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m61.jpg (ref. 37). Using the relation Δd<<Δs, we can write An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m62.jpg near the d-wave resonance, where An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m63.jpg and An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m64.jpg. This has the form of an isolated resonance and we can apply the findings of ref. 24, namely that An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m65.jpg, to determine An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m66.jpg. We note that ref. 17 predicts An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m67.jpg, which is larger than our determination of An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m68.jpg.

In this determination, we ignore the coupling between the centre of mass and relative motion that arises from the fact that K and Rb experience different trapping potentials in the optical lattice. We use an effective trap frequency An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m69.jpg that governs the dynamics in the relative coordinate. Here, mRb and mK are the masses of the Rb and K atom, respectively. The trap frequency for Rb is given by An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m70.jpg, and for the 1,064 nm lattice the trap frequency for K is An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m71.jpg. For Vlatt=35ER, An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m72.jpg.

Density distribution of doublons

The dashed lines in Fig. 3a,b have been obtained by random sampling of initial doublon positions according to a Gaussian probability distribution of the filling fraction. A peak filling of 10% and widths σx=σy=6.5σz=21 sites have been used, corresponding to N≈2,000 sites occupied with a doublon. The experimentally determined cloud sizes are slightly larger (σx=25-42 sites), but we confirmed that the resulting doublon fraction is converged with respect to the cloud size. In-situ absorption images of the cloud are consistent with a Gaussian distribution of 5–10% peak filling. Tunnelling of Rb is neglected in the model, where initially each K atom is localized on a site containing a Rb atom and the doublon fraction is defined as the probability to find the K atom on a site with Rb after the evolution time An external file that holds a picture, illustration, etc.
Object name is ncomms11279-m73.jpg.

Band mapping

To measure the excited-band fraction of the initial K atoms, we use a band-mapping technique (Fig. 3d). Starting with the K atoms in the 3D lattice plus optical dipole trap potential, we turn off the lattice in 1 ms and allow the the K gas to expand in the optical dipole trap for a quarter trap period38. We image the cloud with a probe beam that propagates along the vertical direction.

Additional information

How to cite this article: Covey, J. P. et al. Doublon dynamics and polar molecule production in an optical lattice. Nat. Commun. 7:11279 doi: 10.1038/ncomms11279 (2016).

Acknowledgments

We thank M. Wall for useful discussions. We acknowledge funding for this work from NIST, NSF grant number 1125844, NSF-PIF-1211914, AFOSR-MURI, ARO-MURI, and ARO. J.P.C. acknowledges funding from the NDSEG fellowship. Part of the computation for this work was performed at the University of Oklahoma Supercomputing Center for Education and Research (OSCER).

Footnotes

Author contributions The experimental work and data analysis were carried out by J.P.C., S.A.M., M.T.M., Z.F., D.S.J. and J.Y. Theoretical modelling and calculations are provided by M.G., A.S.-N., J.S., P.S.J. and A.M.R. All authors discussed the results and contributed to the preparation of the manuscript.

References

  • Baranov M. Theoretical progress in many-body physics with ultracold dipolar gases. Phys. Rep. 464, 71–111 (2008).
  • Pupillo G., Micheli A., Büchler H. P. & Zoller P. in Cold Molecules: Theory, Experiment, and Applications (eds Krems R., Friedrich B., Stwalley W. C. CRC Press (2009).
  • Carr L. D., DeMille D., Krems R. V. & Ye J. Cold and ultracold molecules: science, technology and applications. New J. Phys. 11, 055049 (2009).
  • Lahaye T., Menotti C., Santos L., Lewenstein M. & Pfau T. The physics of dipolar bosonic quantum gases. Rep. Prog. Phys. 72, 126401 (2009).
  • Lemeshko M., Krems R. V., Doyle J. M. & Kais S. Manipulation of molecules with electromagnetic fields. Mol. Phys. 111, 1648–1682 (2013).
  • Yao N. Y. et al. . Topological flat bands from dipolar spin systems. Phys. Rev. Lett. 109, 266804 (2012). [PubMed]
  • Gorshkov A. V. et al. . Tunable superfluidity and quantum magnetism with ultracold polar molecules. Phys. Rev. Lett. 107, 115301 (2011). [PubMed]
  • Syzranov S. V., Wall M. L., Gurarie V. & Rey A. M. Spin-orbital dynamics in a system of polar molecules. Nat. Commun. 5, 5391 (2014). [PubMed]
  • Yan B. et al. . Observation of dipolar spin-exchange interactions with lattice-confined polar molecules. Nature 501, 521–525 (2013). [PubMed]
  • Hazzard K. R. A. et al. . Many-body dynamics of dipolar molecules in an optical lattice. Phys. Rev. Lett. 113, 195302 (2014). [PubMed]
  • Moses S. A. et al. . Creation of a low-entropy quantum gas of polar molecules in an optical lattice. Science 350, 659–662 (2015). [PubMed]
  • Safavi-Naini A., Wall M. L. & Rey A. M. Role of interspecies interactions in the preparation of a low-entropy gas of polar molecules in a lattice. Phys. Rev. A 92, 063416 (2015).
  • Winkler K. et al. . Repulsively bound atom pairs in an optical lattice. Nature 441, 853–856 (2006). [PubMed]
  • Zaccanti M. et al. . Control of the interaction in a Fermi-Bose mixture. Phys. Rev. A 74, 041605 (2006).
  • Julienne P. S. Ultracold molecules from ultracold atoms: a case study with the KRb molecule. Faraday Discuss. 142, 361–388 (2009). [PubMed]
  • Bloom R. S., Hu M.-G., Cumby T. D. & Jin D. S. Tests of universal three-body physics in an ultracold Bose-Fermi mixture. Phys. Rev. Lett. 111, 105301 (2013). [PubMed]
  • Ruzic B. P., Greene C. H. & Bohn J. L. Quantum defect theory for high-partial-wave cold collisions. Phys. Rev. A 87, 032706 (2013).
  • Busch T., Englert B.-G., Rzewski K. & Wilkens M. Two cold atoms in a harmonic trap. Found. Phys. 28, 549–559 (1998).
  • Köhl M., Moritz H., Stöferle T., Günter K. & Esslinger T. Fermionic atoms in a three dimensional optical lattice: observing Fermi surfaces, dynamics, and interactions. Phys. Rev. Lett. 94, 080403 (2005). [PubMed]
  • Thalhammer G. et al. . Long-lived feshbach molecules in a three-dimensional optical lattice. Phys. Rev. Lett. 96, 050402 (2006). [PubMed]
  • Greif D., Tarruell L., Uehlinger T., Jördens R. & Esslinger T. Probing nearest-neighbor correlations of ultracold fermions in an optical lattice. Phys. Rev. Lett. 106, 145302 (2011). [PubMed]
  • Chotia A. et al. . Long-lived dipolar molecules and Feshbach molecules in a 3d optical lattice. Phys. Rev. Lett. 108, 080405 (2012). [PubMed]
  • Mies F. H., Tiesinga E. & Julienne P. S. Manipulation of Feshbach resonances in ultracold atomic collisions using time-dependent magnetic fields. Phys. Rev. A 61, 022721 (2000).
  • Julienne P. S., Tiesinga E. & Köhler T. Making cold molecules by time-dependent feshbach resonances. J. Mod. Opt. 51, 1787–1806 (2004).
  • Klempt C. et al. . Radio-frequency association of heteronuclear Feshbach molecules. Phys. Rev. A 78, 061602 (2008).
  • Jürgensen O., Meinert F., Mark M. J., Nägerl H.-C. & Lühmann D.-S. Observation of density-induced tunneling. Phys. Rev. Lett. 113, 193003 (2014). [PubMed]
  • Greiner M., Mandel O., Esslinger T., Hänsch T. W. & Bloch I. Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms. Nature 415, 39–44 (2002). [PubMed]
  • Heinze J. et al. . Multiband spectroscopy of ultracold fermions: observation of reduced tunneling in attractive Bose-Fermi mixtures. Phys. Rev. Lett. 107, 135303 (2011). [PubMed]
  • Jachymski K., Idziaszek Z. & Calarco T. Feshbach resonances in a nonseparable trap. Phys. Rev. A 87, 042701 (2013).
  • Strohmaier N. et al. . Observation of elastic doublon decay in the Fermi-Hubbard model. Phys. Rev. Lett. 104, 080401 (2010). [PubMed]
  • Ospelkaus S. et al. . Quantum-state controlled chemical reactions of ultracold potassium-rubidium molecules. Science 327, 853–857 (2010). [PubMed]
  • Best T. et al. . Role of interactions in 87Rb-40K Bose-Fermi mixtures in a 3d optical lattice. Phys. Rev. Lett. 102, 030408 (2009). [PubMed]
  • Gopalakrishnan S., Martin I. & Demler E. A. Quantum quasicrystals of spin-orbit-coupled dipolar bosons. Phys. Rev. Lett. 111, 185304 (2013). [PubMed]
  • Martin I., Gopalakrishnan S. & Demler E. A. Weak crystallization theory of metallic alloys. Preprint at http://arxiv.org/abs/1506.03077 (2015).
  • Schreiber M. et al. . Observation of many-body localization of interacting fermions in a quasirandom optical lattice. Science 349, 842–845 (2015). [PubMed]
  • Smith J. et al. . Many-body localization in a quantum simulator with programmable random disorder. Preprint at http://arxiv.org/abs/1508.07026 (2015).
  • Jachymski K. & Julienne P. S. Analytical model of overlapping Feshbach resonances. Phys. Rev. A 88, 052701 (2013).
  • Murthy P. A. et al. . Matter-wave Fourier optics with a strongly interacting two-dimensional Fermi gas. Phys. Rev. A 90, 043611 (2014).

Articles from Nature Communications are provided here courtesy of Nature Publishing Group