Next Article in Journal
Skin Sensitization Testing—What’s Next?
Next Article in Special Issue
Influence of Structure-modifying Agents in the Synthesis of Zr-doped SBA-15 Silica and Their Use as Catalysts in the Furfural Hydrogenation to Obtain High Value-added Products through the Meerwein-Ponndorf-Verley Reduction
Previous Article in Journal
Correction: Huang, F.-C., et al. Arabidopsis RETICULON-LIKE3 (RTNLB3) and RTNLB8 Participate in Agrobacterium-Mediated Plant Transformation. Int. J. Mol. Sci. 2018, 19, 638
Previous Article in Special Issue
Steam Reforming of Model Bio-Oil Aqueous Fraction Using Ni-(Cu, Co, Cr)/SBA-15 Catalysts
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 3
10.3390/ijms20030665
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Simulation of Naphthalene, Phenanthrene, and Pyrene Adsorption on MCM-41

by
Xiong Yang
1,2,†,
Chuanzhao Zhang
3,†,
Lijun Jiang
1,†,
Ziyi Li
1,2,*,
Yingshu Liu
1,2,
Haoyu Wang
3,
Yi Xing
1,2 and
Ralph T. Yang
4
1
School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
Beijing Higher Institution Engineering Research Center of Energy Conservation and Environmental Protection, Beijing 100083, China
3
College of Biochemical Engineering, Beijing Union University, Beijing 100023, China
4
Department of Chemical Engineering, University of Michigan, Ann Arbor, MI 48109, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2019, 20(3), 665; https://doi.org/10.3390/ijms20030665
Submission received: 28 December 2018 / Revised: 25 January 2019 / Accepted: 31 January 2019 / Published: 3 February 2019
(This article belongs to the Special Issue Mesoporous Materials: From Synthesis to Applications)
Browse Figures
Versions Notes

Abstract

:
The adsorption of three typical polycyclic aromatic hydrocarbons (PAHs), naphthalene, phenanthrene, and pyrene with different ring numbers, on a common mesoporous material (MCM-41) was simulated based on a well-validated model. The adsorption equilibriums (isotherms), states (angle distributions and density profiles), and interactions (radial distribution functions) of three PAHs within the mesopores were studied in detail. The results show that the simulated isotherms agreed with previous experimental results. Each of the PAHs with flat molecules showed an adsorption configuration that was parallel to the surface of the pore, in the following order according to the degree of arrangement: pyrene (Pyr) > phenanthrene (Phe) > naphthalene (Nap). In terms of the interaction forces, there were no hydrogen bonds or other strong polar forces between the PAHs and MCM-41, and the O–H bond on the adsorbent surface had a unique angle in relation to the PAH molecular plane. The polarities of different H atoms on the PAHs were roughly the same, while those of the C atoms on the PAHs decreased from the molecular centers to the edges. The increasing area of the π-electron plane on the PAHs with the increasing ring number could lead to stronger adsorption interactions, and thus a shorter distance between the adsorbate and the adsorbent.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) generally refer to hydrocarbons that have two or more benzene rings [1,2]. PAHs, commonly known as carcinogens, mutagens, teratogens, and bioaccumulators, can cause great damage in humans and the environment [3,4,5,6,7,8]. They are also a part of the atmospheric particulate pollutant PM2.5 [9]. In 1976, the United States Environmental Protection Agency (EPA)’s list of 129 “priority pollutants” included 16 types of PAH compounds [10]. Low-ring PAHs such as naphthalene (Nap), phenanthrene (Phe) and pyrene (Pyr) are the main components in the atmosphere [11]. The adsorption method has become one of the important ways to purify PAHs due to its high efficiency and low energy consumption [12,13,14]. Ordered mesoporous adsorbents have shown high adsorption capacities for their tailorable pore systems, especially in the adsorption of molecules such as PAHs [15,16]. Mastering the adsorption mechanism of PAHs on mesoporous materials is critical to optimize the material’s performance. Molecular simulation has obvious benefits in the adsorption isotherm calculation and adsorption state analysis of adsorption molecules. It can analyze the adsorption mechanism [17,18] and the nature of action in depth, and it has been widely used in adsorbent selection and design [19,20].
MCM-41 [21,22,23] is a typically ordered mesoporous adsorbent from a family of silicate and alumosilicate solids that has been employed to investigate the adsorption phenomena on various gases [24,25]. Since the emergence of MCM-41, a large number of adsorption-related material structures and physical and chemical properties have been analyzed [23,26]. Molecular simulations with MCM-41 have been conducted with the Monte Carlo (MC) method to calculate the various gas adsorption conditions [27,28,29,30,31]. The MCM-41 model evolved from earlier simple models [32,33,34,35,36] (e.g., single-layer cylindrical atomic wall, and the rough model of simple potential function), to later all-atom models [37,38] (e.g., carving out models in amorphous silica, and ensuring the bonding relationship between silicon and oxygen atoms), and recent templated synthesis models for restoring real synthetic conditions [39,40]. Among them, the carving-out modeling method is easier to implement, and calculations based on this model have been widely used for gas adsorption because of the good agreement with experimental results. Jing et al. [41] discussed the adsorption effects of two structural forms of unit cell models, and provided ideas for improving computational efficiency. Ugliengo et al. [42] optimized the MCM-41 structure by the density functional method. Pajzderska et al. [27] and Coasne et al. [43] studied the surface roughness, pore shape, and surface undulation of mesoporous channels, and made corresponding structural improvements. The above works have important significance in the improvement of the overall structure, surface morphology, and adsorption sites of the MCM-41 model.
With the development of molecular simulation technology, MCM-41 has made great progress in the field of adsorption separation as a typical mesoporous molecular sieve [28,29,44,45]. With regard to small molecule gas adsorption, Bonnaud et al. [46] characterized and analyzed the adsorption state and the diffusion phenomenon of water molecules in the pores of MCM-41 [46]. In the adsorption simulation of low-volatility compounds and macromolecular organics such as drugs, the adsorption characteristics of ibuprofen [47], tert-butanol [48], and neopentane [49] have been studied, as well as the influence of water molecules on adsorption [50,51]. Coasne et al. [30,31,52] studied the role of organic molecules such as benzene and toluene in the pores based on the combination of the angle of the benzene ring plane and the density profile. Compared with common small molecular gases, macromolecular organic compounds have more obvious molecular structures such as prominent geometric skeletons, polypolar sites, and deformations. It is beneficial to grasp the distribution of these properties in a spatial capacity to judge the state of the molecules in the pores and analyze the adsorption mechanism. At present, there are few simulation studies on the adsorption of PAHs. Analyses of the adsorption states and the adsorption mechanism are relatively lacking. Studying the adsorption of typical PAHs such as Nap, Phe, and Pyr on mesoporous adsorbents is helpful to understand the adsorption mechanism and its application in the removal of PAHs. Analysis of the state of molecules inside the adsorbent is extremely important for research and the improvement of adsorption efficiency.
In this paper, adsorption of the three PAHs Nap, Phe, and Pyr will be simulated based on the grand canonical Monte Carlo (GCMC) method. The potential adsorption mechanism will be discussed based on the adsorption states of the three molecules in the pores.

