Efficient Parameterization of Density Functional Tight-Binding for 5f-Elements: A Th-O Case Study.

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th-O system, which was subsequently applied to the study of ThO2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO2. Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

Unveiling Hidden Shake-Up Features in the Uranyl M4-Edge Spectrum.

The M4,5-edge high energy resolution X-ray absorption near-edge structure (HR-XANES) spectra of actinyls offer valuable insights into the electronic structure and bonding properties of heavy-element complexes. To conduct a comprehensive spectral analysis, it is essential to employ computational methods that accurately account for relativistic effects and electron correlation. In this work, we utilize variational relativistic multireference configurational interaction methods to compute and analyze the X-ray M4-edge absorption spectrum of uranyl. By employing these advanced computational techniques, we achieve excellent agreement between the calculated spectral features and experimental observations. Moreover, the calculations unveil significant shake-up features, which arise from the intricate interplay between strongly correlated 3d core-electron and ligand excitations. This research provides important theoretical insights into the spectral characteristics of heavy-element complexes. Furthermore, it establishes the foundation for utilizing M4,5-edge spectroscopy as a means to investigate the chemical activities of such complexes. By leveraging this technique, we can gain a deeper understanding of the bonding behavior and reactivity of heavy-element compounds.

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman-Rubin Convergence Diagnostics.

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. In the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman-Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

Insights into the Mechanism of Neptunium Oxidation to the Heptavalent State.

Neptunium can exist in multiple oxidation states, including the rare and poorly understood heptavalent form. In this work, we monitored the formation of heptavalent neptunium [Np(VII)O4(OH)2]3- during ozonolysis of aqueous MOH (M=Li, Na, K) solutions using a combined experimental and theoretical approach. All experimental reactions were closely monitored via absorption and vibrational spectroscopy to follow both the oxidation state and the speciation of neptunium guided by the calculated vibrational frequencies for various neptunium species. The mechanism of the reaction partly involves oxidative dissolution of transient Np(VI) oxide/hydroxide solid phases, the identity of which are dependent on the co-precipitating counter-cation Li+/Na+/K+. Additional calculations suggest that the most favorable energetic pathway occurs through the reaction of a [Np(V)O2(OH)4]3- with the hydroxide radical to form [Np(VI)O2(OH)4]2-, followed by an additional oxidation with HO⋅ to create [Np(VII)O4(OH)2]3-.

Advances in heavy alkaline earth chemistry provide insight into complexation of weakly polarizing Ra2+, Ba2+, and Sr2+ cations.

Numerous technologies-with catalytic, therapeutic, and diagnostic applications-would benefit from improved chelation strategies for heavy alkaline earth elements: Ra2+, Ba2+, and Sr2+. Unfortunately, chelating these metals is challenging because of their large size and weak polarizing power. We found 18-crown-6-tetracarboxylic acid (H4COCO) bound Ra2+, Ba2+, and Sr2+ to form M(HxCOCO)x-2. Upon isolating radioactive 223Ra from its parent radionuclides (227Ac and 227Th), 223Ra2+ reacted with the fully deprotonated COCO4- chelator to generate Ra(COCO)2-(aq) (log KRa(COCO)2- = 5.97 ± 0.01), a rare example of a molecular radium complex. Comparative analyses with Sr2+ and Ba2+ congeners informed on what attributes engendered success in heavy alkaline earth complexation. Chelators with high negative charge [-4 for Ra(COCO)2-(aq)] and many donor atoms [≥11 in Ra(COCO)2-(aq)] provided a framework for stable complex formation. These conditions achieved steric saturation and overcame the weak polarization powers associated with these large dicationic metals.

Impact of Surface Oxidation on the Morphology of Uranium Dioxide Nanoparticles.

The undeniable importance of nanoparticles has led to vast efforts, in many fields of science, to understand their chemical and physical properties. In this paper, the morphology dependence of f-element nanoparticles is correlated to the oxygen environment and the type and coverage of capping ligands. This dependence was evaluated by first-principles calculations of the surface energies of different crystallographic planes (001, 110, and 111) as a function of the relative oxygen chemical potential and under the influence of different ligands. Uranium dioxide nanoparticles were the focus of this study due to their high sensitivity to oxidation compared to thorium dioxide nanoparticles, a homoleptic material but insensitive to oxidation. To fully explain the experimental observations of uranium dioxide nanocrystals, theoretical modeling shows that the consideration of surfaces with different oxidation conditions is necessary. It is shown that, for materials with low oxidation potential, such as uranium dioxide, the oxygen environment and capping ligand concentration are competing factors in determining the nanoparticle morphology.

Unusual Actinyl Complexes with a Redox-Active N,S-Donor Ligand.

