A potential energy function with harmonic intramolecular and Lennard-Jones plus Coulombic intermolecular
terms was tested in molecular dynamics simulations of liquid nitromethane. Parameter values were adjusted
iteratively until satisfactory agreement with density functional pair calculations and experimental data was
achieved. The properties computed using the NVT and NPT ensembles were the heat of vaporization, dielectric
constant, self-diffusion coefficient, density, heat capacity at constant pressure, pair correlation functions, single
molecule and collective dipole moment reorientation times, the vibrational spectrum, and the effect of increasing
pressure upon the C−N stretching frequency. Overall, the results were in reasonable accord with experimental
results, the greatest discrepancy being for the dielectric constant. It was concluded, on the basis of the
reorientation times and the calculated molecular surface electrostatic potential, that the intermolecular
interactions in liquid nitromethane at 1 atm are not highly directional and site-specific.