Electric Field Gradient in Chiral and Tetrahedral Molecules within High-Order LRESC Formalism.

In this work, we present the electric field gradient (EFG) given by the linear response elimination of the small component (LRESC) scheme up to the 1/c4 order (c is the speed of light in vacuum) in CHFClX (X = Br, I, At) chiral molecules, together with CHF2Br and CH2FX (X = Br, I, At) tetrahedral systems. The former could be good candidates for further parity violation studies, especially when heavy atoms are surrounding. In this context, the LRESC scheme demonstrates effective applicability to large tetrahedral and chiral molecules that incorporate heavy elements, with relativistic effects playing a crucial role. The LRESC results of EFG exhibit an excellent agreement with those calculated at the four-component level, giving differences of only hundredths order in a.u. (atomic units) for the bromine nucleus and less than 0.1 a.u. for the iodine nucleus. Regarding the other nuclei, for the chiral molecules, there is a heavy atom effect on the light atom (HALA) for chlorine and fluorine atoms as the substituent halogen atom becomes heavier. Furthermore, the electronic part of the EFG for the central carbon and the fluorine nuclei presents an important dependence with the environment in the molecules under study. With accurate calculations of the EFG and tabulated nuclear quadrupole moment, the nuclear quadrupole coupling constant is obtained within the LRESC scheme, including for the first time correlation effects on the spin-dependent corrections with this methodology, providing results close to the experimental ones for Cl, Br, and I atoms. At the Hartree-Fock level, the differences are around 6% for Br and I nuclei, and at the density functional theory level with the LDA and PBE0 functionals, the differences are no more than 2%.

DOCSIC: A Mean-Field Method for Orbital-by-Orbital Self-Interaction Correction.

We introduce a new method to remove the one-electron self-interaction error in approximate density functional calculations on an orbital-by-orbital basis, as originally proposed by Perdew and Zunger [Phys. Rev. B 1981, 23, 5048]. This method is motivated by a recent proposal by Pederson et al. [J. Chem. Phys. 2014, 140, 121103] to remove self-interaction that employs orbitals derived from the real-space density matrix, known as FLOSIC (Fermi Löwdin orbitals self-interaction correction). However, instead of Fermi Löwdin orbitals, our scheme utilizes columns of the density matrix to determine localized orbitals, like the localization procedure proposed by Fuemmeler et al. [J. Chem. Theory Comput. 2023, 19, 8572]. The new method, dubbed DOCSIC for density matrix as orbital coefficients self-interaction correction, contrasts with traditional Perdew-Zunger or FLOSIC in that it does not incorporate additional optimization parameters, and, unlike the average density self-interaction correction of Ciofini et al. [Chem. Phys. Lett. 2003, 380, 12], it makes use of localized orbitals. Another advantage of DOCSIC is that it can be implemented as a mean-field formalism. We show details of the self-consistent generalized Kohn-Sham implementation, some illustrative results, and we finally highlight its advantages and limitations.

Generalized Spin in the Variance-Based Wave Function Optimization Method within the Doubly Occupied Configuration Interaction Framework.

In this work, we implement a generalized spin formulation of the doubly occupied configuration interaction methodology using the energy variance of the N-electron Hamiltonian. We perform the optimization of the N-electron wave functions and calculate their corresponding energies, using a unified variational treatment for ground and excited states based on the energy variance, which allows us to describe the entire energy spectra on an equal footing. We analyze the effects produced by the breakdown of the S2 and Sz symmetries in the spectra of model hydrogenic clusters in terms of energies and spin-related quantities, arising from the restricted, unrestricted, and generalized spin methods. The results are compared with other related methods as well as full configuration interaction.

Density Matrix Implementation of the Fermi-Löwdin Orbital Self-Interaction Correction Method.

The Fermi-Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi-Löwdin orbitals which are used to remove the self-interaction error in the Perdew-Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi-Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn-Sham with one additional term in the Kohn-Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn-Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. We analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville-von Neumann equation to calculate excitation energies.

Generalized spin σ-SCF method.