Understanding the fundamental chemistry of soft N,S-donor ligands with actinides across the series is critical for separation science toward sustainable nuclear energy. This task is particularly challenging when the ligands are redox active. We herein report a series of actinyl complexes with a N,S-donor redox-active ligand that stabilizes different oxidation states across the actinide series. These complexes are isolated and characterized in the gas phase, along with high-level electronic structure studies. The redox-active N,S-donor ligand in the products, C5H4NS, acts as a monoanion in [UVIO2(C5H4NS-)]+ but as a neutral radical with unpaired electrons localized on the sulfur atom in [NpVO2(C5H4NS•)]+ and [PuVO2(C5H4NS•)]+, resulting in different oxidation states for uranium and transuranic elements. This is rationalized by considering the relative energy levels of actinyl(VI) 5f orbitals and S 3p lone pair orbitals of the C5H4NS- ligand and the cooperativity between An-N and An-S bonds that provides additional stability for the transuranic elements.

Structure and Dynamics of NaCl/KCl/CaCl2-EuCln (n = 2, 3) Molten Salts.

Modern molten salt reactor design and the techniques of electrorefining spent nuclear fuels require a better understanding of the chemical and physical behavior of lanthanide/actinide ions with different oxidation states dissolved in various solvent salts. The molecular structures and dynamics that are driven by the short-range interactions between solute cations and anions and long-range solute and solvent cations are still unclear. In order to study the structural change of solute cations caused by different solvent salts, we performed first-principles molecular dynamics simulations in molten salts and extended X-ray absorption fine structure (EXAFS) measurements for the cooled molten salt samples to identify the local coordination environment of Eu2+ and Eu3+ ions in CaCl2, NaCl, and KCl. The simulations reveal that with the increasing polarizing the outer sphere cations from K+ to Na+ to Ca2+, the coordination number (CN) of Cl- in the first solvation shell increases from 5.6 (Eu2+) and 5.9 (Eu3+) in KCl to 6.9 (Eu2+) and 7.0 (Eu3+) in CaCl2. This coordination change is validated by the EXAFS measurements, in which the CN of Cl- around Eu increases from 5 in KCl to 7 in CaCl2. Our simulation shows that the fewer Cl- ions coordinated to Eu leads to a more rigid first coordination shell with longer lifetime. Furthermore, the diffusivities of Eu2+/Eu3+ are related to the rigidity of their first coordination shell of Cl-: the more rigid the first coordination shell is, the slower the solute cations diffuse.

Architector for high-throughput cross-periodic table 3D complex building.

Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector, a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector, which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

Uranium-Mediated Peroxide Activation and a Precursor toward an Elusive Uranium cis-Dioxo Fleeting Intermediate.

The activation of chalcogen-chalcogen bonds using organometallic uranium complexes has been well documented for S-S, Se-Se, and Te-Te bonds. In stark contrast, reports concerning the ability of a uranium complex to activate the O-O bond of an organic peroxide are exceedingly rare. Herein, we describe the peroxide O-O bond cleavage of 9,10-diphenylanthracene-9,10-endoperoxide in nonaqueous media, mediated by a uranium(III) precursor [((Me,AdArO)3N)UIII(dme)] to generate a stable uranium(V) bis-alkoxide complex, namely, [((Me,AdArO)3N)UV(DPAP)]. This reaction proceeds via an isolable, alkoxide-bridged diuranium(IV/IV) species, implying that the oxidative addition occurs in two sequential, single-electron oxidations of the metal center, including rebound of a terminal oxygen radical. This uranium(V) bis-alkoxide can then be reduced with KC8 to form a uranium(IV) complex, which upon exposure to UV light, in solution, releases 9,10-diphenylanthracene to generate a cyclic uranyl trimer through formal two-electron photooxidation. Analysis of the mechanism of this photochemical oxidation via density functional theory (DFT) calculations indicates that the formation of this uranyl trimer occurs through a fleeting uranium cis-dioxo intermediate. At room temperature, this cis-configured dioxo species rapidly isomerizes to a more stable trans configuration through the release of one of the alkoxide ligands from the complex, which then goes on to form the isolated uranyl trimer complex.

Influencing Bonding Interactions of the Neptunyl (V, VI) Cations with Electron-Donating and -Withdrawing Groups.

Neptunium makes up the largest percentage of minor actinides found in spent nuclear fuel, yet separations of this element have proven difficult due to its rich redox chemistry. Developing new reprocessing techniques should rely on understanding how to control the Np oxidation state and its interactions with different ligands. Designing new ligands for separations requires understanding how to properly tune a system toward a desired trait through functionalization. Emerging technologies for minor actinide separations focus on ligands containing carboxylate or pyridine functional groups, which are desirable due to their high degree of functionalization. Here, we use DFT calculations to study the interactions of carboxylate and polypyridine ligands with the neptunyl cation [Np(V/VI)O2]+/2+. A systematic study is performed by varying the electronic properties of the carboxylate and polypyridine ligands through the inclusion of different electron-withdrawing and electron-donating R groups. We focus on how these groups can affect geometric properties, electronic structure, and bonding characterization as a function of the metal oxidation state and ligand character and discuss how these factors can play a role in neptunium ligand design principles.

