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

Can Twisted Double Bonds Facilitate Stepwise [2 + 2] Cycloadditions?

Computational studies for a series of low to high strain anti-Bredt alkenes suggest that those with highly twisted bridgehead double bonds and a small singlet-triplet energy gap may undergo facile stepwise [2 + 2] cycloadditions to furnish four membered rings. A selection of reaction substrates, including ethylene, acetylene, perfluoroethylene, and cyclooctyne are considered.

An unexpected synthesis of azepinone derivatives through a metal-free photochemical cascade reaction

Azepinone derivatives are privileged in organic synthesis and pharmaceuticals. Synthetic approaches to these frameworks are limited to complex substrates, strong bases, high power UV light or noble metal catalysis. We herein report a mild synthesis of azepinone derivatives by a photochemical generation of 2-aryloxyaryl nitrene, [2 + 1] annulation, ring expansion/water addition cascade reaction without using any metal catalyst. Among the different nitrene precursors tested, 2-aryloxyaryl azides performed best under blue light irradiation and Brønsted acid catalysis. The reaction scope is broad and the obtained products underwent divergent transformations to afford other related compounds. A computational study suggests a pathway involving a step-wise aziridine formation, followed by a ring-expansion to the seven-membered heterocycle. Finally, water is added in a regio-selective manner, this is accelerated by the added TsOH.

Machine learning classification of disrotatory IRC and conrotatory non-IRC trajectory motion for cyclopropyl radical ring opening

Quasiclassical trajectory analysis is now a standard tool to analyze non-minimum energy pathway motion of organic reactions. However, due to the large amount of information associated with trajectories, quantitative analysis of the dynamic origin of reaction selectivity is complex. For the electrocyclic ring opening of cyclopropyl radical, more than 4000 trajectories were run showing that allyl radicals are formed through a mixture of disrotatory intrinsic reaction coordinate (IRC) motion as well as conrotatory non-IRC motion. Geometric, vibrational mode, and atomic velocity transition-state features from these trajectories were used for supervised machine learning analysis with classification algorithms. Accuracy >80% with a random forest model enabled quantitative and qualitative assessment of transition-state trajectory features controlling disrotatory IRC versus conrotatory non-IRC motion. This analysis revealed that there are two key vibrational modes where their directional combination provides prediction of IRC versus non-IRC motion.

Spin-projected QM/MM Free Energy Simulations for Oxidation Reaction of Guanine in B-DNA by Singlet Oxygen

Guanine is the most susceptible base to oxidation damage induced by reactive oxygen species including singlet oxygen (¹ O₂ , ¹ Δ_g ). We clarify whether the first step of guanine oxidation in B-DNA proceeds via either a zwitterionic or a diradical intermediate. The free energy profiles are calculated by means of a combined quantum mechanical and molecular mechanical (QM/MM) method coupled with the adaptive biasing force (ABF) method. To describe the open-shell electronic structure of ¹ O₂ correctly, the broken-symmetry spin-unrestricted density functional theory (BS-UDFT) with an approximate spin projection (AP) correction is applied to the QM region. We find that the effect of spin contamination on the activation and reaction free energies is up to ∼8 kcal mol^-1 , which is too large to be neglected. The QM(AP-ULC-BLYP)/MM-based free energy calculations also reveal that the reaction proceeds through a diradical transition state, followed by a conversion to a zwitterionic intermediate. Our computed activation energy of 5.2 kcal mol^-1 matches experimentally observed range (0∼6 kcal mol^-1 ).

Density Functional Theory Study of Mechanisms of [8 + 2] Cycloadditions of Dienylfurans/Dienylisobenzofurans with DMAD

The mechanisms of [8 + 2] cycloaddition reactions between dienylfurans/dienylisobenzofurans and the activated alkyne, DMAD (dimethyl acetylenedicarboxylate), have been investigated by DFT calculations. The former [8 + 2] reaction is stepwise, starting from attack of the diene substituent on furan, not the furyl moiety in dienylfurans, to DMAD to give a diradical intermediate, which then undergoes ring closure to form the second bond between DMAD and the furan moiety, generating the final [8 + 2] cycloadducts. In contrast, the latter [8 + 2] reaction starts from [4 + 2] cycloaddition of the diene in the furan ring of dienylisobenzofurans toward DMAD, followed by the rate-determining stepwise [1,5]-vinyl shift, forming the [8 + 2] products. The different mechanisms of [8 + 2] reactions are attributed to the facts that for dienylfurans, the reactive diene part is the diene substituent on furan, but in the case of dienylisobenzofurans, it is the diene in the furan ring (its reaction with DMAD to generate an aromatic benzene ring is the driving force for this regiochemistry). Consequently, the [8 + 2] reactions begin with the reaction of the most reactive part of tetraene (either the diene substituent on furan for dienylfurans or the diene in the furan ring for dienylisobenzofurans) with DMAD. FMO analysis and kinetic study have been carried out to gain more information of the reaction mechanisms. Two [8 + 2] reactions of dienylisobenzofurans with different substituents toward DMAD have also been further analyzed by DFT calculations in this paper.

