Density-Functional Theory Study of Iron(III) Hydrolysis in Aqueous Solution

Characterization of iron(III) in aqueous and alkaline environments with ab initio and ReaxFF potentials

The iron(III) complexes [Fe(H2O)n(OH)m]3-m (n + m = 5, 6, m ≤ 3) and corresponding proton transfer reactions are studied with total energy calculations, the nudged elastic band (NEB) method, and molecular dynamics (MD) simulations using ab initio and a modification of reactive force field potentials, the ReaxFF-AQ potentials, based on the implementation according to Böhm et al. [J. Phys. Chem. C 120, 10849-10856 (2016)]. Applying ab initio potentials, the energies for the reactions [Fe(H2O)n(OH)m]3-m + H2O → [Fe(H2O)n-1(OH)m+1]2-m + H3O+ in a gaseous environment are in good agreement with comparable theoretical results. In an aqueous (aq) or alkaline environment, with the aid of NEB computations, respective minimum energy paths with energy barriers of up to 14.6 kcal/mol and a collective transfer of protons are modeled. Within MD simulations at room temperature, a permanent transfer of protons around the iron(III) ion is observed. The information gained concerning the geometrical and energetic properties of water and the [Fe(H2O)n(OH)m]3-m complexes from the ab initio computations has been used as reference data to optimize parameters for the O-H-Fe interaction within the ReaxFF-AQ approach. For the optimized ReaxFF-AQ parameter set, the statistical properties of the basic water model, such as the radial distribution functions and the proton hopping functions, are evaluated. For the [Fe(H2O)n(OH)m]3-m complexes, it was found that while geometrical and energetic properties are in good agreement with the ab initio data for gaseous environment, the statistical properties as obtained from the MD simulations are only partly in accordance with the ab initio results for the iron(III) complexes in aqueous or alkaline environments.

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.

On The Nature of FeIV =Oaq in Aqueous Media: A DFT analysis

FeIV =Oaq is a key intermediate in many advanced oxidation processes and probably in biological systems. It is usually referred to as FeIV =O2+ . The pKa's of FeIV =Oaq as derived by DFT are: pKa1=2.37 M06 L/6-311++G(d,p) (SDD for Fe) and pKa2=7.79 M06 L/6-311++G(d,p) (SDD for Fe). This means that in neutral solutions, FeIV =Oaq is a mixture of (H2 O)4 (OH)FeIV =O+ and (H2 O)2 (OH)2 FeIV =O. The oxidation potential of FeIV =Oaq in an acidic solution, E0 {(H2 O)5 FeIV =O2+ /FeIII (H2 O)6 3+ , pH 0.0} is calculated with and without a second solvation sphere and the recommended value is between 2.86 V (B3LYP/Def2-TZVP, with a second solvation sphere) and 2.23 V (M06 L/Def2-TZVP without a second solvation sphere). This means that FeIV =Oaq is the strongest oxidizing agent formed in systems involving FeVI O4 2- even in neutral media.

Curious Behavior of Fe3+-As3+ Chemical Interactions and Nucleation of Clusters in Aqueous Medium

The simultaneous presence of Fe³⁺ and As³⁺ ions in groundwater (higher ppb or lower ppm level concentrations at circumneutral pH) as well as in acid mine drainages (AMDs)/industrial wastewater (up to few thousand ppm concentration at strongly acidic pH) are quite common. Therefore, understanding the chemical interactions prevalent between Fe³⁺ and As³⁺ ions in aqueous medium leading to nucleation of ionic clusters/solids, followed by aggregation and growth, is of great environmental significance. In the present work, we attempt to probe the nucleation process of Fe³⁺-As³⁺ clusters in solutions of various concentrations and pHs (from AMD to groundwater-like) using a combination of experimental and theoretical techniques. Interestingly, our study reveals nucleation of primary FeAs clusters in nearly all of them independent of concentration or pH. Theoretical studies employed density functional theory (DFT) to predict the primary clusters as stable Fe₄As₄ units. The surprising resemblance of these clusters with known Fe³⁺-As³⁺ minerals at the local level was observed experimentally, which provides an important clue about solid-phase growth from a range of Fe³⁺-As³⁺ solutions. Our experimental findings are further supported by a stepwise reaction mechanism established from detailed DFT studies.

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.

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.

Study on the grain size control of metatitanic acid in a mixture acid system based on Arrhenius and Boltzmann fitting

Herein, to control the particle size of metatitanic acid produced via titanium thermal hydrolysis in sulfuric–chloric mixture acid (SCMA) solutions, the relationship between its grain size and hydrolysis parameters is discussed, and the corresponding mathematical model was established using the experimental data. Firstly, Ti(OH)(SO₄)(Cl)(H₂O)₃ was selected as the most likely initial structure in the SCMA solution by comparing the experimental and corresponding simulated Raman spectra by density functional theory (DFT). Secondly, according to the predicted initial structure of TiO²⁺ and the experimental data for the hydrolysis process, with an increase in the concentration of TiO²⁺ and reaction temperature, the hydrolysis rate and grain size increased, while the agglomerate particle size decreased. Finally, a mathematic model was established and fitted by the Arrhenius equation and the Boltzmann distribution to describe the relationship between the grain size and hydrolysis parameters, as follows:

Herein, to control the particle size of metatitanic acid produced via titanium thermal hydrolysis in sulfuric–chloric acid (SCMA) solutions, the relationship between its grain size and hydrolysis parameters are discussed, and the corresponding methematical model is discussed using the experimental data.

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.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

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.