Numerical Density Matrix Renormalization Group Research Articles

Proximate Tomonaga-Luttinger liquid in an anisotropic Kitaev-Gamma model

The characterization of quantum spin liquid phases in Kitaev materials has been a subject of intensive studies over the recent years, both theoretically and experimentally. Most theoretical studies have focused on an isotropically interacting model with its coupling strength being equivalent on each bond in an attempt to simplify the problem. In this Letter, we study an extended spin-1/2 Kitaev-Γ model on a honeycomb lattice with an additional tuning parameter that controls the coupling strength on one of the bonds: We connect the limit of isolated Kitaev-Γ chains, which is known to exhibit an emergent SU(2)1 Tomonaga-Luttinger liquid phase [Yang , ], to the two-dimensional model. We report on an instance, in which properties of the Tomonaga-Luttinger liquid persist for finite interchain coupling forming a phase we call the proximate Tomonaga-Luttinger liquid (pTLL). This pTLL is strongly anisotropic in character analogous to , and features spinon-like excitations similar to those of the antiferromagnetic Heisenberg chain. We use numerical exact diagonalization and density matrix renormalization group on various cluster geometries in a complementary way to overcome finite-size limitations. Published by the American Physical Society 2024

Detecting hidden order in fractional Chern insulators

Topological phase transitions go beyond Ginzburg and Landau's paradigm of spontaneous symmetry breaking and occur without an associated local order parameter. Instead, such transitions can be characterized by the emergence of nonlocal order parameters, which require measurements on extensively many particles simultaneously—an impossible venture in real materials. On the other hand, quantum simulators have demonstrated such measurements, making them prime candidates for experimental confirmation of nonlocal topological order. Here, building upon the recent advances in preparing few-particle fractional Chern insulators using ultracold atoms and photons, we propose a realistic scheme for detecting the hidden off-diagonal long-range order (HODLRO) characterizing Laughlin states. Furthermore, we demonstrate the existence of this hidden order in fractional Chern insulators, specifically for the ν=1/2-Laughlin state in the isotropic Hofstadter-Bose-Hubbard model. This is achieved by large-scale numerical density matrix renormalization group (DMRG) simulations based on matrix product states, for which we formulate an efficient sampling procedure providing direct access to HODLRO in close analogy to the proposed experimental scheme. We confirm the characteristic power-law scaling of HODLRO, with an exponent 1/ν=2, and show that its detection requires only a few thousand snapshots. This makes our scheme realistically achievable with current technology and paves the way for further analysis of nonlocal topological orders, e.g., in topological states with non-Abelian anyonic excitations. Published by the American Physical Society 2024

CFT and lattice correlators near an RG domain wall between minimal models

Conformal interfaces separating two conformal field theories (CFTs) provide maps between different CFTs, and naturally exist in nature as domain walls between different phases. One particularly interesting construction of a conformal interface is the renormalization group (RG) domain wall between CFTs. For a given Virasoro minimal model \mathcal{M}_{k+3,k+2}ℳk+3,k+2, an RG domain wall can be generated by a specific deformation which triggers an RG flow towards its adjacent Virasoro minimal model \mathcal{M}_{k+2,k+1}ℳk+2,k+1 with the deformation turned on over part of the space. An algebraic construction of this domain wall was proposed by Gaiotto in [J. High Energy Phys. 12, 103 (2012)]. In this paper, we will provide a study of this RG domain wall for the minimal case k=2k=2, which can be thought of as a nonperturbative check of the construction. In this case the wall is separating the Tricritical Ising Model (TIM) CFT and the Ising Model (IM) CFT. We will check the analytical results of correlation functions from the RG brane construction with the numerical density matrix renormalization group (DMRG) calculation using a lattice model proposed in [arXiv.1206.1332, Science 344, 280 (2014)], and find a perfect agreement. We comment on possible experimental realizations of this RG domain wall.

