Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

AnB8-: borozene complexes with monovalent actinide

In light of recently reported monovalent lanthanide in borozene complexes LnB8- (Ln = La, Pr, Tb, Tm, Yb), the corresponding AnB8- (An = Ac, Pa, Bk, Md, No) actinide species within the same group were theoretically investigated in respect of oxidation state, stability, electronic structure and chemical bonding pattern. Our investigations reveal the feasibility of actinides, especially for the late actinide borozene compounds (BkB8-, MdB8-, NoB8-) adopting a monovalent oxidation state of +I, a phenomenon fine-tuned by the doubly aromatic borozene B82-. Early actinides (AcB8-, PaB8-) however exhibit a tendency towards higher trivalent oxidation states.

Low-Overpotential Rechargeable Na-CO2 Batteries Enabled by an Oxygen-Vacancy-Rich Cobalt Oxide Catalyst

Rechargeable sodium-carbon dioxide (Na-CO2) batteries have been proposed as a promising CO2 utilization technique, which could realize CO2 reduction and generate electricity at the same time. They suffer, however, from several daunting problems, including sluggish CO2 reduction and evolution kinetics, large polarization, and poor cycling stability. In this study, a rambutan-like Co3O4 hollow sphere catalyst with abundant oxygen vacancies was synthesized and employed as an air cathode for Na-CO2 batteries. Density functional theory calculations reveal that the abundant oxygen vacancies on Co3O4 possess superior CO2 binding capability, accelerating CO2 electroreduction, and thereby improving the discharge capacity. In addition, the oxygen vacancies also contribute to decrease the CO2 decomposition free energy barrier, which is beneficial for reducing the overpotential further and improving round-trip efficiency. Benefiting from the excellent catalytic ability of rambutan-like Co3O4 hollow spheres with abundant oxygen vacancies, the fabricated Na-CO2 batteries exhibit extraordinary electrochemical performance with a large discharge capacity of 8371.3 mA h g-1, a small overpotential of 1.53 V at a current density of 50 mA g-1, and good cycling stability over 85 cycles. These results provide new insights into the rational design of air cathode catalysts to accelerate practical applications of rechargeable Na-CO2 batteries and potentially Na-air batteries.

Ultrastrong Electron-Phonon Coupling in Uranium-Organic Frameworks Leading to Inverse Luminescence Temperature Dependence

Electron-phonon interactions, crucial in condensed matter, are rarely seen in Metal-Organic Frameworks (MOFs). Detecting these interactions typically involves analyzing luminescence in lanthanide- or actinide-based compounds. Prior studies on Ln- and Ac-based MOFs at high temperatures revealed additional peaks, but these were too faint for thorough analysis. In our research, we fabricated a high-quality, crystalline uranium-based MOF (KIT-U-1) thin film using a layer-by-layer method. Under UV light, this film showed two distinct "hot bands," indicating a strong electron-phonon interaction. At 77 K, these bands were absent, but at 300 K, a new emission band appeared with half the intensity of the main luminescence. Surprisingly, a second hot band emerged above 320 K, deviating from previous findings in rare-earth compounds. We conducted a detailed ab-initio analysis employing time-dependent density functional theory to understand this unusual behaviour and to identify the lattice vibration responsible for the strong electron-phonon coupling. The KIT-U-1 film's hot-band emission was then utilized to create a highly sensitive, single-compound optical thermometer. This underscores the potential of high-quality MOF thin films in exploiting the unique luminescence of lanthanides and actinides for advanced applications.

Dynamic dissolution of Cm3+ ions incorporated at the calcite-water interface: an ab initio molecular dynamics simulation study

