Next Article in Journal
Mineral Montmorillonite Valorization by Developing Ni and Mo–Ni Catalysts for Third-Generation Green Diesel Production
Next Article in Special Issue
Periodic Mesoporous Organosilica Nanoparticles for CO2 Adsorption at Standard Temperature and Pressure
Previous Article in Journal
Stability of Cocaine, Opiates, and Metabolites in Dried Saliva Spots
Previous Article in Special Issue
Brownmillerites CaFeO2.5 and SrFeO2.5 as Catalyst Support for CO Oxidation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure and Population of Complex Ionic Species in FeCl2 Aqueous Solution by X-ray Absorption Spectroscopy

1
Materials Research Laboratory, University of Nova Gorica, Vipavska 13, SI-5001 Nova Gorica, Slovenia
2
Department of Low and Medium Energy Physics, J. Stefan Institute, Jamova 39, SI-1001 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(3), 642; https://doi.org/10.3390/molecules27030642
Submission received: 16 December 2021 / Revised: 13 January 2022 / Accepted: 18 January 2022 / Published: 19 January 2022
(This article belongs to the Special Issue Sol-Gel Functional Materials)

Abstract

:
Technologies for mass production require cheap and abundant materials such as ferrous chloride (FeCl2). The literature survey shows the lack of experimental studies to validate theoretical conclusions related to the population of ionic Fe-species in the aqueous FeCl2 solution. Here, we present an in situ X-ray absorption study of the structure of the ionic species in the FeCl2 aqueous solution at different concentrations (1–4 molL−1) and temperatures (25–80 °C). We found that at low temperature and low FeCl2 concentration, the octahedral first coordination sphere around Fe is occupied by one Cl ion at a distance of 2.33 (±0.02) Å and five water molecules at a distance of 2.095 (±0.005) Å. The structure of the ionic complex gradually changes with an increase in temperature and/or concentration. The apical water molecule is substituted by a chlorine ion to yield a neutral Fe[Cl2(H2O)4]0. The observed substitutional mechanism is facilitated by the presence of the intramolecular hydrogen bonds as well as entropic reasons. The transition from the single charged Fe[Cl(H2O)5]+ to the neutral Fe[Cl2(H2O)4]0 causes a significant drop in the solution conductivity, which well correlates with the existing conductivity models.

Graphical Abstract

1. Introduction

Understanding electrolyte solution properties, such as a population of ionic species, local structure, coordination, hydration, and ion association as a function of concentration and temperature, is important for numerous emerging applications, especially for the growing field of energy storage [1,2] The development of new technologies for mass production relies on cheap and abundant materials [3,4,5]. This supports the use of iron in the form of ferrous chloride for different applications. In the field of energy storage application, it is used as an aqueous electrolyte for iron flow batteries [3,6,7,8,9,10], and Power-to-Solid energy storage technology [11]. Ferrous chloride is used in a manufacturing process of the lithium iron orthosilicate (Li2FeSiO4) [12] and FeSe2 cathode [13], FeCl2-graphite composite anode used for lithium-ion batteries [14], and Na-FeCl2 batteries [15]. Additionally, the use of ferrous chloride is reported for wearable thermoelectric cells [16,17], in environmental applications for arsenic stabilization [18], as well in spintronics for the growth of single-layer magnetic material [19,20].
The existing commercial flow battery technologies have limited energy density due to the solubility limits of the electroactive species [21]. The energy losses of the electrochemical cell are caused by electrodes resistance, electrolyte resistance, and the steric impediment of the ion exchange membranes [22,23,24,25]. For polymeric separators, such as ion exchange membranes used for electrodialysis, reverse osmosis, flow batteries, etc., there is required a fundamental understanding of the influence of ion size and ion charge on salt transport properties [26]. Consequently, the optimization of electrodes, electrolytes, membranes, and other crucial parts of the electrochemical devices is recognized as the important key requirement for the future deployment of new energy storage systems [27].
Molecular dynamic simulations [28,29,30] and experimental analyses [31,32,33] showed that the first hydrated shell of Fe2+ in a pure aqueous solution is characterized by a regular octahedron with a Fe-O distance of about 2.10 Å. A density functional theory analysis of the optimized structures predicted a Fe(H2O)62+ octahedral coordination with inequivalent Fe-O distances suggesting the Jahn–Teller distortion [34]. Opposite to that, four studies ruled out the occurrence of the Jahn–Teller distortion in an Fe2+ aqueous solution [31,33,35,36]. The first one was done by combined ab initio quantum mechanical/molecular mechanical molecular dynamics simulations [35], while the other studies involved both the molecular dynamics simulations and X-ray absorption structure method [31,33,36]. In the second hydration shell, 12.4 water molecules were found [35].
In the chloride solutions, the chloride ions participate in the structure of Fe ionic species. Lee et al. [37] calculated, based on estimated equilibrium constants, that in the FeCl2 solution, the predominant ionic species are [FeCl]+ ions. In the entire concentration range, their share does not drop below 90%. According to them, the other two significantly represented species are Fe2+ and [FeCl2]0. With the increase in FeCl2 concentration, the share of Fe2+ decreases on account of an increase in the neutral [FeCl2]0. This was confirmed for high temperatures (>200 °C), high pressures (>300 bar), and high Cl concentrations [38,39]. The calculated ion association constant for the association of Fe2+ and Cl to FeCl+ is −0.88(5) [40].
Zhao et al. [41] calculated the stability constant of different Fe-chloride complexes from UV and near IR spectra. The results are not consistent with that of Lee et al. [37]. Their calculations showed that Fe2+ ions are predominantly present in highly concentrated FeCl2 solutions (up to about 3 molL−1). Only then did the other complex Fe-chlorine species start to prevail. They also claim that with an increase in concentration and temperature, the Fe-complexes undergo a configurational transformation from octahedral to tetrahedral coordination. The transformation is driven by an increase in entropy change due to a replacement of water molecules by the chloride ions [41].
The literature survey shows that the experimental confirmation of the theoretical studies, related to the population of different ionic Fe-species in ferrous chloride (FeCl2) aqueous solution (hereafter FeCl2 (aq)), have not been performed. This is critical because the theoretical studies have not offered to same conclusions and, consequently, we cannot be sure what real ionic structures in this solution are. For that reason, we focused on the experimental investigation of ionic species present in the FeCl2 (aq) at different concentrations and temperatures. X-ray absorption spectroscopy (XAS) is the method of choice to study the valence state and local structure of the ionic species in solutions and disordered materials [42,43,44,45,46]. We used in situ Fe K-edge XANES (X-ray Absorption Near Edge Structure) and EXAFS (Extended X-ray Absorption Fine Structure) analysis to monitor the Fe valence state [32,45,47] and local structure around the Fe cations in the aqueous solutions [31,32,33]. The target concentration range was from 1 to 4 molL−1 and the temperature range was from 25 °C to 80 °C. Using the conductivity model developed by Zhang et al. [48], we showed that the conclusions of the XAS studies well correlate with the conductivity measurements.