Z3 and (×Z3)3 symmetry protected topological paramagnets

We identify two-dimensional three-state Potts paramagnets with gapless edge modes on a triangular lattice protected by (×Z3)3 ≡ Z3 × Z3 × Z3 symmetry and smaller Z3 symmetry. We derive microscopic models for the gapless edge, uncover their symmetries and analyze the conformal properties. We study the properties of the gapless edge by employing the numerical density-matrix renormalization group (DMRG) simulation and exact diagonalization. We discuss the corresponding conformal field theory, its central charge, and the scaling dimension of the corresponding primary field. We argue, that the low energy limit of our edge modes defined by the SUk(3)/SUk(2) coset conformal field theory with the level k = 2. The discussed two-dimensional models realize a variety of symmetry-protected topological phases, opening a window for studies of the unconventional quantum criticalities between them.

Effective Bi-Layer Model Hamiltonian and Density-Matrix Renormalization Group Study for the High-T c Superconductivity in La3Ni2O7 under High Pressure

High-T c superconductivity with possible T c ≈ 80 K has been reported in the single crystal of La3Ni2O7 under high pressure. Based on the electronic structure given by the density functional theory calculations, we propose an effective bi-layer model Hamiltonian including both 3d z 2 and 3d x 2–y 2 orbital electrons of the nickel cations. The main feature of the model is that the 3d z 2 electrons form inter-layer σ-bonding and anti-bonding bands via the apical oxygen anions between the two layers, while the 3d x 2–y 2 electrons hybridize with the 3d z 2 electrons within each NiO2 plane. The chemical potential difference of these two orbital electrons ensures that the 3d z 2 orbitals are close to half-filling and the 3d x 2–y 2 orbitals are near quarter-filling. The strong on-site Hubbard repulsion of the 3d z 2 orbital electrons gives rise to an effective inter-layer antiferromagnetic spin super-exchange J. Applying pressure can self dope holes on the 3d z 2 orbitals with the same amount of electrons doped on the 3d x 2–y 2 orbitals. By performing numerical density-matrix renormalization group calculations on a minimum setup and focusing on the limit of large J and small doping of 3d z 2 orbitals, we find the superconducting instability on both the 3d z 2 and 3d x 2–y 2 orbitals by calculating the equal-time spin singlet pair–pair correlation function. Our numerical results may provide useful insights in the high-T c superconductivity in single crystal La3Ni2O7 under high pressure.

Quantum correlations in the frustrated XY model on the honeycomb lattice

We have investigated the spin-1/2 XY frustrated antiferromagnetic Heisenberg honeycomb model, which features an intermediate region in its ground state phase diagram that is not well understood. The two dominant phases in the diagram are the quantum spin-liquid (QSL) and the antiferromagnetic Ising order. Quantum correlations suggest that the QSL phase is likely to exhibit entanglement. To explore this possibility, we utilized numerical Lanczos and density matrix renormalization group (DMRG) methods to calculate concurrence, quantum discord (QD), and entanglement entropy. The results of our study indicate the existence of quantum entanglement within the intermediate region, implying a greater probability for the dominance of the quantum spin-liquid (QSL) phase over the antiferromagnetic Ising order. This discovery underscores the importance of considering quantum correlations in comprehending the model’s behavior and provides insight into the complex nature of quantum systems.

Density Matrix Renormalization Group for Continuous Quantum Systems.

We introduce a versatile and practical framework for applying matrix product state techniques to continuous quantum systems. We divide space into multiple segments and generate continuous basis functions for the many-body state in each segment. By combining this mapping with existing numerical density matrix renormalization group routines, we show how one can accurately obtain the ground-state wave function, spatial correlations, and spatial entanglement entropy directly in the continuum. For a prototypical mesoscopic system of strongly interacting bosons we demonstrate faster convergence than standard grid-based discretization. We illustrate the power of our approach by studying a superfluid-insulator transition in an external potential. We outline how one can directly apply or generalize this technique to a wide variety of experimentally relevant problems across condensed matter physics and quantum field theory.

