One of the major challenges in realizing antiferromagnetic and superfluid phases in optical lattices is the ability to cool fermions. We determine constraints on the entropy for observing these phases in two-dimensional Hubbard models using determinantal quantum Monte Carlo simulations. We find that an entropy per particle approximately = ln2 is sufficient to observe the insulating gap in the repulsive Hubbard model at half-filling, or the pairing pseudogap in the attractive case. Observing antiferromagnetic correlations or superfluidity in 2D systems requires a further reduction in entropy by a factor of 3 or more. In contrast with higher dimensions, we find that adiabatic cooling is not useful to achieve the required low temperatures. We also show that double-occupancy measurements are useful for thermometry for temperatures greater than the nearest-neighbor hopping energy.
One of the major challenges in realizing antiferromagnetic and superfluid phases in optical lattices is the ability to cool fermions. We determine constraints on the entropy for observing these phases in two-dimensional Hubbard models using determinantal quantum Monte Carlo simulations. We find that an entropy per particle approximately = ln2 is sufficient to observe the insulating gap in the repulsive Hubbard model at half-filling, or the pairing pseudogap in the attractive case. Observing antiferromagnetic correlations or superfluidity in 2D systems requires a further reduction in entropy by a factor of 3 or more. In contrast with higher dimensions, we find that adiabatic cooling is not useful to achieve the required low temperatures. We also show that double-occupancy measurements are useful for thermometry for temperatures greater than the nearest-neighbor hopping energy.
A major challenge in realizing antiferromagnetic and superfluid phases in optical lattices is the ability to cool fermions. We determine the equation of state for the 3D repulsive Fermi-Hubbard model as a function of the chemical potential, temperature, and repulsion using unbiased determinantal quantum Monte Carlo methods, and we then use the local density approximation to model a harmonic trap. We show that increasing repulsion leads to cooling but only in a trap, due to the redistribution of entropy from the center to the metallic wings. Thus, even when the average entropy per particle is larger than that required for antiferromagnetism in the homogeneous system, the trap enables the formation of an antiferromagnetic Mott phase.
Abstract The celebrated antiferromagnetic (AFM) phase transition was realized in a most recent optical lattice experiment for the 3D fermionic Hubbard model [ Shao et al . Nature 632 267 (2024) ]. Despite this important progress, it was observed that the AFM structure factor (and also the critical entropy) reaches its maximum at an interaction strength U / t ≃ 11.75, which is significantly larger than the theoretical prediction of U / t ≃ 8. Here, we resolve this discrepancy by studying the interplay between the thermal entropy, density disorder, and antiferromagnetism in the half-filled 3D Hubbard model, using numerically exact auxiliary-field quantum Monte Carlo simulations. We have achieved an accurate entropy phase diagram, enabling us to simulate arbitrary entropy path on the temperature-interaction plane and track experimental parameters effectively. We find that above the discrepancy can be quantitatively explained by the entropy increase associated with increasing interaction strength in experiments, and together with the lattice density disorder present in the experimental setup. We further investigate the entropy dependence of double occupancy and predict universal behaviors that could serve as valuable probes in future optical lattice experiments.
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.
In the experimental context of cold-fermion optical lattices, we discuss the possibilities to approach the pseudogap or ordered phases by manipulating the scattering length or the strength of the laser-induced lattice potential. Using the two-particle self-consistent approach, as well as quantum Monte Carlo simulations, we provide isentropic curves for the two- and three-dimensional Hubbard models at half-filling. These quantitative results are important for practical attempts to reach the ordered antiferromagnetic phase in experiments on optical lattices of two-component fermions. We find that adiabatically turning on the interaction in two dimensions to cool the system is not very effective. In three dimensions, adiabatic cooling to the antiferromagnetic phase can be achieved in such a manner, although the cooling efficiency is not as high as initially suggested by dynamical mean-field theory. Adiabatic cooling by turning off the repulsion beginning at strong coupling is possible in certain cases.
Antiferromagnetic, Mott-insulator, d-wave, and gossamer superfluid phases are calculated for two-dimensional (2D) square lattices from the extended Hubbard $(t\text{\ensuremath{-}}J\text{\ensuremath{-}}U)$ model using the Gutzwiller projection method and renormalized mean-field theory. The phase separation between antiferromagnetic and d-wave superfluid phases is found near half filling when the on-site repulsion exceeds $U\ensuremath{\gtrsim}7.3t$ and coincides with a first-order transition in the double occupancy. The phase separation is thus predicted for 2D optical lattices with ultracold Fermi atoms whereas it is inhibited in cuprates by Coulomb frustration which instead may lead to stripes. In a confined optical lattice the resulting density distribution is discontinuous and with extended Mott plateau which enhances the antiferromagnetic phase but suppresses the superfluid phase. Observation of Mott-insulator, antiferromagnetic, stripe, and superfluid phases in density and momentum distributions and correlations is discussed.
We study finite-temperature properties of strongly interacting fermions in the honeycomb lattice using numerical linked-cluster expansions and determinantal quantum Monte Carlo simulations. We analyze a number of thermodynamic quantities, including the entropy, the specific heat, uniform and staggered spin susceptibilities, short-range spin correlations, and the double occupancy at and away from half filling. We examine the viability of adiabatic cooling by increasing the interaction strength for homogeneous as well as for trapped systems. For the homogeneous case, this process is found to be more efficient at finite doping than at half filling. That, in turn, leads to an efficient adiabatic cooling in the presence of a trap, which, starting with even relatively high entropies, can drive the system to have a Mott insulating phase with substantial antiferromagnetic correlations.