Multiconfiguration Self-consistent Field Method Research Articles

Fully Polarizable Multiconfigurational Self-Consistent Field/Fluctuating Charges Approach.

A multiscale model based on the coupling of the multiconfigurational self-consistent field (MCSCF) method and the classical atomistic polarizable fluctuating charges (FQ) force field is presented. The resulting MCSCF/FQ approach is validated by exploiting the CASSCF scheme through application to compute vertical excitation energies of formaldehyde and para-nitroaniline in aqueous solution. The procedure is integrated with molecular dynamics simulations to capture the solute's conformational changes and the dynamic aspects of solvation. Comparative analysis with alternative solvent models, gas-phase calculations, and experimental data provides insights into the model's accuracy in reproducing solute-solvent molecular interactions and spectral signals.

Semiclassical Nonadiabatic Molecular Dynamics Using Linearized Pair-Density Functional Theory.

Nonadiabatic molecular dynamics is an effective method for modeling nonradiative decay in electronically excited molecules. Its accuracy depends strongly on the quality of the potential energy surfaces, and its affordability for long direct-dynamic simulations with adequate ensemble averaging depends strongly on the cost of the required electronic structure calculations. Linearized pair-density functional theory (L-PDFT) is a recently developed post-self-consistent-field multireference method that can model potential energy surfaces with an accuracy similar to expensive multireference perturbation theories but at a computational cost similar to the underlying multiconfiguration self-consistent field method. Here, we integrate the SHARC dynamics and PySCF electronic structure code to utilize L-PDFT for electronically nonadiabatic calculations and use the combined programs to study the photoisomerization reaction of cis-azomethane. We show that L-PDFT is able to successfully simulate the photoisomerization without crashes, and it yields results similar to the more expensive extended multistate complete active space second-order perturbation theory. This shows that L-PDFT can model internal conversion, and it demonstrates its promise for broader photodynamics applications.

Linear Response Equations Revisited: A Simple and Efficient Iterative Algorithm.

We present an algorithm to solve the linear response equations for Hartree-Fock, Density Functional Theory, and the Multiconfigurational Self-Consistent Field method that is both simple and efficient. The algorithm makes use of the well-established symmetric and antisymmetric combinations of trial vectors but further orthogonalizes them with respect to the scalar product induced by the response matrix. This leads to a standard, symmetric block eigenvalue problem in the expansion subspace that can be solved by diagonalizing a symmetric, positive definite matrix half the size of the expansion space. Numerical tests showed that the algorithm is robust and stable.

Potential energy curves of molecular nitrogen up to N

The potential energy curves for molecular ions up to N are calculated in an ab initio manner using the multi configurational self-consistent field method. Specifically, we implement in an automatic way a previously used double loop optimisation scheme within the multi configurational self-consisted field method. We obtain the potential energy curves up to N ions with any combination of core, inner valence, and outer valence holes. Finally, we provide the code used to generate these potential energy curves.

Understanding the Excited-State Relaxation Mechanisms of Xanthophyll Lutein by Multi-configurational Electronic Structure Calculations.

The contradictory behaviors in light harvesting and non-photochemical quenching make xanthophyll lutein the most attractive functional molecule in photosynthesis. Despite several theoretical simulations on the spectral properties and excited-state dynamics, the atomic-level photophysical mechanisms need to be further studied and established, especially for an accurate description of geometric and electronic structures of conical intersections for the lowest several electronic states of lutein. In the present work, semiempirical OM2/MRCI and multi-configurational restricted active space self-consistent field methods were performed to optimize the minima and conical intersections in and between the 1Ag-, 2Ag-, 1Bu+, and 1Bu- states. Meanwhile, the relative energies were refined by MS-CASPT2(10,8)/6-31G*, which can reproduce correct electronic state properties as those in the spectroscopic experiments. Based on the above calculation results, we proposed a possible excited-state relaxation mechanism for lutein from its initially populated 1Bu+ state. Once excited to the optically bright 1Bu+ state, the system will propagate along the key reaction coordinate, i.e., the stretching vibration of the conjugated carbon chain. During this period of time, the 1Bu- state will participate in and forms a resonance state between the 1Bu- and 1Bu+ states. Later, the system will rapidly hop to the 2Ag- state via the 1Bu+/2Ag- conical intersection. Finally, the lutein molecule will survive in the 2Ag- state for a relatively long time before it internally converts to the ground state directly or via a twisted S1/S0 conical intersection. Notably, though the photophysical picture may be very different in solvents and proteins, the current theoretical study proposed a promising calculation protocol and also provided many valuable mechanistic insights for lutein and similar carotenoids.