2. Results

2.1. Fe K-Edge XANES Results

The Fe K-edge XANES analysis is used to determine the valence state and local symmetry of Fe cations in the 1 molL−1 FeCl2 (aq) in a temperature range from RT to 80 °C and in the FeCl2 (aq) with higher concentrations up to 4 molL−1. Normalized Fe K-edge XANES spectra are shown in Figure 1, together with the spectra of selected Fe reference compounds with known Fe valence states and local symmetry of Fe atom neighborhood (1 molL−1 FeCl3 (aq) as a reference for Fe3+ and crystalline FeSO4·7H2O as a reference for octahedrally coordinated Fe2+ cations).
Different local environments of the Fe cations result in different Fe K-edge profiles and pre-edge lines in the XANES spectra. The energy position of the Fe absorption edge and the pre-edge features are correlated with the valence state of the absorbing atom in the sample. With increasing oxidation, each absorption feature in the XANES spectrum is shifted to higher energies. The Fe K-edge shift of about 4.5 eV was found between the spectra of the two- and three-valent iron compounds [32,45,47,49].
The Fe K-edge XANES spectra of all FeCl2 (aq) (Figure 1) exhibit the same edge profile, characteristic for octahedrally coordinated Fe2+ cations, as in the case of reference FeSO4·7H2O sample or [Fe(H2O)6]2+ complexes in water solution [32]. The energy position of the Fe K-edge in all FeCl2 (aq) solutions is identical, coinciding with the edge position of the reference FeSO4·7H2O compound. Therefore, the XANES results show that all Fe cations in the FeCl2 (aq) at all investigated concentrations up to 4 molL−1 are in divalent form with octahedral coordination in the nearest coordination shell. There are also no changes of Fe oxidation state in the solution during heating up to 80 °C. The results give no indications for the presence of tetrahedral Fe(II)-chloride complex [FeCl4] with a characteristic pre-edge peak in XANES spectrum at 7112 eV [33,49], which is in agreement with previous findings that indicated the formation of tetrahedral FeCl4 species in the aqueous solution at RT only at high Cl concentrations above 8 molL−1 or at high temperatures above 150 °C at 6 molL−1 [33,49].

2.2. Fe K-Edge EXAFS Results

The Fe K-edge EXAFS analysis was used to directly probe the local structure around Fe cations in the selected temperature and concentration range. Fourier transforms (FT) of the k3-weighted Fe K-edge EXAFS spectra (Figure 2, Figure 3 and Figure 4) exhibit a strong peak at about 2 Å, which can be ascribed to photoelectron backscattering on neighbor atoms in the first Fe coordination shell, while the signal from second and more distant Fe hydration shells in the solution is negligible. The weak signal in the R range between 3 and 5 Å can be attributed to the multiple scatterings of the photoelectron within the nearest coordination shell of neighbors. Qualitative comparisons of the Fe FT EXAFS spectra in the R range between 1 and 2.5 Å show that the average Fe neighborhood in the solutions is very similar but not the same. A progressive decrease of the nearest neighbor peak is observed with an increase in temperature and concentration of the FeCl2 (aq).
Structural parameters of the average local Fe neighborhood (type and the average number of neighbors, the radii, and Debye–Waller factor of neighbor shells) are quantitatively resolved from the Fe EXAFS spectra by comparing the measured EXAFS signal with model signal, constructed ab initio with the FEFF6 program code [50], in which the photoelectron scattering paths are calculated ab initio from a tentative spatial distribution of Fe neighbor atoms. The atomic species of neighbors are identified in the fit by their specific scattering factor and phase shift. The FEFF model of Fe cations in the solution was based on previous structural studies of FeCl2 (aq), which indicate that different Fe complexes ([FeClx(H2O)6−x]2−x for x = 0–2) with octahedral geometry may be expected in FeCl2 (aq) at concentrations up to 4 molL−1 and temperatures below 100 °C [31,32,33,49]. The model structures of Fe-ionic species in FeCl2 (aq) are schematically presented in Figure 5.
The FEFF model of the nearest Fe coordination shell included two single scattering paths, one for Cl neighbors at 2.3 Å and the second for oxygen neighbors at 2.1 Å, and all multiple (triangular and linear) scattering paths with a total length up to 5 Å. Three variable parameters for each type of neighbor in the first coordination sphere were varied: the coordination number (N), the distance to the neighbor atoms (R), and the Debye–Waller factors (σ2). The amplitude reduction factor (S02) was fixed at the value of 0.85. A shift of energy origin ΔEo, common to all scattering paths, was also varied. The total number of O and Cl neighbor atoms in the nearest Fe coordination shell was constrained to six and arranged in a distorted octahedron, as suggested by XANES analysis. The parameters of multiple scattering paths within the distorted octahedra of nearest neighbors were constrained with the values of structural parameters obtained by modeling the single scattering paths on the same neighbors.
The spectra were fitted in the k range from 3.4 Å−1 to 13 Å−1. A very good fit is obtained for all spectra in the R range of 1 Å up to 5 Å. (Figure 2, Figure 3 and Figure 4). A complete list of the best-fit parameters is listed in Table 1. The parameters of multiple scattering paths are listed in Table S1 in the Supplementary Materials.
In all cases, oxygen and chlorine atoms were identified in the nearest Fe coordination shell, with the Fe–O distance of 2.095 ± 0.005 Å and Fe–Cl distance of 2.33 ± 0.02 Å, in agreement with previous findings [31,32,33]. In the 1 molL−1 FeCl2 (aq) at RT, we found 5.0 oxygen and 1.0 chlorine neighbors on average, clearly indicating that the dominant species are [FeCl(H2O)5]+ complexes. A systematic increase in the number of Cl neighbors and a decrease in the number of O neighbors is observed at higher temperatures (Table 1).
The Debye–Waller factors of Fe-Cl coordination are significantly larger, indicating larger static disorder for Cl neighbors present in the first hydration shell of the octahedral Fe2+ ion. A smaller increase of the Debye–Waller factor of both species is detected with increased temperature as expected due to larger thermal disorder in the first coordination shell at higher temperatures. The results indicate that at higher temperatures the [FeCl2(H2O)4]0 complexes predominate over [FeCl(H2O)5]+.
Similar results are obtained also for solutions with higher concentrations at RT (Table 1). A systematic increase in the number of Cl neighbors and a decrease in the number of O neighbors is observed at higher concentrations, indicating that for the 3 molL−1 and 4 molL−1 solutions [FeCl2(H2O)4]0 are dominant complexes already at RT.
The EXAFS results indicate similar local structure around Fe2+ cations in FeCl2 (aq) as observed for Zn2+ cations in zinc(II) chloride in aqueous solution [44]. For a concentration of 1 molL−1, the existence of a six-fold Zn coordination with about 4.8 water molecules and 1.2 Cl ions were determined. With an increase in the Cl concentration, the water molecules bound to the Zn2+ ion are gradually exchanged with the chloride ions.

