Benchmark ab Initio Energy Profiles for the Gas-Phase SN2 Reactions Y- + CH3X → CH3Y + X- (X,Y = F,Cl,Br). Validation of Hybrid DFT Methods

On Tautomerism and Amphoterism: An In-Depth Structural and Physicochemical Characterization of Ammeline and Some of Its Salts

Ammeline is a simple, readily available, molecular compound, which has been known for nearly 200 years. Despite that, no proper structural characterization of ammeline has been conducted so far. For this reason, the prevalent tautomeric form of ammeline in the solid remained unknown to this date. In the course of this study, its crystal structure was finally established by single-crystal X-ray diffraction. In this structure, ammeline is exclusively found as its 4,6-diamino-1,3,5-triazin-2(1H)-one tautomer and adopts layered structure with an exceptionally high hydrogen bond density. Ammeline shows an interesting amphoteric behavior. Therefore, the synthesis and structural characterization of some of its salts were carried out to investigate the influence of the protonation degree on its molecular structure. In particular, the crystal structure of silver ammelinate monohydrate was solved as the first reported structure containing deprotonated ammeline. Moreover, the crystal structures of three different modifications of ammelinium perchlorate were elucidated and the transformation conditions between them were studied. Lastly, the crystal structure of ammelinediium diperchlorate monohydrate, containing unprecedented doubly protonated ammeline, was determined. The products' thermal behavior was studied by differential thermal analysis and thermogravimetric analysis. The perchlorate salts were additionally examined for their potential as insensitive high-energy-density materials.

Stereo-electronic factors influencing the stability of hydroperoxyalkyl radicals: transferability of chemical trends across hydrocarbons and ab initio methods

The hydroperoxyalkyl radicals (˙QOOH) are known to play a significant role in combustion and tropospheric processes, yet their direct spectroscopic detection remains challenging. In this study, we investigate molecular stereo-electronic effects influencing the kinetic and thermodynamic stability of a ˙QOOH along its formation path from the precursor, alkylperoxyl radical (ROO˙), and the depletion path resulting in the formation of cyclic ether + ˙OH. We focus on reactive intermediates encountered in the oxidation of acyclic hydrocarbon radicals: ethyl, isopropyl, isobutyl, tert-butyl, neopentyl, and their alicyclic counterparts: cyclohexyl, cyclohexenyl, and cyclohexadienyl. We report reaction energies and barriers calculated with the highly accurate method Weizmann-1 (W1) for the channels: ROO˙ ⇌ ˙QOOH, ROO˙ ⇌ alkene + ˙OOH, ˙QOOH ⇌ alkene + ˙OOH, and ˙QOOH ⇌ cyclic ether + ˙OH. Using W1 results as a reference, we have systematically benchmarked the accuracy of popular density functional theory (DFT), composite thermochemistry methods, and an explicitly correlated coupled-cluster method. We ascertain inductive, resonance, and steric effects on the overall stability of ˙QOOH and computationally investigate the possibility of forming more stable species. With new reactions as test cases, we probe the capacity of various ab initio methods to yield quantitative insights on the elementary steps of combustion.

How do density functionals affect the Hirshfeld atom refinement?

In this work, the effect of mixing different amounts of Hartree-Fock (HF) exchange with hybrid density functionals applied to the Hirshfeld atom refinement (HAR) of urea and oxalic acid dihydrate is explored. Together, the influence of using different basis sets, methods (including MP2 and HF) and cluster sizes (to model bulk effects) is studied. The results show that changing the amount of HF exchange, no matter the level of theory, has an impact almost exclusively on the H atom refinement parameters. Contrary to pure quantum mechanical calculations where good geometries are obtained with intermediate HF exchange mixtures, in the HAR the best match with neutron diffraction reference values is not necessarily found for these admixtures. While the non-hydrogen covalent bond lengths are insensitive to the combination of method or basis set employed, the X-H bond lengths always increase proportionally to the HF exchange for the analysed systems. This outcome is opposite to what is normally observed from geometry optimisations, i.e., shorter bonds are obtained with greater HF exchange. Additionally, the thermal ellipsoids tend to shrink with larger HF exchange, especially for the H atoms involved in strong hydrogen bonding. Thus, it may be the case that the development of density functionals or basis sets suitable for quantum crystallography should take a different path than those fitted for quantum chemistry calculations.

Hydrogen Bonding versus Halogen Bonding: Spectroscopic Investigation of Gas-Phase Complexes Involving Bromide and Chloromethanes

