Next Article in Journal
Phytochemistry, Bioactivities and Traditional Uses of Michelia × alba
Next Article in Special Issue
Synthesis and Biomedical Applications of Highly Porous Metal–Organic Frameworks
Previous Article in Journal
Diclofenac Diminished the Unfolded Protein Response (UPR) Induced by Tunicamycin in Human Endothelial Cells
Previous Article in Special Issue
Linker Functionalization Strategy for Water Adsorption in Metal–Organic Frameworks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Surface Modification Strategy for Enhanced NO2 Capture in Metal–Organic Frameworks

by
Dionysios Raptis
1,
Charalampos Livas
1,
George Stavroglou
1,
Rafaela Maria Giappa
1,2,
Emmanuel Tylianakis
2,
Taxiarchis Stergiannakos
1 and
George E. Froudakis
1,*
1
Department of Chemistry, University of Crete, Voutes Campus, GR-71003 Heraklion, Crete, Greece
2
Department of Materials Science and Technology, University of Crete, Voutes Campus, GR-71003 Heraklion, Crete, Greece
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(11), 3448; https://doi.org/10.3390/molecules27113448
Submission received: 19 April 2022 / Revised: 23 May 2022 / Accepted: 24 May 2022 / Published: 26 May 2022

Abstract

:
The interaction strength of nitrogen dioxide (NO2) with a set of 43 functionalized benzene molecules was investigated by performing density functional theory (DFT) calculations. The functional groups under study were strategically selected as potential modifications of the organic linker of existing metal–organic frameworks (MOFs) in order to enhance their uptake of NO2 molecules. Among the functional groups considered, the highest interaction energy with NO2 (5.4 kcal/mol) was found for phenyl hydrogen sulfate (-OSO3H) at the RI-DSD-BLYP/def2-TZVPP level of theory—an interaction almost three times larger than the corresponding binding energy for non-functionalized benzene (2.0 kcal/mol). The groups with the strongest NO2 interactions (-OSO3H, -PO3H2, -OPO3H2) were selected for functionalizing the linker of IRMOF-8 and investigating the trend in their NO2 uptake capacities with grand canonical Monte Carlo (GCMC) simulations at ambient temperature for a wide pressure range. The predicted isotherms show a profound enhancement of the NO2 uptake with the introduction of the strongly-binding functional groups in the framework, rendering them promising modification candidates for improving the NO2 uptake performance not only in MOFs but also in various other porous materials.

1. Introduction

