Density Matrix Renormalization Group Technique Research Articles

Entanglement properties of the Haldane phases: A finite system-size approach

We study the bond-alternating Heisenberg model using the finite-size density-matrix renormalization group (DMRG) technique and analytical arguments based on the matrix product state, where we pay particular attention to the boundary-condition dependence on the entanglement spectrum of the system. We show that, in the antiperiodic boundary condition (APBC), the parity quantum numbers are equivalent to the topological invariants characterizing the topological phases protected by the bond-centered inversion and $\pi$ rotation about $z$ axis. We also show that the odd parity in the APBC, which characterizes topologically nontrivial phases, can be extracted as a two-fold degeneracy in the entanglement spectrum even with finite system size. We then determine the phase diagram of the model with the uniaxial single-ion anisotropy using the level spectroscopy method in the DMRG technique. These results not only suggest the detectability of the symmetry protected topological (SPT) phases via general twisted boundary conditions but also provide a useful and precise numerical tool for discussing the SPT phases in the exact diagonalization and DMRG techniques.

Kitaev exchange and field-induced quantum spin-liquid states in honeycomb α-RuCl3.

Large anisotropic exchange in 5d and 4d oxides and halides open the door to new types of magnetic ground states and excitations, inconceivable a decade ago. A prominent case is the Kitaev spin liquid, host of remarkable properties such as protection of quantum information and the emergence of Majorana fermions. Here we discuss the promise for spin-liquid behavior in the 4d5 honeycomb halide α-RuCl3. From advanced electronic-structure calculations, we find that the Kitaev interaction is ferromagnetic, as in 5d5 iridium honeycomb oxides, and indeed defines the largest superexchange energy scale. A ferromagnetic Kitaev coupling is also supported by a detailed analysis of the field-dependent magnetization. Using exact diagonalization and density-matrix renormalization group techniques for extended Kitaev-Heisenberg spin Hamiltonians, we find indications for a transition from zigzag order to a gapped spin liquid when applying magnetic field. Our results offer a unified picture on recent magnetic and spectroscopic measurements on this material and open new perspectives on the prospect of realizing quantum spin liquids in d5 halides and oxides in general.

Synthetic-gauge-field-induced resonances and Fulde-Ferrell-Larkin-Ovchinnikov states in a one-dimensional optical lattice

Coherent coupling generated by laser light between the hyperfine states of atoms, loaded in a 1D optical lattice, gives rise to the "synthetic dimension" system which is equivalent to a Hofstadter model in a finite strip of square lattice. An SU(M) symmetric attractive interaction in conjunction with the synthetic gauge field present in this system gives rise to unusual effects. We study the two- body problem of the system using the T-matrix formalism. We show that the two-body ground states pick up a finite momentum and can transform into two-body resonance like features in the scattering continuum with a large change in the phase shift. As a result, even for this 1D system, a critical amount of attraction is needed to form bound states. These phenomena have spectacular effects on the many body physics of the system analyzed using the numerical density matrix renormalization group technique. We show that the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states form in the system even for a "balanced" gas and the FFLO momentum of the pairs scales linearly with flux. Considering suitable measures, we investigate interesting properties of these states. We also discuss a possibility of realization of a generalized interesting topological model, called the Creutz ladder.

Diffusive and Subdiffusive Spin Transport in the Ergodic Phase of a Many-Body Localizable System.

We study high temperature spin transport in a disordered Heisenberg chain in the ergodic regime. By employing a density matrix renormalization group technique for the study of the stationary states of the boundary-driven Lindblad equation we are able to study extremely large systems (400 spins). We find both a diffusive and a subdiffusive phase depending on the strength of the disorder and on the anisotropy parameter of the Heisenberg chain. Studying finite-size effects, we show numerically and theoretically that a very large crossover length exists that controls the passage of a clean-system dominated dynamics to one observed in the thermodynamic limit. Such a large length scale, being larger than the sizes studied before, explains previous conflicting results. We also predict spatial profiles of magnetization in steady states of generic nondiffusive systems.

Interaction-driven fractional quantum Hall state of hard-core bosons on kagome lattice at one-third filling

