Introduction
Molybdenum (Mo) occurs as a trace element in most soils but it is essential for nearly all known life forms on Earth (Neunhäuser et al. Reference Neunhäuser, Berreck and Insam2001; Mendel and Bittner Reference Mendel and Bittner2006). Both in the global biogeochemical oxidation-reduction cycles of carbon, nitrogen, and sulfur (Stiefel Reference Stiefel2002; Williams and Frausto da Silva Reference Williams and Frausto da Silva2002) and in the metabolism of living organisms (Bortels Reference Bortels1930) Mo holds a key position in the active sites of a wide range of metalloenzymes (Norrish Reference Norrish, Nicholas and Egan1975; Mendel and Bittner Reference Mendel and Bittner2006).
The most abundant and only bioavailable Mo species in environmental systems is Mo6+ , in the form of the molybdate anion (MoO4 2−) (Mendel and Bittner Reference Mendel and Bittner2006). In contrast to the oceans, where Mo concentrations are up to 10.3 mg/kg (Wang et al. Reference Wang, Aller and Sañudo-Wilhelmy2011), the soil contains only 1–5 mg of Mo kg–1 (He et al. Reference He, Yang and Stoffella2005). Most soil-bound Mo is considered to not be bioavailable (Žemberyová et al. Reference Žemberyová, Hagarová, Zimová, Barteková and Kuss2010) because of surface adsorption (Goldberg et al. Reference Goldberg, Forster and Godfrey1996) or structural incorporation (Richmond et al. Reference Richmond, Hockridge, Loan and Parkinson2004; Carroll and Richmond Reference Richmond, Hockridge, Loan and Parkinson2008; Das et al. Reference Das, Essilfe-Dughan and Hendry2016) in the soil minerals. Various studies have emphasized the adsorption mechanisms of Mo to different soil constituents (e.g. Goldberg et al. Reference Goldberg, Forster and Godfrey1996; Brinza et al. Reference Brinza, Benning and Statham2008), whereas information about the structural incorporation of Mo in minerals remains scarce (Richmond et al. Reference Richmond, Hockridge, Loan and Parkinson2004; Carroll and Richmond Reference Carroll and Richmond2008; Žemberyová et al. Reference Žemberyová, Hagarová, Zimová, Barteková and Kuss2010; Wang et al. Reference Žemberyová, Hagarová, Zimová, Barteková and Kuss2011; Brinza et al. Reference Brinza, Vu, Shaw, Mosselmans and Benning2015; Das et al. Reference Das, Essilfe-Dughan and Hendry2016).
Fixation of Mo in German forest soils by ferric oxides was documented by Lang and Kaupenjohann (Reference Lang and Kaupenjohann1999) who showed that most Mo in these soils is not adsorbed but rather is associated with crystalline iron oxides (135–1071 µg/kg). The easily exchangeable Mo fraction amounts to only 5–28 µg/kg and the poorly crystalline iron oxides retain 44–407 µg/kg. Furthermore, shortly after Mo fertilization, forest soils retain molybdenum by surface adsorption to iron oxides, while, with increasing time, Mo is incorporated into the iron oxide structure and thus becomes inaccessible for further biogeochemical processes (Lang Reference Lang2000).
Common iron oxide mineral phases in soils are goethite (α-FeOOH), lepidocrocite (γ-FeOOH), feroxyhyte (δ-FeOOH), hematite (α-Fe2O3), maghemite (γ-Fe2O3), magnetite (Fe3O4), and ferrihydrite (FeOOH·nH2O). Not only do they have large adsorption capacities for trace elements in solution (Peacock and Sherman Reference Peacock and Sherman2004; Rout et al. Reference Rout, Mohapatra and Anand2012) but also for element uptake via coprecipitation (Martínez and McBride Reference Martínez and McBride1998, Reference Martínez and McBride2001; Manceau et al. Reference Manceau, Schlegel, Musso, Sole, Gauthier, Petit and Trolard2000). These properties make them well known scavengers for metals (Dardenne et al. Reference Dardenne, Schäfer, Lindqvist-Reis, Denecke, Plaschke, Rothe and Kim2002; Tarassov et al. Reference Tarassov, Mihailova, Tarassova and Konstantinov2002; Mullet et al. Reference Mullet, Demoission, Humbert, Michot and Vantelon2007) and metalloids (e.g. Sieber et al. Reference Sieber, Sanchez, Turner and Somorjai1985; Scheinost et al. Reference Scheinost, Rossberg, Vantelon, Xifra, Kretzschmar, Leuz and Johnson2006; Bolanz et al. Reference Bolanz, Bläss, Ackermann, Ciobota, Rösch, Tarcea, Popp and Majzlan2013b) in both natural and contaminated environments (Hochella et al. Reference Hochella, Kasama, Putnis, Putnis and Moore2005). Iron oxide phases can, therefore, control the bioavailability of nutrients as well as contaminants such as arsenic (Sherman and Randall Reference Sherman and Randall2003; Bolanz et al. Reference Bolanz, Wierzbicka-Wieczorek, Čaplovičová, Uhlík, Göttlicher, Steininger and Majzlan2013a, Reference Bolanz, Bläss, Ackermann, Ciobota, Rösch, Tarcea, Popp and Majzlanb) and lead (Vu et al. Reference Vu, Shaw, Brinza and Benning2010) in soils.
Previous studies focused predominantly on homovalent and homotopic substitutions, where Fe3+ is replaced by other trivalent cations such as Al3+, Cr3+, or Mn3+ (Cornell and Giovanoli Reference Cornell and Giovanoli1987; Wells et al. Reference Wells, Gilkes and Fitzpatrick2001; Mullet et al. Reference Mullet, Demoission, Humbert, Michot and Vantelon2007). The substitution by aliovalent cations such as Mo6+, however, depends on intrinsic charge-balancing reactions. Several mechanisms were described or assumed in previous investigations, such as the formation of Fe3+ vacancies (Ciobanu et al. Reference Ciobanu, Wade, Cook, Schmidt Mumm and Giles2013), reduction of surrounding Fe3+ to Fe2+ (Balko and Clarkson Reference Balko and Clarkson2001), the formation of O2– vacancies (Khan et al. Reference Khan, Chen, Boolchand and Smirniotis2008), and the formation of structural clusters intergrown within the iron oxide host phase (Bolanz et al. Reference Bolanz, Bläss, Ackermann, Ciobota, Rösch, Tarcea, Popp and Majzlan2013b). For Mo6+, the most abundant Mo species in oxic systems (Barron et al. Reference Barron, Wurzburger, Bellenger, Wright, Kraepiel and Hedin2009), almost no convincing data for structural incorporation in iron oxides exist, however.
The objective of the current study was to investigate the mechanisms of structural incorporation of Mo in ferrihydrite, goethite, and hematite. Synthesis protocols for these Mo-bearing iron oxides needed to be developed, therefore, and the resulting products needed to be analyzed by XAS. The research focused on: (1) developing reproducible methods to synthesize Mo-bearing iron oxides; (2) evaluating the mechanisms of structural incorporation of Mo6+ in these oxides; and (3) assessing the effects of Mo6+ incorporation on the structures of iron oxides.
METHODS
Synthesis of Pure and Mo-bearing Iron Oxides
The investigated iron oxides were synthesized in the absence or presence of Mo6+ according to modified procedures described by Schwertmann and Cornell (Reference Schwertmann and Cornell2000). All of the chemicals (which have the brand name Emsure®) were supplied by Merck (Germany) and met ACS, Reag. Ph. Eur., and ISO guidelines for analytical-grade chemicals.
For the synthesis of pure six-line ferrihydrite, 20 g of iron(III) nitrate nonahydrate was dissolved in 1.5 L deionized water preheated to 75°C, under rapid stirring. Immediately afterward, the solution was placed in an oven for 10 min at 75°C. During this time, the solution changed from amber-colored to dark reddish brown. The dark brown solution was plunged into ice water for rapid cooling. For Mo-bearing six-line ferrihydrite, selected volumes of a Na2MoO4⋅2H2O solution adjusted to pH 4 with diluted HNO3 were added to previously-made iron(III) nitrate nonahydrate solution. The Mo:Fe molar ratios used were 0 (blank), 0.01, 0.02, 0.05, and 0.10.
The formation of hematite was based on the transformation of two-line ferrihydrite and started with the synthesis of pure and Mo-bearing ferrihydrite. 20 g of iron(III) nitrate nonahydrate (Fe(NO3)3⋅9H2O) was dissolved in 300 mL of deionized water under vigorous stirring. For Mo-bearing two-line ferrihydrite, selected volumes of a Na2MoO4⋅2H2O solution, previously adjusted to pH 4 with diluted HNO3, were added, so that the following initial Mo:Fe molar ratios applied: 0 (blank), 0.01, 0.02, 0.05, 0.10, 0.20, 0.50, and 1.00. Subsequent addition of a 1 M KOH solution raised the pH to 7–8 and the solution changed from amber to a dark brown suspension of two-line ferrihydrite. These ferrihydrite suspensions were filtered and the still moist, thick filtrate was transferred into 100 mL pressure reaction vessels and placed in a preheated oven at 200°C for 7 days.
For the synthesis of goethite, a method based on the transformation of ferrihydrite was suggested by Schwertmann and Cornell (Reference Schwertmann and Cornell2000). There, 50 mL of a 1 M iron(III) nitrate nonahydrate solution was poured into a 1 L polyethylene flask and 90 mL of 5 M KOH solution was added quickly under constant stirring. Reddish brown ferrihydrite precipitated immediately and the suspension was diluted to 1 L with deionized water. Aging of this suspension in a closed flask at 70°C for 3 days led to the transformation of ferrihydrite to ochreous goethite. Molybdenum in various Mo:Fe molar ratios was introduced by adding 50 mL of a pH 4-adjusted Na2MoO4⋅2H2O solution to the 50 mL of previously made 1 M iron(III) nitrate nonahydrate solution. After addition of 90 mL of 5 M KOH solution, the resulting Mo-bearing ferrihydrite suspension was left to age in a preheated oven at 70°C for 3 days.
All resulting suspensions and solids were treated with 200 mL of a 0.2 M solution of K2HPO4 and KH2PO4 (according to Liang and Zhu Reference Liang and Zhu2016) in order to desorb molybdate ions attached to the outer surface of the iron oxide particles. Subsequently, the suspensions were washed by dialysis for about a week until the conductivity of the dialysis solution, deionized water, decreased to <10 µS/cm, then filtered and air-dried. After drying, the precipitates were ground for ~10 min with an agate mortar and pestle and sealed in glass vials. In the case of six-line ferrihydrite, before starting the desorption step, the resulting suspension was washed by dialysis, dried in an oven at 30°C, and the dry product then ground. The powder obtained was added to 200 mL of the 0.2 M K2HPO4 and KH2PO4 solution and subsequently treated like the above-mentioned suspensions. If not implemented in this order, six-line ferrihydrite in suspension would transform to two-line ferrihydrite during desorption.
Powder X-ray Diffraction Analysis (XRD)
Powder XRD patterns of all synthesized phases and reference compounds were collected with a Bruker D8 AXS Advance diffractometer (Bruker, Karlsruhe, Germany) using CuKα radiation (λ = 1.5406 Å) and equipped with a LYNXEYE detector. Samples of the pure iron oxides and the Mo-bearing samples were measured for full-profile refinement at 20°C from 5 to 120°2θ with a dwell time of 1 s and a step width of 0.01°2θ. The powder XRD patterns were processed using the software JANA2006 (Petříček et al. Reference Petříček, Dušek and Palatinus2014), and for the fitting, models of six-line ferrihydrite (Michel et al. Reference Michel, Ehm, Antao, Lee, Chupas, Liu, Strongin, Schoonen, Phillips and Parise2007), hematite (Blake et al. Reference Blake, Hessevick, Zoltai and Finger1966), goethite (Gualtieri and Venturelli Reference Gualtieri and Venturelli1999), and akaganéite (Post et al. Reference Post, Heaney, Von Dreele and Hanson2003) were used. The statistical measures of deviation of the fit from the data obtained are given in Table 1, including the residuals for the weighted profile (Rwp) and the goodness-of-fit (GOF).
Values below the detection limit are marked n.v. (no value). Values marked with * represent the estimated uncertainty for the Fe loss per mole of Mo. The numbers in parentheses indicate the uncertainty for the last digit. For each fit, the weighted R-factor (Rwp) and goodness of fit (GOF) are given
Inductively-coupled Plasma Optical Emission Spectrometry (ICP-OES)
A simultaneous radial ICP-OES spectrometer 725ES (Agilent Technologies, Waldbronn, Germany) with a charge-coupled device (CCD) detector was used to determine the chemical composition of all samples (Fe, Mo). Before measurements, each sample was dried at 105°C for 24 h. Subsequently, ~100 mg of each sample was dissolved at room temperature in 5 mL of concentrated hydrochloric acid (HCl), adjusted to a defined volume of 10 mL with deionized water and stored in corrosion-resistant sample tubes.
X-ray Absorption Near-edge Structure (XANES) and Extended X-ray Absorption Fine-structure (EXAFS) Spectroscopy
X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine-structure (EXAFS) spectra of the synthesized Mo-bearing iron oxides and reference materials (MoO2, MoO3, and Na2MoO4⋅2H2O) were collected at the Mo L 3 -edge (2520 eV) at the SUL-X beamline of the Karlsruhe Research Accelerator (KARA), Germany, and at the Mo K-edge (20,020 eV) at the B18 beamline of the Diamond Light Source (Didcot, UK). A silicon (111) crystal pair with a fixed beam exit as monochromator was used for XANES and EXAFS measurements at the SUL-X and B18 beamlines.
Using the software XAFSmass (Klementiev Reference Klementiev2012), a defined mass of each sample was calculated for optimum absorption and subsequently mixed with cellulose powder, homogenized, and pressed into pellets. This resulted in a Mo concentration of 2.710 wt.% for the six-line ferrihydrite, 2.733 wt.% for the hematite, and 0.013 wt.% for the goethite pellet.
The XANES spectra at the Mo L 3 -edge were collected at the SUL-X beamline with a step width of 5 eV in the region of –80 to –50 eV, 1 eV in the region –50 to –20 eV below the edge, 0.2 eV in the region –20 to + 20 eV relative to the edge, and 0.05 Å−1 up to 5.3 k (~105 eV) above the absorption edge. EXAFS spectra collected at the Diamond Light Source were recorded with a step width of 0.3 eV in the region –200 to 1000 eV relative to the edge (20,020 eV).
At KARA, all spectra were measured simultaneously in transmission and fluorescence mode, while spectra collected at Diamond were measured in transmission mode. A Mo0 foil was mounted between ionization chambers 2 and 3 at the beamlines. The Mo K- or L 3 -edge was recorded parallel to each sample spectrum for energy calibration. Five spectra were measured for each sample as a minimum and two spectra for each reference compound.
All EXAFS spectra were processed and fit by the Demeter software package (Ravel and Newville Reference Ravel and Newville2005). For the shell-by-shell fitting of the k 3 -weighted EXAFS spectra, paths generated from models of six-line ferrihydrite (Michel et al. Reference Michel, Ehm, Antao, Lee, Chupas, Liu, Strongin, Schoonen, Phillips and Parise2007), a hematite supercell based on the structure of Blake et al. (Reference Blake, Hessevick, Zoltai and Finger1966), and a goethite supercell based on the structure of Gualtieri and Venturelli (Reference Gualtieri and Venturelli1999) were used. In these models, Mo replaced one Fe atom. Parameters for the paths were calculated by FEFF6 (de Leon et al. Reference De Leon, Rehr, Zabinsky and Albers1991) from the Demeter software package (Ravel and Newville Reference Ravel and Newville2005).
The statistical measures of the deviation of the fit from data obtained are given, as well as the number of variables allowed (to float) in the fit (N var), the number of independent data points (N idp), and the R factor. While N idp is equal to 2ΔkR/π, where k and R are the fitting ranges, the R factor is defined as:
where χ is the magnitude of the EXAFS oscillations and x is the set of variables to be refined.
The EXAFS spectra of the references were fitted with the respective structures. Their spectra were used to identify problems that could be encountered when fitting a model to the structurally unknown samples.
Results
Chemical Composition of Pure and Mo-bearing Iron Oxides
The chemical compositions of all synthesized iron oxide phases, obtained by ICP-OES, showed an increase in the Mo content in the solid phase with the increasing initial Mo:Fe molar ratio in solution (Table 1).
Goethite takes up less Mo than any of the other phases studied. Synthesized at an initial Mo:Fe molar ratio of 1.0, goethite contained only 0.03 wt.% Mo. In comparison, six-line ferrihydrite or hematite synthesized at an initial Mo:Fe ratio of 0.01 contained much more Mo, 2.61 or 1.12 wt.%, respectively. Despite their low Mo concentrations, the goethite samples were included in this study because crystallographic changes in the samples could still be observed. For hematite with larger Mo concentrations, the molar loss of Fe for each mole of Mo (Fe loss) was calculated by (Bolanz et al. Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017):
where Fe wt.% Mo(0) is the Fe content of the Mo-free six-line ferrihydrite or hematite in wt.% and Fe wt.% Mo(x) stands for the Fe content of Mo-bearing six-line ferrihydrite or hematite, M Mo and M Fe are the molecular weights of Mo and Fe, respectively, and Mo wt.% Mo(x) is the Mo content in wt.% of the Mo-bearing hematite.
The results of the calculations for the synthesized hematite showed that the Fe loss per mole of Mo incorporated (Table 1) increased from ~0.76 at low initial Mo:Fe molar ratios in hem_0.01 to ~2.10 in hem_1.0, where the initial Mo:Fe molar ratio used for synthesizing this sample was 1.0.
Because no Mo-bearing six-line ferrihydrite could be synthesized as a pure single phase and in spite of their large Mo concentrations, no calculations were done for the loss of Fe per Mo atom incorporated.
Powder X-ray Diffraction
Powder XRD patterns of six-line ferrihydrites at various initial Mo:Fe molar ratios (0, 0.01, 0.05, 0.1) (Fig. 1) revealed that akaganéite (β-FeOOH) was present as an impurity in the samples, documented by pXRD peaks at ~17 (Bragg indices 002), ~27 (103), and ~57 (215)°2θ. The intensities of the akaganéite peaks increased with increasing initial Mo:Fe molar ratios. Low intensities and very broad peaks testified to the poor crystallinity of both products and prevented quantitative evaluation of lattice parameters or phase fractions. Only syntheses up to an initial Mo:Fe molar ratio of 0.1 were considered further. The pXRD data showed that higher Mo:Fe ratios resulted in the precipitation of two-line ferrihydrite instead of six-line ferrihydrite.
Powder XRD patterns of the synthesized hematite (Fig. 2) and goethite (Fig. 3) samples displayed only XRD maxima characteristic of each individual phase. The refined lattice parameters and cell volumes of the full-profile fits are given in Table 1. For hematite, in particular, the unit-cell dimensions correlated strongly with the Mo content. The lattice parameter c and the cell volume increased from 13.7615(3) Å and 302.60 Å3 in the absence of Mo to 13.8177(2) Å and 302.84 Å3 at an initial Mo:Fe molar ratio of 1.0. On the contrary, the lattice parameter a decreased from 5.03890(6) Å for pure hematite to 5.03065(4) Å for hem_1.0 (Table 1, Fig. 4a). Expansion of the lattice parameter c agreed with observations for hematite grown in the presence of P5+ (Gálvez et al. Reference Gálvez, Barron and Torrent1999) and As5+ (Bolanz et al. Reference Bolanz, Bläss, Ackermann, Ciobota, Rösch, Tarcea, Popp and Majzlan2013b).
For goethite, the lattice parameter a increased from 9.9551(4) Å in the absence of Mo to 9.9622(3) Å at an initial Mo:Fe molar ratio of 1; parameter c showed negative correlation with increasing Mo content (Fig. 4b). Unit-cell parameter b and the cell volume did not correlate with the Mo concentration.
X-ray Absorption Near-edge Structure (XANES) Spectroscopy
Normalized Mo K-edge XANES profiles of the synthesized phases (Fig. 5) showed almost the same pattern with a distinct pre-edge feature at ~19,996.0 eV, the absorption edge at ~20,006.4 eV, and the crest of the white line at ~20,028.4 eV. Based on the intensity of their pre-edge feature, the samples were divided into two groups. The first group comprised only the spectrum of sample goe_1.0. This spectrum showed a more distinct pre-edge feature than the spectra of hematite and six-line ferrihydrite (Fig. 5) and a slightly less pronounced feature directly before the crest of the white line. In these features, it resembled the spectrum of Na2MoO4⋅2H2O where hexavalent Mo is in tetrahedral coordination with oxygen. The second group includes the spectra of hematite and six-line ferrihydrite. Their pre-edge features resembled more those in MoO3, which consists of layers of MoO6 octahedra in an orthorhombic crystal. None of the samples investigated resembled MoO2 which crystallizes in a monoclinic cell and has a rutile-like structure consisting of distorted Mo(IV) octahedra.
The XANES spectra at the Mo L 3 -edge showed splitting of the white line which is a function of Mo coordination (Fig. 6). The first peak was always more intense than the second one, but the energy differences between the peaks (i.e. peak-splitting values) differed considerably. Hematite spectra showed the greatest values with 3.21–3.27 eV, followed by six-line ferrihydrite with 2.96–3.01 eV. The values for the goethite spectra were considerably smaller than for the other two samples with 2.60–2.82 eV, respectively.
Extended X-ray Absorption Fine Structure (EXAFS) Spectroscopy at the Mo K-edge
Local Mo structure in six-line ferrihydrite. As seen in the pXRD analyses (Fig. 1), six-line ferrihydrite was contaminated by a small amount of akaganéite. All X-ray absorption spectra of these samples were affected by this two-component system. Detailed evaluation of the local structure of Mo in six-line ferrihydrite was complicated. Fitting of the recorded EXAFS spectra sufficed for a general statement concerning the coordination of Mo in this structure and its surrounding short-range order.
Comparison of the EXAFS spectra of pure Mo-bearing akaganéite with the six-line ferrihydrite/akaganéite mixture (Fig. 7) showed that the first shell of both spectra matched well, while the second-shell features of both samples located between 2.2 and 3.4 Å (uncorrected for phase shift) were significantly different. This difference enabled analysis of the position of Mo in the structure of six-line ferrihydrite.
Although the XANES measurements suggested octahedral Mo coordination, fitting of Mo in both octahedral and tetrahedral coordination with the structural model of ferrihydrite from Michel et al. (Reference Michel, Ehm, Antao, Lee, Chupas, Liu, Strongin, Schoonen, Phillips and Parise2007) was implemented. In this model, Fe exists in two different octahedrally coordinated sites (Fe(1) and Fe(2)) and one tetrahedral site Fe(3).
The fitting results for octahedral Mo6+ at the Fe(1) site suggested Mo–O distances of between 1.788(2) and 2.293(3) Å (corrected for phase shift, Table 2), which are consistent with the strongly distorted MoO6 octahedron in sardignaite (Orlandi et al. Reference Orlandi, Pasero and Bigi2010) used by Bolanz et al. (Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017) to fit the first shell of Mo-substituted akaganéite. While the next Fe neighbors coordinated with Mo within the akaganéite structure of Bolanz et al. (Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017) are at 3.311, 3.510, and 3.640 Å (two atoms for each path), the fitting results of the ferrihydrite model showed two additional groups of two Fe atoms at 2.860(2) Å and 3.068(2) Å. To compensate for the strongly distorted MoO6 octahedron, the six remaining Fe atoms split into four Mo-Fe paths with different interatomic distances ranging from 3.204(5) to 3.663(7) Å. All paths of the Fe shell shared one value for the Debye–Waller factor (σ2), but had different Δr values except two groups of Fe atoms at distances of 2.860(2), 3.068(2), and 3.224(2) Å, as well as 3.569(7) and 3.663(7) Å, respectively. The formation of Fe vacancies as described for akaganéite by Bolanz et al. (Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017) was not observed.
ΔE 0 = 0, N var = 14, N ind = 28, R = 0.007, k-range = 3–12, R-range = 1.000–6.000
With Mo–O distances of between 1.764(2) and 2.362(8) Å (corrected for phase shift, Table 3), the fitting results for octahedral Mo6+ at the Fe(2) site were in good agreement with the observed values at the Fe(1) site. In contrast to the fit at the Fe(1) site, no Fe atoms were observed between 2.8 and 3.2 Å. The first two Fe atoms of the second shell were located at 3.209(1) Å. The incorporation of the strongly distorted MoO6 octahedron at the Fe(2) site resulted in a considerable shift of the last seven Mo-Fe paths. These paths were longer than the corresponding calculated Fe–Fe distances. The three furthermost Mo-Fe paths, which were split, showed distances of 3.667(9) and 3.823(1) Å, respectively. The calculated Fe–Fe distances, however, were located at 3.598 Å. All paths of the Fe shell shared the same Debye-Waller factor, which was smaller than that for Mo6+ fit at the Fe(1) site. This could imply that the incorporation of Mo6+ at the Fe(2) site leads to lower structural disorder than for Fe(1).
ΔE 0 = 0, N var = 14, N ind = 28, R = 0.003, k-range = 3–12, R-range = 1.000–6.000
When considering suitable sites for the incorporation of Mo6+ into the structure of six-line ferrihydrite, the coordination number of Mo6+ is essential, either tetrahedral or octahedral (Shannon Reference Shannon1979). Short distances of ~1.7–1.8 Å were considered to be typical for tetrahedral coordination and longer distances of 1.9–2.1 Å indicated octahedral coordination.
For six-line ferrihydrite, fitting with tetrahedral Mo (Fig. 8) resulted in significant distortion of the MoO4 tetrahedron. The distances of three Mo–O paths (1.783(2) Å) were shorter than calculated Fe–O distances (1.953 Å) for a Mo-free model; the fourth O atom at 2.052(9) Å was far more distant than O atoms typical for tetrahedral Mo6+ coordination. The Mo–Fe distances with bridging O atoms with short Mo–O distances were also relatively short, and vice versa (Table 4).
ΔE 0 = 0, N var = 19, N ind = 31, R = 0.027, k-range = 3–13, R-range = 1.000–6.000
In summary, a model with tetrahedral Mo6+ in six-line ferrihydrite led to an exceptional distortion of the tetrahedron and its surroundings. The goodness of fit (RFe(3) = 0.027) was worse than the fit with octahedral Mo6+ (RFe(1) = 0.007 and RFe(2) = 0.003). Hence, fitting with octahedral Mo6+ yielded significantly better results.
Local Mo structure in hematite. The Mo XANES spectra suggested octahedral Mo6+ coordination in the synthetic hematite. The EXAFS spectrum of the sample hem_1.0 was fitted with paths calculated from a hematite supercell (based on the structure of Blake et al. Reference Blake, Hessevick, Zoltai and Finger1966), in which one Mo replaces one Fe (Fig. 9). Due to the significant distortion of the MoO6 octahedron, the six O atoms of the first shell were fitted with three different Mo–O paths and distances of 1.790(3), 2.067(6), and 2.31(2) Å. All three paths were assigned the same amplitude and Debye-Waller factor, but different individual values for Δr. The path calculated for the face-sharing Fe3+ was excluded because such a close position of highly charged cations (Mo6+-Mo6+ or Fe3+-Mo6+) is unlikely. The second shell consisted of Fe atoms forming edge-sharing octahedra with distances of 3.169(8) Å and corner-sharing octahedra at 3.412(9) Å. They shared the same σ2, but had different values for Δr. This approach improved the fit significantly and created an Fe vacancy for charge balance. The deviation of the interatomic distances in the third Fe shell with 5.30(2) Å, compared to 5.44 Å for Mo-free hematite, were significant despite the overall distance. These differences may indicate different degrees of Mo substitution in the hematite structure (Brinza et al. Reference Brinza, Vu, Shaw, Mosselmans and Benning2015).
Local Mo structure in goethite. For goethite, Mo6+ was modeled in both tetrahedral and octahedral coordination (Fig. 9). For Mo in tetrahedral coordination, the EXAFS spectrum of sample goe_1.0 was fitted with the calculated paths of a combined model of a goethite (Gualtieri and Venturelli Reference Gualtieri and Venturelli1999) supercell and the structure of FeMoO4 (Rapposch et al. Reference Rapposch, Anderson and Kostiner1980). The fitting shows that Mo6+ was surrounded by four oxygen atoms at a distance of 1.780(2) Å (Table 5), at the higher end of the range of Mo–O distances in FeMoO4 (1.633–1.779 Å). Fitting of the second shell was performed using the structure of a goethite supercell. The addition of the eight nearest Fe atoms (at Mo–Fe distances of 3.023 to 3.584 Å) and all attempts to refine this model worsened the goodness of fit significantly. Other alternatives were sought, therefore, to explain the second shell features at 2.2–2.5 and 2.7–3.4 Å (uncorrected for phase shift) of the EXAFS spectrum (Fig. 9). The notion of tetrahedral incorporation of Mo6+ in the goethite structure was abandoned and replaced by a model (Table 5) in which tetrahedral MoO4 2– was adsorbed onto goethite. Fitting of Mo6+ as tetrahedrally adsorbed MoO4 2− onto goethite resulted in Mo–Fe distances of 2.90(4) and 3.58(2) Å, both with different values for σ2 (Table 5). The goodness of fit was R = 0.079.
ΔE 0 = 0, N var = 8, N ind = 36, R = 0.079, k-range = 2.5–14.0, R-range = 1.000–6.000
For comparison, Mo6+ was also modeled in octahedral coordination (Fig. 9). Paths created by a Mo-bearing goethite supercell to fit this model resulted in five different Mo–O paths (Table 6) for octahedrally coordinated Mo6+. The Mo–O distances in the first shell ranged from 1.783(5) Å to 2.29(2) Å. The second shell consisted of four edge- and four corner-sharing Fe octahedra with Mo–Fe distances of 2.93(2) and 3.32(3) Å, both for the edge-sharing octahedra, and 3.59(1) Å for the corner-sharing ones, respectively. Removing one of the closest Fe atoms from the fit (located at 2.93(2) Å) reduced the amplitude of the second shell feature at 2.2–2.5 Å and, therefore, improved the fit significantly. The goodness of this fit was with R = 0.036, less than half that of the fit for tetrahedral MoO4 2− adsorbed to goethite (R = 0.079) and, therefore, in better agreement with the measured EXAFS spectrum of goe_1.0.
ΔE 0 = 0, N var = 13, N ind = 36, R = 0.036, k-range = 2.5–14.0, R-range = 1.000–6.000
Discussion
Six-line Ferrihydrite
Powder XRD measurements of the synthesized Mo-bearing six-line ferrihydrites indicated the coprecipitation of chlorine-free akaganéite, as described previously by Linehan et al. (Reference Linehan, Darab, Matson, Chen and Amonette1997), Richmond et al. (Reference Richmond, Hockridge, Loan and Parkinson2004), and Carroll and Richmond (Reference Carroll and Richmond2008). All further compositional and structural investigations were, therefore, affected by this two-component system.
Comparing the Mo L 3 -edge XANES spectra of the six-line ferrihydrite phases with spectra of Mo-bearing akaganéite of Bolanz et al. (Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017) showed that the peak-splitting values of the synthesized samples (2.96 to 3.01 eV) were considerably greater than the 2.7 eV measured for akaganéite (Bolanz et al. Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017). Mo6+ was, therefore, considered to be octahedrally coordinated in six-line ferrihydrite. Nevertheless, because of the generally low peak-splitting values, the possibility of tetrahedrally incorporated Mo6+ was still checked.
Fitting of the EXAFS spectra with Mo in octahedral and tetrahedral coordination suggested that octahedral Mo6+ led to a better fit, applying the structure of Michel et al. (Reference Michel, Ehm, Antao, Lee, Chupas, Liu, Strongin, Schoonen, Phillips and Parise2007). The formation of Fe3+ vacancies as described by Bolanz et al. (Reference Bolanz, Grauer, Cooper, Göttlicher, Steininger, Perry and Küsel2017) for akaganéite or as for hematite (this study) was not observed. Possible ways to balance the excess charge of Mo6+ could be the release of three H+ atoms or the incorporation of three additional hydroxide ions (OH−) in the structure of ferrihydrite, which still is not fully understood. Compared to listed values of Mo–O distances for MoO4 tetrahedra in other phases (Table 7), the incorporation of tetrahedrally coordinated Mo6+ in six-line ferrihydrite resulted in significant distortion of the MoO4 tetrahedron. The longest Mo–O distance, 1.87(2) Å for (NH4)8Mo10O34 (Benchrifa and de Pape Reference Benchrifa and de Pape1990), was still much shorter than 2.050(9) Å, the greatest distance for the MoO4 tetrahedron modeled for six-line ferrihydrite. Distances of 1.9–2.1 Å are characteristic of octahedral coordination of Mo (Shannon Reference Shannon1979) and the incorporation of Mo in tetrahedral coordination would, therefore, result in the destabilization of said tetrahedron and its surroundings.
The authors admit, at this point, that neither the structure of Mo-bearing six-line ferrihydrite as a whole nor the corresponding charge-balance reaction could be resolved. Further experiments and a working method to synthesize pure Mo-bearing six-line ferrihydrite are needed, therefore.
Hematite
Refining of the hematite pXRD spectra resulted in notable changes in the unit-cell parameters a and c together with the increase in unit-cell volume with increasing the initial Mo:Fe molar ratio (Table 1). Similar changes were observed in studies dealing with incorporation of As5+ (Bolanz et al. Reference Bolanz, Bläss, Ackermann, Ciobota, Rösch, Tarcea, Popp and Majzlan2013b) and P5+ (Gálvez et al. Reference Gálvez, Barron and Torrent1999) in hematite. Comparable ionic radii of octahedral Fe3+ (0.65 Å, high-spin state) and Mo6+ (0.59 Å) (Shannon Reference Shannon1979) suggest that Mo6+ could substitute for Fe3+ in a crystal structure. Systematic changes of the cell parameters and the absence of impurities support this assumption.
The XANES spectra at the Mo L 3 -edge are sensitive indicators for the Mo coordination in various structures (Bare et al. Reference Bare, Mitchell, Maj, Edwin, Vrieland and Gland1993; Aritani et al. Reference Aritani, Tanaka, Funabiki and Yoshida1996). The ligand field splitting of the 4d orbitals causes characteristic splitting of the white line (Hedman et al. Reference Hedman, Penner-Hahn, Hodgson, Hodgson, Hedman and Penner-Hahn1984). The split white line is attributed to the electron transition from 2p3/2 to split 4d states, i.e. t2 (dxy, dxz, dyz) and e (dx 2-y 2, dz 2) (Aritani et al. Reference Aritani, Tanaka, Funabiki and Yoshida1996). According to ligand-field theory, for tetrahedral coordination, the intensity of the first peak should be lower than the second (e:t2 = 2:3) and vice versa for octahedral coordination (t2g:eg = 3:2) (Bare et al. Reference Bare, Mitchell, Maj, Edwin, Vrieland and Gland1993). Typical peak-splitting values for tetrahedral compounds lie between 1.8 and 2.4 eV, whereas octahedral values are greater at 3.1–4.5 eV (Bare et al. Reference Bare, Mitchell, Maj, Edwin, Vrieland and Gland1993). Values of 3.21–3.27 eV for the peak splitting of hematite, therefore, could be attributed to the incorporation of a distorted MoO6 octahedron (Aritani et al. Reference Aritani, Tanaka, Funabiki and Yoshida1996).
Replacing one Fe3+ of a hematite supercell by Mo6+ and subsequent fitting of this structure supported the formation of a distorted MoO6 octahedron (for Mo–O distances see Table 8) in the first shell. Several factors could contribute to the local distortion of this polyhedron, such as differences in Fe (0.65 Å) and Mo (0.59 Å) ionic radii (Shannon Reference Shannon1979) and the greater electrostatic charge of Mo6+ compared to Fe3+.
ΔE 0 = 0, N var = 15, N ind = 24, R = 0.027, k-range = 2.5–12, R-range = 1.000–5.000
The Fe vacancy resulting from the exclusion of the closest face-sharing Fe3 + atom in this fit enabled the structure to equilibrate the surplus charge of the incorporated Mo6 + . The calculations for Fe loss suggested a loss of ~2.10 Fe atoms per Mo atom incorporated (Table 1) and are in good agreement, therefore, with this model and support this valence-balancing process:
The decisive step for Mo incorporation by hematite and goethite seems to be the assimilation of Mo into the poorly crystalline iron oxide, ferrihydrite, where coordination of Mo6+ changes from tetrahedral to octahedral. Considering that ferrihydrite is a metastable phase which transforms over time to more stable iron oxides, previously incorporated Mo6+ is then adopted by the newly formed phases. In the next step, ferrihydrite transforms either to hematite by a topotactic solid-state transformation (Cornell and Schwertmann Reference Cornell and Schwertmann2003; Cudennec and Lecerf Reference Cudennec and Lecerf2006), which is probably why larger amounts of Mo6+ are incorporated in the structure of hematite than in goethite, or to goethite by a dissolution-reprecipitation process, where significant amounts of the previously bound Mo6+ are released.
Goethite
Analyzing the refined data of the goethite pXRD spectra, small but noticeable variation of the unit-cell parameters a and c with increasing initial Mo:Fe molar ratio (Fig. 4b) were determined. Observations of similar changes for the incorporation of W into goethite by Kreissl et al. (Reference Kreissl, Bolanz, Göttlicher, Steininger, Tarassov and Markl2016) and the analogous behavior of the hematite samples suggested incorporation of Mo6+ into the crystal structure of the synthesized goethites. No coprecipitated phases were detected in the pXRD spectra. Changes in the chemical composition support this assumption, therefore.
Evaluation of the Mo K- and L 3 -edge spectra of the synthesized goethites led to only a vague impression of the coordination of Mo in goethite. Peak-splitting values of between 2.60 and 2.82 eV in its L 3 -edge XANES spectra were distinctly smaller than for hematite and six-line ferrihydrite. These values are situated between the listed values for tetrahedrally and octahedrally coordinated Mo6+. The intensity ratio (~3:2) of the split peaks in its L 3 -edge XANES spectrum clearly points to the octahedral coordination of Mo6+, however.
In order to resolve the coordination of Mo6+ in the sample, both tetrahedral as well as octahedral coordination for the fitting of the goethite EXAFS spectra were considered. For tetrahedrally coordinated Mo6+, one could assume that the incorporation of the fitted MoO4 2− tetrahedron into the used supercell structure in place of an Fe octahedron would greatly distort the goethite structure. Such distortion would be reflected by significant changes in the unit-cell parameters and especially of the unit-cell volume. Given the small amount of Mo in the solids (0.0274 wt.% at the most), only slight changes in the lattice parameters can be observed. The underlying cause of the observed variations in these parameters cannot be stated with certainty. Not only could they stem from the incorporation of Mo in the goethites, but also from the measured particle sizes, from internal strain of the probes, and from adsorbed ions. The observed changes in the lattice parameters could, at most, serve as a possible hint that no MoO4 2− tetrahedra were incorporated, therefore.
To pursue the approach of tetrahedrally coordinated Mo, the adsorption of tetrahedrally coordinated MoO4 2− on the surface of goethite was assumed. The fitting resulted in a model where MoO4 2− tetrahedra were adsorbed as inner-sphere complexes on the goethite surface. Adsorption would explain the split second-shell feature between 2.2 and 3.4 Å (not corrected for phase shift) in the EXAFS spectrum of goethite (Fig. 9) fairly well. The first peak of this feature between 2.2 and 2.5 Å (not corrected for phase shift) could be explained by an edge-sharing complex via bidentate mono-nuclear linkage. The second, more intense peak (2.7–3.4 Å) could be explained by corner-sharing, inner-sphere surface complexes bound by monodentate mononuclear linkage. The Mo–Fe distances for both complexes, 2.90(4) Å for the bidentate one and 3.58(2) Å for the monodentate complex, were in good agreement with the values suggested by Arai (Reference Arai2010), ~2.81 Å for edge-sharing bidentate Mo-Fe complexes, as well as ~3.59 Å for As-Fe and ~3.63 Å for Cr-Fe linkages postulated by Fendorf et al. (Reference Fendorf, Eick, Gossl and Sparks1997).
As expected, considering the peak-splitting values of the Mo L 3 -edge XANES spectrum of goethite, fitting of the structure where one Fe3+ of the goethite supercell is replaced by Mo6+, results in the formation of a strongly distorted MoO6 octahedron (for Mo–O distances see Table 6) in the first shell. Three positive excess charges created by the substitution of Fe3+ by Mo6+ must be compensated. One possible way to compensate the excess charge would be the release of three H+ atoms. This could lead to a significant reduction of the hydrogen content in goethite, especially for greater Mo concentrations. This deprotonation mechanism was shown by Bolanz et al. (Reference Bolanz, Bläss, Ackermann, Ciobota, Rösch, Tarcea, Popp and Majzlan2013b) to work for low concentrations of highly charged ions such as Sb5+ only. If the concentration of such ions increases, the goethite structure would probably collapse due to missing hydrogen bonds. In dehydroxylation experiments, the break-down of the goethite structure was shown by Gualtieri and Venturelli (Reference Gualtieri and Venturelli1999) to have occurred long before the depletion of the hydroxyl groups. Based on spectroscopic analyses, the formation of Fe3+ vacancies as previously assumed for Mo-bearing hematite seems to be the more likely charge-balancing reaction. This would reduce the four edge-sharing Fe octahedra related to the Mo octahedron incorporated to just three and further improve the goodness of the fit.
As seen by comparing the results of this model with those in which tetrahedral MoO4 2− was adsorbed to goethite, the fit of octahedrally coordinated Mo is in better agreement with the measured EXAFS spectrum of goe_1.0. This model goes much better with the observations of the pXRD measurements. Incorporation of the distorted MoO6 octahedron caused the structural disorder of the synthesized goethites to increase and explains the observed alteration of the cell parameters, therefore (Table 1, Fig. 4b) (Cornell and Schwertmann Reference Cornell and Schwertmann2003). No additional phases containing MoO4 2− were observed in the pXRD spectra.
SUMMARY AND CONCLUSIONS
Iron oxides hematite, six-line ferrihydrite, and goethite have the ability to control the bioavailability of micro-nutrients such as Mo and contaminants in soils (Cornell and Schwertmann Reference Cornell and Schwertmann2003). Their large adsorption capacities for trace elements in solutions (Peacock and Sherman Reference Peacock and Sherman2004; Rout et al. Reference Rout, Mohapatra and Anand2012) but also for element uptake via coprecipitation (Martínez and McBride Reference Martínez and McBride1998, Reference Martínez and McBride2001; Manceau et al. Reference Manceau, Schlegel, Musso, Sole, Gauthier, Petit and Trolard2000) make them well known scavengers for metals and metalloids (e.g. Sieber et al. Reference Sieber, Sanchez, Turner and Somorjai1985; Dardenne et al. Reference Dardenne, Schäfer, Lindqvist-Reis, Denecke, Plaschke, Rothe and Kim2002). Although surface adsorption might play a significant role in short-term retention of such substances, previous studies proposed that structural incorporation might likewise be an important process to contain these substances but for prolonged periods.
Based on chemical and spectroscopic data presented in the current study, the incorporation mechanisms of Mo6+ in the structure of the most common iron oxide phases (i.e. hematite, six-line ferrihydrite, and goethite) were investigated.
Clear indications for changes in the structure of the investigated iron phases are provided by significantly correlated changes for the unit-cell parameters of hematite and goethite with increasing Mo contents. XANES measurements at the Mo L 3 -edge sustain the structural incorporation of Mo6+ in the structure of all three phases. The splitting of the white line suggests the transformation of the originally tetrahedrally coordinated Mo6+ in Na2MoO4⋅2H2O which was used for the synthesis to octahedrally coordinated Mo6+ in the iron oxide phases obtained. Fitting of the EXAFS spectra with structural models of hematite, ferrihydrite, and goethite, with Mo residing in the Fe position, led to reasonable Mo–O distances which matched a strongly distorted MoO6 octahedron incorporated in the structure of all three phases and excluded theoretical models of MoO4 2− tetrahedra adsorbed on the surface of goethite or incorporated in the crystal structures of ferrihydrite and goethite.
These findings support the structural incorporation of Mo6+ in hematite, six-line ferrihydrite, and goethite by heterovalent substitution of one octahedrally coordinated Fe3+ by Mo6+. In the case of hematite and goethite, this substitution led to the formation of an additional Fe3+ vacancy, which compensated the three surplus charges.
Characterizing the incorporation mechanisms of Mo6+ into the crystal structure of six-line ferrihydrite, goethite, and hematite, some of the most common iron oxides, improves understanding of its long-time retention in soils and may represent a long-term control on Mo, as well as other similar metalloids with high coordination numbers in environmental settings. This is relevant to a number of systems, including mining environments, where potentially hazardous metals can be released to the surrounding ecosystem, and agricultural areas, where the bioavailability of MoO4 2− affects the growth of crops and the health of grazing ruminants. In general, the retention of Mo6+ by incorporation into iron oxides makes MoO4 2− less bioavailable for plants and microorganisms (Lang and Kaupenjohann Reference Lang and Kaupenjohann1999), therefore affecting the global nitrogen cycle, as nitrogen fixation and nitrate assimilation can be restricted.
One of the most environmentally consequential polyvalent metals is uranium. Its importance arises as it is the most significant radioactive radionuclide by mass of all environmental contaminants. Because of its long half-life (4.5 Ga), its behavior, and its decay chain, studies dealing with the limitation of uranium migration and its long-term retention are of great environmental importance. In particular, uranium could be incorporated into hematite during ferrihydrite aggregation and the early stages of crystallization, according to Marshall et al. (Reference Marshall, Morris, Law, Livens, Mosselmans, Bots and Shaw2014). The XAS data by those authors indicated the incorporation of U6+ into a distorted octahedral coordination replacing Fe3+. Atomistic simulations by Kerisit et al. (Reference Kerisit, Felmy and Ilton2011, Reference Kerisit, Bylaska, Massey, McBriarty and Ilton2016) support these findings. The replacement of Fe3+ by U6+ in greatly distorted octahedral coordination bears considerable resemblance to the behavior of Mo6+ as seen in the present study. As no suitable charge-compensation mechanism within the constraints of the EXAFS measurements of Marshall et al. (Reference Marshall, Morris, Law, Livens, Mosselmans, Bots and Shaw2014) could be identified, the possibility exists that U6+ follows not only the same incorporation mechanism but also the same charge-compensation mechanism as Mo6+ during its incorporation into hematite. Additional simulations (Kerisit et al. Reference Kerisit, Felmy and Ilton2011, Reference Kerisit, Bylaska, Massey, McBriarty and Ilton2016) support the octahedral incorporation of U6+ in place of Fe3+ in goethite, with the result that one can assume that both Mo6+ and U6+ incorporation into hematite and goethite are subjected to the same mechanisms.
Incorporation of Mo in various iron oxides can be considered, therefore, as a reasonable indication of the uranium scavenger behavior of such minerals; in order to confirm this hypothesis, more comparative analyses are needed. If this assumption proves to be true, molybdenum can be used as a safe surrogate for uranium in some experiments concerning its migration in natural engineered environments as well as for long-term immobilization.
ACKNOWLEDGMENTS
Open Access funding has been enabled and organized by Projekt DEAL, (Germany). The authors thank Dirk Merten and Ines Kamp (Friedrich-Schiller-University Jena) for their support during the ICP-OES analyses. They are grateful to the Karlsruhe Institute of Technology for the beamtime on the SUL-X beamline at the Karlsruhe Research Accelerator (KARA) and acknowledge the Diamond Light Source for time on beamline B18 under proposal SP13606. The Friedrich-Schiller-University, Germany, where most of the experiments were conducted, is also acknowledged. This project was funded by the Deutsche Forschungsgemeinschaft (DFG) through project BO 4534/4-1.
Funding
Funding sources are as stated in the Acknowledgments.
Declarations
Conflict of Interest
The authors declare that they have no conflict of interest.