Self-consistent Field Wave Functions Research Articles

Geometric Optimization of Restricted-Open and Complete Active Space Self-Consistent Field Wave Functions.

We explore Riemannian optimization methods for Restricted-Open-shell Hartree-Fock (ROHF) and Complete Active Space Self-Consistent Field (CASSCF) methods. After showing that ROHF and CASSCF can be reformulated as optimization problems on so-called "flag manifolds", we review Riemannian optimization basics and their application to these specific problems. We compare these methods to traditional ones and find robust convergence properties without fine-tuning of numerical parameters. Our study suggests that Riemannian optimization is a valuable addition to orbital optimization for ROHF and CASSCF, warranting further investigation.

Restricted Open-Shell Hartree-Fock Method for a General Configuration State Function Featuring Arbitrarily Complex Spin-Couplings.

In this work, we present a general spin restricted open-shell Hartree-Fock (ROHF) implementation that is able to generate self-consistent field (SCF) wave functions for an arbitrary configuration state function (CSF). These CSFs can contain an arbitrary number of unpaired electrons in arbitrary spin-couplings. The resulting method is named CSF-ROHF. We demonstrate that starting from the ROHF energy expression, for example, the one given by Edwards and Zerner, it is possible to obtain the values of the ROHF vector-coupling coefficients by setting up an open-shell for each group of consecutive parallel-coupled spins dictated by the unique spin-coupling pattern of any given CSF. To achieve this important and nontrivial goal, we employ the machinery of the iterative configuration expansion configuration interaction (ICE-CI) method, which is able to tackle general CI problems on the basis of spin-adapted CSFs. This development allows for the efficient generation of SCF spin-eigenfunctions for systems with complex spin-coupling patterns, such as polymetallic chains and metal clusters, while maintaining SCF scaling with system size (quadratic or less, depending on the specific algorithm and approximations chosen).

Potential-Energy Curves for the Ground and Several Electronic States of NdO and NdS.

Potential energy curves (PECs) were calculated for 21 and 18 electronic states of NdO and NdS molecules, respectively. In each case, static electron correlation effects were described by incomplete model space multiconfiguration self-consistent field wave functions based on an active space that included the most important valence orbitals. Dynamic electron correlation was included by the multireference second-order generalized Van Vleck perturbation theory method. Scalar-relativistic contributions were included by the effective core potential approach, using def2-TZVPP basis sets. Spin-dependent relativistic corrections were determined to be small and negligible for the Nd atom and so were not included in the calculations. The 21 and 18 electronic states of NdO and NdS were predicted to be in the excitation energy range of ∼3.2 and ∼2.7 eV, respectively. The ground electronic states of NdO and NdS were determined as 15H (6s4fσ4fϕ4fδ) and 15H (4fϕ4fπ4fπ6s), with spectroscopic constants: bond length Re = 1.780 and 2.325 Å, and harmonic frequency ωe = 891 and 538 cm-1, respectively.

Theoretical Prediction of Core-Level Binding Energies: Analysis of Unexpected Errors.

The analysis of the C(1s) and O(1s) core-level binding energies (CLBEs) of selected molecules computed by means of total energy Hartree-Fock (ΔSCF-HF) differences shows that in some cases, the calculated values for the C(1s) are larger than the experiment, which is unexpected. The origin of these unexpected errors of the Hartree-Fock ΔSCF BEs is shown to arise from static, nondynamical, electron correlation effects which are larger for the ion than for the neutral system. Once these static correlation effects are included by using complete active space self-consistent field (CASSCF) wave functions that include internal correlation terms, the resulting ΔSCF BEs are, as expected, smaller than measured values.

Spinless formulation of linearized adiabatic connection approximation and its comparison with the second order N-electron valence state perturbation theory.

The adiabatic connection (AC) approximation, along with its linearized variant AC0, was introduced as a method of obtaining dynamic correlation energy. When using a complete active space self-consistent field (CASSCF) wave function as a reference, the AC0 approximation is considered one of the most efficient multi-reference perturbation theories. It only involves the use of 1st- and 2nd-order reduced density matrices. However, some numerical results have indicated that the excitation energies predicted by AC0 are not as reliable as those from the second-order N-electron valence state perturbation theory (NEVPT2). In this study, we develop a spinless formulation of AC0 based on the Dyall Hamiltonian and provide a detailed comparison between AC0 and NEVPT2 approaches. We demonstrate the components within the correlation energy expressions that are common to both methods and those unique to either AC0 or NEVPT2. We investigate the role of the terms exclusive to NEVPT2 and explore the possibility of enhancing AC0's performance in this regard.

Complex Linear Response Functions for a Multiconfigurational Self-Consistent Field Wave Function in a High Performance Computing Environment.

