Multireference Perturbation Theory Methods Research Articles

Ab Initio Calculation of UV-vis Absorption of Parent Mg, Fe, Co, Ni, Cu, and Zn Metalloporphyrins.

Relativistic restricted active space (RAS) second-order multireference perturbation theory (MRPT2) methods, incorporating spin-orbit (SO) coupling perturbatively via state interaction (SO-MRPT2/RASSCF), were used to reproduce the absorption spectra of parent metalloporphyrins containing the Mg2+, Zn2+, Co2+, Ni2+, Cu2+, or FeCl2+ ions in the 12,500-40,000 cm-1 region. Particular attention was paid to the interaction between the porphyrin ring and the metal 3d electrons in states of different multiplicities (we used metal 3d and double d-shell or 3d' orbitals). For this class of compounds, the N-electron valence state perturbation theory (NEVPT2) method is superior to the complete active space perturbation theory (CASPT2) and successfully reproduces the energies of all four characteristic transitions (Q, B, N, and L) of closed-shell metalloporphyrins. Inclusion of SO coupling was found to have very little effect on excitation energies and oscillator strengths. For FeCl2+ porphyrin, we treated ligand-to-metal charge-transfer (LMCT; π,d), metal ligand field (d,d), and metal-to-ligand charge-transfer (MLCT; d,π*) transitions within the same framework. The broad and intense spectral features associated with its B (Soret) band are attributed to multiconfigurational LMCT (d,π*) bands involving strong metal-ligand orbital mixing.

Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory.

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.

Molecular Geometry and Electronic Structure of Copper Corroles.

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree-Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S0) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T0), which is approximately 5 kcal/mol higher in energy compared to the S0 state for both the bare and substituted copper corroles. Notably, the planarity of the T0 geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S0 and T0 geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

Photodissociation dynamics and UV absorption spectrum of acetone oxide (CH3)2COO.

The photodynamics and B1A' ← X1A' absorption spectrum of acetone oxide, (CH3)2COO, are studied theoretically from first principles. The underlying adiabatic potential energy curves (and surfaces) are computed by a second-order multireference perturbation theory method and diabatized using a diabatization by ansatz scheme. To confirm the results, for selected geometries EOM-CCSD and XMS-RS2C calculations were also performed. The dynamical calculation rests on the multi-configuration time-dependent Hartree wavepacket propagation method. The experimental absorption spectrum is reproduced satisfactorily. This result serves to validate the Hamiltonian model built within the quasi-diabatic representation. Contrary to the smallest Criegee intermediate, CH2OO, it is found that the vibronic coupling between the B and C states of (CH3)2COO plays an essential role in reproducing the experimental absorption spectrum. Time-dependent electronic populations reveal a faster decay than for the smaller system CH2OO. This is interpreted in terms of the stronger coupling between the B and C states in the larger system leading to a shorter lifetime for the B state than in CH2OO.

Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking.

In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) approach from a computational quantum chemistry perspective and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQMC for general main-group chemistry applications. For systems where single determinant wave functions are qualitatively accurate, we expect the accuracy of ph-AFQMC in conjunction with a single-determinant trial wave function to be between that of coupled-cluster with singles and doubles (CCSD) and CCSD with perturbative triples (CCSD(T)). For these applications, ph-AFQMC should be a method of choice when canonical CCSD(T) is too expensive to run. For systems where multireference (MR) wave functions are needed for qualitative accuracy, ph-AFQMC is far more accurate than MR perturbation theory methods and competitive with MR configuration interaction (MRCI) methods. Due to the computational efficiency of ph-AFQMC compared to MRCI, we recommended ph-AFQMC as a method of choice for handling dynamic correlation in MR problems. We conclude with a discussion of important directions for future development of the ph-AFQMC approach.

ΔNO and the complexities of electron correlation in simple hydrogen clusters.

The Δ natural orbital (ΔNO) two-electron density matrix (2-RDM) and energy expression are derived from a multideterminantal wave function. The approximate ΔNO 2-RDM is combined with an on-top density functional and a double-counting correction to capture electron correlation. A trust-region Newton's method optimization algorithm for the simultaneous optimization of ΔNO orbitals and occupancies is introduced and compared to the previous iterative diagonalization algorithm. The combination of ΔNO and two different on-top density functionals, Colle-Salvetti (CS) and Opposite-spin exponential cusp and Fermi-hole correction (OF), is assessed on small hydrogen clusters and compared to density functional, single-reference coupled-cluster, and multireference perturbation theory (MRMP2) methods. The ΔNO-CS and ΔNO-OF methods outperform the single-reference methods and are comparable to MRMP2. However, there is a distinct qualitative error in the ΔNO potential energy surface for H4 compared to the exact. This discrepancy is explained through analysis of the ΔNO orbitals, occupancies, and the two-electron density.

