Skip to main content

A density functional theory study of Fe(II)/Fe(III) distribution in single layer green rust: a cluster approach


Green rust (GR) is a potentially important compound for the reduction of heavy metal and organic pollutants in subsurface environment because of its high Fe(II) content, but many details of the actual reaction mechanism are lacking. The reductive capacity distribution within GR is a key to understand how and where the redox reaction occurs and computational chemistry can provide more details about the electronic properties of green rust. We constructed three sizes of cluster models of single layer GR (i.e., without interlayer molecules or ions) and calculated the charge distribution of these structures using density functional theory. We found that the Fe(II) and Fe(III) are distributed unevenly in the single layer GR. Within a certain range of Fe(II)/Fe(III) ratios, the outer iron atoms behave more like Fe(III) and the inner iron atoms behave more like Fe(II). These findings indicate that the interior of GR is more reductive than the outer parts and will provide new information to understand the GR redox interactions.


Green rusts (GR) are a family of Fe(II), Fe(III) layered double hydroxides (LDH) that frequently form in oxygen-poor, Fe(II)-rich soils and waters [1,2,3,4,5]. GRs are composed of positively charged Fe(II), Fe(III) hydroxide layers that alternate with hydrated interlayers containing anions (e.g. SO42−, CO32−, Cl, etc.) and occasionally cations for charge compensation. GRs are represented by the general formula [FeII(1-x)FeIIIx(OH)2]x+[(x/n)An,mH2O]x, where x represents the molar fraction of the ferric ion that usually ranges from 0.25 to 0.33 and An− denotes intercalated anions [2, 6,7,8]. GRs are classified into two types based on the anion they intercalate: [1] GR type 1 has a rhombohedral unit cell and intercalates planar or monatomic anions (e.g., Cl, CO32−) [9, 10]. GR type 2 has a hexagonal unit cell and intercalates three-dimensional anions (e.g., SO42−) [11].

GRs have been widely investigated for removal of organic and inorganic contaminants from waters and soils [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27] (e.g., SeO42− [13], U6+ [15], TcO4 [17], Ag+, Au3+, Cu2+ and Hg2+ [18], CCl4 [16, 22, 23], NO3 [14], CrO42− [24,25,26,27]) due to their excellent reducing capacity. In these studies, researchers have proposed several mechanisms to explain redox reactions by GRs. [28] For example, Hansen et al. [14] suggested that nitrate (NO3) reduction by chloride GR is faster compared to other GRs because chloride interlayer exchange by nitrate, which is then reduced in the interlayer. Similarly, several studies argued that the reduction of chromate (CrO42−) occurs in the GR interlayer following anion exchange with chromate [25, 29,30,31]. In contrast, Thomas et al. [26] proposed that chromate is directly reduced at sulphate GR particle surface sites by electrons shuttled from the particle’s interior via electron hopping [32]. Also, a similar surface process was propsed by Choi et al. [21] and O’Loughlin et al. [15] for the reduction of perchloroethene (C2Cl4) and U(VI) by GR, respectively. Indeed, C2Cl4 and U(VI) are unlikely to intercalate into the GR interlayer. Also, the TEM images in O’Loughlin et al. [15] show UO2 nanoparticles decorating GR particle edge surfaces after reaction further supporting reduction reaction at GR particle edges [26, 32]. Lastly, recent studies have shown that GR cannot reduce chlorinate ethenes [33], however, if a catalyst such as bone char is added to reaction, chlorinated ethenes can be rapidly reduced by GR, which is contrast to previous study that shows pure GRs cannot reduce chlorinated ethenes [33,34,35], suggesting that electron transfer can occur in GRs.

These examples demonstrate that despite many experimental studies on GR reductive capacity with contaminants, there is still much confusion about where on GR particles reduction reactions occur, and whether GR hydroxide sheets can shuttle electrons from their interior to the exterior. To unravel the active redox sites on GR hydroxide sheets, we must gain a better understanding of GR electronic properties at a molecular level and evaluate the distribution of Fe(II) and Fe(III) atoms within single GR layers. In this study, we constructed 3 different sized cluster models of a single GR layer and then analyzed its electronic properties using density functional theory (DFT). Specifically, we evaluated the relationship between spin state and total electronic energy and used an implicit solvent model to take into account the interlayer water between the GR hydroxide sheets. The obtained results give new insights into the distribution of Fe(II) and Fe(III) species in GR hydroxide sheets, allowing to make suggestions of GR active redox sites and their role in reactions with contaminants.

Calculation methods

Calculation methods

All geometry optimizations were performed with the Turbomole program, v6.5 [36]. The COSMO implicit solvent model [37] with an infinite dielectric constant was combined with Becke–Perdew (BP) functional [38, 39] and triple-ζvalence plus polarization (TZVP) [40] basis set for all calculations. To solve the convergence difficulties of the calculation that result from the high degree of freedom in the models, a higher orbital shifting parameter (0.3) was applied to all calculations.