Fractional spin excitations and conductance in the spiral-staircase Heisenberg ladder

We investigate theoretically the spiral staircase Heisenberg spin-$1/2$ ladder in the presence of antiferromagnetic long-range spin interactions and a uniform magnetic field. As a special case we also consider the Kondo necklace model. If the magnetizations of the two chains forming the ladder satisfy a certain resonance condition, involving interchain couplings as perturbations, the system is in a partially gapped magnetic phase hosting excitations characterized by fractional spins, whose values can be changed by the magnetic field. We show that these fractional spin excitations can be probed by spin currents in a transport setup with a spin conductance that reveals the fractionalized spin. In some special cases, the spin conductance reaches universal values in units of $(g\mu_B)^2/h$, where $g$ is the $g$-factor, $\mu_B$ the Bohr magneton, and $h$ the Planck constant. We obtain our results with the help of bosonization and numerical density matrix renormalization group methods.

Two-Hole Ground State: Dichotomy in Pairing Symmetry

A single-hole ground state Ansatz for the two-dimensional t-J model has been recently studied by the variational Monte Carlo (VMC) method. Such a doped hole behaves like a "twisted" non-Landau quasiparticle characterized by an emergent quantum number in agreement with exact numerics. In this work, we further investigate the ground state of two holes by VMC. It is found that the two holes strongly attract each other to form a pairing state with a new quantum number the same as obtained by the numerical exact diagonalization and density matrix renormalization group (DMRG) calculations. A unique feature of this pairing state is a dichotomy in the pairing symmetry, i.e., a d-wave in terms of the electron c operators and an s-wave in terms of the new quasiparticles, as explicitly illustrated in the ground state wave function. A similar VMC study of a two-hole wave function for the t-J two-leg ladder also yields a good agreement with the DMRG result. We demonstrate that the pairing mechanism responsible for the strong binding here is not due to the long-range antiferromagnetic order nor the resonating-valence-bound pairing in the spin background but is the consequence of the quantum phase-strings created by the hopping of holes. The resulting spin-current pattern mediating the pairing force is explicitly illustrated in the VMC calculation. Physical implications to superconductivity at finite doping are also discussed.

Quantum phases of two-dimensional Z2 gauge theory coupled to single-component fermion matter

We investigate the rich quantum phase diagram of Wegner's theory of discrete Ising gauge fields interacting with $U(1)$ symmetric single-component fermion matter hopping on a two-dimensional square lattice. In particular limits the model reduces to (i) pure $\mathbb{Z}_2$ even and odd gauge theories, (ii) free fermions in a static background of deconfined $\mathbb{Z}_2$ gauge fields, (iii) the kinetic Rokhsar-Kivelson quantum dimer model at a generic dimer filling. We develop a local transformation that maps the lattice gauge theory onto a model of $\mathbb{Z}_2$ gauge-invariant spin $1/2$ degrees of freedom. Using the mapping, we perform numerical density matrix renormalization group calculations that corroborate our understanding of the limits identified above. Moreover, in the absence of the magnetic plaquette term, we reveal signatures of topologically ordered Dirac semimetal and staggered Mott insulator phases at half-filling. At strong coupling, the lattice gauge theory displays fracton phenomenology with isolated fermions being completely frozen and dimers exhibiting restricted mobility. In that limit, we predict that in the ground state dimers form compact clusters, whose hopping is suppressed exponentially in their size. We determine the band structure of the smallest clusters numerically using exact diagonalization. The rich phenomenology discussed in this paper can be probed in analog and digital quantum simulators of discrete gauge theories and in Kitaev spin-orbital liquids.

Prediction of Toric Code Topological Order from Rydberg Blockade

