Dynamics study of the OH + NH3 hydrogen abstraction reaction using QCT calculations based on an analytical potential energy surface

A Hessian-Free Method to Prevent Zero-Point Energy Leakage in Classical Trajectories

The problem associated with the zero-point energy (ZPE) leak in classical trajectory calculations is well known. Since ZPE is a manifestation of the quantum uncertainty principle, there are no restrictions on energy during the classical propagation of nuclei. This phenomenon can lead to unphysical results, such as forming products without the ZPE in the internal vibrational degrees of freedom (DOFs). The ZPE leakage also permits reactions below the quantum threshold for the reaction. We have developed a new Hessian-free method, inspired by the Lowe-Andersen thermostat model, to prevent energy dipping below a threshold in the local-pair (LP) vibrational DOFs. The idea is to pump the leaked energy to the corresponding local vibrational mode taken from the other vibrational DOFs. We have applied the new correction protocol on the ab-initio ground-state molecular dynamics simulation of the water dimer (H₂O)₂, which dissociates due to unphysical ZPE spilling from high-frequency OH modes. The LP-ZPE method has been able to prevent the ZPE spilling of the OH stretching modes by pumping back the leaked energy into the corresponding modes, while this energy is taken from the other modes of the dimer itself, keeping the system as a microcanonical ensemble.

Early Chemistry of Nicotine Degradation in Heat-Not-Burn Smoking Devices and Conventional Cigarettes: Implications for Users and Second- and Third-Hand Smokers

Nicotine exposure results in health risks not only for smokers but also for second- and third-hand smokers. Unraveling nicotine's degradation mechanism and the harmful chemicals that are produced under different conditions is vital to assess exposure risks. We performed a theoretical study to describe the early chemistry of nicotine degradation by investigating two important reactions that nicotine can undergo: hydrogen abstraction by hydroxyl radicals and unimolecular dissociation. The former contributes to the control of the degradation mechanism below 800 K due to a non-Arrhenius kinetics, which implies an enhancement of reactivity as temperature decreases. The latter becomes important at higher temperatures due to its larger activation energy. This change in the degradation mechanism is expected to affect the composition of vapors inhaled by smokers and room occupants. Conventional cigarettes, which operate at temperatures higher than 1000 K, are more prone to yield harmful pyridinyl radicals via nicotine dissociation, while nicotine in electronic cigarettes and vaporizers, with operating temperatures below 600 K, will be more likely degraded by hydroxyl radicals, resulting in a vapor with a different composition. Although low-temperature nicotine delivery devices have been claimed to be less harmful due to their nonburning operating conditions, the non-Arrhenius kinetics that we observed for the degradation mechanism below 873 K suggests that nicotine degradation may be more rapidly initiated as temperature is reduced, indicating that these devices may be more harmful than it is commonly assumed.

Quantum calculations of the photoelectron spectra of the OH-·NH3 anion: implications for OH + NH3→ H2O + NH2 reaction dynamics

We present the results of quantum dynamics calculations for analyzing the experimentally measured photoelectron spectra of the OH^-·NH₃ anion complex. Detachment of an excess electron of OH^-·NH₃ initially produces a molecular arrangement, which is close to the transition-state structure of the neutral OH + NH₃→ H₂O + NH₂ hydrogen abstraction reaction due to the Franck-Condon principle, and thus finally leads to the OH + NH₃ or H₂O + NH₂ asymptotic channel. We used both the path integral method and the reduced-dimensionality quantum wave packet method to simulate the photoelectron spectra of the OH^-·NH₃ anion. The calculated spectra were found to be in qualitative agreement with the experimental spectra. It was found that the photodetached complex mainly dissociates into the OH + NH₃ channel; however, we found that the hydrogen exchange process also contributes to the photodetachment spectra.

Quasi-Classical Trajectory Study of the CN + NH3 Reaction Based on a Global Potential Energy Surface

Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.

