À) with water is modelled in ONONO 2 Á (H 2 O) 4 clusters. Molecular Dynamics simulations using second-order Møller-Plesset perturbation (MP2) theory support the feasibility of the reaction of a charge-separated species to produce HONO and nitric acid.It is well known that nitrous acid (HONO) plays a major role in the chemistry of the atmosphere, e.g. as a source of hydroxyl (OH) radicals. 1,2 While heterogeneous reactions of NO 2 with water on tropospheric surfaces are believed to be important in generating HONO, the actual mechanisms involved are not known with certainty. On a chemical basis, the involvement of the NO 2 dimer seems reasonable, 3 although other mechanisms have been proposed for production of HONO both in the dark and under irradiation. [4][5][6][7][8][9][10][11] Initial steps of one proposed mechanism 3 include NO 2 dimerization to form N 2 O 4 and autoionization of asymmetric N 2 O 4 , i.e. ONONO 2 , to form (NO + )(NO 3 À ).ONONO 2 -(NO + )(NO 3 À )Reaction of the (NO + )(NO 3 À ) ion pair with water to form HONO and nitric acid (HONO 2 ) is the final step. A number of studies have addressed the formation of NO 2 dimers in both the symmetric (with a central NN bond) and asymmetric (ONONO 2 ) forms, and the interconversion between the two. [12][13][14][15][16][17][18] While the symmetric form of the dimer is favoured, there are pathways that can generate the asymmetric trans-ONONO 2 form that is the most likely precursor to the ion pair and subsequently HONO. 12,13,15,16,18 de Jesus Medeiros and Pimentel, in a study of the dimerization of NO 2 and the symmetric to asymmetric isomerization of the dimer in water clusters, predicted direct HONO formation from the NO 2 Á(H 2 O) n + NO 2 reaction and from the transition state for symmetric to asymmetric isomerization. 15 These modifications of the proposed mechanism described above bypass formation of (NO + )(NO 3 À ) through autoionization or, in the case of the NO 2 Á(H 2 O) n + NO 2 reaction, formation of the NO 2 dimer. Further exploration of the (NO + )(NO 3 À ) path will allow for comparison between potential routes and provide additional insight into the detailed mechanism that converts NO 2 in the presence of water to HONO. This is particularly useful for understanding other important reactions of NO 2 on surfaces in the presence of water, such as the reaction with HCl to form ClNO, since such chemistry likely shares common intermediates with NO 2 hydrolysis. [19][20][21] In ab initio Molecular Dynamics [AIMD], the reaction trajectory is propagated classically on the ab initio potential surface providing a means for studying reactive processes and the influence of the immediate environment. 22 The system wave function or electron density is calculated at each time step, i.e. on-the-fly, for determination of the energy, the interatomic forces and other properties. Partial charges are an obvious metric to investigate charge separation in the reaction of an ion pair as charge characteristics can indicate a shift in the nature of the interactio...