The physical realization of Z2 topological order as encountered in the paradigmatic toric code has proven to be an elusive goal. We predict that this phase of matter can be realized in a two-dimensional array of Rydberg atoms placed on the ruby lattice, at specific values of the Rydberg blockade radius. First, we show that the blockade model—also known as a “PXP” model—realizes a monomer-dimer model on the kagome lattice with a single-site kinetic term. This model can be interpreted as a Z2 gauge theory whose dynamics is generated by monomer fluctuations. We obtain its phase diagram using the numerical density matrix renormalization group method and find a topological quantum liquid (TQL) as evidenced by multiple measures including (i) a continuous transition between two featureless phases, (ii) a topological entanglement entropy of ln2 as measured in various geometries, (iii) degenerate topological ground states, and (iv) the expected modular matrix from ground state overlap. Next, we show that the TQL persists upon including realistic, algebraically decaying van der Waals interactions V(r)∼1/r6 for a choice of lattice parameters. Moreover, we can directly access topological loop operators, including the Fredenhagen-Marcu order parameter. We show how these can be measured experimentally using a dynamic protocol, providing a “smoking gun” experimental signature of the TQL phase. Finally, we show how to trap an emergent anyon and realize different topological boundary conditions, and we discuss the implications for exploring fault-tolerant quantum memories.12 MoreReceived 24 November 2020Accepted 19 May 2021DOI:https://doi.org/10.1103/PhysRevX.11.031005Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasTopological phases of matterTopological quantum computingPhysical SystemsRydberg atoms & moleculesCondensed Matter, Materials & Applied Physics

Modified strong-coupling treatment of a spin- 12 Heisenberg trimerized chain developed from the exactly solved Ising-Heisenberg diamond chain

Quantum spin-1/2 antiferromagnetic Heisenberg trimerized chain with strong intradimer and weak monomer-dimer coupling constants is studied using the novel many-body perturbation expansion, which is developed from the exactly solved spin-1/2 Ising-Heisenberg diamond chain preserving correlations between all interacting spins of the trimerized chain unlike the standard perturbation scheme developed on the grounds of noninteracting spin monomers and dimers. The Heisenberg trimerized chain shows the intermediate one-third plateau, which was also observed in the magnetization curve of the polymeric compound Cu$_3$(P$_2$O$_6$OH)$_2$ affording its experimental realization. Within the modified strong-coupling method we have obtained the effective Hamiltonians for the magnetic-field range from zero to one-third plateau, and from one-third plateau to the saturation magnetization. The second-order perturbation theory shows reliable results even for the moderate ratio between weaker dimer-monomer and stronger intradimer coupling constants. We have also examined thermodynamic properties and recovered the low-temperature peak in the specific heat. The accuracy of the developed method is tested through a comparison with numerical density-matrix renormalization group and quantum Monte Carlo simulations. Using the results for the effective Hamiltonian we suggest straightforward procedure for finding the microscopic parameters of one-dimensional trimerized magnetic compounds with strong intradimer and weak monomer-dimer couplings. We found the refined values for the coupling constants of Cu$_3$(P$_2$O$_6$OH)$_2$ by matching the theoretical results with the available experimental data for the magnetization and magnetic susceptibility in a wide range of temperatures and magnetic fields.

Current jumps in flat-band ladders with Dzyaloshinskii-Moriya interactions

Localized magnons states, due to flat bands in the spectrum, is an intensely studied phenomenon and can be found in many frustrated magnets of different spatial dimensionality. The presence of Dzyaloshinskii-Moriya (DM) interactions may change radically the behavior in such systems. In this context, we study a paradigmatic example of a one-dimensional frustrated antiferromagnet, the sawtooth chain in the presence of DM interactions. Using both path integrals methods and numerical Density Matrix Renormalization Group, we revisit the physics of localized magnons and determine the consequences of the DM interaction on the ground state. We have studied the spin current behavior, finding three different regimes. First, a Luttinger-liquid regime where the spin current shows a step behavior as a function of parameter $D$, at a low magnetic field. Increasing the magnetic field, the system is in the Meissner phase at the $m = 1/2$ plateau, where the spin current is proportional to the DM parameter. Finally, further increasing the magnetic field and for finite $D$ there is a small stiffness regime where the spin current shows, at fixed magnetization, a jump to large values at $D = 0$, a phenomenon also due to the flat band.