State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives.

Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au25(SR)18 - is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.

Advancing the Am Extractant Design through the Interplay among Planarity, Preorganization, and Substitution Effects.

Advancing the field of chemical separations is important for nearly every area of science and technology. Some of the most challenging separations are associated with the americium ion Am(III) for its extraction in the nuclear fuel cycle, 241Am production for industrial usage, and environmental cleanup efforts. Herein, we study a series of extractants, using first-principle calculations, to identify the electronic properties that preferentially influence Am(III) binding in separations. As the most used extractant family and because it affords a high degree of functionalization, the polypyridyl family of extractants is chosen to study the effects of the planarity of the structure, preorganization of coordinating atoms, and substitution of various functional groups. The actinyl ions are used as a structurally simplified surrogate model to quickly screen the most promising candidates that can separate these metal ions. The down-selected extractants are then tested for the Am(III)/Eu(III) system. Our results show that π interactions, especially those between the central terpyridine ring and Am(III), play a crucial role in separation. Adding an electron-donating group onto the terpyridine backbone increases the binding energies to Am(III) and stabilizes Am-terpyridine coordination. Increasing the planarity of the extractant increases the binding strength as well, although this effect is found to be rather weak. Preorganizing the coordinating atoms of an extractant to their binding configuration as in the bound metal complex speeds up the binding process and significantly improves the kinetics of the separation process. This conclusion is validated by the synthesized 1,2-dihydrodipyrido[4,3-b;5,6-b]acridine (13) extractant, a preorganized derivative of the terpyridine extractant, which we experimentally showed was four times more effective than terpyridine at separating Am3+ from Eu3+ (SFAm/Eu ∼ 23 ± 1).

Hydrazine Energy Storage: Displacing N2 H4 from the Metal Coordination Sphere.

Hydrogen carriers, such as hydrazine (N2 H4 ), may facilitate long duration energy storage, a vital component for resilient grids by enabling more renewable energy generation. Lanthanide coordination chemistry with N2 H4 as well as efforts to displace N2 H4 from the metal coordination sphere to develop an efficient catalytic production cycle were detailed. Modeling the equilibrium of different ligand coordination, it was predicted that strong sigma donor molecules would be required to displace N2 H4 . Monitoring competition experiments with nuclear magnetic resonance confirmed that trimethyl phosphine oxide, dimethylformamide, and dimethyl sulfoxide displaced N2 H4 in large or small lanthanide complexes.

Synthesis, solid-state, solution, and theoretical characterization of an "in-cage" scandium-NOTA complex.

Developing chelators that strongly and selectively bind rare-earth elements (Sc, Y, La, and lanthanides) represents a longstanding fundamental challenge in inorganic chemistry. Solving these challenges is becoming more important because of increasing use of rare-earth elements in numerous technologies, ranging from paramagnets to luminescent materials. Within this context, we interrogated the complexation chemistry of the scandium(III) (Sc3+) trication with the hexadentate 1,4,7-triazacyclononane-1,4,7-triacetic acid (H3NOTA) chelator. This H3NOTA chelator is often regarded as an underperformer for complexing Sc3+. A common assumption is that metalation does not fully encapsulate Sc3+ within the NOTA3- macrocycle, leaving Sc3+ on the periphery of the chelate and susceptible to demetalation. Herein, we developed a synthetic approach that contradicted those assumptions. We confirmed that our procedure forced Sc3+ into the NOTA3- binding pocket by using single crystal X-ray diffraction to determine the Na[Sc(NOTA)(OOCCH3)] structure. Density functional theory (DFT) and 45Sc nuclear magnetic resonance (NMR) spectroscopy showed Sc3+ encapsulation was retained when the crystals were dissolved. Solution-phase and DFT studies revealed that [Sc(NOTA)(OOCCH3)]1- could accommodate an additional H2O capping ligand. Thermodynamic properties associated with the Sc-OOCCH3 and Sc-H2O capping ligand interactions demonstrated that these capping ligands occupied critical roles in stabilizing the [Sc(NOTA)] chelation complex.

Spectroscopic and electrochemical characterization of a Pr4+ imidophosphorane complex and the redox chemistry of Nd3+ and Dy3+ complexes.

