No Arabic abstract
We study thermodynamics of the 3D Hubbard model at half filling on approach to the Neel transition by means of large-scale unbiased Diagrammatic Determinant Monte Carlo simulations. We obtain the transition temperature in the strongly correlated regime, as well as temperature dependence of energy, entropy, double occupancy, and the nearest-neighbor spin correlation function. Our results improve the accuracy of previous unbiased studies and present accurate benchmarks in the ongoing effort to realize the antiferromagnetic state of matter with ultracold atoms in optical lattices.
We study the 3D Hubbard model with anisotropic nearest neighbor tunneling amplitudes using the dynamical cluster approximation and compare the results with a quantum simulation experiment using ultracold fermions in an optical lattice, focussing on magnetic correlations. We find that the short-range spin correlations are significantly enhanced in the direction with stronger tunneling amplitudes. Our results agree with the experimental observations and show that the experimental temperature is lower than the strong tunneling amplitude. We characterize the system by examining the spin correlations beyond neighboring sites and determine the distribution of density, entropy and spin correlation in the trapped system. We furthermore investigate the dependence of the critical entropy at the Neel transition on anisotropy.
The realization of antiferromagnetic (AF) correlations in ultracold fermionic atoms on an optical lattice is a significant achievement. Experiments have been carried out in one, two, and three dimensions, and have also studied anisotropic configurations with stronger tunneling in some lattice directions. Such anisotropy is relevant to the physics of cuprate superconductors and other strongly correlated materials. Moreover, this anisotropy might be harnessed to enhance AF order. Here we numerically investigate, using Determinant Quantum Monte Carlo, a simple realization of anisotropy in the 3D Hubbard model in which the tunneling between planes, $t_perp$, is unequal to the intraplane tunneling $t$. This model interpolates between the three-dimensional isotropic ($t_perp = t$) and two-dimensional ($t_perp =0$) systems. We show that at fixed interaction strength to tunneling ratio ($U/t$), anisotropy can enhance the magnetic structure factor relative to both 2D and 3D results. However, this enhancement occurs at interaction strengths below those for which the Neel temperature $T_{rm Nacute{e}el}$ is largest, in such a way that the structure factor cannot be made to exceed its value in isotropic 3D systems at the optimal $U/t$. We characterize the 2D-3D crossover in terms of the magnetic structure factor, real space spin correlations, number of doubly-occupied sites, and thermodynamic observables. An interesting implication of our results stems from the entropys dependence on anisotropy. As the system evolves from 3D to 2D, the entropy at a fixed temperature increases. Correspondingly, at fixed entropy, the temperature will decrease going from 3D to 2D. This suggests a cooling protocol in which the dimensionality is adiabatically changed from 3D to 2D.
We present a simple trick that allows to consider the sum of all connected Feynman diagrams at fixed position of interaction vertices for general fermionic models. With our approach one achieves superior performance compared to Diagrammatic Monte Carlo, while rendering the algorithmic part dramatically simpler. As we consider the sum of all connected diagrams at once, we allow for cancellations between diagrams with different signs, alleviating the sign problem. Moreover, the complexity of the calculation grows exponentially with the order of the expansion, which should be constrasted with the factorial growth of the standard diagrammatic technique. We illustrate the efficiency of the technique for the two-dimensional Fermi-Hubbard model.
We investigate the existence of quantum disentangled liquid (QDL) states in the half-filled Hubbard model on bipartite lattices. In the one dimensional case we employ a combination of integrability and strong coupling expansion methods to argue that there are indeed finite energy-density eigenstates that exhibit QDL behaviour in the sense of J. Stat. Mech. P10010 (2014). The states exhibiting the QDL property are atypical in the sense that while their entropy density is non-zero, it is smaller than that of thermal states at the same energy density. We argue that for U >> t these latter thermal states exhibit a weaker form of the QDL property, which carries over to the higher dimensional case.
The SU(2) symmetric Fermi-Hubbard model (FHM) plays an essential role in strongly correlated fermionic many-body systems. In the one particle per site and strongly interacting limit ${U/t gg 1}$, it is effectively described by the Heisenberg Hamiltonian. In this limit, enlarging the spin and extending the typical SU(2) symmetry to SU($N$) has been predicted to give exotic phases of matter in the ground state, with a complicated dependence on $N$. This raises the question of what --- if any --- are the finite-temperature signatures of these phases, especially in the currently experimentally relevant regime near or above the superexchange energy. We explore this question for thermodynamic observables by numerically calculating the thermodynamics of the SU($N$) FHM in the two-dimensional square lattice near densities of one particle per site, using determinant Quantum Monte Carlo and Numerical Linked Cluster Expansion. Interestingly, we find that for temperatures above the superexchange energy, where the correlation length is short, the energy, number of on-site pairs, and kinetic energy are universal functions of $N$. Although the physics in the regime studied is well beyond what can be captured by low-order high-temperature series, we show that an analytic description of the scaling is possible in terms of only one- and two-site calculations.