Hydrogen bonding and halogen bonding are important non-covalent interactions that are known to occur in large molecular systems, such as in proteins and crystal structures. Although these interactions are important on a large scale, studying hydrogen and halogen bonding in small, gas-phase chemical species allows for the binding strengths to be determined and compared at a fundamental level. In this study, anion photoelectron spectra are presented for the gas-phase complexes involving bromide and the four chloromethanes, CH₃ Cl, CH₂ Cl₂ , CHCl₃ , and CCl₄ . The stabilisation energy and electron binding energy associated with each complex are determined experimentally, and the spectra are rationalised by high-level CCSD(T) calculations to determine the non-covalent interactions binding the complexes. These calculations involve nucleophilic bromide and electrophilic bromine interactions with chloromethanes, where the binding motifs, dissociation energies and vertical detachment energies are compared in terms of hydrogen bonding and halogen bonding.

Spectroscopic Study of the Br- +CH3 I→I- +CH3 Br SN 2 Reaction

Mass spectrometry and anion photoelectron spectroscopy have been used to study the gas-phaseS N 2 ${{{\rm S}}_{{\rm N}}2}$ reaction involvingB r - ${{{\rm B}{\rm r}}^{-}}$ andC H 3 I ${{{\rm C}{\rm H}}_{3}{\rm I}}$ . The anion photoelectron spectra associated with the reaction intermediates of thisS N 2 ${{{\rm S}}_{{\rm N}}2}$ reaction are presented. High-level CCSD(T) calculations have been utilised to investigate the reaction intermediates that may form as a result of theS N 2 ${{{\rm S}}_{{\rm N}}2}$ reaction along various different reaction pathways, including back-side attack and front-side attack. In addition, simulated vertical detachment energies of each reaction intermediate have been calculated to rationalise the photoelectron spectra.

Accurate Quantum Chemical Prediction of Gas-Phase Anion Binding Affinities and Their Structure-Binding Relationships

This paper systematically examines the performance of contemporary wavefunction and density functional theory methods to identify robust and cost-efficient methods for predicting gas-phase anion binding energies. This includes the local coupled cluster LNO-CCSD(T) and DLPNO-CCSD(T), as well as double-hybrid DSD-PBEP86-D3(BJ) and various hybrid functionals M06-2X, B3LYP-D3(BJ), ωB97M-V, and ωB97X-V. The focus is on dual-hydrogen-bonding anion receptors that are commonly found in supramolecular chemistry and organocatalysis, namely, (thio)ureas, deltamides, (thio)squaramides, and croconamides as well as the yet-to-be-explored rhodizonamides. Of the methods examined, M06-2X emerged as the overall best performing method as the other functionals including DSD-PBEP86-D3(BJ) and the local coupled cluster DLPNO-CCSD(T) method displayed systematic errors that increase with the degree of carbonylation of the receptors. Hybrid ONIOM models that employed semiempirical methods (PM7, GFN1-xTB, and GFN2-xTB) and "threefold"-corrected small-basis set potentials (HF-3c, B97-3c, and PBEh-3c) were explored, and the best models resulted in 50- to 500-fold reduction in CPU time compared to W1-local. These calculations provide important insight into the structure-binding relationships where there is a direct correlation between Brønsted acidity and anion binding affinity, though the strength of the correlation also depends on other factors such as hydrogen-bonding geometry and the geometrical distortion that the receptor needs to undergo to bind the anion.

The SN2 reaction and its relationship with the Walden inversion, the Finkelstein and Menshutkin reactions together with theoretical calculations for the Finkelstein reaction

This communication gives an overview of the relationships between four reactions that although related were not always perceived as such: SN2, Walden, Finkelstein, and Menshutkin. Binary interactions (SN2 & Walden, SN2 & Menshutkin, SN2 & Finkelstein, Walden & Menshutkin, Walden & Finkelstein, Menshutkin & Finkelstein) were reported. Carbon, silicon, nitrogen, and phosphorus as central atoms and fluorides, chlorides, bromides, and iodides as lateral atoms were considered. Theoretical calculations provide Gibbs free energies that were analyzed with linear models to obtain the halide contributions. The M06-2x DFT computational method and the 6-311++G(d,p) basis set have been used for all atoms except for iodine where the effective core potential def2-TZVP basis set was used. Concerning the central atom pairs, carbon/silicon vs. nitrogen/phosphorus, we reported here for the first time that the effect of valence expansion was known for Si but not for P. Concerning the lateral halogen atoms, some empirical models including the interaction between F and I as entering and leaving groups explain the Gibbs free energies.

TSNet: predicting transition state structures with tensor field networks and transfer learning

