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:1905.01558 · Journal of Chemical Theory and Computation

Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method

Leon Freitag, Yingjin Ma, Alberto Baiardi, Stefan Knecht, Markus Reiher

Published Oct 31, 2019• 24 citations
arXivHTMLPDF

Abstract

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.

DOI
10.1021/acs.jctc.9b00969

Related papers

Algorithm for analytic nuclear energy gradients of state averaged DMRG-CASSCF theory with newly derived coupled-perturbed equations

Tsubasa Iino, Toru Shiozaki, Takeshi Yanai

We present an algorithm for evaluating analytic nuclear energy gradients of the state-averaged density matrix renormalization group complete-active-space self-consistent field (SA-DMRG-CASSCF) theory based on the newly derived coupled-perturbed (CP) DMRG-CASSCF equations. The Lagrangian for the conventional SA-CASSCF analytic gradient theory is extended to the SA-DMRG-CASSCF variant that can fully consider a whole set of constraints on the parameters of multi-root canonical matrix product states formed at all the DMRG block configurations. An efficient algorithm to solve the CP-DMRG-CASSCF equations for determining the multipliers was developed. The complexity of the resultant analytic gradient algorithm is overall the same as that of the unperturbed SA-DMRG-CASSCF algorithm. In addition, a reduced-scaling approach was developed to directly compute the SA reduced density matrices (SA-RDMs) and their perturbed ones without calculating separate state-specific RDMs. As part of our implementation scheme, we neglect the term associated with the constraint on the active orbitals in terms of the active-active rotation in the Lagrangian. Thus, errors from the true analytic gradients may be caused in this scheme. The proposed gradient algorithm was tested with the spin-adapted implementation by checking how accurately the computed analytic energy gradients reproduce numerical gradients of the SA-DMRG-CASSCF energies using a common number of renormalized bases. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized bases, which is required to attain adequate accuracy in DMRG's total energies.

The Journal of Chemical Physics•Dec 16, 2022•5 citations

The density matrix renormalization group self-consistent field method: Orbital optimization with the density matrix renormalization group method in the active space

Dominika Zgid, Marcel Nooijen

We present the density matrix renormalization group self-consistent field (DMRG-SCF) approach that is analogous to the complete active space self-consisted field (CASSCF) method but instead of using for the description of the active space the full configuration interaction (FCI) method, the DMRG-SCF uses the density matrix renormalization group (DMRG) method. The DMRG-SCF approach, similarly to CASSCF, properly describes the multiconfigurational character of the wave function but avoids the exponential scaling of the FCI method and replaces it with a polynomial scaling. Hence, calculations for a larger number of orbitals and electrons in the active space are possible since the DMRG method provides an efficient tool to automatically select from the full Hilbert space the many-body contracted basis states that are the most important for the description of the wave function.

The Journal of Chemical Physics•Apr 11, 2008•202 citations

Density matrix renormalization group (DMRG) method as a common tool for large active-space CASSCF/CASPT2 calculations

Naoki Nakatani, Sheng Guo

This paper describes an interface between the density matrix renormalization group (DMRG) method and the complete active-space self-consistent field (CASSCF) method and its analytical gradient, as well as an extension to the second-order perturbation theory (CASPT2) method. This interfacing allows large active-space multi-reference computations to be easily performed. The interface and its extension are both implemented in terms of reduced density matrices (RDMs) which can be efficiently computed via the DMRG sweep algorithm. We also present benchmark results showing that, in practice, the DMRG-CASSCF calculations scale with active-space size in a polynomial manner in the case of quasi-1D systems. Geometry optimization of a binuclear iron-sulfur cluster using the DMRG-CASSCF analytical gradient is demonstrated, indicating that the inclusion of the valence p-orbitals of sulfur and double-shell d-orbitals of iron lead to non-negligible changes in the geometry compared to the results of small active-space calculations. With the exception of the selection of M values, many computational settings in these practical DMRG calculations have been tuned and black-boxed in our interface, and so the resulting DMRG-CASSCF and DMRG-CASPT2 calculations are now available to novice users as a common tool to compute strongly correlated electronic wavefunctions.

The Journal of Chemical Physics•Mar 1, 2017•91 citations