The molecular tetravalent oxidation state for praseodymium is observed in solution via oxidation of the anionic trivalent precursor [K][Pr3+(NP(1,2-bis-tBu-diamidoethane)(NEt2))4] (1-Pr(NP*)) with AgI at -35 °C. The Pr4+ complex is characterized in solution via cyclic voltammetry, UV-vis-NIR electronic absorption spectroscopy, and EPR spectroscopy. Electrochemical analyses of [K][Ln3+(NP(1,2-bis-tBu-diamidoethane)(NEt2))4] (Ln = Nd and Dy) by cyclic voltammetry are reported and, in conjunction with theoretical modeling of electronic structure and oxidation potential, are indicative of principal ligand oxidations in contrast to the metal-centered oxidation observed for 1-Pr(NP*). The identification of a tetravalent praseodymium complex in in situ UV-vis and EPR experiments is further validated by theoretical modeling of the redox chemistry and the UV-vis spectrum. The latter study was performed by extended multistate pair-density functional theory (XMS-PDFT) and implicates a multiconfigurational ground state for the tetravalent praseodymium complex.

Relativistic Effects in Modeling the Ligand K-Edge X-ray Absorption Near-Edge Structure of Uranium Complexes.

Accurate modeling of the complex electronic structure of actinide complexes requires full inclusion of relativistic effects. In this study, we examine the effect of explicit inclusion of spin-orbit coupling (SOC) versus scalar relativistic effects on the predicted spectra for heavy-element complexes. In this study, we employ a relativistic two-component Hamiltonian in the X2C form with all of the electrons in the system being considered explicitly to compare and contrast with previous studies that included the relativistic effects by means of relativistic effective core potentials (RECPs). A few uranium complexes are chosen as model systems. Comparison of the computed Cl K-edge X-ray absorption spectra with experimental data shows significantly improved agreement when a variational relativistic treatment of SOC is performed. In particular, we note the importance of SOC terms to obtain not only correct transition energies but also correct intensities for these heavy-element complexes because of the redistribution of ligand bonding character among the valence MOs. While RECPs generally agree well with all-electron scalar relativistic calculations, there are some differences in the predicted spectra of open-shell systems. These methods are still suitable for broad application to analyze the qualitative nature of transitions in X-ray absorption spectra, but caution is recommended for quantitative analysis, as SOC can be non-negligible for both open- and closed-shell heavy-element systems.

An Allyl Uranium(IV) Sandwich Complex: Are ϕ Bonding Interactions Possible?

A method to explore head-to-head ϕ back-bonding from uranium f-orbitals into allyl π* orbitals has been pursued. Anionic allyl groups were coordinated to uranium with tethered anilide ligands, then the products were investigated by using NMR spectroscopy, single-crystal XRD, and theoretical methods. The (allyl)silylanilide ligand, N-((dimethyl)prop-2-enylsilyl)-2,6-diisopropylaniline (LH), was used as either the fully protonated, singly deprotonated, or doubly deprotonated form, thereby highlighting the stability and versatility of the silylanilide motif. A free, neutral allyl group was observed in UI2 (L1)2 (1), which was synthesized by using the mono-deprotonated ligand [K][N-((dimethyl)prop-2-enyl)silyl)-2,6-diisopropylanilide] (L1). The desired homoleptic sandwich complex U[L2]2 (2) was prepared from all three ligand precursors, but the most consistent results came from using the dipotassium salt of the doubly deprotonated ligand [K]2 [N-((dimethyl)propenidesilyl)-2,6-diisopropylanilide] (L2). This allyl-based sandwich complex was studied by using theoretical techniques with supporting experimental spectroscopy to investigate the potential for phi (ϕ) back-bonding. The bonding between UIV and the allyl fragments is best described as ligand-to-metal electron donation from a two carbon fragment-localized electron density into empty f-orbitals.

N2-to-NH3 conversion by excess electrons trapped in point vacancies on 5f-element dioxide surfaces.

Ammonia (NH3) is one of the basic chemicals in artificial fertilizers and a promising carbon-free energy storage carrier. Its industrial synthesis is typically realized via the Haber-Bosch process using traditional iron-based catalysts. Developing advanced catalysts that can reduce the N2 activation barrier and make NH3 synthesis more efficient is a long-term goal in the field. Most heterogeneous catalysts for N2-to-NH3 conversion are multicomponent systems with singly dispersed metal clusters on supporting materials to activate N2 and H2 molecules. Herein, we report single-component heterogeneous catalysts based on 5f actinide dioxide surfaces (ThO2 and UO2) with oxygen vacancies for N2-to-NH3 conversion. The reaction cycle we propose is enabled by a dual-site mechanism, where N2 and H2 can be activated at different vacancy sites on the same surface; NH3 is subsequently formed by H- migration on the surface via associative pathways. Oxygen vacancies recover to their initial states after the release of two molecules of NH3, making it possible for the catalytic cycle to continue. Our work demonstrates the catalytic activities of oxygen vacancies on 5f actinide dioxide surfaces for N2 activation, which may inspire the search for highly efficient, single-component catalysts that are easy to synthesize and control for NH3 conversion.