Incorporation of Mo⁶⁺ in Ferrihydrite, Goethite, and Hematite

Abstract

Among all iron oxides, hematite (α-Fe₂O₃), goethite (α-FeOOH), and ferrihydrite (FeOOH⋅nH₂O) are the most common mineral species. While immobilization of Mo⁶⁺ by surface adsorption on ferric oxides has been studied extensively, the mechanisms of incorporation in their structure have been researched little. The objective of this study was to investigate the relation between Mo content and its structural incorporation in hematite, goethite, and six-line ferrihydrite by a combination of X-ray absorption spectroscopy (XAS), powder X-ray diffraction (pXRD), and inductively-coupled plasma optical emission spectrometry (ICP-OES). Synthesized in the presence of Mo, the hematite, goethite, and six-line ferrihydrite phases incorporated up to 8.52, 0.03, and 17.49 wt. % Mo, respectively. For hematite and goethite, pXRD analyses did not indicate the presence of separate Mo phases. Refined unit-cell parameters correlated with increasing Mo concentration in hematite and goethite. The unit-cell parameters indicated an increase in structural disorder within both phases and, therefore, supported the structural incorporation of Mo in hematite and goethite. Analysis of pXRD measurements of Mo-bearing six-line ferrihydrites revealed small amounts of coprecipitated akaganéite. X-ray absorption near edge structure (XANES) measurements at the Mo L₃-edge indicated a strong distortion of the MoO₆ octahedra in all three phases. Fitting of extended X-ray absorption fine structure (EXAFS) spectra of the Mo K-edge supported the presence of such distorted octahedra in a coordination environment similar to the Fe position in the investigated specimen. Incorporation of Mo⁶⁺ at the Fe³⁺-position for both hematite and goethite resulted in the formation of one Fe vacancy in close proximity to the newly incorporated Mo⁶⁺ and, therefore, charge balance within the hematite and goethite structures.

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. 2001; Mendel and Bittner 2006). Both in the global biogeochemical oxidation-reduction cycles of carbon, nitrogen, and sulfur (Stiefel 2002; Williams and Frausto da Silva 2002) and in the metabolism of living organisms (Bortels 1930) Mo holds a key position in the active sites of a wide range of metalloenzymes (Norrish 1975; Mendel and Bittner 2006).

The most abundant and only bioavailable Mo species in environmental systems is Mo⁶⁺ , in the form of the molybdate anion (MoO₄ ²⁻) (Mendel and Bittner 2006). In contrast to the oceans, where Mo concentrations are up to 10.3 mg/kg (Wang et al. 2011), the soil contains only 1–5 mg of Mo kg^–1 (He et al. 2005). Most soil-bound Mo is considered to not be bioavailable (Žemberyová et al. 2010) because of surface adsorption (Goldberg et al. 1996) or structural incorporation (Richmond et al. 2004; Carroll and Richmond 2008; Das et al. 2016) in the soil minerals. Various studies have emphasized the adsorption mechanisms of Mo to different soil constituents (e.g. Goldberg et al. 1996; Brinza et al. 2008), whereas information about the structural incorporation of Mo in minerals remains scarce (Richmond et al. 2004; Carroll and Richmond 2008; Žemberyová et al. 2010; Wang et al. 2011; Brinza et al. 2015; Das et al. 2016).

Fixation of Mo in German forest soils by ferric oxides was documented by Lang and Kaupenjohann (1999) 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 2000).

