Helium Droplet Infrared Spectroscopy of the Butyl Radicals

Butyl radicals ( n-, s-, i-,} and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo- tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700-3125 cm-1, allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800-3000 cm-1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on CCSD(T) quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

Demarcating Noncovalent and Covalent Bond Territories: Imine-CO2 Complexes and Cooperative CO2 Capture

Chemical bond territory is rich with covalently bonded molecules wherein a strong bond is formed by equal or unequal sharing of a quantum of electrons. The noncovalent version of the bonding scenarios expands the chemical bonding territory to a weak domain wherein the interplay of electrostatic and π-effects, dipole-dipole, dipole-induced dipole, and induced dipole-induced dipole interactions, and hydrophobic effects occur. Here we study both the covalent and noncovalent interactive behavior of cyclic and acyclic imine-based functional molecules (XN) with CO₂. All parent XN systems preferred the formation of noncovalent (nc) complex XN···CO₂, while more saturated such systems (XN') produced both nc and covalent (c) complexes XN'⁺-(CO₂)^-. In all such cases, crossover from an nc to c complex is clearly demarcated with the identification of a transition state (ts). The complexes XN'···CO₂ and XN'⁺-(CO₂)^- are bond stretch isomers, and they define the weak and strong bonding territories, respectively, while the ts appears as the demarcation point of the two territories. Cluster formation of XN with CO₂ reinforces the interaction between them, and all become covalent clusters of general formula (XN⁺-(CO₂)^-)_n. The positive cooperativity associated with the NH···OC hydrogen bond formation between any two XN'⁺-(CO₂)^- units strengthened the N-C coordinate covalent bond and led to massive stabilization of the cluster. For instance, the stabilizing interaction between the XN unit with CO₂ is increased from 2-7 kcal/mol range in a monomer complex to 14-31 kcal/mol range for the octamer cluster (XN'⁺-(CO₂)^-)₈. The cooperativity effect compensates for the large reduction in the entropy of cluster formation. Several imine systems showed the exergonic formation of the cluster and are predicted as potential candidates for CO₂ capture and conversion.

Mini-Review on Structure-Reactivity Relationship for Aromatic Molecules: Recent Advances

Recent advances in quantifying nucleophilic reactivities in chemical reactions and intermolecular interactions of aromatic molecules are reviewed. This survey covers experimental (IR frequency shifts induced by hydrogen bonding) and theoretical (modeling of potential energy surfaces, atomic charges, molecular electrostatic potential) approaches in characterizing chemical reactivity. Recent advances in software developments assisting the evaluation of the reactive sites for electrophilic aromatic substitution are briefly discussed.

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

Herein, we report a comprehensive study of CO₂ 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 NaO^tBu. 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 CO₂ molecule. Based on these discoveries, we designed a pyridine-based Mn pincer catalyst system, which could achieve CO₂ hydroboration in low-temperature and base-free conditions through a metal-ligand cooperation mechanism.

Substituent, Solvent, and Dispersion Effects on the Zwitterionic Character and Dimerization Thermochemistry of the Group 6 Fulvene Metal Tricarbonyl Complexes

Dimerizations of fulvene metal tricarbonyl complexes of the type (C₅H₄CRR')M(CO)₃ (R, R' = MeO, Me, H; M = Cr, Mo, W) to form a metal-metal bond and a new carbon-carbon bond, thereby giving binuclear cyclopentadienyl metal carbonyl derivatives, are predicted to be thermochemically favored but to have significant activation energies ranging from ΔE = 19 to 42 kcal/mol. However, the introduction of dimethylamino but not methoxy substituents onto the exocyclic carbon atom changes the situation drastically so that the monomers [C₅H₄CH(NMe₂)]M(CO)₃ and [C₅H₄C(NMe₂)₂]M(CO)₃ become strongly thermochemically favored, lying ΔE = 43 kcal/mol (M = W) to 63 kcal/mol (M = Cr) below their corresponding dimers. In such dimethylamino-substituted (fulvene)M(CO)₃ derivatives, the M-C distance to the exocyclic fulvene carbon is lengthened beyond the bonding distance to give a zwitterionic structure with a pentahapto fulvene ligand. Such M-C distances in (fulvene)M(CO)₃ complexes, which have preferred zwitterionic structures, increase with increasing solvent polarity (i.e., dielectric constant) until a saturation point is reached.

Fermi.jl: A Modern Design for Quantum Chemistry