Low-Lying Electronic States of the Nickel Dimer.

The generalized Van Vleck second order multireference perturbation theory (GVVPT2) method was used to investigate the low-lying electronic states of Ni2. Because the nickel atom has an excitation energy of only 0.025 eV to its first excited state (the least in the first row of transition elements), Ni2 has a particularly large number of low-lying states. Full potential energy curves (PECs) of more than a dozen low-lying electronic states of Ni2, resulting from the atomic combinations 3F4 + 3F4 and 3D3 + 3D3, were computed. In agreement with previous theoretical studies, we found the lowest lying states of Ni2 to correlate with the 3D3 + 3D3 dissociation limit, and the holes in the d-subshells were in the subspace of delta orbitals (i.e., the so-dubbed δδ-states). In particular, the ground state was determined as X 1Γg and had spectroscopic constants: bond length (R e) = 2.26 Å, harmonic frequency (ωe) = 276.0 cm−1, and binding energy (D e) = 1.75 eV; whereas the 1 1Σg + excited state (with spectroscopic constants: R e = 2.26 Å, ωe = 276.8 cm−1, and D e = 1.75) of the 3D3 + 3D3 dissociation channel lay at only 16.4 cm−1 (0.002 eV) above the ground state at the equilibrium geometry. Inclusion of scalar relativistic effects through the spin-free exact two component (sf-X2C) method reduced the bond lengths of both of these two states to 2.20 Å, and increased their binding energies to 1.95 eV and harmonic frequencies to 296.0 cm−1 for X 1Γg and 297.0 cm−1 for 1 1Σg +. These values are in good agreement with experimental values of R e = 2.1545 ± 0.0004 Å, ωe = 280 ± 20 cm−1, and D 0 = 2.042 ± 0.002 eV for the ground state. All states considered within the 3F4 + 3F4 dissociation channel proved to be energetically high-lying and van der Waals-like in nature. In contrast to most previous theoretical studies of Ni2, full PECs of all considered electronic states of the molecule were produced.

Multi-reference perturbation theory study on the RaCl molecule promising for the laser cooling

The potential energy curves of the ground and eight low-lying excited atomic and ionic terms of the RaCl molecule are calculated for the first time using the multi-reference perturbation theory method at the CASSCF/XMCQDPT2 level of theory including the spin-orbit coupling. The electronic term energies, equilibrium internuclear distances, dissociation energies, transition dipole moments, sequences of vibrational energies, harmonic vibrational frequencies, Franck–Condon factors, and radiative lifetimes are predicted. The data obtained in this study suggests the possibility of the direct laser cooling of the RaCl molecules.

Ab initio study on the excited states of pyrene and its derivatives using multi-reference perturbation theory methods

Low-lying singlet excited states of pyrene derivatives originated from the 1La and 1Lb states of pyrene have decisive influences on their absorption and fluorescence emission behaviors. Calculation of these excited states with quantitative accuracy is required for the theoretical design of pyrene derivatives tailored to target applications; this has been a long-standing challenge for ab initio quantum chemical calculations. In this study, we explore an adequate computational scheme through calculations of pyrene and its phenyl-substituted derivatives using multi-reference perturbation theory (MRPT) methods. All valence π orbitals on the pyrene moiety were assigned to the active orbitals. Computational load was reduced by restricting the electron excitations within the active orbitals in the preparation of reference configuration space. A generalized multi-configuration quasi-degenerate perturbation theory (GMCQDPT) was adopted to treat the reference space other than the complete active space. The calculated 1La and 1Lb excitation energies of pyrene are in good agreement with the experimental values. Calculations of 1,3,6,8-tetraphenyl pyrene suggest that the energetic ordering of 1La and 1Lb is inverted through tetraphenyl substitution and its lowest singlet excited state is the 1La parentage of pyrene, which is consistent with the experimentally deduced scheme. These results are not readily obtained by MRPT calculations with a limited number of active orbitals and single-reference theory calculations. Diphenyl pyrenes (DPPy) were also calculated at the same level of theory to investigate the dependence on the substitution positions of phenyl groups.

Analytic first-order derivatives of partially contracted n-electron valence state second-order perturbation theory (PC-NEVPT2)

