Abstract
This study is focused on the adsorption of hexavalent chromium ions Cr(vi) using graphene oxide (GO). The GO was prepared by chemical oxidation (Hummers method) of graphite particles. The synthesized GO adsorbent was characterized by Fourier transform infrared spectroscopy and UV-Vis spectroscopy. It was used for the adsorption of Cr(vi) ions. The theoretical calculations based on density functional theory and Monte Carlo calculations were used to explore the preferable adsorption site, interaction type, and adsorption energy of GO toward the Cr(vi) ions. Moreover, the most stable adsorption sites were used to calculate and plot noncovalent interactions. The obtained results are important as they give molecular insights regarding the nature of the interaction between GO surface and the adsorbent Cr(vi) ions. The found adsorption energy of −143.80 kcal/mol is indicative of the high adsorptive tendency of this material. The adsorption capacity value of GO toward these ions is q = 240.361 mg/g.
1 Introduction
As a result of human activity (mining, chemical industry, manufacturing, agriculture, etc.), water resources are partially contaminated to some extent, and this problem continues to be present. Heavy metals also contribute to this variety of chemical pollutants, as effluent products or by-products of the industry [1,2,3]. These pollutants are of great and long-lasting concern for environmental and health safety due to their properties such as nonbiodegradability and bioaccumulation ability in living organisms [1,4]. A common heavy metal polluting species found frequently in the environment is hexavalent chromium ions Cr(vi) [5]. Generally, chromium exists in several stable oxidation states such as Cr(0), the trivalent Cr(iii), and the hexavalent Cr(vi) species. Cr(vi) ions in the forms of chromate, dichromate, and chromium trioxide are considered the most toxic forms of chromium. In this form, they possess an increased oxidizing potential, solubility, and mobility in the environment and living organisms. Cr(vi) is considered a carcinogen [6] and reflects toxicity toward many living organisms (plants, microorganisms, animals, etc.) [3,4,5,6]. The source of Cr(vi) in wastewaters is the effluents discharged from the industrial treatment [7]. The most cost-effective process for the removal of Cr(vi) from wastewater is the adsorption process [8].
Recently, the 2D carbon nanomaterials such as graphene oxide (GO) due to their extraordinary properties are also explored as adsorbents. GO is a product obtained from the chemical exfoliation and oxidation of layered graphite [9,10,11,12,13]. The adsorptive properties of GO are due to its huge surface area and the oxygenated functional groups on its structure [10,11]. This material is ideal as an adsorbent in aqueous media as it is highly hydrophilic, forming stable aqueous dispersions. Theoretical studies are valuable and often applied to study carbon nanomaterials to better understand their interaction with different types of molecules or ions [14,15].
To experimentally put in evidence the adsorptive properties of GO toward the dichromate ions, AAS measurements were performed. After evaluating experimentally the adsorptive ability of GO, theoretical calculations based on density functional theory (DFT) and Monte Carlo (MC) were used to obtain molecular-level insights for better understanding this process.
2 Methods
2.1 Synthesis of GO
The GO synthesis is based on the procedure published by Abdolhosseinzadeh et al. [16]. One gram of graphite flakes was added to 50 mL of cooled concentrated sulfuric acid (stirred in an ice-water bath). Then, 3 g of potassium permanganate was gradually added carefully by keeping the temperature below 10°C. The formed suspension was stirred at room temperature for 25 min followed by 5 min sonication in the ultrasonic bath. After repeating the stirring–sonication process for 12 times, the reaction content is quenched by the addition of 200 mL of DI water. This solution is left for 2 h under ultrasonication. A 20 mL 30% hydrogen solution was added to the graphite oxide suspension and stirred until no gas evolution is observed. The GO was washed with 1 M hydrochloric acid solution and deionized water several times for 30 min each by centrifugation. The obtained GO precipitates are dried at room temperature.
2.2 Adsorption studies
For the evaluation of adsorptive properties of GO, 100 mL of dichromate solution (50 mg L−1) was added into a 200 mL vessel with a fixed GO mass of adsorbent (0.015 g) and was agitated at a speed of 500 rpm at 25°C. The dichromate content of the solution is analyzed (after the adsorbent was separated from the samples by centrifugation and filtering) after 1–60 min by atomic absorption spectrometry (AAS).
The amount of adsorption q (mg kg−1) is calculated as follows:
where Cinit and Ceq are the initial concentration and equilibrium concentration of Cr(vi) ions, respectively (mg L−1), V is the volume of aqueous phase (mL), and m is the mass (g) of GO.
The percentage of adsorption (%) can be calculated using the following equation:
2.3 AAS measurements
AAS measurements were done using a High-Resolution Continuum Atomic Absorption Spectrometer with continuum light source – AAS Model novAA® 300 from Analytik Jena AG.
2.4 Computational details
Ab initio calculations were carried out with the DFT using the Dmol3 software. A previously employed GO model served for the adsorption study consisted of a 5 × 5 graphene layer (with epoxy, carboxy, and hydroxyl functional groups on its surface) [18]. The calculations were based on the triple-numeric quality with the polarization functions (TNP) basis set [19]. The exchange–correlation energy is described by the Perdew–Burke–Ernzerhof (PBE) functional within the generalized gradient approximation (GGA) [18,20]. The van der Waals interaction on the calculations was accounted for using the Tkatchenko–Scheffler method [21]. Self-consistent-field (SCF) convergence criterion was set to the root-mean-square change in the electronic density to be less than 10−6 eV. All energy minima were characterized by performing a vibrational analysis to ensure the lack of imaginary frequencies [22,23].
The adsorption energy [14,18,23] was calculated as follows:
where EGO/Cr(vi) is the total energy of the adsorption system. EGO and ECr(vi) are the energies of the isolated GO and Cr(vi) ions. The noncovalent interaction (NCI) was computed using Multiwfn software [24]. The NCI surface is plotted using Visual Molecular Dynamics software [25].
The electron density difference (EDD) was calculated using Cambridge Sequential Total Energy Package software [26]. The geometry optimization is performed by using the GGA-PBE [18,20]. A norm-conserving pseudopotential was used for the calculations. The Broyden–Fletcher–Goldfarb–Shanno algorithm was applied for the optimization and EDD of the model, with the SCF tolerance set at 2.0 × 10−6 eV/atom.
2.5 MC calculations
For MC, the interaction is between dichromate ions and the surface of GO in the simulated adsorption aqueous media. The calculation was performed under periodic boundary conditions. For the simulation, one potassium dichromate + 100 water molecules (as a solvent) were used (the system charges were assigned using the charge equilibration [QEq] method). The MC is performed by using several cycles of simulated annealing. During the annealing process, the temperature was set automatically at 10,000–100 K for each annealing cycle. The adsorption configurations were acquired as the temperature was slowly decreased. The Metropolis MC method was used to evaluate the adsorption configurations (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies – COMPASS II force field [27] [energy convergence tolerance of 2 × 10−5 kcal/mol and force convergence tolerance of 0.001 kcal/mol/A]) of the interaction between the molecules and the substrates.
Ethical approval: The conducted research is not related to either human or animal use.
3 Results and discussion
Figure 1 presents the Fourier transform infrared spectroscopy (FTIR) spectra of GO.
FTIR spectra of the synthesized GO.
The spectra of GO show the presence of bands related to C–O (νC–O at 1,050 cm−1), C–O–C (νC–O–C at 1,230 cm−1), C–OH (νC–O–H at 1,380 cm−1), and CO in carboxylic acid and carbonyl moieties that are present mostly along sheet edges but also on the basal plane of graphene sheets (νCO at 1,715 cm−1) and a broad peak between 3,000 and 3,550 cm−1 corresponding to O–H vibration [13,28,29]. The UV-Vis spectra of GO are presented in Figure 2.
UV-Vis absorption spectra of GO in DI water.
In water, the GO dispersion has a peak at around 232 nm, attributed to π → π* transitions of aromatic C═C bonds, and a peak at around 301 nm due to n → π* transition of aromatic C═O bonds [13,30]. As evaluated by AAS, the adsorption of the Cr(vi) onto the GO is a fast process. The adsorption percentage after a short contact time (t = 15 min) reaches 73.15% (q = 233.35 mg/g), and this value almost remains the same regardless of the time interval up to 24 h, where this value is 74.12%. The same holds also for the adsorption capacity. In this case, the adsorption capacity value is q = 240.361 mg/g. A value that is in the same range as the GO/polyamidoamine dendrimer composites for the Cr(vi) adsorption [17].
From Figure 3, the adsorption energies for the Cr(vi) ion onto the GO surface are in the range from −133.7 to −143.8 kcal/mol depending on the adsorption site. The highest adsorption energy is obtained for the Pose_3, in which the dichromate ion (where its oxygen atom bonded between two Cr atoms) is located near the surface hydroxyl group of the GO.
![Figure 3 (a) Four different explored adsorption sites for the dichromate ion on the GO surface (the GO model is approximated by using a 5 × 5 graphene model that contains epoxy, hydroxyl, and carboxy groups on its surface) (level of theory: GGA/PBE/DFT-D = TS [Tkatchenko and Scheffler]//TNP – triple numerical polarized basis set; basis file: 4.4) and (b) the corresponding adsorption energy values.](/document/doi/10.1515/chem-2020-0148/asset/graphic/j_chem-2020-0148_fig_003.jpg)
(a) Four different explored adsorption sites for the dichromate ion on the GO surface (the GO model is approximated by using a 5 × 5 graphene model that contains epoxy, hydroxyl, and carboxy groups on its surface) (level of theory: GGA/PBE/DFT-D = TS [Tkatchenko and Scheffler]//TNP – triple numerical polarized basis set; basis file: 4.4) and (b) the corresponding adsorption energy values.
To correctly assess the interaction type [14,22,23] between the chromate ions and GO surface, the NCI plot and the reduced density gradient (RDG) vs. sign(λ)ρ were computed [31] for the highest adsorption site, in this case, Pose_3. The visualization of the NCI surface (Figure 4(a)) and the RDG vs sign(λ)ρ (Figure 4(b)) supports that the adsorption interaction is a van der Waals one.
(a) Noncovalent interaction surfaces and (b) the plot of RDG vs sign(λ)ρ for the interaction of GO/dichromate ion (Pose_3).
MC calculations are frequently used to understand the adsorptive properties of different molecular entities at the surface of materials [32,33], although this calculation is semi-quantitative. The adsorption geometry and energy distributions for the GO/Cr(vi)/100 H2O system obtained by the adsorption locator module are presented in Figure 5.
(a) Side and top views of the equilibrium adsorption configurations of the dichromate ions (one dichromate ion + 100 water molecules) onto the GO model (5 × 5 graphene layer containing: an epoxy, carboxy, and hydroxyl functional groups on its surface obtained using MC simulations) and (b) adsorption energy distributions for dichromate ions.
As shown in Figure 5, the adsorption energy of Cr(vi) is 290.32 kcal mol−1, which indicates the presence of strong adsorptive interaction between the GO surface and the Cr(vi). The adsorption energy values are comparable with those found in the previous studies (where a different and a smaller GO model is used) [34]. It is important to observe that as computed by DFT calculation, the preferred adsorption site of the Cr(vi) is also predicted by MC. This is associated with the fact that electrostatic interactions can take place at this adsorption site as found previously by Mondal et al. [34].
The EDD map of the adsorbed dichromate ion on the GO surface (Pose_3) is presented in Figure 6. The areas presented in yellow depict the increase (gain) of the electron density and those in blue depict the decrease (lost) among the interacting entities: GO and the dichromate ion. The electrons from the GO surface are transferred to the dichromate ion during the adsorption process. Electron increase and a decrease among oxygen atoms of dichromate ion and the hydroxyl group support the electrostatic type of interaction among GO and dichromate ion.
EDD for the adsorbed dichromate ion on the GO surface (Pose_3).
4 Conclusions
The adsorption of the Cr(vi) ions onto GO is explored both experimentally using AAS measurements and theoretically using DFT and MC calculations. The experimental results show that the GO has fast adsorption and high adsorption capacity toward the Cr(vi) ions. The adsorption capacity is 240.3 mg/g, making GO a possible candidate for the application in the removal of this toxic chemical species. Furthermore, the high adsorption energy values as evidenced by theoretical calculation show that these ions strongly interact with the groups in the graphene surface, particularly with the hydroxyl ones.
The computed NCI surfaces point out that in the adsorption of these ions, an important role is played by van der Waals interaction. The EDD shows that the electrostatic interaction takes place where the dichromate ions attract the electron density from the GO surface. The present study not only puts in evidence the practical applicability of GO as an effective adsorbent for Cr(vi) ions but also gives molecular insights regarding the adsorption sites of GO, interaction type, and adsorption energetics of this material toward the toxic hexavalent chromium ions.
Acknowledgments
The authors gratefully acknowledge the support from the Ministry of Education, Science and Technology of Kosovo (No. 2-5069) for providing him with the computing resources. The Turkish Chemistry Society support is highly appreciated.
Conflict of interest: The authors declare no conflict of interest.
References
[1] Nagajyoti PC, Lee KD, Sreekanth TVM. Heavy metals, occurrence and toxicity for plants: a review. Environ Chem Lett. 2010;8:199–216.10.1007/s10311-010-0297-8Search in Google Scholar
[2] McGrath SP, Zhao FJ, Lombi E. Plant and rhizosphere processes involved in phytoremediation of metal-contaminated soils. Plant Soil. 2001;232:207–14.10.1007/978-94-010-0566-1_20Search in Google Scholar
[3] Ali H, Khan E, Ilahi I. Environmental chemistry and ecotoxicology of hazardous heavy metals: Environmental persistence, toxicity, and bioaccumulation. J Chem. 2019;2019:14.10.1155/2019/6730305Search in Google Scholar
[4] Tchounwou PB, Yedjou CG, Patlolla AK, Sutton DJ. Heavy metal toxicity and the environment. EXS. 2012;101:133–64.10.1007/978-3-7643-8340-4_6Search in Google Scholar PubMed PubMed Central
[5] Oliveira H. Chromium as an environmental pollutant: insights on induced plant toxicity. J Bot. 2012;2012:1–8.10.1155/2012/375843Search in Google Scholar
[6] Park RM, Bena JF, Stayner LT, Smith RJ, Gibb HJ, Lees PSJ. Hexavalent chromium and lung cancer in the chromate industry: a quantitative risk assessment. Risk Anal. 2004;24:1099–108.10.1111/j.0272-4332.2004.00512.xSearch in Google Scholar PubMed
[7] Aydin YA, Aksoy ND. Adsorption of chromium on chitosan: optimization, kinetics and thermodynamics. Chem Eng J. 2009;151:188–94.10.1016/j.cej.2009.02.010Search in Google Scholar
[8] Fernández-López JA, Angosto JM, Avilés MD. Biosorption of hexavalent chromium from aqueous medium with opuntia biomass. Sci World J. 2014;2014:670249.10.1155/2014/670249Search in Google Scholar PubMed PubMed Central
[9] Park S, Ruoff RS. Chemical methods for the production of graphenes. Nat Nanotechnol. 2009;4:217–24.10.1038/nnano.2009.58Search in Google Scholar PubMed
[10] Chen D, Feng H, Li J. Graphene oxide: preparation, functionalization, and electrochemical applications. Chem Rev. 2012;112:6027–53.10.1021/cr300115gSearch in Google Scholar PubMed
[11] Dideikin AT, Vul’ AY. Graphene oxide and derivatives: the place in graphene family. Front Phys. 2019;6:149.10.3389/fphy.2018.00149Search in Google Scholar
[12] Marcano DC, Kosynkin DV, Berlin JM, Sinitskii A, Sun Z, Slesarev A, et al. Improved synthesis of graphene oxide. ACS Nano. 2010;4:4806–14.10.1021/nn1006368Search in Google Scholar PubMed
[13] Dreyer DR, Park S, Bielawski CW, Ruoff RS. The chemistry of graphene oxide. Chem Soc Rev. 2010;39:228–40.10.1039/B917103GSearch in Google Scholar
[14] Gao H, Feng W, Li X, Li N, Du Y, Wu Y, et al. Insights into the non-covalent interaction between modified nucleobases and graphene nanoflake from first-principles. Phys E. 2019;107:73–79.10.1016/j.physe.2018.11.002Search in Google Scholar
[15] Bai H, Gao H, Feng W, Zhao Y, Wu Y. Interaction in Li@fullerenes and Li+@fullerenes: first principle insights to Li-based endohedral fullerenes. Nanomaterials. 2019;9:630.10.3390/nano9040630Search in Google Scholar PubMed PubMed Central
[16] Abdolhosseinzadeh S, Asgharzadeh H, Kim HS. Fast and fully-scalable synthesis of reduced graphene oxide. Sci Rep. 2015;5:1–7.10.1038/srep10160Search in Google Scholar PubMed PubMed Central
[17] Liu H, Zhang F, Peng Z. Adsorption mechanism of Cr(vi) onto GO/PAMAMs composites. Sci Rep. 2019;9:1–12.10.1038/s41598-019-40344-9Search in Google Scholar PubMed PubMed Central
[18] Berisha A. Interactions between the aryldiazonium cations and graphene oxide: a DFT study. J Chem. 2019;2019:5.10.1155/2019/5126071Search in Google Scholar
[19] Delley B. Ground-state enthalpies: evaluation of electronic structure approaches with emphasis on the density functional method. J Phys Chem A. 2006;110:13632–9.10.1021/jp0653611Search in Google Scholar PubMed
[20] Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys Rev Lett. 1996;77:3865–8.10.1103/PhysRevLett.77.3865Search in Google Scholar PubMed
[21] Tkatchenko A, Scheffler M. Accurate molecular van der Waals interactions from ground-state electron density and free-atom reference data. Phys Rev Lett. 2009;102:073005.10.1103/PhysRevLett.102.073005Search in Google Scholar PubMed
[22] Berisha A. The influence of the grafted aryl groups on the solvation properties of the graphyne and graphdiyne – a MD study. Open Chem. 2019;17:703–10.10.1515/chem-2019-0083Search in Google Scholar
[23] Berisha A, Combellas C, Kanoufi F, Médard J, Decorse P, Mangeney C, et al. Alkyl-modified gold surfaces: characterization of the Au–C bond. Langmuir. 2018;34:11264–71.10.1021/acs.langmuir.8b01584Search in Google Scholar PubMed
[24] Lu T, Chen F. Multiwfn: a multifunctional wavefunction analyzer. J Comput Chem. 2012;33:580–92.10.1002/jcc.22885Search in Google Scholar PubMed
[25] Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14:33–38.10.1016/0263-7855(96)00018-5Search in Google Scholar PubMed
[26] Clark SJ, Segall MD, Pickard CJ, Hasnip PJ, Probert MJ, Refson K, et al. First principles methods using CASTEP. Z Kristallogr. 2005;220:567–70.10.1524/zkri.220.5.567.65075Search in Google Scholar
[27] Sun H, Jin Z, Yang C, Akkermans RLC, Robertson SH, Spenley NA, et al. COMPASS II: extended coverage for polymer and drug-like molecule databases. J Mol Model. 2016;22:1–10.10.1007/s00894-016-2909-0Search in Google Scholar PubMed
[28] Ossonon BD, Bélanger D. Synthesis and characterization of sulfophenyl-functionalized reduced graphene oxide sheets. RSC Adv. 2017;7:27224–34.10.1039/C6RA28311JSearch in Google Scholar
[29] Lomeda JR, Doyle CD, Kosynkin DV, Hwang WF, Tour JM. Diazonium functionalization of surfactant-wrapped chemically converted graphene sheets. J Am Chem Soc. 2008;130:16201–6.10.1021/ja806499wSearch in Google Scholar PubMed
[30] Rommozzi E, Zannotti M, Giovannetti R, D’Amato C, Ferraro S, Minicucci M, et al. Reduced graphene oxide/TiO2 nanocomposite: from synthesis to characterization for efficient visible light photocatalytic applications. Catalysts. 2018;8:598.10.3390/catal8120598Search in Google Scholar
[31] Řezáč J, Hobza P. Benchmark calculations of interaction energies in noncovalent complexes and their applications. Chem Rev. 2016;116:5038–71.10.1021/acs.chemrev.5b00526Search in Google Scholar PubMed
[32] Dagdag O, Berisha A, Safi Z, Hamed O, Jodeh S, Verma C, et al. DGEBA-polyaminoamide as effective anti-corrosive material for 15CDV6 steel in NaCl medium: Computational and experimental studies. J Appl Polym Sci. 2020;137:48402.10.1002/app.48402Search in Google Scholar
[33] Berisha A. Experimental, Monte Carlo and molecular dynamic study on corrosion inhibition of mild steel by pyridine derivatives in aqueous perchloric acid. Electrochem. 2020;1:188–99.10.3390/electrochem1020013Search in Google Scholar
[34] Mondal MK, Roy D, Chowdhury P. Designed functionalization of reduced graphene oxide for sorption of Cr(vi) over a wide pH range: a theoretical and experimental perspective. New J Chem. 2018;42:16960–71.10.1039/C8NJ03794ASearch in Google Scholar
© 2020 Ardhmeri Alija et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.