Analytic energy gradients for the driven similarity renormalization group multireference second-order perturbation theory

Shuhe Wang, Chenyang Li, Francesco A. Evangelista

We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.

arXiv:2109.14548•Sep 29, 2021•2 citations

Advances in the density matrix renormalization group method for use in quantum chemistry

Dominika Zgid

Despite the success of modern quantum chemistry in predicting properties of organic molecules, the treatment of inorganic systems, which have many close lying states, remains out of quantitative reach for current methods. To treat non-dynamic correlation, we take advantage of the density matrix renormalization group (DMRG) method that has become very successful in the field of solid state physics. We present a detailed study of the DMRG method, and we pay special attention to the evolution of the understanding behind the mathematical structure of the DMRG wave function. Our primary goal is to develop a density matrix renormalization group self--consistent--field (DMRG-SCF) approach, analogous to the complete active space self--consistent field (CASSCF) method, but dealing with large active spaces that are too demanding for the full configuration interaction (FCI) method. \nAs a first step towards such a DMRG-SCF procedure, we present a spin-adapted DMRG algorithm designed to target spin- and spatial-symmetry states that are hard to obtain while using an unrestricted algorithm. \nOur next step is a modification of the DMRG algorithm to obtain decreasing energy at every step during the sweep. This monotonically convergent DMRG scheme lets us obtain the two-body density matrix as a by--product of the existing procedure without any additional cost in storage. Additionally, the two-body density matrix produced at convergence using this scheme is free from the N-representability problem that is present when the two--body density matrix is produced with the two-site DMRG scheme without additional storage cost. Finally, taking advantage of the modifications developed herein, we present results obtained from our DMRG-SCF method. Lastly, we discuss possible ways of incorporating dynamical correlation into the DMRG scheme, in order to obtain a modern multireference approach.

•May 1, 2008•1 citation

Second-Order Self-Consistent-Field Density-Matrix Renormalization Group

Yingjin Ma, Stefan Knecht, Sebastian F. Keller +1 more

We present a matrix-product state (MPS)-based quadratically convergent density-matrix renormalization group self-consistent-field (DMRG-SCF) approach. Following a proposal by Werner and Knowles (J. Chem. Phys. 1985, 82, 5053), our DMRG-SCF algorithm is based on a direct minimization of an energy expression which is correct to second order with respect to changes in the molecular orbital basis. We exploit a simultaneous optimization of the MPS wave function and molecular orbitals in order to achieve quadratic convergence. In contrast to previously reported (augmented Hessian) Newton-Raphson and superconfiguration-interaction algorithms for DMRG-SCF, energy convergence beyond a quadratic scaling is possible in our ansatz. Discarding the set of redundant active-active orbital rotations, the DMRG-SCF energy converges typically within two to four cycles of the self-consistent procedure.

