In this work, we develop molecular models, or “force fields” for the explosives to predict their octanol-water partition coefficients, Henry's law constants and also vapor-liquid equilibria, vapor pressure and critical points. For the non-bonded interactions, Lennard-Jones parameters for each interaction site are obtained from the TraPPE-UA [2,3,4] and TraPPE-EH force fields [5] from similar functional groups. The partial charges are fitted to reproduce partition coefficients and Henry's law constants from the experiment. NPT molecular dynamics simulations coupled with the free energy perturbation technique is used to compute the physicochemical properties and Gibbs-Duhem integration is used to determine the vapor-liquid equilibria, vapor pressure and critical parameters. The log Kow values predicted for DNAN and MNA are 1.648 and 1.934, respectively and are within 3 % of experiment.
[1] J.W.Gillett, Environmental Toxicology and Chemistry. 2, 463 (1983).
[2] C.D.Wick, J.M.Stubbs, N.Rai, J.I.Siepmann, J.Phys.Chem.B, 109 (18974).
[3] M.G.Martin, J.I.Siepmann, J.Phys.Chem.B, 104 (8008).
[4] J.M.Stubbs, J.J.Potoff, J.I.Siepmann, J.Phys.Chem.B, 108 (17596).
[5] N.Rai, J.I.Siepmann, J.Phys.Chem.B, 111 (2007).