Common iron oxide mineral phases in soils are goethite (α-FeOOH), lepidocrocite (γ-FeOOH), feroxyhyte (δ-FeOOH), hematite (α-Fe₂O₃), maghemite (γ-Fe₂O₃), magnetite (Fe₃O₄), and ferrihydrite (FeOOH·nH₂O). Not only do they have large adsorption capacities for trace elements in solution (Peacock and Sherman 2004; Rout et al. 2012) but also for element uptake via coprecipitation (Martínez and McBride 1998, 2001; Manceau et al. 2000). These properties make them well known scavengers for metals (Dardenne et al. 2002; Tarassov et al. 2002; Mullet et al. 2007) and metalloids (e.g. Sieber et al. 1985; Scheinost et al. 2006; Bolanz et al. 2013b) in both natural and contaminated environments (Hochella et al. 2005). Iron oxide phases can, therefore, control the bioavailability of nutrients as well as contaminants such as arsenic (Sherman and Randall 2003; Bolanz et al. 2013a, b) and lead (Vu et al. 2010) in soils.

Previous studies focused predominantly on homovalent and homotopic substitutions, where Fe³⁺ is replaced by other trivalent cations such as Al³⁺, Cr³⁺, or Mn³⁺ (Cornell and Giovanoli 1987; Wells et al. 2001; Mullet et al. 2007). The substitution by aliovalent cations such as Mo⁶⁺, however, depends on intrinsic charge-balancing reactions. Several mechanisms were described or assumed in previous investigations, such as the formation of Fe³⁺ vacancies (Ciobanu et al. 2013), reduction of surrounding Fe³⁺ to Fe²⁺ (Balko and Clarkson 2001), the formation of O^2– vacancies (Khan et al. 2008), and the formation of structural clusters intergrown within the iron oxide host phase (Bolanz et al. 2013b). For Mo⁶⁺, the most abundant Mo species in oxic systems (Barron et al. 2009), 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 Mo⁶⁺ in these oxides; and (3) assessing the effects of Mo⁶⁺ 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 Mo⁶⁺ according to modified procedures described by Schwertmann and Cornell (2000). 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 Na₂MoO₄⋅2H₂O solution adjusted to pH 4 with diluted HNO₃ 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(NO₃)₃⋅9H₂O) was dissolved in 300 mL of deionized water under vigorous stirring. For Mo-bearing two-line ferrihydrite, selected volumes of a Na₂MoO₄⋅2H₂O solution, previously adjusted to pH 4 with diluted HNO₃, 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 (2000). 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 Na₂MoO₄⋅2H₂O 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 K₂HPO₄ and KH₂PO₄ (according to Liang and Zhu 2016) 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 K₂HPO₄ and KH₂PO₄ 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. 2014), and for the fitting, models of six-line ferrihydrite (Michel et al. 2007), hematite (Blake et al. 1966), goethite (Gualtieri and Venturelli 1999), and akaganéite (Post et al. 2003) 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 (R_wp) and the goodness-of-fit (GOF).

Table 1 Chemistry, Fe loss per mole of Mo, and unit-cell parameters of pure and Mo-bearing phases as determined by ICP-OES and full profile fitting with the goethite structure by Gualtieri and Venturelli (1999) and the hematite structure by Blake et al. (1966)

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 (R_wp) 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 (MoO₂, MoO₃, and Na₂MoO₄⋅2H₂O) were collected at the Mo L ₃ -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 2012), 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 ₃ -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 Å⁻¹ 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 Mo⁰ foil was mounted between ionization chambers 2 and 3 at the beamlines. The Mo K- or L ₃ -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 2005). For the shell-by-shell fitting of the k ³ -weighted EXAFS spectra, paths generated from models of six-line ferrihydrite (Michel et al. 2007), a hematite supercell based on the structure of Blake et al. (1966), and a goethite supercell based on the structure of Gualtieri and Venturelli (1999) were used. In these models, Mo replaced one Fe atom. Parameters for the paths were calculated by FEFF6 (de Leon et al. 1991) from the Demeter software package (Ravel and Newville 2005).

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:

(1) $R=\sum {\left({\chi }_i^{\mathrm{data}},-,{\chi }_i^{\mathrm{model}},\left(x\right)\right)}^2/\sum {\left({\chi }_i^{\mathrm{data}}\right)}^2$

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. 2017):

