Molecular dynamics simulation of liquid N2O4<-->2NO2 by orientation-sensitive pairwise potential. III. Reaction dynamics

The dissociation and association dynamics of N2O4 [see text] 2NO2 in liquid state are studied by classical molecular dynamics simulations of reactive liquid NO2. An OSPP+LJ potential between NO2 molecules, which is a sum of an orientation-sensitive pairwise potential (OSPP) between N-N atoms proposed in Paper I [J. Chem. Phys. 115, 10852 (2001)] and Lennard-Jones potentials between N-O and O-O atoms, has been used in the simulation. The reaction dynamics is studied as a function of well depth De and anisotropy factors of the OSPP potential: Atheta (0< or =Atheta< or =1) for the rocking angle and Atau (0< or =Atau< or =0.5) for the torsional angle of relative NO2-NO2 orientation. The lifetime tauD of initially prepared NO2 dimers is found to increase as De increases, Atheta increases, and Atau decreases. Dissociation and association dynamics are studied in detail around the extreme limit of pure NO2-dimer liquid: De=0.12 x 10(-18) J, Atheta=0.5, and Atau=0.1, which has been found to reproduce both the observed liquid phase equilibrium properties and Raman band shapes of the dissociation mode very well. The dissociation dynamics from microscopic reaction trajectories is compared with the potential of the mean force (PMF) as a function of the N-N distance R. The PMF of reactive liquid NO2 shows a transition state barrier at R=2.3-2.5 A, and NO2-trimer structure is found to be formed at the barrier. Two types of dissociation of the NO2 dimer-the dissociation by collisional activation of the reactive mode to cross the dissociation limit and the NO2-mediated dissociation via bond transfer-are studied. The latter needs less free energy and is found to be much more probable. The dissociation trajectories and PMF in reactive liquid NO2 are compared with those of a reactive NO2 pair in inert solvent N2O4.