INCORPORATION OF Mo6+ IN FERRIHYDRITE, GOETHITE, AND HEMATITE

Among all iron oxides, hematite (α-Fe2O3), goethite (α-FeOOH), and ferrihydrite (FeOOH⋅nH2O) are the most common mineral species. While immobilization of Mo6+ 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 L3-edge indicated a strong distortion of the MoO6 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 Mo6+ at the Fe3+-position for both hematite and goethite resulted in the formation of one Fe vacancy in close proximity to the newly incorporated Mo6+ 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 6+ , in the form of the molybdate anion (MoO 4 2− ) (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).
Previous studies focused predominantly on homovalent and homotopic substitutions, where Fe 3+ is replaced by other trivalent cations such as Al 3+ , Cr 3+ , or Mn 3+ (Cornell and Giovanoli 1987;Wells et al. 2001;Mullet et al. 2007).The substitution by aliovalent cations such as Mo 6+ , however, depends on intrinsic charge-balancing reactions.Several mechanisms were described or assumed in previous investigations, such as the formation of Fe 3+ vacancies (Ciobanu et al. 2013), reduction of surrounding Fe 3+ to Fe 2+ (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 6+ , 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 6+ in these oxides; and (3) assessing the effects of Mo 6+ incorporation on the structures of iron oxides.

Synthesis of Pure and Mo-bearing Iron Oxides
The investigated iron oxides were synthesized in the absence or presence of Mo 6+ 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 2 MoO 4 ⋅2H 2 O solution adjusted to pH 4 with diluted HNO 3 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 3 ) 3 ⋅9H 2 O) was dissolved in 300 mL of deionized water under vigorous stirring.For Mo-bearing two-line ferrihydrite, selected volumes of a Na 2 MoO 4 ⋅2H 2 O solution, previously adjusted to pH 4 with diluted HNO 3 , 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 2 MoO 4 ⋅2H 2 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 2 HPO 4 and KH 2 PO 4 (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 2 HPO 4 and KH 2 PO 4 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).

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 2 , MoO 3 , and Na 2 MoO 4 ⋅2H 2 O) 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 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 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 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 57.33 ± 0.34 n.v.6L-Feh_0.0145.14 ± 0.34 2.6072 ± 0.0074 6L-Feh_0.0244.25 ± 0.28 3.7133 ± 0.0398 6L-Feh_0.0537.07 ± 0.25 9.5864 ± 0.0997 6L-Feh_0.1 31.97 ± 0.19 17.4853 ± 0.1546 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 Mo 0 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 2005).For the shell-by-shell fitting of the k 3 -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: 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.

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): 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 (2) Mo(x)  wt.% × M Fe 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 P 5+ (Gálvez et al. 1999) and As 5+ (Bolanz et al. 2013b).
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 2 MoO 4 ⋅2H 2 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 3 , which consists of layers of MoO 6 octahedra in an orthorhombic crystal.None of the samples investigated resembled MoO 2 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.27eV, followed by six-line ferrihydrite with 2.96-3.01eV.The values for the goethite spectra were considerably smaller than for the other two samples with 2.60-2.82eV, 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 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 6+ 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 6 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 6+ fit at the Fe(1) site.This could imply that the incorporation of Mo 6+ at the Fe(2) site leads to lower structural disorder than for Fe(1).
When considering suitable sites for the incorporation of Mo 6+ into the structure of six-line ferrihydrite, the coordination number of Mo 6+ 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.4).
In summary, a model with tetrahedral Mo 6+ in sixline 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 6+ (R Fe(1) = 0.007 and R Fe(2) = 0.003).Hence, fitting with octahedral Mo 6+ yielded significantly better results.
Local Mo structure in hematite.The Mo XANES spectra suggested octahedral Mo 6+ 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 6 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 3+ was excluded because such a close position of highly charged cations (Mo 6+ -Mo 6+ or Fe 3+ -Mo 6+ ) 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. 2015).
Local Mo structure in goethite.For goethite, Mo 6+ 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 4 (Rapposch et al. 1980).The fitting shows that Mo 6+ 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 4 (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 6+ in the goethite structure was abandoned and replaced by a model (Table 5) in which tetrahedral MoO 4 2-was adsorbed onto goethite.Fitting of Mo 6+ as tetrahedrally adsorbed MoO 4 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.
For comparison, Mo 6+ was also modeled in octahedral coordination (Fig. 9).Paths created by a  6) for octahedrally coordinated Mo 6+ .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 4 2− adsorbed to goethite (R = 0.079) and, therefore, in better agreement with the measured EXAFS spectrum of goe_1.0.

Six-line Ferrihydrite
Powder XRD measurements of the synthesized Mo-bearing six-line ferrihydrites indicated the coprecipitation of chlorinefree 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.

