Fulminic acid: a quasibent spectacle.

Fulminic acid (HCNO) played a critical role in the early development of organic chemistry, and chemists have sought to discern the structure and characteristics of this molecule and its isomers for over 200 years. The mercurial nature of the extremely flat H-C-N bending potential of fulminic acid, with a nearly vanishing harmonic vibrational frequency at linearity, remains enigmatic and refractory to electronic structure theory, as dramatic variation with both orbital basis set and electron correlation method is witnessed. To solve this problem using rigorous electronic wavefunction methods, we have employed focal point analyses (FPA) to ascertain the ab initio limit of optimized linear and bent geometries, corresponding vibrational frequencies, and the HCN + O(3P) → HCNO reaction energy. Electron correlation treatments as extensive as CCSDT(Q), CCSDTQ(P), and even CCSDTQP(H) were employed, and complete basis set (CBS) extrapolations were performed using the cc-pCVXZ (X = 4-6) basis sets. Core electron correlation, scalar relativistic effects (MVD1), and diagonal Born-Oppenheimer corrections (DBOC) were all included and found to contribute significantly in determining whether vibrationless HCNO is linear or bent. At the all-electron (AE) CCSD(T)/CBS level, HCNO is a linear molecule with ω5(H-C-N bend) = 120 cm-1. However, composite AE-CCSDT(Q)/CBS computations give an imaginary frequency (51i cm-1) at the linear optimized geometry, leading to an equilibrium structure with an H-C-N angle of 173.9°. Finally, at the AE-CCSDTQ(P)/CBS level, HCNO is once again linear with ω5 = 45 cm-1, and inclusion of both MVD1 and DBOC effects yields ω5 = 32 cm-1. A host of other topics has also been investigated for fulminic acid, including the dependence of re and ωi predictions on a variety of CBS extrapolation formulas, the question of multireference character, the N-O bond energy and enthalpy of formation, and issues that give rise to the quasibent appellation.

Trinuclear and Tetranuclear Ruthenium Carbonyl Nitrosyls: Oxidation of a Carbonyl Ligand by an Adjacent Nitrosyl Ligand.

Trinuclear and tetranuclear ruthenium carbonyls of the types Ru3(CO)n(NO)2, Ru3(N)(CO)n(NO), Ru3(N)2(CO)n, Ru3(N)(CO)n(NCO), Ru3(CO)n(NCO)(NO), Ru4(N)(CO)n(NO), Ru4(N)(CO)n(NCO), and Ru4(N)2(CO)n related to species observed experimentally in the chemistry of Ru3(CO)10(µ-NO)2 have been investigated using density functional theory. In all cases, the experimentally observed structures have been found to be low-energy structures. The low-energy trinuclear structures typically have a central strongly bent Ru-Ru-Ru chain with terminal CO groups and bridging nitrosyl, isocyanate, and/or nitride ligands across the end of the chain. The low-energy tetranuclear structures typically have a central Ru4N unit with terminal CO groups and a non-bonded pair of ruthenium atoms bridged by a nitrosyl or isocyanate group.

Automatic Differentiation for Explicitly Correlated MP2.

Automatic differentiation (AD) offers a route to achieve arbitrary-order derivatives of challenging wave function methods without the use of analytic gradients or response theory. Currently, AD has been predominantly used in methods where first- and/or second-order derivatives are available, but it has not been applied to methods lacking available derivatives. The most robust approximation of explicitly correlated MP2, MP2-F12/3C(FIX)+CABS, is one such method. By comparing the results of MP2-F12 computed with AD versus finite-differences, it is shown that (a) optimized geometries match to about 10-3 Å for bond lengths and a 10-6 degree for angles, and (b) dipole moments match to about 10-6 D. Hessians were observed to have poorer agreement with numerical results (10-5), which is attributed to deficiencies in AD implementations currently. However, it is notable that vibrational frequencies match within 10-2 cm-1. The use of AD also allowed the prediction of MP2-F12/3C(FIX)+CABS IR intensities for the first time.

Hydrogen bonding in 2,2,2-trifluoroethanol.

