Density Matrix Renormalization Group Algorithm Research Articles

Terminable Transitions in a Topological Fermionic Ladder.

Interacting fermionic ladders are versatile platforms to study quantum phases of matter, such as different types of Mott insulators. In particular, there are D-Mott and S-Mott states that hold preformed fermion pairs and become paired-fermion liquids upon doping (d wave and s wave, respectively). We show that the D-Mott and S-Mott phases are in fact two facets of the same topological phase and that the transition between them is terminable. These results provide a quantum analog of the well-known terminable liquid-to-gas transition. However, the phenomenology we uncover is even richer, as the order of the transition may alternate between continuous and first order, depending on the interaction details. Most importantly, the terminable transition is robust in the sense that it is guaranteed to appear for weak, but arbitrary couplings. We discuss a minimal model where some analytical insights can be obtained, a generic model where the effect persists; and a model-independent field-theoretical study demonstrating the general phenomenon. The role of symmetry and the edge states is briefly discussed. The numerical results are obtained using the variational uniform matrix-product state (VUMPS) formalism for infinite systems, as well as the density-matrix renormalization group (DMRG) algorithm for finite systems.

Zero-Field Splitting Tensor of the Triplet Excited States of Aromatic Molecules: A Valence Full-π Complete Active Space Self-Consistent Field Study.

A method to predict the D tensor in the molecular frame with multiconfigurational wave functions in large active space was proposed, and the spin properties of the lowest triplets of aromatic molecules were examined with full-π active space; such calculations were challenging because the size of active space grows exponentially with the number of π electrons. In this method, the exponential growth of complexity is resolved by the density matrix renormalization group (DMRG) algorithm. From the D tensor, we can directly determine the direction of the magnetic axes and the ZFS parameters, D- and E-values, of the phenomenological spin Hamiltonian with their signs, which are not usually obtained in ESR experiments. The method using the DMRG-CASSCF wave function can give correct results even when the sign of D- and E-values is sensitive to the accuracy of the prediction of the D tensor and existing methods fail to predict the correct magnetic axes.

Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase.

The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).

Block2: A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.

Flexible DMRG-Based Framework for Anharmonic Vibrational Calculations.

We present a novel formulation of the vibrational density matrix renormalization group (vDMRG) algorithm tailored to strongly anharmonic molecules described by general, high-dimensional model representations of potential energy surfaces. For this purpose, we extend the vDMRG framework to support vibrational Hamiltonians expressed in the so-called n-mode second-quantization formalism. The resulting n-mode vDMRG method offers full flexibility with respect to both the functional form of the PES and the choice of the single-particle basis set. We leverage this framework to apply, for the first time, vDMRG based on an anharmonic modal basis set optimized with the vibrational self-consistent field algorithm on an on-the-fly constructed PES. We also extend the n-mode vDMRG framework to include excited-state-targeting algorithms in order to efficiently calculate anharmonic transition frequencies. We demonstrate the capabilities of our novel n-mode vDMRG framework for methyloxirane, a challenging molecule with 24 coupled vibrational modes.

Chebyshev pseudosite matrix product state approach for the spectral functions of electron-phonon coupling systems

The electron-phonon ($e$-ph) coupling system often has a large number of phonon degrees of freedom, whose spectral functions are numerically difficult to compute using matrix product state (MPS) formalisms. To solve this problem, we propose a new and practical method that combines the Chebyshev MPS and the pseudosite density matrix renormalization group (DMRG) algorithm. The Chebyshev vector is represented by a pseudosite MPS with global $U(1)$ fermion symmetry, which maps $2^{N_p}$ bosonic degrees of freedom onto $N_p$ pseudosites, each with two states. This approach can handle arbitrary $e$-ph coupling Hamiltonians where pseudosite DMRG performs efficiently. We use this method to study the spectral functions of the doped extended Hubbard-Holstein model in a regime of strong Coulomb repulsion, which has not been studied extensively before. Key features of the excitation spectra are captured at a modest computational cost. Our results show that weak extended $e$-ph couplings can increase the spectral weight of the holon-folding branch at low phonon frequencies, in agreement with angle-resolved photoemission observations on one-dimensional (1D) cuprates.

Density-Matrix Renormalization Group Algorithm for Simulating Quantum Circuits with a Finite Fidelity