Magnetic Behavior of Trigonal (Bi-)pyramidal 3d8 Mononuclear Nanomagnets: The Case of [Ni(MDABCO)2Cl3]ClO4.

This paper attempts to shed light on the origin of the magnetic behavior specific to trigonal bi- and pyramidal 3d8 mono- and polynuclear nanomagnets. The focus lies on entirely unraveling the system's intrinsic microscopic mechanisms and fundamental quantum mechanical relations governing the underlying electron dynamics. To this end, we develop a self-consistent approach to characterize, in great detail, all electron correlations and the ensuing fine structure of the energy spectra of a broad class of 3d8 systems. The mathematical framework is based on the multiconfigurational self-consistent field method and is devised to account for prospective quantum mechanical constraints that may confine the electron orbital dynamics while preserving the properties of all measurable quantities. We successfully characterize the experimentally observed magnetic anisotropy properties of a slightly distorted trigonal bipyramidal Ni2+ coordination complex, demonstrating that such compounds do not exhibit intrinsic huge zero-field splitting and inherent giant magnetic anisotropy. We reproduce qualitatively and quantitatively the behavior of the low-field magnetic susceptibility, magnetization, low-, and high-field electron paramagnetic resonance spectroscopy measurements and provide an in-depth analysis of the obtained results.

Electron and ion spectroscopy of azobenzene in the valence and core shells.

Azobenzene is a prototype and a building block of a class of molecules of extreme technological interest as molecular photo-switches. We present a joint experimental and theoretical study of its response to irradiation with light across the UV to x-ray spectrum. The study of valence and inner shell photo-ionization and excitation processes combined with measurement of valence photoelectron-photoion coincidence and mass spectra across the core thresholds provides a detailed insight into the site- and state-selected photo-induced processes. Photo-ionization and excitation measurements are interpreted via the multi-configurational restricted active space self-consistent field method corrected by second order perturbation theory. Using static modeling, we demonstrate that the carbon and nitrogen K edges of azobenzene are suitable candidates for exploring its photoinduced dynamics thanks to the transient signals appearing in background-free regions of the NEXAFS and XPS.

Intermolecular-Type Conical Intersections in Benzene Dimer.

The equilibrium and conical intersection geometries of the benzene dimer were computed in the framework of the conventional, linear-response time-dependent and spin-flipped time-dependent density functional theories (known as DFT, TDDFT and SF-TDDFT) as well as using the multiconfigurational complete active space self-consistent field (CASSCF) method considering the minimally augmented def2-TZVPP and the 6-31G(d,p) basis sets. It was found that the stacking distance between the benzene monomers decreases by about 0.5 Å in the first electronic excited state, due to the stronger intermolecular interaction energy, bringing the two monomers closer together. Intermolecular-type conical intersection (CI) geometries can be formed between the two benzene molecules, when (i) both monomer rings show planar deformation and (ii) weaker (approximately 1.6-1.8 Å long) C-C bonds are formed between the two monomers, with parallel and antiparallel orientation with respect to the monomer. These intermolecular-type CIs look energetically more favorable than dimeric CIs containing only one deformed monomer. The validity of the dimer-type CI geometries obtained by SF-TDDFT was confirmed by the CASSCF method. The nudged elastic band method used for finding the optimal relaxation path has confirmed both the accessibility of these intermolecular-type CIs and the possibility of the radiationless deactivation of the electronic excited states through these CI geometries. Although not as energetically favorable as the previous two CI geometries, there are other CI geometries characterized by the relative rotation of monomers at different angles around a vertical C-C axis.

Multiconfigurational SCF and Short-Range DFT Combined with Polarizable Density Embedding: Comparison of Linear-Response and State-Specific Solvatochromic Shifts of Acrolein and Para-nitrophenolate in Water.