Approximating molecular wave functions involves heavy numerical effort; therefore, codes for such tasks are written completely or partially in efficient languages such as C, C++, and Fortran. While these tools are dominant throughout quantum chemistry packages, the efficient development of new methods is often hindered by the complexity associated with code development. In order to ameliorate this scenario, some software packages take a dual approach where a simpler, higher-level language, such as Python, substitutes the traditional ones wherever performance is not critical. Julia is a novel, dynamically typed, programming language that aims to solve this two-language problem. It gained attention because of its modern and intuitive design, while still being highly optimized to compete with "low-level" languages. Recently, some chemistry-related projects have emerged exploring the capabilities of Julia. Herein, we introduce the quantum chemistry package Fermi.jl, which contains the first implementations of post-Hartree-Fock methods written in Julia. Its design makes use of many Julia core features, including multiple dispatch, metaprogramming, and interactive usage. Fermi.jl is a modular package, where new methods and implementations can be easily added to the existing code. Furthermore, it is designed to maximize code reusability by relying on general functions with specialized methods for particular cases. The feasibility of the project is explored through evaluating the performance of popular ab initio methods. It is our hope that this project motivates the usage of Julia within the community and brings new contributions into Fermi.jl.

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

The Br + (H2O)3 → HBr + (H2O)2OH reaction has been investigated using the CCSD(T) method with the basis sets as large as cc-pVQZ(-PP). The Br + (H2O)3 reaction is also compared with related Br + H2O/(H2O)2 and F/Cl + (H2O)3 reactions.

Cumulants as the variables of density cumulant theory: A path to Hermitian triples

We study the combination of orbital-optimized density cumulant theory and a new parameterization of reduced density matrices in which the variables are the particle-hole cumulant elements. We call this combination OλDCT. We find that this new Ansatz solves problems identified in the previous unitary coupled cluster Ansatz for density cumulant theory: the theory is now free of near-zero denominators between occupied and virtual blocks, can correctly describe the dissociation of H₂, and is rigorously size-extensive. In addition, the new Ansatz has fewer terms than the previous unitary Ansatz, and the optimal orbitals delivered by the exact theory are the natural orbitals. Numerical studies on systems amenable to full configuration interaction show that the amplitudes from the previous ODC-12 method approximate the exact amplitudes predicted by this Ansatz. Studies on equilibrium properties of diatomic molecules show that even with the new Ansatz, it is necessary to include triples to improve the accuracy of the method compared to orbital-optimized linearized coupled cluster doubles. With a simple iterative triples correction, OλDCT outperforms other orbital-optimized methods truncated at comparable levels in the amplitudes, as well as coupled cluster single and doubles with perturbative triples [CCSD(T)]. By adding four more terms to the cumulant parameterization, OλDCT outperforms CCSDT while having the same O(V⁵O³) scaling.

Substituent Effects on Aluminyl Anions and Derived Systems: A High-Level Theory

Aluminyl anions are low-valent aluminum species bearing a lone pair of electrons and a negative charge. These systems have drawn recent synthetic interest for their nucleophilic nature, allowing for the activation of σ-bonds, and have been proposed as a pathway to hydrogen energy storage. In this research, we provide high-level ab initio geometries and energies for both the simplest aluminyl anion (AlH₂^-) and several substituted derivatives. Geometries are reported using the gold-standard CCSD(T)/aug-cc-pV(T+d)Z level of theory. Energies were extrapolated to the complete basis set limit through the focal point approach, utilizing coupled-cluster methods through perturbative quadruples and basis sets up to five-ζ quality. Geometries were rationalized using electrostatic, steric, and orbital donation effects. The donation from substituents to Al is accompanied by back-donation effects, a property traditionally thought of in transition-metal systems. Stereoelectronic effects through the secondary orbital interaction play a fundamental role in stabilizing these low-valent aluminum compounds and would likely also affect the feasibility of their use within several industrial applications. The energetic analysis of the formation of each substituted anion is rationalized as the result of three energetic schemes. The effectiveness of these schemes for determining the relative formation energies is discussed.

Potential energy profile for the Cl + (H2O)3 → HCl + (H2O)2OH reaction. A CCSD(T) study

