Abstract

Molecular simulations (energy minimizations and molecular dynamics) of an n-hexane soot model developed by Smith and co-workers (M. S. Akhter, A. R. Chughtai and D. M. Smith, Appl. Spectrosc., 1985, 39, 143; ) were performed. The MM+ (N. L. Allinger, J. Am. Chem. Soc., 1977, 395, 157; ) and COMPASS (H. Sun, J. Phys. Chem., 1998, 102, 7338; ) force fields were tested for their ability to produce realistic soot nanoparticle structure. The interaction of pyrene with the model soot was simulated. Quantum mechanical calculations on smaller soot fragments were carried out. Starting from an initial 2D structure, energy minimizations are not able to produce the observed layering within soot with either force field. Results of molecular dynamics simulations indicate that the COMPASS force field does a reasonably accurate job of reproducing observations of soot structure. Increasing the system size from a 683 to a 2732 atom soot model does not have a significant effect on predicted structures. Neither does the addition of water molecules surrounding the soot model. Pyrene fits within the soot structure without disrupting the interlayer spacing. Polycyclic aromatic hydrocarbons (PAH), such as pyrene, may strongly partition into soot and have slow desorption kinetics because the PAH-soot bonding is similar to soot–soot interactions. Diffusion of PAH into soot micropores may allow the PAH to be irreversibly adsorbed and sequestered so that they partition slowly back into an aqueous phase causing dis-equilibrium between soil organic matter and porewater.