Ion Association in AlCl3 Aqueous Solutions from Constrained First-Principles Molecular Dynamics

Speciation of La3+-Cl- Complexes in Hydrothermal Fluids from Deep Potential Molecular Dynamics

The stability of rare earth element (REE) complexes plays a crucial role in quantitatively assessing their hydrothermal migration and transformation. However, reliable data are lacking under high-temperature hydrothermal conditions, which hampers our understanding of the association behavior of REE. Here a deep learning potential model for the LaCl3-H2O system in hydrothermal fluids is developed based on the first-principles density functional theory calculations. The model accurately predicts the radial distribution functions compared to ab initio molecular dynamics (AIMD) simulations. Furthermore, species of La-Cl complexes, the dissociation pathway of the La-Cl complexes dissociation process, and the potential of mean forces and corresponding association constants (logK) for LaCln3-n (n = 1-4) are extensively investigated under a wide range of temperatures and pressures. Empirical density models for logK calculation are fitted with these data and can accurately predict logK data from both experimental results and AIMD simulations. The distribution of La-Cl species is also evaluated across a wide range of temperatures, pressures, and initial chloride concentration conditions. The results show that La-Cl complexes are prone to forming in a low-density solution, and the number of bonded Cl- ions increases with rising temperature. In contrast, in a high-density solution, La3+ dominates and becomes the more prevalent species.

Probing aqueous ions with non-local Auger relaxation

Non-local analogues of Auger decay are increasingly recognized as important relaxation processes in the condensed phase. Here, we explore non-local autoionization, specifically Intermolecular Coulombic Decay (ICD), of a series of aqueous-phase isoelectronic cations following 1s core-level ionization. In particular, we focus on Na+, Mg2+, and Al3+ ions. We unambiguously identify the ICD contribution to the K-edge Auger spectrum. The different strength of the ion-water interactions is manifested by varying intensities of the respective signals: the ICD signal intensity is greatest for the Al3+ case, weaker for Mg2+, and absent for weakly-solvent-bound Na+. With the assistance of ab initio calculations and molecular dynamics simulations, we provide a microscopic understanding of the non-local decay processes. We assign the ICD signals to decay processes ending in two-hole states, delocalized between the central ion and neighbouring water. Importantly, these processes are shown to be highly selective with respect to the promoted water solvent ionization channels. Furthermore, using a core-hole-clock analysis, the associated ICD timescales are estimated to be around 76 fs for Mg2+ and 34 fs for Al3+. Building on these results, we argue that Auger and ICD spectroscopy represents a unique tool for the exploration of intra- and inter-molecular structure in the liquid phase, simultaneously providing both structural and electronic information.

High-Rate Aqueous Aluminum-Ion Batteries Enabled by Confined Iodine Conversion Chemistry

Most reported cathode materials for rechargeable aqueous Al metal batteries are based on an intercalative-type chemistry mechanism. Herein, iodine embedded in MOF-derived N-doped microporous carbon polyhedrons (I₂ @ZIF-8-C) is proposed to be a conversion-type cathode material for aqueous aluminum-ion batteries based on "water-in-salt" electrolytes. Compared with the conventional Al-I₂ battery using ionic liquid electrolyte, the proposed aqueous Al-I₂ battery delivers much enhanced electrochemical performance in terms of specific capacity and voltage plateaus. Benefitting from the confined liquid-solid conversion of iodine in hierarchical N-doped microporous carbon polyhedrons and enhanced reaction kinetics of aqueous electrolytes, the I₂ @ZIF-8-C electrode delivers high reversibility, superior specific capacity (≈219.8 mAh g^-1 at 2 A g^-1 ), and high rate performance (≈102.6 mAh g^-1 at 8 A g^-1 ). The reversible reaction between I₂ and I^- , with I₃ ^- and I₅ ^- as intermediates, is confirmed via ex situ Raman spectra and X-ray photoelectron spectroscopy. Furthermore, solid-state hydrogel electrolyte is employed to fabricate a flexible Al-I₂ battery, which shows performance comparable to batteries using liquid electrolyte and can be integrated to power wearable devices as a reliable energy supply.

Cation Effect of Chloride Salting Agents on Transition Metal Ion Hydration and Solvent Extraction by the Basic Extractant Methyltrioctylammonium Chloride