Single scattering paths
Fitting of the EXAFS spectra with Mo in octahedral and tetrahedral coordination suggested that octahedral Mo 6+ led to a better fit, applying the structure of Michel et al. (2007).The formation of Fe 3+ 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 6+ 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 4 tetrahedra in other phases (Table 7), the incorporation of tetrahedrally coordinated Mo 6+ in six-line ferrihydrite resulted in significant distortion of the MoO 4 tetrahedron.The longest Mo-O distance, 1.87(2) Å for (NH 4 ) 8 Mo 10 O 34 (Benchrifa and de Pape 1990), was still much shorter than 2.050(9) Å, the greatest distance for the MoO 4 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.
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 5+ (Bolanz et al. 2013b) and P 5+ (Gálvez et al. 1999) in hematite.Comparable ionic radii of octahedral Fe 3+ (0.65 Å, high-spin state) and Mo 6+ (0.59 Å) (Shannon 1979) suggest that Mo 6+ could substitute for Fe 3+ 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. 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. ) (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 = 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.27eV for the peak splitting of hematite, therefore, could be attributed to the incorporation of a distorted MoO 6 octahedron (Aritani et al. 1996).
Replacing one Fe 3+ of a hematite supercell by Mo 6+ and subsequent fitting of this structure supported the formation of a distorted MoO 6 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 6+ compared to Fe 3+ .
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: 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 6+ changes from tetrahedral to octahedral.Considering that ferrihydrite is a metastable phase which transforms over time to more stable iron oxides, previously incorporated Mo 6+ is then adopted by the newly formed phases.In the next step, ferrihydrite transforms either to hematite by a topotactic solidstate transformation (Cornell and Schwertmann 2003;Cudennec and Lecerf 2006), which is probably why larger amounts of Mo 6+ 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 6+ 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 6+ 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 Mo 6+ .The intensity ratio (~3:2) of the split peaks in its L 3 -edge XANES spectrum clearly points to the octahedral coordination of Mo 6+ , however.
In order to resolve the coordination of Mo 6+ in the sample, both tetrahedral as well as octahedral coordination for the fitting of the goethite EXAFS spectra were considered.For tetrahedrally coordinated Mo 6+ , one could assume that the incorporation of the fitted MoO 4 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 (3) 2 Fe 3+ → Mo 6+ + vacancy Table 4 Scattering paths for tetrahedrally coordinated Mo 6+ in six-line ferrihydrite (6L-Feh_0.1)at the Fe(3) site ΔE 0 = 0, N var = 19, N ind = 31, R = 0.027, k-range = 3-13, R-range = 1.000-6.000strain 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 4 2− tetrahedra were incorporated, therefore.
To pursue the approach of tetrahedrally coordinated Mo, the adsorption of tetrahedrally coordinated MoO 4 2− on the surface of goethite was assumed.The fitting resulted in a model where MoO 4 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   (1997).
As expected, considering the peak-splitting values of the Mo L 3 -edge XANES spectrum of goethite, fitting of the structure where one Fe 3+ of the goethite supercell is replaced by Mo 6+ , results in the formation of a strongly distorted MoO 6 octahedron (for Mo-O distances see Table 6) in the first shell.Three positive excess charges created by the substitution of Fe 3+ by Mo 6+ 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 5+ 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 3+ 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 4 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 MoO 6 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 4 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 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 andMcBride 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 6+ 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 Mo 6+ in the structure of all three phases.The splitting of the white line suggests the transformation of the originally tetrahedrally  (Hazen et al. 1985) 1.771 CdMoO 4 (Daturi et al. 1997) 1.751 FeMoO 4 (Rapposch et al. 1980) 1.633-1.780FeVMoO 7 (Le Bail et al. 1995) 1.631-1.741Na 2 MoO 4 (Fortes 2015) 1.772 K 2 MoO 4 (Gatehouse and Leverett 1969) 1.742-1.785(NH 4 ) 8 Mo 10 O 34 (Benchrifa and de Pape 1990) 1.701-1.867which matched a strongly distorted MoO 6 octahedron incorporated in the structure of all three phases and excluded theoretical models of MoO 4 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 Mo 6+ in hematite, six-line ferrihydrite, and goethite by heterovalent substitution of one octahedrally coordinated Fe 3+ by Mo 6+ .In the case of hematite and goethite, this substitution led to the formation of an additional Fe 3+ vacancy, which compensated the three surplus charges.
Characterizing the incorporation mechanisms of Mo 6+ 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 4 2− affects the growth of crops and the health of grazing ruminants.In general, the retention of Mo 6+ by incorporation into iron oxides makes MoO 4 2− 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 longterm 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 6+ into a distorted octahedral coordination replacing Fe 3+ .Atomistic simulations by Kerisit et al. (2011Kerisit et al. ( , 2016) ) support these findings.The replacement of Fe 3+ by U 6+ in greatly distorted octahedral coordination bears considerable resemblance to the behavior of Mo 6+ 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 6+ follows not only the same incorporation mechanism but also the same charge-compensation mechanism as Mo 6+ during its incorporation into hematite.Additional simulations (Kerisit et al. 2011(Kerisit et al. , 2016) ) support the octahedral incorporation of U 6+ in place of Fe 3+ in goethite, with the result that one can assume that both Mo 6+ and U 6+ 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.

Fig. 1
Fig.1Powder 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 sixline ferrihydrite (black) were taken fromJansen et al. (2002), and for akaganéite (blue) fromPost et al. (2003)

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

Fig. 5
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 2 , MoO 3 , and Na 2 MoO 4 .A more detailed section showing the pre-edge feature of the spectra is inset.Each spectrum represents the average of at least three spectra

Fig. 7
Fig. 7 Normalized k 3 -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

Fig. 8
Fig. 8 Normalized k 3 -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 e. t 2 (d xy , d xz , d yz ) and e (d

Fig. 9
Fig. 9 Normalized k 3 -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

Table 6 .
Scattering paths for octahedrally coordinated Mo 6+ Fe complexes, as well as ~3.59 Å for As-Fe and ~3.63 Å for Cr-Fe linkages postulated byFendorf et al.

Table 7 .
Interatomic distances of Mo 6+ and O 2− (d Mo-O ) in selected oxysalts of molybdenum