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

Autor: Tsubasa Iino, Toru Shiozaki, Takeshi Yanai
Rok vydání: 2023
Předmět:
Zdroj: The Journal of Chemical Physics. 158:054107
ISSN: 1089-7690
0021-9606
DOI: 10.1063/5.0130636
Popis: 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 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 the whole set of constraints on the parameters of multi-root canonical matrix product states (MPSs) 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-active (AA) orbital rotation in the Lagrangian. This rotation, usually realized by the orbital localization, regulates the noninvariance of the DMRG-CASSCF wavefunction to the AA rotation. 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 the numerical gradients of SA-DMRG-CASSCF energies using a common number of renormalized basis. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized basis, which is required to attain adequate accuracy in the DMRG's total energies.
Databáze: OpenAIRE