GR structure model

The GR structure was based on crystallographic data for sulphate GR (GRSO4) provided in Christiansen et al. [3] with following formula: NaFe(II)6Fe(III)3(SO4)2(OH)18.12H2O. GRSO4 particles form hexagonal platelets, which consist of hydroxide layers where all octahedral sites are occupied, and the interlayer spaces are filled with octahedrally hydrated sodium and sulphate ions, along with additional water. Here, we only focus on the hydroxide sheet structure (i.e., Fe(II), Fe(III) and OH ions), for which we constructed three different sized hexagonal clusters representing single GR hydroxide layers, shown in Fig. 1. We modified the edges to ensure that every Fe atom was coordinated to O atoms of six hydroxyl groups. The small (GR2 × 6), medium (GR3 × 6) and large (GR4 × 6) cluster have 2, 3 and 4 Fe atoms, respectively, located on each of the 6 edges (Fig. 1a).

Fig. 1

a Schematic representation of single layer GR clusters with different sizes: (S) small (GR2 × 6), (M) medium (GR3 × 6) and (L) large (GR4 × 6). The brown dots represent Fe atoms. b Top and c side view of single layer GR clusters, where Fe atoms are colored brown, O atoms red and H atoms white

We define the ratio of Fe(II) to Fe(III) atoms in the GR cluster by applying a specific charge to the clusters. This in turn defines the number of unpaired electrons (Fe atoms are assumed to be in high spin states, which means Fe(II) has 4 unpaired electrons and Fe(III) has 5 unpaired electrons; discussed in Sect. 2.3). The tested cluster charges, Fe(II)/Fe(III) ratios and numbers of unpaired electrons are shown for each cluster in Table 1.

Table 1 Range of applied Fe(II)/Fe(III) ratios, charges and numbers of unpaired electron (NUE) for GR cluster model

Multiplicities assumption

To verify our assumption of Fe electronic occupation, the energy differences between a set of possible spin states for the built GR cluster models were calculated and compared. Specifically, we calculated all possible spin states for GR2 × 6(− 5) (i.e., GR2 × 6 with applied cluster charge of − 5), while for GR2 × 6(− 6), GR3 × 6(− 10), and GR4 × 6(− 9), we only calculated a low, intermediate and high spin state as shown in Fig. 2. Across all performed calculations, the high spin states yield the lowest energy compared to the intermediate and low spin states. Looking at GR2 × 6(− 5), where all possible spin states were calculated, the trends are not perfectly linear between energy and spin state, however, the general trend that the highest spin state have the lowest energy is still evident. These results indicate that for all three GR cluster models, the high spin state is thermodynamically the most favorable state, which is also consistent with Hund’s rule for single atoms.

Fig. 2

Cluster free energy differences between spin states for the GR clusters. GR2 × 6(− 5): from left to right are spin states with unpaired electrons from 1 to 33 (i.e., 1, 3, 5…). GR2 × 6(− 6), GR3 × 6(− 10) and GR4 × 6(− 9): from left to right are low, intermediate and high spin states respectively. For all models, the energies of the highest spin states are set to be 0

A previous study showed that standard exchange–correlation functionals disfavour high spin state in iron complexes, whereas hybrid and some meta-GGA functionals are generally better at predicting the correct spin state [41]. Despite this, our calculations still show that our models prefer to be in high spin states. However, we still performed a comparison study of different functionals’ preference of spin state for our models. We compared BP86 with 6 other functionals, which represent three basically different approaches of DFT functionals: B-LYP [39, 42] and PBE [43] represent standard exchange–correlation functionals, B3-LYP [42, 44,45,46], PBE0 [47] and TPSSh [48] represent hybrid functionals, and TPSS [49] represents meta-GGA functional. We used GR2 × 6(− 8) as the example, calculated its fully optimized minimum energy at different spin states using different functionals with same TZVP basis set. Figure 3 shows the energy differences for GR2 × 6(− 8) when iron atoms are at low, intermediate and high spin states, with different functionals. For B-LYP and PBE, the intermediate spin state has the minimum energy, where the energy differences between intermediate and high spin state are very small. All hybrid functionals gave good support for the high spin state, especially PBE0, for which the high spin state has an energy that is 5.11 eV (i.e., 117.84 kcal/mol) lower than the intermediate spin state. B-P only slightly disfavour high spin state compared to hybrid functionals, which is consistent with the previous study [41], however, the energy of high spin state is still significantly lower than the energy of low and intermediate spin state, qualitatively reproducing the spin behaviour of more accurate hybrid functionals. Compared to B-P, TPSS disfavours high spin state considerably more. In our case, the B-P functional show good quality for predicting iron spin state, in qualitative agreement with the more accurate hybrid methods. The use of the BP functional also allows us to apply solvation treatment using COSMO-RS theory. [50]