2. Results and Discussion

2.1. Adsorption Equlibrium

Figure 1 provides the GCMC simulation and experimental results [53,54] of the adsorption isotherms of Nap, Phe, and Pyr on MCM-41 at 398 K (normal flue gas temperature). The isotherms from the simulation calculation were in good agreement with the experimental values, with the following order according to the adsorption capacity: Pyr > Phe   > Nap, corresponding to the inverse sequence of the volatility of the PAH molecules. The results indicate that the current MCM-41 unit cell model with the calculation condition was reasonable and can truly reflect the adsorption of PAHs at low concentrations. The three isotherms of the PAHs belonged to Type I according to the IUPIC classification [55]. Because the PAH concentrations were low, all of the isotherms did not reach the saturated adsorption capacity. The adsorption capacity Q (mmmol·kg−1) can be determined with the formula:
Q = 1000 ∗ N P A H / M a d s o r b e n t ,
where N P A H is the number of PAHs in the pores of MCM-41 and M a d s o r b e n t (g·mol−1) is the molar mass of the adsorbent.
Figure 2 presents a diagram of PAH adsorption as an instantaneous map under different pressures, and it reflects the overall distribution state of PAH molecule adsorption in the cylindrical channel. PAH pore adsorption was mainly based on single-layer adsorption, and this also indicated that the isotherm belonged to Type I. It can also be observed that the PAHs were adsorbed on the surface of the mesoporous channel in a tiled manner. This will be discussed in the following section.

2.2. Adsorption State

2.2.1. Angle Distribution

It was important to explore the angular distribution of PAH adsorption in the pores since the PAHs exhibited a planar molecular structure. Thus, the adsorption state and configuration of PAHs on the pore surfaces was statistically obtained. Here, the angle between the plane normal vector of the adsorbed PAHs and the axis direction of the cylindrical channel was defined as θ (theta). Figure 3 provides the angle distributions of Nap, Phe, and Pyr adsorption on MCM-41. The θ of the three PAH molecules as they adsorbed in the pore was mainly distributed between 80 and 90 degrees, and the probability of θ between 0 and 15 degrees was close to zero. This indicates that there are no PAHs perpendicular to the axial direction of the channel; that is to say, the PAHs are mainly adsorbed on the surface of the mesoporous channel in a tiled manner, which is consistent with Figure 3. The degree of the parallel distribution (θ tends to 90°) of adsorbate over the channel surface increased with increasing PAH weight, indicating a more stable interaction between the adsorbate and adsorbent.
In order to explain the orientation of the adsorbate molecules in the pores of MCM-41, an order index, S, was defined as follows [30]:
S = ( 3 2 cos 2 θ − 1 2 )
where S tended to be −0.5 when the adsorbate molecules tended to be parallel to the surface, and S tended to be 1 when the adsorbate molecules tended to be perpendicular to the surface. The statistical results of the S values of the three PAHs are shown in Table 1. The S values of Phe and Pyr were close to −0.4, indicating that the two molecular planes were close to being parallel to the axial direction of the channel. The S value of Nap was only −0.2, indicating that the tiling property of the adsorption configuration on the surface was weak, and there was less molecular arrangement order compared to that of Phe and Pyr. The S values of the three PAHs did not reach the ideal value of −0.5, probably because the thermal motion of the molecule at a higher temperature of 398 K caused local perturbations of the molecule at the adsorbed position. In this study, the adsorbates were at low gas concentrations, and monolayer adsorption dominated. Three non-polar PAH molecules exhibited relatively weak physical adsorption properties on the surface. Therefore, the angle and position of their distribution may also be disturbed by the interaction of the surrounding adsorbate molecules.

2.2.2. Density Profile

