A Stable Aluminum Tris(dithiolene) Triradical.

A stable aluminum tris(dithiolene) triradical (3) was experimentally realized through a low-temperature reaction of the sterically demanding lithium dithiolene radical (2) with aluminum iodide. Compound 3 was characterized by single-crystal X-ray diffraction, UV-vis and EPR spectroscopy, SQUID magnetometry, and theoretical computations. The quartet ground state of triradical 3 has been unambiguously confirmed by variable-temperature continuous wave EPR experiments and SQUID magnetometry. Both SQUID magnetometry and broken-symmetry DFT computations reveal a small doublet-quartet energy gap [ΔEDQ = 0.18 kcal mol-1 (SQUID); ΔEDQ = 0.14 kcal mol-1 (DFT)]. The pulsed EPR experiment (electron spin echo envelop modulation) provides further evidence for the interaction of these dithiolene-based radicals with the central aluminum nucleus of 3.

Correction: Pericyclic reaction benchmarks: hierarchical computations targeting CCSDT(Q)/CBS and analysis of DFT performance.

Correction for 'Pericyclic reaction benchmarks: hierarchical computations targeting CCSDT(Q)/CBS and analysis of DFT performance' by Pascal Vermeeren et al., Phys. Chem. Chem. Phys., 2022, 24, 18028-18042, https://doi.org/10.1039/D2CP02234F.

Structures and Energetics of E2H3+ (E = As, Sb, and Bi) Cations.

E2H2 (E = As, Sb, Bi) structures involving multiple bonds have attracted much attention recently. The E2H3+ cations (protonated E2H2) are predicted to be viable with substantial proton affinities (>180 kcal/mol). Herein, the bonding characters and energetics of a number of E2H3+ isomers are explored through CCSD(T) and DFT methods. For the As2H3+ system, the CCSD(T)/cc-pVQZ-PP method predicts that the vinylidene-like structure lies lowest in energy, with the trans and cis isomers higher by 6.7 and 9.3 kcal/mol, respectively. However, for Sb2H3+ and Bi2H3+ systems, the trans isomer is the global minimum, while the energies of the cis and vinylidene-like structures are higher, respectively, by 2.0 and 2.4 kcal/mol for Sb2H3+ and 1.6 and 15.0 kcal/mol for Bi2H3+. Thus, the vinyledene-like structure is the lowest energy for the arsenic system but only a transition state of the bismuth system. With permanent dipole moments, all minima may be observable in microwave experiments. Besides, we have also obtained transition states and planar-cis structures with higher energies. The current results should provide new insights into the various isomers and provide a number of predictions for future experiments.

Convergent ab initio analysis of the multi-channel HOBr + H reaction.

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br-O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies -0.2 kcal mol-1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol-1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol-1 barrier. The final product energies lie -37.2, -67.9, and -7.3 kcal mol-1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin-orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol-1.

Germanium(II) Dithiolene Complexes.

The 1 : 2 reaction of the imidazole-based dithiolate (2) with GeCl2 • dioxane in THF/TMEDA gives 3, a TMEDA-complexed dithiolene-based germylene. Compound 3 is converted to monothiolate-complexed (5) and N-heterocyclic carbene-complexed (7) germanium(II) dithiolene complexes via Lewis base ligand exchange. A bis-dithiolene-based germylene (8), involving a 3c-4e S-Ge-S bond, has also been synthesized through controlled hydrolysis of 7. The bonding nature of 3, 5, and 8 was investigated by both experimental and theoretical methods.