Nitrogen dioxide (NO2) belongs to a group of highly reactive gases known as nitrogen oxides (NOx) and primarily gets in the air from the burning of fuels [1]. The NO2 generated by the exhaust gases of the industries as well as by our private cars become an important air pollutant, the toxicity of which has an impact on the environment and human health [2]. NO2 is formed in the combustion processes of heating systems. The main pollutant emitted directly from hydrocarbon combustion is nitric oxide (NO), along with a small proportion of nitrogen dioxide (NO2). Nitrogen oxide is oxidized by ozone (O3) in the atmosphere, on a 10 min time scale, to give NO2 [3]. According to the EPA (Environmental Protection Agency), when NO2 is present in the air, it can be harmful to human health due to the irritation created in the airways of the human respiratory system [4]. In addition, NO2 and other NOx molecules can interact with water, oxygen, and other chemicals in the atmosphere to form acid rain [5,6], whose harmful ecological effects are detrimental and most prominent in aquatic environments. Due to the serious impacts of nitrogen pollution [2], there are several separation methods of NO2 from industrial gas mixtures [7] and various adsorbents which have been used for the removal of NO2; some of them including zeolites [8], wood-based activated carbon [9], graphite oxides and iron composites [10]. Zhu et al., in an attempt to adsorb NO2 in N-doping activated carbon, observed that N-doping leads to easier adsorption of NO2 molecules, thus increasing NO2 physisorption energies [11]. After theoretically investigating NO2 adsorption on the surface of a silicon carbide (SiC) nanotube, Iranimanesh et al. proposed it for NO2 gas pollutant sensing and removal [12]. Lately, the separation of NO2 from industrial gas mixtures through its adsorption has been clearly identified as a possible application of metal–organic frameworks (MOFs) [13,14].
In the last two decades, there has been a rapid growth of metal–organic frameworks (MOFs) in the field of porous materials, and their applications vary from adsorption and separation of gases [15] to catalysis [16] and drug delivery [17]. MOFs are crystalline materials with extremely high porosity (up to 90% free volume) and a huge internal surface area extending beyond 6000 m2/g [18]. Due to these properties and considering the plethora of both inorganic and organic building blocks, MOFs have attracted great research interest as high-capacity adsorbents to meet various separation demands [18,19,20]. One of the important structure-to-property flexibilities of the MOF structures is the potential of their organic linkers to be modified with the incorporation of various functional groups in order to tune their interaction with selected molecules [18].
Several previous studies have shown that the introduction of functional groups into MOFs leads to enhanced gas uptake performance. In the work of Frysali et al. [21], the introduction of a sulfate anion in the phenyl ring has the highest interaction energy (−5.4 kcal/mol) with CO2, a value almost two times larger than the corresponding binding energy for benzene (−2.9 kcal/mol). Klontzas et al. [22] showed that the gravimetric capacity of the Li modified IRMOF-8 was calculated to 10 wt% at 77 K and 100 bar, while the corresponding values show great promise also at room temperature with an uptake of 4.5 wt%, performances significantly enhanced with respect to the unmodified counterpart (up to three times stronger). For NO2, the theoretical studies made by Fioretos et al. [23] have shown that NO2 interacts stronger with functionalized benzenes such as aniline, phenol, and toluene (with binding energies of −2.26 kcal/mol, −1.72 kcal/mol and −2.02 kcal/mol respectively) than with benzene (−1.67 kcal/mol). The amino (-NH2) substituent can be particularly beneficial as the interaction of strongly polar molecules, such as NO2, with the amino-substituted aromatic rings is characterized by the contribution of electrostatic dipole–dipole forces resulting in enhanced adsorption.
Taking into account that many MOF frameworks have a phenyl group in their organic linker, together with the effectiveness of the organic linker functionalization strategy for tuning their interaction with guest molecules, in this work, we investigate the interaction of NO2 molecules with a series of 43 carefully selected benzene molecules by means of density functional theory (DFT) calculations. The selection of the functional groups was based on chemical intuition and findings of previous similar studies [22,24,25,26,27]. Subsequently, in order to verify the effectiveness of the functionalization on the enhancement of NO2 uptake in MOFs, grand canonical Monte Carlo (GCMC) simulations were performed for the IRMOF-8, modified by the three functional groups that showed the strongest interaction with NO2.

2. Computational Methods

2.1. Density Functional Theory

To investigate the interaction of NO2 molecules with the organic linkers of MOFs, we start with the simplified model of a benzene ring. A large set of 43 functional groups was examined for their binding strength towards the NO2 molecule. The conformations of the functionalized benzene molecules were optimized using Gershom Martin’s double-hybrid density functional DSD BLYP in the resolution of identity (RI) approximation [28] along with the def2-TZVPP basis set and with the corresponding auxiliary basis set for the RI approximation [29,30,31], including the D3BJ (Becke–Johnson damping version) empirical correction term for the dispersion interactions as proposed by Grimme [32,33,34,35]. All geometries were optimized without any symmetry constraints, and the optimized minimum-energy structures were verified as stationary points on the potential energy surface by performing numerical harmonic vibrational frequency calculations. All calculations were performed with the Orca 4.2 program package [36,37]. The dimer energies were corrected for the basis set superposition error (BSSE) using the counterpoise (CP) method as proposed by Boys and Bernardi [38].
The electron density redistribution was calculated as the difference between the electron density of the functionalized benzene–NO2 (ΔD) complex and the electron densities of the functionalized benzene (D(functionalized benzene)) and NO2 (D(NO2)) molecule according to the formula:
ΔD = D(functionalized benzene-NO2) − D(functionalized benzene) − D(NO2)
All electron densities were calculated at the RI-DSD-BLYP/def2-TZVPP level of theory. Mathematical operations on the electron densities along with the visualization of the electron density difference, were done using gOpenMol graphics program [39,40].

2.2. Grand Canonical Monte Carlo