Block-Localized Excitation for Excimer Complex and Diabatic Coupling

We describe a block-localized excitation (BLE) method to carry out constrained optimization of block-localized orbitals for constructing valence bond-like, diabatic excited configurations using multistate density functional theory (MSDFT). The method is an extension of the previous block-localized wave function method through a fragment-based ΔSCF approach to optimize excited determinants within a molecular complex. In BLE, both the number of electrons and the electronic spin of different fragments in a whole system can be constrained, whereas electrostatic, exchange, and polarization interactions among different blocks can be fully taken into account of. To avoid optimization collapse to unwanted states, a ΔSCF projection scheme and a maximum overlap of wave function approach have been presented. The method is illustrated by the excimer complex of two naphthalene molecules. With a minimum of eight spin-adapted configurational state functions, it was found that the inversion of L_a^- and L_b^- states near the optimal structure of the excimer complex is correctly produced, which is in quantitative agreement with DMRG-CASPT2 calculations and experiments. Trends in the computed transfer integrals associated with excited-state energy transfer both in the singlet and triplet states are discussed. The results suggest that MSDFT may be used as an efficient approach to treat intermolecular interactions in excited states with a minimal active space (MAS) for interpretation of the results and for dynamic simulations, although the selection of a small active space is often system dependent.

Deciphering the origin of million-fold reactivity observed for the open core diiron [HO-FeIII-O-FeIV[double bond, length as m-dash]O]2+ species towards C-H bond activation: role of spin-states, spin-coupling, and spin-cooperation

High-valent metal-oxo species have been characterised as key intermediates in both heme and non-heme enzymes that are found to perform efficient aliphatic hydroxylation, epoxidation, halogenation, and dehydrogenation reactions. Several biomimetic model complexes have been synthesised over the years to mimic both the structure and function of metalloenzymes. The diamond-core [Fe2(μ-O)2] is one of the celebrated models in this context as this has been proposed as the catalytically active species in soluble methane monooxygenase enzymes (sMMO), which perform the challenging chemical conversion of methane to methanol at ease. In this context, a report of open core [HO(L)FeIII-O-FeIV(O)(L)]2+ (1) gains attention as this activates C-H bonds a million-fold faster compared to the diamond-core structure and has the dual catalytic ability to perform hydroxylation as well as desaturation with organic substrates. In this study, we have employed density functional methods to probe the origin of the very high reactivity observed for this complex and also to shed light on how this complex performs efficient hydroxylation and desaturation of alkanes. By modelling fifteen possible spin-states for 1 that could potentially participate in the reaction mechanism, our calculations reveal a doublet ground state for 1 arising from antiferromagnetic coupling between the quartet FeIV centre and the sextet FeIII centre, which regulates the reactivity of this species. The unusual stabilisation of the high-spin ground state for FeIV[double bond, length as m-dash]O is due to the strong overlap of with the orbital, reducing the antibonding interactions via spin-cooperation. The electronic structure features computed for 1 are consistent with experiments offering confidence in the methodology chosen. Further, we have probed various mechanistic pathways for the C-H bond activation as well as -OH rebound/desaturation of alkanes. An extremely small barrier height computed for the first hydrogen atom abstraction by the terminal FeIV[double bond, length as m-dash]O unit was found to be responsible for the million-fold activation observed in the experiments. The barrier height computed for -OH rebound by the FeIII-OH unit is also smaller suggesting a facile hydroxylation of organic substrates by 1. A strong spin-cooperation between the two iron centres also reduces the barrier for second hydrogen atom abstraction, thus making the desaturation pathway competitive. Both the spin-state as well as spin-coupling between the two metal centres play a crucial role in dictating the reactivity for species 1. By exploring various mechanistic pathways, our study unveils the fact that the bridged μ-oxo group is a poor electrophile for both C-H activation as well for -OH rebound. As more and more evidence is gathered in recent years for the open core geometry of sMMO enzymes, the idea of enhancing the reactivity via an open-core motif has far-reaching consequences.

Prediction on the origin of chemoselectivity in Lewis base-mediated competition cyclizations between allenoates and chalcones: a computational study

The HOMO_TSs and p-orbital contributions of the center atoms were used to predict the origin of chemoselectivity in this work.

Carbene derived diradicaloids - building blocks for singlet fission?

