Nguyen Van Hao1, Nguyen Van Dang1, Do Hoang Tung2, Pham The Tan3, Nguyen Van Tu4, Pham Van Trinh4,5. 1. Faculty of Physics and Technology, TNU - Thai Nguyen University of Science Thai Nguyen City Vietnam dangnv@tnus.edu.vn. 2. Institute of Physics, Vietnam Academy of Science and Technology 18 Hoang Quoc Viet Hanoi Vietnam. 3. Hung Yen University of Technical Education Khoai Chau Hung Yen Province Vietnam. 4. Institute of Materials Science, Vietnam Academy of Science and Technology 18 Hoang Quoc Viet Hanoi Vietnam trinhpv@ims.vast.vn. 5. Graduate University of Science and Technology, Vietnam Academy of Science and Technology 18 Hoang Quoc Viet Hanoi Vietnam.
Toxic metal pollution is a serious environmental problem today. Toxic metal ions are often difficult to biodegrade and easily accumulate in both plants and animals, thus causing danger to humans and organisms.[1-3] Therefore, they are necessary to be removed from wastewater before it goes into the environment. There are many methods to remove toxic metals including chemical reduction, co-precipitation, filtration, ion exchange, ion flotation, adsorption, etc. which have been used to treat toxic metal ion-containing wastewater.[4-9] In particular, the adsorption method is considered a more attractive method due to its simple operation, low cost, easy to recycle adsorbents and high efficiency.[10-15] Many previous reports have shown that there are many materials capable of adsorbing toxic metals such as activated carbon, fly ash, sawdust, rice husks, zeolite, iron oxide and manganese.[16-21] However, these adsorbents have a low effective removal of toxic metals due to their low affinity.[22] Therefore, new materials need to be searched for more effective removal of toxic metals.In recent years, graphene-based materials have attracted a lot of attention due to their unique properties.[23] In which, graphene oxide (GO) is a type of graphene containing oxygen-rich functional groups such as –OH, –CO, COOH… on the surface thereby increasing their reactivity and hydrophilic properties. Therefore, GO could be used in many different fields such as materials technology, energy, health and environmental treatment.[24] Up to now, many techniques have been developed for synthesizing GO including mechanical exfoliation, the liquid-phase exfoliation. The preparation of GO from inexpensive graphite powder is a promising method for low-cost and large-scale manufacturing. However, it is well-known that graphite is currently used as an electrode material for most of commercial batteries and thus could become a critical material in the near future. The reuse of the graphite of the recycled batteries could be a potential solution for the future demand of graphite.Thus, the goal of this work is to prepare GO from the graphite rod of recycled battery by solution plasma exfoliation at atmospheric pressure. The prepared GO having a large surface area and high porosity has been applied for removing Pb(ii) ions from the water. The effects of several factors such as pH, contact time, temperature and initial concentration of Pb(ii) on the adsorption of GO were also analyzed in detail. Besides, the adsorption kinetics, isothermal and thermodynamic modeling studies were performed to examine the adsorption properties of GO for Pb(ii).
Experimental
Materials
A graphite rod is taken from the recycled battery core (AA SIZE, Panasonics – Malaysia). KOH, (NH4)2SO4, NaOH, HCl chemicals were purchased from the Xinyang chemical (China). Pb(NO3)2 in analytical adsorption is supplied by Sigma Aldrich. The experimental solution is double distilled water.
GO synthesis experiments
The solution plasma exfoliated processing of the synthesis of GO using the graphite rod of recycled battery is shown in Fig. 1. The system consists of a 150 ml beaker used as a reaction vessel, two electrodes are placed parallel and deeply immersed in the electrolyte solution. The positive electrode is a graphite rod collected from old AA batteries with a diameter of 5.0 mm and a length of 50 mm. The negative electrode is a pure platinum rod (99.9%). The electrolyte solution is a mixture of 500 ml of 0.5 mM KOH and 0.2 mM (NH4)2SO4. The two electrodes are connected by a home-made DC high voltage supply. It is a modulated DC power supplier with an average voltage of 2.45 kV and a current of 22.4 mA. Plasma is created between the two electrodes at a distance of 0.2 mm (Fig. 1a). After the plasma process, the obtained solution is filtered and washed with PVDF film (0.2 μm) to remove impurities and organic compounds from the battery core. Finally, the samples were dried at 60 °C in a vacuum for 14 hours.
Fig. 1
The solution plasma exfoliated processing of graphene oxide synthesis using graphite rod of recycled battery core (a) experimental setup and (b) possible mechanism.
Characterization
Samples are characterized by field emission scanning electron microscopy (FESEM, Hitachi S4800), atomic force microscopy (AFM, PARK XN100), X-ray diffractometer (Miniflex Rigaku), absorption spectra (UV-Vis, Jasco-V770), Raman spectroscopy by using 532 nm wavelength at room temperature (XploRA, Horiba Jobin-Yvon), and spectrum Fourier transform infrared (FTIR) with 4 cm−1 resolutions (Spectrum Two, PerkinElmer, USA). The BET specific surface areas (SBET) of GO were determined by N2 adsorption/desorption isotherms at 77 K (BET, Builder, SSA-4300).
Batch adsorption experiments
The adsorption efficiency of GO for Pb(ii) was evaluated and investigated via parameters including contact time (5–180 minutes), pH (3–9), contact temperature (293–313 K) and initial concentration (60–200 mg L−1). A Pb(ii) solution of 1000 mg L−1 was prepared using Pb(NO3)2 dispersed in distilled water. The initial pH of Pb(ii) solution is adjusted from 3 to 9 using NaOH and HCl solutions. The pH of the solution is measured by a pH meter (PHS-3C, Hinotek, China).The Pb(ii) adsorption experiments were performed by adding 50 mg of GO to a 100 ml conical glass flask containing 50 ml of Pb(ii) solution and shaking with a shaker (model: MaxQ 4000 Benchtop, Thermo Scientific, USA) over time (at a speed of 120 rpm). The solution is then filtered, centrifuged to separate the supernatant. The Pb(ii) concentration at various contact time was analyzed using AAS (Hitachi, Japan). Adsorption efficiency R (%) and adsorption capacity qe are determined by the following equation:where C0 is the initial concentration of Pb(ii) (mg L−1), Ce is the equilibrium concentration of Pb(ii) after adsorption (mg L−1), qe is the equilibrium adsorption capacity of Pb(ii) (mg g−1), m is the mass of GO adsorbent (g), and V is the volume of Pb(ii) solution (L).
Results and discussions
The characteristics of GO
Fig. 2a and b shows the morphology of GO produced by plasma technology in solution from the recycled battery core. The average thickness and diameter of the prepared GO were determined by the statically counting of about more than 100 sheets from AFM images as shown in Fig. 2c and d. The results indicated that the average thickness of nanosheet size of 3.3 nm corresponds to the number of layers from 8–10 layers, the average diameter of about 1.5 μm. Fig. 2e presents the N2 adsorption/desorption isotherm curves at 77 K of graphite rod and GO. The results show BET surface areas and the total pore volume of graphite and GO are 15.78 m2 g−1, 0.074 cm3 g−1 and 128.37 m2 g−1, 0.36 cm3 g−1, respectively. The characteristics of porous structure and BET surface area of graphite and GO are detailed in Table 1, it demonstrated that the surface area of GO significantly increases compared to graphite.
Fig. 2
(a) SEM and (b) AFM images of GO, the statistical distribution of thickness and diameter (c and d) and (e) N2 adsorption/desorption isotherms exfoliation of graphite and GO at 77 K.
Surface characteristics of graphite and GO
Materials
SBET (m2 g−1)
Average pore diameter (nm)
Total pore volume (cm3 g−1)
Graphite
15.78
9.56
0.074
GO
128.37
19.73
0.36
Fig. 3a shows the UV-Vis spectra of aqueous graphite, GO and GO dispersion after Pb(ii) adsorption process. The results showed that the UV-Vis spectra of GO have two observed features, the first is a shoulder at ∼310 nm, attributed to the plasmon peak related to the n–π* transition due to the presence of the oxygen-containing functional groups and the second feature is a peak appearing near 235 nm are related to the plasmon peak due to the transition of the π–π* electron in C–C the aromatic ring bond of the graphene layers.[25-28] While the UV-Vis spectra of graphite show only one wide band with a slight transition peak of ∼275 nm as a characteristic of this material.[28-33] It is interesting noted that the absorption spectrum of GO after Pb(ii) adsorption has a significant change both in curve form and in spectral intensity, the first is the plasmon peak at 235 nm is shifted to 231 nm, while the shoulder at about 310 nm is shifted to the position 305 nm with a steeper slope. This result shows that Pb(ii) ions were adsorbed on the surface of GO. As a result, the optical absorption of GO is dominated by the π–π* plasmon peak near 230 nm, this plasmon peak depends on the effects related to the nanometer-scale sp2 clusters and the CC, CO and C–O bonds, respectively, with deconvoluted peaks at the binding energy ranges of 284.8–285.0, 288.0–288.9 and 286.0–286.9 eV.[34-40]
Fig. 3
(a) UV-Vis spectra (b) XRD patterns, (c) Raman spectra and (d) FTIR spectra of graphite rod, GO and GO after adsorption of Pb(ii).
Fig. 3b shows the XRD patterns of graphite and GO. The results show that in addition to the peak with strong intensity at the angle of 2θ = 26.24° (002), some peaks appearing at 43.24° (100) and 54.92° (004) which are typical peaks of graphite. As for GO, besides typical peaks of graphite, a new peak with high intensity appears at 2θ = 10.71°. It is noted that the typical peak of graphite (002) is slightly shifted to 26.34° with a lower intensity. This is attributed to the graphite do not completely oxidize during the plasma process and still retain part of its graphite structure. The d-distance of GO was determined to be 0.8745 nm, implying that some oxygen functional groups were attached to the surface of graphene layers. When comparing the XRD patterns of GO before and after Pb(ii) adsorption, the intensity of the peak at 2θ = 10.71° is decreased and the peak position is slightly shifted. In addition, the presence of new peaks at 2θ = 21.46°, 26.11°, 29.53°, 30.70°, 35.50°, 42.33°, 43.65°, 54.42° indicating that Pb(ii) was adsorbed on the surface of GO.[41,42]Fig. 3c shows the Raman spectrum of graphite, GO and GO/Pb. Raman spectra of graphite showed the presence of D-band at 1350 cm−1 with very low intensity representing defects related to vacancy and grain boundaries and G-band (graphite) at 1579 cm−1 with high intensity relative to sp2 carbon atoms due to scattering of the E2g phonons.[43,44] The 2D-band at 2717 cm−1 with relatively high intensity is attributed to a two-phonon double resonance.[45,46] For GO, the Raman spectra all showed a slight shift in the bands (D-band at 1350.2 cm−1, G-band at 1580 cm−1 and 2D-band at 2691 cm−1). The ID/IG ratios of graphite and GO were 0.04 and 0.58, respectively. The increase in the ID/IG ratio can be explained by oxidation of graphite during the plasma exfoliation process in solution and thereby also confirming the degree of surface functionalization of GO (e.g., introduces a few of containing oxygen functional groups).[46] In addition, the 2D band of GO has been slightly shifted to low frequencies with a significantly higher intensity than graphite, which indicated the formation of the graphene sheets from the graphite rod during the plasma generation process.[47,48] After Pb(ii) adsorption, the D band and G band are shifted to 1363 and 1572 cm−1, respectively and the ID/IG ratio increased up to 0.96. The peak shifts and the increase of the ID/IG ratio are attributed to the adsorption of Pb(ii) on the surface of GO.[42]Fig. 3d shows the FTIR spectrum of graphite, GO and GO after adsorption of Pb(ii). For graphite, the FTIR spectrum is a wide range covering the entire spectrum without any peak containing the functional groups such as CO, C–O, COO–.[49] While the FTIR spectrum of GO is prepared by plasma process in a solution from graphite rods taken from the recycled battery core showing signs of oxygen-derived species, the functional groups are obvious in the structure, such as a strong band at 3310 cm−1 attributable to the stretch oscillation of the OH group, a narrower peak at 1637 cm−1 related to the stretching oscillations from carbonyl and carboxylic groups (CO), the peak that appears at 1399 cm−1 corresponds to the stretching oscillations of COO− (ref. 50) and the peak at 1058 cm−1 indicates the stretching oscillations associated with the C–O group on the surface of the GO.[51] This indicates that functional groups containing oxygen already exist on the surface of the GO plates. Compared with the FTIR spectrum of GO before adsorption, the peak intensity at 3310 cm−1 and 1637 cm−1 has a slight decrease and shift after Pb(ii) adsorption. This is explained by the interaction of Pb(ii) ions with oxygen-containing functional groups on the surface of GO. The peak at 1058 cm−1 also has decreased intensity and the red shift. This is probably due to the electron delocalization in the surface structure of GO when Pb(ii) is adsorbed, causing the density to decrease and the bond length increasing, leading to a red shift.[42,52] In addition, a new peak at 1215 cm−1 appeared in the spectrum, attributed to the stretching vibration of CC.
Removal Pb(ii) in water using GO
The adsorption is a process that depends on many factors such as the pH of the solution, temperature, the initial concentration of the toxic metal and the properties of the adsorbent.PH is a parameter that greatly influences the adsorption capacity of the adsorbent. To investigate the effect of PH on the removal of Pb(ii) ions, the pH of the solution was changed in the range from 3 to 9 at the fixed initial concentration of Pb(ii) of 100 mg L−1. Fig. 4a shows the effect of pH on the adsorption capacity of GO for Pb(ii). The results showed that adsorption capacity increased when the pH increased from 3 to 6 and reached a stable value at pH about 6–7 and then decreased slightly. The increase in adsorption capacity can be explained by the competition between H+ and Pb2+ ions decreased on GO or due to the reduction of electrostatic repulsion between GO and Pb2+ surfaces because of positive charge on the surface decreases,[53,54] therefore it is stable at pH values of 6–7. The reduced adsorption at pH > 7 may be due to the precipitation of Pb(ii) on the GO surface.[55,56] As a result, a pH value of 6 was chosen for the subsequent experiments.
Fig. 4
Effect of (a) pH, (b) contact time, (c) initial concentration and (d) temperature on the Pb(ii) adsorption of GO.
Fig. 4b shows the effect of contact time on the adsorption for Pb(ii) of GO. Experimental conditions at pH 6, initial concentration of Pb(ii) 100 mg l−1, the GO concentration is 1 mg ml−1 in 50 ml of solution at room temperature and the contact time is changed from 5 to 180 minutes. The obtained results showed that, in the first 20 minutes, the adsorption capacity increased rapidly. As increasing the adsorption time, the adsorption increased to a stable value of about 40 minutes and this is considered to be the adsorption equilibrium time. Fig. 4c shows the effect of the initial concentration of Pb(ii) on the adsorption of GO. The results indicated that, when the initial concentration of Pb(ii) increased from 60 mg l−1 to 200 mg l−1, the adsorption capacity increased accordingly, from 59.5 mg g−1 to 113.7 mg g−1. This result can be explained by the fact that under optimal pH conditions, the surface of GO will contain H+ ions that help elevate ions interactions with the functional groups on the GO surface or because GO has highly porous, thus it could adsorb more Pb(ii) ions at higher concentrations.[56,57]The effect of temperature on Pb(ii) adsorption by GO was performed by changing the temperature in the range of 293–313 K, during 60 minutes contact time and pH 6 and as illustrated in Fig. 4d. The results showed that, when the temperature increased from 293 to 313 K, the adsorption capacity increased slightly from 96.84 to 97.24 mg g−1. This indicated that the Pb(ii) adsorption process is an endothermic process because it is possible that the porosity and total pore volume of the adsorbent increase, leading to increases in the number of available active sites on the adsorbent as the temperature increases, however, it does not increase much. The same tendency has been obtained in some previous studies.[56-62]To better understand the effect of contact time onto the Pb(ii) adsorption process, the first-order second Pseudo (eqn (3)) and second-order Pseudo models (eqn (4)) were used to fit data:where qe and q are equilibrium adsorption and at time t for Pb(ii) (mg g−1) of GO, respectively and k1, k2 are the rate constants of first-order Pseudo kinetic adsorption and second-order Pseudo adsorption kinetic models, respectively.Fig. 5a shows the linear fit curve of the second-order Pseudo adsorption kinetic model data. The calculated results show that the adsorption kinetics of Pb(ii) on GO fits well with both kinetic models. Table 2 shows the parameters obtained from fit data following the first-order Pseudo and second-order Pseudo models. As a result, the value of R2 for the second-order Pseudo kinetic model (0.9983) is higher than the first-order Pseudo kinetic model (0.9871), which shows the action time in the adsorption process is mainly based on chemical adsorption. These results are consistent with the previous publication.[42,57,63]
Fig. 5
(a) Linear fit curve of the second-order Pseudo adsorption kinetic model data, (b) adsorption isotherms of Pb(ii) onto the GO (1 mg ml−1) at different initial concentrations (60–200 mg l−1) after 60 min and (c) linear plot of ln Kdvs. 1/T for Pb(ii) adsorption onto GO at T = 293, 303, 313 K.
Kinetic parameters for Pb(ii) adsorption models
Kinetic models
Pseudo-first-order model
Pseudo-second-order model
Parameters
R2
qe (mg g−1)
k1 (min−1)
R2
qe (mg g−1)
k2 (g mg−1 min−1)
0.9871
94.98
0.1197
0.9983
102.98
0.0018
To study precisely the equilibrium relationship between GO and Pb(ii), the two most common adsorption isotherm models, the Langmuir and Freundlich isotherm models were utilized to describe the Pb(ii) adsorption onto GO. The Langmuir and Freundlich isotherm models are expressed as.[64,65]where qe (mg g−1) is the equilibrium adsorption capacity of Pb(ii) onto GO, qm (mg g−1) is the maximum adsorption capacity of Pb(ii) onto GO, Ce (mg L−1) is the equilibrium concentration of Pb(ii), kL is the Langmuir constant, kF and 1/n are Freundlich constants. Fig. 5b presents the adsorption isotherms of Pb(ii) on GO at different initial concentrations. The isothermal parameters of those models were given in Table 3. The fit result of the model shows that the value of R2 of the Langmuir model is higher than that of the Freundlich model, which indicates the suitability in Langmuir isotherms of GO for Pb(ii). The adsorption capacity of GO calculated as 180.1 mg g−1 is quite high compared to the list of other adsorbents (Table 4).
Isothermal parameters of Pb(ii) adsorption
Isotherm models
Langmuir
Freundlich
Parameters
R2
qm (mg g−1)
kL (L mg−1)
R2
1/n
kF (g mg−1)
0.9913
180.13
0.008
0.9692
0.4908
8.673
q
max of Pb(ii) by different adsorbents
Adsorbents
qmax (mg g−1)
Year of publication
Ref.
GO
81.3
2019
42
Nickel ferrite–rGO
121.6
2018
66
SV–rGO
95
2019
67
CNTs
38.76
2011
70
GO
111
2016
56
GO-SH1
142.8
2016
56
GO-SH2
200
2016
56
GOA
158.7
2016
68
GO–Zr–P
384.6
2016
65
GO/chitosan
99
2011
69
GO
180.1
—
This work
Thermodynamic parameters including Gibbs free energy ΔG°, entropy change ΔS°, enthalpy change ΔH° are investigated through adsorption of Pb(ii) onto GO at different temperatures and is described as by the following equations:[71]where T is the temperature (K), R is the gas constant (8.314 J mol−1 K−1), Kd is the adsorption equilibrium coefficient, ΔS° and ΔH° are determined by the intercept and slope of the linear fitting of ln(Kd) to 1/T plot, respectively.The thermodynamic parameters determined from Fig. 5c are presented in Table 5. Positive ΔH° values indicate the adsorption of Pb(ii) onto GO as an endothermic process. This result is in good agreement with the previous discussion. Moreover, the free energy change Gibbs ΔG° < 0 and when the environment temperature increases, the value of ΔG° decreases, which shows the absorption of Pb(ii) on the GO surface is more favorable when the temperature increases, this is because at higher temperatures the ions are more prone to dehydration. With the value ΔH° = 5.9 kJ mol−1 < 40 kJ mol−1, so this is a physisorption process.[63,71] At the same time, the value of ΔS° > 0 also confirmed that the adsorption was spontaneous with high affinity.
Thermodynamic parameters of Pb(ii) adsorption onto GO
T (K)
ΔG° (kJ mol−1)
ΔS° (J mol−1 K−1)
ΔH° (kJ mol−1)
293
−8.58
49.45
5.90
303
−9.08
313
−9.57
To investigate the recyclability of the GO material, the regeneration experiment was repeated 7 times. Fig. 6 shows the adsorption efficiency of GO on Pb(ii) when using regenerated GO for 7 repeated uses. The regeneration of the GO adsorbent was evaluated for up to 7 consecutive cycles and after each cycle, the adsorption efficiency was reduced. However, after 6 cycles, the Pb(ii) adsorption efficiency still reached 72.94%. This can be explained by a variation in the surface activity position of GO. This is because of the incomplete separation of the Pb(ii) ions involved in ion-exchange interactions. Besides, another possible reason for performance degradation per cycle could be partial erosion of the GO surface in a strongly acidic environment. This obtained result is consistent with previous studies on Pb(ii) adsorption[42,72-74] and shows the recyclability and stability of the prepared materials.
Fig. 6
Regeneration experiment for Pb(ii) adsorption by GO.
The possible mechanisms of adsorption
Based on the analysis results of UV-Vis absorption spectra, FTIR, Raman spectra and XRD of GO before and after Pb(ii) adsorption as well as the discussion of the results above, Pb(ii) removal mechanism by GO can be proposed as in Fig. 7. From the optical absorption analysis of Pb(ii)-loaded GO (Fig. 3a), it can be seen that the change of plasmon absorption peak is near 230 nm compared to the peak at 235 nm of GO, this indicates that there is Pb(ii) attachment to the surface of GO. This is due to the complexation interaction of Pb(ii) ions and the π–π* bonding electrons involving bonds such as CC (at energy range 284.8–285.0 eV), CO (in the range energy 288.0–288.9 eV), O–CO (at the energy range 286.0–289.6 eV), –COOH (at the energy range 530–531.5 eV) and –CO (at the energy range 531.8–532.2 eV).[75-78]Fig. 3d illustrates the FTIR spectra of the loaded Pb(ii) GO after adsorption showing peaks with a decrease in intensity and transition from 3310 cm−1 to 3406 cm−1, from 1637 cm−1 to 1627 cm−1, which could due to the complexation of Pb(ii) with the –OH group and the –CO group in carboxylic acid, respectively.[79,80] The peak shift to 1075 cm−1 could be attributable to hydroxyl group re-formation on lead oxide or to the C–O stretching vibration of –COOH due to the hydroxyl group complexation of Pb(ii) with sites.[79,81] The appearance of a new peak at 1215 cm−1, is attributed to the stretching vibration of CC group in the nanometer-scale sp2 carbon skeletal lattice. Previous studies of Pb(ii) binding with oxygen-containing functional groups on the surface of oxidized poly-walled CNTs, GO and rGO have been reported[42,79,82,83] and the results from spectral analysis UV-Vis, FTIR, the surface complexing reactions of GO can be represented by the following equations:[83]
Fig. 7
Possible mechanism for Pb(ii) adsorption onto GO.
The schematic representation of the proposed Pb(ii) adsorption at positions on the surface of the GO plate as shown in Fig. 7. The Pb(ii) is mainly attached to the positions of the hydroxyl and carbonyl groups (as given in eqn (11)–(15)). This can be considered as an approximate representation of the adsorption capacity of Pb(ii) on the surface of GO.[79]
Conclusions
We have investigated the method basing on a plasma exfoliated process to prepare GO from a recycled battery core in solution was performed. This is considered to be a simple, inexpensive, novel and environmentally friendly technique. The prepared GO has an average thickness of 3.5 nm and a length of 1.5 μm. The properties of GO are suitable for the adsorption of toxic metals. Experimental results have demonstrated that the adsorption capacity of GO for Pb(ii) can reach a maximum value of 180.1 mg g−1 at room temperature. Studies of isotherms and adsorption kinetic models have shown that adsorption of GO onto Pb(ii) fits well with the Langmuir model. Calculations of thermodynamic parameters (ΔG°, ΔH°, và ΔS°) show that the adsorption of Pb(ii) ions on the GO surface is spontaneous and intrinsically heat-absorbed. The main Pb(ii) adsorption mechanism of GO is the formation of complexation of Pb(ii) in the π–π* bond electrons and the electrostatic interaction between Pb(ii) ions and oxygen-containing functional groups are also dominant mechanisms for GO. All these properties make the solution plasma exfoliated GO as a new carbon material and attract a lot of attention, a promising adsorbent for many environmental treatment applications with a high recyclability and efficiency removal contaminants in aqueous solution.
Conflicts of interest
The authors declare no possible conflict of interests.
Authors: S J Olusegun; G L S Rodrigues; E T F Freitas; L R S Lara; W R Rocha; N D S Mohallem Journal: J Hazard Mater Date: 2019-07-10 Impact factor: 10.588