Fig. 3

Energy differences between low, intermediate and high spin states of GR2 × 6(− 8) calculated using different DFT functionals

Results and discussion

Electronic properties of the single layer GR3 × 6 model

For the GR3 × 6 cluster, 20 different cluster charge models were calculated. Models with cluster charge from -12 to 0, i.e., Fe(II)/Fe(III) ratios are from 15/4 to 3/16, were obtained by carrying out BP/TZVP optimizations, while models with charges from -13 to -16 yielded distorted structures, and models with charge from + 1 to + 3 exhibited electronic occupations that were not in the ground state. GR3 × 6 with charges from -13 to -16 and from + 1 to + 3 were thus not analyzed here, but are further discussed in the supporting information, Additional file 1: Text S1 and Figure S1. In the following section, we will use the GR3 × 6 model with applied -10 charge, GR3 × 6(− 10), as an example to discuss the structure and electronic property of the single layer GR model. For this cluster, the Fe(II)/Fe(III) ratio is approximately 2.2 which matches the Fe(II)/Fe(III) ratios typically observed for experimentally produced sulphate GR [3].

The optimized GR3 × 6(− 10) structure is hexagonal, with all Fe atoms approximately in the same plane (Fig. 4). Given that Fe(III) has more unpaired electrons than Fe(II), we can identify the location of Fe(II) and Fe(III) ions in the structure by assessing the number of unpaired electrons (spin states) in every Fe atom using Mulliken population analysis. The results of this analysis for GR3 × 6(− 10) (Fig. 5) show that Fe atoms at the edges of the simulated GR layer have a higher number of unpaired electrons (i.e., higher spin states) compare to Fe atoms in the interior of the GR layer, indicated by a gap of 0.1. Additionally, the central Fe atom seems to have the lowest number of unpaired electron (i.e., lowest spin state). Thus, the reducing capacity of outer Fe atoms is lower compared to inner Fe atoms, with the central Fe atom having the highest reducing capacity. This means that the reductive capacity across the GR sheet structure is quite uneven, with the inner Fe atoms having Fe(II) like character, while the outer Fe atoms are more Fe(III) like.

Fig. 4

a Top and b side view of GR3 × 6 with a charge of -10

Fig. 5

Number of unpaired electrons for Fe atoms in the GR3 × 6(− 10) model structure, ordered according to decreasing number of unpaired electrons. Fe atoms labeled with O, I and C refer to outer, inner and central atoms, respectively

To verify the observations made for GR3 × 6(− 10), the distribution of unpaired electrons was also calculated for all other stable GR3 × 6 models (i.e., with charge from -12 to 0 i.e., Fe(II)/Fe(III) ratio from 15/4 to 3/16) shown in Fig. 6. Identical to the results of GR3 × 6(− 10), the number of unpaired electrons is considerably higher (i.e., spin states are higher) for outer Fe atoms compared to inner Fe atoms in all 13 model structures. Note that the clear drop in spin state between outer and inner Fe atoms is also still visible when averaging the number of unpaired electrons across all outer and inner Fe atoms, respectively, or by comparing minimum and maximum values for outer and inner Fe atoms, respectively (Table 2). Also, in most models (except for GR3 × 6(− 11) and GR3 × 6(− 12)), the central Fe atom still exhibits the lowest spin state as initially observed for GR3 × 6(− 10) (Fig. 6, Table 2). Overall, all 13 GR3 × 6 models reaffirmed the uneven distribution of reducing capacity (i.e., Fe(II) and Fe(III) atoms) across the single GR sheet. We note that this trend was consistent across the wide range of cluster charge (Fe(II)/Fe(III) ratio) we investigated.

Fig. 6

Numbers of unpaired electrons for Fe irons atoms in all obtained GR3 × 6 models. Fe atoms labeled with O, I and C refers to outer, inner and central atoms, respectively

Table 2 Relative differences in the number of unpaired electrons (NUE) between outer, inner and central Fe atoms for the simulated GR3 × 6 models

Cluster size effects

To ensure results observed for the GR3 × 6 models are not dependent on cluster size, we performed similar electronic property analyses for a smaller (GR2 × 6) and a larger (GR4 × 6) GR single layer structure, shown in Fig. 7a, b. For those two cluster sizes, 8 and 38 Fe(II)/Fe(III) ratios were examined, and 5 and 23 yielded thermodynamically stable structures after geometry optimization, respectively.

Fig. 7

Number of unpaired electrons for the Fe atoms in the stable a GR2 × 6 models (all 5 shown) and b GR4 × 6 models (only 5 shown for clarity, all 23 models are shown in Additional file 1: Figure S2). Fe atoms labeled with O, I and C refer to outer, inner and central atoms, respectively

