The extended multireference quasi-degenerate perturbation theory, proposed by Granovsky [J. Chem. Phys. 134, 214113 (2011)], is combined with internally contracted multi-state complete active space second-order perturbation theory (XMS-CASPT2). The first-order wavefunction is expanded in terms of the union of internally contracted basis functions generated from all the reference functions, which guarantees invariance of the theory with respect to unitary rotations of the reference functions. The method yields improved potentials in the vicinity of avoided crossings and conical intersections. The theory for computing nuclear energy gradients for MS-CASPT2 and XMS-CASPT2 is also presented and the first implementation of these gradient methods is reported. A number of illustrative applications of the new methods are presented.

On behalf of the development team, I review the capabilities of the Brilliantly Advanced General Electronic-structure Library (BAGEL) program package in this article. BAGEL is a newly developed full-fledged program package for electronic-structure computation in quantum chemistry, which is released under the GNU General Public License with many contributions from the developers. The unique features include analytical CASPT2 nuclear energy gradients and derivative couplings, relativistic multireference wave functions based on the Dirac equation, and implementations of novel electronic structure theories. All of the programs are efficiently parallelized using both threads and MPI processes. We also discuss the code generator SMITH3, which has been used to implement some of the programs in BAGEL. The developers' contributions are listed at the end of the main text.

We report the development of the theory and computer program for analytical nuclear energy gradients for (extended) multi-state complete active space perturbation theory (CASPT2) with full internal contraction. The vertical shifts are also considered in this work. This is an extension of the fully internally contracted CASPT2 nuclear gradient program, recently developed for a state-specific variant by us [MacLeod and Shiozaki, J. Chem. Phys. 142, 051103 (2015)]; in this extension, the so-called λ equation is solved to account for the variation of the multi-state CASPT2 energies with respect to the change in the amplitudes obtained in the preceding statespecific CASPT2 calculations, and the Z-vector equations are modified accordingly. The program is parallelized using the MPI3 remote memory access protocol that allows us to perform efficient one-sided communication.The optimized geometries of the ground and excited states of a copper corrole and benzophenone are presented as numerical examples. The code is publicly available under the GNU General Public License.

An internally contracted multireference configuration interaction is developed which employs wave functions that explicitly depend on the electron-electron distance (MRCI-F12). This MRCI-F12 method has the same applicability as the MRCI method, while having much improved basis-set convergence with little extra computational cost. The F12b approximation is used to arrive at a computationally efficient implementation. The MRCI-F12 method is applied to the singlet-triplet separation of methylene, the dissociation energy of ozone, properties of diatomic molecules, and the reaction barrier and exothermicity of the F + H(2) reaction. These examples demonstrate that already with basis sets of moderate size the method provides near complete basis set MRCI accuracy, and hence quantitative agreement with the experimental data. As a side product, we have also implemented the explicitly correlated multireference averaged coupled pair functional method (MRACPF-F12).

The probability of non-radiative transitions in photochemical dynamics is determined by the derivative couplings, the couplings between different electronic states through the nuclear degrees of freedom. Efficient and accurate evaluation of the derivative couplings is, therefore, of central importance to realize reliable computer simulations of photochemical reactions. In this work, the derivative couplings for multistate multireference second-order perturbation theory (MS-CASPT2) and its 'extended' variant (XMS-CASPT2) are studied, in which we present an algorithm for their analytical evaluation. The computational costs for evaluating the derivative couplings are essentially the same as those for calculating the nuclear energy gradients. The geometries and energies calculated with XMS-CASPT2 for small molecules at minimum energy conical intersections (MECIs) are in good agreement with those computed by multireference configuration interaction. As numerical examples, MECIs are optimized using XMS-CASPT2 for stilbene and a GFP model chromophore (the 4-para-hydroxybenzylidene-1,2-dimethyl-imidazolin-5-one anion).

Analytical nuclear gradients for fully internally contracted complete active space second-order perturbation theory (CASPT2) are reported. This implementation has been realized by an automated code generator that can handle spin-free formulas for the CASPT2 energy and its derivatives with respect to variations of molecular orbitals and reference coefficients. The underlying complete active space self-consistent field and the so-called Z-vector equations are solved using density fitting. The implementation has been applied to the vertical and adiabatic ionization potentials of the porphin molecule to illustrate its capability.

We report the development of programs for on-the-fly surface-hopping dynamics simulations in the gas and condensed phases on the potential energy surfaces computed by multistate multireference perturbation theory (XMS-CASPT2) with full internal contraction. On-the-fly nonadiabatic dynamics simulations are made possible by improving the algorithm for XMS-CASPT2 nuclear energy gradient and derivative coupling evaluation. The program is interfaced to a surface-hopping dynamics program, Newton-X, and a classical molecular dynamics package, tinker, to realize such simulations. On-the-fly XMS-CASPT2 surface-hopping dynamics simulations of 9H-adenine and an anionic GFP model chromophore (para-hydroxybenzilideneimidazolin-5-one) in water are presented to demonstrate the applicability of our program to sizable systems. Our program is implemented in the bagel package, which is publicly available under the GNU General Public License.

We present an approach to accurately construct the few-state model Hamiltonians for singlet fission processes on the basis of an ab initio electronic structure method tailored to dimer wave functions, called an active space decomposition strategy. In this method, the electronic structure of molecular dimers is expressed in terms of a linear combination of products of monomer states. We apply this method to tetracene and pentacene, using monomer wave functions computed by the restricted active space (RAS) method. Near-exact wave functions are computed for π-electrons of dimers that contain up to 7 × 1012 electronic configurations. Our product ansatz preserves the diabatic picture of the minimal dimer model, allowing us to accurately identify model Hamiltonians. The wave functions obtained from the model Hamiltonians account for more than 99% of the total wave functions. The resulting model Hamiltonians are shown to be converged with respect to all the parameters in the model, and corroborate previously reported coupling strengths.