CONTEXT: 2,2,2-Trifluoroethanol (TFE) is known as a membrane mimetic solvent. The IR spectrum, 1H NMR spectrum, 13C NMR spin‒lattice relaxation times (T1), and nuclear Overhauser effect (NOE) data are consistent with extensive hydrogen bonding in TFE, but do not lead to structural features of the hydrogen bonding. Hence, DFT computations were carried out. The results predict the existence of a set of H-bonded dimers and trimers. The bond lengths and dihedral angles in these complexes are obtained, together with their dissociation energies. Computations were also performed for the geometry of the two conformers of the isolated monomer. The structure of one of the dimers consists of a 7-member cyclic fragment with a free CF3CH2 side chain. One set of the trimer structures involves the OH of a third monomer H-bonding to one of the F atoms in the CF3 group of the side chain of this dimer, thereby creating three trimer isomers. A fourth trimer cluster is formed from three monomers in which three OH∙∙∙O bonds create a cyclic fragment with three CF3CH2 side chains. The high dissociation energy (with respect to three monomers) indicates the high stability of the trimer complexes. The structural features of the trimer complexes resemble the structure of a conventional liquid crystal molecule and are postulated to resemble the latter in properties and function in solution, but at a much shorter timescale because of the noncovalent bonding. This hydrogen bonding phenomenon of TFE may be related to its function as a membrane memetic solvent. METHODS: Initially, IR and NMR spectroscopic methods were used. Standard procedures were followed. For the computations, a hybrid DFT method with empirical dispersion, ωB97X-D, was used. The basis set, 6-311++G**, is of triple-ζ quality, in which polarization functions and diffuse functions were added for all atoms.

Computationally Driven Discovery of a BCR-ABL1 Kinase Inhibitor with Activity in Multidrug-Resistant Chronic Myeloid Leukemia.

The permeability glycoprotein, encoded by the ABCB1 gene, is widely implicated in multidrug resistance (MDR), as it has been shown to reduce the intracellular concentration of most small molecule therapeutics, including the majority of the breakpoint cluster region Abelson proto-oncogene 1 (BCR-ABL1) kinase inhibitors used in the treatment of Philadelphia chromosome positive (Ph+) leukemias. With this in mind, we describe an integrated theoretical and experimental approach to shed light on substituent effects in the pendant anilino moiety of 4-anilinoquinazolines and 4-anilinoquinoline-3-carbonitrile-based kinase inhibitors and their influence on P-gp-mediated efflux. This analysis culminated in the identification of a hydroxylamine-bearing, dual cSRC/BCR-ABL1 kinase inhibitor 16a that exhibits a marked reduction in P-gp-mediated efflux ratio and potent activity in a Ph+ patient-derived cell line (K562) and an MDR-Ph+ patient-derived cell line (K562/Dox) overexpressing P-gp. Overall, we demonstrate that the P-gp-mediated efflux ratio can be minimized by computationally driven optimization of the molecular dipole and/or cpKa without recourse to intramolecular hydrogen bonds.

Comparative Study of Neutral and Cationic Sn2H2: Toward Laboratory Detection of the Cation.

Group 14 M2H2 isomers (M = Si, Ge, Sn, and Pb) have attracted interest due to their radically differing electronic structures from acetylene. To better understand the Sn-H interactions of the neutral and cationic Sn2H2 structures, we present the most rigorous study of these systems to date. CCSD(T)/cc-pwCVTZ harmonic frequencies are presented as the first predictions for the neutral and cationic species to date. CCSDT(Q)/CBS relative energies are reported using the focal point approach, confirming the butterfly isomer as the global minimum on the potential energy surface for both the neutral and cationic species. In all, there exist 7 minima and 15 transition states. NBO analysis is also performed to elucidate the changes in bond order going from neutral to cation across all isomers of Sn2H2. Our results provide insights into the important Sn-H interaction and provide guidance for future work that may detect Sn2H2+ in the laboratory for the first time.

Accurate and efficient open-source implementation of domain-based local pair natural orbital (DLPNO) coupled-cluster theory using a t1-transformed Hamiltonian.

We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol-1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol-1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms).

Ethynyl Radical Hydrogen Abstraction Energetics and Kinetics Utilizing High-Level Theory.