To verify the effectiveness of the strongest interacting functional group candidates obtained from the DFT calculations to enhance the NO2 uptake in MOFs, we employed Monte Carlo simulations in the grand canonical ensemble. IRMOF-8 was selected as the platform for the uptake calculation and was functionalized with the best performing functional groups, as shown in Figure 1.
The ΝO2 adsorption was calculated for a pressure range up to 1.2 bar at 298 K. Fugacity coefficients for the different thermodynamic states were defined using the Penge–Robinson equation of state [41]. The framework coordinates were taken from the crystallographic data [42] and a cubic periodic box of size 30.1 × 30.1 × 30.1 Å3 was used for all the frameworks, functionalized or parent. Simulations were performed in supercells incorporating enough repeat units such that all edge lengths were greater than 25.6 Å, i.e., twice the Lennard-Jones (LJ) cut-off radius. For each simulation point, 50,000 cycles were performed for system equilibration, followed by additional 100,000 cycles for sampling over the ensemble averages. For the description of the interactions between the IRMOF-8 and the NO2 atoms, LJ + Coulomb potentials were used and each atom of the host or the guest was treated explicitly [43]. The framework of IRMOF-8 was kept rigid during the simulations, while NO2 molecules were allowed to translate and rotate. Nitrogen dioxides were treated as three rigid center molecules with bond lengths between nitrogen and oxygen atoms held fixed at 1.19 Å. For the electrostatic interactions between NO2 molecules and the host material, point charges equal to qO = −0.073 and qN = 0.146 were placed at oxygen and nitrogen sites of NO2, respectively. The framework atoms charges were defined by employing the CHELPG method [44]. For the van der Waals interactions, potential parameters according to COMPASS Force Field [45,46] model were used, with ε = 50.36 and σ = 3.24 Å for Nitrogen atom and ε = 62.51 and σ = 2.93 Å for the oxygen center.
For each MOF framework, the potential parameters were taken from the UFF force field [47], except for organic linker atoms that were treated separately. Lorenz–Berthelot mixing rules were used to describe the NO2−IRMOF-8 interactions. For the functionalized organic linkers, the UFF parameters were found to be inconsistent with DFT interactions. To correct this, the parameters of the classical potential were fitted to reproduce the quantum chemical data (SI). All GCMC calculations were carried out with the RASPA software package [48].

3. Results and Discussion

The interaction of the NO2 molecule with the full set of the 43 strategically functionalized benzenes was calculated at the RI-DSD-BLYP/def2-TZVPP level of theory and can be seen in Table S1 of the SI section. In Table 1, we present the nine functional group candidates with the highest binding energy together with the non-functionalized benzene for comparison. In Figure 2, the DFT optimized dimer geometries of the NO2…C6H5-X systems are shown.
The highest interaction energy (5.4 kcal/mol) with NO2 was found for phenyl hydrogen sulfate (-OSO3H), which is almost three times stronger than the corresponding binding energy of the unfunctionalized benzene (2.0 kcal/mol).
The energetically most favorable structures are characterized by a weak interaction between NO2 and organic linkers. More specifically, the acidic protons of the substituents, especially (-OSO3H, -OPO3H2, -PO3H2, -SOOH, -C(OH)3), tend to interact with NO2′s oxygen site with binding distances between 1.89 Å and 2.22 Å. The trend in the best binding energies is also confirmed by the electron density redistribution plots in Figure 3. Due to the electrostatic nature of NO2′s interactions, these plots can serve as a rule of thumb for predicting the most stable NO2 complexes with organic molecules. The red regions of the electron redistribution plots correspond to rich electron areas located around oxygen atoms, where the blue regions correspond to poor electron areas located mainly around hydrogen atoms. The electron-rich region of the nitrogen atom of NO2 interacts with the electron-poor regions of the functionalized benzene, and the electron-rich regions around oxygen atoms of NO2 interact with the electron-poor regions (around hydrogen atoms) of the functionalized benzenes.
To test the effect of the proposed surface modification strategy on the enhancement of the NO2 capture in MOF structures, we selected IRMOF-08 as the platform and modified its framework by adding one functional group per linker, as shown in Figure 1. We calculated the NO2 adsorption isotherms by performing GCMC simulations for the three functional groups (-OSO3H, -PO3H2, -OPO3H2) with the stronger interactions. From the excess volumetric and gravimetric isotherms at 298 K and pressures up to 1.2 bar (Figure 4a,b), there is a significant enhancement of the uptake due to functionalization. At 1.0 bar, the corresponding volumetric and gravimetric uptake for the IRMOF-8 with the OSO3H group was found at 193 cm3 (STP)/cm3 and 15 mmol/g, respectively, a value 12.8% larger than that of the unmodified IRMOF-8. The corresponding volumetric and gravimetric uptake for the unmodified IRMOF-8 was found at 7.5 cm3 (STP)/cm3 and 0.8 mmol/g.
Figure 5 shows representative snapshots taken at 0.1 and 1.0 bar for the parent and the three functionalized frameworks. In both pressures, the modified material hosts considerably more NO2 molecules than the parent structure due to the stronger binding sites introduced to the structure by functionalization. This is also verified by the fact that the NO2 molecules are located closer to the functional groups.

