Fe(II) hydrolysis in aqueous solution: A DFT study

Insights into rare earth element speciation: unraveling sulfate and hydrolysis complexes through DFT calculations

CONTEXT

Rare earth elements (REE) are indispensable in numerous green technologies owing to their exceptional physical and chemical attributes. Separating REE is a pivotal process to meet the increasing demands of the high-tech industry. Understanding the hydrolysis of REE in aqueous environments marks the initial stride in comprehending their separation mechanisms. Sulfate commonly coexists in high-concentration solutions alongside REE, stemming from mineral processing. We analyzed the hydrolysis of REE and their complexes with sulfate using DFT methods. We present and discuss on the structural characteristics of hydrolysis species and sulfate complexes in alignment with existing experimental data. Estimates of Gibbs free energies for hydrolysis and sulfate complex formation were compared against literature values. REE pose challenges owing to the labile nature of aqua complexes and the pivotal role of system dynamics. We showed that hydrolysis reactions could be suitably modeled, yielding an error margin of approximately 5 kcal mol-1 concerning experimental values, employing the M06 exchange-correlation functional with the SMD implicit solvation model. However, sulfate chemical species proved to be more challenging, exhibiting larger error margins with substantial variations across the REE series. The Raman spectrum analysis of lanthanum sulfate complexes demonstrated excellent agreement with experimental values.

METHOD

We applied the M06, PBE, and PBE0 exchange-correlation functionals combined with def2-TZVP basis sets and SMD to obtain the Gibbs free energies of hydrolysis and sulfate complexation with lanthanides in aqueous solution. The calculations were performed using the ORCA program.

Exploring the electrochemical properties and lithium insertion mechanisms in akaganeite (β-FeOOH) - a combined DFT/experimental study

Akaganeite (β-FeOOH) has been intensely investigated to be used in different electrochemical applications as a cathode material in Li-ion batteries owing to its unique structural characteristics, including channels capable of accommodating and reversibly extracting charged species such as lithium or sodium ions. We revisited the synthesis, and its electrochemical properties based on a combined experimental/theoretical approach aiming to understand the mechanism of the electron transfer in this material. Electrochemical investigations, employing Li2SO4 aqueous electrolyte, unveiled notable alterations in the charge/discharge profiles. The initial discharge curve revealed distinct plateaus at 3.4 V and 2.9 V, with the absence of the former in subsequent cycles, indicating irreversible reactions in the initial cycle. Furthermore, density functional theory (DFT) calculations were employed to elucidate the impact of lithium atom insertion on the electronic and structural properties of akaganeite. We gained insights into the underlying electrochemical processes calculating band structures, density of states, and topological analysis based on Bader's theory. The calculated oxidation potentials (3.2 V) closely matched experimental observations, attributing the 3.2 V plateau to lithium insertion into the akaganeite structure.

Solvation energies of the ferrous ion in water and in ammonia at various temperatures

CONTEXT

The solvation of metal ions is crucial to understanding relevant properties in physics, chemistry, or biology. Therefore, we present solvation enthalpies and solvation free energies of the ferrous ion in water and ammonia. Our results agree well with the experimental reports for the hydration free energy and hydration enthalpy. We obtained [Formula: see text] kJ mol[Formula: see text] for the hydration free energy and [Formula: see text] kJ mol[Formula: see text] for the hydration enthalpy of ferrous ion in water at room temperature. At ambient temperature, we obtained [Formula: see text] kJ mol[Formula: see text] as the [Formula: see text] ammoniation free energy and [Formula: see text] kJ mol[Formula: see text] for the ammoniation enthalpy. In addition, the free energy of solvation is deeply affected when the temperature increases. This pattern can be attributed to the rise of entropy when the temperature rises. Besides, the temperature does not affect the ammoniation enthalpies and the hydration enthalpy of the [Formula: see text] ion.

METHOD

All the geometry optimizations are performed at the MP2 methods associated with the 6-31++g(d,p) basis set of Pople. solvated phase structures of [Formula: see text] ion in water or in ammonia are performed using the PCM model. The [Formula: see text] program suite was used to perform all the calculations. The program TEMPO was also used to evaluate the temperature sensitivity of the different obtained geometries.

Iron species activating chlorite: Neglected selective oxidation for water treatment