The ethynyl radical, C2H, is found in a variety of different environments ranging from interstellar space and planetary atmospheres to playing an important role in the combustion of various alkynes under fuel-rich conditions. Hydrogen-atom abstraction reactions are common for the ethynyl radical in these contrasting environments. In this study, the C2H + HX → C2H2 + X, where HX = HNCO, trans-HONO, cis-HONO, C2H4, and CH3OH, reactions have been investigated at rigorously high levels of theory, including CCSD(T)-F12a/cc-pVTZ-F12. For the stationary points thus located, much higher levels of theory have been used, with basis sets as large as aug-cc-pV5Z and methods up to CCSDT(Q), and core correlation was also included. These molecules were chosen because they can be found in either interstellar or combustion environments. Various additive energy corrections have been included to converge the relative enthalpies of the stationary points to subchemical accuracy (≤0.5 kcal mol-1). Barriers predicted here (2.19 kcal mol-1 for the HNCO reaction and 0.47 kcal mol-1 for C2H4) are significantly lower than previous predictions. Reliable kinetics were acquired over a wide range of temperatures (50-5000 K), which may be useful for future experimental studies of these reactions.

Binuclear Lantern Complexes of All First-Row 3d Metals Scandium to Zinc: Metal-Metal Bond Lengths and Bond Orders, with Measures of Intrinsic Metal-Metal Bond Strength.

The B3LYP and M06-L functionals with the cc-pVTZ basis set are used to study lantern-type binuclear complexes of all the first-row (3d block) metals scandium to zinc in various low-energy spin states, out of which the ground states are predicted. These complexes are studied as models using mostly the unsubstituted formamidinate ligand. For each metal, metal-metal (MM) bond lengths are related to the formal MM bond orders (zero to five), derived by MO analysis and by electron counting. The predicted ground-state spin multiplicities and MM bond lengths of the model complexes generally agree fairly well with available experimental results on substituted analogues. Finally, values of the formal shortness ratio and Wiberg index for the MM bonds in all of these complexes in all spin states studied are categorized into ranges according to the MM bond orders 0 to 5 in steps of 0.5. The trends shown validate their use in estimating intrinsic metal-metal bond strength regardless of the metal.

MolSym: A Python package for handling symmetry in molecular quantum chemistry.

A consideration of the point group symmetry of molecules is often advantageous from a computational efficiency standpoint and sometimes necessary for the correct treatment of chemical physics problems. Many modern electronic structure software packages include a treatment of symmetry, but these are sometimes incomplete or unusable outside of that program's environment. Therefore, we have developed the MolSym package for handling molecular symmetry and its associated functionalities to provide a platform for including symmetry in the implementation and development of other methods. Features include point group detection, molecule symmetrization, arbitrary generation of symmetry element sets and character tables, and symmetry adapted linear combinations of real spherical harmonic basis functions, Cartesian displacement coordinates, and internal coordinates. We present some of the advantages of using molecular symmetry as achieved by MolSym, particularly with respect to Hartree-Fock theory, and the reduction of finite difference displacements in gradient/Hessian computations. This package is designed to be easily integrated into other software development efforts and may be extended to further symmetry applications.

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.

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.

Biomimetic Frustrated Lewis Pair Catalysts for Hydrogenation of CO to Methanol at Low Temperatures.

The industrial production of methanol through CO hydrogenation using the Cu/ZnO/Al2O3 catalyst requires harsh conditions, and the development of new catalysts with low operating temperatures is highly desirable. In this study, organic biomimetic FLP catalysts with good tolerance to CO poison are theoretically designed. The base-free catalytic reaction contains the 1,1-addition of CO into a formic acid intermediate and the hydrogenation of the formic acid intermediate into methanol. Low-energy spans (25.6, 22.1, and 20.6 kcal/mol) are achieved, indicating that CO can be hydrogenated into methanol at low temperatures. The new extended aromatization-dearomatization effect involving multiple rings is proposed to effectively facilitate the rate-determining CO 1,1-addition step, and a new CO activation model is proposed for organic catalysts.

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.

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.

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.

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.

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.

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.

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.