Molecular nitrogen-N2 properties: The intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].