(2) ${\mathrm{Fe}}_{\mathrm{loss}}=\frac{({\mathrm{Fe}}_{\mathrm{wt}.\%}^{\mathrm{Mo}\left(0\right)}-{\mathrm{Fe}}_{\mathrm{wt}.\%}^{\mathrm{Mo}\left(x\right)})\times M_{\mathrm{Mo}}}{{\mathrm{Mo}}_{\mathrm{wt}.\%}^{\mathrm{Mo}\left(x\right)}\times M_{\mathrm{Fe}}}$

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.

Fig. 1 Powder XRD patterns of six-line ferrihydrite synthesized in the absence of Mo^{6 +} and selected samples grown at initial Mo:Fe molar ratios of 0.01, 0.05, and 0.1 display the increase of akaganéite with increasing Mo:Fe ratios. Miller indices for six-line ferrihydrite (black) were taken from Jansen et al. (2002), and for akaganéite (blue) from Post et al. (2003)

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 ų in the absence of Mo to 13.8177(2) Å and 302.84 ų 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 P⁵⁺ (Gálvez et al. 1999) and As⁵⁺ (Bolanz et al. 2013b).

Fig. 2 Powder XRD patterns of hematite synthesized in the absence of Mo^{6 +} and samples grown at initial Mo:Fe molar ratios of 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, and 1.0. Miller indices were taken from Blake et al. (1966)

Fig. 3 Powder XRD patterns of goethite synthesized in the absence of Mo^{6 +} and samples grown at initial Mo:Fe molar ratios of 0.05, 0.1, 0.25, 0.5, 0.75, and 1.0. Miller indices were taken from Gualtieri and Venturelli (1999)

Fig. 4 Comparison of unit-cell parameters of a hematite and b goethite. Error bars of cell parameters are smaller than the symbols

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 Na₂MoO₄⋅2H₂O 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 MoO₃, which consists of layers of MoO₆ octahedra in an orthorhombic crystal. None of the samples investigated resembled MoO₂ which crystallizes in a monoclinic cell and has a rutile-like structure consisting of distorted Mo(IV) octahedra.

Fig. 5 Normalized Mo K-edge XANES spectra of the synthesized phases goe_1.0, hem_1.0, and 6L-Feh_0.1 and reference substances MoO₂, MoO₃, and Na₂MoO₄. A more detailed section showing the pre-edge feature of the spectra is inset. Each spectrum represents the average of at least three spectra

The XANES spectra at the Mo L ₃ -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.

Fig. 6 Normalized XANES spectra at the L ₃ -edge of samples goe_1.0, hem_1.0, and 6L-Feh_0.1 and reference substances MoO₂, MoO₃, and Na₂MoO₄. While MoO₂ shows no distinctive peak-splitting, the splitting and the ratios between the intensities of both peaks for MoO₃ (3:2) and Na₂MoO₄ (2:3) are characteristic of Mo in octahedral and tetrahedral coordination, 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.

Fig. 7 Normalized k³-weighted Mo K-edge spectra (left), magnitude of the Fourier transform (middle), and part of the Fourier transform of the Mo K-edge spectra (right). The thick black lines represent the measured and processed data of 6L-Feh_0.1, the blue lines represent the akaganéite data. All spectra are uncorrected for phase shift

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. (2007) 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 Mo⁶⁺ 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 MoO₆ octahedron in sardignaite (Orlandi et al. 2010) used by Bolanz et al. (2017) 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. (2017) 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 MoO₆ 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 (σ²), 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. (2017) was not observed.

Table 2 Scattering paths for octahedrally coordinated Mo⁶⁺ in six-line ferrihydrite (6L-Feh_0.1) at the Fe(1) site

ΔE ₀ = 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 Mo⁶⁺ 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 MoO₆ 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 Mo⁶⁺ fit at the Fe(1) site. This could imply that the incorporation of Mo⁶⁺ at the Fe(2) site leads to lower structural disorder than for Fe(1).