Transition states are among the most important molecular structures in chemistry, critical to a variety of fields such as reaction kinetics, catalyst design, and the study of protein function. However, transition states are very unstable, typically only existing on the order of femtoseconds. The transient nature of these structures makes them incredibly difficult to study, thus chemists often turn to simulation. Unfortunately, computer simulation of transition states is also challenging, as they are first-order saddle points on highly dimensional mathematical surfaces. Locating these points is resource intensive and unreliable, resulting in methods which can take very long to converge. Machine learning, a relatively novel class of algorithm, has led to radical changes in several fields of computation, including computer vision and natural language processing due to its aptitude for highly accurate function approximation. While machine learning has been widely adopted throughout computational chemistry as a lightweight alternative to costly quantum mechanical calculations, little research has been pursued which utilizes machine learning for transition state structure optimization. In this paper TSNet is presented, a new end-to-end Siamese message-passing neural network based on tensor field networks shown to be capable of predicting transition state geometries. Also presented is a small dataset of S_N2 reactions which includes transition state structures - the first of its kind built specifically for machine learning. Finally, transfer learning, a low data remedial technique, is explored to understand the viability of pretraining TSNet on widely available chemical data may provide better starting points during training, faster convergence, and lower loss values. Aspects of the new dataset and model shall be discussed in detail, along with motivations and general outlook on the future of machine learning-based transition state prediction.

Thermal activation of methane by MgO+: temperature dependent kinetics, reactive molecular dynamics simulations and statistical modeling

The kinetics methane activation (MgO⁺ + CH₄) was studied experimentally and computationally by running and analyzing reactive atomistic simulations.

Nucleophilic Influences and Origin of the SN 2 Allylic Effect

The potential energy surfaces for the S_N 2 reactions of allyl and propyl chlorides with 21 anionic and neutral nucleophiles was studied by using ωB97X-D/6-311++G(3df,2pd) computations. The "allylic effect" on S_N 2 barriers was observed for all reactions, and compared with propyl substrates, the energy barriers differed by -0.2 to -4.5 kcal mol^-1 in the gas phase. Strong correlations of the S_N 2 net activation barriers with cation affinities, proton affinities, and electrostatic potentials at nuclei demonstrated the powerful influence of electrostatic interactions on these reactions. For the reactions of anionic (but not neutral) nucleophiles with allyl chloride, some of the incoming negative charge (0.2-18 %) migrated into the carbon chains, which would provide secondary stabilization of the S_N 2 transition states. Activation strain analysis provided additional insight into the allylic effect by showing that the energy of geometric distortion for the reactants to reach the S_N 2 transition state was smaller for each allylic reaction than for its propyl analogue. In many cases, the interaction energies between the substrate and nucleophile in this analysis were more favorable for propyl chloride reactions, but this compensation did not overcome the predominant strain energy effect.

Reaction mechanism of the dengue virus serine protease: a QM/MM study

The dengue virus (DENV) is the causative agent of the viral infection dengue fever. In spite of all the efforts made to prevent the spread of the disease, once it is contracted, there is no specific treatment for dengue and the WHO guidelines are limited to rest and symptomatic treatment. In its reproductive cycle, DENV utilizes the NS2B-NS3pro, a serine protease, to cleave the viral polyprotein into its constituents. This enzyme is essential for the virus lifecycle, and presents an attractive target for the development of specific dengue treatments. Here we used a hybrid Quantum Mechanics and Molecular Mechanics (QM/MM) Molecular Dynamics approach and Umbrella Sampling to study the first step (acylation) of the reaction catalyzed by NS2B-NS3pro, using the Pairwise Distance Directed Gaussian PM3 (PDDG/PM3) semi-empirical Hamiltonian for the QM subsystem, and Amber ff99SB for the MM subsystem. Our results indicate that the nucleophilic attack on the substrate by Ser135 occurs in a stepwise manner, in which a proton transfer to His51 first activates Ser135, which only later attacks the substrate. The rate-determining step is the Ser135 activation, with a barrier of 24.1 kcal mol^-1. Water molecules completing the oxyanion hole stabilize the negative charge formed on the carbonyl oxygen of the substrate. The final step in the process is a proton transfer from His51 to the substrate's nitrogen, which happens with a lower barrier of 5.1 kcal mol^-1, and leads directly to the breakage of the peptide bond.

Computational investigation of LiF containing hypersalts

This study explores the design of possible hypersalts starting from the hyperhalogen Li₃F₄ plus a Li atom and the hyperalkali Li₄F₃ plus a F atom.

Ab initio analysis on potential superbases of several hyperlithiated species: Li3F2O and Li3F2OHn (n = 1, 2)

A systematic investigation on the potential basicity of the novel hyperlithiated species Li₃F₂O and Li₃F₂(OH)_n (n = 1, 2) based upon the superalkali cluster Li₃F₂ was conducted using high-level ab initio techniques.