Interacting Rice-Mele model: Bulk and boundaries

We investigate the interacting, one-dimensional Rice-Mele model, a prototypical fermionic model of topological properties. To set the stage, we firstly compute the single-particle spectral function, the local density, and the boundary charge in the absence of interactions. The boundary charge is fully determined by bulk properties indicating a bulk-boundary correspondence. In a large parameter regime it agrees with the one obtained from an effective low-energy theory (arXiv:2004.00463). Secondly, we investigate the robustness of our results towards two-particle interactions. To resum the series of leading logarithms for small gaps, which dismantle plain perturbation theory in the interaction, we use an essentially analytical renormalization group approach. It is controlled for small interactions and can directly be applied to the microscopic lattice model. We benchmark the results against numerical density matrix renormalization group data. The main interaction effect in the bulk is a power-law renormalization of the gap with an interaction dependent exponent. The important characteristics of the boundary charge are unaltered and can be understood from the renormalized bulk properties, elevating the bulk-boundary correspondence to the interacting regime. This requires a consistent treatment not only of the low-energy gap renormalization but also of the high-energy band width one. In contrast to low-energy field theories our renormalization group approach also provides the latter. We show that the interaction spoils the relation between the bulk properties and the number of edge states, consistent with the observation that the Rice-Mele model with finite potential modulation does not reveal any zero-energy edge states.

Time-dependent spintronic anisotropy in magnetic molecules

We theoretically study the quench dynamics of induced anisotropy of a large-spin magnetic molecule coupled to spin-polarized ferromagnetic leads. The real-time evolution is calculated by means of the time-dependent density-matrix numerical renormalization group method implemented within the matrix product states framework, which takes into account all correlations in very accurate manner. We determine the system's response to a quench in the spin-dependent coupling to ferromagnetic leads. In particular, we focus on the transient dynamics associated with crossing from the weak to the strong coupling regime, where the Kondo correlations become important. The dynamics is examined by calculating the time-dependent expectation values of the spin-quadrupole moment and the associated spin operators. We identify the relevant time scales describing the quench dynamics and determine the influence of the molecule's effective exchange coupling and leads spin-polarization on the dynamical behavior of the system. Furthermore, the generalization of our predictions for large values of molecule's spin is considered. Finally, we analyze the effect of finite temperature and show that it gives rise to a reduction of magnetic anisotropy by strong suppression of the time-dependent spin-quadrupole moment due to thermal fluctuations.

Dipole oscillations of fermionic quantum gases along the BEC-BCS crossover in disordered potentials

We investigate dipole oscillations of ultracold Fermi gases along the BEC-BCS crossover through disordered potentials. We observe a disorder-induced damping of oscillations as well as a change of the fundamental Kohn-mode frequency. The measurement results are compared to numerical density matrix renormalization group calculations as well as to a three-dimensional simulation of non-interacting fermions. Experimentally, we find a disorder-dependent damping, which grows approximately with the second power of the disorder strength. Moreover, we observe experimentally a change of oscillation frequency which deviates from the expected behavior of a damped harmonic oscillator on a percent level. While this behavior is qualitatively expected from the theoretical models used, quantitatively the experimental observations show a significantly stronger effect than predicted by theory. Furthermore, while the frequency shift seems to scale differently with interaction strength in the BEC versus BCS regime, the damping coefficient apparently decreases with the strength of interaction, but not with the sign, which changes for BEC and BCS type Fermi gases. This is surprising, as the dominant damping mechanisms are expected to be different in the two regimes.

