Introduction
Reactions at a solid–liquid interphase play a vital role in understanding a dissolution phenomenon. Recent research efforts have helped to improve dissolution-based processes such as the regeneration of ceramic filter media (Salmimies et al., Reference Salmimies, Kallas, Ekberg, Görres, Andreassen, Beck and Häkkinen2013; Smith et al., Reference Smith, Sheridan, van Dyk, Naik, Plint and Turrer2018) and leaching iron from clays (Mandal & Banerjee, Reference Mandal and Banerjee2004). Sulfuric and nitric acid are commonly used leaching agents for iron oxides, but more effective dissolution can be achieved using a solution that contains complexing ligands. Complexing agents such as ethylenediaminetetraacetic acid (EDTA), trisodium nitrilotriacetate (NTA), hydroxyethylethylenediaminetriacetic acid (HEDTA), pentetic acid (DTPA), and methanesulfonic acid (MSA) have been used to enhance the dissolution of iron oxides (Chang & Matijevic, Reference Chang and Matijevic1983; Blesa et al., Reference Blesa, Borghi, Maroto and Regazzoni1984; Wang et al., Reference Wang, Fu, Zhang, Song and Chen2017). Among these ligands, oxalic acid has also received a lot of attention as it can dissociate into strong complexing oxalate ions for the dissolution (Cornell & Schindler, Reference Cornell and Schindler1987; Panias et al., Reference Panias, Taxiarchou, Paspaliaris and Kontopoulos1996; Lee et al., Reference Lee, Tran, Park and Kim2007; Salmimies et al., Reference Salmimies, Vehmaanperä and Häkkinen2016). Oxalic acid was chosen for the present study because it is environmentally friendly. Moreover, it can also dissolve aluminum and manganese, which are often also present in clay minerals (Rennert et al., Reference Rennert, Dietel, Heilek, Dohrmann and Mansfeldt2021).
The reaction scheme for dissolution of hematite in oxalic acid can be described via Eqs. 1–10 (Stumm & Furrer, Reference Stumm, Furrer and Stumm1987; Siffert & Sulzenberger, Reference Siffert and Sulzberger1991; Panias et al., Reference Panias, Taxiarchou, Paspaliaris and Kontopoulos1996).
Dissociation of oxalic acid:
Surface protonation:
Surface complexation:
Electron transfer from ox to Fe(III) and detachment of Fe(II)
Autocatalytic dissolution:
Symbol ‘>’ describes the solid surface, L is the ligand, II or III is the oxidation state of iron in the solid phase, 2+ or 3+ is the oxidation state of iron in the liquid, ‘...’ is the adsorbed species on the solid surface, and ox refers to species formed from oxalic acid, i.e. HC2O4 - or C2O4 -
In previous studies (Lee et al., Reference Lee, Tran, Park and Kim2007; Salmimies et al., Reference Salmimies, Vehmaanperä and Häkkinen2016), the formation of a solid precipitate was observed during the dissolution of hematite in oxalic acid. This reaction product, however, was not considered in the original reaction scheme (Eqs. 1–10). Often, solid precipitation is a drawback that hinders the dissolution process, blocks the filter media, or reduces the whiteness of the clay mineral, which cannot then be used again as a raw material in the ceramics industry. This precipitate is humboldtine, named after Friedrich Wilhelm Heinrich Alexander von Humboldt (1769–1859), naturalist, explorer, and scientist. The formation of humboldtine can be described by Eq. 11, which can take place after the reduction of surface Fe(III) to Fe(II) (Eq. 5).
The thermal decomposition of humboldtine reverses this reaction by restoring both magnetite (Angermann & Töpfer, Reference Angermann and Töpfer2008) and hematite (Ogasawara & Koga, Reference Ogasawara and Koga2014). Such thermal decomposition may also occur on Mars (Applin et al., Reference Applin, Izawa, Cloutis, Goltz and Johnson2015). The formation of humboldtine is of great interest because, though it can have drawbacks in terms of dissolution, its use has recently been demonstrated for high conductivity (Yamada et al., Reference Yamada, Sadakiyo and Kitagawa2009; Huskić et al., Reference Huskić, Pekov, Krivovichev and Frščić2016; Yao et al., Reference Yao, Armstrong, Zhou, Sougrati, Kidkhunthod, Tunmee, Sun, Sattayaporn, Lightfoot, Ji, Jiang, Wu, Tang and Cheng2019) and for photocatalysis (Fan et al., Reference Fan, Zhang, Li, Wang, Zhou, Cheng, Zhou and Shi2016), and thus its facile synthesis has become a relevant research topic (Müller et al., Reference Müller, Bourcet and Hanfland2021). The conditions under which humboldtine is formed in nature remain unknown (Robinson, Reference Robinson2004), but formation has been observed when pyrite rocks are in contact with oxalic acid (Green, Reference Green2004). Moreover, the formation of humboldtine has also been associated with oxalic acid produced by microorganisms near iron surfaces (Aramendia et al., Reference Aramendia, Gomez-Nubla, Bellot-Gurlet, Castro, Arana and Madariaga2015). Clearly, the formation of humboldtine in nature is dependent on the iron composition of the rock substrate (Gadd et al., Reference Gadd, Bahri-Esfahani, Li, Rhee, Wei, Fomina and Liang2014). Very often humboldtine is synthesized via several steps (Müller et al., Reference Müller, Bourcet and Hanfland2021) whereas the formation of humboldtine via dissolution-precipitation reaction under mild conditions (at the temperature of 35°C) could offer a simpler route. Understanding the mechanism causing humboldtine formation under dissolution conditions is, therefore, important. The reaction steps shown in Eqs. 1–10 are, thus, hypothetical. DFT calculations offer a unique window to visualize the actual mechanisms taking place during the dissolution, and such calculations have already been deployed successfully to study the hematite surface (Huang et al., Reference Huang, Kumar Ramadugu and Mason2016; Si et al., Reference Si, Li, Zhou, Liu and Prezhdo2020) and a hematite and benzene system (Dzade et al., Reference Dzade, Roldan and de Leeuw2014).
The aim of the current study, therefore, was to examine the conditions and mechanisms leading to the formation of humboldtine via dissolution of hematite in oxalic acid, with emphasis on the DFT calculations combined with experimental verification. A further objective was to verify the formation of humboldtine in the specific acid system being studied here. The hypothesis was that the DFT calculations would define the actual dissolution mechanisms taking place in the system and provide understanding of the adsorption, electron-transfer, and precipitation processes.
Materials and Methods
Computational Methodology
α-Fe2O3 with lattice parameters a = b = 5.036 Å, c = 13.75 Å, α = β = 90°, and γ = 120° was used here, consistent with previous experimental work (Blake et al., Reference Blake, Hessevick, Zoltai and Finger1966). The (0001) surface with a 2×2 supercell and the nine-layer periodic slab model was constructed. A >12 Å vacuum length was set in the z direction to avoid interactions between the periodic images. The surface was optimized with the atoms on the uppermost four layers allowed to relax. The adsorption energies (E ads) were calculated with the equation: E ads = E complex – (E molecule + E surface), where E complex, E molecule, and E surface represent the total energy of the complex after adsorption, the energy of the free molecule (e.g. H2C2O4, H2O), and energy of the free surface, respectively. The free molecules were also optimized with the same parameters as used for the surface. The Bader charges for all configurations were calculated for the analysis of the electron-transfer processes (Henkelman et al., Reference Henkelman, Arnaldsson and Jónsson2006). In the Bader charge-analysis scheme, the simulation system was partitioned into different volumes by the surfaces of minimum electronic density. The electronic density within one volume is then assigned to the atom in that volume. The Bader charge analysis scheme is a powerful tool which allows the change in the atomic charges during the reaction to be analyzed. Furthermore, the oxidation state-determination method to examine the changes in the oxidation states of Fe during the reaction (Sit et al., Reference Sit, Car, Cohen and Selloni2011) was also adopted. In this method, a 5×5 occupation matrix was constructed for each electron spin by projecting the Kohn-Sham wave functions from the DFT calculations onto the atomic d orbitals of the transition metal ion. The eigenvalues of the matrix are called the occupation numbers. As shown earlier (Sit et al., Reference Sit, Car, Cohen and Selloni2011), the occupation numbers close to unity correspond to the d electrons belonging to the ion. The oxidation state can then be assigned with knowledge of the number of d electrons in the ion. The method has been applied successfully to study the oxidation states of transition metal ions in various problems (Sit et al., Reference Sit, Cohen and Selloni2012; Majumdar et al., Reference Majumdar, Tse, Wu and Yao2017; Wei et al., Reference Wei, Huang, Wang, Zhang, Ge, Ma, Wen, Zhang, Li, Lei, Zhang, Irawan, Liu and Wu2017; Ku and Sit, Reference Ku and Sit2019; Ricca et al., Reference Ricca, Timrov, Cococcioni, Marzari and Aschauer2019; Wang et al., Reference Wang, Qiu, Long, Lu, Hu, Liu and Zhang2019; Jensen et al., Reference Jensen, Tan, Feng, Li, Zhou and Thonhauser2020). Further details of this method can be found in Sit et al. (Reference Sit, Car, Cohen and Selloni2011). All the DFT calculations were performed using the PWscf module of the Quantum Espresso (QE) package (Giannozzi et al., Reference Giannozzi2009). The generalized gradient approximation (GGA) proposed by Perdew, Burke, and Ernzerh (PBE) was used to describe the exchange correlation energy (Perdew et al., Reference Perdew, Burke and Ernzerhof1996). The ultrasoft pseudopotentials were used to simulate electron−ionic core interactions (Vanderbilt, Reference Vanderbilt1990). The van der Waals interaction correction with the DFT-D2 scheme proposed by Grimme (Reference Grimme2006) was used to simulate the weak interactions. The cutoff energies were set equally for all the calculations, which were 30 Rydberg (408 eV) and 240 Rydberg (3265 eV) for the wavefunctions and the augmented part of the charge density, respectively. The 2×2×1 Monkhorst-Pack k-point mesh was found to be adequate in this study after testing. Spin-polarized calculations were carried out for the systems involving α-Fe2O3 because α-Fe2O3 has unpaired electrons. Also, the Hubbard U value was set to 5 eV to consider accurately the strongly correlated d electrons (Dzade et al., Reference Dzade, Roldan and de Leeuw2014). Constrained optimization simulations were carried out to calculate the activation barrier of the C–C bond of the H2C2O4.
Materials
Analytical-grade oxalic dihydrate powder from Merck (99%, Darmstadt, Germany) and Millipore-water were used to prepare 0.33 M oxalic acid solution for the dissolution experiment. The solid synthetic hematite powder was from Alfa Aesar (98%, Kandel, Germany). Analytical-grade, concentrated nitric acid from Merck (65%, Darmstadt, Germany) was used for the preparation of the nitric acid solution for the Fe concentration analysis.
Dissolution Experiment
The dissolution experiment followed a similar set-up as in the previous studies of the authors (Salmimies et al., Reference Salmimies, Vehmaanperä and Häkkinen2016; Vehmaanperä et al., Reference Vehmaanperä, Salmimies and Häkkinen2021). The experiment was done in a 1 L water-jacketed glass reactor, which was equipped with a Lauda Proline RP855 thermostat to keep the reaction temperature constant. The system was equilibrated with air and the experiment was done under light. A pitched-blade turbine with a mixing speed of 800 rpm and four baffles ensured effective mixing. The experiment was run in the following manner: (1) oxalic acid solution was heated to 35°C; (2) 12 g of hematite powder was added to the reactor; and (3) mixing was switched on. Solution samples were collected from the reactor with a syringe, filtered using a 0.22 μm syringe filter, and diluted with 14 wt.% nitric acid to avoid any changes in the samples before analysis.
The slurry was collected after the experiment and filtered using a Büchner funnel equipped with a Whatman® 42 filter paper. After this, the solids were dried in an oven at 40°C. According to Schwertmann & Cornell (Reference Schwertmann and Cornell2003), this temperature is low enough to avoid any phase modification or transformation caused by oven drying. After drying, the solids were ground gently in a mortar before XRD and SEM-EDS analysis.
Equilibrium was reached within 413 h, but the solubility data are not reported here because the present investigation was of the formation of humboldtine using DFT with experimental verification. Experimental data and discussion of the kinetics and equilibrium of dissolution of iron oxides collected with the described set-up have been reported elsewhere (Salmimies et al., Reference Salmimies, Vehmaanperä and Häkkinen2016; Vehmaanperä et al. Reference Vehmaanperä, Salmimies and Häkkinen2021).
Analysis
The pH of solutions from the reactor was measured using a WTW pH 401 i-meter with a WTW SenTix 41 electrode (Xylem Analytics, Weinheim, Germany). The pH ranged from 0.8 to 1.0 during the dissolution, which was favorable for the formation of humboldtine (Lee et al., Reference Lee, Tran, Park and Kim2007). The concentration of total dissolved Fe was determined using a flame atomic absorption spectrometer (Thermo Scientific iCE 3000 AAS with a detection limit of 0.02 mg/L for Fe, China).
The mineralogical composition of the residual solid was analyzed using X-ray diffraction (XRD, Bruker D8 Focus with CuKα radiation, Karlsruhe, Germany). The XRD data were recorded from 10 to 70°2θ with a step size of 0.02°2θ and a scanning speed of 0.02°2θ s–1. The phase identification was done using the DIFFRAC.EVA software and the PDF 4 database. A scanning electron microscope (SEM, Hitachi SU 3500 with Thermo Scientific, UltraDry SDD EDS, Madison, Wisconsin, USA) was used to estimate the elemental composition and morphology of the solids.
Results
DFT Calculations
As in the paper by Dzade et al. (Reference Dzade, Roldan and de Leeuw2014), the present study considered the (0001) surface of hematite because it is the dominant surface obtained under most conditions (Parkinson, Reference Parkinson2016). The adsorption of oxalic acid onto the Fe2O3 surface and the charge transfer were studied using Bader charge analysis and oxidation state determination. For oxalic acid adsorption, spontaneous proton transfer from the oxalic acid to the surface took place. Therefore, only stable structures for singly de-protonated and oxalate adsorption could be obtained. For both cases, electron transfer from the molecule to the surface was observed (Table 1), suggesting that the molecule was oxidized partially upon adsorption. Oxalate adsorption was stronger than singly de-protonated adsorption, and it also showed a larger electron transfer.
A main objective of the DFT study was to understand the adsorption and the electron-transfer processes. In the charge analysis, although oxalate was oxidized, the reduction of the surface Fe was not very significant. A possible reason is that the electron from the oxalate was shared by the whole Fe2O3 slab instead of reducing a particular Fe. This fact was also confirmed by Stone and Morgan (Reference Stone, Morgan and Stumm1987). Clearly, the reduction of Fe did not occur during adsorption as supported by the oxidation states of the Fe (Table 2).
When the C–C bond broke, reduction of two surface Fe atoms from Fe(III) to Fe(II) was noted. The oxidation state calculations method through DFT calculations allowed assignment of d electrons to various Fe ions, and thus to extract their oxidation states. The present analysis indicated that the reaction proceeded first as the reduction of surface Fe(III) to Fe(II) by the C2O4 2– → 2CO2 reaction. Then the remaining oxalate reacted with the surface Fe(II) to form humboldtine. The top view and side view of the structure after the C–C bond was broken with four representative Fe positions are shown in Fig. 1. The DFT results revealed specific locations on the surface, namely, Fe(1) and Fe(2), that were the most favorable sites for initiating the formation of the precipitate.
Experimental Findings
To establish the formation of humboldtine, the solid phase before and after dissolution was analyzed using SEM-EDS and XRD. The SEM-EDS data were as follows: 14 wt.% of carbon, 41 wt.% of iron, and ~1 wt.% of other impurities such as aluminum and silicon. The calculated oxygen concentration was 44 wt.%. The XRD analysis showed that a transformation of the solid phase took place as a result of dissolution (Fig. 2). The XRD pattern was in good agreement with hematite before the dissolution. After dissolution, the pattern showed humboldtine in the α monoclinic phase, a small amount of hematite, and a small fraction of humboldtine in the β phase. The XRD peaks were narrow and sharp, indicating a well-formed crystal structure in both cases.
A remarkable change in color was observed between the samples, indicating the formation of humboldtine which has a typical yellow to amber-yellow color (Lide, Reference Lide1999). The original hematite powder was red whereas the precipitate appeared yellow after the dissolution (Fig. 3a). The SEM images showed very large crystals and chain structures (Fig. 3c,d), in good agreement with the crystallographic structure of humboldtine (Fig. 3b). The DFT results (described above) showed that the iron atoms of the hematite surface can convert to Fe(II) and can precipitate in the presence of oxalates into large quantities of humboldtine crystals given the low solubility (K sp(ferrous oxalate) = 3.2∙10−7) (Speight et al., Reference Speight, Lange and Dean2005).
SUMMARY AND CONCLUSIONS
The present study showed that the adsorption of oxalate ion to the surface of hematite and the electron-transfer of Fe(III) to Fe(II) facilitated the formation of humboldtine. Considering the phase diagram (Lee et al., Reference Lee, Tran, Park and Kim2007) and the XRD results, excess amounts of oxalic acid are seen to favor the formation of crystalline humboldtine. The 1:2 molar ratio of Fe/C2O4 was required to obtain humboldtine as 1 mole of oxalate was required for the electron transfer of Fe(III) to Fe(II), and another 1 mole was needed to form humboldtine. In theory, ~25 g of humboldtine could be formed using 12 g of synthetic hematite and 40 g of oxalic acid dihydrate powder.
If the target was to dissolve hematite without forming humboldtine, this could be achieved by adding the oxalic acid gradually. In addition, humboldtine formation was unlikely to be observed if the reaction scheme was shifted toward carbon dioxide formation.
The most important result of the current study was the discovery of a new pathway to produce large crystals of humboldtine of good quality. The previous method by which humboldtine was synthesized involved several steps, whereas the present method requires only mixing of oxalic acid with hematite powder at the appropriate temperature of 35°C and for the appropriate experimental time of < 400 h. The XRD analysis showed that the crystal structure was in the monoclinic α phase suitable for applications requiring significant ion mobility (Yamada et al., Reference Yamada, Sadakiyo and Kitagawa2009). The samples could be used for anode material in lithium-ion batteries, therefore (Zhang et al., Reference Zhang, Li, Wang, Zhao, Chen, Dai and Yao2020). Moreover, the DFT calculations provided important atomistic details about the structures after oxalic acid adsorption, and information about the direction and magnitude of the charge-transfer process, which were exhibited by both the Bader charge and oxidation state analysis. Oxalic acid was doubly de-protonated and adsorbed on the Fe2O3 (0001) surface. The adsorption process alone did not lead to surface reduction. The charge transfer and oxidation state analysis showed that cleavage of the adsorbed oxalate C–C bond was key to reducing two surface Fe(III) to Fe(II), providing the sites for subsequent reactions of the oxalate to form humboldtine.
Acknowledgments
P.V., R.S., and A.H. designed and performed the experiment, and analyzed the data. B.G., P.H.-L.S., and B.B. were responsible for the DFT calculations and the theoretical interpretation of the results. Toni Väkiparta is kindly acknowledged for providing the SEM images. B.G. and P.H.-L.S. acknowledge the computer facility managed by the School of Energy and Environment (SEE) of the City University of Hong Kong. This research was supported by the Ministry of Education and Culture (Finland).
FUNDING
Open access funding has been provided by LUT University (previously Lappeenranta University of Technology (LUT)). Funding sources are as stated in the Acknowledgments.
Declarations
Conflict of Interest
The authors declare that they have no conflict of interest.