We present novel developments for the highly efficient evaluation of complex linear response functions of a multiconfigurational self-consistent field (MCSCF) wave function as implemented in MultiPsi. Specifically, expressions for the direct evaluation of linear response properties at given frequencies using the complex polarization propagator (CPP) approach have been implemented, within both the Tamm-Dancoff approximation (TDA) and the random phase approximation (RPA). Purely real algebra with symmetric and antisymmetric trial vectors in a shared subspace is used wherein the linear response equations are solved. Two bottlenecks of large scale MC-CPP calculations, namely, the memory footprint and computational time, are addressed. The former is addressed by limiting the size of the subspace of trial vectors by using singular value decomposition (SVD) on either orbital or CI subspaces. The latter is addressed using an efficient parallel implementation as well as the strategy of dynamically adding linear response equations at near-convergence to neighboring roots. Furthermore, a novel methodology for decomposing MC-CPP spectra in terms of intuitive orbital excitations in an approximate fashion is presented. The performance of the code is illustrated with several numerical examples, including the X-ray spectrum of a molecule with nearly one hundred atoms. Additionally, for X-ray spectroscopy, the effect of including or excluding the core orbital in the active space on small covalent metal complexes is discussed.

Excited-State Aromaticity Reversals in Naphthalene and Anthracene.

Aromaticity reversals between the electronic ground (S0) and low-lying singlet (S1, S2) and triplet (T1, T2, T3) states of naphthalene and anthracene are investigated by calculating the respective off-nucleus isotropic magnetic shielding distributions using complete-active-space self-consistent field (CASSCF) wavefunctions involving gauge-including atomic orbitals (GIAOs). The shielding distributions around the aromatic S0, antiaromatic S1 (1Lb), and aromatic S2 (1La) states in naphthalene are found to resemble the outcomes of fusing together the respective S0, S1, and S2 shielding distributions of two benzene rings. In anthracene, 1La is lower in energy than 1Lb, and as a result, the S1 state becomes aromatic, and the S2 state becomes antiaromatic; the corresponding shielding distributions are found to resemble extensions by one ring of those around the S2 and S1 states in naphthalene. The lowest antiaromatic singlet state of either molecule is found to be significantly more antiaromatic than the respective T1 state, which shows that it would be incorrect to assume that the similarity between the (anti)aromaticities of the S1 and T1 states in benzene, cyclobutadiene, and cyclooctatetraene would be maintained in polycyclic aromatic hydrocarbons.

How Different are the Diamagnetic and Paramagnetic Contributions to Off-Nucleus Shielding in Aromatic and Antiaromatic Rings?

The spatial variations in the diamagnetic and paramagnetic contributions to the off-nucleus isotropic shielding, , and to the zz component of the off-nucleus shielding tensor, , around benzene (C6 H6 ) and cyclobutadiene (C4 H4 ) are investigated using complete-active-space self-consistent field wavefunctions. Despite the substantial differences between and around the aromatic C6 H6 and the antiaromatic C4 H4 , the diamagnetic and paramagnetic contributions to these quantities, and , and and , are found to behave similarly in the two molecules, shielding and deshielding, respectively, each ring and its surroundings. The different signs of the most popular aromaticity criterion, the nucleus-independent chemical shift (NICS), in C6 H6 and C4 H4 are shown to follow from a change in the balance between the respective diamagnetic and paramagnetic contributions. Thus, the different NICS values for antiaromatic and antiaromatic molecules cannot be attributed to differences in the ease of access to excited states only; differences in the electron density, which determines the overall bonding picture, also play an important role.

Near-exact nuclear gradients of complete active space self-consistent field wave functions.

In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules. We show that HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant errors. The ability of the HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems, which require a larger active space than possible with a complete active space self-consistent field. Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT). We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional. We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.

DENSITY FUNCTIONAL AND GREEN’S FUNCTIONS METHOD TO COMPUTING SPECTRAL PARAMETERS OF DIATOMIC MOLECULES

It is presented a new effective method of calculating the energy and spectral parameters of diatomic molecules based on the hybrid theory of the quasi-particle density functional and the theory of Green's functions. As an illustration, the data of the calculation of the vertical ionization potentials and the coupling constants (vibrational structure) of the photoelectron spectra of a number of diatomic molecules, in particular, N2, are received. It is presented a detailed comparison of the received results with the data of standard theories of the Hartree-Fock type, the multi-configuration electron propagator method, and the extended theory based on Koopmans' theorem using multi-configurational self-consistent field wave functions with different sets of basis functions. Iit is shown that the consistent, maximally precise consideration of exchange-correlation effects and reorganization effects within the framework of the combined theory leads to a rather significant improvement in the agreement of theoretical and experimental data both in terms of ionization potentials and photoelectron spectra in general.