A balanced treatment of dynamic and static electron correlation is important in computational chemistry, and multireference perturbation theory (MRPT) is able to do this at a reasonable computational cost. In this paper, analytic first-order derivatives, specifically gradients and dipole moments, are developed for a particular MRPT method, state-specific partially contracted n-electron valence state second-order perturbation theory (PC-NEVPT2). Only one linear equation needs to be solved for the derivative calculation if the Z-vector method is employed, which facilitates the practical application of this approach. A comparison of the calculated results with experimental geometrical parameters of O3 indicates excellent agreement although the calculated results for O3 - are slightly outside the experimental error bars. The 0-0 transition energies of various methylpyrimidines and trans-polyacetylene are calculated by performing geometry optimizations and seminumerical second-order geometrical derivative calculations. In particular, the deviations of 0-0 transition energies of trans-polyacetylene from experimental values are consistently less than 0.1 eV with PC-NEVPT2, indicating the reliability of the method. These results demonstrate the importance of adding dynamic electron correlation on top of methods dominated by static electron correlation and of developing analytic derivatives for highly accurate methods.

Quantum dynamical simulation of intramolecular singlet fission in covalently coupled pentacene dimers.

We analyze the dynamics of intramolecular singlet fission in a series of pentacene-based dimers consisting of two pentacene-like chromophores covalently bonded to a phenylene linker in ortho, meta, and para positions. The study uses a quantum dynamical approach that employs a model vibronic Hamiltonian whose parameters are obtained using multireference perturbation theory methods. The results highlight the different role of the direct and mediated mechanism in these systems, showing that the population of the multiexcitonic state, corresponding to the first step of the intramolecular singlet fission process, occurs mainly through a superexchange-like mechanism involving doubly excited or charge transfer states that participate in the process in a virtual way. In addition, the systems investigated provide insight into the roles that built-in geometrical constraints and the electronic structure of the spacer play in the intramolecular singlet fission process.

State-Interaction Pair-Density Functional Theory Can Accurately Describe a Spiro Mixed Valence Compound.

Mixed-valence compounds with strong couplings between electronic states constitute one of the most challenging types of multireference systems for electronic structure theory. Previous work on a model mixed-valence compound, the 2,2',6,6'-tetrahydro-4 H,4' H-5,5'-spirobi[cyclopenta[ c]pyrrole] cation, showed that multireference perturbation theory (MRPT) can give a physical energy surface for the mixed-valence compound only by going to the third order or by using a scheme involving averaging orbital energies in a way specific to mixed-valence systems. In this study, we show that second-order MRPT methods (CASPT2, MS-CASPT2, and XMS-CASPT2) can give good results by calculating the Fock operator for the zeroth-order Hamiltonian using the state-averaged density matrix. We also show that state-interaction pair-density functional theory (SI-PDFT) is free from the unphysical behavior of previously tested second-order MRPT methods for this prototype mixed-valence compound near the avoided crossing. This is very encouraging because of the much lower cost in applying SI-PDFT to large or complex systems.

Intramolecular Singlet Fission: Insights from Quantum Dynamical Simulations.

We investigate the dynamics of intramolecular singlet fission in a dimer consisting of two pentacene-based chromophores covalently bonded to a phenylene spacer using an approach that combines high-level ab initio multireference perturbation theory methods and quantum dynamical simulations. The results show that the population of the multiexcitonic state, corresponding to the first step of singlet fission, is facilitated by the existence of higher-lying doubly excited and charge transfer states that participate in a superexchange-like way. The important role played by high-frequency ring-breathing molecular vibrations in the process is also discussed.

An Isolated Molecule of Iron(II) Phthalocyanin Exhibits Quintet Ground-State: A Nexus between Theory and Experiment.

Iron(II) phthalocyanine (FePc) is an important member of the phthalocyanines family with potential applications in the fields of electrocatalysis, magnetic switching, electrochemical sensing, and phototheranostics. Despite the importance of electronic properties of FePc in these applications, a reliable determination of its ground-state is still challenging. Here we present combined state of the art computational methods and experimental approaches, that is, Mössbauer spectroscopy and Superconducting Quantum Interference Device (SQUID) magnetic measurements to identify the ground state of FePc. While the nature of the ground state obtained with density functional theory (DFT) depends on the functional, giving mostly the triplet state, multi-reference complete active space second-order perturbation theory (CASPT2) and density matrix renormalization group (DMRG) methods assign quintet as the FePc ground-state in gas-phase. This has been confirmed by the hyperfine parameters obtained from 57 Fe Mössbauer spectroscopy performed in frozen monochlorobenzene. The use of monochlorobenzene guarantees an isolated nature of the FePc as indicated by a zero Weiss temperature. The results open doors for exploring the ground state of other metal porphyrin molecules and their controlled spin transitions via external stimuli.

