We use determinantal quantum Monte Carlo simulations and numerical linked-cluster expansions to study thermodynamic properties and short-range spin correlations of fermions in the honeycomb lattice. We find that, at half filling and finite temperatures, nearest-neighbor spin correlations can be stronger in this lattice than in the square lattice, even in regimes where the ground state in the former is a semimetal or a spin liquid. The honeycomb lattice also exhibits a more pronounced anomalous region in the double occupancy that leads to stronger adiabatic cooling than in the square lattice. We discuss the implications of these findings for optical lattice experiments.
We use determinantal quantum Monte Carlo simulations and numerical linked-cluster expansions to study thermodynamic properties and short-range spin correlations of fermions in the honeycomb lattice. We find that, at half filling and finite temperatures, nearest-neighbor spin correlations can be stronger in this lattice than in the square lattice, even in regimes where the ground state in the former is a semimetal or a spin liquid. The honeycomb lattice also exhibits a more pronounced anomalous region in the double occupancy that leads to stronger adiabatic cooling than in the square lattice. We discuss the implications of these findings for optical lattice experiments.
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.
We use numerical linked cluster expansions (NLCEs) to study correlated lattice systems in one and two dimensions.For fermions in the honeycomb lattice, we study its finite-temperature properties and short-range spin correlations using NLCEs and determinantal quantum Monte Carlo (DQMC) simulations.For the homogeneous system, 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.Employing a local density approximation (LDA), we examine the viability of adiabatic cooling by increasing the interaction strength for homogeneous as well as for trapped systems.We also use NLCEs to study thermodynamic properties of the two-dimensional spin-1/2 Ising, XY, and Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattices.In all cases, the nearest-neighbor coupling between the spins takes values ±J with equal probability.We obtain the disorder averaged (over all disorder configurations) energy, entropy, specific heat, and uniform magnetic susceptibility in each case.These results are compared with the corresponding ones in the clean models.Analytic expressions are obtained for low orders in the expansion of these thermodynamic quantities in inverse temperature.For many-body localization in disordered isolated systems, we study quantum quenches in the thermodynamic limit.By a quantum quench it is meant that the initial state is stationary with respect iii to an initial Hamiltonian, which is suddenly changed to a new (time-independent) Hamiltonian.The latter then drives the (unitary) dynamics of the system.We are interested in the time average of observables after relaxation following the quench, which can be computed under the so-called diagonal ensemble.
We study finite-temperature properties of strongly correlated fermions in two-dimensional optical lattices by means of numerical linked cluster expansions, a computational technique that allows one to obtain exact results in the thermodynamic limit. We focus our analysis on the strongly interacting regime, where the on-site repulsion is of the order of or greater than the band width. We compute the equation of state, double occupancy, entropy, uniform susceptibility, and spin correlations for temperatures that are similar to or below the ones achieved in current optical lattice experiments. We provide a quantitative analysis of adiabatic cooling of trapped fermions in two dimensions, by means of both flattening the trapping potential and increasing the interaction strength.
We investigate the thermodynamic properties of a half-filled SU(2N) Hubbard model in the two-dimensional square lattice by the method of the determinant quantum Monte Carlo simulation, which is free of the fermion "sign problem." The large number of hyperfine-spin components enhances spin fluctuations, which facilitates the Pomeranchuk cooling to temperatures comparable to the superexchange energy scale in the case of SU(6). Various physical quantities including entropy, charge fluctuations, and spin correlations are calculated.
Optical lattice experiments with ultracold fermion atoms and quantum gas microscopy have recently realized direct measurements of magnetic correlations at the site-resolved level. We calculate the short-range spin-correlation functions in the ground state of the two-dimensional repulsive Hubbard model with the auxiliary-field quantum Monte Carlo (AFQMC) method. The results are numerically exact at half filling where the fermion sign problem is absent. Away from half filling, we employ the constrained path AFQMC approach to eliminate the exponential computational scaling from the sign problem. The constraint employs unrestricted Hartree-Fock trial wave functions with an effective interaction strength U, which is optimized self-consistently within AFQMC. Large supercells are studied, with twist averaged boundary conditions as needed, to reach the thermodynamic limit. We find that the nearest-neighbor spin correlation always increases with the interaction strength U, contrary to the finite-temperature behavior where a maximum is reached at a finite U value. We also observe a change of sign in the next-nearest-neighbor spin correlation with increasing density, which is a consequence of the buildup of the long-range antiferromagnetic correlation. We expect the results presented in this paper to serve as a benchmark as lower temperatures are reached in ultracold atom experiments.
The half-filled Hubbard model on the honeycomb lattice is investigated by numerically exact large-scale quantum Monte Carlo simulations for lattice sizes up to 2592 sites. By performing careful finite-size scaling for the spin-spin correlation functions, calculated with a high degree of accuracy, we find that the ground state is antiferromagnetically long-range ordered at U / t = 4, where U is the on-site Hubbard interaction and t is the nearest neighbor hopping. Our result is in sharp contrast to a recent report [Meng et al. , Nature 464, 847 (2010)], where instead strong evidence of spin liquid behavior is found.