Adiabatic state preparation of correlated wave functions with nonlinear scheduling functions and broken-symmetry wave functions

Adiabatic state preparation (ASP) can generate the correlated wave function by simulating the time evolution of wave function under the time-dependent Hamiltonian that interpolates the Fock operator and the full electronic Hamiltonian. However, ASP is inherently unsuitable for studying strongly correlated systems, and furthermore practical computational conditions for ASP are unknown. In quest for the suitable computational conditions for practical applications of ASP, we performed numerical simulations of ASP in the potential energy curves of N2, BeH2, and in the C2v quasi-reaction pathway of the Be atom insertion to the H2 molecule, examining the effect of nonlinear scheduling functions and the ASP with broken-symmetry wave functions with the S2 operator as the penalty term, contributing to practical applications of quantum computing to quantum chemistry. Eventually, computational guidelines to generate the correlated wave functions having the square overlap with the complete-active space self-consistent field wave function close to unity are discussed.

Ab initio composite strategies and multireference approaches for lanthanide sulfides and selenides.

The f-block ab initio correlation consistent composite approach was used to predict the dissociation energies of lanthanide sulfides and selenides. Geometry optimizations were carried out using density functional theory and coupled cluster singles, doubles, and perturbative triples with one- and two-component Hamiltonians. For the two-component calculations, relativistic effects were accounted for by utilizing a third-order Douglas-Kroll-Hess Hamiltonian. Spin-orbit coupling was addressed with the Breit-Pauli Hamiltonian within a multireference configuration interaction approach. The state averaged complete active space self-consistent field wavefunctions obtained for the spin-orbit coupling energies were used to assign the ground states of diatomics, and several diagnostics were used to ascertain the multireference character of the molecules.

Wavefunction-based quantum-chemical ab initio calculations for core electron binding energies of small open shell molecules

Core electron binding energies (CEBEs), i.e. ionization energies of 1s core orbitals, are calculated by means of wavefunction-based quantum-chemical ab initio methods for a series of small open-shell molecules containing first-row atoms. The calculations are performed in three steps: (a) Koopmans’ theorem, where the orbitals of the electronic ground state are used unchanged also for the ions, (b) Hartree–Fock or self consistent field (SCF) approximation in which the orbitals are allowed to relax after 1s ionization (ΔSCF), (c) dynamic correlation effects on top of SCF. For open-shell molecules 1s ionization leads to ions in several spin states, mostly to a pair of a triplet and a singlet state. In several cases one or both of these ionic states are only poorly described by a single-reference SCF wavefunction, therefore a multi-reference complete active space self consistent field (CAS-SCF) wavefunction is used instead. The correlation effects are evaluated by means of our multi-reference coupled electron pair approximation program. The accuracy of the calculated CEBEs is in the order of 0.1–0.4 eV. This is in agreement with experimental results for NO and O2. But there exist only very few gas phase data for CEBEs of open-shell molecules.

A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions.

In this work, we present a one-step second-order converger for state-specific (SS) and state-averaged (SA) complete active space self-consistent field (CASSCF) wave functions. Robust convergence is achieved through step restrictions using a trust-region augmented Hessian (TRAH) algorithm. To avoid numerical instabilities, an exponential parameterization of variational configuration parameters is employed, which works with a nonredundant orthogonal complement basis. This is a common approach for SS-CASSCF and is extended to SA-CASSCF wave functions in this work. Our implementation is integral direct and based on intermediates that are formulated in either the sparse atomic-orbital or small active molecular-orbital basis. Thus, it benefits from a combination with efficient integral decomposition techniques, such as the resolution-of-the-identity or the chain-of-spheres for exchange approximations. This facilitates calculations on large molecules, such as a Ni(II) complex with 231 atoms and 5154 basis functions. The runtime performance of TRAH-CASSCF is competitive with the other state-of-the-art implementations of approximate and full second-order algorithms. In comparison with a sophisticated first-order converger, TRAH-CASSCF calculations usually take more iterations to reach convergence and, thus, have longer runtimes. However, TRAH-CASSCF calculations still converge reliably to a true minimum even if the first-order algorithm fails.

DQC: A Python program package for differentiable quantum chemistry.