Dynamics and Novel Mechanisms of SN2 Reactions on ab Initio Analytical Potential Energy Surfaces

We describe a novel theoretical approach to the bimolecular nucleophilic substitution (S_N2) reactions that is based on analytical potential energy surfaces (PESs) obtained by fitting a few tens of thousands high-level ab initio energy points. These PESs allow computing millions of quasi-classical trajectories thereby providing unprecedented statistical accuracy for S_N2 reactions, as well as performing high-dimensional quantum dynamics computations. We developed full-dimensional ab initio PESs for the F^- + CH₃Y [Y = F, Cl, I] systems, which describe the direct and indirect, complex-forming Walden-inversion, the frontside attack, and the new double-inversion pathways as well as the proton-transfer channels. Reaction dynamics simulations on the new PESs revealed (a) a novel double-inversion S_N2 mechanism, (b) frontside complex formation,

Assessing How Correlated Molecular Orbital Calculations Can Perform versus Kohn-Sham DFT: Barrier Heights/Isomerizations

To assess the title issue, 38 hydrogen transfer barrier heights and 38 non-hydrogen transfer barrier heights/isomerizations extracted from extensive databases have been considered, in addition to 4 2 p-isomerization reactions and 6 others for large organic molecules. All Kohn-Sham DFT calculations have employed the popular M06-2X functional, whereas the correlated molecular orbital (MO)-based ones are from single-reference MP2 and CCSD(T) methods. They have all utilized the same basis sets, with raw MO energies subsequently extrapolated to the complete basis set limit without additional cost. MP2 calculations are found to be as cost-effective as DFT ones and often slightly more, while showing a satisfactory accuracy when compared with the reference data. Although the focus is on barrier heights, the results may bear broader implications, in that one may see successes and difficulties of DFT when compared with traditional MO theories for the same data.

Revisiting the Dielectric Constant Effect on the Nucleophile and Leaving Group of Prototypical Backside SN2 Reactions: A Reaction Force and Atomic Contribution Analysis

The solvent effect on the nucleophile and leaving group atoms of the prototypical F^- + CH₃Cl → CH₃F + Cl^- backside bimolecular nucleophilic substitution reaction (S_N2) is analyzed employing the reaction force and the atomic contributions methods on the intrinsic reaction coordinate (IRC). Solvent effects were accounted for using the polarizable continuum solvent model. Calculations were performed employing 11 dielectric constants, ε, ranging from 1.0 to 78.5, to cover a wide spectrum of solvents. The reaction force data reveal that the solvent mainly influences the region of the IRC preceding the energy barrier, where the structural rearrangement to reach the transition state occurs. A detailed analysis of the atomic role in the reaction as a function of ε reveals that the nucleophile and the carbon atom are the ones that contribute the most to the energy barrier. In addition, we investigated the effect of the choice of nucleophile and leaving group on the ΔE₀ and ΔE^‡ of Y^- + CH₃X → YCH₃ + X^- (X, Y = F, Cl, Br, I) in aqueous solution. Our analysis allowed us to find relationships between the atomic contributions to the activation energy and leaving group ability and nucleophilicity.

Polysulfide-1-oxides react with peroxyl radicals as quickly as hindered phenolic antioxidants and do so by a surprising concerted homolytic substitution

Polysulfides are important additives to a wide variety of industrial and consumer products and figure prominently in the chemistry and biology of garlic and related medicinal plants. Although their antioxidant activity in biological contexts has received only recent attention, they have long been ascribed 'secondary antioxidant' activity in the chemical industry, where they are believed to react with the hydroperoxide products of autoxidation to slow the auto-initiation of new autoxidative chain reactions. Herein we demonstrate that the initial products of trisulfide oxidation, trisulfide-1-oxides, are surprisingly reactive 'primary antioxidants', which slow autoxidation by trapping chain-carrying peroxyl radicals. In fact, they do so with rate constants (k = 1-2 × 104 M-1 s-1 at 37 °C) that are indistinguishable from those of the most common primary antioxidants, i.e. hindered phenols, such as BHT. Experimental and computational studies demonstrate that the reaction occurs by a concerted bimolecular homolytic substitution (SH2), liberating a perthiyl radical - which is ca. 16 kcal mol-1 more stable than a peroxyl radical. Interestingly, the (electrophilic) peroxyl radical nominally reacts as a nucleophile - attacking the of the trisulfide-1-oxide - a role hitherto suspected only for its reactions at metal atoms. The analogous reactions of trisulfides are readily reversible and not kinetically competent to inhibit hydrocarbon autoxidation, consistent with the longstanding view that organosulfur compounds must be oxidized to afford significant antioxidant activity. The reactivity of trisulfides and their oxides are contrasted with what is known of their shorter cousins and predictions are made and tested with regards to the reactivity of higher polysulfides and their 1-oxides - the insights from which may be exploited in the design of future antioxidants.