The addition of a nonextractable salt has an important influence on the solvent extraction of metal ions, but the underlying principles are not completely understood yet. However, relating solute hydration mechanisms to solvent extraction equilibria is key to understanding the mechanism of solvent extraction of metal ions as a whole. We have studied the speciation of Co(II), Zn(II), and Cu(II) in aqueous solutions containing different chloride salts to understand their extraction to the basic extractant methyltrioctylammonium chloride (TOMAC). This includes the first speciation profile of Zn(II) in chloride media with the three Zn(II) species [Zn(H₂O)₆]²⁺, [ZnCl₃H₂O]⁻, and [ZnCl₄]^2–. The observed differences in extraction efficiency for a given transition metal ion can be explained by transition metal ion hydration due to ion–solvent interactions, rather than by ion–solute interactions or by differences in speciation. Chloride salting agents bearing a cation with a larger hydration Gibbs free energy reduce the free water content more, resulting in a lower hydration for the transition metal ion. This destabilizes the transition metal chloro complex in the aqueous phase and increases the extraction efficiency. Salting agents with di- and trivalent cations reduce the transition metal chloro complex hydration less than expected, resulting in a lower extraction efficiency. The cations of these salting agents have a very large hydration Gibbs free energy, but the overall hydration of these salts is reduced due to significant salt ion pair formation. The general order of salting-out strength for the extraction of metal ions from chloride salt solutions is Cs⁺ < Rb⁺ < NH₄⁺ ≈ K⁺ < Al³⁺ ≈ Mg²⁺ ≈ Ca²⁺ ≈ Na⁺ < Li⁺. These findings can help in predicting the optimal conditions for metal separation by solvent extraction and also contribute to a broader understanding of the effects of dissolved salts on solutes.

Addition of a nonextractable salt influences the stability and solvent extraction efficiency of metal complexes. Cations of different chloride salts reduce the solution free water content as a function of their increasing hydration energy and decreasing tendency for ion pair formation with chloride anions. These ion−solvent interactions reduce the hydration of metal complexes, increasing their distribution ratios. These effects influence aqueous transition metal complexes more than direct ion−solute interactions and changes in complex speciation.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Insights into water-mediated ion clustering in aqueous CaSO4 solutions: pre-nucleation cluster characteristics studied by ab initio calculations and molecular dynamics simulations

The molecular structure of growth units building crystals is a fundamental issue in the crystallization processes from aqueous solutions. In this work, a systematic investigation of pre-nucleation clusters and their hydration characteristics in aqueous CaSO₄ solutions was performed using ab initio calculations and molecular dynamics (MD) simulations. The results of ab initio calculations and MD simulations indicate that the dominant species in aqueous CaSO₄ solutions are monodentate ion-associated structures. Compared with charged ion clusters, neutral clusters are more likely to be present in an aqueous CaSO₄ solution. Neutral (CaSO₄)_m clusters are probably the growth units involved in the pre-nucleation or crystallization processes. Meanwhile, hydration behavior around ion associated species in aqueous CaSO₄ solutions plays an important role in related phase/polymorphism selections. Upon ion clustering, the residence of some water molecules around Ca²⁺ in ion-associated species is weakened while that of some bridging waters is enhanced due to dual interaction by Ca²⁺ and SO₄^2-. Some phase/polymorphism selections can be achieved in aqueous CaSO₄ solutions by controlling the hydration around pre-nucleation clusters. Moreover, the association trend between calcium and sulfate is found to be relatively strong, which hints at the low solubility of calcium sulfate in water.

Exploring Ion-Ion Interactions in Aqueous Solutions by a Combination of Molecular Dynamics and Neutron Scattering

Recent advances in computational and experimental techniques have allowed for accurate description of ion pairing in aqueous solutions. Free energy methods based on ab initio molecular dynamics, as well as on force fields accounting effectively for electronic polarization, can provide quantitative information about the structures and occurrences of individual types of ion pairs. When properly benchmarked against electronic structure calculations for model systems and against structural experiments, in particular neutron scattering, such force field simulations represent a powerful tool for elucidating interactions of salt ions in complex biological aqueous environments.

Coordination and hydrolysis of plutonium ions in aqueous solution using Car-Parrinello molecular dynamics free energy simulations

Car-Parrinello molecular dynamics (CPMD) simulations have been used to examine the hydration structures, coordination energetics, and the first hydrolysis constants of Pu(3+), Pu(4+), PuO2(+), and PuO2(2+) ions in aqueous solution at 300 K. The coordination numbers and structural properties of the first shell of these ions are in good agreement with available experimental estimates. The hexavalent PuO2(2+) species is coordinated to five aquo ligands while the pentavalent PuO2(+) complex is coordinated to four aquo ligands. The Pu(3+) and Pu(4+) ions are both coordinated to eight water molecules. The first hydrolysis constants obtained for Pu(3+) and PuO2(2+) are 6.65 and 5.70, respectively, all within 0.3 pH unit of the experimental values (6.90 and 5.50, respectively). The hydrolysis constant of Pu(4+), 0.17, disagrees with the value of -0.60 in the most recent update of the Nuclear Energy Agency Thermochemical Database (NEA-TDB) but supports recent experimental findings. The hydrolysis constant of PuO2(+), 9.51, supports the experimental results of Bennett et al. [Radiochim. Acta 1992, 56, 15]. A correlation between the pKa of the first hydrolysis reaction and the effective charge of the plutonium center was found.