We develop a density-matrix renormalization group (DMRG) algorithm for the simulation of quantum circuits. This algorithm can be seen as the extension of time-dependent DMRG from the usual situation of hermitian Hamiltonian matrices to quantum circuits defined by unitary matrices. For small circuit depths, the technique is exact and equivalent to other matrix product state (MPS) based techniques. For larger depths, it becomes approximate in exchange for an exponential speed up in computational time. Like an actual quantum computer, the quality of the DMRG results is characterized by a finite fidelity. However, unlike a quantum computer, the fidelity depends strongly on the quantum circuit considered. For the most difficult possible circuit for this technique, the so-called "quantum supremacy" benchmark of Google Inc. , we find that the DMRG algorithm can generate bit strings of the same quality as the seminal Google experiment on a single computing core. For a more structured circuit used for combinatorial optimization (Quantum Approximate Optimization Algorithm or QAOA), we find a drastic improvement of the DMRG results with error rates dropping by a factor of 100 compared with random quantum circuits. Our results suggest that the current bottleneck of quantum computers is their fidelities rather than the number of qubits.

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules.

We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.

Time-reversal symmetry adaptation in relativistic density matrix renormalization group algorithm.

In the nonrelativistic Schrödinger equation, the total spin S and spin projection M are good quantum numbers. In contrast, spin symmetry is lost in the presence of spin-dependent interactions, such as spin-orbit couplings in relativistic Hamiltonians. Therefore, the relativistic density matrix renormalization group algorithm (R-DMRG) only employing particle number symmetry is much more expensive than nonrelativistic DMRG. In addition, artificial breaking of Kramers degeneracy can happen in the treatment of systems with an odd number of electrons. To overcome these issues, we propose time-reversal symmetry adaptation for R-DMRG. Since the time-reversal operator is antiunitary, this cannot be simply achieved in the usual way. We introduce a time-reversal symmetry-adapted renormalized basis and present strategies to maintain the structure of basis functions during the sweep optimization. With time-reversal symmetry adaptation, only half of the renormalized operators are needed, and the computational costs of Hamiltonian-wavefunction multiplication and renormalization are reduced by half. The present construction of the time-reversal symmetry-adapted basis also directly applies to other tensor network states without loops.

Root- N Krylov-space correction vectors for spectral functions with the density matrix renormalization group

We propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov space decomposition to compute the correction-vectors. Our approach entails the calculation of the root-N (N=2 is the standard square root) of the Hamiltonian propagator using Krylov space decomposition, and repeating this procedure N times to obtain the actual correction-vector. We show that our method greatly alleviates the burden of keeping a large bond dimension at large target frequencies, a problem found with conventional correction-vector DMRG, while achieving better computational performance at large N. We apply our method to spin and charge spectral functions of t-J and Hubbard models in the challenging two-leg ladder geometry, and provide evidence that the root-N approach reaches a much improved resolution compared to conventional correction-vector.

Bond order via cavity-mediated interactions

We numerically study the phase diagram of bosons tightly trapped in the lowest band of an optical lattice and dispersively coupled to a single-mode cavity field. The dynamics is encompassed by an extended Bose-Hubbard model. Here, the cavity-mediated interactions are described by a two-body potential term with a global range and by a correlated tunneling term where the hopping amplitude depends on a global observable. We determine the ground state properties in one dimension by means of the density matrix renormalization group algorithm, focusing on the effects due to the correlated tunneling. The latter is responsible for the onset of bond orders, manifesting in one insulating and two gapless bond ordered phases. We discuss the resulting phases for different geometries that correspond to different relative strengths of the correlated tunneling coefficient. We finally analyze the scaling of the entanglement entropy in the gapless bond ordered phases that appear entirely due to global interactions and determine the corresponding central charges.

Computing time-periodic steady-state currents via the time evolution of tensor network states.

We present an approach based upon binary tree tensor network (BTTN) states for computing steady-state current statistics for a many-particle 1D ratchet subject to volume exclusion interactions. The ratcheted particles, which move on a lattice with periodic boundary conditions subject to a time-periodic drive, can be stochastically evolved in time to sample representative trajectories via a Gillespie method. In lieu of generating realizations of trajectories, a BTTN state can variationally approximate a distribution over the vast number of many-body configurations. We apply the density matrix renormalization group algorithm to initialize BTTN states, which are then propagated in time via the time-dependent variational principle (TDVP) algorithm to yield the steady-state behavior, including the effects of both typical and rare trajectories. The application of the methods to ratchet currents is highlighted, but the approach extends naturally to other interacting lattice models with time-dependent driving. Although trajectory sampling is conceptually and computationally simpler, we discuss situations for which the BTTN TDVP strategy can be beneficial.

Quantum phases of dipolar bosons in one-dimensional optical lattices

We theoretically analyze the phase diagram of a quantum gas of bosons that interact via repulsive dipolar interactions. The bosons are tightly confined by an optical lattice in a quasi-one-dimensional geometry. In the single-band approximation, their dynamics is described by an extended Bose-Hubbard model where the relevant contributions of the dipolar interactions consist of density-density repulsion and correlated tunneling terms. We evaluate the phase diagram for unit density using numerical techniques based on the density-matrix renormalization group algorithm. Our results predict that correlated tunneling can significantly modify the parameter range of the topological insulator phase. At vanishing values of the onsite interactions, moreover, correlated tunneling promotes the onset of a phase with a large number of low-energy metastable configurations.