The density profile can be used to infer the average spatial position of the molecules in the pores by overlapping the profiles of the single-peak type, symmetric double-peak type, and irregular type, which can be used to further analyze the adsorption. Figure 4 shows plots of the central density profile of the three PAHs adsorbed on MCM-41. The abscissa indicated the radial direction of the mesoporous cylindrical channel, and the origin represented the center of the channel. The maximum distance was almost equal to the radius of the pores in the MCM-41 model (around 19 Å) in our simulation; the further away the molecule was from the center of the pores, the closer it was to the edge of the pores. For the special structure of the three PAHs (two or more benzene ring structures), the density distributions of different geometric centers on Nap, Phe, and Pyr, as shown in Figure 5, were refined to further gauge the state of adsorbed molecules.
The molecules of Nap, Phe, and Pyr are stable planar molecules, and there are three special states in the adsorption process [49], which are in parallel with the molecular plane, perpendicular to the surface of the molecular sieve, and exist at a unique angle to the surface. These three different states correspond to three different interactions: the weak interaction between the surface of the adsorbent and the molecular plane, the strong interaction between the surface of the adsorbent and the polar edges of the molecules, and the combination of the two effects.
Figure 4 shows that the density distributions of all centers did not exhibit obvious symmetrical bimodal phenomena. The larger the molecular radius, the farther the distance of the density distribution peak corresponding to the center of the pore (the corresponding distances from the center of benzene rings to the central axis of the pore were 14.35 Å, 14.65 Å, and 14.75 Å), and the stronger the interaction between the molecule and the surface. The peak position relationship of the density profile between the adsorbate molecules Core and CoreX illustrates the following three points. First, with increasing adsorbate molecular size, there was a greater coincidence that the peak position of the density was located between the centers of the benzene ring and the center of the molecule. This corresponded to an increase in the stability of the molecular plane with the edge of the PAHs likely to be deflected by the surface polarity interaction. Second, the two sub-centers of the Nap molecule were partially offset from the peak position of the density at the center, indicating that there was a significant tilt angle between the Nap molecule and the surface of the pore. Third, the difference between the density peak positions of the centers of the benzene ring and the center of the Phe and Pyr molecules was small. It can be concluded from the statistical results of the previous angle distributions that the parallelism between the molecular planes and the surface of the pore was satisfactory.

2.3. Adsorption Interactions

The radial distribution function (RDF, g(r)) is a function of the distance between different particles (molecules or atoms). It is the ratio of the local density to the average density, which indicates the density variation of the particle in space. The RDF method was employed to further analyze the distribution distance between the adsorption sites of the PAHs (rings or atoms) and the active atoms on the adsorbent surface, such as O and H. The RDF between the PAH molecules and different sites of the mesoporous surface is shown in Figure 6. The distances corresponding to the RDF peaks are summarized in Table 2.
Figure 6(1) shows that the H atoms on the PAHs at different positions were the same distance from the O atoms on the surface of MCM-41. This indicated that the H atoms in the same PAH molecule had similar weak activities. The peaks of the RDF distances corresponding to Nap, Phe, and Pyr were 3.43 Å, 3.63 Å, and 3.53 Å, respectively (see Table 2), which were all larger than the common hydrogen bond length (≤3.10 Å). This indicated that there could be limited hydrogen bonds produced between PAH molecules and surface silanol groups (Si–OH).
Figure 6(2) shows there were large differences in the distances between the H atoms on MCM-41 and C atoms (HMCM41–CorePAHs) at different positions on the PAHs. The C atoms of the C–H bonds on the PAH were close to H atoms on MCM-41, and C atoms near the center of the PAH were spaced apart from the H atoms on MCM-41. This distribution and action was a result of the positive and negative charges of the molecule. The H and O atoms on the PAHs or the surface of MCM-41 had positive potential and negative potential electrical properties, respectively, while the electrical property of C atoms on the PAHs was dependent on the positions. The central C atom was positively charged while the outer C atom was negatively charged.
Table 2 shows that the difference in the distance between H and O atoms on MCM-41 and the RDF peak of the PAH centers was 0.44 Å (Nap), 0.4 Å (Phe), and 0.45 Å (Pyr). They were smaller than the average distance between H and O atoms in the surface silanol group (0.94 Å), indicating that the molecular planes of the O–H bond and the PAHs were not perpendicular, but existed at a unique angle. The spacing between the central PAH ring and the H and O atoms of the PAHs decreased with increasing PAH molecular size. The C atoms near the center of the PAHs, which were negatively charged, attracted Si and H atoms that were positively charged, and this led to a decrease in the distance between the PAHs and the inner surface of the adsorbent. This result was consistent with the density peak position results shown in Figure 4.
It can be concluded from the angular distribution and density peak position analysis that the three PAH molecules were adsorbed in a flat position on the inner surface of the adsorbent tunnel. The adsorbates possessed a stable planar structure with H atoms on the molecular edge, and there was no strong polar force generated during the adsorption process that can result in the severe tilt or deflection of the molecular plane. These findings were proved by RDF results, where the hydrogen bonding was weak at the edge of the PAHs from the large distance between H atoms on PAHs and O atoms on MCM-41. Additionally, there was no unique H atom, and the surface generated a prominent force that tilted the PAH plane for the same activity of H atoms at the edge of the PAHs. As we know, the current PAHs were formed by the combination of benzene rings, and large π bonds easily form throughout the molecule. The potential and charge density distributions of Nap, Phe, and Pyr obtained by the DFT simulations, shown in Figure 7, illustrated that the center of the ring on both sides of the PAH molecular plane was negatively charged (red part in Figure 7), which is a neighboring π bond, and the area of action of the π-electron plane increased when the number of rings increased. This led to enhanced, or stronger, adsorption when the adsorption process was governed by dispersion forces [56]. Hence, the heavier the molecular weight of the PAHs, the stronger the interaction with the surface of the pore and the shorter the distance between the adsorbate and adsorbent.

3. Methods and Parameters

3.1. Model Structures