Chlorite (ClO₂ ^-) is the by-product of the water treatment process carried out using chlorine dioxide (ClO₂) as an effective disinfectant and oxidant; however, the reactivation of ClO₂ ^- has commonly been overlooked. Herein, it was unprecedentedly found that ClO₂ ^- could be activated by iron species (Fe_b: Fe⁰, Fe^II, or Fe^III), which contributed to the synchronous removal of ClO₂ ^- and selective oxidative treatment of organic contaminants. However, the above-mentioned activation process presented intensive H⁺-dependent reactivity. The introduction of Fe_b significantly shortened the autocatalysis process via the accumulation of Cl^- or ClO^- during the protonation of ClO₂ ^- driven by ultrasonic field. Furthermore, it was found that the interdependent high-valent-Fe-oxo and ClO₂, after identification, were the dominant active species for accelerating the oxidation process. Accordingly, the unified mechanisms based on coordination catalysis ([Fe ^N (H₂O) _a (ClO _x ^m-) _b ] ^{n +}-P) were putative, and this process was thus used to account for the pollutant removal by the Fe_b-activated protonated ClO₂ ^-. This study pioneers the activation of ClO₂ ^- for water treatment and provides a novel strategy for "waste treating waste". Derivatively, this activation process further provides the preparation methods for sulfones and ClO₂, including the oriented oxidation of sulfoxides to sulfones and the production of ClO₂ for on-site use.

A DFT study of the adsorption of O2 and [Fe(H2O)2(OH)3] on the (001) and (112) surfaces of chalcopyrite

The oxidation of chalcopyrite, CuFeS₂, is still not well understood and relevant in the context of the hydrometallurgical extraction of copper. Herein, we used DFT calculations within the periodic boundary conditions formalism to study the adsorption of O₂ and [Fe(H₂O)₂(OH)₃] molecules on the (001) and (112) surfaces of CuFeS₂. The O₂ molecule adsorbs strongly by a dissociative pathway at sulfur atoms on the (001) surface with an adsorption energy of - 76.5 kcal mol^-1. The surface is chemically modified forming SO₂ groups, in which the S-O bond length is calculated to be 1.47 and 1.54 Å. PDOS and Löwdin charges analyses indicate the oxidation of the sulfur atoms on the surface. We tested different adsorption modes of [Fe(H₂O)₂(OH)₃], and a bidantade coordination with the O_ads-Fe_sur and Fe_ads-S_sur bond lengths of 2.02 and 2.47 Å is the most favorable with an adsorption energy of - 18.8 kcal mol^-1 on the (001) surface. Adsorptions of each species are also observed on the (112) surface, but they are weaker than those observed on the (001) surface.

Exploring the kinetics of actinyl-EDTA reduction by ferrous iron using quantum-mechanical calculations

The reduction of An(vi) (An = U, Np, and Pu) to An(iv) significantly decreases its solubility and mobility. This reaction can be hindered by complexation with inorganic (e.g., carbonate) or organic ligands. Ethylenediaminetetraacetic acid (EDTA) is one such organic ligand that forms stable complexes with actinides. Therefore, it may enhance the mobility of actinides. However, the redox kinetics and mechanisms of actinyl (An(v/vi)O2+/2+)-EDTA are not well characterized yet and are thus studied here using quantum-mechanical calculations. The principle is to approach the actinyl-EDTA and Fe2+ (reductant) in small incremental steps and calculate the system energy at each distance. The overall reaction is then delineated into sub-processes (encounter frequency in bulk solution, formation of outer-sphere complex, transition from outer- to inner-sphere complex, and electron transfer), and reaction rates are determined for each sub-process. The formation of outer-sphere complexes occurs rapidly in microseconds to seconds over a wide range of actinyl concentrations (pM to μM); in contrast, the transition to the inner-sphere complex is relatively slow (milliseconds to a few seconds). Immediate electron transfer to form the pentavalent actinide is observed along the reaction path for Np(vi) and Pu(vi), but not for U(vi). Surprisingly, in acidic conditions, one of the carboxylic groups gets protonated in EDTA of [UO2(edta)]2- rather than one of the amino groups. This process-based series of calculations can be applied to any redox reaction and allows the prediction of changes to the rate law and rate-limiting step in a more fundamental way for different environments.

Molecular Design of Novel Chemicals for Iron Sulfide Scale Removal