4. Conclusions

In this work, we studied by means of density functional theory the interaction of a chemically diverse set of 43 functionalized benzenes with the NO2 molecule.
The highest interaction energy with NO2 (5.4 kcal/mol at the RI-DSD-BLYP/def2-TZVPP level of theory) was found for phenyl hydrogen sulfate (-OSO3H)—an interaction almost three times larger than the corresponding binding energy for non-functionalized benzene (2.0 kcal/mol).
The groups with the highest NO2 binding (-OSO3H, -PO3H2, -OPO3H2) were selected for functionalizing the linker of IRMOF-8 and investigating the trend in their NO2 uptake capacities with grand canonical Monte Carlo (GCMC) simulations. GCMC simulations showed a clear enhancement of the NO2 uptake both gravimetrically and volumetrically at 298 K and pressures up to 1.2 bar for the functionalized MOF, an enhancement even more pronounced at low pressures; at 0.1 bar, the volumetric uptake becomes 40 or 60 or 110 times larger than the unmodified IRMOF-08 by introducing -OPO3H2, -OSO3H, -PO3H2 functional groups, respectively. Based on this significant enhancement, we propose our surface functionalization as a general strategy for improving the NO2 adsorption uptake not only in MOFs, but also in various other porous materials. Our theoretical results can serve as high accuracy reference calculations and guide synthesis towards materials with high NO2 adsorption capacity.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27113448/s1, Table S1: Sorted binding energies (kcal/mol) of the NO2…C6H5-X systems under study, calculated at the RI-DSD-BLYP D3(BJ)/def2-TZVPP level of theory. All interaction energy values have been corrected for the basis set superposition error (BSSE) by the full counterpoise method [38]. Percentage of binding energy enhancement with the introduction of the FG compared to benzene; Figure S1: Global minima geometries and binding energy values (in kcal/mol) of all the systems in this study; Figure S2: Electron density redistribution plots of the optimized geometries of the NO2…C6H5-X complexes. With red and blue being the regions that gain and lose electron density upon the formation of the complex, respectively; Figure S3: Fitting of the (ε, σ) parameters of the UFF [47] potential on the QM data obtained from the ab-initio scan of NO2 over benzene; Figure S4: Fitting of the (ε, σ) parameters for the NO2…C6H5-OPO3H2 interaction; Figure S5: Fitting of the (ε, σ) parameters for the NO2…C6H5-OSO3H interaction; Figure S6: Fitting of the (ε, σ) parameters for the NO2…C6H5-PO3H2 interaction.

Author Contributions

Conceptualization and methodology, G.E.F.; investigation, D.R., C.L., G.S., R.M.G., E.T. and T.S.; writing, editing, and formal analysis, D.R., R.M.G., T.S. and G.E.F.; All authors have read and agreed to the published version of the manuscript.

Funding