An all-atom model of MCM-41 was established with the carving-out method [36,37] based on the rhombohedral cell model [57], as shown in Figure 8. The main mesoporous structure of the model consisted of excavated cylindrical holes in amorphous silica. Each unit contained four cylindrical channels. The atomic bonds in the whole unit cell were saturated by the connection of the oxygen (O) and hydrogen (H) atoms to the unsaturated bonds from silicon (Si) and O atoms (the Si atoms with three OH atoms were removed). The model was optimized and relaxed based on interactions described by the Dreiding force field. The ESP charges were calculated according to the method from Zhuo et al. [44] to finalize the structure of MCM-41. The model structure parameters are shown in Table S1, and comparisons of the characterization results with experimental values [52] are shown in Figure S1 in the Supplementary Material.
Nap, Phe, and Pyr are typical benzene-ring structures, as shown in Figure 9. The molecular models of the three PAHs were determined with the density function theory (DFT) method, with all-electron atomic charges obtained by the ESP method. In the DFT calculation, the Perdew–Wang method was employed to exchange correlation functions in the calculation, where the base group used was double numerical plus polarization (DNP). The distribution of the electrostatic potentials on the electron isopycnic surface of the three PAHs is shown in Figure 7.

3.2. Calculation Method

In the GCMC simulation, the adsorption temperature was fixed at 398 K, and the adsorbate molecules were regarded as rigid molecules. The Dreiding force field was used to describe the energy forms of the adsorbent–adsorbate and adsorbate–adsorbate interactions. The atomic force field form and interaction parameters are shown in Table S2. Electrostatic interactions are described in Ewald format with an accuracy of 10−4 kcal mol−1. The cutoff radii for the van der Waals interactions and hydrogen bonds were 15.5 Å and 4.5 Å, respectively. The total step number was 1 × 107, in which the first 5 × 106 steps were used to balance the adsorption state, and the next 5 × 106 steps were used to calculate the state average of the entire ensemble. The PAHs’ gas partial pressure in the simulation was set at a lower level (~100 ppm) according to the common concentration ranges of PAHs in actual situations [58]. The fugacity value required in the calculation process was obtained by converting the pressure into the physical parameters of Nap, Phe, and Pyr by the Peng–Robinson equation. The parameters such as critical temperature, critical pressure, and the eccentricity factor are shown in Table S3.
The adsorption model results under the pressure of 0.8 Pmax were used to further analyze the adsorption state of the three PAH molecules in the main channel of MCM-41. The adsorption states including the angle distribution and the density peak were analyzed. The radius radial distribution function representative of adsorbate–adsorbent mutual force was finally discussed.

4. Conclusions

Based on the GCMC simulation, the MCM-41 model was established, and the adsorption of Nap, Phe, and Pyr on MCM-41 was simulated. The adsorption model was validated by use of experimental adsorption isotherms. During the adsorption, no perpendicular state between the PAH molecular plane and the central axis of the adsorbent channel was observed. The order according to the degree of adsorption was: Pyr > Phe > Nap. The three PAHs were adsorbed in a flat position on the surface of the pore, and the configuration planes of Phe and Pyr were more parallel to the surface compared to Nap. The RDF curves indicated that there was no hydrogen bonding produced during the adsorption process. The H atoms in the same PAHs had the same activity, but the adsorption activity of the C atoms decreased from the central position, and the O–H bond on the adsorbent surface had a unique angle with the plane of the PAHs. The area of action of the π-electron plane increased when the number of rings increased, leading to enhanced, stronger adsorption, and thus, a shorter distance between the adsorbate and adsorbent.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/3/665/s1.

Author Contributions

Conceptualization, X.Y., Z.L. and R.T.Y.; methodology, C.Z., L.J. and Z.L.; writing, X.Y., L.J. and H.W.; project administration, Y.L. and Y.X.

Funding

This research was funded by National Natural Science Foundation of China (51604016, 51478083, 21676025), Beijing Natural Science Foundation (No. 8182019, 8174064), National Key R&D Program of China (No. 2017YFC0210302), the National Postdoctoral Program for Innovative Talents (No. BX201700029), the China Postdoctoral Science Foundation (No. 2018M630075), the Fundamental Research Funds for the Central Universities (No. FRF-TP-17-063A1, FRF-BD-18-015A).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

NapNaphthalene
PhePhenanthrene
PyrPyrene
GCMCGrand Canonical Monte Carlo