We introduce a generalization of the σ-SCF method to approximate noncollinear spin ground and excited single-reference electronic states by minimizing the Hamiltonian variance. The new method is based on the σ-SCF method, originally proposed by Ye et al. [J. Chem. Phys. 147, 214104 (2017)], and provides a prescription to determine ground and excited noncollinear spin states on an equal footing. Our implementation was carried out utilizing an initial simulated annealing stage followed by a mean-field iterative self-consistent approach to simplify the cumbersome search introduced by generalizing the spin degrees of freedom. The simulated annealing stage ensures a broad exploration of the Hilbert space spanned by the generalized spin single-reference states with random complex element-wise rotations of the generalized density matrix elements in the simulated annealing stage. The mean-field iterative self-consistent stage employs an effective Fockian derived from the variance, which is utilized to converge tightly to the solutions. This process helps us to easily find complex spin structures, avoiding manipulating the initial guess. As proof-of-concept tests, we present results for Hn (n = 3-7) planar rings and polyhedral clusters with geometrical spin frustration. We show that most of these systems have noncollinear spin excited states that can be interpreted in terms of geometric spin frustration. These states are not directly targeted by energy minimization methods, which are meant to converge to the ground state. This stresses the capability of the σ-SCF methodology to find approximate noncollinear spin structures as mean-field excited states.

High order relativistic corrections on the electric field gradient within the LRESC formalism.

In this work, we present relativistic corrections to the electric field gradient (EFG) given by the Linear Response Elimination of the Small Component (LRESC) scheme at 1/c2 order and including for the first time spin-dependent (SD) corrections at 1/c4 order. We show that these new terms improve the performance of LRESC as results with this methodology are very close to those calculated at the four-component Dirac-Hartree-Fock (4c-DHF) level. We assess the new corrections in BrY and AtY di-halogen (Y = F, Cl, Br, I, and At) and XZY bi-linear molecules (Z = Zn, Cd, and Hg; X, Y = F, Cl, Br, I, and At). At the 4c-DHF level, we analyze the contributions coming from the large and small components of the relativistic 4c wave function to the electronic part of EFG and compare them with the LRESC corrections to find their electronic origin. For the HgX2 (X = Cl, Br, and I) subset, when the SD correcting terms are included, LRESC calculations match very well with 4c-DHF ones and those from the literature, with differences less than 1% for molecules containing up to three heavy atoms. We show that LRESC gives accurate values of EFG, allowing the analysis of the electronic origin of relativistic effects in terms of well-known nonrelativistic operators.

Performance of the LRESC Model on top of DFT Functionals for Relativistic NMR Shielding Calculations.

The linear response within the elimination of the small component model (LRESC) is an insightful and computationally efficient method for including relativistic effects on molecular properties like the nuclear magnetic shielding constants, spin-rotation constant, g-tensor, and electric field gradient of heavy atom containing molecules with atoms belonging up to the sixth row of the periodic table. One of its main advantages is its capacity to analyze the electronic origin of the different relativistic correcting terms. Until now, it was always applied on top of Hartree-Fock ground-state wave functions (LRESC/HF) to calculate and analyze NMR shieldings. In this work, we show the performance of the LRESC formalism on top of some density functional theory (DFT) functionals to compute tin shielding constants in SnX4 (X = H, F, Cl, Br, I) molecular systems. We analyze the performance of each LRESC/DFT scheme on reproducing the electronic mechanisms of the shieldings, taking as a benchmark the results of relativistic calculations at the RPA level of approach (4c/RPA). As in previous works, we divide the LRESC relativistic correcting terms into two groups: core-dependent and ligand-dependent contributions. It is shown here that core-dependent corrections are well-reproduced for the selected DFT functionals, but some differences arise in the ligand-dependent ones. We focus on the performance of different functionals, including the same electron correlation part but containing different amounts of HF exchange. The best results are obtained for the BHandHLYP functional (50% of HF exchange) and the worst for BLYP (0%). When the percentage of HF exchange increases, ligand-dependent contributions are better described, and the final LRESC/DFT results are closer to those obtained with LRESC/HF and 4c/RPA methods. The spin-orbit correction to the shielding constant is one of the main ligand-dependent contributions (there are two more) with total value depending on the amount of HF exchange included in the functional. When the amount of HF exchange decreases, the spin-orbit contribution becomes larger, overestimating the shielding constant even when nonrelativisitc DFT values are much smaller than the nonrelativistic HF ones, as it happens for the heaviest molecular system studied here (SnI4).