A multi-level quantum mechanics and molecular mechanics study of SN2 reaction at nitrogen: NH2Cl + OH(-) in aqueous solution

We employed a multi-level quantum mechanics and molecular mechanics approach to study the reaction NH2Cl + OH(-) in aqueous solution. The multi-level quantum method (including the DFT method with both the B3LYP and M06-2X exchange-correlation functionals and the CCSD(T) method, and both methods with the aug-cc-pVDZ basis set) was used to treat the quantum reaction region in different stages of the calculation in order to obtain an accurate potential of mean force. The obtained free energy activation barriers at the DFT/MM level of theory yielded a big difference of 21.8 kcal mol(-1) with the B3LYP functional and 27.4 kcal mol(-1) with the M06-2X functional respectively. Nonetheless, the barrier heights become very close when shifted from DFT to CCSD(T): 22.4 kcal mol(-1) and 22.9 kcal mol(-1) at CCSD(T)(B3LYP)/MM and CCSD(T)(M06-2X)/MM levels of theory, respectively. The free reaction energy obtained using CCSD(T)(M06-2X)/MM shows an excellent agreement with the one calculated using the available gas-phase data. Aqueous solution plays a significant role in shaping the reaction profile. In total, the water solution contributes 13.3 kcal mol(-1) and 14.6 kcal mol(-1) to the free energy barrier heights at CCSD(T)(B3LYP)/MM and CCSD(T)(M06-2X)/MM respectively. The title reaction at nitrogen is a faster reaction than the corresponding reaction at carbon, CH3Cl + OH(-).

Benchmark Databases for Nonbonded Interactions and Their Use To Test Density Functional Theory

We present four benchmark databases of binding energies for nonbonded complexes. Four types of nonbonded interactions are considered:  hydrogen bonding, charge transfer, dipole interactions, and weak interactions. We tested 44 DFT methods and 1 WFT method against the new databases; one of the DFT methods (PBE1KCIS) is new, and all of the other methods are from the literature. Among the tested methods, the PBE, PBE1PBE, B3P86, MPW1K, B97-1, and BHandHLYP functionals give the best performance for hydrogen bonding. MPWB1K, MP2, MPW1B95, MPW1K, and BHandHLYP give the best performances for charge-transfer interactions, and MPW3LYP, B97-1, PBE1KCIS, B98, and PBE1PBE give the best performance for dipole interactions. Finally, MP2, B97-1, MPWB1K, PBE1KCIS, and MPW1B95 give the best performance for weak interactions. Overall, MPWB1K is the best of all the tested DFT methods, with a relative error (highly averaged) of only 11%, and MPW1K, PBE1PBE, and B98 are the best of the tested DFT methods that do not contain kinetic energy density. Moving up the rungs of Jacob's ladder for nonempirical DFT, PBE improves significantly over the LSDA, and TPSS improve slightly over PBE (on average) for nonbonded interactions.

Evaluation of the SCF Combination of KS-DFT and 3D-RISM-KH; Solvation Effect on Conformational Equilibria, Tautomerization Energies, and Activation Barriers

The effect of solvation on conformational equilibria, tautomerization energies, and activation barriers in simple SN2 reactions is reproduced by using the self-consistent field coupling of the Kohn-Sham density functional theory (KS-DFT) for electronic structure and the three-dimensional reference interaction site model with the closure approximation of Kovalenko and Hirata (3D-RISM-KH) for molecular solvation structure. These examples are used in order to validate the implementation of the 3D-RISM-KH method in the Amsterdam Density Functional package. The computations of the free energy difference in the trans/gauche conformational equilibrium for 1,2-dichloroethane in different solvents; the relative tautomerization free energy for cytosine, isocytosine, and guanine; and the free energy activation barrier for a CH3X-type (X = F, Cl, Br) SN2 reaction exhibit agreement with the experimental data. The method is also applied to the electronic and hydration structure of carbon single-wall nanotubes of different diameters, including the effect of water located in the inner space of the nanotubes. A comparison with continuum models of solvation (including COSMO) as well as with other more precise and computationally more expensive calculations is made to demonstrate the accuracy and predictive capability of the new KS-DFT/3D-RISM-KH method.

Searching for Saddle Points by Using the Nudged Elastic Band Method: An Implementation for Gas-Phase Systems