References

  1. Casarett, L.J.; Doull, J. Toxicology: The Basic Science of Poisons; McGraw-Hill: New York, NY, USA, 1996. [Google Scholar]
  2. Li, Z.; Liu, Y.; Yang, X.; Xing, Y.; Wang, Z.; Yang, Q.; Yang, R.T. Desorption kinetics of naphthalene and acenaphthene over two activated carbons via thermogravimetric analysis. Energy Fuels 2015, 29, 5303–5310. [Google Scholar] [CrossRef]
  3. Cheng, X.; Kan, A.T.; Tomson, M.B. Naphthalene adsorption and desorption from aqueous C60 fullerene. J. Chem. Eng. Data 2004, 49, 675–683. [Google Scholar] [CrossRef]
  4. Yang, K.U.N.; Wu, W.; Jing, Q.; Zhu, L. Aqueous adsorption of aniline, phenol, and their substitutes by multi-walled carbon nanotubes. Environ. Sci. Technol. 2008, 42, 7931–7936. [Google Scholar] [CrossRef] [PubMed]
  5. Owabor, C.N.; Ogbeide, S.E.; Susu, A.A. Adsorption and desorption kinetics of naphthalene, anthracene, and pyrene in soil matrix. Pet. Sci. Technol. 2010, 28, 504. [Google Scholar] [CrossRef]
  6. Long, A.S.; Lemieux, C.L.; Gagné, R.; Lambert, I.B.; White, P.A. Genetic Toxicity of Complex Mixtures of Polycyclic Aromatic Hydrocarbons: Evaluating Dose-Additivity in a Transgenic Mouse Model. Environ. Sci. Technol. 2017, 51, 8138–8148. [Google Scholar] [CrossRef] [PubMed]
  7. Liamin, M.; Le Mentec, H.; Evrard, B.; Huc, L.; Chalmel, F.; Boutet-Robinet, E.; Le Ferrec, E.; Sparfel, L. Genome-Wide Transcriptional and Functional Analysis of Human T Lymphocytes Treated with Benzo [α] pyrene. Int. J. Mol. Sci. 2018, 19, 3626. [Google Scholar] [CrossRef]
  8. Libalova, H.; Rossner, P.; Vrbova, K.; Brzicova, T.; Sikorova, J.; Vojtisek-Lom, M.; Beranek, V.; Klema, J.; Ciganek, M.; Neca, J.; et al. Comparative analysis of toxic responses of organic extracts from diesel and selected alternative fuels engine emissions in human lung BEAS-2B cells. Int. J. Mol. Sci. 2016, 17, 1833. [Google Scholar] [CrossRef]
  9. Mastral, A.M.; Callen, M.S. A review on polycyclic aromatic hydrocarbon (PAH) emissions from energy generation. Environ. Sci. Technol. 2000, 34, 3051–3057. [Google Scholar] [CrossRef]
  10. Keith, L.; Telliard, W. ES&T special report: Priority pollutants: Ia perspective view. Environ. Sci. Technol. 1979, 13, 416–423. [Google Scholar]
  11. Albuquerque, M.; Coutinho, M.; Borrego, C. Long-term monitoring and seasonal analysis of polycyclic aromatic hydrocarbons (PAHs) measured over a decade in the ambient air of Porto, Portugal. Sci. Total Environ. 2016, 543, 439–448. [Google Scholar] [CrossRef]
  12. Dowaidar, A.M.; El-Shahawi, M.S.; Ashour, I. Adsorption of Polycyclic Aromatic Hydrocarbons onto Activated Carbon from Non-Aqueous Media: 1. The Influence of the Organic Solvent Polarity. Sep. Sci. Technol. 2007, 42, 3609–3622. [Google Scholar] [CrossRef]
  13. Mastral, A.M.; Callen, M.; Murillo, R. Assessment of PAH emissions as a function of coal combustion variables. Fuel 1996, 75, 1533–1536. [Google Scholar] [CrossRef]
  14. Zhou, H.C.; Zhong, Z.P.; Jin, B.S.; Huang, Y.J.; Xiao, R. Experimental study on the removal of PAHs using in-duct activated carbon injection. Chemosphere 2005, 59, 861–869. [Google Scholar] [CrossRef] [PubMed]
  15. Li, Z.; Liu, Y.; Yang, X.; Xing, Y.; Tsai, C.; Wang, Z.; Yang, Q.; Yang, R.T. Desorption of polycyclic aromatic hydrocarbons on mesoporous sorbents: Thermogravimetric experiments and kinetics study. Ind. Eng. Chem. Res. 2016, 55, 1183–1191. [Google Scholar] [CrossRef]
  16. Araújo, R.S.; Azevedo, D.C.S.; Cavalcante, C.L., Jr.; Jiménez-López, A.; Rodríguez-Castellón, E. Adsorption of polycyclic aromatic hydrocarbons (PAHs) from isooctane solutions by mesoporous molecular sieves: Influence of the surface acidity. Microporous Mesoporous Mater. 2008, 108, 213–222. [Google Scholar] [CrossRef]
  17. Shang, J.; Li, G.; Singh, R.; Gu, Q.; Nairn, K.M.; Bastow, T.J.; Medhekar, N.; Doherty, C.M.; Hill, A.J.; Liu, J.Z.; et al. Discriminative separation of gases by a “molecular trapdoor” mechanism in chabazite zeolites. J. Am. Chem. Soc. 2012, 134, 19246–19253. [Google Scholar] [CrossRef] [PubMed]
  18. Shang, J.; Li, G.; Webley, P.A.; Liu, J.Z. A density functional theory study for the adsorption of various gases on a caesium-exchanged trapdoor chabazite. Comput. Mater. Sci. 2016, 122, 307–313. [Google Scholar] [CrossRef]
  19. Li, Z.; Liu, Y.; Zhang, C.; Yang, X.; Ren, J.; Jiang, L. Methane recovery from coal bed gas using modified activated carbons: A combined method for assessing the role of functional groups. Energy Fuels 2015, 29, 6858–6865. [Google Scholar] [CrossRef]
  20. Li, G.K.; Shang, J.; Gu, Q.; Awati, R.V.; Jensen, N.; Grant, A.; Zhang, X.; Sholl, D.S.; Liu, J.Z.; Webley, P.A.; et al. Temperature-regulated guest admission and release in microporous materials. Nat. Commun. 2017, 8, 15777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Beck, J.S.; Vartuli, J.C.; Roth, W.J.; Leonowicz, M.E.; Kresge, C.T.; Schmitt, K.D.; Chu, C.T.W.; Olson, D.H.; Sheppard, E.W.; McCullen, S.B.; et al. A new family of mesoporous molecular sieves prepared with liquid crystal templates. J. Am. Chem. Soc. 1992, 114, 10834–10843. [Google Scholar] [CrossRef]
  22. Kresge, C.T.; Leonowicz, M.E.; Roth, W.J.; Vartuli, J.C.; Beck, J.S. Ordered mesoporous molecular sieves synthesized by a liquid-crystal template mechanism. Nature 1992, 359, 710. [Google Scholar] [CrossRef]
  23. Selvam, P.; Bhatia, S.K.; Sonwane, C.G. Recent advances in processing and characterization of periodic mesoporous MCM-41 silicate molecular sieves. Ind. Eng. Chem. Res. 2001, 40, 3237–3261. [Google Scholar] [CrossRef]
  24. Corma, A. From microporous to mesoporous molecular sieve materials and their use in catalysis. Chem. Rev. 1997, 97, 2373–2420. [Google Scholar] [CrossRef] [PubMed]
  25. Ciesla, U.; Schüth, F. Ordered mesoporous materials. Microporous Mesoporous Mater. 1999, 27, 131–149. [Google Scholar] [CrossRef]
  26. Rimola, A.; Costa, D.; Sodupe, M.; Lambert, J.F.; Ugliengo, P. Silica surface features and their role in the adsorption of biomolecules: Computational modeling and experiments. Chem. Rev. 2013, 113, 4216–4313. [Google Scholar] [CrossRef]
  27. Pajzderska, A.; Gonzalez, M.A.; Mielcarek, J.; Wąsicki, J. Water behavior in MCM-41 as a function of pore filling and temperature studied by NMR and molecular dynamics simulations. J. Phys. Chem. C 2014, 118, 23701–23710. [Google Scholar] [CrossRef]
  28. Coasne, B.; Galarneau, A.; Di Renzo, F.; Pellenq, R.J. Molecular simulation of nitrogen adsorption in nanoporous silica. Langmuir 2010, 26, 10872–10881. [Google Scholar] [CrossRef]
  29. Chang, S.C.; Chien, S.Y.; Chen, C.L. Analyzing adsorption characteristics of CO2, N2 and H2O in MCM-41 silica by molecular simulation. Appl. Surf. Sci. 2015, 331, 225–233. [Google Scholar] [CrossRef]
  30. Coasne, B.; Galarneau, A.; Pellenq, R.J.M.; di Renzo, F. Adsorption, intrusion and freezing in porous silica: The view from the nanoscale. Chem. Soc. Rev. 2013, 42, 4141–4171. [Google Scholar] [CrossRef]
  31. Coasne, B.; Alba-Simionesco, C.; Audonnet, F.; Dosseh, G.; Gubbins, K.E. Adsorption and structure of benzene on silica surfaces and in nanopores. Langmuir 2009, 25, 10648–10659. [Google Scholar] [CrossRef]
  32. Ravikovitch, P.I.; Domhnaill, S.C.Ó.; Neimark, A.V.; Schueth, F.; Unger, K.K. Capillary hysteresis in nanopores: Theoretical and experimental studies of nitrogen adsorption on MCM-41. Langmuir 1995, 11, 4765–4772. [Google Scholar] [CrossRef]
  33. Maddox, M.W.; Gubbins, K.E. Molecular simulation of fluid adsorption in buckytubes and MCM-41. Int. J. Thermophys. 1994, 15, 1115–1123. [Google Scholar] [CrossRef]
  34. Maddox, M.W.; Olivier, J.P.; Gubbins, K.E. Characterization of MCM-41 using molecular simulation: Heterogeneity effects. Langmuir 1997, 13, 1737–1745. [Google Scholar] [CrossRef]
  35. Koh, C.A.; Montanari, T.; Nooney, R.I.; Tahir, S.F.; Westacott, R.E. Experimental and computer simulation studies of the removal of carbon dioxide from mixtures with methane using AlPO4-5 and MCM-41. Langmuir 1999, 15, 6043–6049. [Google Scholar] [CrossRef]
  36. He, Y.; Seaton, N.A. Heats of adsorption and adsorption heterogeneity for methane, ethane, and carbon dioxide in MCM-41. Langmuir 2006, 22, 1150–1155. [Google Scholar] [CrossRef] [PubMed]
  37. Brodka, A.; Zerda, T.W. Properties of liquid acetone in silica pores: Molecular dynamics simulation. J. Chem. Phys. 1996, 104, 6319–6326. [Google Scholar] [CrossRef]
  38. Pellenq RJ, M.; Levitz, P.E. Capillary condensation in a disordered mesoporous medium: A grand canonical Monte Carlo study. Mol. Phys. 2002, 100, 2059–2077. [Google Scholar] [CrossRef]
  39. Gelb, L.D.; Gubbins, K.E. Characterization of porous glasses: Simulation models, adsorption isotherms, and the Brunauer–Emmett–Teller analysis method. Langmuir 1998, 14, 2097–2111. [Google Scholar] [CrossRef]
  40. Pérez-Sánchez, G.; Gomes, J.R.B.; Jorge, M. Modeling self-assembly of silica/surfactant mesostructures in the templated synthesis of nanoporous solids. Langmuir 2013, 29, 2387–2396. [Google Scholar] [CrossRef]
  41. Jing, Y.; Wei, L.; Wang, Y.; Yu, Y. Molecular simulation of MCM-41: Structural properties and adsorption of CO2, N2 and flue gas. Chem. Eng. J. 2013, 220, 264–275. [Google Scholar] [CrossRef]
  42. Ugliengo, P.; Sodupe, M.; Musso, F.; Bush, I.J.; Orlando, R.; Dovesi, R. Realistic models of hydroxylated amorphous silica surfaces and MCM-41 mesoporous material simulated by large-scale periodic B3LYP calculations. Adv. Mater. 2008, 20, 4579–4583. [Google Scholar] [CrossRef]
  43. Coasne, B.; Pellenq, R.J.M. Grand canonical Monte Carlo simulation of argon adsorption at the surface of silica nanopores: Effect of pore size, pore morphology, and surface roughness. J. Chem. Phys. 2004, 120, 2913–2922. [Google Scholar] [CrossRef] [PubMed]
  44. Zhuo, S.; Huang, Y.; Hu, J.; Liu, H.; Hu, Y.; Jiang, J. Computer simulation for adsorption of CO2, N2 and flue gas in a mimetic MCM-41. J. Phys. Chem. C 2008, 112, 11295–11300. [Google Scholar] [CrossRef]
  45. Zhang, X.; Shao, X.; Wang, W.; Cao, D. Molecular modeling of selectivity of single-walled carbon nanotube and MCM-41 for separation of methane and carbon dioxide. Sep. Purif. Technol. 2010, 74, 280–287. [Google Scholar] [CrossRef]
  46. Bonnaud, P.A.; Coasne, B.; Pellenq, R.J.M. Molecular simulation of water confined in nanoporous silica. J. Phys. Condens. Matter 2010, 22, 284110. [Google Scholar] [CrossRef] [PubMed]
  47. Brás, A.R.; Fonseca, I.M.; Dionísio, M.; Schönhals, A.; Affouard, F.; Correia, N.T. Influence of nanoscale confinement on the molecular mobility of ibuprofen. J. Phys. Chem. C 2014, 118, 13857–13868. [Google Scholar] [CrossRef]
  48. Ghoufi, A.; Hureau, I.; Lefort, R.; Morineau, D. Hydrogen-bond-induced supermolecular assemblies in a nanoconfined tertiary alcohol. J. Phys. Chem. C 2011, 115, 17761–17767. [Google Scholar] [CrossRef]
  49. Herdes, C.; Ferreiro-Rangel, C.A.; Duren, T. Predicting neopentane isosteric enthalpy of adsorption at zero coverage in MCM-41. Langmuir 2011, 27, 6738–6743. [Google Scholar] [CrossRef]
  50. Lerbret, A.; Lelong, G.; Mason, P.E.; Saboungi, M.L.; Brady, J.W. Molecular dynamics and neutron scattering study of glucose solutions confined in MCM-41. J. Phys. Chem. B 2011, 115, 910–918. [Google Scholar] [CrossRef]
  51. Huang, J.; Meagher, M.M. Pervaporative recovery of n-butanol from aqueous solutions and ABE fermentation broth using thin-film silicalite-filled silicone composite membranes. J. Membr. Sci. 2001, 192, 231–242. [Google Scholar] [CrossRef]
  52. Coasne, B.; Fourkas, J.T. Structure and dynamics of benzene confined in silica nanopores. J. Phys. Chem. C 2011, 115, 15471–15479. [Google Scholar] [CrossRef]
  53. Li, Z.; Liu, Y.; Yang, X.; Xing, Y.; Yang, Q.; Yang, R.T. Adsorption thermodynamics and desorption properties of gaseous polycyclic aromatic hydrocarbons on mesoporous adsorbents. Adsorption 2017, 23, 361–371. [Google Scholar] [CrossRef]
  54. Li, Z.; Liu, Y.; Yang, X.; Xing, Y.; Tsai, C.J.; Meng, M.; Yang, R.T. Performance of mesoporous silicas and carbon in adsorptive removal of phenanthrene as a typical gaseous polycyclic aromatic hydrocarbon. Microporous Mesoporous Mater. 2017, 239, 9–18. [Google Scholar] [CrossRef]
  55. Yang, R.T. Gas Separation by Adsorption Processes; Butterworth-Heinemann: Oxford, UK, 2013. [Google Scholar]
  56. Moreno-Castilla, C. Adsorption of organic molecules from aqueous solutions on carbon materials. Carbon 2004, 42, 83–94. [Google Scholar] [CrossRef] [Green Version]
  57. Siperstein, F.R.; Gubbins, K.E. Synthesis and characterization of templated mesoporous materials using molecular simulation. Mol. Simul. 2001, 27, 339–352. [Google Scholar] [CrossRef]
  58. Liu, Y.; Li, Z.; Yang, X.; Xing, Y.; Tsai, C.; Yang, Q.; Wang, Z.; Yang, R.T. Performance of mesoporous silicas (MCM-41 and SBA-15) and carbon (CMK-3) in the removal of gas-phase naphthalene: Adsorption capacity, rate and regenerability. RSC Adv. 2016, 6, 21193–21203. [Google Scholar] [CrossRef]