The stability of actinide-mineral solid solution in a water environment is critical for assessing the safety of nuclear-waste geological repositories and studying actinide migration in natural systems. However, the dissolution behavior of actinide ions incorporated at the mineral-water interface is still unclear at the atomic level. Herein, we present metadynamics simulations of the reaction pathways, thermodynamics and kinetics of trivalent curium ions (Cm3+) dissolving from calcite surfaces. Cm3+ ions incorporated in different calcite surfaces (i.e., terrace and stepped surfaces) with distinct coordination environments have different reaction pathways, free energy barriers and free energy changes. We found that Cm dissolution from a stepped surface is more favorable than that from a terrace surface, both thermodynamically and kinetically. In addition, water molecules seem to promote the detachment of curium ions from the surface by exerting a pulling force via water coordination with Cm3+ and a pushing force via proton migration to the surface layer and water diffusion in the vacant Cm site. Thus, the findings from this work prove to be a milestone in revealing the dynamic dissolution mechanism of trivalent actinides from minerals and would also help predict the dissolution behaviors of other metal ions at the solid-water interface in chemical and environmental sciences.

Greater transferability and accuracy of norm-conserving pseudopotentials using nonlinear core corrections

We present an investigation into the transferability of pseudopotentials (PPs) with a nonlinear core correction (NLCC) using the Goedecker, Teter, and Hutter (GTH) protocol across a range of pure GGA, meta-GGA and hybrid functionals, and their impact on thermochemical and non-thermochemical properties. The GTH-NLCC PP for the PBE density functional demonstrates remarkable transferability to the PBE0 and ωB97X-V exchange-correlation functionals, and relative to no NLCC, improves agreement significantly for thermochemical benchmarks compared to all-electron calculations. On the other hand, the B97M-rV meta-GGA functional performs poorly with the PBE-derived GTH-NLCC PP, which is mitigated by reoptimizing the NLCC parameters for this specific functional. The findings reveal that atomization energies exhibit the greatest improvements from use of the NLCC, which thus provides an important correction needed for covalent interactions relevant to applications involving chemical reactivity. Finally we test the NLCC-GTH PPs when combined with medium-size TZV2P molecularly optimized (MOLOPT) basis sets which are typically utilized in condensed phase simulations, and show that they lead to consistently good results when compared to all-electron calculations for atomization energies, ionization potentials, barrier heights, and non-covalent interactions, but lead to somewhat larger errors for electron affinities.

A computational study on the coordination modes and electron absorption spectra of the complexes U(iv) with N,N,N',N'-tetramethyl-diglycolamide and anions

Lipophilic N,N,N',N'-tetraalkyl-diglycolamides (TRDGAs) are promising extractants for actinides separation in spent nuclear fuel reprocessing. Usually, in the extracted complexes of actinide and lanthanide ions of various oxidation states, the metal ions are completely surrounded by 2 or 3 TRDGA molecules, and the counter anions do not directly coordinate with them. In contrast, the extracted complexes of U(iv) from different media presenting different absorption spectra indicate that the anions (Cl- and NO3-) are directly involved in the coordination with U(iv) in the first inner sphere. Based on this exceptional observation in solvent extraction, taking the coordination of U(iv) with N,N,N',N'-tetramethyl-diglycolamide (TMDGA, the smallest analogue of TRDGA) as the research object, we mimic the behaviours of counterions (Cl- and NO3-) and the water molecule during coordination of TMDGA with U(iv), especially combining with the simulation of the absorption spectra. We demonstrate that during the complexing of TMDGA to U(iv), the counterion Cl- will occupy one coordination number in the inner coordination sphere, and NO3- will occupy two by bidentate type; however, the ubiquitous water cannot squeeze in the inner coordination sphere. In addition, the coordination of Cl- and NO3- is proved to favour the extraction with the lower binding energy. Moreover, the simulation of absorption spectra is in good agreement with the observation from experiments, further verifying the aforementioned conclusion. This work in some way will provide guidance to improve the computation methods in research of actinides by mimicking the absorption spectra of actinide ions in different complexes.

Enhanced Hydrogen Bonds of the (H2O)n (n = 4-8) Clusters Confined in Uranyl Peroxide Cluster Na20(UO2)20(O2)30