Table 3 Scattering paths for octahedrally coordinated Mo⁶⁺ in six-line ferrihydrite (6L-Feh_0.1) at the Fe(2) site

ΔE ₀ = 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 Mo⁶⁺ into the structure of six-line ferrihydrite, the coordination number of Mo⁶⁺ is essential, either tetrahedral or octahedral (Shannon 1979). 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 MoO₄ 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 Mo⁶⁺ coordination. The Mo–Fe distances with bridging O atoms with short Mo–O distances were also relatively short, and vice versa (Table 4).

Fig. 8 Normalized k³-weighted Mo K-edge spectra (left), magnitude of the Fourier transform (middle), and part of the Fourier transform of the Mo K-edge spectra (right). The thick black lines represent the measured and processed data of 6L-Feh_0.1, the dashed blue lines represent the fits. All spectra are uncorrected for phase shift

Table 4 Scattering paths for tetrahedrally coordinated Mo⁶⁺ in six-line ferrihydrite (6L-Feh_0.1) at the Fe(3) site

ΔE ₀ = 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 Mo⁶⁺ in six-line ferrihydrite led to an exceptional distortion of the tetrahedron and its surroundings. The goodness of fit (R_{Fe(3) =} 0.027) was worse than the fit with octahedral Mo⁶⁺ (R_Fe(1) = 0.007 and R_Fe(2) = 0.003). Hence, fitting with octahedral Mo⁶⁺ yielded significantly better results.

Local Mo structure in hematite. The Mo XANES spectra suggested octahedral Mo⁶⁺ 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. 1966), in which one Mo replaces one Fe (Fig. 9). Due to the significant distortion of the MoO₆ 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 Fe³⁺ was excluded because such a close position of highly charged cations (Mo⁶⁺-Mo⁶⁺ or Fe³⁺-Mo⁶⁺) 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 σ², 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. 2015).

Fig. 9 Normalized k³-weighted Mo K-edge spectra (left), magnitude of the Fourier transform (middle), and part of the Fourier transform of the Mo K-edge spectra (right). The thick black lines represent the measured and processed data of hem_1.0 and goe_1.0, the dashed blue lines represent the fits. All spectra are uncorrected for phase shift

Local Mo structure in goethite. For goethite, Mo⁶⁺ 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 1999) supercell and the structure of FeMoO₄ (Rapposch et al. 1980). The fitting shows that Mo⁶⁺ 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 FeMoO₄ (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 Mo⁶⁺ in the goethite structure was abandoned and replaced by a model (Table 5) in which tetrahedral MoO₄ ^2– was adsorbed onto goethite. Fitting of Mo⁶⁺ as tetrahedrally adsorbed MoO₄ ²⁻ onto goethite resulted in Mo–Fe distances of 2.90(4) and 3.58(2) Å, both with different values for σ² (Table 5). The goodness of fit was R = 0.079.

Table 5 Scattering paths for tetrahedrally coordinated MoO₄ ²⁻ adsorbed on goethite (goe_1.0)

ΔE ₀ = 0, N _var = 8, N _ind = 36, R = 0.079, k-range = 2.5–14.0, R-range = 1.000–6.000

For comparison, Mo⁶⁺ 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 Mo⁶⁺. 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 MoO₄ ²⁻ adsorbed to goethite (R = 0.079) and, therefore, in better agreement with the measured EXAFS spectrum of goe_1.0.

Table 6. Scattering paths for octahedrally coordinated Mo⁶⁺ in goethite (goe_1.0)

ΔE ₀ = 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. (1997), Richmond et al. (2004), and Carroll and Richmond (2008). All further compositional and structural investigations were, therefore, affected by this two-component system.

Comparing the Mo L ₃ -edge XANES spectra of the six-line ferrihydrite phases with spectra of Mo-bearing akaganéite of Bolanz et al. (2017) 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. 2017). Mo⁶⁺ 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 Mo⁶⁺ was still checked.

