Recent experiments on ultracold atomic alkali gases in a one-dimensional optical lattice have demonstrated the transition from a gas of soft-core bosons to a Tonks-Girardeau gas in the hard-core limit, where one-dimensional bosons behave like fermions in many respects. We have studied the underlying many-body physics through numerical simulations which accommodate both the soft-core and hard-core limits in one single framework. We find that the Tonks-Girardeau gas is reached only at the strongest optical lattice potentials. Results for slightly higher densities, where the gas develops a Mott-like phase already at weaker optical lattice potentials, show that these Mott-like short-range correlations do not enhance the convergence to the hard-core limit.
This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems.
We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the
1/k4 momentum tail.
An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling.
In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.
Bosonic atoms confined in optical lattices are described by the Bose-Hubbard model and can exist in two different phases, Mott insulator or superfluid, depending on the strength of the system parameters. In the vicinity of the phase boundary, there are degeneracies that occur between every two adjacent Mott lobes. Because of this, nondegenerate perturbation theory fails to give meaningful results for the condensate density: it predicts a phase transition in a point of the phase diagram where no transition occurs. Motivated by this, we develop two different degenerate perturbative methods to solve the degeneracy-related problems. Moreover, we study the one-dimensional repulsively interacting Bose gas under harmonic confinement, with special attention to the asymptotic behavior of the momentum distribution, which is a universal
k−4 decay characterized by the Tan's contact. The latter constitutes a direct signature of the short-range correlations in such an interacting system and provides valuable insights about the role of the interparticle interactions. We investigate the system constituted of
N interacting particles in the strongly interacting limit. In such a regime, the strong interparticle interaction makes the bosons behave similarly to the ideal Fermi gas. Because of the difficulty in analytically solving the system for
N particles at finite interaction, the Tonks-Girardeau regime provides a favorable scenario to probe the contact. Therefore, we are able to provide an analytical formula for the Tan's contact. Furthermore, we analyze the scaling properties of the Tan's contact in terms of
N in the high-temperature regime as well as in the strongly interacting regime. Finally, we compare our analytical calculations of the Tan's contact to quantum Monte Carlo simulations and discuss some fundamental differences between the canonical and the grand-canonical ensembles.
We report the observation of a one-dimensional (1D) Tonks-Girardeau (TG) gas of bosons moving freely in 1D. Although TG gas bosons are strongly interacting, they behave very much like noninteracting fermions. We enter the TG regime with cold rubidium-87 atoms by trapping them with a combination of two light traps. By changing the trap intensities, and hence the atomic interaction strength, the atoms can be made to act either like a Bose-Einstein condensate or like a TG gas. We measure the total 1D energy and the length of the gas. With no free parameters and over a wide range of coupling strengths, our data fit the exact solution for the ground state of a 1D Bose gas.
We summarize in the present work exact results obtained for Tonks-Girardeau gases on one-dimensional optical lattices both for the ground state and nonequilibrium dynamics. On the theoretical side, impenetrable bosons offer the opportunity to study strongly interacting systems in one-dimensional lattices exactly, by means of the Jordan-Wigner transformation, and hence contribute to the topic of strong correlations at the center of interest in both condensed matter physics and quantum gases. This motivation is further enhanced by recent experimental realizations of such systems with ultracold atoms. After having shown their universal properties in equilibrium, we concentrate on their nonequilibrium dynamics. It will be shown that, starting from a pure Fock state, quasi-long-range correlations develop dynamically and lead to the formation of quasicondensates with a momentum determined by the underlying lattice. We expect this effect to be relevant for atom lasers with full control of the wavelength. Then, we will show that the free evolution of an initially confined Tonks-Girardeau gas leads to a momentum distribution that approaches at long times that of the equivalent fermionic system, giving rise to a bosonic gas with a Fermi edge, and hence a fermionization that can only be obtained out of equilibrium. Remarkably, although the momentum distribution function of the Tonks-Girardeau gas becomes equal to the one of the fermions, no loss in coherence is observed in the system, as reflected by a large occupation of eigenstates of the one-particle density matrix.
We consider a homogeneous 1D Bose gas with contact interactions and a large attractive coupling constant. This system can be realized in tight waveguides by exploiting a confinement induced resonance of the effective 1D scattering amplitude. By using the diffusion Monte Carlo method we show that, for small densities, the gaslike state is well described by a gas of hard rods. The critical density for cluster formation is estimated using the variational Monte Carlo method. The behavior of the correlation functions and of the frequency of the lowest breathing mode for harmonically trapped systems shows that the gas is more strongly correlated than in the Tonks-Girardeau regime.