Research

Research

We develop new methods and algorithms for electronic structure theory. Click the links below to find out more!

Two-electron reduced-density-matrix-based methods

The standard approach for capturing nondynamical correlation effects in quantum chemical computations is the complete active space self-consistent field (CASSCF) method. In most implementations of CASSCF, the electronic structure of the active space is described by a configuration interaction (CI) wave function. Because the size of the CI wave function grows exponentially with the size of the active space, the application of CASSCF to large actives spaces is nontrivial. The treatment of large active spaces requires one to either (i) abandon CI in favor of some other wave function expansion that scales polynomially with system size or (ii) abandon the wave function altogether. Methods that employ the two-electron reduced-density matrix (2-RDM) as the central variable, such as the variational 2-RDM method (v2RDM), facilitate the design of polynomially-scaling CASSCF. The challenge in working directly with the 2-RDM is in ensuring that that an optimized 2-RDM actually corresponds to an antisymmetric

We have developed a free and open-source implementation of a v2RDM-driven CASSCF method as a plugin to the Psi4 electronic structure package. The active-space 2-RDM is obtained from a semidefinite optimization procedure, without the use of the N-electron wave function. Our CASSCF implementation is applicable to systems with active spaces as large or larger than 50 electrons in 50 orbitals [a (50,50) active space] and thousands of external orbitals (see Fig. 1).

Figure 1. The singlet/triplet state splitting at the v2RDM CASSCF/cc-pVXZ (X=D,T) level of theory. Experimental and DMRG values are taken from J. Hachmann, J. J. Dorando, M. Aviles, and G. K.-L. Chan, J. Chem. Phys.

You can read more about v2RDM methods and our implementation here:

The direct optimization of the 2-RDM in the v2RDM approach is a semidefinite programming (SDP) problem, and the SDP solver we employ is built upon an augmented Lagrangian formalism. Because the Lagrangian is stationary with respect to variations in both the primal and dual solutions to the SDP problem, it can be shown that expressions for v2RDM-CASSCF analytic energy gradients are identical to those that arise in conventional, CI-based CASSCF methods. We have implemented analytic energy gradients for v2RDM-CASSCF which can be used for geometry optimizations and, in principle

Figure 2. Symmetry-unique carbon carbon bond lengths (Å) obtained at the v2RDM-CASSCF/cc-pVDZ level of theory (red) and derived from experiment [

You can read more about our derivation and implementation of analytic energy gradients for v2RDM-CASSCF here:

The ground-state

So, to model excited states in strongly-correlated systems, we employ and approach known as the extended random phase approximation (ERPA). In the simplest form of the ERPA, excited states are approximated as single excitations out of the ground state. In this case, one can determine excited-state energies and wave functions with knowledge of only ground-state 1- and 2-RDMs, the natural variables of the v2RDM method!

More details coming soon!

In the meantime, you can read more about pure-state

We have developed a method to extract excitation energies from the time-evolution of a one-electron reduced-density matrix (1-RDM) obtained from variational 2-RDM methods. The equation of motion (EOM) for the 1-RDM depends on the 2-RDM. In our method, the 2-RDM at a given time is obtained by a ground-state-like minimization of the energy with respect to the elements of the 2-RDM, subject to a subset of known

The method provides a good description of excitations to states of a different symmetry than the ground state. For H

The method also predicts accurate excitation energies to the first excited state of different symmetry than the ground state in 4-electron systems like the set of ions that are isoelectronic with beryllium and the molecule, LiH.

You can read more about time-dependent RDM methods here:

We are also developing a new implementation of the v2RDM-CASSCF approach within the Q-Chem electronic structure package. The Q-Chem implementation will boast a number of useful features, including (i) effective core potentials, (ii) analytic energy gradients within the density fitting approximation, and (iii) compatability graphical processing units for faster v2RDM-CASSCF computations. The v2RDM-CASSCF code will be available in the 5.0.2 release of Q-Chem.