One of the challenging features of studying model Hamiltonians with cold atoms in optical lattices is the presence of spatial inhomogeneities induced by the confining potential, which results in the coexistence of different phases.This paper presents quantum Monte Carlo results comparing methods for confining fermions in two dimensions, including conventional diagonal confinement, a recently proposed ``off-diagonal confinement'', as well as a trap which produces uniform density in the lattice. At constant entropy and for currently accessible temperatures, we show that (1) diagonal confinement results in the strongest signature of magnetic order, primarily because of its judicious use of entropy sinks at the trap edge and that (2) for d-wave pairing, a trap with uniform density is optimal and can be effectively implemented via off-diagonal confinement. This feature is important to any prospective search for superconductivity in optical lattices.
One of the challenging features of studying model Hamiltonians with cold atoms in optical lattices is the presence of spatial inhomogeneities induced by the confining potential, which results in the coexistence of different phases.This paper presents quantum Monte Carlo results comparing methods for confining fermions in two dimensions, including conventional diagonal confinement, a recently proposed ``off-diagonal confinement'', as well as a trap which produces uniform density in the lattice. At constant entropy and for currently accessible temperatures, we show that (1) diagonal confinement results in the strongest signature of magnetic order, primarily because of its judicious use of entropy sinks at the trap edge and that (2) for d-wave pairing, a trap with uniform density is optimal and can be effectively implemented via off-diagonal confinement. This feature is important to any prospective search for superconductivity in optical lattices.
In a recent letter [Phys. Rev. Lett. 104, 167201 (2010)] we proposed a new confining method for ultracold atoms on optical lattices, which is based on off-diagonal confinement (ODC). This method was shown to have distinct advantages over the conventional diagonal confinement (DC), which makes use of a trapping potential, such as the existence of pure Mott phases and highly populated condensates. In this manuscript we show that the ODC method can also lead to lower temperatures than the DC method for a wide range of control parameters. Using exact diagonalization we determine this range of parameters for the hard-core case and then we extend our results to the soft-core case by performing quantum Monte Carlo (QMC) simulations for both DC and ODC systems at fixed temperature and analyzing the corresponding entropies. We also propose a method for measuring the entropy in QMC simulations.
The low-temperature properties of the two-dimensional attractive Hubbard model are strongly influenced by the fermion density. Away from half-filling, there is a finite-temperature transition to a phase with s-wave pairing order. However, Tc is suppressed to zero at half-filling, where long-range charge-density-wave order also appears, degenerate with superconductivity. This paper presents determinant quantum Monte Carlo simulations of the attractive Hubbard model in the presence of a confining potential Vtrap which makes the fermion density $\ensuremath{\rho}$ inhomogeneous across the lattice. Pair correlations are shown to be large at low temperatures in regions of the trapped system with incommensurate filling, and to exhibit a minimum as the local density $\ensuremath{\rho}(\mathbit{i})$ passes through one fermion per site. In this ring of $\ensuremath{\rho}(\mathbit{i})=1$, charge order is enhanced. A comparison is made between treating Vtrap within the local-density approximation (LDA) and in an ab initio manner. It is argued that certain sharp features of the LDA result at integer filling do not survive the proximity of doped sites. The critical temperature of confined systems of fixed characteristic density is estimated.
A major focus in condensed matter physics is to study the origin of exotic quantum phases such as coexistent and inhomogeneous phases, quantum criticality, and secondary ordered phases close to quantum critical points. Exotic phenomena in strongly correlated systems occur due to competing complex interactions of spin, charge, lattice, and orbital degrees of freedom. Complex quantum phases in strongly correlated systems are challenging, but they might be very useful due to their possible functionality to make advance devices. In order to understand, utilize, and optimize such behaviors, we need to improve our understanding of these systems. Studies in cold atom systems are of interest since cold atom experiments provide a control on model parameters. In this thesis, we use novel analytical and computational techniques to treat strongly interacting bosonic systems. Taking advantage of the very versatile quantum Monte Carlo Stochastic Green Function algorithm, we studied several interesting problems. First, we explore a recently developed new confining method for cold atoms on optical lattices. Atoms are confined via a hopping integral that decreases as a function of the distance from the center of the lattice. This method might lead to lower temperatures than existing diagonal confinement methods. Next, we study the ground state phase diagram of interacting bosons on a ring-shape lattice with a region of weak hopping integrals. The model, an extension of the well known Bose-Hubbard model, develops a novel local Mott phase in addition to the usual Mott and superfluid phases in the homogeneous system. This might provide a new insight to the description of atomtronics applications. Finally, we study the two species Bose Hubbard model in a two-dimensional lattice. This model presents novel phases due to the complexity associated with multiple species. Its phase diagram shows ordered and coexistence phases including a ferromagnetic phase separated phase with high entropy. This phase might be accessible experimentally. The novel phases found from our studies are linked to experiments on ultra-cold atoms trapped by laser beams.
The emergence of local phases in a trapped two-component Fermi gas in an optical lattice is studied using quantum Monte Carlo simulations. We treat temperatures that are comparable to or lower than those presently achievable in experiments and large enough systems that both magnetic and paired phases can be detected by inspection of the behavior of suitable short-range correlations. We use the latter to suggest the interaction strength and temperature range at which experimental observation of incipient magnetism and d-wave pairing are more likely and evaluate the relation between entropy and temperature in two-dimensional confined fermionic systems.
The ability to cool fermions in optical lattices to ultra cold temperatures has led to an interdisciplinary area of research, that has attracted a lot of attention in recent years. An interesting development in this area is the possibility to realize models for strongly correlated fermions in the laboratory, such as the fermionic Hubbard Model. Determinantal Quantum Monte Carlo simulations have proven to be an important tool in the study of fermionic atoms. Nonetheless, it is important to compare the results and efficiency of different methods. Here comparisons with Numerical Linked Cluster Expansion and Dynamical Mean Field Theory data for double occupation and short range correlations, both relevant to current optical lattice experiments, will be presented and discussed. Another topic relevant in the context of optical lattice experiments is the study of metal insulator transitions. Indeed, the Mott insulating phase has been realized and observed in two-flavor mixtures of fermionic atoms loaded on optical lattices, being characterized both by the double occupation and the compressibility. An interesting point that has been addressed in the literature over the years is whether the same fermion-fermion interaction, responsible for the Mott insulating state, could drive an insulating system metallic. Here we show that, when fermions are loaded in optical lattices with spatially varying interactions a correlation induced Mott insulator to metal transition can take place. The spatial modulation of the interactions was recently demonstrated and opens the possibility for the experimental realization of such exotic phases.
We study properties of ultracold quantum gases trapped in 1D and under the influence of an underlying optical lattice. We analyze both fermionic and bosonic systems following different numerically exact approaches, like quantum Monte Carlo for soft-core bosons and two-species fermions, and a recently developed exact approach for HCB. The exposition is organized as follows. In Chap. 2, we study single species degenerated fermions on 1D lattices. We show that a splitting of the system takes place in a region of the spectrum with eigenstates that have a non-vanishing weight only in a fraction of the system. We also show that if on top of the lattice an alternating potential is introduced, doubling the original periodicity, an additional slicing' of the system can be achieved. The width and number of such regions can be controlled in a given energy range by the amplitude of the new modulation. By filling these systems the ground state may exhibit insulating regions due to the presence of local gaps, which can be observed in the local density of states. As a complement to these 1D results we present in Appendix A results obtained in 2D lattices. In Chap. 3 we review in some detail the two quantum Monte Carlo (QMC) techniques that are relevant to our study of confined soft-core bosons and two-species fermions on optical lattices, the worldline QMC and the zero temperature projector QMC, respectively. The basis of the classical Monte Carlo method are also presented, included a proof of the Central Limit theorem (Appendix B) that establishes its foundations. Properties of soft-core bosons confined on optical lattices are presented in Chap. 4. QMC simulations show that due to the presence of the slowly varying confining potential the concept of commensurability, typical for periodic systems, loses its meaning. Mott domains appears for a continuous range of fillings and always coexist with superfluid phases, as it was mentioned before. The latter feature is reflected by the behavior of the global compressibility, which never vanishes. A local compressibility, proportional to the variance of the local density, is defined in order to characterize the local Mott-insulating phases. Finally, a phase diagram for trapped systems is also presented. In Chap. 5 we study the ground state properties of two-components fermions confined on optical lattices. Like for the bosonic case, Mott domains also appear for a continuous range of fillings and always coexist with compressible phases. We define a local order parameter (that we denominate local compressibility), which characterizes the Mott insulating phases in an unambiguous way. By means of this local compressibility, we study in detail the interphase between the metallic and insulating region finding that critical behavior sets in revealing a new critical exponent. Furthermore, the behavior of the local compressibility and the variance of the density are found to be universal in this case, independently of the confining potential and the strength of the interaction, as usual in critical phenomena. We present in Chap. 6 a recently developed exact numerical approach that allows to study ground state properties of HCB confined on 1-D lattices. This exact treatment is applied to study the off-diagonal behavior of the one-particle density matrix (OPDM) and related quantities in the equilibrium case. We find that the OPDM decays as a power-law x^-1/2 for large-x, irrespective of the confining potential chosen. The power-law above is shown to determine the scaling of the occupation of the lowest natural orbital in the thermodynamic limit. This scaling and its finite size corrections are also studied for arbitrary powers of the confining potential. In addition, we find a power-law decay of the NO occupations at low densities, and show that its exponent is also universal. The low density limit in the lattice, equivalent to continuous systems, is also analyzed in detail. The approach above is generalized in Chap. 7 to study the non-equilibrium dynamics of HCB in 1-D configurations with an underlying lattice. The presence of the lattice enables the creation of pure Fock states of HCB (a HCB per lattice site) where there is no coherence in the system. We show that quasi-long range correlations develop in the equal-time-one-particle density matrix when such states are allowed to evolve freely, and that they lead to the formation of quasi-condensates of HCB at finite momentum. In addition,we obtain an universal power-law describing the population of the quasi-condensate as a function of time. Finally, we discuss how such systems can be used to create atom lasers with a wave-length that can be controlled through the lattice parameter.