Four different reaction pathways are initially located for the reaction of Cl atom plus water trimer Cl + (H₂O)₃ → HCl + (H₂O)₂OH using a standard DFT method. As found for the analogous fluorine reaction, the geometrical and energetic results for the four chlorine pathways are closely related. However, the energetics for the Cl reaction are very different from those for fluorine. In the present paper, we investigate the lowest-energy chlorine pathway using the "gold standard" CCSD(T) method in conjunction with correlation-consistent basis sets up to cc-pVQZ. Structurally, the stationary points for the water trimer reaction Cl + (H₂O)₃ may be compared to those for the water monomer reaction Cl + H₂O and water dimer reaction Cl + (H₂O)₂. Based on the CCSD(T) energies, the title reaction is endothermic by 19.3 kcal mol^-1, with a classical barrier height of 16.7 kcal mol^-1 between the reactants and the exit complex. There is no barrier for the reverse reaction. The Cl⋯(H₂O)₃ entrance complex lies 5.3 kcal mol^-1 below the separated reactants. The HCl⋯(H₂O)₂OH exit complex is bound by 8.6 kcal mol^-1 relative to the separated products. The Cl + (H₂O)₃ reaction is somewhat similar to the analogous Cl + (H₂O)₂ reaction, but qualitatively different from the Cl + H₂O reaction. It is reasonable to expect that the reactions between the chlorine atom and larger water clusters may be similar to the Cl + (H₂O)₃ reaction. The potential energy profile for the Cl + (H₂O)₃ reaction is radically different from that for the valence isoelectronic F + (H₂O)₃ system, which may be related to the different bond energies between HCl and HF.

Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

Kinetic Stability of Pentazole

The utility of high energy density materials (HEDMs) comes from their thermodynamic properties which arise from specific structural features that contribute to energy storage. Studies of such structural features seek to increase our understanding of these energy storage mechanisms in order to further enhance their properties. High-nitrogen-containing HEDMs are of particular interest because they are less toxic than traditional HEDMs. Pentazole is the largest of the nitrogen rings which has been synthesized and considered for an HEDM; however, few experimental studies exist due to the difficulty involved in the synthesis, and most previous theoretical studies employed composite methods where lower level geometries were used with higher level methods. Here, the decomposition reaction of pentazole is studied. Geometries, fundamental frequencies, and energies for each of the stationary points of the decomposition pathway are computed using ab initio methods up to CCSDT(Q). Decomposition rates are calculated over a range of temperatures using canonical transition state theory in order to determine the kinetic stability of pentazole. Based on the present results, it would be difficult for pentazole to act as an HEDM, requiring temperatures close to 200 K to achieve a suitable level of stability.

Group 15 and 16 Nitrene-Like Pnictinidenes

Pnictinidenes are an increasingly relevant species in main group chemistry and generally exhibit proclivity for the triplet electronic ground state. However, the elusive singlet electronic states are often desired for chemical applications. We predict the singlet-triplet energy differences (ΔE_ST =E_Singlet -E_Triplet ) of simple group 15 and 16 substituted pnictinidenes (Pn-R; Pn=P, As, Sb, or Bi) with highly reliable focal-point analyses targeting the CCSDTQ/CBS level of theory. The only cases we predict to have favorable singlet states are P-PH₂ (-3.2 kcal mol^-1 ) and P-NH₂ (-0.2 kcal mol^-1 ). ΔE_ST trends are discussed in light of the geometric predictions as well as qualitative natural bond order analysis to elucidate some of the important electronic structure features. Our work provides a rigorous benchmark for the ΔE_ST of fundamental Pn-R moieties and provides a firm foundation for the continued study of heavier pnictinidenes.

Carbene-Stabilized Dithiolene (L0 ) Zwitterions

A series of reactions between Lewis bases and an imidazole-based dithione dimer (1) has been investigated. Both cyclic(alkyl)(amino)carbene (CAAC) (2) and N-heterocyclic carbene (NHC) (4), in addition to N-heterocyclic silylene (NHSi) (6), demonstrate the capability to cleave the sulphur-sulphur bonds in 1, giving carbene-stabilized dithiolene (L⁰ ) zwitterions (3 and 5) and a spirocyclic silicon-dithiolene compound (7), respectively. The bonding nature of 3, 5, and 7 are probed by both experimental and theoretical methods.

Contrasting the Mechanism of H2 Activation by Monomeric and Potassium-Stabilized Dimeric AlI Complexes: Do Potassium Atoms Exert any Cooperative Effect?