3. Discussion

The increase in the population of Cl ions in the first coordination sphere of the Fe-ionic species with temperature and concentration can be rationalized based on the obtained Fe K-edge EXAFS results (Table 2). At the low temperature and low FeCl2 concentration, the Fe2+ ions in the FeCl2 (aq) are octahedrally coordinated with one Cl and five H2O, which can be written as [FeCl(H2O)5]+. The polar water molecules are oriented with their negative oxygen side towards Fe2+ and all five Fe-O distances are the same (within error bars)—2.095 ± 0.005 Å (Figure 5). The Cl ion is bonded to Fe at a distance of 2.33 ± 0.02 Å, which is a characteristic Fe-Cl bond distance for chloride ions in the first hydration shell of the octahedral Fe2+ ion [33]. The formation of intramolecular hydrogen bonds between the hydrogen atoms of the basal water molecules and the highly electronegative apical chlorine can be expected in such a tightly bound hydration complex. This restrains the orientation of the water molecules within the complex ionic structure and does not allow for a significant dynamic structure disorder as noticed by Roscioni et al. [46] for the hydration shell around Zn2+. In addition, the electronegativity of chlorine shifts the electron density away from the apical oxygen, which results in the weakening and destabilizing of its electrostatic bond with Fe. As a consequence, with an increase in temperature, the breaking of the bond becomes more probable, and the substitution of the apical water molecule by the Cl ion, which due to its charge can form a stronger electrostatic bond, possible (Figure 5).
The values of the Fe-O Debye–Waller factors of about 0.007 Å2 are typical for the hydration shell of divalent ions in solutions [31,33], where the metal cation forms a quite stable hydration complex and fluctuations of the first-shell water molecules are due to the thermal motion only [46]. The Debye–Waller factors of Fe-Cl coordination are significantly larger, indicating larger static disorder for the Cl neighbors present in the first hydration shell of the octahedral Fe2+ ionic complex.
The small increase in the Debye–Waller factors with temperature and concentration indicates that the Fe2+ cation forms a quite stable and tightly bound hydration complex. This again supports the presence of the intramolecular hydrogen bonds between the water molecules and two apical chlorine ions, which restrain the thermal oscillations in the high-temperature Fe[Cl2(H2O)4]0 even more than in [FeCl(H2O)5]+.
Apparently, the described substitution mechanism decreases the entropy of the system because the initially free Cl ions get bound to Fe2+ to form the neutral Fe[Cl2(H2O)4]0 complex. This is why the system at low temperature and low chloride concentration prefers one and not two Cl in the first coordination shell. However, with the increase in temperature, the increased thermal oscillation of the apical water molecule further destabilizes its bond with Fe. The decrease in the substitutional enthalpy overcomes the entropy drop (-T∆S) and makes the change of free energy for this process negative, thus enabling the substitution.
The observed change in the Fe coordination shell with the increase in the FeCl2 concentration can also be explained considering the thermodynamic state of the system. Here, we have to emphasize that we are working with highly concentrated solutions up to 4 molL−1. In such solution, the nominal species ratio Fe2+:Cl:H2O is about 1:2:14. Taking into account the structure of the monovalent ion complex, the ratio [FeCl(H2O)5]+:Cl:H2O becomes 1:1:9. This causes the shortage of the water molecules available for hydration of the Cl ions. The studies show that a typical first hydration sphere of Cl ion consists of about 6 water molecules [51,52]. Not accounting for the outer hydration spheres, this makes the species ratio [FeCl(H2O)5]+:[Cl(H2O)6]:H2O to only be 1:1:3. So, the majority of the water molecules are used up for the first hydration sphere of the present ionic species, and not enough are available for outer hydration layers. This is thermodynamically very unfavorable. The free energy of such a system highly increases and destabilizes the system. As a consequence, the formation of the neutral Fe[Cl2(H2O)4]0 species becomes thermodynamically favorable. For each transformation of [FeCl(H2O)5]+ into Fe[Cl2(H2O)4]0, one water molecule is released in addition to the water molecules that are released from the Cl hydration sphere.
With the increase in concentration, the association of the charged Cl and [FeCl(H2O)5]+ into neutral Fe[Cl2(H2O)4]0 must influence the conductivity of the solutions due to the decrease in the concentration of the charged ionic species. Figure 6 shows that the conductivity in the dilute regime increases with the concentration. This is expected as the ionic mobility and the population of species are not affected much by the increase in the concentration whereas the number of charge carriers (i.e., charged ionic species) significantly increases. With a further increase in the concentration, the ionic mobility is getting hindered due to ionic interactions, and the population of the ionic species changes in favor of the neutral species. This makes the initial increase in conductivity to slow down and reaches the maximum value at around 2.2 mol·kg−1. Eventually, at even higher concentrations, the conductivity decreases sharply as a result of an ever lower concentration of charged ionic species and a hindered ionic mobility due to a lack of free water molecules. The decrease is much sharper than in the case of, e.g., divalent CaCl2. Up to the concentration of 4.5 mol·kg−1, the conductivity of CaCl2 (aq) decreases from its maximum for only about 10%, while it decreases by 35% for FeCl2 (aq).
The observed association of the ionic species that causes the nonlinear relationship between a concentration of the electrolyte and a concentration of free ions can also be deduced from the model developed by Zhang et al. [48]:
κ = P 1 T + P 2 m n exp P 3 m T P 4
where κ, T, and m are conductivity, temperature, and molality. P1, P2, P3, P4, and n are constants independent of concentration and temperature.
The measurements from Figure 6 were fitted by Equation (1) to extract the model parameters. If the exponential value n is equal to 1, the relationship is linear while the deviation from 1 reflects nonlinearity. For NaCl (aq) and CaCl2 (aq), the n value is 1 ± 0.005, which indicates the negligible association of ions. For FeCl2 (aq) the n value is 0.94, which supports the conclusions of XAS studies on the ionic association.

