Potential-energy surface of the CH3CO + O2 reaction has been calculated by ab initio quantum chemistry methods. The geometries were optimized using the second-order Moller-Plesset theory (MP2) with the 6-311G(d,p) basis set and the coupled-cluster theory with single and double excitations (CCSD) with the correlation consistent polarized valence double zeta (cc-pVDZ) basis set. The relative energies were calculated using the Gaussian-3 second-order Moller-Plesset theory with the CCSD/cc-pVDZ geometries. Multireference self-consistent-field and MP2 methods were also employed using the 6-311G(d,p) and 6-311++G(3df,2p) basis sets. Both addition/elimination and direct abstraction mechanisms have been investigated. It was revealed that acetylperoxy radical [CH3C(O)OO] is the initial adduct and the formation of OH and alpha-lactone [CH2CO2(1A')] is the only energetically accessible decomposition channel. The other channels, e.g., abstraction, HO2 + CH2CO, O + CH3CO2, CO + CH3O2, and CO2 + CH3O, are negligible. Multichannel Rice-Ramsperger-Kassel-Marcus theory and transition state theory (E-resolved) were employed to calculate the overall and individual rate coefficients and the temperature and pressure dependences. Fairly good agreement between theory and experiments has been obtained without any adjustable parameters. It was concluded that at pressures below 3 Torr, OH and CH2CO2(1A') are the major nascent products of the oxidation of acetyl radicals, although CH2CO2(1A') might either undergo unimolecular decomposition to form the final products of CH2O + CO or react with OH and Cl to generate H2O and HCl. The acetylperoxy radicals formed by collisional stabilization are the major products at the elevated pressures. In atmosphere, the yield of acetylperoxy is nearly unity and the contribution of OH is only marginal.

In this paper, the surface tension and critical properties for the TIP4P/2005 and BLYPSP-4F models are reported. A clear dependence of surface tension on the van der Waals cutoff radius (r vdw ) is shown when van der Waals interactions are modeled with a simple cutoff scheme. A linear extrapolation formula is proposed that can be used to determine the infinite r vdw surface tension through a few simulations with finite r vdw . A procedure for determining liquid and vapor densities is proposed that does not require fitting to a profile function. Although the critical temperature of water is also found to depend on the choice of r vdw , the dependence is weaker. We argue that a r vdw of 1.75 nm is a good compromise for water simulations when long-range van der Waals correction is not applied. Since the majority of computational programs do not support rigorous treatment of long-range dispersion, the establishment of a minimal acceptable r vdw is important for the simulation of a variety of inhomogeneous systems, such as water bubbles, and water in confined environments. The BLYPSP-4F model predicts room temperature surface tension marginally better than TIP4P/2005 but overestimates the critical temperature. This is expected since only liquid configurations were fit during the development of the BLYPSP-4F potential. The potential is expected to underestimate the stability of vapor and thus overestimate the region of stability for the liquid. C 2015 AIP Publishing LLC. [http://dx

Relative rate studies were carried out to determine the temperature dependent rate constant ratio k1/k2a: CH3O· + O2 → HCHO + HO2· and CH3O· + NO2 (+M) → CH3ONO2 (+M) over the temperature range 250–333 K in an environmental chamber at 700 Torr using Fourier transform infrared detection. Absolute rate constants k2 were determined using laser flash photolysis/laser-induced fluorescence under the same conditions. The analogous experiments were carried out for the reactions of the perdeuterated methoxy radical (CD3O·). Absolute rate constants k2 were in excellent agreement with the recommendations of the JPL Data Evaluation panel. The combined data (i.e., k1/k2 and k2) allow the determination of k1 as 1.3(–0.5)(+0.9) × 10(–14) exp[−(663 ± 144)/T] cm(3) s(–1), corresponding to 1.4 × 10(–15) cm(3) s(–1) at 298 K. The rate constant at 298 K is in excellent agreement with previous work, but the observed temperature dependence is less than was previously reported. The deuterium isotope effect, kH/kD, can be expressed in the Arrhenius form as k1/k3 = (1.7(–0.4)(+0.5)) exp((306 ± 70)/T). The deuterium isotope effect does not appear to be greatly influenced by tunneling, which is consistent with a previous theoretical work by Hu and Dibble. (Hu, H.; Dibble, T. S., J. Phys. Chem. A 2013, 117, 14230–14242.)

Water is perhaps the most simulated liquid. Recently three water models have been developed following the adaptive force matching (AFM) method that provides excellent predictions of water properties with only electronic structure information as a reference.Compared to many other electronic structure based force fields that rely on fairly sophisticated energy expressions, the AFM water models use point-charge based energy expressions that are