A Theoretical Study on the Inclusion of Fe, Cu, and Zn in Illite Clays

Research Article A Theoretical Study on the Inclusion of Fe, Cu, and Zn in Illite Clays Antonio Sánchez-Coronilla , 1 Elisa I. Martín , 2 Francisco José Fernández-de-Cordova, 3 Francisco Javier Santos, 1 and José Hidalgo Toledo 1 1 Departamento de Química Física, Facultad de Farmacia, Universidad de Sevilla, E-41012 Sevilla, Spain 2 Departamento de Ingeniería Química, Facultad de Química, Universidad de Sevilla, E-41012 Sevilla, Spain 3 Instituto de Investigaciones Químicas (IIQ), CSIC-Universidad de Sevilla, Avenida Américo Vespucio 49, E-41092 Sevilla, Spain Correspondence should be addressed to Antonio Sánchez-Coronilla; [email protected] and Elisa I. Martín; [email protected] Received 14 December 2018; Accepted 28 March 2019; Published 7 May 2019 Academic Editor: Ruibing Wang Copyright © 2019 Antonio Sánchez-Coronilla et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The inclusion of Fe, Cu, and Zn in (1, 0, 0), (0, 0.5, 0), and (0.5, 0.5, 0.5) sites of an illite with the KAl 2 Si 4 O 12 H 2 structure has been studied. For the inclusion of the metals, their common oxidation states were chosen, that is, 0, +2, +3 and 0, +1, +2, for Fe and Cu, respectively, while 0 and +2 for Zn. Periodic DFT calculations were performed to know the most favourable site of incorporation of the ions. Energetically the most favourable site for the inclusion corresponds to the (1, 0, 0) coordinate for all the ions independently of their oxidation state. However, the highest oxidation state of the metals (Fe 3+ ,Cu 2+ , and Zn 2+ ) was the most favoured for being incorporated into the illite structure and was the selected ion for the discussion. In those structures, metal oxygen interaction plays an important role in stabilizing the systems. Structural and energetic results indicate that illite presents good adsorption characteristic of those Fe 3+ ,Cu 2+ , and Zn 2+ in the (1, 0, 0) site. Thus, those ions may be available for plants for its extraction by phytoextraction techniques and the consequent soil regeneration. The inclusion of a second metallic ion revealed the most favourable inclusion corresponding to the inclusion of Fe 3+ ion. The inclusion of this ion modifies the coordination sphere around the first metal being available for subsequent extraction by phytoremediation or other techniques for clean-up of the soil and its regeneration. 1. Introduction Clay minerals are major components in soils and sedimentary rocks [1]. They play an important role in environmental processes such as nutrient cycling, plant growth, organic matter maturation, contaminant migration, and petroleum production and could coexist with methanogens in anoxic environments [2, 3]. Of relevance is the role played in methane adsorption and radionuclide adsorption on clay rocks [4, 5]. The latter is related to their use as sealing barriers in nuclear waste and spent nuclear fuel repositories [5]. Regarding the protective effect of mineral clays, illite is one of the clays that best counteract the disruption of the intestinal barrier permeability caused by mycotoxins [6]. Illite is a common soil component [1]. Illite is a natural and available mineral in the world with a variety of physicochemical characteristics such as porosity, small particle size, good bonding ability, great electrical insulating property, and thermal stability [7]. Illite is found as abundant phase in natural clays as a consequence of the transformation of smectite with general formula K y Al 4 (Si 8-y ,Al y )O 20 (OH) 4 usually with 1<y<15, but always with y<2where potassium is commonly localized in interlayer sites [8]. Silicon is replaced by aluminum isomorphically in the tetrahedral layer to maintain global charge [9]. It is of note that soils where illite is the predominant phase are characterized by high adsorption and low desorption capacities [10]. This clay has been used as a low-cost adsorbent for industrial waste treatment and as adsorbent for fast and efficient removal of phosphate from aqueous solution [10]. The phosphate adsorption is pH-dependent, and illite shows a great potential as adsorbent [11]. Hindawi Journal of Nanomaterials Volume 2019, Article ID 4546350, 14 pages https://doi.org/10.1155/2019/4546350 Soil contamination by heavy metals (As, Cd, Cr, Cu, Pb, and Zn) is one of the major environmental problems raising critical concerns for both human health and ecosystems due to their carcinogenic and mutagenic effects on animals and humans. The clean-up of most of these soils is mandatory for the area to be reclaimed and to minimize the entry of potentially toxic elements into the food chain [12–15]. Thus, polluted soils by metals are a source of negative environmental impact. These metals can pollute groundwater and aquifers with consequent environmental disaster. Therefore, both the extraction of heavy metals from polluted soils and a consequent regeneration of the contaminated soils have attracted the interest of researchers for the last years [16–18]. In fact, a greenhouse technique for recovering K Al Si O H Cu (a) K Al Si O H Cu (b) K Al Si O H Cu (c) K Al Si O H Cu (d) Figure 1: The local geometry for illite, KAl 2 Si 4 O 12 H 2 (a). Coordinates for the inclusion of Cu into the illite structure: (b) (0, 0.5, 0), (c) (0.5, 0.5, 0.5), and (d) (1, 0, 0). 2 Journal of Nanomaterials polluted soils based on the use of plant is phytoremedation [12, 18]. Phytoremediation is both an efficient and an economically feasible method that can be used for the removal of heavy metals and radionuclides as well as for organic pollutants such as polynuclear aromatic hydrocarbons, polychlorinated biphenyls, and pesticides [12, 18, 19]. Terrestrial plants have the capacity of adsorbing heavy metals present in the soil [14, 20]. This capacity depends on chemical factors such as the nature of the different metals and their availability within the structure of the soils. Therefore, the sorption/desorption of the metals depends on their nature, that is, the oxidation state of those metals [13, 14, 17, 21]. Thus, the earlier review of literature reveals the key role that plays the sorption and fixation of the metals into the soils for being extracted and the next soil regeneration. In this sense, it is of interest to know both how the metals are distributed into the clay and their oxidation state. In this study, an illite clay with composition KAl 2- Si 4 O 12 H 2 has been selected. Periodic calculations based on the density functional theory (DFT-periodic) were performed to know how the presence of some metals affects the illite structure. The chosen metals were Fe, Cu, and Zn due to their usual presence in soils near mining wastes. We have studied the inclusions of these metals individually and their combination in pairs. The individual atoms have been examined with different oxidation states and placed in the initial structure in different positions. Both the energetically more stable oxidation state and the metal position inside illite were selected. The positions of the individual atoms served as the initial point for the combination of the two metals. Finally, the electronic effect of the metals into the illite structure was explored through the electron localization function (ELF) and density of states (DOS) analysis. 2. Materials and Methods 2.1. Computational Framework. DFT-periodic calculations were performed using the Vienna Ab Initio Simulation Package (VASP) with the projector-augmented wave (PAW) method[22,23].ThenumberofplanewavesinVASPwascontrolled by a cut-offenergy, set in our calculations to Ecut = 500 eV to satisfactorily describe the system. The electron exchange and correlation were treated within the generalized gradient approximation (GGA). In the case of GGA, Perdew-Burke-Ernzerhof (PBE) functional was used. Both the cell shape and atomic positions were optimized using a conjugate-gradient algorithm, where the iterative relaxation of atomic positions was stopped when the forces on the atoms were less than 0.01eV/Å. Gaussian smearing with kBT=01 eV was applied [22–25]. Illite (KAl 2 Si 4 O 12 H 2 ) was characterized by the C2/m space group [26]. The tag KSPACING determined the number of k-points in the mesh. A value of 0.3 for sampling the Brillouin zone for the illite bulk was enough to obtain negligible changes in the optimized cell parameters and energy. The resulting cell parameters were a=5312 Å,b=8868 Å, and c=10401 Å with α=γ=90 ° ;β= 100 55 ° , in agreement with the experimental values (a=5202 Å,b=8979 Å,c= 10 226 Å;α=γ=90 ° ;β= 101 57 ° ) [26]. The metals Fe, Cu, and Zn were chosen, and the effect of its inclusion in the illite structure was studied from three positions that correspond to the coordinates (0, 0.5, 0), (0.5, 0.5, 0.5), and (1, 0, 0) in the unit cell. The density of states (DOS) and projected density of states (PDOS) for the relaxed structures were obtained using the tetrahedron method with Blöchl corrections and a KSPACING of 0.2. Electron localization function (ELF) images and structure images were obtained using the Vaspview software and Vesta, respectively [27, 28]. 3. Results and Discussion 3.1. Inclusion of Fe, Cu, and Zn in the Illite Structure. The inclusion of Fe, Cu, and Zn in the illite structure has been performed in three different coordinates, which are (0, 0.5, 0), (0.5, 0.5, 0.5), and (1, 0, 0) in the illite cell. Figure 1 shows the local geometry for the most stable configuration of the optimized illite structure. Also in this figure, the local geometry inclusion of Cu is shown as an example of the positions of the metals within the cell. The common oxidation states of Fe, Cu, and Zn were selected for their inclusion into the illite structure. Thus, for the inclusion of Fe and Cu, three oxidation states were chosen: 0, +2, +3 and 0, +1, +2, respectively, with only two states for Zn, that is, 0 and +2. To gain insight into the energetics of the effect of the inclusion of the metals in the illite structure, the difference in energy as consequence of the inclusion of the metals in the three coordinates was calculated as follows: ΔE=Emetal/illite −Eillite , 1 where Emetal/illite and Eillite are the total energies of the illite structure with the metal and free, respectively. The differences in energy for the different metals are included in Table 1. According to the energetic results shown in Table 1, the most favoured site for the adsorption of the metals corresponds to the coordinate (1, 0, 0). In the case of Zn 0 , similar results were obtained for the inclusion in coordinates (0.5, 0.5, 0.5) and (1, 0, 0), with the inclusion in the latter coordinate being slightly more stable. Of note is that in all positions, Table 1: Difference in energy for the inclusion of the metals in different sites of the illite structure. ΔE(eV) (0, 0.5, 0) (0.5, 0.5, 0.5) (1, 0, 0) Fe 3+ -28.35 -28.84 -31.59 Fe 2+ -22.62 -24.46 -28.08 Fe 0 -11.80 -15.97 -16.02 Cu 2+ -17.13 -20.71 -22.40 Cu + -16.79 -20.05 -21.13 Cu 0 -7.03 -10.97 -12.54 Zn 2+ -19.91 -20.39 -23.84 Zn 0 -4.49 -10.87 -10.90 3Journal of Nanomaterials the highest oxidation state of the metals (Fe 3+ ,Cu 2+ , and Zn 2+ ) is the most favoured for being incorporated into the illite structure. This result is of interest because it implies a major availability of those elements in the illite structure for being extracted, which agrees with experimental results of desorption of those metals [5, 14, 19, 21, 29]. In this way, from now on the discussion of the results will be mainly focused on the inclusion of the highest oxidation states of the metals. Figure 2 shows the local geometry in the coordination shell around Fe 3+ ,Cu 2+ , and Zn 2+ for the most stable configuration of the optimized structures in the coordinate (1, 0, 0) of the illite structure. For the purpose of comparison, the figure also includes the illite structure. As an example, the HO-M-OH angles and M-OH distances (M=Fe 3+ ,Cu 2+ , and Zn 2+ ) of the cut corresponding with the middle face of the structures as indicated in the dotted circle in Figure 2 are shown in Table 2. As Table 2 shows, αand βangles hardly change and are around 180 ° . In addition, it is also observed that, generally, the M-OH distances are slightly shorter for the inclusion of Fe 3+ than for Cu 2+ and Zn 2+ . Moreover, the (001) basal distance for the systems after the incorporation of Fe 3+ ,Cu 2+ , and Zn 2+ ions is 10.10053, 10.15051, and 10.15604Å, respectively. Such results suggest a strengthening of the interactions between the metals and the oxygen, which may be thought to be exclusively due to the smaller ionic radii of Fe 3+ than Cu 2+ and Zn 2+ , but it is also necessary to consider electronic effects produced in the structures, as will be discussed below from the ELF and PDOS. Thus, the ELF plot in 2D corresponding to the M-OH (M=Fe 3+ ,Cu 2+ , and Zn 2+ ) interactions for the inclusion in the (1, 0, 0) site of illite is shown in Figure 2. In this figure, the contour plot around the oxygen atoms is different for the systems with Fe 3+ ,Cu 2+ , and Zn 2+ when it is compared with those of illite. There is a modification of the halos of the contour plot Fe Zn Al Si O H Cu 01 Al Si O H (a) Fe Zn Al Si O H Cu A l Si O H a b c  d Central plane view 01 (b) Fe Zn Al Si O H Cu Al Si O H 01 (c) Fe Zn Al Si O H Cu Al S i O H 01 (d) Figure 2: The local geometry around the (1, 0, 0) inclusion site with the ELF contour plots from the central plane view of (a) illite, (b) Fe 3+ , (c) Cu 2+ , and (d) Zn 2+ . Table 2: Geometrical parameters, identified in Figure 2, of the structures studied. Structure Distance (Å) Angle ( ° ) abcd αβ Fe 3+ 2.27301 2.27290 3.43240 3.43234 179.9952 179.9964 Cu 2+ 2.27344 2.27328 3.43260 3.43154 179.9071 179.9410 Zn 2+ 2.27353 2.27332 3.43314 3.43300 179.9970 179.9966 4 Journal of Nanomaterials E (eV) Illite −2 0 20 40 60 120 110 100 90 80 70 60 50 40 30 20 10 0 −202468 80 100 120 02468 DOS (1/eV) O p-states PDOS (1/eV) E (eV) (a) Fe3+ −2 −100 −80 −60 −40 −20 0 20 40 60 80 100 0246 DOS (1/eV) E (eV) (b) Cu2+ −120 −100 −80 −60 −40 −20 0 20 40 60 80 10 120 −20246 E (eV) DOS (1/eV) (c) Figure 3: Continued. 5Journal of Nanomaterials towards Fe 3+ and Cu 2+ by means of slight electron localization (sky blue colour) which provides stability to the metaloxygen interaction. This highlights an attraction between the oxygen and the metals which explains the geometrical results expounded above (Table 2). The contour plot ELF for the system with Zn 2+ resembles that of illite. However, there is also a slight directionality in the outlines of the electron localization of oxygen towards Zn 2+ . ELF analysis indicates the strengthening of the interaction between the oxygen atoms from illite and metals. The inclusion of the metals in other coordinates causes distortion of the structure as compared to the parent structure. As an example, the structural distortion is clearly reflected in the variation in the Fe-OH distances (ca. 1.66 and 1.67Å, for aand bin Figure 2, respectively) and HO-Fe-OH angles (ca. 172.33 ° ) of the central plane view for the inclusion of Fe 3+ in the (0, 0.5, 0) site. That distortion of the structure is related to structural tensions and thereby with instability. Thus, it is understandable that the inclusion of the metals in the (1, 0, 0) site leads to the most stable structures with the minimal geometrical differences when compared to the illite structure. To further understand the electronic properties that underlie the metal-oxygen interactions, the electronic structure of the inclusion of Fe 3+ ,Cu 2+ , and Zn 2+ in the (1, 0, 0) site will be discussed by means of DOS and PDOS (Figure 3). The illite band gap is around 4.2eV (Figure 3(a)). The analysis of the projected density of states (PDOS) on the oxygen atoms of illite indicates that the valence and conduction bands are mostly composed of the O p-states (inset of Figure 3(a)). With the inclusion of the metals appear new bands in the gap for the inclusion of Fe 3+ and Cu 2+ , while for the inclusion of Zn 2+ only a new band near the valence band maximum appears (Figure 3). The presence of bands in the gaps for the inclusion of Fe 3+ led to the existence of new lowerenergy electronic transition with the lowest band gap value around 2.3eV, followed by the systems with Zn 2+ and Cu 2+ , ca. 3.5 and 3.9eV, respectively (Figure 3). This energy electronic transition should lead to changes in the experimental UV-Vis absorption spectra as compared to that of illite. In this sense, experimental band gap values can be obtained both by using UV-Vis spectra in reflectance diffuse mode, determined as reported elsewhere [30], and finally by comparing with those obtained in the calculations. The distribution and the chemical interaction of the metals within the structure are of interest because a rich transfer of the metals during water movement along the soil pore may exist, which has been previously shown by optical measurements in surface and subsurface soil horizons [31]. The analysis of PDOS for the inclusion of Fe 3+ and Cu 2+ reveals that those bands in the gap belong mostly to the contribution of Fe 3+ and Cu 2+ d-states that overlap with the bands from the O p-states (Figures 4(a) and 4(b)). This overlapping is indicative of the oxygen metal interaction which agrees with the above discussion from ELF analysis. A different situation is obtained from the analysis of the PDOS for the inclusion of Zn 2+ (Figure 4(c)). The presence of Zn 2+ in the illite structure E (eV) Zn2+ 100 80 60 40 20 0 −20246 DOS (1/eV) (d) Figure 3: Total DOS plot for (a) illite and the inclusion of (b) Fe 3+ , (c) Cu 2+ , and (d) Zn 2+ in the (1, 0, 0) site. In the inset: PDOS of the O p-states from illite. 6 Journal of Nanomaterials O p-statesFe d-states E (eV) E (eV) PDOS (1/eV) PDOS (1/eV) −100 −20 2 4 6 −80 −60 −40 −20 0 0 20 40 60 80 100 −60 −40 −20 0 0 20 40 60 −20 246 60 50 40 30 20 10 0 0 −20 −40 −60 (a) O p-statesCu d-states E (eV) E (eV) PDOS (1/eV) PDOS (1/eV) −60 −40 −20 0 0 20 40 60 −20246 −60 −40 −20 0 0 20 40 60 −2 0246 (b) Figure 4: Continued. 7Journal of Nanomaterials indicates the main contribution of the Zn p-states to the edge of the valence band that overlap with the bands from the O p-states of the structure (Figure 4(c)), revealing the interaction between the Zn and oxygen atoms in the structure, as stated above. Finally, based on this structural and energetic discussion, it seems reasonable to consider that illite presents a good adsorption character of those metals (Fe 3+ ,Cu 2+ , and Zn 2+ ) in the (1, 0, 0) site. Therefore, for instance, those metals may be available for plants for its extraction by using phytoextraction techniques and the consequent soil regeneration. 3.2. Inclusion of a Second Fe, Cu, and Zn Atom in the Illite Structure. As far as metallic compounds are hard to be found alone, it is of interest to explore the effect that the presence of two metallic atoms has on the illite structure. From the results expounded above, the structures with Fe 3+ ,Cu 2+ , and Zn 2+ in the (1, 0, 0) site for exploring theoretically the changes in the structural and electronic properties of the illite when a second metallic (Fe, Cu, or Zn) atom is included were selected. The oxidation state selected for the inclusion of the second element was Fe 3+ ,Cu 2+ , and Zn 2+ based on the energetic analysis shown above. In order to evaluate the structural changes ought to the presence of two metals, site (0, 0.5, 0) was selected for the inclusion of the second metal. For the study, we selected the (0, 0.5, 0) site instead of the (0.5, 0.5, 0.5) site because the potassium atoms in the latter will be placed around the metals and they may hinder the electronic interactions between the two metals. The local geometry for the most stable configuration of the optimized structures for the incorporation of the second Fe 3+ ,Cu 2+ , and Zn 2+ ions is shown in Figure 5. In Table 3 is included the difference in energy when a second metal is included in the illite structure. This difference is defined as ΔE=Emetal/metal/illite −Emetal/illite , 2 where Emetal/metal/illite and Emetal/illite are the total energies of the illite structure with the two metals and with only one, respectively. According to the results shown in Table 3, the most stable interactions correspond to those systems with the inclusion of the second Fe 3+ ion (structures 1,5, and 8 in Figure 5). The inclusion of the second atom modifies slightly the illite structure. The most striking changes are produced around the metals, so attention will be drawn to the Al atoms around the metals, as indicated in the central plane view of Figure 6. The geometrical changes for structures 1, 5, and 8are shown in Table 4. The inclusion of the second Fe 3+ decreases the Al-Fe distances (aand bin Figure 6) in the following order 1>8>5, with angle αfrom Figure 6 next to 180 ° for structure 5, with this structure being the most favoured energetically, in agreement with the energetic data from Table 3. The shortening of distances aand bin Figure 6 is accompanied by an increase in the values in the distances between Al-M (M=Fe, Cu, and Zn), namely, c and din Figure 6. Although the second Fe 3+ ion is surrounded by Al cations, Fe 3+ is stabilized by the two hydroxyl groups placed around this ion, as the distances around 1.7-1.8Å for all systems are attested (e and f in Table 4), and this will be corroborated below by the PDOS analysis. The ELF analysis reveals the existence of electron localization (green and sky blue colour) between the Al and Fe atoms in structures 1,5, and 8(dotted square in Figure 6). This result implies a certain degree of stabilizing interactions in the system by a redistribution of electronic localization focusing from Al towards Fe which enables its interaction with the two hydroxyl groups from structures 1,5, and 8. Thus, the inclusion of the second metal performs a kind of synergistic effect in the illite structure with regards the accommodation of the other metals in the structure for being available for a subsequent extraction by phytoremediation or other techniques. Finally, the possible electron changes originated by the presence of a second metal in the illite structure were analysed by means of DOS (Figure 7). The incorporation of Fe 3+ modifies the gap of illite with the appearance of new bands near the edge of the valence band for the three structures (1,5, and 8). Those bands are mainly from the Fe and O p-statesZn p-states E (eV) E (eV) PDOS (1/eV) PDOS (1/eV) −2 1.0 0.8 0.6 0.4 0.2 0.0 0246 −2 100 80 60 40 20 0.0 0246 (c) Figure 4: PDOS plot for the inclusion of (a) Fe 3+ (b) Cu 2+ , and (c) Zn 2+ in the (1, 0, 0) site. 8 Journal of Nanomaterials O d- and p-states, respectively (Figure 8). The overlaps of those bands corroborate with the Fe-O interaction for the inclusion of the second metal and the stabilization of the system. 4. Conclusions This paper presents the theoretical study of the stability and electronic effect of the inclusion of Fe, Cu, and Zn in different sites of an illite with the KAl 2 Si 4 O 12 H 2 structure. Three sites were selected in the illite structure, namely, (1, 0, 0), (0.5, 0.5, 0.5), and (0, 0.5, 0). For the inclusion of the metals, their common oxidation states were chosen: 0, +2, +3 and 0, +1, +2, for Fe and Cu, respectively, while for Zn, only two states 132 4 65 7 98 Fe ZnAl Si O H KCu Figure 5: The local geometry for the second inclusion of Fe 3+ ,Cu 2+ , and Zn 2+ ions in the (0, 0.5, 0) site for the structures with Fe 3+ (from 1 to 3), Cu 2+ (from 4 to 6), and Zn 2+ (from 7 to 9) in the (1, 0, 0) site. Table 3: Difference in energy when a second metal is included in the illite structure as shown in Figure 5. Structure ΔE(eV) 1-10.78 2-3.46 3-2.20 4-6.80 5-13.49 6-1.73 7-3.71 8-12.23 9-1.72 9Journal of Nanomaterials