The polarizable density embedding model is combined with the multiconfigurational self-consistent field and MC-srDFT electronic structure methods to calculate solvatochromic shifts of the n-π* absorption of acrolein and the π-π* absorption of the para-nitrophenolate anion in aqueous solution. Differences between linear-response (LR) and state-specific (SS) solvent shifts are analyzed by assessing the contributions of different terms in the solvent potential. This comparison shows that the differences are not only due to the intrinsically different response of LR and SS excitation energies to the polarizability of the environment but also due to a different response to the static part of the environment potential. These observations show that even in nonpolarizable environments, LR and SS calculations based on SCF (orbital optimization) methods do not necessarily agree on the spectral shift. The difference can be as large as, or even dominate, the difference due to dynamical polarization.

Ab initio approaches for [formula omitted] and [formula omitted] ions towards modeling of the [formula omitted] ion in cold helium plasma

As a first step towards realistic modeling of transport properties of the N2+ ion in helium gas, detailed investigations of the electronic structure of the ion as well as the N2+/He collision complex have been performed with the main focus on computational efficiency and accuracy. A broad range of correlation consistent basis sets and representative orbital spaces have been considered for both computationally cheap multi-configuration self-consistent field (MCSCF) methods and benchmark multi-reference configuration interaction (MRCI) approaches. It has been found that the computationally advantageous MCSCF approach, even if combined with basis sets of moderate sizes, leads to an acceptable accuracy and is thus suitable for direct dynamics simulations of N2+/He collisions to be addressed in subsequent works.

Analytical Energy Gradient for State-Averaged Orbital-Optimized Variational Quantum Eigensolvers and Its Application to a Photochemical Reaction.

Elucidating photochemical reactions is vital to understanding various biochemical phenomena and developing functional materials such as artificial photosynthesis and organic solar cells, albeit with notorious difficulty in both experiments and theories. The best theoretical way so far to analyze photochemical reactions at the level of ab initio electronic structure is the state-averaged multiconfigurational self-consistent field (SA-MCSCF) method. However, the exponential computational cost of classical computers with the increasing number of molecular orbitals hinders applications of SA-MCSCF for large systems we are interested in. Utilizing quantum computers was recently proposed as a promising approach to overcome such computational cost, dubbed as state-averaged orbital-optimized variational quantum eigensolver (SA-OO-VQE). Here, we extend a theory of SA-OO-VQE so that analytical gradients of energy can be evaluated by standard techniques that are feasible with near-term quantum computers. The analytical gradients, known only for the state-specific OO-VQE in previous studies, allow us to determine various characteristics of photochemical reactions such as the conical intersection (CI) points. We perform a proof-of-principle calculation of our methods by applying it to the photochemical cis-trans isomerization of 1,3,3,3-tetrafluoropropene. Numerical simulations of quantum circuits and measurements can correctly capture the photochemical reaction pathway of this model system, including the CI points. Our results illustrate the possibility of leveraging quantum computers for studying photochemical reactions.

Potential Energy Curves of Molecular Nitrogen for Singly and Doubly Ionized States with Core and Valence Holes.

Theoretical description of potential energy curves (PECs) of molecular ions is essential for interpretation and prediction of coupled electron-nuclear dynamics following ionization of parent molecule. However, an accurate representation of these PECs for core or inner valence ionized state is nontrivial, especially at stretched geometries for double- or triple-bonded systems. In this work, we report PECs of singly and doubly ionized states of molecular nitrogen using state-of-the-art quantum chemical methods. The valence, inner valence, and core ionized states have been computed. A double-loop optimization scheme that separates the treatment of the core and the valence orbitals during the orbital optimization step of the multiconfiguration self-consistent field method has been implemented. This technique allows the energy to be converged to any desired ionized state with any number of core or inner-shell holes. The present work also compares the PECs obtained using both delocalized and localized sets of orbitals for the core hole states. The PECs of a number of singly and doubly ionized valence states have also been computed and compared with previous studies. The computed PECs reported here are expected to be of importance for future studies to understand the interplay between photoionization and Auger spectra during the breakup of molecular nitrogen when interacting with intense free electron lasers.

Multiconfiguration Pair-Density Functional Theory for Transition Metal Silicide Bond Dissociation Energies, Bond Lengths, and State Orderings.