Magnetic Properties of Co(II) Complexes with Polyhedral Carborane Ligands.

In this work we present a computational analysis of a new family of magnetic Co(II) single-ion complexes with large magnetic anisotropy based on icosahedral and octahedral carborane ligands. In particular, we extend our previous computational work on mononuclear Co(II) complexes with 1,2-(HS)2-1,2-C2B10H10 and 9,12-(HS)2-1,2-C2B10H10 icosahedral o-carborane ligands to a larger set of complexes where the Co(II) ion is doubly chelated by those ligands and by other two positional isomers belonging to the 1,2-dicarba- closo-dodecaborane family. We also describe Co(II) complexes with octahedral ligands derived from 1,2-dicarba- closo-hexaborane and study the effects of replacing a thiol group by a hydroxy group in both polyhedral geometries, as well as the influence of the position of the carbon atoms. On analysis of the results for a total of 20 complexes, our results show that carborane-based Co(II) single-ion compounds present a distorted-tetrahedral geometry, high-spin ground states, and high values for the magnetic anisotropy parameters. We point out which of these would be suitable candidates to be synthesized and used as molecular magnets.

Magnetic Properties of Mononuclear Co(II) Complexes with Carborane Ligands.

We analyze the magnetic properties of three mononuclear Co(II) coordination complexes using quantum chemical complete active space self-consistent field and N-electron valence perturbation theory approaches. The complexes are characterized by a distorted tetrahedral geometry in which the central ion is doubly chelated by the icosahedral ligands derived from 1,2-(HS)2-1,2-C2B10H10 (complex I), from 1,2-(HS)2-1,2-C2B10H10 and 9,12-(HS)2-1,2-C2B10H10 (complex II), and from 9,12-(HS)2-1,2-C2B10H10 (complex III), which are two positional isomers of dithiolated 1,2-dicarba- closo-dodecaborane (complex I). Complex I was realized experimentally recently (Tu, D.; Shao, D.; Yan, H.; Lu, C. Chem. Commun. 2016, 52, 14326) and served to validate the computational protocol employed in this work, while the remaining two proposed complexes can be considered positional isomers of I. Our calculations show that these complexes present different axial and rhombic zero-field splitting anisotropy parameters and different values of the most significant components of the g tensor. The predicted axial anisotropy D = -147.2 cm-1 for complex II is twice that observed experimentally for complex I, D = -72.8 cm-1, suggesting that this complex may be of interest for practical applications. We also analyze the temperature dependence of the magnetic susceptibility and molar magnetization for these complexes when subject to an external magnetic field. Overall, our results suggest that o-carborane-incorporated Co(II) complexes are worthwhile candidates for experimental exploration as single-ion molecular magnets.

Theoretical analysis of NMR shieldings of group-11 metal halides on MX (M = Cu, Ag, Au; X = H, F, Cl, Br, I) molecular systems, and the appearance of quasi instabilities on AuF.

Accurate calculations of nuclear magnetic shieldings of group-11 metal halides, σ(M; MX) (M = Cu, Ag, Au; X = H, F, Cl, Br, I), were performed with relativistic and nonrelativistic theoretical schemes in order to learn more about the importance of the involved electronic mechanisms that underlie such shieldings. We applied state of the art schemes: polarization propagators at a random phase level of approach (PP-RPA); spin-free Hamiltonian (SF); linear response elimination of small component (LRESC) and density functional theory (DFT) with two different functionals: B3LYP and PBE0. The results from DFT calculations are not close to those from the relativistic polarization propagator calculations at the RPA level of approach (RelPP-RPA), in line with previous results. The spin-orbit (SO) contribution to a shielding constant is important only for MF molecules (M = Cu, Ag, Au). Different electronic mechanisms are considered within the LRESC method, bunched into two groups: core- and ligand-dependent. For the analysed shieldings the core-dependent electronic mechanisms are the most important ones; the ligand-dependent being only important for MF molecules. An out of range value for σ(Au) is found in AuF. It was previously reported in the literature, either originated in the large fluorine electronegativity together with large spin-orbit coupling contributions; or, due to Fermi-contact contributions. We argue here that such an unexpected large value is an artifact originated in the appearance of quasi instabilities, and show how to handle this apparent problem.