There has been a growing interest in realizing topologically nontrivial states of matter in band insulators, where a quantum Hall effect can appear as an intrinsic property of the band structure. While the on-going progress is under way with a number of directions, the possibility of realizing novel interaction-generated topological phases, without the requirement of a nontrivial invariant encoded in single-particle wavefunction or band structure, can significantly extend the class of topological materials and is thus of great importance. Here, we show an interaction-driven topological phase emerging in an extended Bose-Hubbard model on kagome lattice, where the non-interacting band structure is topological trivial with zero Berry curvature in the Brillouin zone. By means of an unbiased state-of-the-art density-matrix renormalization group technique, we identify that the groundstate in a broad parameter region is equivalent to a bosonic fractional quantum Hall Laughlin state, based on the characterization of universal properties including groundstate degeneracy, edge excitations and anyonic quasiparticle statistics. Our work paves a way of finding interaction induced topological phase at the phase boundary of conventionally ordered solid phases.

Non-Markovian dynamics in chiral quantum networks with spins and photons

We study the dynamics of chiral quantum networks consisting of nodes coupled by unidirectional or asymmetric bidirectional quantum channels. In contrast to familiar photonic networks where driven two-level atoms exchange photons via 1D photonic nanostructures, we propose and study a setup where interactions between the atoms are mediated by spin excitations (magnons) in 1D $XX$ spin chains representing spin waveguides. While Markovian quantum network theory eliminates quantum channels as structureless reservoirs in a Born-Markov approximation to obtain a master equation for the nodes, we are interested in non-Markovian dynamics. This arises from the nonlinear character of the dispersion with band-edge effects, and from finite spin propagation velocities leading to time delays in interactions. To account for the non-Markovian dynamics we treat the quantum degrees of freedom of the nodes and connecting channels as a composite spin system with the surrounding of the quantum network as a Markovian bath, allowing for an efficient solution with time-dependent density matrix renormalization group techniques. We illustrate our approach showing non-Markovian effects in the driven-dissipative formation of quantum dimers, and we present examples for quantum information protocols involving quantum state transfer with engineered elements as basic building blocks of quantum spintronic circuits.

Magnetic phase transition in coherently coupled Bose gases in optical lattices

We describe the ground state of a gas of bosonic atoms with two coherently coupled internal levels in a deep optical lattice in a one dimensional geometry. In the single-band approximation this system is described by a Bose-Hubbard Hamiltonian. The system has a superfluid and a Mott insulating phase which can be either paramagnetic or ferromagnetic. We characterize the quantum phase transitions at unit filling by means of a density matrix renormalization group technique and compare it with a mean-field approach. The presence of the ferromagnetic Ising-like transition modifies the Mott lobes. In the Mott insulating region the system maps to the ferromagnetic spin-1/2 XXZ model in a transverse field and the numerical results compare very well with the analytical results obtained from the spin model. In the superfluid regime quantum fluctuations strongly modify the phase transition with respect to the well established mean-field three dimensional classical bifurcation.

Symmetry-conserving purification of quantum states within the density matrix renormalization group

The density matrix renormalization group (DMRG) algorithm was originally designed to efficiently compute the zero temperature or ground-state properties of one dimensional strongly correlated quantum systems. The development of the algorithm at finite temperature has been a topic of much interest, because of the usefulness of thermodynamics quantities in understanding the physics of condensed matter systems, and because of the increased complexity associated with efficiently computing temperature-dependent properties. The ancilla method is a DMRG technique that enables the computation of these thermodynamic quantities. In this paper, we review the ancilla method, and improve its performance by working on reduced Hilbert spaces and using canonical approaches. We furthermore explore its applicability beyond spins systems to t-J and Hubbard models.

Entanglement entropy and the Berezinskii–Kosterlitz–Thouless phase transition in the J1–J2 Heisenberg chain

Using the density matrix renormalization group (DMRG) technique, we study the Berezinskii–Kosterlitz–Thouless (BKT) quantum phase transition (QPT) in the J1–J2 Heisenberg chain model from the quantum entanglement point of view. It is found that the gap behavior between two neighboring two-site entanglement entropies as well as the first derivative of both the two-site entropy and the block entropy can be used as indicators for the BKT phase transition in this model. The corresponding size dependent scaling behaviors are analyzed, respectively. Our numerical results give direct evidence for the effectiveness of the entanglement in the BKT-type QPT indicating from different aspects.

Magnetic Excitations of Spin Nematic State in Frustrated Ferromagnetic Chain

