We characterize the stationary points along the Walden inversion, front-side attack, and double-inversion pathways of the X– + CH3Y and X– + SiH3Y [X, Y = F, Cl, Br, I] SN2 reactions using chemically accurate CCSD­(T)-F12b/aug-cc-pVnZ [n = D, T, Q] levels of theory. At the carbon center, Walden inversion dominates and proceeds via prereaction (X–···H3CY) and postreaction (XCH3···Y–) ion-dipole wells separated by a usually submerged transition state (X–H3C–Y)−, front-side attack occurs over high barriers, double inversion is the lowest-energy retention pathway for X = F, and hydrogen- (F–···HCH2Y) and halogen-bonded (X–···YCH3) complexes exist in the entrance channel. At the silicon center, Walden inversion proceeds through a single minimum (X–SiH3–Y)−, the front-side attack is competitive via a usually submerged transition state separating pre- and postreaction minima having X–Si–Y angles close to 90°, double inversion occurs over positive, often high barriers, and hydrogen- and halogen-bonded complexes are not found. In addition to the SN2 channels (Y– + CH3X/SiH3X), we report reaction enthalpies for proton abstraction (HX + CH2Y–/SiH2Y–), hydride substitution (H– + CH2XY/SiH2XY), XH···Y– complex formation (XH···Y– + 1CH2/1SiH2), and halogen abstraction (XY + CH3 –/SiH3 – and XY– + CH3/SiH3).

We develop a full-dimensional analytical potential energy surface (PES) for the F− + SiH3Cl reaction using Robosurfer for automatically sampling the configuration space, the robust [CCSD-F12b + BCCD(T) − BCCD]/aug-cc-pVTZ composite level of theory for computing the energy points, and the permutationally invariant polynomial method for fitting. Evolution of the fitting error and the percentage of the unphysical trajectories are monitored as a function of the iteration steps/number of energy points and polynomial order. Quasi-classical trajectory simulations on the new PES reveal rich dynamics resulting in high-probability SN2 (SiH3F + Cl−) and proton-transfer (SiH2Cl− + HF) products as well as several lower-probability channels, such as SiH2F− + HCl, SiH2FCl + H−, SiH2 + FHCl−, SiHFCl− + H2, SiHF + H2 + Cl−, and SiH2 + HF + Cl−. The Walden-inversion and front-side-attack-retention SN2 pathways are found to be competitive, producing nearly racemic products at high collision energies. The detailed atomic-level mechanisms of the various reaction pathways and channels as well as the accuracy of the analytical PES are analyzed along representative trajectories.

We determine the proton affinity (PA) and gas-phase basicity (GB) of amino acid α-alanine at a chemically accurate level by performing explicitly-correlated CCSD(T)-F12b/aug-cc-pVDZ geometry optimizations and normal mode vibrational frequency calculations as well as CCSD(T)-F12b/aug-cc-pVTZ energy computations at the possible neutral and protonated geometries. Temperature effects at 298.15 K considering translational, rotational, and vibrational enthalpy and entropy corrections are obtained via standard statistical mechanics utilizing the molecular geometries and the harmonic vibrational energy levels. Both the amino nitrogen (N) and the carbonyl oxygen (O) atoms are proven to be potential protonation sites and a systematic conformational search reveals 3 N-and 9 O-protonated conformers in the 0.00-7.88 and 25.43-30.43 kcal/mol energy ranges at 0 K, respectively. The final computed PA and GB values at (0)298.15 K in case of N-protonation are (214.47)216.80 and 207.07 kcal/mol, respectively, whereas the corresponding values for O-protonation are (189.04)190.63 and 182.31 kcal/mol. The results of the benchmark high-level coupled-cluster computations are utilized to assess the accuracy of several lowerlevel cost-effective methods such as MP2 and density functional theory with various functionals (SOGGA11-X, M06-2X, PBE0, B3LYP, M06, TPSS).

No abstract