The kinetics of the reaction of atomic sulfur with acetylene (S (3P) + C2H2) were investigated experimentally via the flash photolysis resonance fluorescence method, and the theoretical potential energy surface for the reaction CN + SO was modeled via the density functional and configuration interaction computational methods. Sulfur is of interest in modern chemistry due to its relevance in combustion and atmospheric chemistry, in the Claus process, in soot and diamond-film formation and in astrochemistry. Experimental conditions ranged from 295 – 1015 K and 10 – 400 Torr of argon. Pressure-dependence was shown at all experimental temperatures. The room temperature high-pressure limit second order rate constant was (2.10 ± 0.08) × 10-13 cm3 molecule-1 s-1. The Arrhenius plot of the high-pressure limit rate constants gave an Ea of (11.34 ± 0.03) kJ mol-1 and a pre-exponential factor of (2.14 ± 0.19) × 10-11 cm3 molecule-1 s-1. S (3P) + C2H2 is likely an adduct forming reaction due to pressure-dependence (also supported by a statistical mechanics analysis) which involves intersystem crossing. The potential energy surface for CN + SO was calculated at the B3LYP/6-311G(d) level and refined at the QCISD/6-311G(d) level. The PES was compared to that of the analogous reaction …

In this work we address two problems in computational chemistry relevant to biomolecular modeling. In the first project, we consider the conformer space of melatonin as a a representative example of “real-life” flexible biomolecules. Geometries for all 52 unique conformers are optimized using spin-component scaled MP2, and then relative energies are obtained at the CCSD (T) level near the complete basis set limit. These are then used to validate a variety of DFT methods with and without empirical dispersion corrections, as well as some lower-level ab initio methods. Basis set convergence is found to be relatively slow due to internal C-H…O and C-H…N contacts. Absent dispersion corrections, many DFT functionals will transpose the two lowest conformers. Dispersion corrections resolve the problem for most functionals. Double hybrids yield particularly good performance, as does MP2.5. In the second project, we propose a simple DFT-based diagnostic for nondynamical correlation effects. Aλ= (1-TAE [ΧλC]/TAE[XC])/λ where TAE is the total atomization energy, XC the “pure” DFT exchange-correlation functional, and ΧλC the corresponding hybrid with 100λ% HF-type exchange. The diagnostic is a good predictor for sensitivity of energetics to the level of theory, unlike most of the wavefunction-based diagnostics. For GGA functionals, Aλ values approaching unity …