A new implementation of the Nudged Elastic Band (NEB) optimization method is presented. This approach uses a global procedure that yields the whole reaction path, and thus it provides an alternative to the sequential optimization of the transition state and consequent calculation of the minimum energy path. Furthermore the algorithm is very useful when one is not sure if a saddle point exists, because it can be used to eliminate the possibility of a saddle point when one does not exist. Three different versions of the NEB algorithm have been implemented. The influences of various parameters and methodological choices on the performance of the method have been studied, and the quality of the results is assessed by comparison with the saddle point and minimum energy path calculations sequential method. Recommendations are made for algorithmic choices and default parameters.

Multireference Model Chemistries for Thermochemical Kinetics

By combining the generalized valence bond ansatz of correlated participating orbitals (CPO) with the complete-active-space prescription for selecting configurations and with the use of multireference second order perturbation theory (MRMP2) for including dynamical correlation, we define three levels of multireference (MR) theoretical model chemistries for electronic structure calculations of chemical reaction energies and barrier heights. The three levels differ in their choice of which orbitals are considered to be participating; the choices are called nominal (nom-CPO), moderate (mod-CPO), and extended (ext-CPO). Combining any of these three choices with a method for treatment of dynamical correlation energy and a one-electron basis set yields a theoretical model chemistry. Unlike the full-valence choice of active orbitals, the CPO choices lead to active spaces that contain the orbitals needed to include important static correlation effects on chemical reactions but do not increase with the size of the nonparticipating portion of the system, and hence they remain viable computational options even for many large and complex reacting systems. The accuracies of the new levels, combined with the MG3S basis set (a partially augmented, multiply polarized valence triple-ζ basis with appropriately tight d functions for 3p-block elements) and with the fully augmented correlation-consistent aug-cc-pVTZ basis set, are assessed against a previously presented database of barrier heights for diverse reaction types. We find that nom-CPO level captures the bulk of the static correlation energy, and MRMP2/nom-CPO calculations have an average error of only 1.4 kcal/mol in barrier heights, which may be compared to 5.0 kcal/mol for single-reference MP2 theory, 2.5 kcal/mol for CCSD, and 4.1 and 1.0 kcal/mol for the B3LYP and M06-2X density functionals, respectively. The accuracy of MRMP2/CPO for transition structure bond lengths and donor-acceptor distances is excellent, with a mean unsigned error of only 0.007 Å as compared to 0.018 Å for CCSD, 0.019 Å for M06-2X, and 0.039 Å for MP2 and B3LYP. We also introduce a new multireference diagnostic, called the M diagnostic, that allows one to measure the importance of static correlation in a given reagent or transition state.

Computation of Accurate Activation Barriers for Methyl-Transfer Reactions of Sulfonium and Ammonium Salts in Aqueous Solution

The energetics of methyl-transfer reactions from dimethylammonium, tetramethylammonium, and trimethylsulfonium to dimethylamine were computed with density functional theory, MP2, CBS-QB3, and quantum mechanics/molecular mechanics (QM/MM) Monte Carlo methods. At the CBS-QB3 level, the gas-phase activation enthalpies are computed to be 9.9, 15.3, and 7.9 kcal/mol, respectively. MP2/6-31+G(d,p) activation enthalpies are in best agreement with the CBS-QB3 results. The effects of aqueous solvation on these reactions were studied with polarizable continuum model, generalized Born/surface area (GB/SA), and QM/MM Monte Carlo simulations utilizing free-energy perturbation theory in which the PDDG/PM3 semiempirical Hamiltonian for the QM and explicit TIP4P water molecules in the MM region were used. In the aqueous phase, all of these reactions proceed more slowly when compared to the gas phase, since the charged reactants are stabilized more than the transition structure geometries with delocalized positive charges. In order to obtain the aqueous-phase activation free energies, the gas-phase activation free energies were corrected with the solvation free energies obtained from single-point conductor-like polarizable continuum model and GB/SA calculations for the stationary points along the reaction coordinate.

Assessment of density-functionals for describing the X(-) + CH3ONO2 gas-phase reactions with X = F, OH, CH2CN

