Pomeranchuk cooling of SU(2N) ultracold fermions in optical lattices

Emergence of Competing Orders and Possible Quantum Spin Liquid in SU(N) Fermions

In the past few decades, tremendous efforts have been made toward understanding the exotic physics emerging from competition between various ordering tendencies in strongly correlated systems. Employing state-of-the-art quantum Monte Carlo simulation, we investigate an interacting SU(N) fermionic model with varying interaction strength and value of N, and we unveil the ground-state phase diagram of the model exhibiting a plethora of exotic phases. For small values of N-namely, N=2, 3-the ground state is an antiferromagnetic (AFM) phase, whereas in the large-N limit, a staggered valence bond solid (VBS) order is dominant. For intermediate values of N such as N=4, 5, remarkably, our study reveals that distinct VBS orders appear in the weak and strong coupling regimes. More fantastically, the competition between staggered and columnar VBS ordering tendencies gives rise to a Mott insulating phase without spontaneous symmetry breaking (SSB), existing in a large interacting parameter regime, which is consistent with a gapped quantum spin liquid. Our study not only provides a platform to investigate the fundamental physics of quantum many-body systems-it also offers a novel route toward searching for exotic states of matter such as quantum spin liquid in realistic quantum materials.

Tangent Space Approach for Thermal Tensor Network Simulations of the 2D Hubbard Model

Accurate simulations of the two-dimensional (2D) Hubbard model constitute one of the most challenging problems in condensed matter and quantum physics. Here we develop a tangent space tensor renormalization group (tanTRG) approach for the calculations of the 2D Hubbard model at finite temperature. An optimal evolution of the density operator is achieved in tanTRG with a mild O(D^{3}) complexity, where the bond dimension D controls the accuracy. With the tanTRG approach we boost the low-temperature calculations of large-scale 2D Hubbard systems on up to a width-8 cylinder and 10×10 square lattice. For the half-filled Hubbard model, the obtained results are in excellent agreement with those of determinant quantum Monte Carlo (DQMC). Moreover, tanTRG can be used to explore the low-temperature, finite-doping regime inaccessible for DQMC. The calculated charge compressibility and Matsubara Green's function are found to reflect the strange metal and pseudogap behaviors, respectively. The superconductive pairing susceptibility is computed down to a low temperature of approximately 1/24 of the hopping energy, where we find d-wave pairing responses are most significant near the optimal doping. Equipped with the tangent-space technique, tanTRG constitutes a well-controlled, highly efficient and accurate tensor network method for strongly correlated 2D lattice models at finite temperature.

Interaction Effects with Varying N in SU(N) Symmetric Fermion Lattice Systems

The interaction effects in ultracold Fermi gases with SU(N) symmetry are studied nonperturbatively in half filled one-dimensional lattices by employing quantum Monte Carlo simulations. We find that, as N increases, weak and strong interacting systems are driven to a crossover region, but from opposite directions as a convergence of itinerancy and Mottness. In the weak interaction region, particles are nearly itinerant, and interparticle collisions are enhanced by N, resulting in the amplification of interaction effects. In contrast, in the strong coupling region, increasing N softens the Mott-insulating background through the enhanced virtual hopping processes. The crossover region exhibits nearly N-independent physical quantities, including the relative bandwidth, Fermi distribution, and the spin structure factor. The difference between even-N and odd-N systems is most prominent at small N's with strong interactions, since the odd case allows local real hopping with an energy scale much larger than the virtual one. The above effects can be experimentally tested in ultracold atom experiments with alkaline-earth(-like) fermions such as ^{87}Sr (^{173}Yb).

Fermion-induced quantum critical points