Theoretical study of intramolecular singlet fission in xanthene-bonded pentacene dimers

Singlet fission is a spin-allowed process by which two non-interacting triplet excitons are generated from a singlet exciton. This type of multiexciton generation mechanism may have the potential to overcome the Shockley-Queisser limit and could enable the development of more efficient solar cells. In this contribution, we analyze intramolecular singlet fission in two modified pentacene dimers covalently bonded to a xanthene spacer employing density functional theory and multireference perturbation theory methods. The systems investigated differ in the electronic structure and bulkiness of the substituents used and in this work we analyze the impact that these aspects have in the relative contributions of the direct and mediated mechanisms to the intramolecular singlet fission.

Initial excited-state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations.

The initial S1 excited-state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight-membered ring (locked-11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS-CASPT2, MS-CASPT2) as well as single-reference coupled-cluster CC2 method. The analysis of XMCQDPT2-based geometries reveals rather weak coupling between in-plane and out-of-plane structural evolution and minor energetical relaxation of three locked-11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out-of-plane deformation resulting in a very steep S1 energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground-state structures, the excited-state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4-0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked-11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.

Electronic Resonances of Nucleobases Using Stabilization Methods.

The interaction of low-energy electrons with nucleobases is important because of their potential to damage nucleic acids. In this work we investigate several low-lying resonances of nucleobases using quantum chemical methods including static and dynamical correlation coupled with orbital stabilization methods. Specifically, the equation of motion for electron affinities via couple cluster singles and doubles and multireference perturbation theory methods were used, and their performance was explored. Low-lying π* resonances were calculated and compared to previous theoretical and experimental results, showing good agreement for positions and widths. Feshbach/Core-excited resonances generated by attachment of an electron to a ππ* excited state of the bases were also calculated, providing for the first time accurate information for these resonances. Mixing between configurations corresponding to shape and Feshbach/core-excited resonances is present in all nucleobases, which complicates the theoretical treatment and necessitates multiconfigurational approaches for a proper description.

Photodissociation dynamics of the pyridinyl radical: Time-dependent quantum wave-packet calculations.

The H-atom photodissociationreaction from the pyridinyl radical (C5H5NH) via the low-lying πσ* excited electronic state is investigated by nonadiabatic time-dependent quantum wave-packet dynamics calculations. A model comprising three electronic states and three nuclear coordinates has been constructed using ab initio multi-configurational self-consistent-field and multi-reference perturbation theory methods. Two conical intersections among the three lowest electronic states have been characterized in the framework of the linear vibronic-coupling model. Time-dependent wave-packet simulations have been performed using the multi-configuration time-dependent Hartree method. The population dynamics of the diabatic and adiabatic electronic states and the time-dependent dissociation behavior are analyzed for various vibrational initial conditions. The results provide detailed mechanistic insight into the photoinduced H-atom dissociation process from a hypervalent aromatic radical and show that an efficient dissociationreaction through two conical intersections is possible.

Ab initio multi-reference perturbation theory calculations of the ground and some excited electronic states of the RbYb molecule

The potential energy curves of the ground and twelve low-lying excited electronic states of the RbYb molecule have been calculated using the multi-reference perturbation theory method at the CASSCF/XMCQDPT2 level of theory including the spin-orbit coupling. The electronic term energies, equilibrium internuclear distances, dissociation energies, transition dipole moments, the sequences of vibrational energies, the harmonic vibrational frequencies, and the Franck–Condon factors have been predicted. The potential energy curves at the asymptotic limits are in a good agreement with the experimental data for rubidium and ytterbium atoms. The obtained data would be useful for spectral experiments with RbYb molecules.

Ab initio multi-reference perturbation theory calculations of the ground and low-lying electronic states of the KRb molecule

The potential energy curves of the low-lying electronic states correlating up to the limit K(4p)+Rb(5s) of KRb molecule have been calculated using the multi-reference perturbation theory method at the CASSCF/XMCQDPT2 level of theory without and with spin–orbit coupling. The calculated parameters of the ground X1Σ+ state are in the best agreement among all previously performed ab initio calculations for the KRb molecule. The calculated vibrational intervals of the ground electronic term of the 39K85Rb molecule describe the experiment with the accuracy within ±1cm−1. The calculated intensities of the 21Σ+ (v′=3, J′=26)→X1Σ+ (v″=0…24, J″=25, 27) transitions satisfactory reflect the experimentally observed intensities distribution.