Determination of the experimental equilibrium structure of solid nitromethane using path-integral molecular dynamics simulations

Abstract

Path-integral molecular dynamics (PIMD) simulations with an empirical interaction potential have been used to determine the experimental equilibrium structure of solid nitromethane at 4.2 and 15 K. By comparing the time-averaged molecular structure determined in a PIMD simulation to the calculated minimum-energy (zero-temperature) molecular structure, we have derived structural corrections that describe the effects of thermal motion. These corrections were subsequently used to determine the equilibrium structure of nitromethane from the experimental time-averaged structure. We find that the corrections to the intramolecular and intermolecular bond distances, as well as to the torsion angles, are quite significant, particularly for those atoms participating in the anharmonic motion of the methyl group. Our results demonstrate that simple harmonic models of thermal motion may not be sufficiently accurate, even at low temperatures, while molecular simulations employing more realistic potential-energy surfaces can provide important insight into the role and magnitude of anharmonic atomic motions.