The multichannel i-propyl + O2 reaction system: A model of secondary alkyl radical oxidation.

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol-1 below the reactants, respectively. Two ß-hydrogen transfer transition states (TS2, TS2') lie above the reactants by (1.4, 2.5) kcal mol-1 and display large Born-Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol-1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol-1 of the corresponding lowest-energy minima.

Concordant Mode Approach for Molecular Vibrations.

The Concordant Mode Approach (CMA) is advanced as a novel hierarchy for increasing the system size and level of theory feasible for quantum chemical computations of harmonic vibrational frequencies. The key concept behind CMA is that transferrable, internal-coordinate normal modes computed at an appropriate lower level of theory (B) comprise a superb basis for converging to vibrational frequencies given by a higher level of theory (A). Accordingly, high-level harmonic frequencies can be evaluated via CMA from a collection of single-point energies that essentially scales linearly in the number of atoms, providing nearly order-of-magnitude CPU time speedups. The accuracy of CMA methods was established by comprehensive tests on over 120 molecules for target Level A = CCSD(T)/cc-pVTZ with auxiliary Level B choices of both CCSD(T)/cc-pVDZ and B3LYP/6-31G(2df,p). Remarkably, the frequency residuals given by the diagonal CMA-0A(nc) scheme exhibit mean absolute deviations (MADs) of only 0.2 cm-1 and standard deviations less than 0.5 cm-1; the corresponding zero-point vibrational energies (ZPVEs) have negligible errors in the vicinity of 0.3 cm-1.

Activation of Ammonia by a Carbene-Stabilized Dithiolene Zwitterion.

A carbene-stabilized dithiolene zwitterion (3) activates ammonia, affording 4• and 5, through both single-electron transfer (SET) and hydrogen atom transfer (HAT). Reaction products were characterized spectroscopically and by single-crystal X-ray diffraction. The mechanism of the formation of 4• and 5 was probed by experimental and computational methods. This discovery is the first example of metal-free ammonia activation via HAT.

Pericyclic reaction benchmarks: hierarchical computations targeting CCSDT(Q)/CBS and analysis of DFT performance.

Hierarchical, convergent ab initio benchmark computations were performed followed by a systematic analysis of DFT performance for five pericyclic reactions comprising Diels-Alder, 1,3-dipolar cycloaddition, electrocyclic rearrangement, sigmatropic rearrangement, and double group transfer prototypes. Focal point analyses (FPA) extrapolating to the ab initio limit were executed via explicit quantum chemical computations with electron correlation treatments through CCSDT(Q) and correlation-consistent Gaussian basis sets up to aug'-cc-pV5Z. Optimized geometric structures and vibrational frequencies of all stationary points were obtained at the CCSD(T)/cc-pVTZ level of theory. The FPA reaction barriers and energies exhibit convergence to within a few tenths of a kcal mol-1. The FPA benchmarks were used to evaluate the performance of 60 density functionals (eight dispersion-corrected), covering the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. The meta-hybrid M06-2X functional provided the best overall performance [mean absolute error (MAE) of 1.1 kcal mol-1] followed closely by the double-hybrids B2K-PLYP, mPW2K-PLYP, and revDSD-PBEP86 [MAE of 1.4-1.5 kcal mol-1]. The regularly used GGA functional BP86 gave a higher MAE of 5.8 kcal mol-1, but it qualitatively described the trends in reaction barriers and energies. Importantly, we established that accurate yet efficient meta-hybrid or double-hybrid DFT potential energy surfaces can be acquired based on geometries from the computationally efficient and robust BP86/DZP level.

Substituted Ortho-Benzynes: Properties of the Triple Bond.

Ortho-benzyne has been well studied by both experiment and theory. Its substituted variants, however, have been less carefully examined. Benchmark data are computed for unsubstituted ortho-benzyne using several density functional theory functionals and basis sets, up to cc-pVQZ. Optimized geometries for the substituted ortho-benzyne as well as harmonic vibrational frequencies and singlet-triplet splittings are computed using the benchmarked functionals. A proximal (syn)OH substitution causes a mean θ1 distortion of +8.1 ± 1.4° from ortho-benzyne. Substituting in the proximal position with F shifts the singlet-triplet splitting by +4.5 ± 0.4 kcal mol-1 from ortho-benzyne. Natural bond orbital analysis, including natural Coulomb electrostatics, elucidates the presence of three influences from the selected substituents: hyperconjugative, resonance, and electrostatic effects.

The bismuth tetramer Bi4: the ν3 key to experimental observation.

The spectroscopic identification of Bi4 has been very elusive. Two constitutional Bi4 isomers of Td and C2v symmetry are investigated and each is found to be a local energetic minimum. The optimized geometries and vibrational frequencies of these two isomers are obtained at the CCSD(T)/cc-pVQZ-PP level of theory, utilizing the Stoll, Metz, and Dolg 60-electron effective core potential. The fundamental frequencies of the Td isomer are obtained at the same level of theory. The focal point analysis method, from a maximum basis set of cc-pV5Z-PP, and proceeding to a maximum correlation method of CCSDTQ, was employed to determine the dissociation energy of Bi4 (Td) into two Bi2 and the adiabatic energy difference between the C2v and Td isomers of Bi4. These quantities are predicted to be +65 kcal mol-1 and +39 kcal mol-1, respectively. Two electron vertical excitation energies between the Td and C2v electronic configurations are computed to be 156 kcal mol-1 for the Td isomer and 9 kcal mol-1 for the C2v isomer. The most probable approach to laboratory spectroscopic identification of Bi4 is via an infrared spectrum. The predicted fundamentals (cm-1) with harmonic IR intensities in parentheses (km mol-1) are 94(0), 123(0.23), and 167(0) for the Td isomer. The moderate IR intensity for the only allowed fundamental may explain why Bi4 has yet to be observed. Through natural bond orbital analysis, the C2v isomer of Bi4 was discovered to exhibit "long-bonding" between the furthest apart 'wing' atoms. This long-bonding is postulated to be facilitated by the σ-bonding orbital between the 'spine' atoms of the C2v isomer.

Unusual nucleophilic reactivity of a dithiolene-based N-heterocyclic silane.

While the dithiolene-based N-heterocyclic silane (4) reacts with two equivalents of BX3 (X = Br, I) to give zwitterionic Lewis adducts 5 and 8, respectively, the parallel reaction of 4 with BCl3 results in 10, a dithiolene-substituted N-heterocyclic silane, via the Si-S bond cleavage. Unlike 5, the labile 8 may be readily converted to 9via BI3-mediated cleavage of the Si-N bond. The formation of 5 and 8 confirms that 4 uniquely possesses dual nucleophilic sites: (a) the terminal sulphur atom of the dithiolene moiety; and (b) the backbone carbon of the N-heterocyclic silane unit.

From Carbene-Dithiolene Zwitterion Mediated B-H Bond Activation to BH3·SMe2-Assisted Boron-Boron Bond Formation.

The 1:1 reaction of the carbene-stabilized dithiolene zwitterion 1 with BH3·SMe2 gave the dithiolene-based hydroborane 2 and the doubly hydrogen-capped CAAC species 3 via hydride-coupled reverse electron transfer processes. The mechanism of this transformation was probed computationally using density functional theory. The subsequent 2:1 reaction of 2 with 1 resulted in 4 and 3, suggesting that 1 can mediate the B-H bond activation not only for BH3 but also for monohydroboranes. In the presence of BH3·SMe2, 2 was unexpectedly converted to the corresponding diborane(4) complex 5 through a dehydrocoupling reaction at an elevated temperature.