A unified theory of quantum critical points beyond the conventional Landau-Ginzburg-Wilson paradigm remains unknown. According to Landau cubic criterion, phase transitions should be first-order when cubic terms of order parameters are allowed by symmetry in the Landau-Ginzburg free energy. Here, from renormalization group analysis, we show that second-order quantum phase transitions can occur at such putatively first-order transitions in interacting two-dimensional Dirac semimetals. As such type of Landau-forbidden quantum critical points are induced by gapless fermions, we call them fermion-induced quantum critical points. We further introduce a microscopic model of SU(N) fermions on the honeycomb lattice featuring a transition between Dirac semimetals and Kekule valence bond solids. Remarkably, our large-scale sign-problem-free Majorana quantum Monte Carlo simulations show convincing evidences of a fermion-induced quantum critical points for N = 2, 3, 4, 5 and 6, consistent with the renormalization group analysis. We finally discuss possible experimental realizations of the fermion-induced quantum critical points in graphene and graphene-like materials.Quantum phase transitions are governed by Landau-Ginzburg theory and the exceptions are rare. Here, Li et al. propose a type of Landau-forbidden quantum critical points induced by gapless fermions in two-dimensional Dirac semimetals.

Majorana Positivity and the Fermion Sign Problem of Quantum Monte Carlo Simulations

The sign problem is a major obstacle in quantum Monte Carlo simulations for many-body fermion systems. We examine this problem with a new perspective based on the Majorana reflection positivity and Majorana Kramers positivity. Two sufficient conditions are proven for the absence of the fermion sign problem. Our proof provides a unified description for all the interacting lattice fermion models previously known to be free of the sign problem based on the auxiliary field quantum Monte Carlo method. It also allows us to identify a number of new sign-problem-free interacting fermion models including, but not limited to, lattice fermion models with repulsive interactions but without particle-hole symmetry, and interacting topological insulators with spin-flip terms.

Ultracold Fermi gases with emergent SU(N) symmetry

We review recent experimental and theoretical progress on ultracold alkaline-earth Fermi gases with emergent SU(N) symmetry. Emphasis is placed on describing the ground-breaking experimental achievements of recent years. The latter include (1) the cooling to below quantum degeneracy of various isotopes of ytterbium and strontium, (2) the demonstration of optical Feshbach resonances and the optical Stern-Gerlach effect, (3) the realization of a Mott insulator of (173)Yb atoms, (4) the creation of various kinds of Fermi-Bose mixtures and (5) the observation of many-body physics in optical lattice clocks. On the theory side, we survey the zoo of phases that have been predicted for both gases in a trap and loaded into an optical lattice, focusing on two and three dimensional systems. We also discuss some of the challenges that lie ahead for the realization of such phases such as reaching the temperature scale required to observe magnetic and more exotic quantum orders. The challenge of dealing with collisional relaxation of excited electronic levels is also discussed.

Competing orders in the 2D half-filled SU(2N) Hubbard model through the pinning-field quantum Monte Carlo simulations

We nonperturbatively investigate the ground state magnetic properties of the 2D half-filled SU(2N) Hubbard model in the square lattice by using the projector determinant quantum Monte Carlo simulations combined with the method of local pinning fields. Long-range Néel orders are found for both the SU(4) and SU(6) cases at small and intermediate values of U. In both cases, the long-range Néel moments exhibit nonmonotonic behavior with respect to U, which first grow and then drop as U increases. This result is fundamentally different from the SU(2) case in which the Néel moments increase monotonically and saturate. In the SU(6) case, a transition to the columnar dimer phase is found in the strong interaction regime.

Dimerized solids and resonating plaquette order in SU(N)-Dirac fermions

We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our numerical results to a saddle-point approximation in the large-N limit. From the large-N regime down to the SU(6) case, the insulating state is found to be a columnar valence bond crystal, with a direct transition to the semimetal at weak, finite coupling, in agreement with the mean-field result in the large-N limit. At SU(4) however, the insulator exhibits a subtly different valence bond crystal structure, stabilized by resonating valence bond plaquettes. In the SU(2) limit, our results support a direct transition between the semimetal and an antiferromagnetic insulator.