The energetics of the ECO2, SN2@C and SN2@N channels of X(-) + CH3ONO2 (X = F, OH, CH2CN) gas-phase reactions were computed using the CCSD(T)/CBS method. This benchmark extends a previous study with X = OH [M. A. F. de Souza et al., J. Am. Chem. Soc., 2012, 134, 19004] and was used to ascertain the accuracy and robustness of nineteen density-functionals for describing these potential energy profiles (PEP) as well as the kinetic product distributions obtained from RRKM calculations. Assessments were based on the mean unsigned error (MUE), the mean signed error (MSE), the #best : #worst (BW) criterion and the statistical confidence interval (CI) for the MSE. In general, double-hybrid (DH) functionals perform better than the range-separated ones, and both are better than the global-hybrid functionals. Based on the MUE and CI criteria the B2GPPLYP, B2PLYP, M08-SO, BMK, ωB97X-D, CAM-B3LYP, M06, M08-HX, ωB97X and B97-K functionals show the best performance in the description of these PEPs. Within this set, the B2GPPLYP functional is the most accurate and robust. The RRKM results indicate that the DHs are the best for describing the selectivities of these reactions. Compared to CCSD(T), the B2PLYP method has a relative error of only ca. 1% for the selectivity and the accuracy to provide the correct conclusion concerning the nonstatistical behavior of these reactions.

Double-inversion mechanisms of the X⁻ + CH₃Y [X,Y = F, Cl, Br, I] SN2 reactions

The double-inversion and front-side attack transition states as well as the proton-abstraction channels of the X(-) + CH3Y [X,Y = F, Cl, Br, I] reactions are characterized by the explicitly correlated CCSD(T)-F12b/aug-cc-pVTZ(-PP) level of theory using small-core relativistic effective core potentials and the corresponding aug-cc-pVTZ-PP bases for Br and I. In the X = F case the double-inversion classical(adiabatic) barrier heights are 28.7(25.6), 15.8(13.4), 13.2(11.0), and 8.6(6.6) kcal mol(-1) for Y = F, Cl, Br, and I, respectively, whereas the barrier heights are in the 40-90 kcal mol(-1) range for the other 12 reactions. The abstraction channels are always above the double-inversion saddle points. For X = F, the front-side attack classical(adiabatic) barrier heights, 45.8(44.8), 31.0(30.3), 24.7(24.2), and 19.5(19.3) kcal mol(-1) for Y = F, Cl, Br, and I, respectively, are higher than the corresponding double-inversion ones, whereas for the other systems the front-side attack saddle points are in the 35-70 kcal mol(-1) range. The double-inversion transition states have XH···CH2Y(-) structures with Cs point-group symmetry, and the front-side attack saddle points have either Cs (X = F or X = Y) or C1 symmetry with XCY angles in the 78-88° range. On the basis of the previous reaction dynamics simulations and the minimum energy path computations along the inversion coordinate of selected XH···CH2Y(-) systems, we suggest that the double inversion may be a general mechanism for SN2 reactions.

Computational comparison of reactions of CS2 with CHX•- (X = F, Cl, and Br): Do F, Cl, and Br substitutions effect differently?

Calculations based on density functional theory (DFT) have been carried out for the gas-phase ion-molecule reactions of CS 2 with CHX •- ( X = F , Cl , and Br ) to investigate the effect of the halogen-substituent on the reaction mechanism. The doublet potential energy surfaces (PESs), involving two striking reaction patterns named by the middle-C attack and the end-S attack in terms of anion attack on the C and S atoms of CS 2, have been explored and characterized in detail. Compared with the results of the end-S attack, reaction with the middle-C attack pattern displays more efficiency on PES. For a given reaction pattern, the reactions of CHCl •- and CHBr •- occur with similar efficiencies and reactivity trends. The CHF •- anion displays remarkably different reactivity, which is traced to its lower electron binding energy and the effect of the electronegative fluorine substituent. This is in good agreement with the experimental observation. According to the property of the product branching ratios in experiment, dynamical effect is used to evaluate the reason that the end-S attack with less energetically favorable obtained in our theoretical studies is comparable to the middle-C attack.

Benchmarking Compound Methods (CBS-QB3, CBS-APNO, G3, G4, W1BD) against the Active Thermochemical Tables: A Litmus Test for Cost-Effective Molecular Formation Enthalpies

