- Open Access
Molecular mechanics and quantum mechanical modeling of hexane soot structure and interactions with pyrene
© The Royal Society of Chemistry and the Division of Geochemistry of the American Chemical Society 2000
- Received: 9 August 2000
- Accepted: 1 September 2000
- Published: 14 September 2000
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; ref. 1) were performed. The MM+ (N. L. Allinger, J. Am. Chem. Soc., 1977, 395, 157; ref. 2) and COMPASS (H. Sun, J. Phys. Chem., 1998, 102, 7338; ref. 3) 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.
- Polycyclic Aromatic Hydrocarbon
- Molecular Dynamic Simulation
- Force Field
- Interlayer Spacing
The importance of black carbon in the environment has been recognized for some time (see ref. 4 and references therein). Soot is one type of black carbon generated from combustion of fossil fuels. When soot particles form an aerosol, soot surfaces can act as a heterogeneous catalyst in a number of environmentally important atmospheric reactions. [5–7] When deposited in soils and sediments, soot has a significant influence on the transport and bioavailability of some organic contaminants due to their strong partitioning into soot. In addition, organic contaminants, such as polycyclic aromatic hydrocarbons (PAH), can form simultaneously with soot in the combustion process. Therefore, trace amounts of soot can play a disproportionate role in the long-term sequestration of otherwise biodegradable compounds. In addition, soot may represent the more hydrophobic end-member in the distributed reactivity model for soil organic matter of Weber and co-workers or the glassy/rubbery model of Pignatello and co-workers. Both of these models suggest a component of soil organic matter that has a high affinity for hydrophobic organic molecules with a rigid structure that can support voids large enough for molecules to diffuse into. Soot has these qualities although it certainly is not the only component of soil organic matter that may play a role in the long-term adsorption of organic contaminants.
Although the role of soot in the environment and the analytical chemistry of soot have been studied extensively, molecular simulations of soot structure and chemistry are rare. Some studies have modeled activated carbon, but the composition and structure of activated carbon are distinct from those of soot (see ref. 7 and references therein). The goal of this study is to discover if molecular modeling techniques with an appropriate force field will produce structural results consistent with observations of soot structure. If so, molecular modeling may provide a method for understanding the strong partitioning and irreversible adsorption of PAH into soot. (Note: 'irreversible' is defined as a process that does not undergo the back reaction via the same pathway as the forward reaction. This does not mean that PAH can not be desorbed once they are sorbed. For instance, if the sorbate changes the structure of the sorbent, then the desorption reaction will likely occur through a different mechanism and the sorption is 'irreversible'.)
The basic model used in this study was taken from the work of Smith and co-workers. This structure was generated from experimental studies on hexane soot. The model was derived from compositional analyses, spectroscopic analyses (e.g., FTIR, Raman, EPR and 13C CP/MAS NMR) and adsorption experiments. Consequently, the initial composition and structure used in this study should be an excellent approximation for hexane soot. The composition of this model is C489H160O34 (683 atoms). The structure used in this study is not meant to be a model for all soots, carbon black or char. Soots with different origins are expected to have different compositions and structures, so the results presented in this paper may not be applicable for a wide variety of substances commonly known as black carbon. The purpose of this paper is to determine which computational methods give reasonable results regarding soot structure and behavior. However, many of the physicochemical principles modeled herein may operate to a greater or lesser degree in various types of black carbon.
Molecular mechanics simulations were performed with two different force fields in two different molecular modeling programs. The MM+ force field was used within the program HyperChem 5.1; and the COMPASS force field was used with the program Cerius (see ref. 17). A major difference between these force fields is that the parameters within COMPASS were optimized to fit data of condensed phases; whereas MM+ parameters were optimized against gas-phase molecules. Thus, COMPASS should be more likely to give accurate results for a soot particle. Both sets of force fields were used within energy minimizations and molecular dynamics (MD) simulations. The main difference between these two methods is that energy minimizations occur at 0 K and result in a local energy minimum structure; whereas MD simulations include temperature and time and allow for a larger sampling of configuration space. This difference is critical when examining possible structures of large, complex molecules. The soot model was first subjected to energy minimization from the 2D structure, then MD simulations were carried out using the minimum energy structure as the starting configuration. MD simulations were performed with a time step of 1 × 10-15 s at 300 K for 10 000 to 100 000 time steps.
The 683 atom soot model represents a small fragment of a soot particle which can be an aggregate 0.1 to 1 μm in size made up of smaller particles on the order of 30 nm in diameter. Hence, the effect of model system size on the calculated results was investigated. A larger soot model was constructed from the energy-minimized structure of four 683 atom soot molecules. These models were randomly placed in proximity to one another and subjected to energy re-minimization and MD simulations with the COMPASS force field in Cerius (see ref. 17). The final diameter of this 2732 atom soot model is approximately 5 nm, closer to the actual soot particle diameter. Calculations on larger models with 30 nm diameters may be necessary in the future, however, in order to better characterize system-size effects on the soot model structure.
Polycyclic aromatic hydrocarbons (PAH) strongly partition into soots and may form simultaneously with soots. Hence, the interaction of the PAH pyrene (C16H10) with soot was modeled. Pyrene molecules were placed in various positions near the energy minimized soot model. Initial positioning of the four pyrene molecules was made to reflect π–π (London dispersion), H–π (dipole–dipole), H–O (H-bonding) and H–H (van der Waals) interactions between the pyrene and soot. Thus, the starting configuration was given the opportunity to find stable arrangements for a variety of bonding arrangements rather than pre-selecting only one type of bond. The pyrene–soot system was then subjected to energy minimizations with both the MM+ and COMPASS force fields. Following energy minimization, MD simulations of 10 000 to 100 000 time steps were carried out. The system size was then quadrupled as explained above. Energy minimizations and MD simulations (100 000 time steps) were then performed on the 2732 atom soot model with 16 pyrene molecules using only the COMPASS force field.
The above calculations were performed with the model system isolated in vacuo. To test the possible effects of solvation on the soot structure and partitioning of pyrene, MD simulations were performed with water molecules present. Runs of 100 000 time steps with the COMPASS force field in Cerius (ref. 17) were carried out on the 683 atom soot model including four pyrene molecules surrounded by 2000 water molecules in a nanodroplet (i.e., no periodic boundary conditions). The soot plus pyrene model system was taken from the previous MD simulation of this system in vacuo using the same force field.
Quantum mechanical calculations on fragments of the soot structure were performed with the programs GAUSSIAN 98 and DMol (ref. 21). Molecular mechanics force fields describe atom–atom interactions in terms of parameters derived for each atom type. Hence, bonding arrangements not included in the original parameterization may not be modeled adequately. Intermolecular bonding, such as H-bonding or π–π bonding, may be particularly problematic for analytical force fields because the parameters are typically optimized to reproduce properties of isolated molecules. For example, neither of the force fields used in this study contain specific H-bonding parameters. Instead, the force fields rely on interactions between partial charges to describe H-bonding even though H-bonding involves electron transfer. Molecular orbital (MO) or density functional theory (DFT) quantum mechanical calculations are generally more reliable than molecular mechanics results because the quantum methods solve for the electron density of a given molecular configuration. Thus, the quantum mechanical calculations are more flexible and accurate, especially when it comes to modeling unusual bonding or intermolecular forces. On the other hand, the quantum methods are much more computationally intensive than molecular mechanics. Consequently, it is impractical to model systems of hundreds or thousands of atoms with quantum mechanics. Similarly, long MD simulations of molecules are not practical with quantum methods. (Exceptionally demanding calculations are possible, however.) The solution to this problem in this study was to isolate fragments of the soot for quantum mechanical modeling. In this manner, long-range configurations are studied with molecular mechanics, and details of short-range bonding are studied with quantum mechanics.
Two tests were carried out. First, intermolecular H-bonding between phenolic and carboxylic anhydride functional groups thought to exist in hexane soot was computed using MP2/6-31G* calculations in GAUSSIAN 98. The 'MP2' indicates that electron correlation is handled with Møller-Plesset second-order perturbation theory. The "6-31G" indicates that 6 Gaussian functions were used to describe the electron density in each atomic orbital, and the valence orbitals were split into two sets – one of 3 Gaussians and one of 1 Gaussian. The "*" indicates that d-orbitals were included on all atoms except H to account for polarization. This energy minimization was performed to test the accuracy of intermolecular H-bonding in soot as predicted by MD simulations using the COMPASS force field. Second, DMol calculations using double-numeric basis sets with p-polarization functions and gradient-corrected exchange and correlation functionals[26, 27] were performed on a soot fragment and a pyrene molecule. The accuracy of this method should be similar to the MP2/6-31G* calculations above but computationally more efficient. These static energy calculations were used to test the reliability of interactions simulated between soot and pyrene with the COMPASS force field.
The initial hexane soot model structure consists of aromatic regions from 2 to 64 fused rings that are terminated with H atoms and carbonyl groups. The aromatic fragments are linked via ether linkages or intermolecular H-bonding. Sergides et al. constructed the 3D structure of hexane soot by layering the fragments of the 2D model. Layering is observed in high-resolution transmission electron microscopy (HRTEM) studies. The aromatic fragments of soot are similar to graphite, so this 3D structure is a reasonable approximation. However, the planar configuration of these soot models[1, 28] requires that the ether linkages have C–O–C angles of approximately 180°, whereas normal C–O–C bonds in ethers are nearly tetrahedral (e.g., 112°).
The soot structures predicted in MD simulations are dramatically different, depending on which force field was used. MD simulation of the 683 atom soot model for 10 000 time steps at 300 K using the MM+ force field in HyperChem results in little change from the energy minimized structure. (Note: MD simulations were terminated after l0 ps because configurational changes were not observed.) The resulting structure does not reproduce the layered structures observed in HRTEM analyses. Thus, the structure resulting from MM+ is unreliable. MD simulations with the COMPASS force field cause the structure to collapse rapidly. Fig. 1b shows the model soot structure generated with COMPASS-based MD simulations after 100 000 time steps. A layered structure is clearly evident between the sections of the soot model.
Density analyses of the energy-minimized soot model gives a density of approximately 1.36 g cm-3. The estimated density of the more compact, layered structure produced using COMPASS (Fig. 1b) is 1.37 g cm -3. Soot densities can be as low as 1.57 g cm-3. The lower density of the model structure may be partially an effect of the small model size because the 2732 atom model has a somewhat higher calculated density of 1.44 g cm -3. Larger model systems may be necessary to achieve the density of real soot particles.
Molecular orbital calculations on H-bonding within soot
To test whether this type of H-bond is realistic or just an artifact of the simulation, quantum mechanical calculations on a phenol–carboxylic anhydride dimer were carried out. An energy minimization using the MP2/6-31G* method in GAUSSIAN 98 was used to determine the H-bonding distance between these two molecules. The MP2/6-31G* method should give a reasonably accurate prediction in this case because the 6-31G* basis set does a reliable job of reproducing structures for these types of organic molecules, and MP2 electron correlation increases the accuracy of the H-bond description. The optimized structure in Fig. 5b shows that the predicted H-bond distance is 2.0 A, somewhat longer than the distance in the soot model. The 0.2 A difference is not dramatic, but somewhat surprising because the edge-on orientation of Fig. 5b should promote stronger H-bonding compared to the parallel orientation of the aromatic rings in Fig. 5a. Thus, one would expect the H-bond distance from the MP2/6-31G* calculations to be smaller rather than larger than that found in the MD simulations. One reason for this increased H-bonding in the MD simulation could be that π–π bonding determines the interlayer spacing such that the H-bonding is optimized within a constrained structure. This interpretation is consistent with the observation that the interlayer spacing does not change significantly in the soot simulation with water present, as discussed above.
Density functional theory calculations on soot fragment and pyrene
The observation that the COMPASS force field performs better for modeling soot than the MM+ force field is not surprising. A force field developed for modeling condensed phase systems should include a better representation of the long-range forces (e.g., van der Waals bonding, π–π bonding, H-bonding) operative in a solid phase like soot compared to a force field developed for gas-phase molecules. The point of this comparison is not to criticize the MM+ force field but to make the point that selection of a force field is of paramount importance in molecular simulations. The availability of user-friendly molecular modeling programs has made MD simulations possible for researchers whose main interest is not molecular simulation. However, the results of simulations will not be accurate unless the force field used to produce the results has been tested on similar systems. Another source of error is that simplifying a system to create a molecular model may eliminate several important physical or chemical factors involved with the process of interest. Thus, in this study, the simulations were geared to increase the level of complexity by increasing the system size and adding waters of solvation. Furthermore, the results were compared against experimental observations and tested by higher-level theory calculations. Perhaps most importantly, the initial model used for the simulations was generated by exhaustive experimental and analytical research.
The soot structure produced and its interactions with pyrene are consistent with experimental observations. This is not to say that more accurate agreement between simulation and experiment is not possible. On the contrary, this study suggests that larger model systems simulated for longer periods of time may produce more accurate soot densities and a more realistic structure. Higher temperature simulations may be necessary to anneal the soot model and overcome long relaxation times in the soot–water simulations. Similarly, longer simulations at elevated temperatures may be useful for studying PAH partitioning and diffusion within soots.
With the above limitations, the current simulation results help explain the observed PAH–soot partitioning and the sluggish desorption of PAH from soot. There are two sources of PAH in soots found in sediments: PAH co-formed with the soot and PAH partitioned into the soot from an aqueous phase. The former group is more strongly sequestered than the latter. Based on the position of pyrene within the model soot, the conclusion of this research is that the pyrene is part of the soot structure and does not expand the interlayer spacing. These structures were generated in MD simulations of soot and pyrene simultaneously, so they illustrate how a PAH could be incorporated with the normal soot structure. Removing the pyrene molecule not only requires that it diffuse through the tortuous micropores of the soot, but the energy lost due to formation of a large "vacancy" in the soot structure must also be regained. If there is no species diffusing into the soot to replace the PAH, then the soot structure must undergo considerable relaxation to anneal out this "vacancy". Under ambient temperature conditions found in soils and surface sediments, large-scale relaxation of the rigid soot structure does not seem likely. Hence, the bioavailability of PAH within soot, especially those PAH that co-formed with the soot, would be limited. Consequently, the environmental risk associated with these sequestered PAH would be decreased.
This research was supported by the Office of Naval Research. Computational resources were also supplied by the Aeronautical Systems Center and Space & Naval Warfare Systems Center through the DoD High Performance Computing initiative. The author would like to thank Louise Criscenti and two anonymous reviewers for improving the manuscript and Jonathan Mathews of the PSU Energy Institute for performing the density analysis calculations.
- Akhter MS, Chughtai AR, Smith DM: The structure of hexane soot I: Spectroscopic studies. Appl Spectrosc. 1985, 39: 143-10.1366/0003702854249114.View ArticleGoogle Scholar
- Allinger NL: Conformational analysis. 130. MM2. A hydro carbon force field utilizing V1 and V2 torsional terms. J Am Chem Soc. 1977, 99: 8127-10.1021/ja00467a001.View ArticleGoogle Scholar
- Sun H: COMPASS: An ab initio force-field optimized for condensed-phase applications – overview with details on alKane and benzene compounds. J Phys Chem. 1998, 102: 7338-View ArticleGoogle Scholar
- Goldberg ED: Black Carbon in the Environment. 1985, Wiley, New York, 198-Google Scholar
- Chughtai AR, Atteya MMO, Kim J, Konowalchuk BK, Smith DM: Adsorption and adsorbate interaction at soot particle surfaces. Carbon. 1998, 36: 1573-10.1016/S0008-6223(98)00116-X.View ArticleGoogle Scholar
- Amman M, Kalberer M, Baltensperger U: Heterogeneous production of nitrous acid on soot in polluted air masses. Nature. 1998, 395: 157-10.1038/25965.View ArticleGoogle Scholar
- Koehler BG, Nicholson VT, Roe HG, Whitney ES: A Fourier transform infrared study of the adsorption of n-hexane soot from -130 to -40°C. J Geophys Res. 1999, 104: 55070-Google Scholar
- Gustafsson , Gschwend PM: Soot as a strong partition medium for polycyclic aromatic hydrocarbons in aquatic systems. ACS Symp Ser. 1997, 671: 365-View ArticleGoogle Scholar
- Dobbins RA, Fletcher RA, Lu W: Laser microprobe analysis of soot precursor particles and carbonaceous soot. Combust Flame. 1995, 100: 301-10.1016/0010-2180(94)00047-V.View ArticleGoogle Scholar
- McGroddy SE, Farrington JW: Sediment porewater partitioning of polycyclic aromatic hydrocarbons in three cores from Boston Harbor, Massachusetts. Environ Sci Technol. 1995, 29: 1542-10.1021/es00006a016.View ArticleGoogle Scholar
- Huang W, Weber WJ: A distributed reactivity model for sorption by soils and sediments. 10. Relationships between desorption hysteresis and the diagenetic profiles of organic domains II and III. Environ Sci Technol. 1997, 31: 2562-10.1021/es960995e.View ArticleGoogle Scholar
- Xing B, Pignatello JJ: Dual-mode sorption of low-polarity compounds in glassy poly(vinyl chloride) and soil organic matter. Environ Sci Technol. 1997, 31: 792-10.1021/es960481f.View ArticleGoogle Scholar
- Muller EA, Gubbins KE: Molecular simulation study of hydrophilic and hydrophobic behavior of activated carbon surfaces. Carbon. 1998, 36: 1433-10.1016/S0008-6223(98)00135-3.View ArticleGoogle Scholar
- Smith DM, Chughtai AR: The surface structure and reactivity of black carbon. Colloids Surf A. 1995, 105: 47-10.1016/0927-7757(95)03337-1.View ArticleGoogle Scholar
- Katrinak KA, Rez P, Perkes PR, Buseck PR: Fractal geometry of carbonaceous aggregates from an urban aerosol. Environ Sci Technol. 1993, 27: 539-10.1021/es00040a013.View ArticleGoogle Scholar
- HyperChem, version 5.1. 1998, Hypercube, Inc., Gainesville, FLGoogle Scholar
- Cerius2, version 3.8. 1999, Molecular Simulations Inc. San Diego, CAGoogle Scholar
- Kubicki JD, Bleam WF: An Introduction to Molecular Modeling, in Molecular Modeling of Clays and Mineral Surfaces. Edited by: ed Kubicki JD, Bleam WF. 2000, CMS lectures, Clay Minerals Society,Google Scholar
- Frenkel D, Smit B: Understanding Molecular Simulation. 1995, Academic Press, San Diego, CA, 444-Google Scholar
- Frisch MJ, et al: GAUSSIAN 98, revision A.7. 1998, Gaussian, Inc., Pittsburgh, PAGoogle Scholar
- DMol3. 1998, Molecular Simulations Inc. San Diego, CAGoogle Scholar
- Reed AE, Curtiss LA, Weinhold F: Intermolecular interactions from a natural bond orbital, donor–acceptor view point. Chem Rev. 1988, 88: 899-10.1021/cr00088a005.View ArticleGoogle Scholar
- Hass KC, Schneider WF, Curioni A, Andreoni W: The chemistry of water on alumina surfaces: Reaction dynamics from first principles. Science. 1998, 282: 265-10.1126/science.282.5387.265.View ArticleGoogle Scholar
- Møller C, Plesset MS: Note on an approximation treatment for many-electron systems. Phys Rev. 1934, 46: 618-10.1103/PhysRev.46.618.View ArticleGoogle Scholar
- Foresman JB, Frisch Æ: Exploring Chemistry with Electronic Structure Methods: A Guide to Using Gaussian. 1996, Gaussian Inc. Pittsburgh, PA, 269-2Google Scholar
- Becke AD: Density functional thermochemistry. 3. The role of exact exchange. J Chem Phys. 1993, 98: 5648-10.1063/1.464913.View ArticleGoogle Scholar
- Lee CT, Yang WT, Parr RG: Development of the Colle–Salvetti correlation energy formula into a functional of the electron density. Phys Rev B. 1988, 37: 785-10.1103/PhysRevB.37.785.View ArticleGoogle Scholar
- Sergides CA, Jassim JA, Chughtai AR, Smith DM: Appl Spectrosc. 1987, 41: 482-10.1366/0003702874448805.View ArticleGoogle Scholar
- Palotás AB, Rainey LC, Feldermann CJ, Sarofim AF, Vander Sande JB: Soot morphology: An application of image analysis in high-resolution transmission electron microscopy. Microsc Res Tech. 1996, 33: 266-10.1002/(SICI)1097-0029(19960215)33:3<266::AID-JEMT4>3.0.CO;2-O.View ArticleGoogle Scholar
- Iino M: Network structure of coals and association behavior of coal-derived materials. Fuel Process Technol. 2000, 62: 89-10.1016/S0378-3820(99)00120-4.View ArticleGoogle Scholar
- Faulon JL, Mathews JP, Carlson GA, Hatcher PG: Correlation between micropore and fractal dimension of bitumi nous coal based on computer generated models. Energy Fuels. 1994, 8: 408-10.1021/ef00044a019.View ArticleGoogle Scholar
- Ruoff RS, Thornton T, Smith D: Density of fullerene containing soot as determined by helium pycnometry. Chem Phys Lett. 1991, 186: 456-10.1016/0009-2614(91)90208-Q.View ArticleGoogle Scholar
- Rockne KJ, Taghon GL, Kosson DS: Pore structure of soot deposits from several combustion sources. Chemosphere. 2000, 41: 1125-10.1016/S0045-6535(00)00040-0.View ArticleGoogle Scholar
- Palmer DC: CrystalMaker 4.0: A Crystal Structures Program for MacOS Computers. CrystalMaker Software, UK. 1998Google Scholar