Parametrization and Molecular Dynamics Simulations of Nitrogen Oxyanions and Oxyacids for Applications in Atmospheric and Biomolecular Sciences

Nitrogen oxyanions
and oxyacids are important agents in atmospheric
chemistry and medical biology. Although their chemical behavior in
solution is relatively well understood, they may behave very differently
at the water/air interface of atmospheric aerosols or at the membrane/water
interface of cells. Here, we developed a fully classical model for
molecular dynamics simulations of NO3–, NO2–, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field
versions. The model successfully accounted for the poorly structured
solvation shell and ion pairing tendency of NO3–. Accurate pure-liquid properties and hydration free energies were
obtained for the oxyacids. Simulations at the water/air interface
showed a local enrichment of HNO3 and depletion of NO3–. The effect was discussed in light of
earlier spectroscopic data and ab initio calculations, suggesting
that HNO3 behaves as a weaker acid at the surface of water.
Our model will hopefully allow for efficient and accurate simulations
of nitrogen oxyanions and oxyacids in solution and at microheterogeneous
interface environments.