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.

Exploring the Tl 2 H 2 potential energy surface: A comparative analysis with group 13 systems and experiment.

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl 2 H 2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug-cc-pwCVQZ-PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero-point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M 2 H 2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

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.

A wealth of structures for the Ge2H2+ radical cation: comparison of theory and experiment.

Five structures of Ge2H2 and Ge2H2+ are investigated in this study. Optimized geometries at the CCSD(T)/cc-pwCVQZ-PP level of theory were obtained. Focal point analyses were performed on these optimized geometries to determine relative energies using the CCSD(T) method with polarized basis sets up to quintuple-zeta. Energy corrections include full T and pertubative (Q) coupled-cluster effects plus anharmonic corrections to the zero-point vibrational energy. Relative ordering in energy from lowest to highest of the five Ge2H2+ structures is butterfly, germylidene, monobridged, trans, then linear. In neutral Ge2H2, the monobridged structure lies lower in energy than the germylidene structure. Fundamental vibrational frequencies and IR intensities were computed for the minima at the CCSD(T)/cc-pwCVTZ-PP level of theory to compare with experimental research. Partial atomic charges and natural bonding orbital analyses indicated that the positive charge of Ge2H2+ is contained in the region of the Ge-Ge bond.

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.

Lantern-type dinickel complexes: An exploration of possibilities for nickel-nickel bonding with bridging bidentate ligands.

Many binuclear nickel complexes have NiNi distances suggesting NiNi covalent bonds, including lantern-type complexes with bridging bidentate ligands. This DFT study treats tetragonal, trigonal, and digonal lantern-type complexes with the formamidinate, guanidinate, and formate ligands, besides some others. Formal bond orders (ranging from zero to two) are assigned to all the NiNi bonds on the basis of MO occupancy considerations. A VB-based electron counting approach assigns plausible resonance structures to the dinickel cores. Model tetragonal complexes with the dimethylformamidinate and the dithioformate ligands have singlet ground states whose non-covalently bonded NiNi distances are close to those in their experimentally known counterparts. Trigonal dinickel complexes are unknown, but are predicted to have quartet ground states with NiNi bonds of order 0.5. The model digonal complexes are predicted to have triplet ground states, but the predicted NiNi bond lengths are longer than those found in their experimentally known counterparts. This could owe to inadequate treatment of electron correlation by DFT in these short NiNi bonds with their multiconfigurational character. All the NiNi bond distances here are categorized into ranges according to the NiNi bond orders of 0, 0.5, 1, 1.5, and 2, no NiNi bonds of order higher than two being identified. The NiNi bonds of given order in these lantern-type complexes are consistently shorter than the corresponding NiNi bonds in dinickel complexes having carbonyl ligands, attributable to the metalmetal bond lengthening effect of CO ligands.

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 highly exothermic hydrogen abstraction reaction H2Te + OH → H2O + TeH: comparison with analogous reactions for H2Se and H2S.

The "gold standard" CCSD(T) method is adopted along with the correlation consistent basis sets up to aug-cc-pV5Z-PP to study the mechanism of the hydrogen abstraction reaction H2Te + OH. The predicted geometries and vibrational frequencies for reactants and products are in good agreement with the available experimental results. With the ZPVE corrections, the transition state in the favorable pathway of this reaction energetically lies 1.2 kcal mol-1 below the reactants, which is lower than the analogous relative energies for the H2Se + OH reaction (-0.7 kcal mol-1), the H2S + OH reaction (+0.8 kcal mol-1) and the H2O + OH reaction (+9.0 kcal mol-1). Accordingly, the exothermic reaction energies for these related reactions are predicted to be 47.8 (H2Te), 37.7 (H2Se), 27.1 (H2S), and 0.0 (H2O) kcal mol-1, respectively. Geometrically, the low-lying reactant complexes for H2Te + OH and H2Se + OH are two-center three-electron hemibonded structures, whereas those for H2S + OH and H2O + OH are hydrogen-bonded. With ZPVE and spin-orbit coupling corrections, the relative energies for the reactant complex, transition state, product complex, and the products for the H2Te + OH reaction are estimated to be -13.1, -1.0, -52.0, and -52.6 kcal mol-1, respectively. Finally, twenty-eight DFT functionals have been tested systematically to assess their ability in describing the potential energy surface of the H2Te + OH reaction. The best of these functionals for the corresponding energtics are -9.9, -1.4, -46.4, and -45.4 kcal mol-1 (MPWB1K), or -13.1, -2.4, -57.1, and -54.6 kcal mol-1 (M06-2X), respectively.

