High-level analytical potential-energy-surface-based dynamics of the OH- + CH3CH2Cl SN2 and E2 reactions in full (24) dimensions.

We develop a coupled-cluster full-dimensional global potential energy surface (PES) for the OH- + CH3CH2Cl reactive system, using the Robosurfer program package, which automatically samples configurations along PES-based trajectories as well as performs ab initio computations with Molpro and fitting with the monomial symmetrization approach. The analytical PES accurately describes both the bimolecular nucleophilic substitution (SN2) and elimination (E2) channels leading to the Cl- + CH3CH2OH and Cl- + H2O + C2H4 products, respectively, and allows efficient quasi-classical trajectory (QCT) simulations. QCT computations on the new PES provide accurate statistically-converged integral and differential cross sections for the OH- + CH3CH2Cl reaction, revealing the competing dynamics and mechanisms of the SN2 and E2 (anti, syn, ß-α transfer) channels as well as various additional pathways leading to induced inversion of the CH3CH2Cl reactant, H-exchange between the reactants, H2O⋯Cl- complex formation, and H2O + CH3CHCl- products via proton abstraction.

The influence of fluorination on the dynamics of the F- + CH3CH2I reaction.

The competition between the bimolecular nucleophilic substitution (SN2) and base-induced elimination (E2) reaction and their intrinsic reactivity is of key interest in organic chemistry. To investigate the effect of suppressing the E2 pathway on SN2 reactivity, we compared the reactions F- + CH3CH2I and F- + CF3CH2I. Differential cross-sections have been measured in a crossed-beam setup combined with velocity map imaging, giving insight into the underlying mechanisms of the individual pathways. Additionally, we employed a selected-ion flow tube to obtain reaction rates and high-level ab initio computations to characterize the different reaction pathways and product channels. The fluorination of the ß-carbon not only suppresses the E2-reaction but opens up additional channels involving the abstraction of fluorine. The overall SN2 reactivity is reduced compared to the non-fluorinated iodoethane. This reduction is presumably due to the competition with the highly reactive channels forming FHF- and CF2CI-.

Dynamics of the Cl- + CH3I reaction on a high-level ab initio analytical potential energy surface.

We have developed a full-dimensional analytical ab initio potential energy surface (PES) for the Cl- + CH3I reaction using the Robosurfer program system. The energy points have been computed using a robust composite method defined as CCSD-F12b + BCCD(T) - BCCD with the aug-cc-pVTZ(-PP) basis set and have been fitted by the permutationally invariant polynomial approach. Quasi-classical trajectory simulations on the new PES reveal that two product channels are open in the collision energy (Ecoll) range of 1-80 kcal/mol, i.e., SN2 leading to I- + CH3Cl and iodine abstraction (above ∼45 kcal/mol) resulting in ICl- + CH3. Scattering angle, initial attack angle, product translational energy, and product internal energy distributions show that the SN2 reaction is indirect at low Ecoll and becomes direct-rebound-back-side (CH3-side) attack-type, as Ecoll increases. Iodine abstraction mainly proceeds with direct stripping mechanism with side-on/back-side attack preference. Comparison with crossed-beam experiments and previous direct dynamics simulations shows quantitative or qualitative agreement and also highlights possible theoretical and/or experimental issues motivating further research.

Protonation of serine: conformers, proton affinities and gas-phase basicities at the "gold standard" and beyond.

The potential energy surfaces (PESs) of serine and its protonated counterparts are investigated to determine the structures of the minima. A total of 95 neutral serine, 15 N-(amino-) and 46 O-(carbonyl-)protonated serine conformers are found. Their relative energies, geometries and harmonic vibrational frequencies are determined at the MP2/aug-cc-pVDZ level of theory. To obtain highly accurate thermodynamic values, further computations are performed: the ten conformers with the lowest relative energies from each molecule type (neutral, N- and O-protonated) are further optimized using the explicitly correlated CCSD(T)-F12a/cc-pVDZ-F12 method (for neutral serine, harmonic vibrational frequencies were also computed). In addition, auxiliary corrections were determined: basis-set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation and second-order Douglas-Kroll relativistic effects along with zero-point energy contributions. Two important thermodynamic parameters (at 298.15 K), proton affinity (PA)/gas-phase basicity (GB) are calculated considering the two different protonation sites: 218.05 ± 0.2/209.86 ± 0.6 kcal mol-1 and 205.87 ± 0.2/196.36 ± 0.3 kcal mol-1 for the amino and carbonyl sites, respectively. The uncertainty of the determined values was approximated based on various sources including auxiliary corrections, basis-set effects, harmonic vibrational frequencies.

