Generalized spin orbital GW theory for spin‐frustrated and spin‐degenerate systems

We proposed the ab initio linear combination of Gaussian type orbital (LCGTO) generalized spin orbital GW (GSO–GW) method and calculated triangular hydrogen molecules as models of the noncolinear magnetic clusters. A remarkable improvement of ionized potentials (IPs) by the GW procedure for GHFS solutions is observed in comparison with calculational results by full CI. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 369–374, 2001     

Spin–orbit coupling of spin‐frustrated systems

We have investigated the effects of spin–orbit (SO) interactions on noncollinear molecular magnetism by combining the classical Dzyaloshinsky–Moriya (DM) model and ab initio generalized spin orbital (GSO) method. We have derived an estimation scheme of the magnetic anisotropy energy (MAE) and the Dzyaloshinsky vector based on the SO first‐order perturbation theory (SOPT1) for GSO Hartree–Fock (GHF) solutions. We found that the fundamental results of GHF‐SOPT1 method can be reproduced by diagonalizing the core Hamiltonian plus SO terms, and that the spin topologies of odd‐ring systems can be determined by the topological indices of the singly occupied molecular orbitals. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio calculation of the Dzyaloshinskii–Moriya parameters: Spin–orbit GSO‐HF,DFT, and CI approaches

We present ab initio methods to determine the Dzyaloshinskii–Moriya (DM) parameter, which provides the anisotropic effects of noncollinear spin systems. For this purpose, we explore various general spin orbital (GSO) approaches, such as Hartree–Fock (HF), density functional theory (DFT), and configuration interaction (CI), with one‐electron spin–orbit coupling (SOC1). As examples, two simple D_3h‐symmetric models, H₃ and B(CH₂)₃, are examined. Implications of the computational results are discussed in relation to as isotropic and anisotropic interactions of molecular‐based magnets. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

High‐level ab initio study of the N+(3P)+SH2 reactions in the gas phase: Role of spin‐forbidden pathways

The singlet and triplet potential energy surfaces involved in N⁺+SH₂ reactions have been explored using high‐level ab initio techniques. The geometries of the stationary points were optimized at the QCISD/6‐311G(df,p) level. The final energies were obtained in CCSD(T)/6‐311+G(3df,2p) single‐point calculations. The results obtained show that, although the N⁺(¹D)+SH₂ entrance channel is higher in energy than the N⁺(³P)+SH₂ one, most of the [H₂, S, N]⁺ singlet state cations are lower in energy than the corresponding triplets, due to their different bonding characteristics. Both singlet and triplet potential energy surfaces are quite close each other, and crossover between them can occur. The minimum energy crossing points were located by means of CASSCF(6,5) calculations. The spin‐orbit couplings show that the transition probability from the triplet to the singlet potential energy surface is significantly large. One of the most important consequences is that some of the products of the reaction, such as SH⁺, can be formed in typical spin‐forbidden processes. Since all the relevant structures along these pathways are much lower in energy than the reactants, this mechanism should be accessible even at low impact energies and therefore could be important in processes taking place in interstellar media. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical calculation of the electronic states with spin‐orbit effects of the molecule NaRb

The potential energy curves of the molecule NaRb have been calculated for the 60 low‐lying electronic states in the Ω‐representation. Using an ab‐initio method the calculation is based on nonempirical pseudo‐potential in the interval 3.0a_o ≤ R ≤ 44.0a_o of the internuclear distance. The spin‐orbit effects have been taken into account through a semiempirical spin‐orbit pseudo‐potential added to the electrostatic Hamiltonian with Gaussian basis sets for both atoms. The spectroscopic constants have been calculated for 42 states and the components of the spin‐orbit splitting have been identified for the states (1, 2, 5)³Π and (1, 2)³Δ. The comparison of the present results with those available in literature shows a good agreement, whereas the other results, to the best of our knowledge, are given here for the first time. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A full‐pivoting algorithm for the Cholesky decomposition of two‐electron repulsion and spin‐orbit coupling integrals

