arXiv:1702.03650 · Physical Review Letters

Sliced Basis Density Matrix Renormalization Group for Electronic Structure

Edwin Miles Stoudenmire, Steven R. White

Published Jul 24, 2017• 40 citations

Abstract

We introduce a hybrid approach to applying the density matrix renormalization group to continuous systems, combining a grid approximation along one direction with a finite Gaussian basis set for the remaining two directions. This approach is especially useful for chainlike molecules, where the grid is used in the long direction. For hydrogen chain systems, the computational time scales approximately linearly with the number of atoms, as we show with near-exact minimal basis set calculations with up to 1000 atoms. The linear scaling comes from both the localization of the basis and a compression method for the long-ranged two-electron interaction. For shorter hydrogen chains, we show results with up to triple-ζ bases.

DOI: 10.1103/physrevlett.119.046401

Related papers

Quantum Chemistry Density Matrix Renormalization Group in the Discrete Variable Representation

Bing Gu, Jiajun Ren, Junzhe Zhang

We present a numerical implementation of the quantum chemistry density matrix renormalization group (DMRG) using the hybrid discrete variable representation (DVR)/Gaussian basis set. The z -axis of real space is discretized by a DVR basis set and each transversal plane is described by the eigenstates of the transversal core Hamiltonian, represented in a set of primitive Gaussian basis functions. Such a hybrid basis can reduce the computation of two-electron repulsion integrals. One main advantage of using the local DVR basis sets over finite-difference methods for real space discretization is that the kinetic energy operator matrix elements can be computed exactly. This method is first applied to a one-dimensional pseudo-hydrogen chain under screened Coulomb potential and then to a realistic hydrogen chain, whereby the DMRG/DVR results are comparable in accuracy to the full configuration interaction results.

Journal of Chemical Theory and ComputationJul 2, 20250 citations

One-dimensional Continuum Electronic Structure Calculations with the Density Matrix Renormalization Group

Edwin Miles Stoudenmire, Lucas O. Wagner, Steven R. White +1 more

We extend the density matrix renormalization group to compute exact ground states of continuum many-electron systems in one dimension with long-range interactions. We find the exact ground state of a chain of 100 strongly correlated artificial hydrogen atoms. The method can be used to simulate 1d cold atom systems and to study density functional theory in an exact setting. To illustrate, we find an interacting, extended system which is an insulator but whose Kohn-Sham system is metallic.

Jul 12, 20111 citation

Multisliced gausslet basis sets for electronic structure

Steven R. White, Edwin Miles Stoudenmire

We introduce highly local basis sets for electronic structure which are very efficient for correlation calculations near the complete basis set limit. Our approach is based on gausslets, recently introduced waveletlike smooth orthogonal functions. We adapt the gausslets to particular systems using one-dimensional (1D) coordinate transformations, putting more basis functions near the nuclei, while maintaining orthogonality. Three-dimensional basis functions are composed out of products of the 1D functions in an efficient way called multislicing. We demonstrate these bases with both Hartree-Fock and density matrix renormalization group calculations on hydrogen chain systems. With both methods, we can go to higher accuracy in the complete basis set limit than is practical for conventional Gaussian basis sets, with errors near 0.1 mH per atom.

arXiv:1809.10258Physical review. B./Physical review. BFeb 11, 201929 citations

One-Dimensional Continuum Electronic Structure with the Density-Matrix Renormalization Group and Its Implications for Density-Functional Theory

Edwin Miles Stoudenmire, Lucas O. Wagner, Steven R. White +1 more

We extend the density matrix renormalization group to compute exact ground states of continuum many-electron systems in one dimension with long-range interactions. We find the exact ground state of a chain of 100 strongly correlated artificial hydrogen atoms. The method can be used to simulate 1D cold atom systems and to study density-functional theory in an exact setting. To illustrate, we find an interacting, extended system which is an insulator but whose Kohn-Sham system is metallic.

arXiv:1107.2394Physical Review LettersAug 1, 201287 citations

Constructing Hubbard models for the hydrogen chain using sliced-basis density matrix renormalization group

Randy C. Sawaya, Steven R. White

Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.

Physical review. B./Physical review. BJan 31, 20229 citations

Direct density matrix renormalization group approaches for strong correlation effects in quantum chemistry

Philipp Bertone Snajberk

An implementation of the single-site density matrix renormalization group (DMRG) approach for the quantum-chemical application to calculate approximate as well as exact electronic ground-state energies and wave functions of atoms and molecules is presented.The algorithm is formulated in terms of the variational and sizeconsistent class of wave functions, matrix product states, and their operator analogues, matrix product operators.A procedure to construct quantum-chemical Hamiltonian operators with long-range Coulomb interactions as matrix product operators is developed.The many-electron Hamiltonian operators are represented as sums and products of one-electron operators.For the efficient matrix product operator construction of quantum-chemical Hamiltonians, different basis spaces including discrete position and momentum space as well as molecular and atomic orbitals are investigated.Discrete position space yields the best matrix product operator size scaling, but introduces artificial long-range couplings in the kinetic energy term of the Hamiltonian which are a consequence of the DMRG algorithm in terms of matrix product states.The huge amount of discrete space points needed for an accurate description poses a problem because of the electron-electron interaction.Therefore, a momentum space or Fourier series expansion of the electron-electron interaction is shown which makes it possible to compress its matrix product operator representation.DMRG applied to quantum-chemical Hamiltonians in matrix product operator form can serve as means to directly (without mean-field treatment in advance) investigate quantum-chemical problems and find the analytically exact ground-state solution of small atoms and molecules.The next basis space analyzed in detail is the atomic orbital basis.The matrix product operator construction and the direct application of the DMRG method in this basis is examined.A physical compression scheme based on a singular value decomposition of the electron-electron interaction in the atomic orbital basis is demonstrated.Furthermore, DMRG calculations in the atomic orbital basis for certain systems where mean-field theory yields qualitatively wrong results and where higher-order approximative correction methods like truncated configuration interaction and coupled cluster cannot fully capture strong electron correlation effects are compared to results obtained in the molecular orbital basis.The comparison shows that atomic orbitals describe the physical situation well and thus lead to a tremendous reduction of the matrix product state dimension and the computational effort.i

Jan 1, 20171 citation