Transition metal silicides are promising materials for improved electronic devices, and this motivates achieving a better understanding of transition metal bonds to silicon. Here we model the ground and excited state bond dissociations of VSi, NbSi, and TaSi using a complete active space (CAS) wave function and a separated-pair (SP) wave function combined with two post-self-consistent field techniques: complete active space with perturbation theory at second order and multiconfiguration pair-density functional theory. The SP approximation is a multiconfiguration self-consistent field method with a selection of configurations based on generalized valence bond theory without the perfect pairing approximation. For both CAS and SP, the active-space composition corresponds to the nominal correlated-participating-orbital scheme. The ground state and low-lying excited states are explored to predict the state ordering for each molecule, and potential energy curves are calculated for the ground state to compare to experiment. The experimental bond dissociation energies of the three diatomic molecules are predicted with eight on-top pair-density functionals with a typical error of 0.2 eV for a CAS wave function and a typical error of 0.3 eV for the SP approximation. We also provide a survey of the accuracy achieved by the SP and extended separated-pair approximations for a broader set of 25 transition metal–ligand bond dissociation energies.

Auxiliary Atomic Relay Center Facilitates Enhanced Magnetic Couplings in Blatter's Radical.

The recent accomplishments in obtaining strong ferromagnetic exchange interactions in organic diradicals have made the field quite fascinating and even more promising toward its technological applications. In this context, herein, we report a unique combination of remarkably strong ferromagnetic exchange interactions coupled with molecular rigidity, utilizing superstable Blatter's radical as a spin source. The planar analogues of the parent Blatter's radical obtained by annulation with a chalcogen coupled to nitronyl nitroxide (NN) are investigated using density functional theory along with the wave function-based multiconfigurational self-consistent field methods, for example, complete active space self-consistent field (CASSCF)-N-electron valence state perturbation theory (NEVPT2). The calculations reveal phenomenal modulation in exchange couplings upon annulation such that remarkably strong ferromagnetic interactions are realized especially for a certain class of the Blatter-NN diradicals. The modulation of spin-spin interactions is rationalized by variation in spin density distribution and molecular torsional angles. We demonstrate that annulation in OMMs opens an additional coupling pathway via auxiliary X-atom acting as the atomic relay center which strongly manipulates the magnitude of exchange couplings.

Molecular magnetism in the multi-configurational self-consistent field method

We develop a structured theoretical framework used in our recent articles (2019 Eur. Phys. J. B 92 93 and 2020 Phys. Rev. B 101 094427) to characterize the unusual behavior of the magnetic spectrum, magnetization and magnetic susceptibility of the molecular magnet Ni4Mo12. The theoretical background is based on the molecular orbital theory in conjunction with the multi-configurational self-consistent field method and results in a post-Hartree–Fock scheme for constructing the corresponding energy spectrum. Furthermore, we construct a bilinear spin-like Hamiltonian involving discrete coupling parameters accounting for the relevant spectroscopic magnetic excitations, magnetization and magnetic susceptibility. The explicit expressions of the eigenenergies of the ensuing Hamiltonian are determined and the physical origin of broadening and splitting of experimentally observed peaks in the magnetic spectra is discussed. To demonstrate the efficiency of our method we compute the spectral properties of a spin-one magnetic dimer. The present approach may be applied to a variety of magnetic units based on transition metals and rare Earth elements.

Revisiting 12C16O[formula omitted]first negative system: An experimental and theoretical study

In this article, a new study of the 12C16O+First Negative system (B2Σ+→X2Σ+), is performed both theoretically and experimentally. From an experimental point of view, high-resolution Fourier transform spectroscopy in the ultraviolet region was used. On the theoretical side, ab initioelectronic energies were calculated employing the state-averaged multiconfiguration self-consistent field (MCSCF) method, followed by the icMRCI/AV6Z level of theory with Davidson modification. Through these two analyses, theoretical and experimental, values for the spectroscopic constants, potential energy curves, and parameters concerning the transition are presented. The results, both theoretical and experimental, obtained in this work are compared with values obtained in previous works, pointing out certain discrepancies with more recent works and providing new and improved values for the band origins of this transition.