Fitting of the EXAFS spectra with Mo in octahedral and tetrahedral coordination suggested that octahedral Mo⁶⁺ led to a better fit, applying the structure of Michel et al. (2007). The formation of Fe³⁺ vacancies as described by Bolanz et al. (2017) for akaganéite or as for hematite (this study) was not observed. Possible ways to balance the excess charge of Mo⁶⁺ 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 MoO₄ tetrahedra in other phases (Table 7), the incorporation of tetrahedrally coordinated Mo⁶⁺ in six-line ferrihydrite resulted in significant distortion of the MoO₄ tetrahedron. The longest Mo–O distance, 1.87(2) Å for (NH₄)₈Mo₁₀O₃₄ (Benchrifa and de Pape 1990), was still much shorter than 2.050(9) Å, the greatest distance for the MoO₄ tetrahedron modeled for six-line ferrihydrite. Distances of 1.9–2.1 Å are characteristic of octahedral coordination of Mo (Shannon 1979) and the incorporation of Mo in tetrahedral coordination would, therefore, result in the destabilization of said tetrahedron and its surroundings.

Table 7. Interatomic distances of Mo⁶⁺ and O²⁻ (d _Mo-O) in selected oxysalts of molybdenum

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 As⁵⁺ (Bolanz et al. 2013b) and P⁵⁺ (Gálvez et al. 1999) in hematite. Comparable ionic radii of octahedral Fe³⁺ (0.65 Å, high-spin state) and Mo⁶⁺ (0.59 Å) (Shannon 1979) suggest that Mo⁶⁺ could substitute for Fe³⁺ 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 ₃ -edge are sensitive indicators for the Mo coordination in various structures (Bare et al. 1993; Aritani et al. 1996). The ligand field splitting of the 4d orbitals causes characteristic splitting of the white line (Hedman et al. 1984). The split white line is attributed to the electron transition from 2p_3/2 to split 4d states, i.e. t₂ (d_xy, d_xz, d_yz) and e (d_x ²_-y ², d_z ²) (Aritani et al. 1996). According to ligand-field theory, for tetrahedral coordination, the intensity of the first peak should be lower than the second (e:t₂ = 2:3) and vice versa for octahedral coordination (t_2g:e_g = 3:2) (Bare et al. 1993). 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. 1993). Values of 3.21–3.27 eV for the peak splitting of hematite, therefore, could be attributed to the incorporation of a distorted MoO₆ octahedron (Aritani et al. 1996).

Replacing one Fe³⁺ of a hematite supercell by Mo⁶⁺ and subsequent fitting of this structure supported the formation of a distorted MoO⁶ 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 1979) and the greater electrostatic charge of Mo⁶⁺ compared to Fe³⁺.

Table 8. Scattering paths for octahedrally coordinated Mo⁶⁺ in hematite (hem_1.0)

ΔE ₀ = 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 Fe^{3 +} atom in this fit enabled the structure to equilibrate the surplus charge of the incorporated Mo^{6 +} . 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:

(3) $2{\mathrm{Fe}}^{3+}\to {\mathrm{Mo}}^{6+}+\mathrm{vacancy}$

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 Mo⁶⁺ changes from tetrahedral to octahedral. Considering that ferrihydrite is a metastable phase which transforms over time to more stable iron oxides, previously incorporated Mo⁶⁺ 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 2003; Cudennec and Lecerf 2006), which is probably why larger amounts of Mo⁶⁺ 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 Mo⁶⁺ 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. (2016) and the analogous behavior of the hematite samples suggested incorporation of Mo⁶⁺ 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 ₃ -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 ₃ -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 Mo⁶⁺. The intensity ratio (~3:2) of the split peaks in its L ₃ -edge XANES spectrum clearly points to the octahedral coordination of Mo⁶⁺, however.

