Jakub Chalupský

Complete active space second-order perturbation theory with cumulant approximation for extended active-space wavefunction from density matrix renormalization group

We report an extension of our previous development that incorporated quantum-chemical density matrix renormalization group (DMRG) into the complete active space second-order perturbation theory (CASPT2) [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. In the previous study, the combined theory, referred to as DMRG-CASPT2, was built upon the use of pseudo-canonical molecular orbitals (PCMOs) for one-electron basis. Within the PCMO basis, the construction of the four-particle reduced density matrix (4-RDM) using DMRG can be greatly facilitated because of simplicity in the multiplication of 4-RDM and diagonal Fock matrix in the CASPT2 equation. In this work, we develop an approach to use more suited orbital basis in DMRG-CASPT2 calculations, e.g., localized molecular orbitals, in order to extend the domain of applicability. Because the multiplication of 4-RDM and generalized Fock matrix is no longer simple in general orbitals, an approximation is made to it using the cumulant reconstruction neglecting higher-particle cumulants. Also, we present the details of the algorithm to compute 3-RDM of the DMRG wavefunction as an extension of the 2-RDM algorithm of Zgid et al. [J. Chem. Phys. 128, 144115 (2008)] and Chan et al. [J. Chem. Phys. 128, 144117 (2008)]. The performance of the extended DMRG-CASPT2 approach was examined for large-scale multireference systems, such as low-lying excited states of long-chain polyenes and isomerization potential of {[Cu(NH3)3]2O2}(2+).

Effect of spin-orbit coupling on reduction potentials of octahedral ruthenium(II/III) and osmium(II/III) complexes.

Reduction potentials of several M(2+/3+) (M = Ru, Os) octahedral complexes, namely, [M(H2O)6](2+/3+), [MCl6](4-/3-), [M(NH3)6](2+/3+), [M(en)3](2+/3+) [M(bipy)3](2+/3+), and [M(CN)6](4-/3-), were calculated using the CASSCF/CASPT2/CASSI and MRCI methods including spin-orbit coupling (SOC) by means of first-order quasi-degenerate perturbation theory. It was shown that the effect of SOC accounts for a systematic shift of approximately -70 mV in the reduction potentials of the studied ruthenium (II/III) complexes and an approximately -300 mV shift for the osmium(II/III) complexes. SOC splits the sixfold-degenerate (2)T(2g) ground electronic state (in ideal octahedral symmetry) of the M(3+) ions into the E((5/2)g) Kramers doublet and G((3/2)g) quartet, which were calculated to split by 1354-1573 cm(-1) in the Ru(3+) complexes and 4155-5061 cm(-1) in the Os(3+) complexes. It was demonstrated that this splitting represents the main contribution to the stabilization of the M(3+) ground state with respect to the closed-shell (1)A(1g) ground state in M(2+) systems. Moreover, it was shown that the accuracy of the calculated reduction potentials depends on the calculated solvation energies of both the oxidized and reduced forms. For smaller ligands, it involves explicit inclusion of the second solvation sphere into the calculations, whereas implicit solvation models yield results of sufficient accuracy for complexes with larger ligands. In such cases (e.g., [M(bipy)3](2+/3+) and its derivatives), very good agreement between the calculated (SOC-corrected) values of the reduction potentials and the available experimental values was obtained. These results led us to the conclusion that especially for Os(2+/3+) complexes, inclusion of SOC is necessary to avoid systematic errors of approximately 300 mV in the calculated reduction potentials.

Multireference Ab Initio Calculations of g tensors for Trinuclear Copper Clusters in Multicopper Oxidases.

EPR spectroscopy has proven to be an indispensable tool in elucidating the structure of metal sites in proteins. In recent years, experimental EPR data have been complemented by theoretical calculations, which have become a standard tool of many quantum chemical packages. However, there have only been a few attempts to calculate EPR g tensors for exchange-coupled systems with more than two spins. In this work, we present a quantum chemical study of structural, electronic, and magnetic properties of intermediates in the reaction cycle of multicopper oxidases and of their inorganic models. All these systems contain three copper(II) ions bridged by hydroxide or O(2-) anions and their ground states are antiferromagnetically coupled doublets. We demonstrate that only multireference methods, such as CASSCF/CASPT2 or MRCI can yield qualitatively correct results (compared to the experimental values) and consider the accuracy of the calculated EPR g tensors as the current benchmark of quantum chemical methods. By decomposing the calculated g tensors into terms arising from interactions of the ground state with the various excited states, the origin of the zero-field splitting is explained. The results of the study demonstrate that a truly quantitative prediction of the g tensors of exchange-coupled systems is a great challenge to contemporary theory. The predictions strongly depend on small energy differences that are difficult to predict with sufficient accuracy by any quantum chemical method that is applicable to systems of the size of our target systems.

Reactivity of the Binuclear Non-Heme Iron Active Site of Δ9 Desaturase Studied by Large-Scale Multireference Ab Initio Calculations

The results of density matrix renormalization group complete active space self-consistent field (DMRG-CASSCF) and second-order perturbation theory (DMRG-CASPT2) calculations are presented on various structural alternatives for the O-O and first C-H activating step of the catalytic cycle of the binuclear nonheme iron enzyme Δ(9) desaturase. This enzyme is capable of inserting a double bond into an alkyl chain by double hydrogen (H) atom abstraction using molecular O2. The reaction step studied here is presumably associated with the highest activation barrier along the full pathway; therefore, its quantitative assessment is of key importance to the understanding of the catalysis. The DMRG approach allows unprecedentedly large active spaces for the explicit correlation of electrons in the large part of the chemically important valence space, which is apparently conditio sine qua non for obtaining well-converged reaction energetics. The derived reaction mechanism involves protonation of the previously characterized 1,2-μ peroxy Fe(III)Fe(III) (P) intermediate to a 1,1-μ hydroperoxy species, which abstracts an H atom from the C10 site of the substrate. An Fe(IV)-oxo unit is generated concomitantly, supposedly capable of the second H atom abstraction from C9. In addition, several popular DFT functionals were compared to the computed DMRG-CASPT2 data. Notably, many of these show a preference for heterolytic C-H cleavage, erroneously predicting substrate hydroxylation. This study shows that, despite its limitations, DMRG-CASPT2 is a significant methodological advancement toward the accurate computational treatment of complex bioinorganic systems, such as those with the highly open-shell diiron active sites.

Mono- and binuclear non-heme iron chemistry from a theoretical perspective.

In this minireview, we provide an account of the current state-of-the-art developments in the area of mono- and binuclear non-heme enzymes (NHFe and NHFe2) and the smaller NHFe(2) synthetic models, mostly from a theoretical and computational perspective. The sheer complexity, and at the same time the beauty, of the NHFe(2) world represents a challenge for experimental as well as theoretical methods. We emphasize that the concerted progress on both theoretical and experimental side is a conditio sine qua non for future understanding, exploration and utilization of the NHFe(2) systems. After briefly discussing the current challenges and advances in the computational methodology, we review the recent spectroscopic and computational studies of NHFe(2) enzymatic and inorganic systems and highlight the correlations between various experimental data (spectroscopic, kinetic, thermodynamic, electrochemical) and computations. Throughout, we attempt to keep in mind the most fascinating and attractive phenomenon in the NHFe(2) chemistry, which is the fact that despite the strong oxidative power of many reactive intermediates, the NHFe(2) enzymes perform catalysis with high selectivity. We conclude with our personal viewpoint and hope that further developments in quantum chemistry and especially in the field of multireference wave function methods are needed to have a solid theoretical basis for the NHFe(2) studies, mostly by providing benchmarking and calibration of the computationally efficient and easy-to-use DFT methods.

Multistate Complete-Active-Space Second-Order Perturbation Theory Based on Density Matrix Renormalization Group Reference States

We present the development of the multistate multireference second-order perturbation theory (CASPT2) with multiroot references, which are described using the density matrix renormalization group (DMRG) method to handle a large active space. The multistate first-order wave functions are expanded into the internally contracted (IC) basis of the single-state single-reference (SS-SR) scheme, which is shown to be the most feasible variant to use DMRG references. The feasibility of the SS-SR scheme comes from two factors: first, it formally does not require the fourth-order transition reduced density matrix (TRDM) and second, the computational complexity scales linearly with the number of the reference states. The extended multistate (XMS) treatment is further incorporated, giving suited treatment of the zeroth-order Hamiltonian despite the fact that the SS-SR based IC basis is not invariant with respect to the XMS rotation. In addition, the state-specific fourth-order reduced density matrix (RDM) is eliminated in an approximate fashion using the cumulant reconstruction formula, as also done in the previous state-specific DMRG-cu(4)-CASPT2 approach. The resultant method, referred to as DMRG-cu(4)-XMS-CASPT2, uses the RDMs and TRDMs of up to third-order provided by the DMRG calculation. The multistate potential energy curves of the photoisomerization of diarylethene derivatives with CAS(26e,24o) are presented to illustrate the applicability of our theoretical approach.

Influence of the choice of projection manifolds in the CASPT2 implementation

ABSTRACT The Complete Active Space Second-Order Perturbation Theory (CASPT2) is well-established as a high-accuracy electronic structure method. It was originally implemented in the early 1990s to an efficient computer code in the molcas program suite, and this implementation has been extensively used as a standard tool. Here, we report a comparison of it against two independent computer-aided implementations of the CASPT2 method, revealing that the CASPT2 energies provided by the original code of molcas (version 8 or earlier) are inconsistent with the predictions of the lately developed computer-aided implementations. It is shown that this error is associated with the projections of the first-order equation onto the fully internally contracted multireference bases which are partially inconsistent between the left- and right-hand sides. The degree of the errors is assessed by performing illustrative CASPT2 calculations. The errors in total CASPT2 energies are demonstrated to be negligible relative to chemical accuracy in many cases, while there is a difficult case where they may substantially alter chemical description. The incorporation of the consistent projections into molcas has been carried out, which is available in the version 8 sp1.

Molecular g-tensors from analytical response theory and quasi-degenerate perturbation theory in the framework of complete active space self-consistent field method

The molecular g-tensor is an important spectroscopic parameter provided by electron para magnetic resonance (EPR) measurement and often needs to be interpreted using computational methods. Here, we present two new implementations based on the first-order and second-order perturbation theories to calculate the g-tensors within the complete-active space self-consistent field (CASSCF) wave function model. In the first-order method, the quasi-degenerate perturbation theory (QDPT) is employed for constructing relativistic CASSCF states perturbed with the spin–orbit coupling operator, which is described effectively in one-electron form with the flexible nuclear screening spin–orbit approximation introduced recently by us. The second-order method is a newly reported approach built upon the linear response theory which accounts for the perturbation with respect to external magnetic field. It is implemented with the coupled–perturbed CASSCF (CP-CASSCF) approach, which provides an equivalent of untruncated sum-over-states expansion. The comparison of the performances between the first-order and second-order methods is shown for various molecules containing light to heavy elements, highlighting their relative strength and weakness. The formulations of QDPT and CP-CASSCF approaches as well as the derivation of the second-order Douglas–Kroll–Hess picture change of Zeeman operators are given in detail.