Organic singlet diradicaloids promise application in non-linear optics, electronic devices and singlet fission. The stabilization of carbon allotropes/cumulenes (C1, C2, C4) by carbenes has been equally an area of high activity. Combining these fields, we showed recently that carbene scaffolds allow as well for the design of diradicaloids. Herein, we report a comprehensive computational investigation (CASSCF/NEVPT2; fractional occupation DFT) on the electronic properties of carbene-bridge-carbene type diradicaloids. We delineate how to adjust the properties of these ensembles through the choice of carbene and bridge and show that already a short C2 bridge results in remarkable diradicaloid character. The choice of the carbene separately tunes the energies of the S1 and T1 excited states, whereas the bridge adjusts the overall energy level of the excited states. Accordingly, we develop guidelines on how to tailor the electronic properties of these molecules. Of particular note, fractional occupation DFT is an excellent tool to predict singlet-triplet gaps.

Coordination and Insertion: Competitive Channels for Borylene Reactions

Monovalent boron, free borylene species of the form B-R are notoriously unstable. Consequently, there are substantial gaps in the literature concerning the potential utility of those species in organic and inorganic synthesis either as ligands or as critical intermediates in reactions. We show that the relative stability of borylene complexes varies widely, depending on the electron donating ability of the R group. We find that borylenes can form, in the gas phase, weak sigma hole type interactions to saturated carbon centers and stronger dative bonds to tetravalent silicon and germanium. An insertion reaction of the form FH₃M + BR → FH₂MBHR competes against dative bonding, however, and the reaction is barrierless in several cases when M = Si and in a few cases when M = Ge. For M = C, the barriers are high enough to stabilize monovalent boron complexes. In each case, the barrier heights to M-H bond activation and BR insertion are very sensitive to the nucleophilicity of BR. We confirm, at the MP2(full) and CCSD(T) levels, a substantial preference in borylenes for the singlet over the triplet state. An account is provided at the B3LYP-D3 and MP2(full) levels for the facile insertion reaction on the singlet surface when M = Si and for the stability of FH₃M·BR type complexes and the higher barriers to insertion when M = C and Ge.

Reaction of Aniline with Singlet Oxygen (O21Δg)

Dissolved organic matter (DOM) acts as an effective photochemical sensitizer that produces the singlet delta state of molecular oxygen (O₂¹Δ_g), a powerful oxidizer that removes aniline from aqueous solutions. However, the exact mode of this reaction, the p- to o-iminobenzoquinone ratio, and the selectivity of one over the other remain largely speculative. This contribution resolves these uncertainties. We report, for the first time, a comprehensive mechanistic and kinetic account of the oxidation of aniline with the singlet delta oxygen using B3LYP and M06 functionals in both gas and aqueous phases. Reaction mechanisms have been mapped out at E, H, and G scales. The 1,4-cycloaddition of O₂¹Δ_g to aniline forms a 1,4-peroxide intermediate (M1), which isomerizes via a closed-shell mechanism to generate a p-iminobenzoquinone molecule. On the other hand, the O₂¹Δ_g ene-type reaction forms an o-iminobenzoquinone product when the hydroperoxyl bond breaks, splitting hydroxyl from the 1,2-hydroperoxide (M3) moiety. The gas-phase model predicts the formation of both p- and o-iminobenzoquinones. In the latter model, the M1 adduct displays the selectivity of up to 96%. A water-solvation model predicts that M1 decomposes further, forming only p-iminobenzoquinone with a rate constant of k = 1.85 × 10⁹ (L/(mol s)) at T = 313 K. These results corroborate the recent experimental findings of product concentration profile in which p-iminobenzoquinonine represents the only detected product.

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.

Potential nonadiabatic reactions: ring-opening 4,6-dimethylidenebicyclo[3.1.0]hex-2-ene derivatives to aromatic reactive intermediates

Potential singlet-triplet surface crossings for the ring opening of 4,6-dimethylidenebicyclo[3.1.0]hex-2-ene derivatives were explored using density functional theory (DFT) and complete active space self-consistent field (CASSCF) methods. Since these ring openings involve relatively high energy species that lead to relatively stable aromatic species, a good scenario for potential nonadiabatic events, we posited that the reaction paths of these ring openings might come close to or cross excited state surfaces. At the DFT level of theory, all reaction paths exhibited characteristics suggestive of singlet-triplet intersections along their paths. 6-Methylidenebicyclo[3.1.0]hex-3-en-2-one and a closely related derivative (4-methylidenebicyclo[3.1.0]hex-2-en-6-one) were explored at the CASSCF level of theory; CASSCF results were qualitatively similar to DFT results and yielded spin-orbit couplings of 1.1-1.4 cm(-1) at the singlet-triplet crossing points.