A significant reduction in the computational effort for the evaluation of the electronic repulsion integrals (ERI) in ab initio quantum chemistry calculations is obtained by using Cholesky decomposition (CD), a numerical procedure that can remove the zero or small eigenvalues of the ERI positive (semi)definite matrix, while avoiding the calculation of the entire matrix. Conversely, due to its antisymmetric character, CD cannot be directly applied to the matrix representation of the spatial part of the two‐electron spin‐orbit coupling (2e‐SOC) integrals. Here, we present a computational strategy to achieve a Cholesky representation of the spatial part of the 2e‐SOC integrals, and propose a new efficient CD algorithm for both ERI and 2e‐SOC integrals. The proposed algorithm differs from previous CD implementations by the extensive use of a full‐pivoting design, which allows a univocal definition of the Cholesky basis, once the CD δ threshold is made explicit. We show that is the upper limit for the errors affecting the reconstructed 2e‐SOC integrals. The proposed strategy was implemented in the ab initio program Computational Emulator of Rare Earth Systems (CERES), and tested for computational performance on both the ERI and 2e‐SOC integrals evaluation. © 2017 Wiley Periodicals, Inc.     

Density functional study of tetrahedral manganese clusters

We have applied the ab initio density functional theory (DFT) based on generalized spin orbital (GSO) for tetrahedral manganese clusters with a noncollinear spin structure. It was found that three-dimensional (3D) spin states become clearly ground-states by all of GSO-localized spin density approximation, GSO-generalized gradient approximation, and GSO-hybrid methods for tetrahedral Mn(II)₄O₄. These results are consistent with the experimental results.     

Automated reaction path searches for spin‐forbidden reactions

Many catalytic and biomolecular reactions containing transition metals involve changes in the electronic spin state. These processes are referred to as “spin‐forbidden” reactions within nonrelativistic quantum mechanics framework. To understand detailed reaction mechanisms of spin‐forbidden reactions, one must characterize reaction pathways on potential energy surfaces with different spin states and then identify crossing points. Here we propose a practical computational scheme, where only the lowest mixed‐spin eigenstate obtained from the diagonalization of the spin‐coupled Hamiltonian matrix is used in reaction path search calculations. We applied this method to the ^6,4FeO⁺ + H₂ → ^6,4Fe⁺ + H₂O, ^6,4FeO⁺ + CH₄ → ^6,4Fe⁺ + CH₃OH, and ⁷Mn⁺ + OCS → ⁵MnS⁺ + CO reactions, for which crossings between the different spin states are known to play essential roles in the overall reaction kinetics. © 2018 Wiley Periodicals, Inc.     

NMR spin–spin coupling constants in hydrogen‐bonded glycine clusters

The influence of the hydrogen bond formation on the NMR spin–spin coupling constants (SSCC), including the Fermi contact (FC), the diamagnetic spin‐orbit, the paramagnetic spin‐orbit, and the spin dipole term, has been investigated systematically for the homogeneous glycine cluster, in gas phase, containing up to three monomers. The one‐bond and two‐bond SSCCs for several intramolecular (through covalent bond) and intermolecular (across the hydrogen‐bond) atomic pairs are calculated employing the density functional theory with B3LYP and KT3 functionals and different types of extended basis sets. The ab initio SOPPA(CCSD) is used as benchmark for the SSCCs of the glycine monomer. The hydrogen bonding is found to cause significant variations in the one‐bond SSCCs, mostly due to contribution from electronic interactions. However, the nature of variation depends on the type of oxygen atom (proton‐acceptor or proton‐donor) present in the interaction. Two‐bond intermolecular coupling constants vary more than the corresponding one‐bond constants when the size of the cluster increases. Among the four Ramsey terms that constitute the total SSCC, the FC term is the most dominant contributor followed by the paramagnetic spin‐orbit term in all one‐bond interaction.     

Electronic transition moment with spin‐orbit coupling of the molecular ion KRb+

By using the electronic wave functions obtained from an ab initio calculation, including the spin‐orbit coupling, the electronic transition moments have been investigated for two bound states of symmetry Ω = 1/2 and Ω = 3/2 of the molecular ion KRb⁺. Based on a canonical functions approach for the determination of the vibrational wave functions, the matrix elements have been calculated for the bound states considered for v = 0, 10, 20 with v′‐ v = 0, 1, 2, …, 6; by using the same canonical approach, the eigenvalues and abscissas of the corresponding turning points (r_min and r_max) have been investigated for these states that obtained from a theoretical ab initio calculation up to v = 105. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Nuclear spin–spin coupling and nuclear motion