This research has been co-financed by the European Regional Development Fund of the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship, and Innovation, under the call RESEARCH—CREATE—INNOVATE (project code: T2EΔK-01976).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hill, S.C.; Douglas Smoot, L. Modeling of Nitrogen Oxides Formation and Destruction in Combustion Systems. Prog. Energy Combust. Sci. 2000, 26, 417–458. [Google Scholar] [CrossRef]
  2. Camargo, J.A.; Alonso, Á. Ecological and Toxicological Effects of Inorganic Nitrogen Pollution in Aquatic Ecosystems: A Global Assessment. Environ. Int. 2006, 32, 831–849. [Google Scholar] [CrossRef] [PubMed]
  3. Han, S.; Bian, H.; Feng, Y.; Liu, A.; Li, X.; Zeng, F.; Zhang, X. Analysis of the Relationship between O3, NO and NO2 in Tianjin, China. Aerosol Air Qual. Res. 2011, 11, 128–139. [Google Scholar] [CrossRef] [Green Version]
  4. EPA. U.S. Integrated Science Assessment for Oxides of Nitrogen–Health Criteria; US Environmental Protection Agency: Washington, DC, USA, 2016.
  5. Samet, J.M.; Utell, M.J. The Risk of Nitrogen Dioxide: What Have We Learned from Epidemiological and Clinical Studies? Toxicol. Ind. Health 1990, 6, 247–262. [Google Scholar] [CrossRef]
  6. Hesterberg, T.W.; Bunn, W.B.; McClellan, R.O.; Hamade, A.K.; Long, C.M.; Valberg, P.A. Critical Review of the Human Data on Short-Term Nitrogen Dioxide (NO2) Exposures: Evidence for NO2 No-Effect Levels. Crit. Rev. Toxicol. 2009, 39, 743–781. [Google Scholar] [CrossRef]
  7. Chokbunpiam, T.; Chanajaree, R.; Caro, J.; Janke, W.; Remsungnen, T.; Hannongbua, S.; Fritzsche, S. Separation of Nitrogen Dioxide from the Gas Mixture with Nitrogen by Use of ZIF Materials; Computer Simulation Studies. Comput. Mater. Sci. 2019, 168, 246–252. [Google Scholar] [CrossRef]
  8. Gal, A.; Kurahashi, M.; Kuzumoto, M. Effect of O3 on NO2 Sorption from Gas over H-Y Zeolite:  Supposition on the Nitrate Anion Formation with NO2 and O3 as Coreactants. J. Phys. Chem. A 2000, 104, 10821–10824. [Google Scholar] [CrossRef]
  9. Muckenhuber, H.; Grothe, H. A DRIFTS Study of the Heterogeneous Reaction of NO2 with Carbonaceous Materials at Elevated Temperature. Carbon N. Y. 2007, 45, 321–329. [Google Scholar] [CrossRef]
  10. Pietrzak, R.; Bandosz, T.J. Reactive Adsorption of NO2 at Dry Conditions on Sewage Sludge-Derived Materials. Environ. Sci. Technol. 2007, 41, 7516–7522. [Google Scholar] [CrossRef]
  11. Zhu, X.; Zhang, L.; Zhang, M.; Ma, C. Effect of N-Doping on NO2 Adsorption and Reduction over Activated Carbon: An Experimental and Computational Study. Fuel 2019, 258, 116109. [Google Scholar] [CrossRef]
  12. Iranimanesh, A.; Yousefi, M.; Mirzaei, M. DFT Approach on SiC Nanotube for NO2 Gas Pollutant Removal. Lab-in-Silico 2021, 2, 38–43. [Google Scholar] [CrossRef]
  13. Li, J.R.; Kuppler, R.J.; Zhou, H.C. Selective Gas Adsorption and Separation in Metal-Organic Frameworks. Chem. Soc. Rev. 2009, 38, 1477–1504. [Google Scholar] [CrossRef] [PubMed]
  14. Ebrahim, A.M.; Levasseur, B.; Bandosz, T.J. Interactions of NO2 with Zr-Based MOF: Effects of the Size of Organic Linkers on NO2 Adsorption at Ambient Conditions. Langmuir 2013, 29, 168–174. [Google Scholar] [CrossRef] [PubMed]
  15. Khan, N.A.; Hasan, Z.; Jhung, S.H. Adsorptive Removal of Hazardous Materials Using Metal-Organic Frameworks (MOFs): A Review. J. Hazard. Mater. 2013, 244–245, 444–456. [Google Scholar] [CrossRef]
  16. Wang, W.; Xu, X.; Zhou, W.; Shao, Z. Recent Progress in Metal-Organic Frameworks for Applications in Electrocatalytic and Photocatalytic Water Splitting. Adv. Sci. 2017, 4, 1600371. [Google Scholar] [CrossRef]
  17. Della Rocca, J.; Liu, D.; Lin, W. Nanoscale Metal-Organic Frameworks for Biomedical Imaging and Drug Delivery. Acc. Chem. Res. 2011, 44, 957–968. [Google Scholar] [CrossRef] [Green Version]
  18. Zhou, H.C.; Long, J.R.; Yaghi, O.M. Introduction to Metal-Organic Frameworks. Chem. Rev. 2012, 112, 673–674. [Google Scholar] [CrossRef]
  19. Zhou, H.C.J.; Kitagawa, S. Metal-Organic Frameworks (MOFs). Chem. Soc. Rev. 2014, 43, 5415–5418. [Google Scholar] [CrossRef] [Green Version]
  20. Bobbitt, N.S.; Mendonca, M.L.; Howarth, A.J.; Islamoglu, T.; Hupp, J.T.; Farha, O.K.; Snurr, R.Q. Metal-Organic Frameworks for the Removal of Toxic Industrial Chemicals and Chemical Warfare Agents. Chem. Soc. Rev. 2017, 46, 3357–3385. [Google Scholar] [CrossRef]
  21. Frysali, M.G.; Klontzas, E.; Tylianakis, E.; Froudakis, G.E. Tuning the Interaction Strength and the Adsorption of CO2 in Metal Organic Frameworks by Functionalization of the Organic Linkers. Microporous Mesoporous Mater. 2016, 227, 144–151. [Google Scholar] [CrossRef]
  22. Klontzas, E.; Mavrandonakis, A.; Tylianakis, E.; Froudakis, G.E. Improving Hydrogen Storage Capacity of MOF by Functionalization of the Organic Linker with Lithium Atoms. Nano Lett. 2008, 8, 1572–1576. [Google Scholar] [CrossRef] [PubMed]
  23. Fioretos, K.A.; Psofogiannakis, G.M.; Froudakis, G.E. Ab-Initio Study of the Adsorption and Separation of NO x and SO x Gases in Functionalized IRMOF Ligands. J. Phys. Chem. C 2011, 115, 24906–24914. [Google Scholar] [CrossRef]
  24. Stergiannakos, T.; Tylianakis, E.; Klontzas, E.; Froudakis, G.E. Enhancement of Hydrogen Adsorption in Metal-Organic Frameworks by Mg 2+ Functionalization: A Multiscale Computational Study. J. Phys. Chem. C 2010, 114, 16855–16858. [Google Scholar] [CrossRef]
  25. Frysali, M.G.; Klontzas, E.; Froudakis, G.E. Ab Initio Study of the Adsorption of CO(2) on Functionalized Benzenes. Chemphyschem 2014, 15, 905–911. [Google Scholar] [CrossRef]
  26. Giappa, R.M.; Tylianakis, E.; Di Gennaro, M.; Gkagkas, K.; Froudakis, G.E. A Combination of Multi-Scale Calculations with Machine Learning for Investigating Hydrogen Storage in Metal Organic Frameworks. Int. J. Hydrogen Energy 2021, 46, 27612–27621. [Google Scholar] [CrossRef]
  27. Giappa, R.M.; Papadopoulos, A.G.; Klontzas, E.; Tylianakis, E.; Froudakis, G.E. Linker Functionalization Strategy for Water Adsorption in Metal–Organic Frameworks. Molecules 2022, 27, 2614. [Google Scholar] [CrossRef]
  28. Skylaris, C.-K.; Gagliardi, L.; Handy, N.C.; Ioannou, A.G.; Spencer, S.; Willetts, A. On the Resolution of Identity Coulomb Energy Approximation in Density Functional Theory. J. Mol. Struct. THEOCHEM 2000, 501–502, 229–239. [Google Scholar] [CrossRef] [Green Version]
  29. Kozuch, S.; Gruzman, D.; Martin, J.M.L. DSD-BLYP: A General Purpose Double Hybrid Density Functional Including Spin Component Scaling and Dispersion Correction. J. Phys. Chem. C 2010, 114, 20801–20808. [Google Scholar] [CrossRef]
  30. Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  31. Weigend, F. Accurate Coulomb-Fitting Basis Sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef]
  32. Caldeweyher, E.; Bannwarth, C.; Grimme, S. Extension of the D3 Dispersion Coefficient Model. J. Chem. Phys. 2017, 147, 034112. [Google Scholar] [CrossRef] [PubMed]
  33. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  35. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef] [PubMed]
  36. Neese, F. The ORCA Program System. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 73–78. [Google Scholar] [CrossRef]
  37. Neese, F. Software Update: The ORCA Program System, Version 4.0. WIREs Comput. Mol. Sci. 2018, 8, e1327. [Google Scholar] [CrossRef]
  38. Boys, S.F.; Bernardi, F. The Calculation of Small Molecular Interactions by the Differences of Separate Total Energies. Some Procedures with Reduced Errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  39. Laaksonen, L. A Graphics Program for the Analysis and Display of Molecular Dynamics Trajectories. J. Mol. Graph. 1992, 10, 33–34. [Google Scholar] [CrossRef]
  40. Bergman, D.L.; Laaksonen, L.; Laaksonen, A. Visualization of Solvation Structures in Liquid Mixtures. J. Mol. Graph. Model. 1997, 15, 301–306. [Google Scholar] [CrossRef]
  41. Peng, D.-Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
  42. Mohamed, E.; Jaheon, K.; Nathaniel, R.; David, V.; Joseph, W.; Michael, O.; Yaghi, O.M. Systematic Design of Pore Size and Functionality in Isoreticular MOFs and Their Application in Methane Storage. Science 2002, 295, 469–472. [Google Scholar] [CrossRef] [Green Version]
  43. Jones, J.E.; Chapman, S. On the Determination of Molecular Fields.—II. From the Equation of State of a Gas. Proc. R. Soc. London. Ser. A Contain. Pap. A Math. Phys. Character 1924, 106, 463–477. [Google Scholar] [CrossRef] [Green Version]
  44. Breneman, C.M.; Wiberg, K.B. Determining Atom-centered Monopoles from Molecular Electrostatic Potentials. The Need for High Sampling Density in Formamide Conformational Analysis. J. Comput. Chem. 1990, 11, 361–373. [Google Scholar] [CrossRef]
  45. Yang, J.; Ren, Y.; Tian, A.; Sun, H. COMPASS Force Field for 14 Inorganic Molecules, He, Ne, Ar, Kr, Xe, H2, O2, N2, NO, CO, CO2, NO2, CS2, and SO2, in Liquid Phases. J. Phys. Chem. B 2000, 104, 4951–4957. [Google Scholar] [CrossRef]
  46. Matito-Martos, I.; Rahbari, A.; Martin-Calvo, A.; Dubbeldam, D.; Vlugt, T.J.H.; Calero, S. Adsorption Equilibrium of Nitrogen Dioxide in Porous Materials. Phys. Chem. Chem. Phys. 2018, 20, 4189–4199. [Google Scholar] [CrossRef] [PubMed]
  47. Rappe, A.K.; Casewit, C.J.; Colwell, K.S.; Goddard, W.A.; Skiff, W.M. UFF, a Full Periodic Table Force Field for Molecular Mechanics and Molecular Dynamics Simulations. J. Am. Chem. Soc. 1992, 114, 10024–10035. [Google Scholar] [CrossRef]
  48. Dubbeldam, D.; Calero, S.; Ellis, D.E.; Snurr, R.Q. RASPA: Molecular Simulation Software for Adsorption and Diffusion in Flexible Nanoporous Materials. Mol. Simul. 2016, 42, 81–101. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The functionalized linker of IRMOF-8 considered in the GCMC simulations; the original IRMOF-8 linker (a), the -OPO3H2 (b), -OSO3H (c), and -PO3H2 (d) functionalized linker. Carbon, hydrogen, oxygen, sulfur, and phosphorus atoms are depicted as gray, white, red, yellow, and orange spheres, respectively.