High-Level Coupled-Cluster Study on Substituent Effects in H2 Activation by Low-Valent Aluminyl Anions.

The synthesis of novel aluminyl anion complexes has been well exploited in recent years. Moreover, the elucidation of the structure and reactivity of these complexes opens the path toward a new understanding of low-valent aluminum complexes and their chemistry. This work computationally treats the substituent effect on aluminyl anions to discover suitable alternatives for H2 activation at a high level of theory utilizing coupled-cluster techniques extrapolated to the complete basis set. The results reveal that the simplest AlH2- system is the most reactive toward the activation of H2, but due to the low steric demand, severe difficulty in the stabilization of this system makes its use nonviable. However, the results indicate that, in principle, aluminyl systems with -C, -CN, -NC, and -N chelating centers would be the best choices of ligand toward the activation of molecular hydrogen by taking care of suitable steric demand to prevent dimerization of the catalysts. Furthermore, computations show that monosubstitution (besides -H) in aluminyl anions is preferred over disubstitution. So our predictions show that bidentate ligands may yield less reactive aluminyl anions to activate H2 than monodentate ones.

Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF4-(H2O)n=1,2,3,4.

This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.

Simulation of the VUV Absorption Spectra of Oxygenates and Hydrocarbons: A Joint Theoretical-Experimental Study.

Vacuum UV absorption spectroscopy is regularly used to provide unambiguous identification of a target species, insight into the electronic structure of molecules, and quantitative species concentrations. As molecules of interest have become more complex, theoretical spectra have been used in tandem with laboratory spectroscopic analysis or as a replacement when experimental data is unavailable. However, it is difficult to determine which theoretical methodologies can best simulate experiment. This study examined the performance of EOM-CCSD and 10 TD-DFT functionals (B3LYP, BH&HLYP, BMK, CAM-B3LYP, HSE, M06-2X, M11, PBE0, ωB97X-D, and X3LYP) to produce reliable vacuum UV absorption spectra for 19 small oxygenates and hydrocarbons using vertical excitation energies. The simulated spectra were analyzed against experiment using both a qualitative analysis and quantitative metrics, including cosine similarity, relative integral change, mean signed error, and mean absolute error. Based on our ranking system, it was determined that M06-2X was consistently the top performing TD-DFT method with BMK, CAM-B3LYP, and ωB97X-D also producing reliable spectra for these small combustion species.

Comparison of multifidelity machine learning models for potential energy surfaces.

Multifidelity modeling is a technique for fusing the information from two or more datasets into one model. It is particularly advantageous when one dataset contains few accurate results and the other contains many less accurate results. Within the context of modeling potential energy surfaces, the low-fidelity dataset can be made up of a large number of inexpensive energy computations that provide adequate coverage of the N-dimensional space spanned by the molecular internal coordinates. The high-fidelity dataset can provide fewer but more accurate electronic energies for the molecule in question. Here, we compare the performance of several neural network-based approaches to multifidelity modeling. We show that the four methods (dual, Δ-learning, weight transfer, and Meng-Karniadakis neural networks) outperform a traditional implementation of a neural network, given the same amount of training data. We also show that the Δ-learning approach is the most practical and tends to provide the most accurate model.

Enthalpies of formation for Criegee intermediates: A correlation energy convergence study.

Criegee intermediates, formed from the ozonolysis of alkenes, are known to have a role in atmospheric chemistry, including the modulation of the oxidizing capacity of the troposphere. Although studies have been conducted since their discovery, the synthesis of these species in the laboratory has ushered in a new wave of investigations of these structures, both theoretically and experimentally. In some of these theoretical studies, high-order corrections for correlation energy are included to account for the mid multi-reference character found in these systems. Many of these studies include a focus on kinetics; therefore, the calculated energies should be accurate (<1 kcal/mol in error). In this research, we compute the enthalpies of formation for a small set of Criegee intermediates, including higher-order coupled cluster corrections for correlation energy up to coupled cluster with perturbative quintuple excitations. The enthalpies of formation for formaldehyde oxide, anti-acetaldehyde oxide, syn-acetaldehyde oxide, and acetone oxide are presented at 0 K as 26.5, 15.6, 12.2, and 0.1 kcal mol-1, respectively. Additionally, we do not recommend the coupled cluster with perturbative quadruple excitations [CCSDT(Q)] energy correction, as it is approximately twice as large as that of the coupled cluster with full quadruple excitations (CCSDTQ). Half of the CCSDT(Q) energy correction may be included as a reliable, cost-effective estimation of CCSDTQ energies for Criegee intermediates.

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.

