Heterogeneous Interaction of N2O5 with HCl Doped H2SO4 under Stratospheric Conditions: ClNO2 and Cl2 Yields

Distinctive Heterogeneous Reaction Mechanism of ClNO2 on the Air-Water Surface Containing Cl

The heterogeneous reaction of nitryl chloride (ClNO2) on the air-water surface plays a significant role in the chloride lifecycle. The air-water surface is ubiquitous on ice surfaces under supercooled conditions, affecting the uptake and heterogeneous reaction processes of trace gases. Previous studies suggest that ClNO2 is formed on Cl-doped ice surfaces following the N2O5 uptake. Herein, a distinctive heterogeneous reaction mechanism of ClNO2 is suggested on an air-water surface containing Cl under supercooled conditions using combined classic molecular dynamics (MD) and Born-Oppenheimer MD simulations. It is found that N2O5 dissociates into a NO2+ and NO3- ionic pair on the top air-water surface. In the top layer of the surface containing barely any Cl-, NO2+ proceeds through hydrolysis and produces H3O+ and HNO3. Thus, surface acidification appears because of H3O+ yields. With NO2+ diffusion to the deep layer of the surface, NO2+ reacts with Cl- and forms ClNO2. Note that ClNO2 formation competes with NO2+ hydrolysis, and the rate of ClNO2 formation is 27.7[Cl-] larger than that of NO2+ hydrolysis. Afterward, the reaction of ClNO2 with Cl- becomes barrierless with the catalysis by H3O+, which is not feasible on a neutral surface. Cl2 is thus generated and escapes into the atmosphere (low solubility of Cl2), contributing to the Cl radical. The proposed mechanism bolsters the current understanding of ClNO2's fate and its role in Cl chemistry in extremely cold environments like the Arctic and other high-latitude regions in wintertime.

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.