Scale deposition is a pertinent challenge in the oil and gas industry. Scales formed from iron sulfide are one of the troublous scales, particularly pyrite. Moreover, the use of biodegradable environmentally friendly chemicals reduces the cost compared to the conventional removal process. In this work, the chelating abilities of four novel chemicals, designed using the in silico technique of density functional theory (DFT), are studied as potential iron sulfide scale removers. Only one of the chemicals containing a hydroxamate functional group had a good chelating ability with Fe2+. The chelating strength and ecotoxicological properties of this chemical were compared to diethylenetriaminepentaacetic acid (DTPA), an already established iron sulfide scale remover. The new promising chemical surpassed DTPA in being a safer chemical and having a greater binding affinity to Fe2+ upon optimization, hence, a better choice. The presence of oxime (-NHOH) and carbonyl (C=O) moieties in the new chemical showed that the bidentate form of chelation is favored. Moreover, the presence of an intramolecular hydrogen bond enhanced its chelating ability.

Development of a New Borax-Based Formulation for the Removal of Pyrite Scales

In the oil and gas industry, pyrite forms one of the most hardened scales in reservoirs, which hinders the flow of fluids. Consequently, this leads to blockage of the downhole tubular, formation damage, and complete shutdown of production and operational processes. Herein, a new green formulation based on borax (K2B4O7) is proposed for pyrite scale removal. The temperature effect, disk rotational speed, and borax concentration have been investigated using a rotating disk apparatus. Also, XPS and SEM-EDX analyses were conducted on the pyrite disk surface before and after the treatment with the green formulation. The new formulation showed the potential ability to dissolve pyrite without generating the toxic hydrogen sulfide (H2S). The dissolution rate of the scale in the new formulation is increased by 16% compared to that in a previous green formulation composed of 20 wt %DTPA+9 wt % K2CO3. Molecular modeling technique using DFT was used to study the solvation energies of Fe2+ and Fe3+. The latter had a higher solvation energy than the former, which confirmed that upon using the borax-based formulation to oxidize Fe2+ to Fe3+. It will aid the dissolution of pyrite scales. The new formulation achieved a corrosion rate that is 25 times lower than that of 15 wt % HCl, which is commercially used in treating scales. Finally, the proposed new formulation does not require the use of corrosion inhibitors; hence, it is expected to result in a more economical scale treatment method.

Thermodynamic Simulation of Carbonate Cements-Water-Carbon Dioxide Equilibrium in Sandstone for Prediction of Precipitation/Dissolution of Carbonate Cements

Carbonate cements, such as calcite, dolomite, ferrocalcite and ankerite, play important roles in the formation of pores in sandstones: precipitation of carbonate cements modifies pores and inhibits compaction, while dissolution creates secondary pores. This work proposed a precipitation-dissolution model for carbonate cements-CO2-H2O system by means of ion equilibrium concentration ([M2+], M = Ca, Mg, Fe or Mn) with different factors, such as temperature, depth, pH, [Formula: see text], variable rock composition and overpressure. Precipitation-dissolution reaction routes were also analyzed by minimization of the total Gibbs free energy (ΔG). Δ[M2+], the variation of [Ca2+], [Fe2+], [Mg2+] or [Mn2+] for every 100 m of burial depths, is used to predict precipitation or dissolution. The calculation results indicate that the increasing temperature results in decrease of equilibrium constant of reactions, while the increasing pressure results in a relatively smaller increase of equilibrium constant; As a result, with increasing burial depth, which brings about increase of both temperature and pressure, carbonate cements dissolve firstly and produces the maximal dissolved amounts, and then precipitation happens with further increasing depth; For example, calcite is dissolving from 0.0 km to 3.0 km with a maximal value of [Ca2+] at depth of 0.8 km, and then precipitates with depth deeper than 3.0 km. Meanwhile, with an increasing CO2 mole fraction in the gaseous phase from 0.1% to 10.0% in carbonate systems, the aqueous concentration of metal ions increases, e.g., dissolved amount of CaFe0.7Mg0.3(CO3)2 increases and reaches maximum of 1.78 mmol·L-1 and 8.26 mmol·L-1 at burial depth of 0.7 km with CO2 mole fraction of 0.1% and 10.0%, respectively. For the influence of overpressure in the calcite system, with overpressure ranging from 36 MPa to 83 MPa, pH reaches a minimum of 6.8 at overpressure of 51 MPa; meanwhile, Δ[Ca2+] increases slightly from -2.24 mmol·L-1 to -2.17 mmol·L-1 and remains negative, indicating it is also a precipitation process at burial depth of 3.9 km where overpressure generated. The method used in this study can be applied in assessing burial precipitation-dissolution processes and predicting possible pores in reservoirs with carbonate cement-water-carbon dioxide.