By exploiting density-matrix renormalization group techniques, we investigate the dynamical spin structure factor of a spin-1/2 Heisenberg chain with ferromagnetic nearest-neighbor and antiferromagnetic next-nearest-neighbor exchange interactions in an applied magnetic field. In a field-induced spin nematic regime, we find gapless longitudinal and gapped transverse spin excitation spectra, in accordance with quasi-long-ranged longitudinal and short-ranged transverse spin correlations, respectively. The gapless point coincides with the dominant longitudinal spin correlation, whereas the gap position exhibits a characteristic field dependence contradicting the dominant transverse spin correlation.

Computational investigation on tunable optical band gap in armchair polyacenes.

Polyacenes in their armchair geometry (phenacenes) have recently been found to possess appealing electronic and optical properties with higher chemical stability and comparatively larger band gap as compared to linear polyacenes. They also behave as high-temperature superconductors upon alkali metal doping. Moreover, the optical properties of crystalline picene can be finely tuned by applying external pressure. We investigated the variation of optical gap as a function of altering the interplanar distances between parallel cofacial phenacene dimers. We employed both time-dependent density functional theory and density matrix renormalization group (DMRG) technique to investigate the lowest singlet excitations in phenacene dimer. Our study showed that the lowest singlet excitation in these systems evolved as a function of interplanar separation. The optical excitation energy gap decreases as a function of inverse interplanar separation of the phenacene dimer. The distant dependent variation of optical absorption at the dimer level may be comparable with experimental observation in picene crystal under pressure. DMRG study also demonstrates that besides picene, electronic properties of higher phenacenes can also be tunable by altering interplanar separation.

Phase diagram of theZ3parafermionic chain with chiral interactions

Parafermions are exotic quasiparticles with non-Abelian fractional statistics that can be realized and stabilized in 1-dimensional models that are generalizations of the Kitaev p-wave wire. We study the simplest generalization, i.e. the $Z_3$ parafermionic chain. Using a Jordan-Wigner transform we focus on the equivalent three-state chiral clock model, and study its rich phase diagram using the density matrix renormalization group technique. We perform our analyses using quantum entanglement diagnostics which allow us to determine phase boundaries, and the nature of the phase transitions. In particular, we study the transition between the topological and trivial phases, as well as to an intervening incommensurate phase which appears in a wide region of the phase diagram. The phase diagram is predicted to contain a Lifshitz type transition which we confirm using entanglement measures. We also attempt to locate and characterize a putative tricritical point in the phase diagram where the three above mentioned phases meet at a single point.

Competing Abelian and non-Abelian topological orders inν=1/3+1/3quantum Hall bilayers

Bilayer quantum Hall systems, realized either in two separated wells or in the lowest two sub-bands of a wide quantum well, provide an experimentally realizable way to tune between competing quantum orders at the same filling fraction. Using newly developed density matrix renormalization group techniques combined with exact diagonalization, we return to the problem of quantum Hall bilayers at filling $\nu = 1/3 + 1/3$. We first consider the Coulomb interaction at bilayer separation $d$, bilayer tunneling energy $\Delta_\textrm{SAS}$, and individual layer width $w$, where we find a phase diagram which includes three competing Abelian phases: a bilayer-Laughlin phase (two nearly decoupled $\nu = 1/3$ layers); a bilayer-spin singlet phase; and a bilayer-symmetric phase. We also study the order of the transitions between these phases. A variety of non-Abelian phases have also been proposed for these systems. While absent in the simplest phase diagram, by slightly modifying the interlayer repulsion we find a robust non-Abelian phase which we identify as the "interlayer-Pfaffian" phase. In addition to non-Abelian statistics similar to the Moore-Read state, it exhibits a novel form of bilayer-spin charge separation. Our results suggest that $\nu = 1/3 + 1/3$ systems merit further experimental study.

Topological superconductivity in the one-dimensional interacting Creutz model

We consider one-dimensional topological insulators characterized by zero energy end states. In presence of proximity induced pairing, those end states can become Majorana states. We study here the fate of those various end states when Hubbard electron-electron repulsive interactions are added, using a combination of mean-field theory and density matrix renormalization group techniques.

Entanglement properties of the nontrivial Haldane insulator in the 1D extended Bose-Hubbard model

We investigate the entanglement properties of a nontrivial topological phase in the one-dimensional (1D) Bose-Hubbard model with additional nearest-neighbor repulsion. Employing the large-scale density-matrix renormalization group technique we show that a gapped insulating phase protected by lattice inversion symmetry, the so-called Haldane insulator, appears between the Mott and density wave phases in the intermediate-coupling regime. The phase boundaries were determined from the central charge via the von Neumann entropy. The Haldane insulator reveals a characteristic degeneracy in the entanglement spectra. Breaking the lattice inversion symmetry strongly affects the distinctive gapped dispersion of the dynamical charge response of the bosonic Haldane insulator.