Bell Correlations at Ising Quantum Critical Points.

When a collection of distant observers share an entangled quantum state, the statistical correlations among their measurements may violate a many-body Bell inequality, demonstrating a nonlocal behavior. Focusing on the Ising model in a transverse field with power-law (1/r^{α}) ferromagnetic interactions, we show that a permutationally invariant Bell inequality based on two-body correlations is violated in the vicinity of the quantum-critical point. This observation, obtained via analytical spin-wave calculations and numerical density-matrix renormalization group computations, is traced back to the squeezing of collective-spin fluctuations generated by quantum-critical correlations. We observe a maximal violation for infinite-range interactions (α=0), namely, when interactions and correlations are themselves permutationally invariant.

Topological Mott insulator with bosonic edge modes in one-dimensional fermionic superlattices

We investigate topological phase transitions driven by interaction and identify a novel topological Mott insulator state in one-dimensional fermionic optical superlattices through numerical density matrix renormalization group (DMRG) method. Remarkably, the low-energy edge excitations change from spin-1/2 fermionic single-particle modes to spin-1 bosonic collective modes across the phase transition. Due to spin-charge separation, the low-energy theory is governed by an effective spin superexchange model, whereas the charge degree of freedom is fully gapped out. Such topological Mott state can be characterized by a spin Chern number and gapless magnon modes protected by a finite spin gap. The proposed experimental setup is simple and may pave the way for the experimental observation of exotic topological Mott states.

Interaction-dependent anisotropy of fractional quantum Hall states

A fractional quantum Hall (FQH) system with broken rotational symmetry exploits its geometric degree of freedom to minimize its ground state energy. The mass anisotropy of bare particles interacting isotropically is partially inherited by the many-body FQH state, and the extent to which it does so depends on the type of interaction, filling fraction and ground state phase. Using numerical infinite density matrix renormalization group simulations, we investigate the transference of elliptical ($C_2$-symmetric) anisotropy from the band mass of the bare particles to the FQH states, for various power law interactions. We map out the response of FQH states to small anisotropy as a function of power law exponent, filling, and statistics (bosonic or fermionic) of the constituents. Interestingly, we find a non-analyticity in the linear response of the FQH state at a special filling-dependent value of the power law exponent, above which the interaction effectively becomes zero-range (point-like). We also investigate the the effect of $C_4$-symmetric band distortions, where we observe a strikingly different dependence on filling.

Dynamics of a mobile impurity in a two-leg bosonic ladder

We have analyzed the behavior of a mobile quantum impurity in a bath formed by a two-leg bosonic ladder by a combination of field theory (Tomonaga-Luttinger liquid) and numerical (Density Matrix Renormalization Group) techniques. Computing the Green's function of the impurity as a function of time at different momenta, we find a power law decay at zero momentum, which signals the breakdown of any quasi-particle description of the impurity motion. We compute the exponent both for the limits of weak and strong impurity-bath interactions. At small impurity-bath interaction, we find that the impurity experiences the ladder as a single channel one-dimensional bath, but effective coupling is reduced by a factor of $\sqrt 2$, thus the impurity is less mobile in the ladder compared to a one dimensional bath. We compared the numerical results for the exponent at zero momentum with a semi-analytical expression that was initially established for the chain and find excellent agreement without adjustable parameters. We analyze the dependence of the exponent in the transverse hopping in the bath and find surprisingly an increase of the exponent at variance with the naive extrapolation of the single channel regime. We study the momentum dependence of the impurity Green's function and find that, as for the single chain, two different regime of motion exist, one dominated by infrared metatrophy and a more conventional polaronic behavior. We compute the critical momentum between these two regimes and compare with prediction based on the structure factor of the bath. In the polaronic regime we also compute numerically the lifetime of the polaron. Finally we discuss how our results could be measured in cold atomic experiments.