Water is a basic resource and an essential component of living organisms. It often exhibits some novel properties under confinement. The water clusters (H₂O)_n (n = 4-8) confined in the cavity of uranyl peroxide cluster Na₂₀(UO₂)₂₀(O₂)₃₀ (U₂₀) have been computationally investigated by using ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) calculations in this study. The results show that the confined water clusters can form hydrogen bonds with the internal oxygen atoms (O_uranyl) of U₂₀, and their conformations changed significantly. The average lengths (2.553-2.645 Å) of hydrogen bonds in confined (H₂O)_n are shorter than those (2.731-2.841 Å) in the corresponding free water clusters. Moreover, these confined hydrogen bonds show better hydrogen bond patterns according to the quantified indices. The natural bond orbital (NBO) calculations determine that there is electron transferring from the U₂₀ to its interior (H₂O)_n. It is the main reason for enhancing hydrogen bond interactions among the confined water molecules because their oxygen atoms are more negatively charged and their hydrogen atoms are more positively charged. The quantum theory of atoms in molecules (QTAIM) and interacting quantum atoms (IQA) analyses indicate that the confined hydrogen bonds are more covalent, based on the significant electron density ρ(r) and local energy density H(r) at the bond critical points (BCPs), and the stronger energies of interatomic exchange interactions (V_xc). These findings may help to promote the communication of confined water clusters and enrich the understating of confined hydrogen bonds.

The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory

The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. We evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol^-1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.

Ultrafiltration separation of Am(VI)-polyoxometalate from lanthanides

Partitioning of americium from lanthanides (Ln) present in used nuclear fuel plays a key role in the sustainable development of nuclear energy1-3. This task is extremely challenging because thermodynamically stable Am(III) and Ln(III) ions have nearly identical ionic radii and coordination chemistry. Oxidization of Am(III) to Am(VI) produces AmO22+ ions distinct with Ln(III) ions, which has the potential to facilitate separations in principle. However, the rapid reduction of Am(VI) back to Am(III) by radiolysis products and organic reagents required for the traditional separation protocols including solvent and solid extractions hampers practical redox-based separations. Herein, we report a nanoscale polyoxometalate (POM) cluster with a vacancy site compatible with the selective coordination of hexavalent actinides (238U, 237Np, 242Pu and 243Am) over trivalent lanthanides in nitric acid media. To our knowledge, this cluster is the most stable Am(VI) species in aqueous media observed so far. Ultrafiltration-based separation of nanoscale Am(VI)-POM clusters from hydrated lanthanide ions by commercially available, fine-pored membranes enables the development of a once-through americium/lanthanide separation strategy that is highly efficient and rapid, does not involve any organic components and requires minimal energy input.

Theoretical insights into the coordination structures, stabilities and electronic spectra of Cm3+ species at the gibbsite-water interface: A computational study combing ab initio molecular dynamics and wave function theory

The information of structure and stability of actinide species is key to understand the sorption mechanism of actinides at mineral-water interface. Such information is approximately derived from experimental spectroscopic measurements and needs to be accurately obtained by a direct atomic-scale modelling. Herein, systematic first-principles calculations and ab initio molecular dynamics (AIMD) simulations are carried out to study the coordination structures and absorption energies of Cm(III) surface complexes at gibbsite-water interface. Eleven representative complexing sites are investigated. The most stable Cm³⁺ sorption species are predicted to be a tridentate surface complex in weakly acidic/neutral solution condition and a bidentate one in the alkaline solution condition. Moreover, luminescence spectra of the Cm³⁺ aqua ion and the two surface complexes are predicted based on the high-accuracy ab initio wave function theory (WFT). The results give a gradually decreasing emission energy in good agreement with experimental observation of a red shift of peak maximum with pH increasing from 5 to 11. This work is a comprehensive computational study involving AIMD and ab initio WFT methods to gain the coordination structures, stabilities, and electronic spectra of actinide sorption species at the mineral-water interface, thus providing important theoretical support for geological disposal of actinide waste.

Norm-Conserving 4f-in-Core Pseudopotentials and Basis Sets Optimized for Trivalent Lanthanides (Ln = Ce-Lu)