Automatic differentiation represents a paradigm shift in scientific programming, where evaluating both functions and their derivatives is required for most applications. By removing the need to explicitly derive expressions for gradients, development times can be shortened and calculations can be simplified. For these reasons, automatic differentiation has fueled the rapid growth of a variety of sophisticated machine learning techniques over the past decade, but is now also increasingly showing its value to support ab initio simulations of quantum systems and enhance computational quantum chemistry. Here, we present an open-source differentiable quantum chemistry simulation code and explore applications facilitated by automatic differentiation: (1) calculating molecular perturbation properties, (2) reoptimizing a basis set for hydrocarbons, (3) checking the stability of self-consistent field wave functions, and (4) predicting molecular properties via alchemical perturbations.

Localized Spin Rotations: A Size-Consistent Approach to Nonorthogonal Configuration Interaction.

Current nonorthogonal configuration interaction (NOCI) methods often use a set of self-consistent field (SCF) states selected based on chemical intuition. However, it may be challenging to track these SCF states across a dissociation profile and the NOCI states recovered may be spin contaminated. In this Article, we propose a method of applying spin rotation on symmetry broken unrestricted Hartree-Fock (sb-UHF) states to generate a basis for NOCI. The dissociation of ethene was examined by localizing spin rotation on each resulting carbene fragment. We show that this gives a size-consistent description of its dissociation and results in spin-pure states at all geometries. The dissociation was also studied with different orbitals, namely, canonical UHF and absolutely localized molecular orbitals (ALMO). Furthermore, we demonstrate that the method can be used to restore spin symmetry of symmetry broken SCF wave functions for molecules of various sizes, marking an improvement over existing NOCI methods.

Electronic structure of the ground and excited states of neutral and charged silicon hydrides, SiHx0/+/−, x = 1–4

The ground and excited electronic states of the titled species are investigated with multi-reference configuration interaction and diffuse basis sets. We found that in addition to the valence orbitals, the inclusion of the 4s, 4p, and especially 3d orbitals (although with minimal population) of silicon in the active space of the reference complete active space self-consistent field wavefunction are necessary for the proper convergence of the calculations. We also demonstrate that the aug-cc-pVTZ basis set provides quite accurate results compared to both larger basis sets and basis set limit results at a lower computational cost. The excited states involve excitations within the 3s and 3p orbitals of silicon (especially for the mono- and di-hydrides), followed by excitations from the Si-H bonding orbitals to either silicon valence or Rydberg (4s, 4p) orbitals. The number of electronic states per energy unit decrease as we add hydrogen atoms, and the first excited state of SiH4 is at 9.0 eV and leads to SiH3 + H. All species have stable ground state structures with all hydrogen atoms bound to silicon, except for SiH4+ and SiH4-. The former dissociates to SiH2+ + H2, while the latter loses an electron or can dissociate forming H2 as well.

Development of anharmonic vibrational structure theory using backflow transformation

A backflow (BF) transformation was applied to a vibrational self-consistent field wavefunction, which is the mean-field approximation in an anharmonic vibrational structure theory. We tested its effectiveness on the H2O molecule using quantum Monte Carlo method. Each normal mode coordinate was transformed into collective coordinates using the correlation function depending on whole coordinates. The vibrational frequencies with the BF transformation reasonably agree with the vibrational configuration interaction (VCI) results. The vibrational probability density to which the BF transformation was applied was also consistent with the VCI results and showed that the BF transformation appropriately take into account the inter-mode couplings.

Spin-Orbit Natural Transition Orbitals and Spin-Forbidden Transitions.

Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. The present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be nonzero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T2-S1 state mixing in pyrazine, spin-forbidden versus spin-allowed 4f-5d transitions in the Tb3+ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy)3].

Covalency of Trivalent Actinide Ions with Different Donor Ligands: Do Density Functional and Multiconfigurational Wavefunction Calculations Corroborate the Observed "Breaks"?

A comprehensive ab initio study of periodic actinide-ligand bonding trends for trivalent actinides is performed. Relativistic density functional theory (DFT) and complete active-space (CAS) self-consistent field wavefunction calculations are used to dissect the chemical bonding in the [AnCl6]3-, [An(CN)6]3-, [An(NCS)6]3-, [An(S2PMe2)3], [An(DPA)3]3-, and [An(HOPO)]- series of actinide (An = U-Es) complexes. Except for some differences for the early actinide complexes with DPA, bond orders and excess 5f-shell populations from donation bonding show qualitatively similar trends in 5f n active-space CAS vs DFT calculations. The influence of spin-orbit coupling on donation bonding is small for the tested systems. Along the actinide series, chemically soft vs chemically harder ligands exhibit clear differences in bonding trends. There are pronounced changes in the 5f populations when moving from Pu to Am or Cm, which correlate with previously noted "breaks" in chemical trends. Bonding involving 5f becomes very weak beyond Cm/Bk. We propose that Cm(III) is a borderline case among the trivalent actinides that can be meaningfully considered to be involved in ground-state 5f covalent bonding.