The spin state distributions calculated for the GR2 × 6 and GR4 × 6 models showcase the same drop in the number of unpaired electrons between outer and central iron atoms as shown for the GR3 × 6 model. This is also verified by comparison of minimum and maximum spin values, or average values (Tables S2 and S4 in Additional file 1) has done for the GR3 × 6 model. Thus all three different sized models support the same electronic configuration for GR hydroxide sheets, with Fe(II) like atoms positioned in the middle of the GR sheet and Fe(III) like atoms on GR sheet edges. The consistency in our observations between different investigated Fe(II)/Fe(III) ratios and cluster sizes strongly suggest that our models provide a reasonable description of the electronic properties of a single GR hydroxide sheet.

Spin contamination

In open-shell system, unrestricted calculations have different set of orbitals for alpha electrons and beta electrons. In this case, spin contamination can occur because the expected value of total spin angular momentum operator, S2, is no longer commensurate with expected spin state of the system [51]. A large spin contamination indicates a failure of the applied computational method. However, because S2 is a two-electron operator, its meaning in density functional theory is diffuse. A previous study argued that S2 can be assigned a diagnostic value also within DFT [52, 53], and we would like to provide the details of spin contamination in our green rust models.

The spin contamination for the three different sized GR models was determined as the difference between the calculated S2 and the ideal value of S(S + 1), i.e., ΔS2 = S2-S(S-1). Calculated ΔS2 values (Table 3) are very small, especially considering the high spin states in the GR models calculated here. It is worth mentioning that the spin contamination increases as the GR model size is increased, but it decreases with as increase in multiplicities within the same sized GR model. Nevertheless, the highest ΔS2 value of 0.25 obtained for GR4 × 6(− 16) is still negligible considering its large number of unpaired electrons.

Table 3 ΔS2 values for all stable GR2 × 6, GR3 × 6 and GR4 × 6 models calculated using B-P functional

We also evaluated the effect of different DFT functionals on spin contamination using the GR2 × 6(− 8) model as an example (Table 4), second order Møller-Plesset (MP2) perturbation theory [54,55,56,57] also applied to determine the spin contamination. All methods show same result of a low degree of spin contamination, with the B-P functional only slightly higher than the hybrid functionals.

Table 4 ΔS2 values of GR2 × 6(− 8) with different DFT functionals and MP2 method. All the results were obtained with fully optimized geometry for the corresponding method

Mulliken and natural bond orbital analysis

Mulliken method is known as the oldest and cheapest way to obtain spin density and atomic charges, but it is also criticised for its lacking polarization effects and basis set dependence [58, 59]. Natural bond orbital (NBO) analysis is considered to be more reliable because it takes electron density and polarization effects into account. Hence, we evaluated the suitability of both the Mulliken and the NBO population analysis. The Mulliken unpaired electrons and NBO spin density for Fe atoms in the three GR models, GR2 × 6(− 8), GR3 × 6(− 10) and GR4 × 6(− 10), is shown in Fig. 8 (at B-P/TZVP level). Although absolute values between the Mulliken and NBO spin density differ, the spin density distributions among iron atoms in our models are basically identical.

Fig. 8

The comparison of spin density of Fe atoms obtained by Mulliken and Natural Bond Orbital (NBO) population analysis for a GR2 × 6(− 8), b GR3 × 6(− 10) and c GR4 × 6(− 10)

Functional effects

The spin density distributions were calculated for GR2 × 6(− 8) and GR3 × 6(− 8) using the different DFT functionals discussed earlier (“Multiplicities assumption” section) and they are shown in the Additional file 1: Figures S3 (GR2 × 6(− 8) with different functionals), S4 (GR2 × 6(− 8) with MP2) and S5 (GR3 × 6(− 10) with PBE0). Two types of spin density distributions are observed: (1) The pure functionals, BLYP, PBE and TPSS, produced similar distributions as the B-P functional in that the center Fe atoms show the lowest spin density, while the outer Fe atoms show relatively high spin states; (2) The hybrid functionals, B3LYP, TPSSh, PBE0 and MP2, include exact exchange, and they show only two discreet values (except for TPSSh) as opposed to the smeared out electronic distribution of pure DFT functionals. These two values can be identified as Fe(III) and Fe(II). All high spin Fe(III) are located on GR sheet edges, meaning that all inner Fe atoms are Fe(II).

Although the pure DFT functionals do not localize the electrons correctly, the distribution trends agree with trends observed for hybrid functionals, which localize the electrons accurately into Fe(II) and Fe(III). Additionally in Additional file 1: Figures S3 and S4, the two Fe atoms that are predicted to be Fe(III) by the hybrid functionals are the two atoms with the highest calculated spin in the pure DFT calculations. Thus, all functionals investigated (pure and hybrid) reveal the same tendency that the outer Fe atoms on average have higher spin state compared to inner Fe atoms. The calculated spin density of the larger cluster GR3 × 6(− 10) with PBE0 also shows a consistent performance, with all six Fe(III) located on the edge of the cluster. If the Fe(III) were randomly distributed, there would only be a 3% likelihood that all Fe(III) would happen to be located at the edge sites, strengthening our conclusions about Fe(III) preferring to be locate at the edge sites.