An Electron Propagator Approach Based on a Multiconfigurational Reference State for the Investigation of Negative-Ion Resonances Using a Complex Absorbing Potential Method.

Negative-ion resonances are important metastable states that result from the collision between an electron and a neutral target. The course of many chemical processes in nature is often dictated by how an intermediate resonance state falls apart. This article reports on the development of an electron propagator (EP) based on a Hamiltonian Ĥ perturbed by a complex absorbing potential (CAP) and a multiconfigurational self-consistent field (MCSCF) initial state to study these resonances. Perturbation of Ĥ by a CAP makes the resonances amenable to a bound-state method like MCSCF. Resonances stand out among the non-resonant states as persistent complex eigenvalues of the perturbed Ĥ when the strength (η) of the CAP is varied. The MCSCF method gives a reliable and accurate description of the target states, especially when the non-dynamical correlations are dominant. The resonance energies are obtained from the poles of the EP. We propose three variants of our EP depending on how the effect of the CAP is introduced. We find that the computationally most efficient variant is the one in which the reference state of the EP is an unperturbed MCSCF wavefunction and a non-zero CAP is defined only on the virtual orbital subspace of the reference state. The onset of the CAP is carefully optimized in order to minimize the artifacts due to reflections from the CAP. An extrapolation method (based on a Padé approximant) and a de-perturbation method are adopted in order to account for the limitations of finite basis sets used and determine the resonance energy in the limit of η → 0. 2P Be-, 2Πg N2-, and 2Π CO- shape resonances are investigated. The position and width of these resonances computed in this study agree well with those reported earlier in the literature.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

The ground and first excited states of HoS studied by four-component relativistic KR-MCSCF and KRCI

The X18.5 ground state and the W17.5 first excited state of the HoS molecule are studied by the four-component Kramers-restricted multiconfiguration self-consistent field (KR-MCSCF) method, with an active space of (4f, 6s)11, and then by the Kramers-restricted configuration interaction method with 27 electrons in Ho(5s), (5p), (4f), (6s), S(3s), and S(3p) correlated. For the X18.5 state, Re, ωe, ΔG1/2, and D0 are calculated as 4.310 au, 477.385 cm−1, 477.72 cm−1, and 3.775 eV, respectively. The values of Re and ΔG1/2 agree well with the observed values of 4.301 au and 463.8811 cm−1, respectively. The dominant configuration state function (CSF) is |(4f5/2)6 (4f7/2,7/2)1 (4f7/2,5/2)1 (4f7/2,3/2)1 (4f7/2,1/2)1 (6s1/2,1/2)1|. This corresponds to the formal charge of Ho2+(5I8)S2−, representing ionic bonding where two electrons, one from [6s] and the other from [4f], are moved from Ho to S. In the doubly occupied S(3p) KR-MCSCF molecular spinors, hybridizations exist between S[3p] atomic spinors and Ho polarization functions [5d*] through which electrons are moved from S to Ho, resulting in covalent bonding. Consequently, ionic bonding is mixed with covalent bonding and 11 electrons of (4f, 6s) move outside the Ho13+([Xe])S2− ion core. The gross atomic charge of Ho is + 0.77 at the CI level. The excitation energy to the W17.5 state is calculated as 0.027 eV, close to the experimental value of 0.031 eV.

MCSCF optimization revisited. II. Combined first- and second-order orbital optimization for large molecules.

A new orbital optimization for the multiconfiguration self-consistent field method is presented. This method combines a second-order (SO) algorithm for the optimization of the active orbitals with the first-order super configuration interaction (SCI) optimization of the remaining closed-virtual rotations and is denoted as the SO-SCI method. The SO-SCI method significantly improves the convergence as compared to the conventional SCI method. In combination with density fitting, the intermediates from the gradient calculation can be reused to evaluate the two-electron integrals required for the active Hessian without introducing a large computational overhead. The orbitals and CI coefficients are optimized alternately, but the CI-orbital coupling is accounted for by the limited memory Broyden-Fletcher-Goldfarb-Shanno quasi-Newton method. This further improves the speed of convergence. The method is applicable to large molecules. The efficiency and robustness of the presented method is demonstrated in benchmark calculations for 21 aromatic molecules as well as for various transition metal complexes with up to 826 electrons and 5154 basis functions.