Benchmark Ab Initio Determination of the Conformers, Proton Affinities, and Gas-Phase Basicities of Cysteine.

A systematic conformational mapping combined with literature data leads to 85 stable neutral cysteine conformers. The implementation of the same mapping process for the protonated counterparts reveals 21 N-(amino-), 64 O-(carbonyl-), and 37 S-(thiol-)protonated cysteine conformers. Their relative energies and harmonic vibrational frequencies are given at the MP2/aug-cc-pVDZ level of theory. Further benchmark ab initio computations are performed for the 10 lowest-lying neutral and protonated amino acid conformers (for each type) such as CCSD(T)-F12a/cc-pVDZ-F12 geometry optimizations (and frequency computations for cysteine) as well as auxiliary correction computations of the basis set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation effects, second-order Douglass-Kroll relativistic effects, and zero-point energy contributions. Boltzmann-averaged 0 (298.15) K proton affinity and [298.15 K gas-phase basicity] values of cysteine are predicted to be 214.96 (216.39) [208.21], 201.83 (203.55) [194.16], and 193.31 (194.74) [186.40] kcal/mol for N-, O-, and S-protonation, respectively, also considering the previously described auxiliary corrections.

Benchmark ab initio proton affinity of glycine.

A systematic conformational search reveals three N- (amino) and eight O- (carbonyl) protonated glycine conformers with benchmark equilibrium(adiabatic) relative energies in the 0.00-7.51(0.00-7.37) and 25.91-31.61(24.45-30.28) kcal mol-1 ranges, respectively. Benchmark ab initio structures of the glycine conformers and its protonated species are obtained at the CCSD(T)-F12b/aug-cc-pVTZ level of theory and the relative energy computations consider basis-set effects up to aug-cc-pVQZ with CCSD(T)-F12b, electron correlation up to CCSDT(Q), core correlation corrections, scalar relativistic effects, and zero-point energy contributions. The best predictions for Boltzmann-averaged 0(298.15) K proton affinities and [298.15 K gas-phase basicities] of glycine are 211.00(212.43)[204.75] and 186.38(187.64)[180.21] kcal mol-1 for N- and O-protonation, respectively, in excellent agreement with experiments.

Conformers of dehydrogenated glycine isomers.

We report a comprehensive ab initio investigation of the conformers of dehydrogenated glycine radicals using the STO-3G, 3-21G, and aug-cc-pVDZ (aVDZ) basis sets and the UHF and UMP2 (H2 N-CH-COOH and HN-CH2 -COOH) as well as MCSCF and MRCI (H2 N-CH2 -COO) methods via two different conformational search strategies generating initial structures for optimizations by (a) removing H atoms from glycine conformers and (b) scanning torsional angles describing internal rotation along the CC, CN, and CO (except for H2 N-CH2 -COO) bonds of the radicals. We find four H2 N-CH-COOH {InCH , IInCH , IIInCH , IVnCH } and seven HN-CH2 -COOH {IpNH , IIpNH , IIInNH , IVpNH VnNH , VIpNH , VIIpNH } conformers with classical(adiabatic) relative energies of {0.00(0.00), 1.57(1.55), 5.25(5.03), 9.85(9.72)} and {0.00(0.00), 0.78(1.06), 1.93(2.08), 3.34(3.16), 3.39(3.29), 5.00(4.86), 9.27(8.87)} kcal/mol, respectively, obtained with UCCSD(T)-F12b/aug-cc-pVTZ(+UCCSD(T)-F12b/aVDZ ZPE correction) and four H2 N-CH2 -COO {IpCOO , IInCOO , IIIpCOO , IVnCOO } conformers with MRCI-F12+Q/aVDZ(+MRCI/aVDZ ZPE correction) energies of {0.00(0.00), 1.65(1.64), 1.78(1.75), 2.21(2.21)} kcal/mol, where n and p denote C1 and Cs symmetry. The MRCI-F12+Q[UCCSD(T)-F12b] InCH → IpNH and InCH → IpCOO classical(adiabatic) isomerization energies are 18.51(17.32)[21.20(20.01)] and 31.88(31.66) kcal/mol, respectively.