Structural and Thermodynamic Analysis of the First Mononuclear Aqueous Aluminum Citrate Complex Using DFT Calculations

Structural and thermodynamic properties of the mononuclear Al/citrate complexes have been theoretically investigated aiming to understand the coordination mechanism at an atomic level. GGA-DFT/PCM calculations have been performed for the different conformations and tautomers arising from the Al(3+) and citric acid (H3L) interaction in aqueous solution. The Gibbs reaction energies were estimated based on the reaction of the trigonal planar Al(OH)3 and H3L to form different Al-citrate complexes. The estimated Gibbs free reaction energies for the [AlL], [AlHL](+), and [Al(OH)L](-) species are in good agreement with the experimental values. In these species, the Al(3+) center is coordinated by two carboxylic and the tertiary hydroxyl groups of the citrate. Conversely to what has been proposed based on the experiments, the present theoretical calculations indicate that the citric acid hydroxyl group remains protonated upon the coordination of Al(3+). In fact, our model turns out to be more consistent with the relative pKa values of citrate protonation groups and with the hydrolysis constant of the H2O bound to Al(3+) leading to better agreement with the available experimental data.

Structures and formation mechanisms of aquo/hydroxo oligomeric beryllium in aqueous solution: a density functional theory study

The structures and formation mechanisms of a wide variety of aquo/hydroxo oligomeric beryllium clusters were investigated using density functional theory. The structural parameters of beryllium clusters were found to vary regularly with the stepwise substitution of bound water molecules in the inner coordination sphere by hydroxyl groups. According to the Gibbs free energies deduced from SMD solvation model computations, unhydrolyzed oligomeric beryllium species are the most favorable products of polymerization, independent of the degrees of hydrolysis of the reactants. Simulation of the formation processes of oligomeric beryllium showed that polymerization, in essence, involves the nucleophilic attack of a terminal hydroxyl group in one BeO4 tetrahedron on the beryllium center in another BeO4 tetrahedron, leading to the bridging of two BeO4 tetrahedrons by a hydroxyl group.

Embedded-cluster calculations in a numeric atomic orbital density-functional theory framework

We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding quantum and molecular mechanical (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level density functional theory calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO2(110).

An effective approach for modifying carbonaceous materials with niobium single sites to improve their catalytic properties

In this paper we show a very simple route for the incorporation of catalytically active niobium species on the surface of carbon materials, such as graphene oxide, carbon nanotubes and activated carbon.

Predicting pKa in Implicit Solvents: Current Status and Future Directions

Computational prediction of condensed phase acidity is a topic of much interest in the field today. We introduce the methods available for predicting gas phase acidity and pKas in aqueous and non-aqueous solvents including high-level electronic structure methods, empirical linear free energy relationships (LFERs), implicit solvent methods, explicit solvent statistical free energy methods, and hybrid implicit–explicit approaches. The focus of this paper is on implicit solvent methods, and we review recent developments including new electronic structure methods, cluster-continuum schemes for calculating ionic solvation free energies, as well as address issues relating to the choice of proton solvation free energy to use with implicit solvation models, and whether thermodynamic cycles are necessary for the computation of pKas. A comparison of the scope and accuracy of implicit solvent methods with ab initio molecular dynamics free energy methods is also presented. The present status of the theory and future directions are outlined.

Antioxidant activity of (+)-bergenin: a phytoconstituent isolated from the bark of Sacoglottis uchi Huber (Humireaceae)

(+)-Bergenin (1) was isolated from Sacoglottis uchi, a species of vegetable found in the Amazon region and popularly used for the treatment of several hepatic problems. The structure of 1 was fully characterized using IR, GC-MS and NMR (1D and 2D) analyses. This phytoconstituent has been used as an oriental folk medicine for the treatment of many diseases and shows antihepatotoxic properties. Tests with beta-carotene, DPPH and a heterogeneous Fenton system were carried out, confirming the antioxidant activity of 1. Theoretical calculations were performed to investigate the formation of the radical derivatives of 1 using *H, *OH, *CH(3), and *CCl(3) as initiator radicals. DFT thermodynamic calculations showed that the methoxyl group (O-6-CH(3)) is the most favorable site for radical attack. Frontier molecular orbital analysis showed that nucleophilic radical attack is favored on the aromatic ring of 1 where the LUMO is localized, with antibonding character with respect to the O-6-CH(3) bond. The possibilities of attack at other sites on 1 were investigated in detail in order to understand the regiospecificity of this reaction.