Figure 1. The functionalized linker of IRMOF-8 considered in the GCMC simulations; the original IRMOF-8 linker (a), the -OPO3H2 (b), -OSO3H (c), and -PO3H2 (d) functionalized linker. Carbon, hydrogen, oxygen, sulfur, and phosphorus atoms are depicted as gray, white, red, yellow, and orange spheres, respectively.
Molecules 27 03448 g001
Figure 2. RI-DSD-BLYP D3(BJ) /def2-TZVPP optimized geometries of benzene and functionalized molecules interacting with NO2. Hydrogen, carbon, oxygen, nitrogen, lithium, sulfur, and phosphorus atoms are depicted as white, gray, red, blue, purple, yellow, and orange spheres, respectively.
Figure 2. RI-DSD-BLYP D3(BJ) /def2-TZVPP optimized geometries of benzene and functionalized molecules interacting with NO2. Hydrogen, carbon, oxygen, nitrogen, lithium, sulfur, and phosphorus atoms are depicted as white, gray, red, blue, purple, yellow, and orange spheres, respectively.
Molecules 27 03448 g002
Figure 3. Electron density redistribution plots of the optimized geometries of the NO2…C6H5-X complexes. With red and blue the regions that gain and lose electron density upon the formation of the complex, respectively.
Figure 3. Electron density redistribution plots of the optimized geometries of the NO2…C6H5-X complexes. With red and blue the regions that gain and lose electron density upon the formation of the complex, respectively.
Molecules 27 03448 g003
Figure 4. Gravimetric (a) and volumetric (b) NO2 uptake isotherms for IRMOF-8 and IRMOF-8-X (X: -OSO3H, -PO3H2, -OPO3H2) at 298 K.
Figure 4. Gravimetric (a) and volumetric (b) NO2 uptake isotherms for IRMOF-8 and IRMOF-8-X (X: -OSO3H, -PO3H2, -OPO3H2) at 298 K.
Molecules 27 03448 g004
Figure 5. Snapshots for the IRMOF-8, IRMOF-8-X (X: -OSO3H, -PO3H2, -OPO3H2) from the GCMC at 298 K-0.1 bar and 298 K-1.0 bar.
Figure 5. Snapshots for the IRMOF-8, IRMOF-8-X (X: -OSO3H, -PO3H2, -OPO3H2) from the GCMC at 298 K-0.1 bar and 298 K-1.0 bar.
Molecules 27 03448 g005
Table 1. Binding energies in kcal/mol of the NO2…C6H5-X systems, calculated at the RI-DSD-BLYP D3(BJ)/def2-TZVPP level of theory and percentage of binding energy enhancement with respect to the introduction of the unfunctionalized benzene. All interaction energies are corrected for the BSSE by the full counterpoise method [39].
Table 1. Binding energies in kcal/mol of the NO2…C6H5-X systems, calculated at the RI-DSD-BLYP D3(BJ)/def2-TZVPP level of theory and percentage of binding energy enhancement with respect to the introduction of the unfunctionalized benzene. All interaction energies are corrected for the BSSE by the full counterpoise method [39].
SystemBinding Energy
(kcal/mol)
Binding Energy Enhancement (%)
NO2…C6H5-OSO3H−5.4170%
NO2…C6H5-OPO3H2−4.6131%
NO2…C6H5-PO3H2−4.2110%
NO2…C6H5-OCONH2−3.470%
NO2…C6H5-C(OH)3−3.470%
NO2…C6H5-SO3H−3.470%
NO2…C6H5-CONH2−3.050%
NO2…C6H5-SOOH−2.944%
NO2…C6H5-COOH−2.944%
NO2…C6H5-H−2.00%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Raptis, D.; Livas, C.; Stavroglou, G.; Giappa, R.M.; Tylianakis, E.; Stergiannakos, T.; Froudakis, G.E. Surface Modification Strategy for Enhanced NO2 Capture in Metal–Organic Frameworks. Molecules 2022, 27, 3448. https://doi.org/10.3390/molecules27113448

AMA Style

Raptis D, Livas C, Stavroglou G, Giappa RM, Tylianakis E, Stergiannakos T, Froudakis GE. Surface Modification Strategy for Enhanced NO2 Capture in Metal–Organic Frameworks. Molecules. 2022; 27(11):3448. https://doi.org/10.3390/molecules27113448

Chicago/Turabian Style

Raptis, Dionysios, Charalampos Livas, George Stavroglou, Rafaela Maria Giappa, Emmanuel Tylianakis, Taxiarchis Stergiannakos, and George E. Froudakis. 2022. "Surface Modification Strategy for Enhanced NO2 Capture in Metal–Organic Frameworks" Molecules 27, no. 11: 3448. https://doi.org/10.3390/molecules27113448

APA Style

Raptis, D., Livas, C., Stavroglou, G., Giappa, R. M., Tylianakis, E., Stergiannakos, T., & Froudakis, G. E. (2022). Surface Modification Strategy for Enhanced NO2 Capture in Metal–Organic Frameworks. Molecules, 27(11), 3448. https://doi.org/10.3390/molecules27113448

Article Metrics

Back to TopTop