In order to resolve the coordination of Mo⁶⁺ in the sample, both tetrahedral as well as octahedral coordination for the fitting of the goethite EXAFS spectra were considered. For tetrahedrally coordinated Mo⁶⁺, one could assume that the incorporation of the fitted MoO₄ ²⁻ 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 MoO₄ ²⁻ tetrahedra were incorporated, therefore.

To pursue the approach of tetrahedrally coordinated Mo, the adsorption of tetrahedrally coordinated MoO₄ ²⁻ on the surface of goethite was assumed. The fitting resulted in a model where MoO₄ ²⁻ 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 (2010), ~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. (1997).

As expected, considering the peak-splitting values of the Mo L ₃ -edge XANES spectrum of goethite, fitting of the structure where one Fe³⁺ of the goethite supercell is replaced by Mo⁶⁺, results in the formation of a strongly distorted MoO₆ octahedron (for Mo–O distances see Table 6) in the first shell. Three positive excess charges created by the substitution of Fe³⁺ by Mo⁶⁺ 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. (2013b) to work for low concentrations of highly charged ions such as Sb⁵⁺ 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 (1999) to have occurred long before the depletion of the hydroxyl groups. Based on spectroscopic analyses, the formation of Fe³⁺ 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 MoO₄ ²⁻ 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 MoO₆ 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 2003). No additional phases containing MoO₄ ²⁻ 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 2003). Their large adsorption capacities for trace elements in solutions (Peacock and Sherman 2004; Rout et al. 2012) but also for element uptake via coprecipitation (Martínez and McBride 1998, 2001; Manceau et al. 2000) make them well known scavengers for metals and metalloids (e.g. Sieber et al. 1985; Dardenne et al. 2002). 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 Mo⁶⁺ 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 ₃ -edge sustain the structural incorporation of Mo⁶⁺ in the structure of all three phases. The splitting of the white line suggests the transformation of the originally tetrahedrally coordinated Mo⁶⁺ in Na₂MoO₄⋅2H₂O which was used for the synthesis to octahedrally coordinated Mo⁶⁺ 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 MoO₆ octahedron incorporated in the structure of all three phases and excluded theoretical models of MoO₄ ²⁻ tetrahedra adsorbed on the surface of goethite or incorporated in the crystal structures of ferrihydrite and goethite.

These findings support the structural incorporation of Mo⁶⁺ in hematite, six-line ferrihydrite, and goethite by heterovalent substitution of one octahedrally coordinated Fe³⁺ by Mo⁶⁺. In the case of hematite and goethite, this substitution led to the formation of an additional Fe³⁺ vacancy, which compensated the three surplus charges.

Characterizing the incorporation mechanisms of Mo⁶⁺ 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 MoO₄ ²⁻ affects the growth of crops and the health of grazing ruminants. In general, the retention of Mo⁶⁺ by incorporation into iron oxides makes MoO₄ ²⁻ less bioavailable for plants and microorganisms (Lang and Kaupenjohann 1999), 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. (2014). The XAS data by those authors indicated the incorporation of U⁶⁺ into a distorted octahedral coordination replacing Fe³⁺. Atomistic simulations by Kerisit et al. (2011, 2016) support these findings. The replacement of Fe³⁺ by U⁶⁺ in greatly distorted octahedral coordination bears considerable resemblance to the behavior of Mo⁶⁺ as seen in the present study. As no suitable charge-compensation mechanism within the constraints of the EXAFS measurements of Marshall et al. (2014) could be identified, the possibility exists that U⁶⁺ follows not only the same incorporation mechanism but also the same charge-compensation mechanism as Mo⁶⁺ during its incorporation into hematite. Additional simulations (Kerisit et al. 2011, 2016) support the octahedral incorporation of U⁶⁺ in place of Fe³⁺ in goethite, with the result that one can assume that both Mo⁶⁺ and U⁶⁺ 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.