Localization of a mobile impurity interacting with an Anderson insulator

Thermalizing and localized many-body quantum systems present two distinct dynamical phases of matter. Recently, the fate of a localized system coupled to a thermalizing system viewed as a quantum bath received significant theoretical and experimental attention. In this work, we study a mobile impurity, representing a small quantum bath, that interacts locally with an Anderson insulator with a finite density of localized particles. Using static Hartree approximation to obtain an effective disorder strength, we formulate an analytic criterion for the perturbative stability of the localization. Next, we use an approximate dynamical Hartree method and the quasi-exact time-evolved block decimation (TEBD) algorithm to study the dynamics of the system. We find that the dynamical Hartree approach which completely ignores entanglement between the impurity and localized particles predicts the delocalization of the system. In contrast, the full numerical simulation of the unitary dynamics with TEBD suggests the stability of localization on numerically accessible timescales. Finally, using an extension of the density matrix renormalization group algorithm to excited states (DMRG-X), we approximate the highly excited eigenstates of the system. We find that the impurity remains localized in the eigenstates and entanglement is enhanced in a finite region around the position of the impurity, confirming the dynamical predictions. Dynamics and the DMRG-X results provide compelling evidence for the stability of localization.

Explicitly Correlated Electronic Structure Calculations with Transcorrelated Matrix Product Operators.

In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron-electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [ J. Chem. Phys. 2020, 153, 164115] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian.

Density matrix renormalization group algorithm for mixed quantum states

Density Matrix Renormalization Group (DMRG) algorithm has been extremely successful for computing the ground states of one-dimensional quantum many-body systems. For problems concerned with mixed quantum states, however, it is less successful in that either such an algorithm does not exist yet or that it may return unphysical solutions. Here we propose a positive matrix product ansatz for mixed quantum states which preserves positivity by construction. More importantly, it allows to build a DMRG algorithm which, the same as the standard DMRG for ground states, iteratively reduces the global optimization problem to local ones of the same type, with the energy converging monotonically in principle. This algorithm is applied for computing both the equilibrium states and the non-equilibrium steady states, and its advantages are numerically demonstrated.

Critical properties of quantum three- and four-state Potts models with boundaries polarized along the transverse field

By computing the low-lying energy excitation spectra with the density matrix renormalization group algorithm we show that boundaries polarized in the direction of the transverse field lead to scale-invariant conformal towers of states at the critical point of the quantum four-state Potts model - a special symmetric case of the Ashkin-Teller model. Furthermore, by direct comparison of the excitation spectra we phenomenologically establish the duality between the transverse-polarized and three-state-mixed boundary conditions at the four-state Potts critical point. Finally, for completeness, we verify that in the quantum three-state Potts model the "new" boundary conditions dual to the mixed ones can be realized by polarizing edge spins along the transverse field.

237Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes.

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. A rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

Long-range order and quantum criticality in antiferromagnetic chains with long-range staggered interactions.

We study quantum phase transitions in Heisenberg antiferromagnetic chains with a staggered power-law decaying long-range interactions. Employing the density-matrix renormalization group (DMRG) algorithm and the fidelity susceptibility as the criticality measure, we establish more accurate values of quantum critical points than the results obtained from the spin-wave approximation, quantum Monte Carlo, and DMRG in literatures. The deviation is especially evident for strong long-range interactions. We extend isotropic long-range interactions to the anisotropic cases and find that kaleidoscope of quantum phases emerge from the interplay of anisotropy of the long-range exchange interaction and symmetry breaking. We demonstrate nonfrustrating long-range interactions induce the true long-range order in Heisenberg antiferromagnetic chains with a continuous symmetry breaking, lifting the restrictions imposed by the Mermin-Wagner theorem.

FFLO state driven by quasiperiodic Zeeman field and its transition to localized states

We numerically study the ground-state properties of a one-dimensional lattice model in the presence of a quasiperiodic Zeeman field, by means of the density matrix renormalization group algorithm. We map out a global phase diagram, which encloses a Bardeen-Cooper-Schrieffer (BCS) phase, a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase where the Cooper pairs obtain a finite center-of-mass momentum, and a localized phase. The FFLO phase is promoted by the quasiperiodic Zeeman field, which demonstrates a power-law pairing correlation with a critical exponent ${\ensuremath{\eta}}_{\mathrm{FFLO}}<1$ decaying much slower than the charge density and spin density correlations. By tuning the interaction and the filling factor, we find the FFLO phase is suppressed in the strong interaction and low filling regions. A large Zeeman field destroys the FFLO state and drives a superconductor-insulator transition. Furthermore, we propose this FFLO phase and the superconductor-insulator transition could be observed in the optical lattice experiment.