High-level theoretical study of the reaction between hydroxyl and ammonia: Accurate rate constants from 200 to 2500 K

Hydrogen abstraction from NH₃ by OH to produce H₂O and NH₂-an important reaction in combustion of NH₃ fuel-was studied with a theoretical approach that combines high level quantum chemistry and advanced chemical kinetics methods. Thermal rate constants calculated from first principles agree well (within 5%-20%) with available experimental data over a temperature range that extends from 200 to 2500 K. Quantum mechanical tunneling effects were found to be important; they lead to a decided curvature and non-Arrhenius behavior for the rate constant.

Mechanisms of the Formation of Adenine, Guanine, and Their Analogues in UV-Irradiated Mixed NH3:H2O Molecular Ices Containing Purine

We investigated the formation mechanisms of the nucleobases adenine and guanine and the nucleobase analogues hypoxanthine, xanthine, isoguanine, and 2,6-diaminopurine in a UV-irradiated mixed 10:1 H₂O:NH₃ ice seeded with precursor purine by using ab initio and density functional theory computations. Our quantum chemical investigations suggest that a multistep reaction mechanism involving purine cation, hydroxyl and amino radicals, together with water and ammonia, explains the experimentally obtained products in an independent study. The relative abundances of these products appear to largely follow from relative thermodynamic stabilities. The key role of the purine cation is likely to be the reason why purine is not functionalized in pure ammonia ice, where cations are promptly neutralized by free electrons from NH₃ ionization. Amine group addition to purine is slightly favored over hydroxyl group attachment based on energetics, but hydroxyl is much more abundant due to higher abundance of H₂O. The amino group is preferentially attached to the 6 position, giving 6-aminopurine, that is, adenine, while the hydroxyl group is preferentially attached to the 2 position, leading to 2-hydroxypurine. A second substitution by hydroxyl or amino group occurs at either the 6 or the 2 position depending on the first substitution. Given that H₂O is far more abundant than NH₃ in the experimentally studied ices (as well as based on interstellar abundances), xanthine and isoguanine are expected to be the most abundant bi-substituted photoproducts. Key Words: Astrophysical ice-Abiotic organic synthesis-Nucleic acids-Origin of life-RNA world. Astrobiology 17, 771-785.

QCT dynamics study of OH/OD + GeH4 reactions. The problem of water bending excitation

The difficulties in the description of the water bending population could be related to the zero-point violation problem in QCT calculations.

Branching Ratios and Vibrational Distributions in Water-Forming Reactions of OH and OD Radicals with Methylamines

Reactions of OH and OD radicals with (CH3)3N, (CH3)2NH, and CH3NH2 were studied by Fourier transform infrared emission spectroscopy (FTIR) of the water product molecules from a fast-flow reactor at 298 K. The rate constants (4.4 ± 0.5) × 10(-11), (5.2 ± 0.8) × 10(-11), and (2.0 ± 0.4) × 10(-11) cm(3) molecule(-1) s(-1) were determined for OD + (CH3)3N, (CH3)2NH, and CH3NH2, respectively, by comparing the HOD emission intensities to the HOD intensity from the OD reaction with H2S. Abstraction from the nitrogen site competes with abstraction from the methyl group, as obtained from an analysis of the HOD and D2O emission intensities from the OD reactions with the deuterated reactants, (CD3)2NH and CD3NH2. After adjustment for the hydrogen-deuterium kinetic isotope effect, the product branching fractions of the hydrogen abstraction from the nitrogen for di- and monomethylamine were found to be 0.34 ± 0.04 and 0.26 ± 0.05, respectively. Vibrational distributions of the H2O, HOD, and D2O molecules are typical for direct hydrogen atom abstraction from polar molecules, even though activation energies are negative because of the formation of pre-transition-state complexes. Comparison is made to the reactions of hydroxyl radicals with ammonia and with other compounds with primary C-H bonds to discuss specific features of disposal of energy to water product.