We present here a set of scalar-relativistic norm-conserving 4f-in-core pseudopotentials, together with complementary valence-shell Gaussian basis sets, for the lanthanide (Ln) series (Ce-Lu). The Goedecker, Teter, and Hutter (GTH) formalism is adopted with the generalized gradient approximation (GGA) for the exchange-correlation Perdew-Burke-Ernzerhof (PBE) functional. The 4f-in-core pseudopotentials are built through attributing 4f-subconfiguration 4fⁿ (n = 1-14) for Ln (Ln = Ce-Lu) into the atomic core region, making it possible to circumvent the difficulty of the description of the open 4fⁿ valence shell. A wide variety of computational benchmarks and tests have been carried out on lanthanide systems including Ln³⁺-containing molecular complexes, aqueous solutions, and bulk solids to validate the accuracy, reliability, and efficiency of the optimized 4f-in-core GTH pseudopotentials and basis sets. The 4f-in-core GTH pseudopotentials successfully replicate the main features of lanthanide structural chemistry and reaction energetics, particularly for nonredox reactions. The chemical bonding features and solvation shells, hydrolysis energetics, acidity constants, and solid-state properties of selected lanthanide systems are also discussed in detail by utilizing these new 4f-in-core GTH pseudopotentials. This work bridges the idea of keeping highly localized 4f electrons in the atomic core and efficient pseudopotential formalism of GTH, thus providing a highly efficient approach for studying lanthanide chemistry in multi-scale modeling of constituent-wise and structurally complicated systems, including electronic structures of the condensed phase and first-principles molecular dynamics simulations.

Ultra-Efficient Americium/Lanthanide Separation through Oxidation State Control

Lanthanide/actinide separation is a worldwide challenge for atomic energy and nuclear waste treatment. Separation of americium (Am), a critical actinide element in the nuclear fuel cycle, from lanthanides (Ln) is highly desirable for minimizing the long-term radiotoxicity of nuclear waste, yet it is extremely challenging given the chemical similarity between trivalent Am(III) and Ln(III). Selective oxidation of Am(III) to a higher oxidation state (OS) could facilitate this separation, but so far, it is far from satisfactory for practical application as a result of the unstable nature of Am in a high OS. Herein, we find a novel strategy to generate stable pentavalent Am (Am(V)) through coordination of Am(III) with a diglycolamide ligand and oxidation with Bi(V) species in the presence of an organic solvent. This strategy leads to efficient stabilization of Am(V) and an extraordinarily high separation factor (>10⁴) of Am from Ln through one single contact in solvent extraction, thereby opening a new avenue to study the high-OS chemistry of Am and fulfill the crucial task of Ln/Am separation in the nuclear fuel cycle. The synergistic coordination and oxidation process is found to occur in the organic solvent, and the mechanism has been well elucidated by quantum-theoretical modeling.

Theoretical Prediction of Graphene-like 2D Uranyl Material with p-Orbital Antiferromagnetism

Versatile graphene-like two-dimensional materials with s-, p- and d-block elements have aroused significant interests because of their extensive applications while there is a lack of f-block one. Herein we report...

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

Actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical for the proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport properties. To tackle this problem, we adopted a machine-learning approach to study ThCl₄-NaCl and UCl₃-NaCl binary systems. The machine-learning potential with the density functional theory accuracy allows us to obtain long molecular dynamics trajectories (ns) for large systems (10³ atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperatures. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination sphere lifetime. Our study also suggests that actinides in molten salts may not follow well-known entropy-scaling laws.

Optimized Pseudopotentials and Basis Sets for Semiempirical Density Functional Theory for Electrocatalysis Applications

A complete set of Goedecker, Teter, and Hutter (GTH) norm-conserving pseudopotentials (PPs) have been optimized, in conjunction with molecular optimized (MOLOPT) basis sets, for both the B97M-rV and ωB97X-V density functionals for members of the main-group elements and 3d and noble metals. The resulting small-core PPs and corresponding DZVP, TZVP, and TZV2P basis sets yield notable improvements compared to the original PBE defaults when validated against all electron calculations for redox reaction energies and geometries, binding energies, and vibrational Stark effects for metal monocarbonyls in vacuum. Further validation of the optimized PP/MOLOPT basis set combinations was performed using ab initio molecular dynamics simulations and shows greatly improved agreement with experimental trends for metal surface relaxations and the adsorption behavior of CO on solid metal surfaces.