Ijms 20 00665 g001 550
Figure 1. Adsorption isotherms for Nap, Phe, and Pyr on the MCM-41 model compared with experiment results at 398K. (a) Nap; (b) Phe; (c) Pyr; Key: red, grand canonical Monte Carlo (GCMC); black, experiment.
Figure 1. Adsorption isotherms for Nap, Phe, and Pyr on the MCM-41 model compared with experiment results at 398K. (a) Nap; (b) Phe; (c) Pyr; Key: red, grand canonical Monte Carlo (GCMC); black, experiment.
Ijms 20 00665 g001
Ijms 20 00665 g002 550
Figure 2. Typical configurations of Nap, Phe and Pyr adsorbed at 398K in the cylindrical nanopore of the model MCM-41. (a) Nap; (b) Phe; (c) Pyr, Key: Si, yellow; O, red; H, white; C, grey.
Figure 2. Typical configurations of Nap, Phe and Pyr adsorbed at 398K in the cylindrical nanopore of the model MCM-41. (a) Nap; (b) Phe; (c) Pyr, Key: Si, yellow; O, red; H, white; C, grey.
Ijms 20 00665 g002
Ijms 20 00665 g003 550
Figure 3. Angle distribution of Nap, Phe, and Pyr on the MCM-41 model. Key: black, Nap; red, Phe; blue, Pyr.
Figure 3. Angle distribution of Nap, Phe, and Pyr on the MCM-41 model. Key: black, Nap; red, Phe; blue, Pyr.
Ijms 20 00665 g003
Ijms 20 00665 g004 550
Figure 4. Density profiles for diffident centers in (a) Nap; (b) Phe; (c) Pyr at 398K on the MCM-41 model.
Figure 4. Density profiles for diffident centers in (a) Nap; (b) Phe; (c) Pyr at 398K on the MCM-41 model.
Ijms 20 00665 g004
Ijms 20 00665 g005 550
Figure 5. Graphical site description of polycyclic aromatic hydrocarbon (PAH) molecules. (a) naphthalene; (b) phenanthrene; (c) pyrene; Key: CoreX, the benzene ring center (X is the number); yellow, symbol of different atoms.
Figure 5. Graphical site description of polycyclic aromatic hydrocarbon (PAH) molecules. (a) naphthalene; (b) phenanthrene; (c) pyrene; Key: CoreX, the benzene ring center (X is the number); yellow, symbol of different atoms.
Ijms 20 00665 g005
Ijms 20 00665 g006 550
Figure 6. Radius distribution function for (a) Nap, (b) Phe, and (c) Pyr. (1) O atoms on MCM-41 and H atoms on PAHs, (2) H atoms on MCM-41 and C atoms on PAHs, (3) H atoms on MCM-41 and centers on PAHs, and (4) O atoms on MCM-41 and centers on PAHs.
Figure 6. Radius distribution function for (a) Nap, (b) Phe, and (c) Pyr. (1) O atoms on MCM-41 and H atoms on PAHs, (2) H atoms on MCM-41 and C atoms on PAHs, (3) H atoms on MCM-41 and centers on PAHs, and (4) O atoms on MCM-41 and centers on PAHs.
Ijms 20 00665 g006
Ijms 20 00665 g007 550
Figure 7. Electrostatic potential and charge density distribution results of Nap, Phe and Pyr. The electrostatic potential distribution ranged from −0.03 to 0.03, with the blue portion representing the positive potential and the red portion representing the negative potential.
Figure 7. Electrostatic potential and charge density distribution results of Nap, Phe and Pyr. The electrostatic potential distribution ranged from −0.03 to 0.03, with the blue portion representing the positive potential and the red portion representing the negative potential.
Ijms 20 00665 g007
Ijms 20 00665 g008 550
Figure 8. The final structure of the MCM-41 model. (a) Ball-and-stick model; (b) CPK model; silicon (Si), yellow; oxygen (O), red; hydrogen (H), white.
Figure 8. The final structure of the MCM-41 model. (a) Ball-and-stick model; (b) CPK model; silicon (Si), yellow; oxygen (O), red; hydrogen (H), white.
Ijms 20 00665 g008
Ijms 20 00665 g009 550
Figure 9. The structures of naphthalene (Nap), phenanthrene (Phe), and pyrene (Pyr) with ESP charges. (a) Nap; (b) Phe; (c) Pyr; Key: carbon (C), grey; H, white; charges, yellow.
Figure 9. The structures of naphthalene (Nap), phenanthrene (Phe), and pyrene (Pyr) with ESP charges. (a) Nap; (b) Phe; (c) Pyr; Key: carbon (C), grey; H, white; charges, yellow.
Ijms 20 00665 g009
Table
Table 1. Order parameter (S) for the molecular direction of Nap, Phe, and Pyr.
Table 1. Order parameter (S) for the molecular direction of Nap, Phe, and Pyr.
AdsorbateSS Values Corresponding States
Nap−0.18430disordered
Phe−0.35681vertical
Pyr−0.3629−0.5parallel
Table
Table 2. Distances between the H or O atom of MCM41 model and CoreX of the three PAHs.
Table 2. Distances between the H or O atom of MCM41 model and CoreX of the three PAHs.
TypeDistance (Å)
NapPhePyr
OMCM41–HPAHs3.433.633.53
HMCM41–CorePAHs4.484.434.38
OMCM41–CorePAHs4.924.834.83
OMCM41–HMCM410.94

