Analysis of wave functions for open-shell molecules

Development of molecular cluster models to probe pyrite surface reactivity

The recent discovery that anaerobic methanogens can reductively dissolve pyrite and utilize dissolution products as a source of iron and sulfur to meet their biosynthetic demands for these elements prompted the development of atomic-scale nanoparticle models, as maquettes of reactive surface sites, for describing the fundamental redox steps that take place at the mineral surface during reduction. The given report describes our computational approach for modeling n(FeS2 ) nanoparticles originated from mineral bulk structure. These maquettes contain a comprehensive set of coordinatively unsaturated Fe(II) sites that are connected via a range of persulfide (S2 2- ) ligation. In addition to the specific maquettes with n = 8, 18, and 32 FeS2 units, we established guidelines for obtaining low-energy structures by considering the pattern of ionic, covalent, and magnetic interactions among the metal and ligand sites. The developed models serve as computational nano-reactors that can be used to describe the reductive dissolution mechanism of pyrite to better understand the reactive sites on the mineral, where microbial extracellular electron-transfer reactions can occur.

Magnetic exchange and valence delocalization in a mixed valence [Fe2+Fe3+Te2]+ complex: insights from theory and interpretations of magnetic and spectroscopic data

A mixed valence binuclear Fe^2.5+-Fe^2.5+ (Robin-Day Class III) transition metal complex, [Fe^2.5+μTe₂Fe^2.5+]^1-, composed of two FeN₂Te₂ pseudo-tetrahedral units with μ-bridging Te^2- ligands was reported to exist in an unprecedented S = 3/2 ground state (Nature Chemistry, https://doi.org/10.1038/s41557-021-00853-5). For this and the homologous complexes containing Se^2- and S^2-, the Anderson-Hasegawa double exchange spin-Hamiltonian was broadly used to interpret the corresponding structural, spectroscopic and magnetic data. First principles multireference ab initio calculations are used here to simulate magnetic and spectroscopic EPR data; analysis of the results affords a rationale for the stabilization of the S = 3/2 ground state of the Fe₂ pair. Complete Active Space Self-Consistent Field (CASSCF) calculations and dynamical correlation accounted for by means of N-Electron Valence Perturbation Theory to Second Order (NEVPT2) reproduce well the g-factors determined from simulations of X-band EPR spectra. A crucial technical tool to achieve these results is: (i) use of a localized orbital formulation of the many-particle problem at the scalar-relativistic CASSCF step; (ii) choice of state averaging over states of a given spin (at the CASCI/NEVPT2 step); and (iii) accounting for spin-orbit coupling within the non-relativistic Born-Oppenheimer (BO) many-particle basis using Quasi-Degenerate Perturbation Theory (QDPT). The inclusion of the S = 5/2 spin manifold reproduced the observed increase in the magnetic susceptibility (χT) in the high temperature range (T > 100 K), which is explained by thermal population of the S = 5/2 excited state at energy 160 cm^-1 above the S = 3/2 ground state. Theoretical values of χT from experimentally reported data points in the temperature range from 3 to 30 K were further computed and analyzed using a model which takes spin-phonon coupling into account. The model considerations and the computational protocols of this study are generally applicable to any Class I/II mixed valence dimer. The work can potentially stimulate further experimental and theoretical work on bi- and oligonuclear transition metal complexes of importance to bioinorganic chemistry and life sciences.

Extraction of local spin-coupled states by second quantized operators

We present a methodology for analyzing chemical bonds embedded in the electronic wave function of molecules, especially in terms of spin correlations or so-called "local spin." In this paper, based on biorthogonal second quantization, the spin correlation functions of molecules are naturally introduced, which enables us to extract local singlet and local triplet elements from the wave function. We also clarify the relationship between these spin correlations and traditional chemical concepts, i.e., resonance structures. Several chemical reactions, including the intramolecular radical cyclization and the formation of preoxetane, are demonstrated to verify the analysis method numerically.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.

Local spin and open quantum systems: clarifying misconceptions, unifying approaches

The theory of open quantum systems (OQSs) is applied to partition the squared spin operator into fragment (local spin) and interfragment (spin-coupling) contributions in a molecular system.

The Right Answer for the Right Reason: My Personal Goal for Quantum Chemistry

A brief history of quantum theory is given to illustrate the barriers to progress caused by preconceived ideas. The biases in my own thinking which I had to overcome to approach the right answer for the right reason are discussed. This is followed by a personal autobiography illustrating how I have led a life of serendipity with no real sense of purpose. Chance events have shaped my life. The algorithms for which I am best known are briefly discussed. Then highlights from the many applications of theory to excited states, bonding in ice, spin properties and magnetism, (e,2e) shake-up spectra, and organic reactions are mentioned. This wide range of applications is mostly due to accidental collaboration with colleagues who sought my help. My real interest was in developing methods which could address these problems.

Selectivity of Cobalt Corrole for CO vs. O2 and N2 in Indoor Pollution

Coal combustion causes indoor pollution of CO. In this work, DFT calculations on cobalt corrole (Co(Cor)) with three most common indoor gas molecules (N₂, O₂ and CO) were performed. The Mulliken spin densities show that the ground states of Co(N₂)(Cor), Co(CO)(Cor) and Co(OC)(Cor) have an anti-ferromagnetic coupling fashion of the electrons on the Co 3d _z² orbital and the π orbital of the corrole ring. However, Co(O₂)(Cor) has a triplet ground state. With the spin contamination corrections, the Co(N₂)(Cor) binding energy was obtained at -50.6 kcal mol^-1 (B3LYP-D3). While CO can interact with Co(Cor) in two different ways, and their binding energies were -22.8 and -10.9 kcal mol^-1 (B3LYP-D3) for Co(CO)(Cor) and Co(OC)(Cor), respectively. The natural bond orbital charges on the axial ligands (NO, CO, OC) are increased upon the chemical bond formation. These are the cause of the shorten metal-ligand bond and the increase of the wavenumber of the metal-ligand bond vibrational transitions. While the charges for O₂ are decreased, leading to bond elongation as well as the decrease of the wavenumber upon complexation. Overall, O₂ was found to be hardly coordinated with Co(Cor). This study provides a detailed molecular understanding of interactions between a gas sensor and gaseous indoor air-pollutants.

Electronic Structures of Anti-Ferromagnetic Tetraradicals: Ab Initio and Semi-Empirical Studies

The energy relationships and electronic structures of the lowest-lying spin states in several anti-ferromagnetic tetraradical model systems are studied with high-level ab initio and semi-empirical methods. The Full-CI method (FCI), the complete active space second-order perturbation theory (CASPT2), and the n-electron valence state perturbation theory (NEVPT2) are employed to obtain reference results. By comparing the energy relationships predicted from the Heisenberg and Hubbard models with ab initio benchmarks, the accuracy of the widely used Heisenberg model for anti-ferromagnetic spin-coupling in low-spin polyradicals is cautiously tested in this work. It is found that the strength of electron correlation (|U/t|) concerning anti-ferromagnetically coupled radical centers could range widely from strong to moderate correlation regimes and could become another degree of freedom besides the spin multiplicity. Accordingly, the Heisenberg-type model works well in the regime of strong correlation, which reproduces well the energy relationships along with the wave functions of all the spin states. In moderately spin-correlated tetraradicals, the results of the prototype Heisenberg model deviate severely from those of multi-reference electron correlation ab initio methods, while the extended Heisenberg model, containing four-body terms, can introduce reasonable corrections and maintains its accuracy in this condition. In the weak correlation regime, both the prototype Heisenberg model and its extended forms containing higher-order correction terms will encounter difficulties. Meanwhile, the Hubbard model shows balanced accuracy from strong to weak correlation cases and can reproduce qualitatively correct electronic structures, which makes it more suitable for the study of anti-ferromagnetic coupling in polyradical systems.

Theory and practice of uncommon molecular electronic configurations

The electronic configuration of the molecule is the foundation of its structure and reactivity. The spin state is one of the key characteristics arising from the ordering of electrons within the molecule's set of orbitals. Organic molecules that have open-shell ground states and interesting physicochemical properties, particularly those influencing their spin alignment, are of immense interest within the up-and-coming field of molecular electronics. In this advanced review, we scrutinize various qualitative rules of orbital occupation and spin alignment, viz., the aufbau principle, Hund's multiplicity rule, and dynamic spin polarization concept, through the prism of quantum mechanics. While such rules hold in selected simple cases, in general the spin state of a system depends on a combination of electronic factors that include Coulomb and Pauli repulsion, nuclear attraction, kinetic energy, orbital relaxation, and static correlation. A number of fascinating chemical systems with spin states that fluctuate between triplet and open-shell singlet, and are responsive to irradiation, pH, and other external stimuli, are highlighted. In addition, we outline a range of organic molecules with intriguing non-aufbau orbital configurations. In such quasi-closed-shell systems, the singly occupied molecular orbital (SOMO) is energetically lower than one or more doubly occupied orbitals. As a result, the SOMO is not affected by electron attachment to or removal from the molecule, and the products of such redox processes are polyradicals. These peculiar species possess attractive conductive and magnetic properties, and a number of them that have already been developed into molecular electronics applications are highlighted in this review. WIREs Comput Mol Sci 2015, 5:440-459. doi: 10.1002/wcms.1233 For further resources related to this article, please visit the WIREs website.

Quantum entanglement in carbon–carbon, carbon–phosphorus and silicon–silicon bonds

We present a quantum entanglement analysis to dissect the bond orders in polyatomic molecules.

Potential energy curves and electronic structure of 3d transition metal hydrides and their cations

We investigate gas-phase neutral and cationic hydrides formed by 3d transition metals from Sc to Cu with density functional theory (DFT) methods. The performance of two exchange-correlation functionals, Boese-Martin for kinetics (BMK) and Tao-Perdew-Staroverov-Scuseria (TPSS), in predicting bond lengths and energetics, electronic structures, dipole moments, and ionization potentials is evaluated in comparison with available experimental data. To ensure a unique self-consistent field (SCF) solution, we use stability analysis, Fermi smearing, and continuity analysis of the potential energy curves. Broken-symmetry approach was adapted in order to get the qualitatively correct description of the bond dissociation. We found that on average BMK predicted values of dissociation energies and ionization potentials are closer to experiment than those obtained with high level wave function theory methods. This agreement deteriorates quickly when the fraction of the Hartree-Fock exchange in DFT functional is decreased. Natural bond orbital (NBO) population analysis was used to describe the details of chemical bonding in the systems studied. The multireference character in the wave function description of the hydrides is reproduced in broken-symmetry DFT description, as evidenced by NBO analysis. We also propose a new scheme to correct for spin contamination arising in broken-symmetry DFT approach. Unlike conventional schemes, our spin correction is introduced for each spin-polarized electron pair individually and therefore is expected to yield more accurate energy values. We derive an expression to extract the energy of the pure singlet state from the energy of the broken-symmetry DFT description of the low spin state and the energies of the high spin states (pentuplet and two spin-contaminated triplets in the case of two spin-polarized electron pairs). The high spin states are build with canonical natural orbitals and do not require SCF convergence.

Comparison of Semiempirical ZILSH and DFT Calculations of Exchange Constants in Fe4 Butterfly Complexes

Magnetic interactions in a series of tetranuclear Fe(3+) complexes with the butterfly core structure have been studied with semiempirical ZILSH and density functional theory (DFT) calculations (B3LYP functional). A theoretical analysis of a previously used method of estimating exchange constants from a restricted number of spin configurations reveals systematic errors arising from asymmetry in the complexes, which cause large variations in results with different choices of spin configurations. Correction factors are derived that yield the correct results obtained from full configuration space (FCS) calculations. Exchange constants obtained from DFT FCS calculations for the "body-body" interaction were large and ferromagnetic, in disagreement with values obtained from empirical fits of magnetic susceptibility data for the complexes, established magnetostructural correlations in polynuclear Fe(3+) complexes, and ZILSH calculations. DFT calculations also gave unreasonably large antiferromagnetic exchange constants for interaction between "wingtip" ions that are not directly bridged, again in disagreement with ZILSH calculations. Estimates of exchange constants for interaction of body and wingtip ions obtained with ZILSH and DFT were similar, with the ZILSH values in slightly better agreement with empirical fits. Considering all interactions, the ZILSH method provides results in better accord with experiment than DFT for these complexes. Additional comparisons of exchange constants obtained with different spin coupling schemes showed that values appropriate for two-center spin eigenfunctions gave consistently better results than values calculated with the local spin operator. The effect of basis set was found to be very small. A brief analysis of these findings is given.

Density functional perturbational orbital theory of spin polarization in electronic systems. II. Transition metal dimer complexes

We present a theoretical scheme for a semiquantitative analysis of electronic structures of magnetic transition metal dimer complexes within spin density functional theory (DFT). Based on the spin polarization perturbational orbital theory [D.-K. Seo, J. Chem. Phys. 125, 154105 (2006)], explicit spin-dependent expressions of the spin orbital energies and coefficients are derived, which allows to understand how spin orbitals form and change their energies and shapes when two magnetic sites are coupled either ferromagnetically or antiferromagnetically. Upon employment of the concept of magnetic orbitals in the active-electron approximation, a general mathematical formula is obtained for the magnetic coupling constant J from the analytical expression for the electronic energy difference between low-spin broken-symmetry and high-spin states. The origin of the potential exchange and kinetic exchange terms based on the one-electron picture is also elucidated. In addition, we provide a general account of the DFT analysis of the magnetic exchange interactions in compounds for which the active-electron approximation is not appropriate.