The vibration and rotation of molecules affects nuclear spin–spin coupling constants. This manifests itself as a temperature dependence of the coupling and also as an isotope effect (after allowing, where necessary, for differing magnetogyric ratios of the two nuclei involved in the isotopic substitution). Within the Born–Oppenheimer approximation, a nuclear spin–spin coupling surface can be defined for each pair of coupled nuclei. This surface is sampled by the nuclei as they undergo the excursions about equilibrium geometry that are governed by the force field. An accurate ab initio carbon–proton spin–spin coupling surface for the methane molecule has been calculated. This was obtained by summing the surfaces for each of the four contributions—Fermi contact, spin–dipolar, orbital paramagnetic, and orbital diamagnetic—expressed as power series in terms of symmetry coordinates. Preliminary calculations for ¹³CH₄ and ¹³CD₄ give a difference of only 6% between the calculated and observed nuclear motion contributions. The observed temperature dependence is also accounted for by the calculations. For these isotopomers, bond stretching plays the dominant role. © 1994 John Wiley & Sons, Inc.     

GRECP/MRD‐CI calculations of spin‐orbit splitting in ground state of Tl and of spectroscopic properties of TlH

The generalized relativistic effective core potential (GRECP) approach is employed in the framework of multireference single‐ and double‐excitation configuration interaction (MRD‐CI) method to calculate the spin‐orbit splitting in the ²P^o ground state of the Tl atom and spectroscopic constants for the 0⁺ ground state of TlH. The 21‐electron GRECP for Tl is used, and the outer core 5s and 5p pseudospinors are frozen with the help of the level shift technique. The spin‐orbit selection scheme with respect to relativistic multireference states and the corresponding code are developed and applied in the calculations. In this procedure both correlation and spin‐orbit interactions are taken into account. A [4,4,4,3,2] basis set is optimized for the Tl atom and employed in the TlH calculations. Very good agreement is found for the equilibrium distance, vibrational frequency, and dissociation energy of the TlH ground state (R_e=1.870 Å, ω_e=1420 cm⁻¹, D_e=2.049 eV) as compared with the experimental data (R_e=1.872 Å, ω_e=1391 cm⁻¹, D_e=2.06 eV). © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 409–421, 2001     

Intramolecular Spin Alignment within Mono‐Oxidized and Photoexcited Anthracene‐Based π Radicals as Prototypical Photomagnetic Molecular Devices: Relationships Between Electrochemical,Photophysical, and Photochemical Control Pathways

AnOV is a π‐conjugated radical built from an anthracene (An) unit linked by a p‐phenylene to an oxoverdazyl (OV) moiety. The mono‐oxidized (cationic) form of AnOV was generated both electrochemically and photochemically (in the presence of an electron acceptor). The triplet nature (S=1) of the electronic ground state of AnOV ⁺ was demonstrated by combining spectroelectrochemistry, electron‐spin resonance (ESR) experiments, and ab initio molecular orbital (MO) calculations. The intramolecular spin alignment (ISA) within AnOV ⁺ results from the ferromagnetic coupling (J_electrochem>0) of the two unpaired electrons located on the oxidized electron donor (An⁺) and on the pendant OV radical. The spin‐density distribution pattern of AnOV ⁺ is akin to that of AnOV when photopromoted ( AnOV *) to its high‐spin (HS) lowest excited quartet (S=3/2) state. This high‐spin state results from the ferromagnetic coupling (J_photophys>0) of the triplet locally excited state of An (³An*) with the doublet ground state of OV. As a shared salient feature, AnOV ⁺ and AnOV * (HS) show a spin delocalization within the domain of activated An in either An⁺ or ³An* (nexus states) forms. The present study essentially contributes to establish and clarify relationships between electrochemical, photophysical, and photochemical pathways to achieve ISA processes within AnOV . In particular, we discuss the impact of the spin polarization of the unpaired electron of OV on electronic features of the An electron‐donating subunit. Close analysis of this polarizing interplay allows one to derive a novel functional paradigm to manipulate electron spins at the intramolecular level with light and under an external magnetic field. Indeed, two original functional elements are identified: light‐triggered donors of spin‐polarized electrons and spin‐selective electron acceptors, which are of potential interest for molecular spintronics.     

Stereochemical behavior of 77Se‐1H spin‐spin coupling constants in pyrazolyl‐1,3‐diselenanes and 1,2‐diselenolane

Conformational study of five derivatives of 2‐(pyrazol‐4‐yl)‐1,3‐diselenane together with related 1,2‐diselenolane in respect to the stereochemical trends of geminal and vicinal ⁷⁷Se‐¹H spin‐spin coupling constants has been carried out by means of high‐level theoretical calculations in combination with experiment. The marked dihedral angle dependences for both types of couplings accounted for the lone pair effect in the case of geminal coupling constants and the Karplus‐type relationship for vicinal couplings have been established, which is of major importance for the stereochemical analysis of saturated selenium containing heterocycles. Copyright © 2012 John Wiley & Sons, Ltd.     