4. Materials and Methods

4.1. Sample Preparation

The FeCl2 (aq) solutions with concentrations from 1 molL−1 to 4 molL−1 were prepared from FeCl2·4H2O (Acros Organics 70 Mesh, <212 µm particle size, 99+% purity) and followed by a reduction procedure to eliminate traces of Fe (III). The reduction procedure consisted of adding 5 mL of HCl per 1 L of solution and Fe powder in excess. Typically, a clear green-colored solution with pH from 3.5 to 4 is obtained in 12–36 h at room temperature. When the target pH is reached the excess of iron was removed. The sample solutions were stored under reductive conditions in air-tight transparent borosilicate glass bottles under a nitrogen gas atmosphere at room temperature during transportation. The correct FeCl2 (aq) concentration was recalculated to account for the reduction procedure.
The measurements of pH, conductivity, and temperature were performed using Phoenix EC45-Multi probes. The conductivity measuring cell, Sentek VPT80C10, was calibrated using Hamilton ZDL 300 mL standard solution with a conductivity of 100 mScm−1 at 25 °C. De-ionized water (Grade 2) was used for all the prepared solutions.

4.2. In Situ XAS of FeCl2 (aq) As a Function of Concentration and Temperature

In situ Fe K-edge X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) measurements of 1 molL−1 FeCl2 aqueous solution in a temperature range from RT to 80 °C and the FeCl2 aqueous solutions at higher concentrations (2, 3, and 4 molL−1) at RT were performed at P65 beamline of PETRA III, at DESY in Hamburg. In addition, 1 molL−1 FeCl3 aqueous solution and crystalline FeSO4·7H2O were measured at room temperature for comparison. The XAS spectra were measured in transmission detection mode. We used our on-purpose designed experimental setup which consists of a PMMA (Polymethyl methacrylate) liquid absorption cell (Figure 7) mounted inside a remotely controlled heating system. The tubular oven allowed the in-situ temperature-controlled measurements at well-defined 10 °C temperature steps. The PMMA cell windows were machined down to 0.2 mm thickness to minimize the absorption. The optimal absorption thickness of the sample solution was adjusted over EPDM rubber spacers thickness providing a solution layer from 0.1 to 0.3 mm. The solution was inserted in the liquid absorption cell by a syringe, in a glove box under a protective atmosphere to prevent oxidation of the Fe and sealed. The sealed cell was inserted in the tubular oven, which was mounted on the beamline between the first and second ionization cells. The oven was closed with aluminum foil (20 µm) windows and filled with He protective atmosphere. The absorption cell was equipped with a thermocouple to read the temperature of the liquid during the in situ XAS experiments.
A Si(111) double crystal monochromator was used with an energy resolution of about 1 eV at 7 keV. The beam size on the sample was 1.5 mm horizontally and 0.2 mm vertically. The absorption spectra were measured in the energy region from −150 eV to +1000 eV relative to the Fe K-edge (7112 eV) in continuous fast (3 min) scans and re-binned to equidistant energy steps of 0.25 eV in the XANES region and equidistant k steps of 0.05 Å−1 in the EXAFS region. The exact energy calibration was established with simultaneous absorption measurement on a 5-micron thick Fe foil placed between the second and third ionization detector. The Fe K-edge XAFS spectra were measured first at RT and then during heating at well-defined 10 °C temperature steps up to 80 °C. At each intermediate temperature state, we performed three to five repetitions of the same scan to check the reproducibility and to improve the signal-to-noise ratio. We measured XAS scans on different spots of the absorption cell window to reduce the eventual effects of radiation damage.
The quantitative analysis of XANES and EXAFS spectra was performed with the Demeter (IFEFFIT) program package [42] in combination with the FEFF6 program code for ab initio calculation of photoelectron scattering paths [50].

Supplementary Materials

The following supporting information can be downloaded online, Table S1: Parameters of photoelectron multiple scattering paths on the neighbors in the nearest coordination shells around Fe cations in the FeCl2 samples.

Author Contributions