arXiv:1611.05972•Journal of Chemical Theory and Computation•May 9, 2017•87 citations
We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.
We present an algorithm for evaluating analytic nuclear energy gradients of the state-averaged density matrix renormalization group complete-active-space self-consistent field (SA-DMRG-CASSCF) theory based on the newly derived coupled-perturbed (CP) DMRG-CASSCF equations. The Lagrangian for the conventional SA-CASSCF analytic gradient theory is extended to the SA-DMRG-CASSCF variant that can fully consider a whole set of constraints on the parameters of multi-root canonical matrix product states formed at all the DMRG block configurations. An efficient algorithm to solve the CP-DMRG-CASSCF equations for determining the multipliers was developed. The complexity of the resultant analytic gradient algorithm is overall the same as that of the unperturbed SA-DMRG-CASSCF algorithm. In addition, a reduced-scaling approach was developed to directly compute the SA reduced density matrices (SA-RDMs) and their perturbed ones without calculating separate state-specific RDMs. As part of our implementation scheme, we neglect the term associated with the constraint on the active orbitals in terms of the active-active rotation in the Lagrangian. Thus, errors from the true analytic gradients may be caused in this scheme. The proposed gradient algorithm was tested with the spin-adapted implementation by checking how accurately the computed analytic energy gradients reproduce numerical gradients of the SA-DMRG-CASSCF energies using a common number of renormalized bases. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized bases, which is required to attain adequate accuracy in DMRG's total energies.
We present the density matrix renormalization group self-consistent field (DMRG-SCF) approach that is analogous to the complete active space self-consisted field (CASSCF) method but instead of using for the description of the active space the full configuration interaction (FCI) method, the DMRG-SCF uses the density matrix renormalization group (DMRG) method. The DMRG-SCF approach, similarly to CASSCF, properly describes the multiconfigurational character of the wave function but avoids the exponential scaling of the FCI method and replaces it with a polynomial scaling. Hence, calculations for a larger number of orbitals and electrons in the active space are possible since the DMRG method provides an efficient tool to automatically select from the full Hilbert space the many-body contracted basis states that are the most important for the description of the wave function.
This paper describes an interface between the density matrix renormalization group (DMRG) method and the complete active-space self-consistent field (CASSCF) method and its analytical gradient, as well as an extension to the second-order perturbation theory (CASPT2) method. This interfacing allows large active-space multi-reference computations to be easily performed. The interface and its extension are both implemented in terms of reduced density matrices (RDMs) which can be efficiently computed via the DMRG sweep algorithm. We also present benchmark results showing that, in practice, the DMRG-CASSCF calculations scale with active-space size in a polynomial manner in the case of quasi-1D systems. Geometry optimization of a binuclear iron-sulfur cluster using the DMRG-CASSCF analytical gradient is demonstrated, indicating that the inclusion of the valence p-orbitals of sulfur and double-shell d-orbitals of iron lead to non-negligible changes in the geometry compared to the results of small active-space calculations. With the exception of the selection of M values, many computational settings in these practical DMRG calculations have been tuned and black-boxed in our interface, and so the resulting DMRG-CASSCF and DMRG-CASPT2 calculations are now available to novice users as a common tool to compute strongly correlated electronic wavefunctions.
We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.
Despite the success of modern quantum chemistry in predicting properties of organic molecules, the treatment of inorganic systems, which have many close lying states, remains out of quantitative reach for current methods. To treat non-dynamic correlation, we take advantage of the density matrix renormalization group (DMRG) method that has become very successful in the field of solid state physics. We present a detailed study of the DMRG method, and we pay special attention to the evolution of the understanding behind the mathematical structure of the DMRG wave function. Our primary goal is to develop a density matrix renormalization group self--consistent--field (DMRG-SCF) approach, analogous to the complete active space self--consistent field (CASSCF) method, but dealing with large active spaces that are too demanding for the full configuration interaction (FCI) method. \nAs a first step towards such a DMRG-SCF procedure, we present a spin-adapted DMRG algorithm designed to target spin- and spatial-symmetry states that are hard to obtain while using an unrestricted algorithm. \nOur next step is a modification of the DMRG algorithm to obtain decreasing energy at every step during the sweep. This monotonically convergent DMRG scheme lets us obtain the two-body density matrix as a by--product of the existing procedure without any additional cost in storage. Additionally, the two-body density matrix produced at convergence using this scheme is free from the N-representability problem that is present when the two--body density matrix is produced with the two-site DMRG scheme without additional storage cost. Finally, taking advantage of the modifications developed herein, we present results obtained from our DMRG-SCF method. Lastly, we discuss possible ways of incorporating dynamical correlation into the DMRG scheme, in order to obtain a modern multireference approach.
We present a matrix-product state (MPS)-based quadratically convergent density-matrix renormalization group self-consistent-field (DMRG-SCF) approach. Following a proposal by Werner and Knowles (J. Chem. Phys. 1985, 82, 5053), our DMRG-SCF algorithm is based on a direct minimization of an energy expression which is correct to second order with respect to changes in the molecular orbital basis. We exploit a simultaneous optimization of the MPS wave function and molecular orbitals in order to achieve quadratic convergence. In contrast to previously reported (augmented Hessian) Newton-Raphson and superconfiguration-interaction algorithms for DMRG-SCF, energy convergence beyond a quadratic scaling is possible in our ansatz. Discarding the set of redundant active-active orbital rotations, the DMRG-SCF energy converges typically within two to four cycles of the self-consistent procedure.