Zero-temperature Phase Diagram Research Articles

Lieb-Schultz-Mattis anomalies and web of dualities induced by gauging in quantum spin chains

Lieb-Schultz-Mattis (LSM) theorems impose non-perturbative constraints on the zero-temperature phase diagrams of quantum lattice Hamiltonians (always assumed to be local in this paper). LSM theorems have recently been interpreted as the lattice counterparts to mixed ’t Hooft anomalies in quantum field theories that arise from a combination of crystalline and global internal symmetry groups. Accordingly, LSM theorems have been reinterpreted as LSM anomalies. In this work, we provide a systematic diagnostic for LSM anomalies in one spatial dimension. We show that gauging subgroups of the global internal symmetry group of a quantum lattice model obeying an LSM anomaly delivers a dual quantum lattice Hamiltonian such that its internal and crystalline symmetries mix non-trivially through a group extension. This mixing of crystalline and internal symmetries after gauging is a direct consequence of the LSM anomaly, i.e., it can be used as a diagnostic of an LSM anomaly. We exemplify this procedure for a quantum spin-1/2 chain obeying an LSM anomaly resulting from combining a global internal \mathbb{Z}^{\,}_{2}×\mathbb{Z}^{\,}_{2}ℤ2×ℤ2 symmetry with translation or reflection symmetry. We establish a triality of models by gauging a \mathbb{Z}^{\,}_{2}\subset\mathbb{Z}^{\,}_{2}×\mathbb{Z}^{\,}_{2}ℤ2⊂ℤ2×ℤ2 symmetry in two ways, one of which amounts to performing a Kramers-Wannier duality, while the other implements a Jordan-Wigner duality. We discuss the mapping of the phase diagram of the quantum spin-1/2 XYZXYZ chains under such a triality. We show that the deconfined quantum critical transitions between Neel and dimer orders are mapped to either topological or conventional Landau-Ginzburg transitions. Finally, we extend our results to \mathbb{Z}^{\,}_{n}ℤn clock models with \mathbb{Z}^{\,}_{n}×\mathbb{Z}^{\,}_{n}ℤn×ℤn global internal symmetry, and provide a reinterpretation of the dual internal symmetries in terms of \mathbb{Z}^{\,}_{n}ℤn charge and dipole symmetries.

Eleven competing phases in the Heisenberg-Gamma ( JΓ ) ladder

The spin-orbit generated Γ interaction is known to induce strong frustration and to be significant in realistic models of materials. To gain an understanding of the possible phases that can arise from this interaction, it is of considerable interest to focus on a limited part of parameter space in a quasi one-dimensional model where high precision numerical results can be obtained. Here we study the Heisenberg–Gamma ( JΓ ) ladder, determining the complete zero temperature phase diagram by analyzing the entanglement spectrum (ES) and energy susceptibility. A total of 11 different phases can be identified, among them the well known rung-singlet (RS) phase and 5 other phases, FM, FM-Z, FM-XY, AF and AF-Z, with conventional long-range magnetic order. The 3 ferromagnetic phases, FM, FM-Z and FM-XY simultaneously have non-zero scalar chirality. Two other phases, the antiferromagnetic Gamma (AΓ) and ferromagnetic Gamma (FΓ) phases, have previously been observed in the Kitaev–Gamma ladder, demonstrating that the AΓ-phase is a symmetry protected topological phase (SPT) protected by TR ×Rb symmetry, the product of time-reversal (TR) and π rotation around the b-axis ( Rb ), while the FΓ-phase is related to the RS phase through a local unitary transformation. The 3 remaining phases, ϒ, Ω and δ show no conventional order, a doubling of the ES and for the ϒ and Ω-phases a gap is clearly present. The δ-phase has a significantly smaller gap and displays incommensurate correlations, with a peak in the static structure factor, S(k) continuously shifting from k/π=2/3 to k=π . In the Ω-phase we find pronounced edge-states consistent with a SPT phase protected by the same TR ×Rb symmetry as the AΓ-phase. The precise nature of the ϒ and δ-phases is less clear.

Quantum Phases from Competing Van der Waals and Dipole-Dipole Interactions of Rydberg Atoms.

Competing short- and long-range interactions represent distinguished ingredients for the formation of complex quantum many-body phases. Their study is hard to realize with conventional quantum simulators. In this regard, Rydberg atoms provide an exception as their excited manifold of states have both density-density and exchange interactions whose strength and range can vary considerably. Focusing on one-dimensional systems, we leverage the Van der Waals and dipole-dipole interactions of the Rydberg atoms to obtain the zero-temperature phase diagram for a uniform chain and a dimer model. For the uniform chain, we can influence the boundaries between ordered phases and a Luttinger liquid phase. For the dimerized case, a new type of bond-order-density-wave phase is identified. This demonstrates the versatility of the Rydberg platform in studying physics involving short- and long-ranged interactions simultaneously.

Phases of cold holographic QCD: Baryons, pions and rho mesons

We improve the holographic description of isospin-asymmetric baryonic matter within the Witten-Sakai-Sugimoto model by accounting for a realistic pion mass, computing the pion condensate dynamically, and including rho meson condensation by allowing the gauge field in the bulk to be anisotropic. This description takes into account the coexistence of baryonic matter with pion and rho meson condensates. Our main result is the zero-temperature phase diagram in the plane of baryon and isospin chemical potentials. We find that the effective pion mass in the baryonic medium increases with baryon density and that, as a consequence, there is no pion condensation in neutron-star matter. Our improved description also predicts that baryons are disfavored at low baryon chemical potentials even for arbitrarily large isospin chemical potential. Instead, rho meson condensation sets in on top of the pion condensate at an isospin chemical potential of about 9.4\, m_\pi9.4mπ. We further observe a highly non-monotonic phase boundary regarding the disappearance of pion condensation.

Mean-field approximation of the O(2) model on 3D lattice Bose systems

New results are presented for bosons defined on the three-dimensional (3D) cubic lattice at very low temperatures. Calculations are made within the mean-filed (or Hartree) theory applied to the so-called O(2) model. The analysis is used to calculate the zero-temperature phase diagram and determine the critical coupling strength λ for a quantum phase transition. Properties of the zero temperature quantum state, such as: the energy per site, the order parameter, the first possible excitation, and the test single-site phase wave function, are displayed and discussed. Their behavior predicts a transition from an ordered (superfluid) to a disordered (normal/insulating) phase at a critical value of the coupling constant λc = 6.

Quench dynamics in the one-dimensional mass-imbalanced ionic Hubbard model

Using the time-dependent Lanczos method, we study the nonequilibrium dynamics of the one-dimensional ionic-mass imbalanced Hubbard chain driven by a quantum quench of the on-site Coulomb interaction, where the system is prepared in the ground state of the Hamiltonian with a different Hubbard interaction. The exact diagonalization method (Lanczos algorithm) is adopted to study the zero temperature phase diagram in equilibrium, which is shown to be in good agreement with previous studies using density matrix renormalization group (DMRG). We then study the nonequilibrium quench dynamics of the spin and charge order parameters by fixing the initial and final Coulomb interactions while changing the quenching time protocols. Our study shows that the time evolution of the charge and spin order parameters strongly depend on the quenching time protocols. In particular, the effective temperature of the system will decrease monotonically as the quenching time is increased. By taking the final Coulomb interaction strength to be in the strong coupling regime, we find that the oscillation frequency of the charge order parameter increases monotonically with the Coulomb interaction. By contrast, the frequency of the spin order parameter decreases monotonically with increasing Coulomb interaction. We explain this result using an effective spin model and a two-site Hubbard model in the strong coupling limit. Finally, we take the final Coulomb interaction strength to be in the weak coupling regime and find that the oscillation frequency of both the charge and spin order parameters increases monotonically with decreasing Coulomb interaction. Our study suggests strategies to engineer the relaxation behavior of interacting quantum many-particle systems.

Even-odd effects in the J1−J2 SU(N) Heisenberg spin chain

The zero-temperature phase diagram of the ${J}_{1}\text{\ensuremath{-}}{J}_{2}\phantom{\rule{4pt}{0ex}}\mathrm{SU}(N)$ antiferromagnetic Heisenberg spin chain is investigated by means of complementary field theory and numerical approaches for general $N$. A fully gapped $\mathrm{SU}(N)$ valence bond solid made of $N$ sites is formed above a critical value of ${J}_{2}/{J}_{1}$ for all $N$. We find that the extension of this $N$-merized phase for larger values of ${J}_{2}$ strongly depends on the parity of $N$. For even $N$, the phase smoothly interpolates to the large ${J}_{2}$ regime where the model can be viewed as a zigzag $\mathrm{SU}(N)$ two-leg spin ladder. The phase exhibits both a $N$-merized ground state and incommensurate spin-spin correlations. In stark contrast to the even case, we show that the $N$-merized phase with odd $N$ has only a finite extent with no incommensuration. A gapless phase in the $\mathrm{SU}{(N)}_{1}$ universality class is stabilized for larger ${J}_{2}$ that stems from the existence of a massless renormalization group flow from $\mathrm{SU}{(N)}_{2}$ to $\mathrm{SU}{(N)}_{1}$ conformal field theories when $N$ is odd.

Quadriexcitons and excitonic condensate in a symmetric electron-hole bilayer with valley degeneracy

Using quantum Monte Carlo simulations we have mapped out the zero-temperature phase diagram of a symmetric electron-hole bilayer with twofold valley degeneracy, as a function of the interlayer distance $d$ and in-layer density $n$. We find that the effect of the valley degeneracy is to shrink the region of stability of the excitonic condensate, in favor of quadriexcitons at small $d$ and of the four-component plasma at large $d$, with minor effects on the value of the excitonic condensate fraction. The enclosure of the condensate in a density window possibly explains why anomalous tunneling conductivity, interpreted as a signature of condensation, is observed only between two finite values of carrier density in graphene bilayers. Our phase diagram may provide directions to select device parameters for future experiments.

Trimer states with Z3 topological order in Rydberg atom arrays

Trimers are defined as two adjacent edges on a graph. We study the quantum states obtained as equal-weight superpositions of all trimer coverings of a lattice, with the constraint of having a trimer on each vertex: the so-called trimer resonating-valence-bond (tRVB) states. Exploiting their tensor network representation, we show that these states can host $\mathbb{Z}_3$ topological order or can be gapless liquids with $\mathrm{U}(1) \times \mathrm{U}(1)$ local symmetry. We prove that this continuous symmetry emerges whenever the lattice can be tripartite such that each trimer covers all the three sublattices. In the gapped case, we demonstrate the stability of topological order against dilution of maximal trimer coverings, which is relevant for realistic models where the density of trimers can fluctuate. Furthermore, we clarify the connection between gapped tRVB states and $\mathbb{Z}_3$ lattice gauge theories by smoothly connecting the former to the $\mathbb{Z}_3$ toric code, and discuss the non-local excitations on top of tRVB states. Finally, we analyze via exact diagonalization the zero-temperature phase diagram of a diluted trimer model on the square lattice and demonstrate that the ground state exhibits topological properties in a narrow region in parameter space. We show that a similar model can be implemented in Rydberg atom arrays exploiting the blockade effect. We investigate dynamical preparation schemes in this setup and provide a viable route for probing experimentally $\mathbb{Z}_3$ quantum spin liquids.

Quantum skyrmion lattices in Heisenberg ferromagnets

Skyrmions are topological magnetic textures that can arise in non-centrosymmetric ferromagnetic materials. In most systems experimentally investigated to date, skyrmions emerge as classical objects. However, the discovery of skyrmions with nanometer length scales has sparked interest in their quantum properties. Here, we simulate the ground states of two-dimensional spin-$1/2$ Heisenberg lattices with Dzyaloshinskii-Moriya interactions and discover a broad region in the zero-temperature phase diagram which hosts quantum skyrmion lattices. We argue that the quantum skyrmion lattice phase can be detected experimentally in the magnetization profile via local magnetic polarization measurements as well as in the spin structure factor measurable via neutron scattering experiments. Finally, we explore the resulting quantum skyrmion state, analyze its real-space polarization profile and show that it is a non-classical state featuring entanglement between quasiparticle and environment mainly localized near the boundary spins of the skyrmion.

Investigation of the temperature phase diagram of the 2D semi-classical ferrimagnetic model

The temperature phase diagram (T over J 2) of the 2D semiclassical Heisenberg model is discussed. The investigation was carried out using the self-consistent Bloch equation and a quantum statistical description of the magnetization. In addition to the three states of the zero-temperature phase diagram, the phase diagram shows two new phases, the paramagnetic phase and the semi-paramagnetic -Néel state. The five occurring ground states, as well as the transitions between these states, are discussed. In particular, attention is paid to the interplay between temperature and frustration.

Benchmark calculations of multiloop pseudofermion fRG