Conceptualization, M.V. and I.A.; Methodology, U.L. and I.A.; Formal Analysis, I.A.; Investigation, U.L. and M.V.; Writing, U.L., M.V. and I.A.; Funding Acquisition, M.V. and I.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency (core programs P1-0112 and P2-0412) and by the project CALIPSOplus under Grant Agreement 730872 from the EU Framework Programme for Research and Innovation HORIZON 2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We acknowledge access to the synchrotron radiation facilities at beamline P65 of PETRA III, DESY, Hamburg (project I-20200036 EC). We would like to thank Edmund Welter of PETRA III for expert advice on beamline operation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. U.S. Department of Energy. Energy Storage Grand Challenge Energy Storage Market Report 2020; U.S. Department of Energy: Washington, DC, USA, 2020; Technical; p. 65.
  2. Comello, S.; Reichelstein, S. The emergence of cost effective battery storage. Nat. Commun. 2019, 10, 2038. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bock, R.; Kleinsteinberg, B.; Selnes-Volseth, B.; Burheim, O.S. A novel iron chloride red-ox concentration flow cell battery (Icfb) concept; power and electrode optimization. Energies 2021, 14, 1109. [Google Scholar] [CrossRef]
  4. Yu, X.; Manthiram, A. Sustainable Battery Materials for Next-Generation Electrical Energy Storage. Adv. Energy Sustain. Res. 2021, 2, 2000102. [Google Scholar] [CrossRef]
  5. Jayathilake, B.S. Electrochemical Pathways for Sustainable Energy Storage and Energy Conversion. Ph.D. Thesis, University of Southern California, Los Angeles, CA, USA, 2019. [Google Scholar]
  6. Qi, Z.; Koenig, G.M. Review Article: Flow battery systems with solid electroactive materials. J. Vac. Sci. Technol. B Nanotechnol. Microelectron. Mater. Process. Meas. Phenom. 2017, 35, 040801. [Google Scholar] [CrossRef] [Green Version]
  7. Gong, K.; Xu, F.; Grunewald, J.B.; Ma, X.; Zhao, Y.; Gu, S.; Yan, Y. All-Soluble All-Iron Aqueous Redox-Flow Battery. ACS Energy Lett. 2016, 1, 89–93. [Google Scholar] [CrossRef]
  8. Yensen, N.; Allen, P.B. Open source all-iron battery for renewable energy storage. HardwareX 2019, 6, e00072. [Google Scholar] [CrossRef]
  9. Petek, T.J.; Hoyt, N.C.; Savinell, R.F.; Wainright, J.S. Slurry electrodes for iron plating in an all-iron flow battery. J. Power Sources 2015, 294, 620–626. [Google Scholar] [CrossRef]
  10. Alotto, P.; Guarnieri, M.; Moro, F. Redox flow batteries for the storage of renewable energy: A review. Renew. Sustain. Energy Rev. 2014, 29, 325–335. [Google Scholar] [CrossRef]
  11. Valant, M. Procedure for Electric Energy Storage in Solid Matter. US Patent 10,906,805 B2, 2 February 2021. [Google Scholar]
  12. Kucuk, K.; Aryal, S.; Moazzen, E.; Timofeeva, E.V.; Segre, C.U. In situ XAS study of the local structure of the nano-Li2FeSiO4/C cathode. J. Phys. Energy 2021, 3, 034015. [Google Scholar] [CrossRef]
  13. Wang, L.; Guo, C.; Zhu, Y.; Zhou, J.; Fan, L.; Qian, Y. A FeCl2-graphite sandwich composite with Cl doping in graphite layers: A new anode material for high-performance Li-ion batteries. Nanoscale 2014, 6, 14174–14179. [Google Scholar] [CrossRef]
  14. Zhang, C.; Zhang, L.; Li, N.; Zhang, X. Studies of FeSe2 cathode materials for Mg-Li hybrid batteries. Energies 2020, 13, 4375. [Google Scholar] [CrossRef]
  15. Zhan, X.; Bowden, M.E.; Lu, X.; Bonnett, J.F.; Lemmon, T.; Reed, D.M.; Sprenkle, V.L.; Li, G. A Low-Cost Durable Na-FeCl2 Battery with Ultrahigh Rate Capability. Adv. Energy Mater. 2020, 10, 1903472. [Google Scholar] [CrossRef]
  16. Yang, P.; Liu, K.; Chen, Q.; Mo, X.; Zhou, Y.; Li, S.; Feng, G.; Zhou, J. Wearable Thermocells Based on Gel Electrolytes for the Utilization of Body Heat. Angew. Chem.-Int. Ed. 2016, 55, 12050–12053. [Google Scholar] [CrossRef] [PubMed]
  17. Bai, C.; Wang, Z.; Yang, S.; Cui, X.; Li, X.; Yin, Y.; Zhang, M.; Wang, T.; Sang, S.; Zhang, W.; et al. Wearable Electronics Based on the Gel Thermogalvanic Electrolyte for Self-Powered Human Health Monitoring. ACS Appl. Mater. Interfaces 2021, 13, 37316–37322. [Google Scholar] [CrossRef] [PubMed]
  18. Lin, Y.; Wu, B.; Ning, P.; Qu, G.; Li, J.; Wang, X.; Xie, R. Stabilization of arsenic in waste slag using FeCl2 or FeCl3 stabilizer. RSC Adv. 2017, 7, 54956–54963. [Google Scholar] [CrossRef] [Green Version]
  19. Zhou, X.; Brzostowski, B.; Durajski, A.; Liu, M.; Xiang, J.; Jiang, T.; Wang, Z.; Chen, S.; Li, P.; Zhong, Z.; et al. Atomically Thin 1T-FeCl 2 Grown by Molecular-Beam Epitaxy. J. Phys. Chem. C 2020, 124, 9416–9423. [Google Scholar] [CrossRef]
  20. Feng, Y.; Wu, X.; Hu, L.; Gao, G. FeCl2/MoS2/FeCl2van der Waals junction for spintronic applications. J. Mater. Chem. C 2020, 8, 14353–14359. [Google Scholar] [CrossRef]
  21. Zhong, F.; Yang, M.; Ding, M.; Jia, C. Organic Electroactive Molecule-Based Electrolytes for Redox Flow Batteries: Status and Challenges of Molecular Design. Front. Chem. 2020, 8, 451. [Google Scholar] [CrossRef] [PubMed]
  22. Kamcev, J.; Sujanani, R.; Jang, E.-S.; Yan, N.; Moe, N.; Paul, D.R.; Freeman, B.D. Salt concentration dependence of ionic conductivity in ion exchange membranes. J. Memb. Sci. 2018, 547, 123–133. [Google Scholar] [CrossRef]
  23. Díaz, J.C.; Kamcev, J. Ionic conductivity of ion-exchange membranes: Measurement techniques and salt concentration dependence. J. Memb. Sci. 2021, 618, 118718. [Google Scholar] [CrossRef]
  24. Zhu, S.; Kingsbury, R.S.; Call, D.F.; Coronell, O. Impact of solution composition on the resistance of ion exchange membranes. J. Memb. Sci. 2018, 554, 39–47. [Google Scholar] [CrossRef]
  25. Bamgbopa, M.O.; Shao-Horn, Y.; Hashaikeh, R.; Almheiri, S. Cyclable membraneless redox flow batteries based on immiscible liquid electrolytes: Demonstration with all-iron redox chemistry. Electrochim. Acta 2018, 267, 41–50. [Google Scholar] [CrossRef]
  26. Geise, G.M.; Paul, D.R.; Freeman, B.D. Progress in Polymer Science Fundamental water and salt transport properties of polymeric materials. Prog. Polym. Sci. 2014, 39, 1–42. [Google Scholar] [CrossRef]
  27. Park, M.; Ryu, J.; Wang, W.; Cho, J. Material design and engineering of next-generation flow-battery technologies. Nat. Rev. Mater. 2016, 2, 16080. [Google Scholar] [CrossRef]
  28. Kuharski, R.A.; Bader, J.S.; Chandler, D.; Sprik, M.; Klein, M.L.; Impey, R.W. Molecular model for aqueous ferrous–ferric electron transfer. J. Chem. Phys. 1988, 89, 3248–3257. [Google Scholar] [CrossRef]
  29. Bader, J.S.; Chandler, D. Computer simulation study of the mean forces between ferrous and ferric ions in water. J. Phys. Chem. 1992, 96, 6423–6427. [Google Scholar] [CrossRef]
  30. Babu, C.S.; Madhusoodanan, M.; Sridhar, G.; Tembe, B.L. Orientations of [Fe(H2O)6]2+ and [Fe(H2O)6]3+ complexes at a reactive separation in water. J. Am. Chem. Soc. 1997, 119, 5679–5681. [Google Scholar] [CrossRef]
  31. D’Angelo, P.; Benfatto, M. Effect of Multielectronic Configurations on the XAFS Analysis at the Fe K Edge. J. Phys. Chem. A 2004, 108, 4505–4514. [Google Scholar] [CrossRef]
  32. Benfatto, M.; Solera, J.A.; García Ruiz, J.; Chaboy, J. Double-channel excitation in the X-ray absorption spectrum of Fe3+ water solutions. Chem. Phys. 2002, 282, 441–450. [Google Scholar] [CrossRef]
  33. Testemale, D.; Brugger, J.; Liu, W.; Etschmann, B.; Hazemann, J.-L. In-situ X-ray absorption study of Iron(II) speciation in brines up to supercritical conditions. Chem. Geol. 2009, 264, 295–310. [Google Scholar] [CrossRef]
  34. Jarzȩcki, A.A.; Anbar, A.D.; Spiro, T.G. DFT Analysis of Fe(H2O) 63+ and Fe(H2O) 62+ Structure and Vibrations; Implications for Isotope Fractionation. J. Phys. Chem. A 2004, 108, 2726–2732. [Google Scholar] [CrossRef]
  35. Remsungnen, T.; Rode, B.M. QM/MM Molecular Dynamics Simulation of the Structure of Hydrated Fe(II) and Fe(III) Ions. J. Phys. Chem. A 2003, 107, 2324–2328. [Google Scholar] [CrossRef]
  36. Ye, Q.; Chen, X.; Zhou, J.; Zhao, H.F.; Chu, W.S.; Zheng, X.S.; Marcelli, A.; Wu, Z.Y. Local hydrated structure of an Fe2+/Fe3+ aqueous solution: An investigation using a combination of molecular dynamics and X-ray absorption fine structure methods. Chin. Phys. C 2013, 37, 038003. [Google Scholar] [CrossRef]
  37. Lee, M.S. Chemical equilibria in ferrous chloride acid solution. Met. Mater. Int. 2004, 10, 387–392. [Google Scholar] [CrossRef]
  38. Ding, K.; Seyfried, W.E. Determination of Fe-Cl complexing in the low pressure supercritical region (NaCl fluid): Iron solubility constraints on pH of subseafloor hydrothermal fluids. Geochim. Cosmochim. Acta 1992, 56, 3681–3692. [Google Scholar] [CrossRef]
  39. Fein, J.B.; Hemley, J.J.; D’Angelo, W.M.; Komninou, A.; Sverjensky, D.A. Experimental study of iron-chloride complexing in hydrothermal fluids. Geochim. Cosmochim. Acta 1992, 56, 3179–3190. [Google Scholar] [CrossRef]
  40. Böhm, F.; Sharma, V.; Schwaab, G.; Havenith, M. The low frequency modes of solvated ions and ion pairs in aqueous electrolyte solutions: Iron(ii) and iron(iii) chloride. Phys. Chem. Chem. Phys. 2015, 17, 19582–19591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Zhao, R.; Pan, P. A spectrophotometric study of Fe(II)-chloride complexes in aqueous solutions from 10 to 100 °C. Can. J. Chem. 2001, 79, 131–144. [Google Scholar] [CrossRef]
  42. Ravel, B.; Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: Data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat. 2005, 12, 537–541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. D’Angelo, P.; Benfatto, M.; Della Longa, S.; Pavel, N.V. Combined XANES and EXAFS analysis of Co2+, Ni2+, and Zn2+ aqueous solutions. Phys. Rev. B 2002, 66, 64209. [Google Scholar] [CrossRef]
  44. D’Angelo, P.; Zitolo, A.; Ceccacci, F.; Caminiti, R.; Aquilanti, G. Structural characterization of zinc(II) chloride in aqueous solution and in the protic ionic liquid ethyl ammonium nitrate by X-ray absorption spectroscopy. J. Chem. Phys. 2011, 135, 154509. [Google Scholar] [CrossRef] [PubMed]
  45. Dominko, R.; Sirisopanaporn, C.; Masquelier, C.; Hanzel, D.; Arčon, I.; Gaberscek, M. On the Origin of the Electrochemical Capacity of Li2Fe0.8Mn0.2{SiO}4. J. Electrochem. Soc. 2010, 157, A1309. [Google Scholar] [CrossRef]
  46. Roscioni, O.M.; D’Angelo, P.; Chillemi, G.; Della Longa, S.; Benfatto, M. Quantitative analysis of XANES spectra of dis{\-}ordered systems based on molecular dynamics. J. Synchrotron Radiat. 2005, 12, 75–79. [Google Scholar] [CrossRef]
  47. Arčon, I.; Kolar, J.; Kodre, A.; Hanžel, D.; Strlič, M. XANES analysis of Fe valence in iron gall inks. X-ray Spectrom. 2007, 36, 199–205. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, W.; Chen, X.; Wang, Y.; Wu, L.; Hu, Y. Experimental and modeling of conductivity for electrolyte solution systems. ACS Omega 2020, 5, 22465–22474. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, W.; Etschmann, B.; Foran, G.J.; Shelley, M.; Brugger, J. Deriving formation constants for aqueous metal complexes from XANES spectra: Zn2+ and Fe2+ chloride complexes in hypersaline solutions. Am. Miner. 2007, 92, 761–770. [Google Scholar] [CrossRef]
  50. Rehr, J.J.; Albers, R.C.; Zabinsky, S.I. High-order multiple-scattering calculations of X-ray-absorption fine structure. Phys. Rev. Lett. 1992, 69, 3397–3400. [Google Scholar] [CrossRef] [PubMed]
  51. Bankura, A.; Santra, B.; DiStasio, R.A., Jr.; Swartz, C.W.; Klein, M.L.; Wu, X. A Systematic Study of Chloride Ion Solvation in Water using van der Waals Inclusive Hybrid Density Functional Theory. Mol. Phys. 2015, 113, 2842–2854. [Google Scholar] [CrossRef] [Green Version]
  52. Ge, L.; Bernasconi, L.; Hunt, P. Linking electronic and molecular structure: Insight into aqueous chloride solvation. Phys. Chem. Chem. Phys. 2013, 15, 13169–13183. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Normalized Fe K-edge XANES spectra of 1 molL−1 FeCl2 (aq) measured in situ from RT to 80 °C, on FeCl2 (aq) at higher concentrations (2, 3, and 4 molL−1) at RT, and on 1 molL−1 FeCl3 (aq) as a reference for Fe3+, and crystalline FeSO4·7H2O as a reference for octahedrally coordinated Fe2+ cations.
