Symmetry and broken symmetry in molecular orbital description of unstable molecules IV: comparison between single- and multi-reference computational results for antiaromtic molecules

Issues on DFT+U calculations of organic diradicals

The diradical state is an important electronic state for understanding molecular functions and should be elucidated for the in silico design of functional molecules and their application to molecular devices. The density functional theory calculation with plane-wave basis and correction of the on-site Coulomb parameter U (DFT+U/plane-wave calculation) is a good candidate of high-throughput calculations of diradical-band interactions. However, it has not been investigated in detail to what extent the DFT+U/plane-wave calculation can be used to calculate organic diradicals with a high degree of accuracy. In the present study, using typical organic diradical molecules (bisphenalenyl molecules) as model systems, the discrepancy in the optimum U values between the two electronic states (open-shell singlet and triplet) that compose the diradical state is detected. The calculated results show that the reason for this U value discrepancy is the difference in electronic delocalisation due to π-conjugation between the open-shell singlet and triplet states, and that the effect of U discrepancy becomes large as diradical character decreases. This indicates that it is necessary to investigate the U value discrepancy with reference to the calculated results by more accurate methods or to experimental values when calculating organic diradicals with low diradical character. For this investigation, the local magnetic moments, unpaired beta electron numbers, and effective magnetic exchange integral values can be used as reference values. For the effective magnetic exchange integral values, the effects of U discrepancy are partially cancelled out. However, because the effects may not be completely offset, care should be taken when using the effective magnetic exchange integral value as a reference. Furthermore, a comparison of DFT+U and hybrid-DFT calculations shows that the DFT+U underestimates the HOMO-LUMO gap of bisphenalenyls, although a qualitative discussion of the gap is possible.

Model calculations for the prediction of the diradical character of physisorbed molecules: p-benzyne/MgO and p-benzyne/SrO

The analysis of the diradical state of functional open-shell molecules is important for understanding their physical properties and chemical reactivity. The diradical character is an important factor in the functional elucidation and design of open-shell molecules. In recent years, attempts have been made to immobilise functional open-shell molecules on surfaces to form devices. However, the influence of surface interactions on the diradical state remains unclear. In this study, the physisorption structures of p-benzyne, which is a typical diradical molecule, on MgO(001) and SrO(001) surfaces are used as models to investigate how the diradical character is affected by physisorption. This is done using approximate spin-projected density functional theory calculations with dispersion correction and plane-wave basis (AP-DFT-D3/plane-wave calculations). The diradical character change (Δy) due to adsorption can be categorised into three factors, namely the change due to the distortion of the diradical molecule (Δydis), the interaction between neighbouring diradical molecules (Δycoh), and molecule-surface interactions (Δysurf). In all the calculated models, physisorption reduced the diradical character (Δy < 0), and the contribution of Δysurf was the largest among the three factors. The calculated results show that adsorption induces electron delocalisation to π-conjugated orbitals and intramolecular charge polarisation, both of which contribute to reducing the occupancy of singly occupied molecular orbitals. This indicates that the diradical character of p-benzyne is reduced by the stabilisation of the resonance structures. Furthermore, geometry optimisation of the surfaces shows that the chemical-soft surface (SrO) varies the diradical character more significantly than the chemical-hard surface (MgO). This study shows that the open-shell electronic state and stack structure of diradical molecules can be controlled through the analysis of the surface diradical state.

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction error (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular, we analyze the effect of scaling down the self-interaction correction (SIC) for three nonempirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and the recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger SIC [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048.], the orbitalwise scaled SIC method [Vydrov, O. A. et al. J. Chem. Phys. 2006, 124, 094108.], and the recent local scaling method [Zope, R. R. et al. J. Chem. Phys. 2019, 151, 214108.]. We compute the magnetic exchange coupling constants using the spin projection and nonprojection approaches for sets of molecules composed of dinuclear and polynuclear H···He models, organic radical molecules, and chlorocuprate and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single-electron regions, PZSIC performs well, but for more complex organic systems and the chlorocuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases, LSIC with kinetic energy density ratio performs better than PZSIC. Analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

Diffusion quantum Monte Carlo method on diradicals using single- and multi-determinant-Jastrow trial wavefunctions and different orbitals

In this work, the diffusion quantum Monte Carlo (DMC) method is employed to calculate the energies of singlet and triplet states for a series of organic diradicals and diatomic diradicals with π2 configuration. Single-determinant-Jastrow (SDJ) trial wavefunctions for triplet states, two-determinant-Jastrow (2DJ) trial wavefunctions for the singlet states, and multi-determinant-Jastrow (MDJ) trial wavefunctions are employed in DMC calculations using restricted open-shell B3LYP (ROB3LYP) orbitals, complete-active-space self-consistent field (CASSCF) orbitals, state-average CASSCF orbitals, or frozen-CASSCF orbitals. Our results show that DMC energies using either SDJ/2DJ or MDJ with ROB3LYP orbitals are close to or lower than those with the other orbitals for organic diradicals, while they are not very sensitive to the employed orbitals for diatomic diradicals. Furthermore, using MDJ can reduce DMC energies to some extent for most of the investigated organic diradicals and some diatomic diradicals. The importance of MDJ on DMC energies can be estimated based on the percentage of main determinants in the CASCI wavefunction. On the other hand, singlet–triplet gaps can be calculated reasonably with DMC using MDJ with a mean absolute error of less than 2 kcal/mol with all these orbitals. CASCI wavefunctions using density functional theory orbitals are preferred in constructing MDJ trial wavefunctions in practical DMC calculations since it is easier to obtain such wavefunctions than CASSCF methods.

Singlet/Triplet State Anti/Aromaticity of CyclopentadienylCation: Sensitivity to Substituent Effect

It is well known that singlet state aromaticity is quite insensitive to substituent effects, in the case of monosubstitution. In this work, we use density functional theory (DFT) calculations to examine the sensitivity of triplet state aromaticity to substituent effects. For this purpose, we chose the singlet state antiaromatic cyclopentadienyl cation, antiaromaticity of which reverses to triplet state aromaticity, conforming to Baird’s rule. The extent of (anti)aromaticity was evaluated by using structural (HOMA), magnetic (NICS), energetic (ISE), and electronic (EDDBp) criteria. We find that the extent of triplet state aromaticity of monosubstituted cyclopentadienyl cations is weaker than the singlet state aromaticity of benzene and is, thus, slightly more sensitive to substituent effects. As an addition to the existing literature data, we also discuss substituent effects on singlet state antiaromaticity of cyclopentadienyl cation.

Accurate singlet-triplet gaps in biradicals via the spin averaged anti-Hermitian contracted Schrödinger equation

The accurate description of biradical systems, and in particular the resolution of their singlet-triplet gaps, has long posed a major challenge to the development of electronic structure theories. Biradicaloid singlet ground states are often marked by strong correlation and, hence, may not be accurately treated by mainstream, single-reference methods such as density functional theory or coupled cluster theory. The anti-Hermitian contracted Schrödinger equation (ACSE), whose fundamental quantity is the two-electron reduced density matrix rather than the N-electron wave function, has previously been shown to account for both dynamic and strong correlations when seeded with a strongly correlated guess from a complete active space (CAS) calculation. Here, we develop a spin-averaged implementation of the ACSE, allowing it to treat higher multiplicity states from the CAS input without additional state preparation. We apply the spin-averaged ACSE to calculate the singlet-triplet gaps in a set of small main group biradicaloids, as well as the organic four-electron biradicals trimethylenemethane and cyclobutadiene, and naphthalene, benchmarking the results against other state-of-the-art methods reported in the literature.

Predicted Reversal in N-Methylazepine/N-Methyl-7-azanorcaradiene Equilibrium Upon Formation of Their N-Oxides

Density functional calculations and up to five different basis sets have been applied to the exploration of the structural, enthalpy and free energy changes upon conversion of the azepine to the corresponding N-oxide. Although it is well known that azepines are typically much more stable than their 7-azanorcaradiene valence isomers, the stabilities are reversed for the corresponding N-oxides. Structural, thermochemical as well as nucleus-independent chemical shift (NICS) criteria are employed to probe the potential aromaticity, antiaromaticity and nonaromaticity of N-methylazepine, its 7-azanorcaradiene valence isomer. For the sake of comparison, analogous studies are performed on N-methylpyrrole and its N-oxide.

Singlet-triplet gaps in diradicals obtained with diffusion quantum Monte Carlo using a Slater-Jastrow trial wavefunction with a minimum number of determinants

Diradicals are essential species in a wide range of chemical processes, whereas the computational study of their electronic structure often remains a challenge due to near-degeneracy of the frontier molecular orbitals. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method is employed to calculate adiabatic energy gaps of some typical diradicals with the Slater-Jastrow trial wavefunction. The antisymmetrized part of the trial wavefunction is taken to be a linear combination of a minimum number of determinants using RB3LYP orbitals from the closed-shell singlet state or ROB3LYP orbitals from the triplet state. Our results show that using the two-determinant-Jastrow trial wavefunction is necessary to achieve reliable energy differences between closed-shell singlet states. The energy of the triplet state with MS = 1 is calculated to be lower than that with MS = 0 with FN-DMC even using trial wavefunctions with spin-pure states as their antisymmetrized parts and this difference is reduced with better orbitals. This indicates that the fixed-node error is smaller for the triplet state with MS = 1. Adiabatic energy gaps obtained from the present FN-DMC calculations are in reasonable agreement with available experimental values. Compared with results of the high level EOM-SF-CC method, energy gaps of FN-DMC with RB3LYP orbitals are slightly better than those using ROB3LYP orbitals and results of EOM-SF-CCSD. The present FN-DMC calculations using the simplest ansatz for the trial wavefunction can achieve reasonable results for these diradicals and they can readily be applied to large diradicals.

Fermi-Löwdin orbital self-interaction correction to magnetic exchange couplings

We analyze the effect of removing self-interaction error on magnetic exchange couplings using the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in the framework of density functional theory (DFT). We compare magnetic exchange couplings obtained from self-interaction-free FLOSIC calculations with the local spin density approximation (LSDA) with several widely used DFT realizations and wave function based methods. To this end, we employ the linear H-He-H model system, six organic radical molecules, and [Cu₂Cl₆]^2- as representatives of different types of magnetic interactions. We show that the simple self-interaction-free version of LSDA improves calculated couplings with respect to LSDA in all cases, even though the nature of the exchange interaction varies across the test set, and in most cases, it yields results comparable to modern hybrids and range-separated approximate functionals.

A Bis(μ-oxido)dinickel(III) Complex with a Triplet Ground State

A bis(μ-oxido)dinickel(III) complex was synthesized and characterized by single crystal X-ray diffraction, resonance Raman, and ESI-mass measurements. Magnetic susceptibility measurements by SQUID and EPR spectroscopy reveal that the complex has a triplet ground state, which is unprecedented for high-valent metal (M) complexes with [M₂ (μ-O)₂ ] diamond core. DFT studies indicate ferromagnetic coupling of the nickel(III) centers. The complex exhibits hydrogen abstraction reactivity and oxygenation reactivity toward external substrates.

Correlated-Participating-Orbitals Pair-Density Functional Method and Application to Multiplet Energy Splittings of Main-Group Divalent Radicals

Predicting the singlet-triplet splittings of divalent radicals is a challenging task for electronic structure theory. In the present work, we investigate the performance of multiconfiguration pair-density functional theory (MC-PDFT) for computing the singlet-triplet splitting for small main-group divalent radicals for which accurate experimental data are available. In order to define theoretical model chemistries that can be assessed consistently, we define three correlated participating orbitals (CPO) schemes (nominal, moderate, and extended, abbreviated as nom, mod, and ext) to define the constitution of complete active spaces, and we test them systematically. Broken-symmetry Kohn-Sham DFT calculations have also been carried out for comparison. We found that the extended CPO-PDFT scheme with translated on-top pair-density functionals have smaller mean unsigned errors than weighted-average broken-symmetry Kohn-Sham DFT with the corresponding exchange-correlation functional. The accuracy of the translated Perdew-Burke-Ernzerhof (tPBE) on-top pair-density functionals with ext-CPO active space is even better than some of the more accurately parametrized exchange-correlation density functionals that we tested; this is very encouraging for MC-PDFT theory.

Unrestricted prescriptions for open-shell singlet diradicals: using economical ab initio and density functional theory to calculate singlet-triplet gaps and bond dissociation curves

Here we present and test several computational prescriptions for calculating singlet-triplet (ST) gap energies and bond dissociation curves for open-shell singlet diradicals using economical unrestricted single reference type calculations. For ST gap energies from Slipchenko and Krylov's atom and molecule test set (C, O, Si, NH, NF, OH(+), O(2), CH(2), and NH(2)(+)) spin unrestricted Hartree-Fock and MP2 energies result in errors greater than 15 kcal/mol. However, spin-projected (SP) Hartree-Fock theory in combination with spin-component-scaled (SCS) or scaled-opposite-spin (SOS) second-order perturbation theory gives ST gap energies with a mean unsigned error (MUE) of less than 2 kcal/mol. Density functionals generally give poor results for unrestricted energies and only the ωB97X-D, the M06, and the M06-2X functionals provide reasonable accuracy after spin-projection with MUE values of 4.7, 4.3, and 3.0 kcal/mol, respectively, with the 6-311++G(2d,2p) basis set. We also present a new one parameter hybrid density functional, diradical-1 (DR-1), based on Adamo and Barone's modified PW exchange functional with the PW91 correlation functional. This DR-1 method gives a mean error (ME) of 0.0 kcal/mol and a MUE value of 1.3 kcal/mol for ST gap energies. As another test of unrestricted methods the bond dissociation curves for methane (CH(4)) and hydrofluoric acid (H-F) were calculated with the M06-2X, DR-1, and ωB97X-D density functionals. All three of these functionals give reasonable results for the methane C-H bond but result in errors greater than 50 kcal/mol for the H-F bond dissociation. Spin-projection is found to significantly degrade bond dissociation curves past ~2.2 Å. Although unrestricted Hartree-Fock theory provides a very poor description of H-F bond dissociation, unrestricted SCS-MP2 and SOS-MP2 methods give accurate results.