Probing the Potential Energy Profile of the I + (H2O)3 → HI + (H2O)2OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included.

Three different pathways for the atomic iodine plus water trimer reaction I + (H2O)3 → HI + (H2O)2OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H2O)3 reactions, three pathways for the I + (H2O)3 reaction are linked in terms of geometry and energetics. To legitimize the results, the "gold standard" CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H2O)3 → HI + (H2O)2OH reaction is predicted to be endothermic by 47.0 kcal mol-1. The submerged transition state is predicted to lie 43.7 kcal mol-1 above the separated reactants. The I···(H2O)3 entrance complex lies below the separated reactants by 4.1 kcal mol-1, and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H2O)2OH exit complex is bound by 4.3 kcal mol-1 in relation to the separated products. Compared with simpler I + (H2O)2 and I + H2O reactions, the I + (H2O)3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H2O)3 reaction. The iodine reaction I + (H2O)3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H2O)3 but much different from the F + (H2O)3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

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.

Acceleration Effect of Bases on Mn Pincer Complex-Catalyzed CO2 Hydroboration.

Herein, we report a comprehensive study of CO2 hydroboration catalyzed by Mn pincer complexes. The traditional metal-ligand cooperation (MLC) mechanism based on the H-Mn-N-Bpin pincer complex is not viable due to the competing abstraction of the Bpin group from the H-Mn-N-Bpin complex by NaOtBu. Instead, we propose an ionic mechanism based on the H-Mn-N-Na species with a low energy span (22.5 kcal/mol) and unveil the acceleration effect of bases. The X groups in the H-Mn-N-X catalyst models are further modulated, and the steric hindrance and H→B donor-acceptor interactions of the X group increase the energy barrier of the hydride transfer. The hydrogen bond and electrostatic interactions of the X group can accelerate the hydride transfer to HCOOBpin and HCHO molecules except for the nonpolar CO2 molecule. Based on these discoveries, we designed a pyridine-based Mn pincer catalyst system, which could achieve CO2 hydroboration in low-temperature and base-free conditions through a metal-ligand cooperation mechanism.

The reaction between the bromine atom and the water trimer: high level theoretical studies.

Three different reaction pathways are found for the reaction of bromine atom (Br) with the lowest-energy structure of the water trimer [uud-(H2O)3], initially using the MPW1K-DFT method. The three bromine pathways have closely related geometries and energetics, analogous to those found for the fluorine and chlorine reactions. The lowest-energy pathway of the Br + uud-(H2O)3 reaction was further investigated using the "gold standard" CCSD(T) method and the correlation-consistent basis sets up to cc-pVQZ(-PP). Based on the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the Br + (H2O)3 reaction is endothermic by 33.3 kcal mol-1. The classical barrier height is 29.0 kcal mol-1 between the reactants and the exit complex, and there is no barrier for the reverse reaction. The Br⋯(H2O)3 entrance complex is found to lie 4.7 kcal mol-1 below the separated reactants, and the HBr⋯(H2O)2OH exit complex is bound by 6.4 kcal mol-1 relative to the separated products. This potential energy profile is further corrected by the zero point energies and spin-orbit coupling effects. Structurally, the Br + (H2O)3 stationary points can be derived from those of the simpler Br + (H2O)2 reaction by judiciously appending a H2O molecule. The Br + (H2O)3 potential energy profile is compared with the Br + (H2O)2 and Br + H2O reactions, as well as to the valence isoelectronic Cl + (H2O)3 and F + (H2O)3 systems. It is reasonable to expect that the reactions between the bromine atom and larger water clusters would be similar to the Br + (H2O)3 reaction.