Figure 1. Normalized Fe K-edge XANES spectra of 1 molL−1 FeCl2 (aq) measured in situ from RT to 80 °C, on FeCl2 (aq) at higher concentrations (2, 3, and 4 molL−1) at RT, and on 1 molL−1 FeCl3 (aq) as a reference for Fe3+, and crystalline FeSO4·7H2O as a reference for octahedrally coordinated Fe2+ cations.
Molecules 27 00642 g001
Figure 2. The k3—weighted Fe K-edge EXAFS spectra of the 1 molL−1 FeCl2 (aq) in the temperature range from RT to 80 °C, and in the FeCl2 (aq) at higher concentrations (2 molL−1, 3 molL−1, and 4 molL−1) at RT, (blue solid line—experiment, red dashed line—best fit EXAFS model in the R range of 1–5 Å). The spectra are displaced vertically for clarity.
Figure 2. The k3—weighted Fe K-edge EXAFS spectra of the 1 molL−1 FeCl2 (aq) in the temperature range from RT to 80 °C, and in the FeCl2 (aq) at higher concentrations (2 molL−1, 3 molL−1, and 4 molL−1) at RT, (blue solid line—experiment, red dashed line—best fit EXAFS model in the R range of 1–5 Å). The spectra are displaced vertically for clarity.
Molecules 27 00642 g002
Figure 3. Fourier transform magnitudes of the k3—weighted Fe K-edge EXAFS spectra of the 1 molL−1 FeCl2 (aq) in the temperature range from RT to 80 °C, and in the FeCl2 (aq) at higher concentrations (2 molL−1, 3 molL−1, and 4 molL−1) at RT, calculated in the k range of 3.4 Å−1 to 13 Å−1 (blue solid line—experiment, red dashed line—best fit EXAFS model in the R range of 1–5 Å).
Figure 3. Fourier transform magnitudes of the k3—weighted Fe K-edge EXAFS spectra of the 1 molL−1 FeCl2 (aq) in the temperature range from RT to 80 °C, and in the FeCl2 (aq) at higher concentrations (2 molL−1, 3 molL−1, and 4 molL−1) at RT, calculated in the k range of 3.4 Å−1 to 13 Å−1 (blue solid line—experiment, red dashed line—best fit EXAFS model in the R range of 1–5 Å).
Molecules 27 00642 g003
Figure 4. Fourier transform magnitude of k3—weighted Fe EXAFS spectrum of the 3 molL−1 FeCl2 (aq) measured at RT, calculated in the k range of 3.4–13 Å−1. Experiment (black solid line), best-fit EXAFS model in the R range of 1–5 Å (magenta dashed line), and single scattering contribution of the Fe-O (blue line) and Fe-Cl neighbors (red line), and contributions of multiple scatterings with a total path length up to 5 Å.
Figure 4. Fourier transform magnitude of k3—weighted Fe EXAFS spectrum of the 3 molL−1 FeCl2 (aq) measured at RT, calculated in the k range of 3.4–13 Å−1. Experiment (black solid line), best-fit EXAFS model in the R range of 1–5 Å (magenta dashed line), and single scattering contribution of the Fe-O (blue line) and Fe-Cl neighbors (red line), and contributions of multiple scatterings with a total path length up to 5 Å.
Molecules 27 00642 g004
Figure 5. Schematic presentation of two major Fe-ionic species in FeCl2 (aq). Fe2+ (yellow) in the monovalent [FeCl(H2O)5]+ is coordinated by five water molecules (red and gray) and an apical Cl ion (green). The intramolecular hydrogen bond restrains movements of the basal water molecules and destabilizes the apical water molecule. The neutral Fe[Cl2(H2O)4]0 is formed by substitution of the apical water molecule by Cl.
Figure 5. Schematic presentation of two major Fe-ionic species in FeCl2 (aq). Fe2+ (yellow) in the monovalent [FeCl(H2O)5]+ is coordinated by five water molecules (red and gray) and an apical Cl ion (green). The intramolecular hydrogen bond restrains movements of the basal water molecules and destabilizes the apical water molecule. The neutral Fe[Cl2(H2O)4]0 is formed by substitution of the apical water molecule by Cl.
Molecules 27 00642 g005
Figure 6. The conductivity of FeCl2 (aq) as a function of concentration (blue dots) and the corresponding fit (red line) of the conductivity model developed by Zhang et al. [48].
Figure 6. The conductivity of FeCl2 (aq) as a function of concentration (blue dots) and the corresponding fit (red line) of the conductivity model developed by Zhang et al. [48].
Molecules 27 00642 g006
Figure 7. CAD drawing and photograph of the experimental PMMA liquid absorption cell, designed and manufactured on-purpose for the reported XAS measurements.
Figure 7. CAD drawing and photograph of the experimental PMMA liquid absorption cell, designed and manufactured on-purpose for the reported XAS measurements.
Molecules 27 00642 g007
Table 1. Parameters of the nearest coordination shells around Fe cations in the FeCl2 samples: average number of neighbor atoms (N), distance (R), and Debye–Waller factor (σ2). Uncertainty of the last digit is given in parentheses. The best fit is obtained with a total number of O and Cl neighbors constrained to 6, the amplitude reduction factor S02 = 0.85, and the shift of the energy origin ΔEo = −4 ± 1 eV. The R-factor (quality of fit parameter) is listed in the last column. Uncertainty of the last digit is given in parentheses.
Table 1. Parameters of the nearest coordination shells around Fe cations in the FeCl2 samples: average number of neighbor atoms (N), distance (R), and Debye–Waller factor (σ2). Uncertainty of the last digit is given in parentheses. The best fit is obtained with a total number of O and Cl neighbors constrained to 6, the amplitude reduction factor S02 = 0.85, and the shift of the energy origin ΔEo = −4 ± 1 eV. The R-factor (quality of fit parameter) is listed in the last column. Uncertainty of the last digit is given in parentheses.
Fe NeighborsNR [Å]σ22]R-Factor
1 molL−1 FeCl2 aq. solution measured at RT
O5.0(3)2.095(5)0.0072(5)0.00030
Cl1.02.31(2)0.015(3)
1 molL−1 FeCl2 aq. solution measured at T = 30 °C
O4.9(3)2.096(5)0.0071(5)0.00020
Cl1.12.34(2)0.015(3)
1 molL−1 FeCl2 aq. solution measured at T = 40 °C
O4.8(3)2.093(5)0.0072(5)0.00016
Cl1.22.33(2)0.015(3)
1 molL−1 FeCl2 aq. solution measured at T = 50 °C
O4.6(3)2.095(5)0.0073(5)0.00013
Cl1.42.33(2)0.017(3)
1 molL−1 FeCl2 aq. solution measured at T = 60 °C
O4.5(3)2.095(5)0.0075(5)0.00009
Cl1.52.33(2)0.018(3)
1 molL−1 FeCl2 aq. solution measured at T = 80 °C
O4.5(3)2.094(5)0.0082(5)0.00012
Cl1.52.33(2)0.017(3)
2 molL−1 FeCl2 aq. solution measured at RT
O4.5(3)2.103(5)0.0065(5)0.00015
Cl1.52.35(2)0.021(3)
3 molL−1 FeCl2 aq. solution measured at RT
O4.0(3)2.097(5)0.0066(5)0.00012
Cl2.02.37(2)0.019(3)
4 molL−1 FeCl2 aq. solution measured at RT
O4.0(3)2.097(5)0.0066(5)0.00012
Cl2.02.37(2)0.019(3)
Table 2. Number of O and Cl neighbors in the first coordination shell around Fe2+ cations in the selected temperature and concentration range, as determined by EXAFS analysis. The average Fe-O and Fe-Cl distances are 2.095 ± 0.005 Å and 2.33 ± 0.02 Å, respectively. Uncertainty in the coordination numbers is ±0.3. The total number of neighbors in the first coordination shell in the EXAFS fit was fixed to 6.
Table 2. Number of O and Cl neighbors in the first coordination shell around Fe2+ cations in the selected temperature and concentration range, as determined by EXAFS analysis. The average Fe-O and Fe-Cl distances are 2.095 ± 0.005 Å and 2.33 ± 0.02 Å, respectively. Uncertainty in the coordination numbers is ±0.3. The total number of neighbors in the first coordination shell in the EXAFS fit was fixed to 6.
Temperature
(°C)
Average Number of Neighbors in 1 molL−1 FeCl2FeCl2 (aq) C. (molL−1)Average Number of Neighbors at RT
ClOClO
201.05.011.05.0
301.14.921.54.5
401.24.832.04.0
501.44.642.04.0
601.54.5
801.54.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luin, U.; Arčon, I.; Valant, M. Structure and Population of Complex Ionic Species in FeCl2 Aqueous Solution by X-ray Absorption Spectroscopy. Molecules 2022, 27, 642. https://doi.org/10.3390/molecules27030642

AMA Style

Luin U, Arčon I, Valant M. Structure and Population of Complex Ionic Species in FeCl2 Aqueous Solution by X-ray Absorption Spectroscopy. Molecules. 2022; 27(3):642. https://doi.org/10.3390/molecules27030642

Chicago/Turabian Style

Luin, Uroš, Iztok Arčon, and Matjaz Valant. 2022. "Structure and Population of Complex Ionic Species in FeCl2 Aqueous Solution by X-ray Absorption Spectroscopy" Molecules 27, no. 3: 642. https://doi.org/10.3390/molecules27030642

APA Style

Luin, U., Arčon, I., & Valant, M. (2022). Structure and Population of Complex Ionic Species in FeCl2 Aqueous Solution by X-ray Absorption Spectroscopy. Molecules, 27(3), 642. https://doi.org/10.3390/molecules27030642

Article Metrics

Back to TopTop