Conclusions and implications

In this work, three different sized cluster models of single GR hydroxide layers were built to study their electronic properties using density functional theory. The calculations showed that a high spin state is favored thermodynamically for the single layer GR model and that a minimum amount of Fe(III) are required, i.e., Fe(II)/Fe(III) ratios between 0.2 and 5, to maintain the hexagonal shape of the GR structure. These ratios are in agreement with ratios measured of synthetic and natural GR samples that typically range between 2 and 3 [3]. All three cluster models showed that the spin states of edge Fe atoms are significantly higher compared to Fe atoms located in the GR interior. This in turn means that Fe(II) and Fe(III) ions are unevenly distributed across the GR hydroxide sheet structure, with edge Fe atoms being more Fe(III)-like, while inner Fe atoms are more Fe(II)-like; thus the GR interior is more reducing compared to the edge. The trends we found were consistent among several pure DFT calculations, hybrid functionals as well as MP2. A reason that may explain why we observed this uneven Fe(II)/Fe(III) distribution is that the chemical environment between edge and inner Fe atoms is different. The edge atoms have larger degree of flexibility compared to the more constrained inner atoms, which are more or less fixed in crystallographic positions. This could allow for a more unconstrained relaxation, including the electronic degrees of freedom. Furthermore, the edge iron atoms bind to three or four hydroxyl ions, instead of bridging oxygen between iron atoms as done by the inner iron atoms. The different chemical environment creates a nonisotropic iron site, which could affect the electronic properties as well, including via a local dipole effect.

Evidently, this GR model excluds the presence of charge-balancing interlayer anions, which potentially could affect the local Fe(II)/Fe(III) distribution, i.e., electronic properties, across the GR hydroxide sheet. However, it is still worth making some initial comparisons to experimental observations here. Taking the example of the reduction of chromate, Cr(VI), by sulphate GR, some studies argued that Cr(VI) is reduced in the interlayer after sulphate exchange with chromate [25, 29,30,31], which would support our observations. In contrast, others argued for Cr(VI) reduction to occur at GRSO4 particle edges [26]. Similarly, the reductions of U(VI) [15] and C2Cl4 [21] by GR are argued to occur on the GR edges which seems to be contrary to our finding that the reduction capacity tends to be higher in the interior of GR crystals. However, the charge hopping mechanism proposed by Wander et al. [32] could explain why reduction can occur on the edge, even if it is less reactive than the central part of GR. Indeed, electron microscopy images of Cr(VI) reacted sulphate GR often show reaction rims with oxidized Fe (oxyhydr) oxide phases, while the GR interiors seem to get fully dissolved [25], which may be another indication that GR interiors are more reactive, i.e., more reducing, and that electrons can readily transfer from GR interior to its edges.

Lastly, our study also suggests that the B-P functional with a TZVP basis set gives reasonable results for the electronic distribution in GR, which opens up possibilities for modelling solvation behavior of interlayer molecules using the implicit solvent method COSMO-RS [50].

Supplementary Information

The full list of Fe(II)/Fe(III) ratios, the number of unpaired electrons in every Fe atoms of GR2 × 6, GR3 × 6 and GR4 × 6, and discussion of Fe(II)/Fe(III) ratio and structure stability are shown in additional file. The discussion about the effect of dielectric constant on the Fe spins (Additional file 1: Text S2, Tables S5–S7 and Figure S6), and magnetic property analysis (Additional file 1: Text S3 and Tables S8–S11) are can also be seen in Additional file 1.



Green rust


Sulphate green rust


Layered double hydroxide


Number of unpaired electrons


Density Functional Theory


Conductor like screening model,


Natural bond orbital


  1. 1.

    Bernal JD, Dasgupta DR, Mackay AL (1959) The oxides and hydroxides of iron and their structural inter-relationships. Clay Miner 4:15–30

    Article  Google Scholar 

  2. 2.

    Abdelmoula M, Trolard F, Bourrié G, Génin JMR (1998) Evidence for the Fe(II)-Fe(III) green rust “Fougerite” mineral occurrence in a hydromorphic soil and its transformation with depth. Hyperfine Interact 112:235–238

    Article  Google Scholar 

  3. 3.

    Christiansen BC, Balic-Zunic T, Petit PO, Frandsen C, Mørup S, Geckeis H, Stipp KA, SLS. (2009) Composition and structure of an iron-bearing, layered double hydroxide (LDH)—Green rust sodium sulphate. Geochim Cosmochim Ac 73:3579–3592

    Article  Google Scholar 

  4. 4.

    Zegeye A, Bonneville S, Benning LG, Sturm A, Fowle DA, Jones C et al (2012) Green rust formation controls nutrient availability in a ferruginous water column. Geology 40:599–602

    Article  Google Scholar 

  5. 5.

    Christiansen BC, Dideriksen K, Katz A, Nedel S, Bovet N, Sorensen HO et al (2014) Incorporation of monovalent cations in sulfate green rust. Inorg Chem 53:8887–8894

    Article  Google Scholar 

  6. 6.

    Ruby C, Aïssa R, Géhin A, Cortot J, Abdelmoula M, Génin J-M (2006) Green rusts synthesis by coprecipitation of FeII–FeIII ions and mass-balance diagram. C R Geosci 338:420–432

    Article  Google Scholar 

  7. 7.

    Refait P, Génin JMR (1997) The mechanisms of oxidation of ferrous hydroxychloride β-Fe2(OH)3Cl in aqueous solution: The formation of akaganeite vs goethite. Corros Sci 39:539–553

    Article  Google Scholar 

  8. 8.

    Génin JMR, Refait P, Simon L, Drissi SH (1998) Preparation and Eh–pH diagrams of Fe(II)–Fe(III) green rust compounds; hyperfine interaction characteristics and stoichiometry of hydroxy-chloride, -sulphate and -carbonate. Hyperfine Interact 111:313–318

    Article  Google Scholar 

  9. 9.

    Refait PH, Abdelmoula M, GÉnin JMR (1998) Mechanisms of formation and structure of green rust one in aqueous corrosion of iron in the presence of chloride ions. Corrosion Sci 40:1547–1560

    Article  Google Scholar 

  10. 10.

    Génin JMR, Christy A, Kuzmann E, Mills S, Ruby C (2014) Structure and occurrences of green rust related new minerals of the fougérite group, trébeurdenite and mössbauerite, belonging to the hydrotalcite supergroup; how Mössbauer spectroscopy helps XRD. Hyperfine Interact 226:459–482

    Article  Google Scholar 

  11. 11.

    Refait P, Bon C, Simon L, Bourrié G, Trolard F, Bessière J et al (2018) Chemical composition and Gibbs standard free energy of formation of Fe(II)-Fe(III) hydroxysulphate green rust and Fe(II) hydroxide. Clay Miner 34:499–510

    Article  Google Scholar 

  12. 12.

    Génin JMR, Olowe AA, Refait P, Simon L (1996) On the stoichiometry and pourbaix diagram of Fe(II)-Fe(III) hydroxy-sulphate or sulphate-containing green rust 2: An electrochemical and Mössbauer spectroscopy study. Corros Sci 38:1751–1762

    Article  Google Scholar 

  13. 13.

    Refait P, Simon L, Génin J-MR (2000) Reduction of SeO42-anions and anoxic formation of Iron(II)−Iron(III) Hydroxy-selenate green rust. Environ Sci Technol 34:819–825

    Article  Google Scholar 

  14. 14.

    Hansen HCB, Guldberg S, Erbs M, Koch CB (2001) Kinetics of nitrate reduction by green rusts—effects of interlayer anion and Fe(II): Fe(III) ratio. Appl Clay Sci 18:81–91

    Article  Google Scholar 

  15. 15.

    O’Loughlin EJ, Kelly SD, Cook RE, Csencsits R, Kemner KM (2003) Reduction of uranium(VI) by mixed iron(II)/iron(III) hydroxide (green rust): formation of UO2 nanoparticles. Environ Sci Technol 37:721–727

    Article  Google Scholar 

  16. 16.

    O’Loughlin EJ, Kemner KM, Burris DR (2003) Effects of Ag(I), Au(III), and Cu(II) on the reductive dechlorination of carbon tetrachloride by green rust. Environ Sci Technol 37:2905–2912

    Article  Google Scholar 

  17. 17.

    Pepper SE, Bunker DJ, Bryan ND, Livens FR, Charnock JM, Pattrick RA et al (2003) Treatment of radioactive wastes: an X-ray absorption spectroscopy study of the reaction of technetium with green rust. J Colloid Interf Sci 268:408–412

    Article  Google Scholar 

  18. 18.

    O’Loughlin EJ, Kelly SD, Kemner KM, Csencsits R, Cook RE (2003) Reduction of AgI, AuIII, CuII, and HgII by FeII/FeIII hydroxysulfate green rust. Chemosphere 53:437–446

    Article  Google Scholar 

  19. 19.

    Trolard F, Bourrié G (2006) Structure of fougerite and green rusts and a thermodynamic model for their stabilities. J Geochem Explor 88:249–251

    Article  Google Scholar 

  20. 20.

    Liang X, Philp RP, Butler EC (2009) Kinetic and isotope analyses of tetrachloroethylene and trichloroethylene degradation by model Fe(II)-bearing minerals. Chemosphere 75:63–69

    Article  Google Scholar 

  21. 21.

    Choi J, Batchelor B, Chung J (2010) Reductive dechlorination of tetrachloroethylene by green rusts modified with copper. Water, Air, & Soil Poll 212:407–417

    Article  Google Scholar 

  22. 22.

    Liang X, Butler EC (2010) Effects of natural organic matter model compounds on the transformation of carbon tetrachloride by chloride green rust. Water Res 44:2125–2132

    Article  Google Scholar 

  23. 23.

    Ayala-Luis KB, Cooper NG, Koch CB, Hansen HC (2012) Efficient dechlorination of carbon tetrachloride by hydrophobic green rust intercalated with dodecanoate anions. Environ Sci Technol 46:3390–3397

    Article  Google Scholar 

  24. 24.

    Williams AG, Scherer MM (2001) Kinetics of Cr(VI) reduction by carbonate green rust. Environ Sci Technol 35:3488–3494

    Article  Google Scholar 

  25. 25.

    Skovbjerg LL, Stipp SLS, Utsunomiya S, Ewing RC (2006) The mechanisms of reduction of hexavalent chromium by green rust sodium sulphate: Formation of Cr-goethite. Geochim Cosmochim Ac 70:3582–3592

    Article  Google Scholar 

  26. 26.

    Thomas AN, Eiche E, Göttlicher J, Steininger R, Benning LG, Freeman HM et al (2018) Products of hexavalent chromium reduction by green rust sodium sulfate and associated reaction mechanisms. Soil Sys 2:58

    Article  Google Scholar 

  27. 27.

    Thomas AN, Eiche E, Gottlicher J, Steininger R, Benning LG, Freeman HM et al (2020) Effects of metal cation substitution on hexavalent chromium reduction by green rust. Geochem Trans 21:2

    Article  Google Scholar 

  28. 28.

    Usman M, Byrne JM, Chaudhary A, Orsetti S, Hanna K, Ruby C et al (2018) Magnetite and green rust: synthesis, properties, and environmental applications of mixed-valent iron minerals. Chem Rev 118:3251–3304

    Article  Google Scholar 

  29. 29.

    Bond DL, Fendorf S (2003) Kinetics and structural constraints of chromate reduction by green rusts. Environ Sci Technol 37:2750–2757

    Article  Google Scholar 

  30. 30.

    Loyaux-Lawniczak S, Refait P, Ehrhardt J-J, Lecomte P, Génin J-MR (2000) Trapping of Cr by formation of ferrihydrite during the reduction of chromate ions by Fe(II)−Fe(III) hydroxysalt green rusts. Environ Sci Technol 34:438–443

    Article  Google Scholar 

  31. 31.

    Loyaux-Lawniczak S, Refait P, Lecomte P, Ehrhardt JJ, Génin JMR (1999) The reduction of chromate ions by Fe(II) layered hydroxides. Hydrol Earth Syst Sci 3:593–599

    Article  Google Scholar 

  32. 32.

    Wander MCF, Rosso KM, Schoonen MAA (2007) Structure and charge hopping dynamics in green rust. J Phys Chem C 111:11414–11423

    Article  Google Scholar 

  33. 33.

    Mangayayam MC, Dideriksen K, Tobler DJ (2018) Can or cannot green rust reduce chlorinated ethenes? Energy Procedia 146:173–178

    Article  Google Scholar 

  34. 34.

    Ai J, Yin W, Hansen HCB (2019) Fast dechlorination of chlorinated ethylenes by green rust in the presence of bone char. Environ Sci Technol Lett 6:191–196

    Article  Google Scholar 

  35. 35.

    Ai J, Ma H, Tobler DJ, Mangayayam MC, Lu C, van den Berg FWJ et al (2020) Bone char mediated dechlorination of trichloroethylene by green rust. Environ Sci Technol 54:3643–3652

    Article  Google Scholar 

  36. 36.

    Ahlrichs R, Bär M, Häser M, Horn H, Kölmel C (1989) Electronic structure calculations on workstation computers: the program system turbomole. Chem Phys Lett 162:165–169

    Article  Google Scholar 

  37. 37.

    Klamt A, Schüürmann G (1993) COSMO: a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J Chem Soc, Perkin Trans 2:799–805

    Article  Google Scholar 

  38. 38.

    Perdew JP (1986) Density-functional approximation for the correlation energy of the inhomogeneous electron gas. Phys Rev B Condens Matter 33:8822–8824

    Article  Google Scholar 

  39. 39.

    Becke AD (1988) Density-functional exchange-energy approximation with correct asymptotic behavior. Phys Rev A Gen Phys 38:3098–3100

    Article  Google Scholar 

  40. 40.

    Schäfer A, Horn H, Ahlrichs R (1992) Fully optimized contracted Gaussian basis sets for atoms Li to Kr. J Chem Phys 97:2571–2577

    Article  Google Scholar 

  41. 41.

    Swart M, Groenhof AR, Ehlers AW, Lammertsma K (2004) Validation of exchange−correlation functionals for spin states of iron complexes. J Phys Chem A 108:5479–5483

    Article  Google Scholar 

  42. 42.

    Lee C, Yang W, Parr RG (1988) Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B Condens Matter 37:785–789

    Article  Google Scholar 

  43. 43.

    Perdew JP, Burke K, Ernzerhof M (1996) Generalized gradient approximation made simple. Phys Rev Lett 77:3865–3868

    Article  Google Scholar 

  44. 44.

    Vosko SH, Wilk L, Nusair M (1980) Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Can J Phys 58:1200–1211

    Article  Google Scholar 

  45. 45.

    Becke AD (1993) Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys 98:5648–5652

    Article  Google Scholar 

  46. 46.

    Stephens PJ, Devlin FJ, Chabalowski CF, Frisch MJ (1994) Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J Phys Chem 98:11623–11627

    Article  Google Scholar 

  47. 47.

    Adamo C, Barone V (1999) Toward reliable density functional methods without adjustable parameters: The PBE0 model. J Chem Phys 110:6158–6170

    Article  Google Scholar 

  48. 48.

    Staroverov VN, Scuseria GE, Tao J, Perdew JP (2003) Comparative assessment of a new nonempirical density functional: molecules and hydrogen-bonded complexes. J Chem Phys 119:12129–12137

    Article  Google Scholar 

  49. 49.

    Tao J, Perdew JP, Staroverov VN, Scuseria GE (2003) Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids. Phys Rev Lett 91:146401

    Article  Google Scholar 

  50. 50.

    Klamt A, Eckert F, Arlt W (2010) COSMO-RS: an alternative to simulation for calculating thermodynamic properties of liquid mixtures. Annu Rev Chem Biomol Eng 1:101–122

    Article  Google Scholar 

  51. 51.

    Szabo A, Ostlund NS (1996) Modern quantum chemistry: introduction to advanced electronic structure theory. Courier Corporation, North Chelmsford

    Google Scholar 

  52. 52.

    Johansson MP, Sundholm D, Gerfen G, Wikstrom M (2002) The spin distribution in low-spin iron porphyrins. J Am Chem Soc 124:11771–11780

    Article  Google Scholar 

  53. 53.

    GrÄFenstein J, Cremer D (2001) On the diagnostic value of (Ŝ2) in Kohn-Sham density functional theory. Molecular Phys 99:981–989

    Article  Google Scholar 

  54. 54.

    Head-Gordon M, Pople JA, Frisch MJ (1988) MP2 energy evaluation by direct methods. Chem Phys Lett 153:503–506

    Article  Google Scholar 

  55. 55.

    Weigend F, Häser M, Patzelt H, Ahlrichs R (1998) RI-MP2: optimized auxiliary basis sets and demonstration of efficiency. Chem Phys Lett 294:143–152

    Article  Google Scholar 

  56. 56.

    Frisch MJ, Head-Gordon M, Pople JA (1990) A direct MP2 gradient method. Chem Phys Lett 166:275–280

    Article  Google Scholar 

  57. 57.

    Frisch MJ, Head-Gordon M, Pople JA (1990) Semi-direct algorithms for the MP2 energy and gradient. Chem Phys Lett 166:281–289

    Article  Google Scholar 

  58. 58.

    Mulliken RS (1955) Electronic population analysis on LCAO–MO molecular wave functions. I J Chem Phys 23:1833–1840

    Article  Google Scholar 

  59. 59.

    Toubal K, Boukabcha N, Tamer Ö, Benhalima N, Altürk S, Avcı D et al (2017) Spectroscopic (FT-IR,1H and 13C NMR) characterization and density functional theory calculations for (Z)-5-(4-nitrobenzyliden)-3-N(2-ethoxyphenyl)-2-thioxo-thiazolidin-4-one (ARNO). J Mol Struc 1147:569–581

    Article  Google Scholar 

Download references


We thank Evren Ataman and Akin Budi for support in calculations and valuable discussions.


This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 675219.

Author information




WS performed the DFT calculations. MPA supervised the computational method and DJT supervised the green rust implications. All authors assisted with writing the manuscript and approved its submission. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Weichao Sun.

Ethics declarations

Competing interests

There are no conflicts to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional tables and figures.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sun, W., Tobler, D.J. & Andersson, M.P. A density functional theory study of Fe(II)/Fe(III) distribution in single layer green rust: a cluster approach. Geochem Trans 22, 3 (2021).

Download citation


  • Green Rust
  • Single Layer
  • Density Functional Theory
  • Fe(II)/Fe(III) Distribution
  • Spin Density
  • Cluster Approach