Share and Cite

MDPI and ACS Style

Yang, X.; Zhang, C.; Jiang, L.; Li, Z.; Liu, Y.; Wang, H.; Xing, Y.; Yang, R.T. Molecular Simulation of Naphthalene, Phenanthrene, and Pyrene Adsorption on MCM-41. Int. J. Mol. Sci. 2019, 20, 665. https://doi.org/10.3390/ijms20030665

AMA Style

Yang X, Zhang C, Jiang L, Li Z, Liu Y, Wang H, Xing Y, Yang RT. Molecular Simulation of Naphthalene, Phenanthrene, and Pyrene Adsorption on MCM-41. International Journal of Molecular Sciences. 2019; 20(3):665. https://doi.org/10.3390/ijms20030665

Chicago/Turabian Style

Yang, Xiong, Chuanzhao Zhang, Lijun Jiang, Ziyi Li, Yingshu Liu, Haoyu Wang, Yi Xing, and Ralph T. Yang. 2019. "Molecular Simulation of Naphthalene, Phenanthrene, and Pyrene Adsorption on MCM-41" International Journal of Molecular Sciences 20, no. 3: 665. https://doi.org/10.3390/ijms20030665

APA Style

Yang, X., Zhang, C., Jiang, L., Li, Z., Liu, Y., Wang, H., Xing, Y., & Yang, R. T. (2019). Molecular Simulation of Naphthalene, Phenanthrene, and Pyrene Adsorption on MCM-41. International Journal of Molecular Sciences, 20(3), 665. https://doi.org/10.3390/ijms20030665

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop