Infinite Density Matrix Renormalization Group Research Articles

Non-Abelian anyons with Rydberg atoms

We study the emergence of topological matter in two-dimensional systems of neutral Rydberg atoms in Ruby lattices. While Abelian anyons have been predicted in such systems, non-Abelian anyons, which would form a substrate for fault-tolerant quantum computing, have not been generated. To generate anyons with non-Abelian braiding statistics, we consider systems with mixed-boundary punctures. We obtain the topologically distinct ground states of the system numerically using the infinite Density Matrix Renormalization Group technique. We discuss how these topological states can be created using ancilla atoms of a different type. We show that a system with $2N+2$ punctures and an equal number of ancilla atoms leads to $N$ logical qubits whose Hilbert space is determined by a set of stabilizing conditions on the ancilla atoms. Quantum gates can be implemented using a set of gates acting on the ancilla atoms that commute with the stabilizers and realize the braiding group of non-Abelian Ising anyons.

Plaquette valence bond state in the spin- 12 J1−J2 XY model on the square lattice

We studied the ground-state phase diagram of spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2} XY$ model on the square lattice with first-neighbor ${J}_{1}$ and second-neighbor ${J}_{2}$ antiferromagnetic interactions using both infinite density matrix renormalization group (iDMRG) and DMRG approaches. We show that a plaquette valence bond phase is realized in an intermediate region $0.50\ensuremath{\le}{J}_{2}/{J}_{1}<0.54$ between a N\'eel magnetic ordered phase at ${J}_{2}/{J}_{1}<0.50$ and a stripy magnetic ordered phase at ${J}_{2}/{J}_{1}\ensuremath{\ge}0.54$. The plaquette valence bond phase is characterized by finite dimer orders in both the horizontal and vertical directions. Contrary to the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model, we do not find numerical evidence for a quantum spin liquid phase in the ${J}_{1}\text{\ensuremath{-}}{J}_{2} XY$ model.

Infinite density matrix renormalization group method for solving the Ising model

Infinite density matrix renormalization group method for solving the Ising model

Robust orbital diamagnetism in correlated Dirac fermions

We study orbital diamagnetism at zero temperature in (2 + 1)-dimensional Dirac fermions with a short-range interaction which exhibits a quantum phase transition to a charge density wave (CDW) phase. We introduce orbital magnetic fields into spinless Dirac fermions on the π-flux square lattice, and analyze them by using infinite density matrix renormalization group. It is found that the diamagnetism remains intact in the Dirac semimetal regime, while it is monotonically suppressed in the CDW regime. Around the quantum critical point of the CDW phase transition, we find a scaling behavior of the diamagnetism characteristic of the chiral Ising universality class. Besides, the scaling analysis implies that the robust orbital diamagnetism at weak magnetic fields in a Dirac semimetal regime would hold not only in our model but also in other interacting Dirac fermion systems as long as scaling regions are wide enough. The scaling behavior may also be regarded as a quantum, magnetic analogue of the critical Casimir effect which has been widely studied for classical phase transitions.

Abelian SU(N)1 chiral spin liquids on the square lattice

In the physics of the fractional quantum Hall (FQH) effect, a zoo of Abelian topological phases can be obtained by varying the magnetic field. Aiming to reach the same phenomenology in spin like systems, we propose a family of $\mathrm{SU}(N)$-symmetric models in the fundamental representation, on the square lattice with short-range interactions restricted to triangular units, a natural generalization for arbitrary $N$ of an SU(3) model studied previously where time-reversal symmetry is broken explicitly. Guided by the recent discovery of ${\mathrm{SU}(2)}_{1}$ and ${\mathrm{SU}(3)}_{1}$ chiral spin liquids (CSL) on similar models we search for topological $\mathrm{SU}({N)}_{1}$ CSL in some range of the Hamiltonian parameters via a combination of complementary numerical methods such as exact diagonalizations (ED), infinite density matrix renormalization group (iDMRG) and infinite Projected Entangled Pair State (iPEPS). Extensive ED on small (periodic and open) clusters up to $N=10$ and an innovative $\mathrm{SU}(N)$-symmetric version of iDMRG to compute entanglement spectra on (infinitely long) cylinders in all topological sectors provide unambiguous signatures of the $\mathrm{SU}({N)}_{1}$ character of the chiral liquids. An SU(4)-symmetric chiral PEPS, constructed in a manner similar to its SU(2) and SU(3) analogs, is shown to give a good variational ansatz of the $N=4$ ground state, with chiral edge modes originating from the PEPS holographic bulk-edge correspondence. Finally, we discuss the possible observation of such Abelian CSL in ultracold atom setups where the possibility of varying $N$ provides a tuning parameter similar to the magnetic field in the physics of the FQH effect.

Towards a topological quantum chemistry description of correlated systems: The case of the Hubbard diamond chain

The recently introduced topological quantum chemistry (TQC) framework has provided a description of universal topological properties of all possible band insulators in all space groups based on crystalline unitary symmetries and time reversal. While this formalism filled the gap between the mathematical classification and the practical diagnosis of topological materials, an obvious limitation is that it only applies to weakly interacting systems-which can be described within band theory. It is an open question to which extent this formalism can be generalized to correlated systems that can exhibit symmetry protected topological phases which are not adiabatically connected to any band insulator. In this work we address the many facettes of this question by considering the specific example of a Hubbard diamond chain. This model features a Mott insulator, a trivial insulating phase and an obstructed atomic limit phase. Here we discuss the nature of the Mott insulator and determine the phase diagram and topology of the interacting model with infinite density matrix renormalization group calculations, variational Monte Carlo simulations and with many-body topological invariants. We then proceed by considering a generalization of the TQC formalism to Green's functions combined with the concept of topological Hamiltonian to identify the topological nature of the phases, using cluster perturbation theory to calculate the Green's functions. The results are benchmarked with the above determined phase diagram and we discuss the applicability and limitations of the approach and its possible extensions.

Bosonic fractional Chern insulating state at integer fillings in a multiband system

The integer quantum Hall state occurs when the Landau levels are fully occupied by the fermions, while the fractional quantum Hall state usually emerges when the Landau level is partially filled by the strongly correlated fermions or bosons. Here, we report two fractional Chern insulating states of the hard-core bosons in a multi-band lattice model hosting topological flat bands with high Chern number. The previously proposed $\nu=1/3$ fractional Chern insulating state inherited from the high Chern number $C=2$ of the lowest topological flat band is revisited by the infinite density matrix renormalization group algorithm. In particular, we numerically identify a bosonic $1/2$-Laughlin-like fractional Chern insulating state at the integer fillings. We show two lower topological flat bands jointly generate an effective $C=1$ Chern band with half-filling. Furthermore, we find a strictly particle-hole-like symmetry between the $\nu$ and $3-\nu$ filling in our model. These findings extend our understanding of quantum Hall states and offer a new route to realize the novel fractional states in the system with multi-bands and high-Chern numbers.

Stability, phase transitions, and numerical breakdown of fractional Chern insulators in higher Chern bands of the Hofstadter model

The Hofstadter model is a popular choice for theorists investigating the fractional quantum Hall effect on lattices, due to its simplicity, infinite selection of topological flat bands, and increasing applicability to real materials. In particular, fractional Chern insulators in bands with Chern number $|C|>1$ can demonstrate richer physical properties than continuum Landau level states and have recently been detected in experiments. Motivated by this, we examine the stability of fractional Chern insulators with higher Chern number in the Hofstadter model, using large-scale infinite density matrix renormalization group (iDMRG) simulations on a thin cylinder. We confirm the existence of fractional states in bands with Chern numbers $C=1,2,3,4,5$ at the filling fractions predicted by the generalized Jain series [Phys. Rev. Lett. 115, 126401 (2015)]. Moreover, we discuss their metal-to-insulator phase transitions, as well as the subtleties in distinguishing between physical and numerical stability. Finally, we comment on the relative suitability of fractional Chern insulators in higher Chern number bands for proposed modern applications.

Four-Spin Terms and the Origin of the Chiral Spin Liquid in Mott Insulators on the Triangular Lattice.

At strong repulsion, the triangular-lattice Hubbard model is described by s=1/2 spins with nearest-neighbor antiferromagnetic Heisenberg interactions and exhibits conventional 120° order. Using the infinite density matrix renormalization group and exact diagonalization, we study the effect of the additional four-spin interactions naturally generated from the underlying Mott-insulator physics of electrons as the repulsion decreases. Although these interactions have historically been connected with a gapless ground state with emergent spinon Fermi surface, we find that, at physically relevant parameters, they stabilize a chiral spin liquid (CSL) of Kalmeyer-Laughlin (KL) type, clarifying observations in recent studies of the Hubbard model. We then present a self-consistent solution based on a mean-field rewriting of the interaction to obtain a Hamiltonian with similarities to the parent Hamiltonian of the KL state, providing a physical understanding for the origin of the CSL.

Abelian topological order of ν=2/5 and 3/7 fractional quantum Hall states in lattice models

Determining the statistics of elementary excitations supported by fractional quantum Hall states is crucial to understanding their properties and potential applications. In this paper, we use the topological entanglement entropy as an indicator of Abelian statistics to investigate the single-component $\nu=2/5$ and $3/7$ states for the Hofstadter model in the band mixing regime. We perform many-body simulations using the infinite cylinder density matrix renormalization group and present an efficient algorithm to construct the area law of entanglement, which accounts for both numerical and statistical errors. Using this algorithm, we show that the $\nu=2/5$ and $3/7$ states exhibit Abelian topological order in the case of two-body nearest-neighbor interactions. Moreover, we discuss the sensitivity of the proposed method and fractional quantum Hall states with respect to interaction range and strength.

Phase transitions in the one-dimensional ionic Hubbard model

We study quantum phase transitions by measuring the bond energy, the number density, and the half-chain entanglement entropy in the one-dimensional ionic Hubbard model. By using the matrix product operator to perform the infinite density matrix renormalization group, we obtain the ground states in the canonical form of matrix product states. Depending on the chemical potential and the staggered potential, the number density and the half-chain entanglement entropy shows clear signatures of a Mott transition. Our results confirm the success of using the matrix product operator method to investigate itinerant fermion systems.

Interplay of local order and topology in the extended Haldane-Hubbard model

We investigate the ground-state phase diagram of the spinful extended Haldane-Hubbard model on the honeycomb lattice using an exact-diagonalization, mean-field variational approach, and further complement it with the infinite density matrix renormalization group, applied to an infinite honeycomb cylinder. This model, governed by both on-site and nearest-neighbor interactions, can result in two types of insulators with finite local order parameters, either with spin or charge ordering. Moreover, a third one, a topologically nontrivial insulator with nonlocal order, is also manifest. We test expectations of previous analyses in spinless versions asserting that once a local order parameter is formed, the topological characteristics of the ground state, associated with a finite Chern number, are no longer present, resulting in a topologically trivial wave function. Our study confirms this overall picture, and highlights how finite-size effects may result in misleading conclusions on the coexistence of finite local order parameters and nontrivial topology in this model.

Many-Body Invariants for Chern and Chiral Hinge Insulators.

We construct new many-body invariants for 2D Chern and 3D chiral hinge insulators characterizing quantized pumping of bulk dipole and quadrupole moments. The many-body invariants are written entirely in terms of many-body ground state wave functions on a torus geometry with twisted boundary conditions and a set of unitary operators. We present a number of supporting arguments for the invariants via topological field theory interpretation, adiabatic pumping argument, and direct mapping to free-fermion band indices. Therefore, the invariants explicitly encircle several different pillars of theoretical descriptions of topological phases. Furthermore, our many-body invariants are written in forms which can be directly employed in various numerics including the exact diagonalization and the density-matrix renormalization group simulations. We finally confirm our invariants by numerical computations including an infinite density-matrix renormalization group on quasi-one-dimensional systems.

Ground-state properties of the one-dimensional Hubbard model with pairing potential

We consider a modification of the one-dimensional Hubbard model by including an external pairing potential. Guided by analytic bosonization results, we quantitatively determine the grand-canonical zero-temperature phase diagram using both finite and infinite density matrix renormalization group algorithm based on the formalism of matrix product states and matrix product operator, respectively. By computing various local quantities as well as the half-system entanglement, we are able to distinguish between Mott, metallic and superconducting phases. We point out the compressible nature of the Mott phase and the fully gapped nature of the many-body spectrum of the superconducting phase, in the presence of explicit U(1)-charge symmetry breaking.

Emergence of nematic paramagnet via quantum order-by-disorder and pseudo-Goldstone modes in Kitaev magnets

The appearance of nontrivial phases in Kitaev materials exposed to an external magnetic field has recently been a subject of intensive studies. Here, we elucidate the relation between the field-induced ground states of the classical and quantum spin models proposed for such materials, by using the infinite density matrix renormalization group (iDMRG) and the linear spin wave theory (LSWT). We consider the $K \Gamma \Gamma'$ model, where $\Gamma$ and $\Gamma'$ are off-diagonal spin exchanges on top of the dominant Kitaev interaction $K$. Focusing on the magnetic field along the $[111]$ direction, we explain the origin of the nematic paramagnet, which breaks the lattice-rotational symmetry and exists in an extended window of magnetic field, in the quantum model. This phenomenon can be understood as the effect of quantum order-by-disorder in the frustrated ferromagnet with a continuous manifold of degenerate ground states discovered in the corresponding classical model. We compute the dynamical spin structure factors using a matrix operator based time evolution and compare them with the predictions from LSWT. We, thus, provide predictions for future inelastic neutron scattering experiments on Kitaev materials in an external magnetic field along the $[111]$ direction. In particular, the nematic paramagnet exhibits a characteristic pseudo-Goldstone mode which results from the lifting of a continuous degeneracy via quantum fluctuations.

Fermi Surfaces of Composite Fermions

The fractional quantum Hall (FQH) effect was discovered in two-dimensional electron systems subject to a large perpendicular magnetic field nearly four decades ago. It helped launch the field of topological phases, and in addition, because of the quenching of the kinetic energy, gave new meaning to the phrase "correlated matter". Most FQH phases are gapped like insulators and superconductors; however, a small subset with even denominator fractional fillings nu of the Landau level, typified by nu = 1/2, are found to be gapless, with a Fermi surface akin to metals. We discuss our results, obtained numerically using the infinite Density Matrix Renormalization Group (iDMRG) scheme, on the effect of non-isotropic distortions with discrete N-fold rotational symmetry of the Fermi surface at zero magnetic field on the Fermi surface of the correlated nu = 1/2 state. We find that while the response for N = 2 (elliptical) distortions is significant (and in agreement with experimental observations with no adjustable parameters), it decreases very rapidly as N is increased. Other anomalies, like resilience to breaking the Fermi surface into disjoint pieces, are also found. This highlights the difference between Fermi surfaces formed from the kinetic energy, and those formed of purely potential energy terms in the Hamiltonian.

Quantum criticality of magnetic catalysis in two-dimensional correlated Dirac fermions

We study quantum criticality of the magnetic field induced charge density wave (CDW) order in correlated spinless Dirac fermions on the $\pi$-flux square lattice at zero temperature as a prototypical example of the magnetic catalysis, by using the infinite density matrix renormalization group. It is found that the CDW order parameter $M(B)$ exhibits an anomalous magnetic field $(B)$ scaling behavior characteristic of the $(2+1)$-dimensional chiral Ising universality class near the quantum critical point, which leads to a strong enhancement of $M(B)$ compared with a mean field result. We also establish a global phase diagram in the interaction-magnetic field plane for the fermionic quantum criticality.

Engineering quantum Hall phases in a synthetic bilayer graphene system

Synthetic quantum Hall bilayer (SQHB), realized by optically driven monolayer graphene in the quantum Hall regime, provides a flexible platform for engineering quantum Hall phases as discussed in [Phys. Rev. Lett. 119, 247403]. The coherent driving which couples two Landau levels mimicks an effective tunneling between synthetic layers. The tunneling strength, the effective Zeeman coupling, and two-body interaction matrix elements are tunable by varying the driving frequency and the driving strength. Using infinite density matrix renormalization group (iDMRG) techniques combined with exact diagonalization (ED), we show that the system exhibits a non-abelian bilayer Fibonacci phase at filling fraction $\nu = 2/3$. Moreover, at integer filling $\nu = 1$, the SQHB exhibits quantum Hall ferromagnetism. Using Hartree-Fock theory and exact diagonalization, we show that excitations of the quantum Hall ferromagnet are topological textures known as skyrmions.

Signatures of a Deconfined Phase Transition on the Shastry-Sutherland Lattice: Applications to Quantum Critical SrCu2(BO3)2

We study a possible deconfined quantum phase transition in a realistic model of a two-dimensional Shastry-Sutherland quantum magnet, using both numerical and field theoretic techniques. Using the infinite density matrix renormalization group (iDMRG) method, we verify the existence of an intermediate plaquette valence bond solid (pVBS) order, with two fold degeneracy, between the dimer and N\'eel ordered phases. We argue that the quantum phase transition between the N\'eel and pVBS orders may be described by a deconfined quantum critical point (DQCP) with an emergent O(4) symmetry. By analyzing the correlation length spectrum obtained from iDMRG, we provide evidence for the DQCP and emergent O(4) symmetry in the lattice model. Such a phase transition has been reported in the recent pressure tuned experiments in the Shastry-Sutherland lattice material $\mathrm{SrCu_2 (BO_3)_2}$. The non-symmorphic lattice structure of the Shastry-Sutherland compound leads to extinction points in the scattering, where we predict sharp signatures of a DQCP in both the phonon and magnon spectra associated with the spinon continuum. The effect of weak interlayer couplings present in the three dimensional material is also discussed. Our results should help guide the experimental study of DQCP in quantum magnets.

Quantum phase transition and criticality in quasi-one-dimensional spinless Dirac fermions

We study quantum criticality of spinless fermions on the quasi one dimensional $\pi$-flux square lattice in cylinder geometry, by using the infinite density matrix renormalization group and abelian bosonization. For a series of the cylinder circumferences $L_y=4n+2=2, 6, \cdots$ with the periodic boundary condition, there are quantum phase transitions from gapped Dirac fermion states to charge density wave (CDW) states. We find that the quantum phase transitions for such circumferences are continuous and belong to the (1+1)-dimensional Ising universality class. On the other hand, when $L_y=4n=4, 8, \cdots$, there are gapless Dirac fermions at the non-interacting point and the phase transition to the CDW state is Gaussian. Both of these two criticalities are described in a unified way by the bosonization. We clarify their intimate relationship and demonstrate that a central charge $c=1/2$ Ising transition line arises as a critical state of an emergent Majorana fermion from the $c=2$ Gaussian transition point.