Generalized population analysis of three‐center two‐electron bonding

Modern valence bond theory, in its spin‐coupled (SC) form, is used to assess the utility of generalized population analysis, using correlated density matrices, for identifying patterns of three‐center two‐electron bonding in B₂H₆, B₃H, and CH₂Li₂. The effect on the relevant three‐center indices of including electron correlation, via the SC method, is assessed by means of comparison with analogous Hartree–Fock values. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Spin‐orbit ab initio and density functional theory investigation of bismuth monoboronyl,BiBO

The molecular properties of bismuth monoboronyl, BiBO, were investigated using high‐level ab initio and density functional theory calculations by including the effect of spin‐orbit coupling (SOC). SOC does not cause any change in the Bi? B bond length of BiBO, by contrast it causes significant elongation of the Bi? B bond of BiBO^?, by ～0.03 Å. The Bi? B bond length of BiBO^? that is calculated by considering SOC is almost identical to that of BiBO; this result is consistent with a recent experimental study. The term values of excited states of BiBO calculated by including SOC are in good agreement with the experimental results. One excited state which was not assigned in the previous experimental study is the Ω = 0⁺ state generated by strong SOC. In the theoretical calculations on molecules containing 6p‐block elements, including SOC is crucial for obtaining results that are consistent with the corresponding experimental results.     

DNA‐backbone radio resistivity induced by spin blockade effect

Coherent control of OH‐free radicals interacting with the spin‐triplet state of a DNA molecule is investigated. A model Hamiltonian for molecular spin singlet‐triplet resonance is developed. We illustrate that the spin‐triplet state in DNA molecules can be efficiently populated, as the spin‐injection rate can be tuned to be orders of magnitudes greater than the decay rate due to small spin‐orbit coupling in organic molecules. Owing to the nano‐second life‐time of OH free radicals, a non‐equilibrium free energy barrier induced by the injected spin triplet state that lasts approximately longer than one‐micro second in room temperature can efficiently block the initial Hydrogen abstraction and DNA damage. For a direct demonstration of the spin‐blockade effect, a molecular simulation based on an ab‐initio Car‐Parrinello molecular dynamics is deployed. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Magnetic Properties and Electronic Structure of Neptunyl(VI) Complexes: Wavefunctions,Orbitals, and Crystal‐Field Models

The electronic structure and magnetic properties of neptunyl(VI), NpO₂²⁺, and two neptunyl complexes, [NpO₂(NO₃)₃]^? and [NpO₂Cl₄]^2?, were studied with a combination of theoretical methods: ab initio relativistic wavefunction methods and density functional theory (DFT), as well as crystal‐field (CF) models with parameters extracted from the ab initio calculations. Natural orbitals for electron density and spin magnetization from wavefunctions including spin–orbit coupling were employed to analyze the connection between the electronic structure and magnetic properties, and to link the results from CF models to the ab initio data. Free complex ions and systems embedded in a crystal environment were studied. Of prime interest were the electron paramagnetic resonance g‐factors and their relation to the complex geometry, ligand coordination, and nature of the nonbonding 5f orbitals. The g‐factors were calculated for the ground and excited states. For [NpO₂Cl₄]^2?, a strong influence of the environment of the complex on its magnetic behavior was demonstrated. Kohn–Sham DFT with standard functionals can produce reasonable g‐factors as long as the calculation converges to a solution resembling the electronic state of interest. However, this is not always straightforward.     

Towards Magnetic Carbo‐meric Molecular Materials

Numerous studies have underlined the putative diradical character of π‐conjugated molecules that can be described by closed‐shell Lewis structures, for instance, p‐dimethylene p–n phenylenes, or long polyacenes. In the latter compounds, the only way to save the aromaticity of the six‐membered rings is to give up the Lewis electron pairing in the singlet biradical ground state. The present work considers the possibility of doing the same by using the basic C₂ units of carbo‐meric architectures. A series of acyclic and cyclic carbo‐meric architectures is studied by using UB3LYP DFT broken‐symmetry calculations, including spin decontaminations and subsequent geometry optimization of the singlet diradical. The C₂ units are shown to stabilize the singlet biradical by spin delocalization, two of them playing approximately the same role as one radical‐insulating 1,4 phenylene moiety. The results are generalized to the investigation of open‐shell polyradical singlet states of rigid hydrocarbon structures, the symmetry and rigidity of which can assist cooperativity and self spin polarization effect. Several synthesis targets with challenging magnetic/spin properties are suggested in the carbo‐mer series.