The theoretical atomization energies of some 45 CxHyOz molecules present in the Active Thermochemical Tables compilation and of particular interest to the combustion chemistry community have been computed using five composite model chemistries as titled. The species contain between 1-8 "heavy" atoms, and a few are conformationally diverse with up to nine conformers. The enthalpies of formation at 0 and 298.15 K are then derived via the atomization method and compared against the recommended values. In general, there is very good agreement between our averaged computed values and those in the ATcT; those for 1,3-cyclopentadiene exceptionally differ considerably, and we show from isodesmic reactions that the true value for 1,3-cyclopentadiene is closer to 134 kJ mol(-1) than the reported 101 kJ mol(-1). If one is restricted to using a single method, statistical measures indicate that the best methods are in the rank order G3 ≈ G4 > W1BD > CBS-APNO > CBS-QB3. The CBS-x methods do on average predict ΔfH(⊖)(298.15 K) within ≈5 kJ mol(-1) but are prone to occasional lapses. There are statistical advantages to be gained from using a number of methods in tandem, and all possible combinations have been tested. We find that the average formation enthalpy coming from using CBS-APNO/G4, CBS-APNO/G3, and G3/G4 show lower mean signed and mean unsigned errors, and lower standard and root-mean-squared deviations, than any of these methods in isolation. Combining these methods also leads to the added benefit of providing an uncertainty rooted in the chemical species under investigation. In general, CBS-APNO and W1BD tend to underestimate the formation enthalpies of target species, whereas CBS-QB3, G3, and G4 have a tendency to overestimate the same. Thus, combining CBS-APNO with a G3/G4 combination leads to an improvement in all statistical measures of accuracy and precision, predicting the ATcT values to within 0.14 ± 4.21 kJ mol(-1), thus rivalling "chemical accuracy" (±4.184 kJ mol(-1)) without the excessive cost associated with higher-level methods such as W1BD.

Quantum mechanical study of solvent effects in a prototype SN2 reaction in solution: Cl- attack on CH3Cl

The nucleophilic attack of a chloride ion on methyl chloride is an important prototype SN2 reaction in organic chemistry that is known to be sensitive to the effects of the surrounding solvent. Herein, we develop a highly accurate Specific Reaction Parameter (SRP) model based on the Austin Model 1 Hamiltonian for chlorine to study the effects of solvation into an aqueous environment on the reaction mechanism. To accomplish this task, we apply high-level quantum mechanical calculations to study the reaction in the gas phase and combined quantum mechanical/molecular mechanical simulations with TIP3P and TIP4P-ew water models and the resulting free energy profiles are compared with those determined from simulations using other fast semi-empirical quantum models. Both gas phase and solution results with the SRP model agree very well with experiment and provide insight into the specific role of solvent on the reaction coordinate. Overall, the newly parameterized SRP Hamiltonian is able to reproduce both the gas phase and solution phase barriers, suggesting it is an accurate and robust model for simulations in the aqueous phase at greatly reduced computational cost relative to comparably accurate ab initio and density functional models.

A concise review of computational studies of the carbon dioxide–epoxide copolymerization reactions

The production of polycarbonates from carbon dioxide and epoxides is an important route by which waste CO₂ is converted into useful products. This review surveys the use of computational chemistry toward understanding this reaction.

Benchmark study on the smallest bimolecular nucleophilic substitution reaction: H⁻+CH₄-->CH₄+H⁻

We report here a benchmark study on the bimolecular nucleophilic substitution (S(N)2) reaction between hydride and methane, for which we have obtained reference energies at the coupled cluster toward full configuration-interaction limit (CC-cf/CBS). Several wavefunction (HF, MP2, coupled cluster) and density functional methods are compared for their reliability regarding these reference data.

Inverse neighboring group participation: explanation of an unusual S→N alkyl migration of isothiuronium salts containing a lactone group

A detailed experimental and DFT study of the S to N alkyl migration of substituted S-(1(3H)-isobenzofuranon-3-yl)isothiuronium bromide to N,N'-dimethyl-N-(3-oxo-1,3-dihydro-2-benzofuran-1-yl)thiourea provided evidence for the existence of an unusual double displacement mechanism involving two consecutive back-side S(N)2 reactions where a carboxylate anion has a crucial role both as a leaving group as well as an internal nucleophile. The thiazetidine zwitterionic species that is involved in this mechanism as an intermediate was characterized by infrared multiphoton dissociation spectroscopy and was trapped with methyl iodide. It was found that the intermediate has a structure of a free ion pair. The double-displacement mechanism can be considered as a new type of inverse lactone neighboring group participation.

DFT study on the reactions of ClO⁻/BrO⁻ with RCl (R = CH₃, C₂H₅, and C₃H₇) in gas phase

Gas-phase reactions of ClO(-)/BrO(-) with RCl (R = CH₃, C₂H₅, and C₃H₇) have been investigated in detail using the popular DFT functional BHandHLYP/aug-cc-pVDZ level of theory. As a result, our findings strongly suggest that the type of reaction is firstly initiated by a typical SN2 fashion. Subsequently, two competitive substitution steps, named as SN2-induced substitution and SN2-induced elimination, respectively, would proceed before the initial SN2 product ion-dipole complex separates, in which the former exhibits less reactivity than the latter. Those are consistent with relevant experimental results. Moreover, we have also explored reactivity difference for the title reactions in term of some factors derived from methyl group, p-π electronic conjugation, ionization energy (IE), as well as molecular orbital (MO) analysis.