Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1809.05808 · Physical Review A

Density-matrix-renormalization-group study of a one-dimensional diatomic molecule beyond the Born-Oppenheimer approximation

Mingru Yang, Steven R. White

Published Feb 21, 2019• 7 citations
arXivHTMLPDF

Abstract

We study one-dimensional (1D) models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate density-matrix renormalization group efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds.

DOI
10.1103/physreva.99.022509

Related papers

Improving matrix product state methods for long-range interactions

Mingru Yang

We study one dimensional models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate DMRG efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin-triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds, given the two electrons in their singlet state in both cases. \nWe propose an improved scheme to do the time dependent variational principle (TDVP) in finite matrix product states (MPS) for two-dimensional systems or one-dimensional systems with long range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state-averaging with global Krylov vectors. We show that the projection error is significantly reduced so that precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.\n \nWe apply our improved TDVP method to investigate the spin squeezing dynamics of the XXZ model of $1/r^\\alpha$ interaction in two dimension. Comparing with the spin squeezing parameter and other observables obtained from discrete truncated Wigner approximation (DTWA), we verify the validity of DTWA and unveil the potential for this method to study dynamics of large-scale spin systems. Our results confirm the existence of large collective regime when $\\alpha > 2$, which can be a guide for future experimental realizations. Combined with the purification method, our improved TDVP method is proved to be also useful to study the thermalization of the long-range interacting system.

•Jan 1, 2020•0 citations

Sliced Basis Density Matrix Renormalization Group for Electronic Structure

Edwin Miles Stoudenmire, Steven R. White

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.

arXiv:1702.03650•Physical Review Letters•Jul 24, 2017•40 citations

Reference electronic structure calculations in one dimension

Lucas O. Wagner, Edwin Miles Stoudenmire, Kieron Burke +1 more

Large strongly correlated systems provide a challenge to modern electronic structure methods, because standard density functionals usually fail and traditional quantum chemical approaches are too demanding. The density-matrix renormalization group method, an extremely powerful tool for solving such systems, has recently been extended to handle long-range interactions on real-space grids, but is most efficient in one dimension where it can provide essentially arbitrary accuracy. Such 1d systems therefore provide a theoretical laboratory for studying strong correlation and developing density functional approximations to handle strong correlation, if they mimic three-dimensional reality sufficiently closely. We demonstrate that this is the case, and provide reference data for exact and standard approximate methods, for future use in this area.

arXiv:1202.4788•Physical Chemistry Chemical Physics•Jan 1, 2012•81 citations

Spin-adapted density matrix renormalization group algorithms for quantum chemistry

Sandeep Sharma, Garnet Kin-Lic Chan

We extend the spin-adapted density matrix renormalization group (DMRG) algorithm of McCulloch and Gulacsi [Europhys. Lett. 57, 852 (2002)] to quantum chemical Hamiltonians. This involves using a quasi-density matrix, to ensure that the renormalized DMRG states are eigenfunctions of Ŝ(2), and the Wigner-Eckart theorem, to reduce overall storage and computational costs. We argue that the spin-adapted DMRG algorithm is most advantageous for low spin states. Consequently, we also implement a singlet-embedding strategy due to Tatsuaki [Phys. Rev. E 61, 3199 (2000)] where we target high spin states as a component of a larger fictitious singlet system. Finally, we present an efficient algorithm to calculate one- and two-body reduced density matrices from the spin-adapted wavefunctions. We evaluate our developments with benchmark calculations on transition metal system active space models. These include the Fe(2)S(2), [Fe(2)S(2)(SCH(3))(4)](2-), and Cr(2) systems. In the case of Fe(2)S(2), the spin-ladder spacing is on the microHartree scale, and here we show that we can target such very closely spaced states. In [Fe(2)S(2)(SCH(3))(4)](2-), we calculate particle and spin correlation functions, to examine the role of sulfur bridging orbitals in the electronic structure. In Cr(2) we demonstrate that spin-adaptation with the Wigner-Eckart theorem and using singlet embedding can yield up to an order of magnitude increase in computational efficiency. Overall, these calculations demonstrate the potential of using spin-adaptation to extend the range of DMRG calculations in complex transition metal problems.

arXiv:1408.5039•The Journal of Chemical Physics•Mar 28, 2012•295 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, 2011•1 citation

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.2394•Physical Review Letters•Aug 1, 2012•87 citations
We study one-dimensional (1D) models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate density-matrix renormalization group efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds.
We study one dimensional models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate DMRG efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin-triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds, given the two electrons in their singlet state in both cases. \nWe propose an improved scheme to do the time dependent variational principle (TDVP) in finite matrix product states (MPS) for two-dimensional systems or one-dimensional systems with long range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state-averaging with global Krylov vectors. We show that the projection error is significantly reduced so that precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.\n \nWe apply our improved TDVP method to investigate the spin squeezing dynamics of the XXZ model of $1/r^\\alpha$ interaction in two dimension. Comparing with the spin squeezing parameter and other observables obtained from discrete truncated Wigner approximation (DTWA), we verify the validity of DTWA and unveil the potential for this method to study dynamics of large-scale spin systems. Our results confirm the existence of large collective regime when alpha>2\\alpha > 2alpha>2, which can be a guide for future experimental realizations. Combined with the purification method, our improved TDVP method is proved to be also useful to study the thermalization of the long-range interacting system.
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.
Large strongly correlated systems provide a challenge to modern electronic structure methods, because standard density functionals usually fail and traditional quantum chemical approaches are too demanding. The density-matrix renormalization group method, an extremely powerful tool for solving such systems, has recently been extended to handle long-range interactions on real-space grids, but is most efficient in one dimension where it can provide essentially arbitrary accuracy. Such 1d systems therefore provide a theoretical laboratory for studying strong correlation and developing density functional approximations to handle strong correlation, if they mimic three-dimensional reality sufficiently closely. We demonstrate that this is the case, and provide reference data for exact and standard approximate methods, for future use in this area.
We extend the spin-adapted density matrix renormalization group (DMRG) algorithm of McCulloch and Gulacsi [Europhys. Lett. 57, 852 (2002)] to quantum chemical Hamiltonians. This involves using a quasi-density matrix, to ensure that the renormalized DMRG states are eigenfunctions of Ŝ(2), and the Wigner-Eckart theorem, to reduce overall storage and computational costs. We argue that the spin-adapted DMRG algorithm is most advantageous for low spin states. Consequently, we also implement a singlet-embedding strategy due to Tatsuaki [Phys. Rev. E 61, 3199 (2000)] where we target high spin states as a component of a larger fictitious singlet system. Finally, we present an efficient algorithm to calculate one- and two-body reduced density matrices from the spin-adapted wavefunctions. We evaluate our developments with benchmark calculations on transition metal system active space models. These include the Fe(2)S(2), [Fe(2)S(2)(SCH(3))(4)](2-), and Cr(2) systems. In the case of Fe(2)S(2), the spin-ladder spacing is on the microHartree scale, and here we show that we can target such very closely spaced states. In [Fe(2)S(2)(SCH(3))(4)](2-), we calculate particle and spin correlation functions, to examine the role of sulfur bridging orbitals in the electronic structure. In Cr(2) we demonstrate that spin-adaptation with the Wigner-Eckart theorem and using singlet embedding can yield up to an order of magnitude increase in computational efficiency. Overall, these calculations demonstrate the potential of using spin-adaptation to extend the range of DMRG calculations in complex transition metal problems.