Effects of magnetic anisotropy on spin dynamics of ferromagnetic frustrated chain

By exploiting density-matrix renormalization group techniques, we investigate the spin dynamics of a spin-1/2 one-dimensional J1-J2 XXZ model with competing ferromagnetic J1 and antiferromagnetic J2 exchange couplings under applied magnetic fields. Numerical results of spin excitation spectra show that in the field-induced spin quadrupole regime, the longitudinal component has a gapless mode and the transverse component has a gapped mode irrespective of the exchange anisotropy. The excitation gap of the transverse spin excitation increases as the exchange anisotropy increases over the XY-like and Ising-like regions, demonstrating that two-magnon bound states are stabilized due to the easy-axis anisotropy.

Comparative density-matrix renormalization group study of symmetry-protected topological phases in spin-1 chain and Bose-Hubbard models

We reexamine the one-dimensional spin-1 $XXZ$ model with on-site uniaxial single-ion anisotropy as to the appearance and characterization of the symmetry-protected topological Haldane phase. By means of large-scale density-matrix renormalization group (DMRG) calculations the central charge can be determined numerically via the von Neumann entropy, from which the ground-sate phase diagram of the model can be derived with high precision. The nontrivial gapped Haldane phase shows up in between the trivial gapped even Haldane and N{\'e}el phases, appearing at large single-ion and spin--exchange interaction anisotropies, respectively. We furthermore carve out a characteristic degeneracy of the lowest entanglement level in the topological Haldane phase, which is determined using a conventional finite-system DMRG technique with both periodic and open boundary conditions. Defining the spin and neutral gaps in analogy to the single-particle and neutral gaps in the intimately connected extended Bose-Hubbard model, we show that the excitation gaps in the spin model qualitatively behave just as for the bosonic system. We finally compute the dynamical spin structure factor in the three different gapped phases and find significant differences in the intensity maximum which might be used to distinguish these phases experimentally.

Critical fluctuations of a confined binary mixture

Exploiting the mapping between a binary mixture and the Ising model we have analyzed the critical fluctuations by means of the density-matrix renormalization group technique. The calculations have been carried out for a two-dimensional Ising strip subject to equally strong surface fields. It was found that the critical Casimir force has significantly different behavior on opposite sides of the capillary condensation line, especially below the critical temperature. It can be concluded that in real binary mixtures t he most attractive force appears at temperatures near �� and at reservoir compositions slightly away from the critical composition.

Destruction of long-range order by quenching of the hopping range in one dimension

We study the dynamics in a one-dimensional hard-core Bose gas with power-law hopping after an abrupt reduction of the hopping range using the time-dependent density-matrix renormalization-group and bosonization techniques. In particular, we focus on the destruction of the Bose-Einstein condensate (BEC), which is present in the initial state in the thermodynamic limit. We argue that this type of quench is akin to a sudden reduction in the effective dimensionality $d$ of the system (from $d>1$ to $d=1)$. We identify two regimes in the evolution of the BEC fraction. For short times the decay of the BEC fraction is Gaussian while for intermediate to long times it is well described by a stretched exponential with an exponent that depends on the initial effective dimensionality of the system. These results are potentially relevant for cold trapped-ion experiments which can simulate an equivalent of hard-core bosons, i.e., spins, with tunable long-range interactions.

Magnetic and nonmagnetic phases in dopedAB2 t−JHubbard chains

We discuss the rich phase diagram of doped AB2 $t-J$ chains by using data from density matrix renormalization group and exact diagonalization techniques. The $J$ vs $\delta$ (hole doping) phase diagram exhibits regions of itinerant ferrimagnetism, incommensurate, resonating valence bond and Nagaoka states, phase separation, and Luttinger liquid (LL) physics. Several features are highlighted, such as the modulated ferrimagnetic structure, the occurrence of Nagaoka spin polarons in the underdoped regime and small values of $J=4t^2/U$, where $t$ is the first-neighbor hopping amplitude and $U$ is the on-site repulsive Coulomb interaction, incommensurate structures with nonzero magnetization, and strong-coupling LL physics in the high-doped regime. We also verify that relevant findings are in agreement with the corresponding findings in square and n-leg ladder lattices. In particular, we mention the instability of Nagaoka ferromagnetism against $J$ and $\delta$.