Aluminyl anions are low-valent, anionic, and carbenoid aluminum species commonly found stabilized with potassium cations from the reaction of Al-halogen precursors and alkali compounds. These systems are very reactive toward the activation of σ-bonds and in reactions with electrophiles. Various research groups have detected that the potassium atoms play a stabilization role via electrostatic and cation ⋯ π interactions with nearby (aromatic)-carbocyclic rings from both the ligand and from the reaction with unsaturated substrates. Since stabilizing K⋯H bonds are witnessed in the activation of this class of molecules, we aim to unveil the role of these metals in the activation of the smaller and less polarizable H₂ molecule, together with a comprehensive characterization of the reaction mechanism. In this work, the activation of H₂ utilizing a NON-xanthene-Al dimer, [K{Al(NON)}]₂ (D) and monomeric, [Al(NON)]^- (M) complexes are studied using density functional theory and high-level coupled-cluster theory to reveal the potential role of K⁺ atoms during the activation of this gas. Furthermore, we aim to reveal whether D is more reactive than M (or vice versa), or if complicity between the two monomer units exits within the D complex toward the activation of H₂ . The results suggest that activation energies using the dimeric and monomeric complexes were found to be very close (around 33 kcal mol^-1 ). However, a partition of activation energies unveiled that the nature of the energy barriers for the monomeric and dimeric complexes are inherently different. The former is dominated by a more substantial distortion of the reactants (and increased interaction energies between them). Interestingly, during the oxidative addition, the distortion of the Al complex is minimal, while H₂ distorts the most, usually over 0.77 Δ E d i s t ≠ . Overall, it is found here that electrostatic and induction energies between the complexes and H₂ are the main stabilizing components up to the respective transition states. The results suggest that the K⁺ atoms act as stabilizers of the dimeric structure, and their cooperative role on the reaction mechanism may be negligible, acting as mere spectators in the activation of H₂ . Cooperation between the two monomers in D is lacking, and therefore the subsequent activation of H₂ is wholly disengaged.

Lantern-Type Divanadium Complexes with Bridging Ligands: Short Metal-Metal Bonds with High Multiple Bond Orders

Vanadium forms binuclear complexes with a variety of ligands often containing V≡V triple bonds. Many tetragonal divanadium paddlewheel complexes with bridging bidentate ligands have been experimentally characterized. This research exhaustively treats model tetragonal, trigonal, and digonal paddlewheel-type divanadium complexes V₂ L_x (L=formamidinate, guanidinate, and carboxylate; x=2, 3, 4), each in the three lowest-energy spin states. The V-V formal bond orders are obtained from metal-metal MO diagrams for representative structures. A number of short V-V multiple bonds of order 3, 3.5, and 4 are found in these model complexes. The short V≡V triple bonds and singlet ground state predicted here for the model tetragonal complexes correspond well with the limited experimental results for the series of known tetragonal paddlewheels. Digonal divanadium lanterns with very short V-V quadruple bonds are predicted as interesting synthetic targets. The V-V bond distances are categorized into distinct ranges according to the formal bond order values from 0.5 to 4. These bond length ranges are compared with the ranges compiled for other divanadium complexes including carbonyl complexes.

Catalyzed reaction of isocyanates (RNCO) with water

The reactions between substituted isocyanates (RNCO) and other small molecules (e.g. water, alcohols, and amines) are of significant industrial importance, particularly for the development of novel polyurethanes and other useful polymers. We present very high-level ab initio computations on the HNCO + H₂O reaction, with results targeting the CCSDT(Q)/CBS//CCSD(T)/cc-pVQZ level of theory. Our results affirm that hydrolysis can occur across both the N[double bond, length as m-dash]C and C[double bond, length as m-dash]O bonds of HNCO via concerted mechanisms to form carbamate or imidic acid with ΔH_0K barrier heights of 38.5 and 47.5 kcal mol^-1. A total of 24 substituted RNCO + H₂O reactions were studied. Geometries obtained with a composite method and refined with CCSD(T)/CBS single point energies determine that substituted RNCO species have a significant influence on these barrier heights, with an extreme case like fluorine lowering both barriers by close to 15 kcal mol^-1 and most common alkyl substituents lowering both by approximately 3 kcal mol^-1. Natural Bond Orbital (NBO) analysis provides evidence that the predicted barrier heights are strongly associated with the occupation of the in-plane C-O* orbital of the RNCO reactant. Key autocatalytic mechanisms are considered in the presence of excess water and RNCO species. Additional waters (one or two) are predicted to lower both barriers significantly at the CCSD(T)/aug-cc-pV(T+d)Z level of theory with strongly electron withdrawing RNCO substituents also increasing these effects, similar to the uncatalyzed case. The 298 K Gibbs energies are only marginally lowered by a second catalyst water molecule, indicating that the decreasing ΔH_0K barriers are offset by loss of translational entropy with more than one catalyst water. Two-step 2RNCO + H₂O mechanisms are characterized for the formation of carbamate and imidic acid. The second step of these two pathways exhibits the largest barrier and presents no clear pattern with respect to substituent choice. Our results indicate that an additional RNCO molecule might catalyze imidic acid formation but have less influence on the efficiency of carbamate formation. We expect that these results lay a firm foundation for the experimental study of substituted isocyanates and their relationship to the energetic pathways of related systems.