The pseudofermion functional renormalization group (pffRG) is a computational method for determining zero-temperature phase diagrams of frustrated quantum magnets. In a recent methodological advance, the commonly employed Katanin truncation of the flow equations was extended to include multiloop corrections, thereby capturing additional contributions from the three-particle vertex (Thoenniss et al. https://arxiv.org/abs/2011.01268; Kiese et al. https://arxiv.org/abs/2011.01269). This development has also stimulated significant progress in the numerical implementation of pffRG, allowing one to track the evolution of pseudofermion vertices under the renormalization group flow with unprecedented accuracy. However, cutting-edge solvers differ in their integration algorithms, heuristics to discretize Matsubara frequency grids, and more. To lend confidence in the numerical robustness of state-of-the-art multiloop pffRG codes, we present and compare results produced with two independently developed and algorithmically distinct solvers for Heisenberg models on three-dimensional lattice geometries. Using the cubic lattice Heisenberg (anti)ferromagnet with nearest and next-nearest neighbor interactions as a generic benchmark model, we find the two codes to quantitatively agree, often up to several orders of magnitude in digital precision, both on the level of spin-spin correlation functions and renormalized fermionic vertices for varying loop orders. These benchmark calculations further substantiate the usage of multiloop pffRG solvers to tackle unconventional forms of quantum magnetism.Graphic abstract

Low-density phase diagram of the three-dimensional electron gas

Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of the three-dimensional homogeneous electron gas at very low density. Fermi fluid and body-centered cubic Wigner crystal ground state energies are determined using Slater-Jastrow-backflow and Slater-Jastrow many-body wave functions at different densities and spin polarizations in finite simulation cells. Finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Unlike previous studies, our results show that the electron gas undergoes a first-order quantum phase transition directly from a paramagnetic fluid to a body-centered cubic crystal at density parameter $r_\text{s} = 86.6(7)$, with no region of stability for an itinerant ferromagnetic fluid. However there is a possible magnetic phase transition from an antiferromagnetic crystal to a ferromagnetic crystal at $r_\text{s}=93(3)$.

Quantum orders in the frustrated Ising model on the bathroom tile lattice

We determine the zero and finite temperature phase diagram of the fully frustrated quantum Ising model on the bathroom tile (4-8) lattice. The phase diagram exhibits a wealth of 2+1d physics, including 1. classical Coulomb dimer liquids of both square and triangular lattice types; 2. quantum order-by-disorder induced phases breaking $\mathbb{Z}_4$, $\mathbb{Z}_6$, and $\mathbb{Z}_8$ symmetries; 3. finite temperature Kosterlitz-Thouless (KT) phases floating over the $\mathbb{Z}_6$ and $\mathbb{Z}_8$ orders; and, 4. staircases of (in)-commensurate symmetry breaking phases at intermediate coupling. We establish this elaborate phase diagram using a combination of dimer model mapping, perturbation theory, Landau analysis and Stochastic Series Expansion Quantum Monte Carlo (QMC-SSE). Our results provide a baseline for studying frustrated magnetism with D-Wave architecture annealers, where the 4-8 lattice can be embedded naturally without `cloning', reducing the number of competing energy scales. Simulations with the D-Wave 2000Q demonstrate qualitative agreement with the high temperature portion of the phase diagram, but are unable to access the low temperature phases.

Frustrated spin-[formula omitted] Heisenberg magnet on an [formula omitted]-stacked honeycomb bilayer: High-order study of the collinear magnetic phases of the [formula omitted]–[formula omitted]–[formula omitted] model

The zero-temperature phase diagram of the frustrated spin-12J1–J2–J1⊥ Heisenberg magnet on an AA-stacked honeycomb bilayer lattice is studied using the coupled cluster method implemented to very high orders. On each monolayer the spins interact via nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic Heisenberg interactions with respective strength parameters J1>0 and J2≡κJ1>0. The two layers are coupled such that NN interlayer pairs of spins also interact via a similar isotropic Heisenberg interaction of strength J1⊥≡δJ1, which may be of either sign. In particular, we locate with high accuracy the complete phase boundaries in the κ–δ half-plane with κ>0 of the two quasiclassical collinear antiferromagnetic phases with Néel or Néel-II magnetic order in each monolayer, and the interlayer NN pairs of spins either aligned (for δ<0) or anti-aligned (for δ>0) to one another. Compared to the two-sublattice Néel order, in which all NN intralayer pairs of spins are antiparallel to one another, the four-sublattice Néel-II order is characterized by NN intralayer pairs of spins on the honeycomb lattice being antiparallel to one another along zigzag (or sawtooth) chains in a specified direction from among the three equivalent honeycomb-lattice directions, and parallel to one another for the corresponding interchain pairs.

Gapless quantum spin liquid and global phase diagram of the spin-1/2 [formula omitted] square antiferromagnetic Heisenberg model

The nature of the zero-temperature phase diagram of the spin-1/2J1-J2 Heisenberg model on a square lattice has been debated in the past three decades, and it remains one of the fundamental problems unsettled in the study of quantum many-body theory. By using the state-of-the-art tensor network method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram of the J1-J2 Heisenberg model up to 24×24 sites, we provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region 0.56≲J2/J1≤0.61 before the system enters the well known collinear antiferromagnetic phase. We stress that we make the first detailed comparison between the results of PEPS and the well-established density matrix renormalization group (DMRG) method through one-to-one direct benchmark for small system sizes, and thus give rise to a very solid PEPS calculation beyond DMRG. Our numerical evidences explicitly demonstrate the huge power of PEPS for highly frustrated spin systems. Finally, an effective field theory is also proposed to understand the physical nature of the discovered gapless QSL and its relation to deconfined quantum critical point (DQCP).

Rise and fall of plaquette order in the Shastry-Sutherland magnet revealed by pseudofermion functional renormalization group

The Shastry-Sutherland (SS) model as a canonical example of frustrated magnetism has been extensively studied. The conventional wisdom has been that the transition from the plaquette valence bond order to the Neel order is direct and potentially realizes a deconfined quantum critical point beyond the Ginzburg-Landau paradigm. This scenario, however, was challenged recently by improved numerics from density matrix renormalization group which offers evidence for a narrow gapless spin liquid between the two phases. Prompted by this controversy and to shed light on this intricate parameter regime from a fresh perspective, we report high-resolution functional renormalization group analysis of the generalized SS model. The flows of over 50 million running couplings provide a detailed picture for the evolution of spin correlations as the frequency/energy scale is dialed from the ultraviolet to the infrared to yield the zero-temperature phase diagram. The singlet dimer phase emerges as a fixed point, the Neel order is characterized by divergence in the vertex function, while the transition into and out of the plaquette order is accompanied by pronounced peaks in the plaquette susceptibility. The plaquette order is suppressed before the onset of the Neel order, lending evidence for a finite spin liquid region for ${J}_{1}/{J}_{2}\ensuremath{\in}(0.77,0.82)$, where the flow is continuous without any indication of divergence.

Phase diagrams of SO( N ) Majorana-Hubbard models: Dimerization, internal symmetry breaking, and fluctuation-induced first-order transitions

We study the zero-temperature phase diagrams of Majorana-Hubbard models with SO($N$) symmetry on two-dimensional honeycomb and $\pi$-flux square lattices, using mean-field and renormalization group approaches. The models can be understood as real counterparts of the SU($N$) Hubbard-Heisenberg models, and may be realized in Abrikosov vortex phases of topological superconductors, or in fractionalized phases of strongly-frustrated spin-orbital magnets. In the weakly-interacting limit, the models feature stable and fully symmetric Majorana semimetal phases. Increasing the interaction strength beyond a finite threshold for large $N$, we find a direct transition towards dimerized phases, which can be understood as staggered valence bond solid orders, in which part of the lattice symmetry is spontaneously broken and the Majorana fermions acquire a mass gap. For small to intermediate $N$, on the other hand, phases with spontaneously broken SO($N$) symmetry, which can be understood as generalized N\'eel antiferromagnets, may be stabilized. These antiferromagnetic phases feature fully gapped fermion spectra for even $N$, but gapless Majorana modes for odd $N$. While the transitions between Majorana semimetal and dimerized phases are strongly first order, the transitions between Majorana semimetal and antiferromagnetic phases are continuous for small $N \leq 3$ and weakly first order for intermediate $N \geq 4$. The weakly-first-order nature of the latter transitions arises from fixed-point annihilation in the corresponding effective field theory, which contains a real symmetric tensorial order parameter coupled to the gapless Majorana degrees of freedom, realizing interesting examples of fluctuation-induced first-order transitions.

Nodeless superconductivity in cuprates with Ba$$_2$$CuO$$_3$$-type structure

In this work zero-temperature phase diagrams of cuprates with Ba$_2$CuO$_3$-type CuO chain structure is investigated. The projective symmetry group analysis is employed in the strong coupling limit, and renormalization group with bosonization analysis is employed in the weak coupling limit. We find that in both of these two limits, large areas of the phase diagrams are filled with nodeless $s$-wave superconducting phases (with weak $d$-wave components), instead of pure $d$-wave phase mostly found in cuprates. This implies that nodeless $s$-wave phase is the dominant superconducting phase in cuprates with Ba$_2$CuO$_3$-type CuO chain structure in low temperature. Other phases are also found, including $(s+d)$-wave superconducting phases and Luttinger liquid phases.

Optimal dense coding and quantum phase transition in Ising-XXZ diamond chain

We theoretically propose a dense coding scheme based on an infinite spin-1/2 Ising-XXZ diamond chain, where the Heisenberg spins dimer can be considered as a quantum channel. Using the transfer-matrix approach, we can obtain the analytical expression of the optimal dense coding capacity χ. The effects of anisotropy, external magnetic field, and temperature on χ in the diamond-like chain are discussed, respectively. It is found that χ is decayed with increasing the temperature, while the valid dense coding (χ>1) can be carried out by tuning the anisotropy parameter. Additionally, a certain external magnetic field can stimulate the enhancement of dense coding capacity. In an infinite spin-1/2 Ising-XXZ diamond chain, it has been shown that there exist two kinds of quantum phase transitions (QPTs) in zero-temperature phase diagram, i.e., one is that the ground state of system changes from the unentangled ferrimagnetic state to the entangled frustrated state and the other is from the entangled frustrated state to the unentangled ferromagnetic state. Here, we propose that optimal dense coding capacity χ can be regarded as a new detector of QPTs in this diamond-like chain, and the relationship between χ and QPTs is well established.