On the Importance of Configuration Search to the Predictivity of Lanthanide Selectivity
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Engineering topics
Yang, Ping
Publications and source records attributed to Yang, Ping.
Advancing Rare-Earth (4 f ) and Actinide (5 f ) Separation through Machine Learning and Automated High-Throughput Experiments
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Light Transfers Through a Koch Shape Cloud
Abstract Modeling radiative transfer in a 3D cloudy atmosphere is critical to climate projections. A recently developed fast 3D radiation parameterization scheme gains some success in quantifying horizontal radiative transfer through cloud sides using cloud area fraction. Based on 3D Monte Carlo simulations of radiative transfer through an idealized single‐layer cloud with Koch‐shaped fractal geometry edges, here we show that radiative energy transport through cloud sides correlates more significantly with cloud area fraction than with cloud perimeter length. The results exemplify the importance of accounting for the horizontal radiative energy exchanges between cloud‐free and cloudy regions with cloud area fraction. Results from additional sensitivity simulations show that increased cloud vertical extent often enhances cloud‐side sunlight leak more significantly than cloud‐side sunlight interception. At low sun elevations, cloud‐side sunlight interception is enhanced more than cloud‐side sunlight leak does with the increase of cloud mass.
58 GEOSCIENCES↗
Efficient Parameterization of Density Functional Tight-Binding for 5 f -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 ThO 2 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 ThO 2 . 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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations
Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.
36 MATERIALS SCIENCE↗
KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models
Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Insights into the Mechanism of Neptunium Oxidation to the Heptavalent State
Abstract 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)O 4 (OH) 2 ] 3− during ozonolysis of aqueous M OH ( 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)O 2 (OH) 4 ] 3− with the hydroxide radical to form [Np(VI)O 2 (OH) 4 ] 2− , followed by an additional oxidation with HO⋅ to create [Np(VII)O 4 (OH) 2 ] 3− .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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. Here, 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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Unveiling Hidden Shake-Up Features in the Uranyl M 4 -Edge Spectrum
The M 4,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 M 4 -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 M 4,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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Advances in heavy alkaline earth chemistry provide insight into complexation of weakly polarizing Ra 2+ , Ba 2+ , and Sr 2+ cations
Numerous technologies—with catalytic, therapeutic, and diagnostic applications—would benefit from improved chelation strategies for heavy alkaline earth elements: Ra 2+ , Ba 2+ , and Sr 2+ . Unfortunately, chelating these metals is challenging because of their large size and weak polarizing power. We found 18-crown-6-tetracarboxylic acid (H 4 COCO) bound Ra 2+ , Ba 2+ , and Sr 2+ to form M(H x COCO) x–2 . Upon isolating radioactive 223 Ra from its parent radionuclides ( 227 Ac and 227 Th), 223 Ra 2+ reacted with the fully deprotonated COCO 4- chelator to generate Ra(COCO) 2- (aq) (log K Ra(COCO)2- = 5.97 ± 0.01), a rare example of a molecular radium complex. Comparative analyses with Sr 2+ and Ba 2+ 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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Linear graphlet models for accurate and interpretable cheminformatics
The surprising effectiveness of topology in the chemical sciences: graphlets in our open-source library, , provide accurate white-box 2D chemical property prediction.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Employing a template synthesis to access diastereopure Np( iv ) and U( iv ) complexes and analysis of their 5f orbitals in bonding
A template methodology was employed as an entryway into diastereopure U( iv ) and Np( iv ) organometallic complexes, and the 5f orbital contributions to bonding were determined.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Sensitivity of Arctic Surface Temperature to Including a Comprehensive Ocean Interior Reflectance to the Ocean Surface Albedo Within the Fully Coupled CESM2
Abstract Almost all current climate models simplify the ocean surface albedo (OSA) by assuming the reflected solar energy without the ocean interior contribution. In this study, an improved ocean surface albedo scheme is incorporated into the Community Earth System Model version 2 (CESM2) to assess the sensitivity of Arctic surface temperature to including ocean interior reflectance to the OSA. Fully coupled CESM2 simulations with and without ocean interior reflectance are subsequently performed, we focus on the analysis of Arctic surface temperature responses. Incorporating ocean interior reflectance increases absorbed solar radiation and warms the ocean, enhancing seasonal heat storage and release across the Arctic Ocean, and increasing sea ice reduction and positive climate feedbacks that elevates Arctic surface temperature. Seasonal variations in air‐surface temperature differences induce changes in turbulent heat flux patterns, concurrently modifying dynamic advection and moisture processes that affect boundary layer humidity and low clouds, especially in winter. Based on partitioning physical processes in the thermodynamic energy equation, surface air warming is induced primarily through positive heating anomalies of vertical advection, latent heat release, and longwave radiative forcing. Through an examination of the surface energy budget, skin temperature warming is driven predominantly by increased downward longwave radiation, positive surface albedo feedback in summer, and increased conductive heat transport from the ocean particularly in winter. Significant effects of ocean interior reflectance on the Arctic Ocean, including sea surface warming and sea ice reduction, justify the importance of ocean interior reflectance in climate models for better understanding of ongoing Arctic climate changes.
Meteorology & Atmospheric Sciences↗
A Refined Understanding of the Ice Cloud Longwave Scattering Effects in Climate Model
Abstract Because longwave (LW) absorption by greenhouse gases and clouds is more significant than the LW scattering effect by clouds, most climate models neglect cloud LW scattering to save computational costs. Ignoring cloud LW scattering directly overestimates outgoing longwave radiation (OLR). This study included ice‐cloud LW scattering treatment in the Exascale Energy Earth System Model (E3SM) version 2 and ran fully‐coupled simulations, prescribed sea surface temperature simulations, and offline radiative transfer calculations to comprehensively assess the impact of ice‐cloud LW scattering on global climate simulation. The instantaneous effect due to ice‐cloud LW scattering reduces the OLR by ∼1 W/m 2 on the global average and 2 W/m 2 on the tropical average. Tropospheric warming and high cloud amount reduction act to partially compensate for such instantaneous OLR reduction caused by the inclusion of LW scattering. When the simulation reaches the equilibrium, the surface warms by 0.66 K on average with respect to the simulation without LW scattering, with the Arctic surface temperature differences more than twice as large as that of the global mean. The impact of including LW scattering on the simulated climate change in response to 4 × CO 2 is also assessed. While including the cloud LW scattering does not significantly modify radiative forcing and total radiative feedback under such a scenario, it results in a 10% more positive cloud feedback.
54 ENVIRONMENTAL SCIENCES↗
Lower-length-scale modeling of chemical additions, corrosion and fission products in select MSR base salts
Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.
36 MATERIALS SCIENCE↗
Ab-Initio molecular dynamics simulations of binary NaCl-ThCl 4 and ternary NaCl-ThCl 4 -UCl 3 molten salts
Molten salt reactor (MSR) with Thorium (Th) fuel cycle has attracted growing attention due to its merits such as safety, low radioactive waste production, and non-proliferation. The design and safe operation of MSR rely on a thorough understanding of the thermophysical properties of molten salts over a wide range of composition and temperature. However, experimental data of Th-based molten salts are limited due to the inherent challenges of dealing with corrosive molten salts and the radioactive nature of actinides. Here, in this work, thermophysical properties including density, heat capacity, thermal expansion coefficient, and mixing energy of binary NaCl-ThCl 4 and ternary NaCl-ThCl 4 -UCl 3 molten salts are explored using ab-initio molecular dynamic simulations (AIMD) with the dDsC dispersion correction. The calculated mixing energy of binary NaCl-ThCl 4 exhibits a minimum close to the eutectic composition. The heat capacity of the mixtures is linearly dependent on the mole fraction of each component. 7-fold and 6-fold coordinated Th complexes are dominant in the mixtures. 8-fold coordinated Th complex increases with ThCl 4 fraction due to the formation of network structures. The average coordination number of Th exhibits a minimum near the eutectic composition. The minimum mixing energy for the ternary mixtures is observed in systems with a composition close to [NaCl] 0.5 [ThCl 4 ] 0.25 [UCl 3 ] 0.25 . The density positively deviates from the ideal solution in mixtures near this composition. These results are important to fill the data and knowledge gap of ThCl 4 molten salts.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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,Ad ArO) 3 N)U III (dme)] to generate a stable uranium(V) bis-alkoxide complex, namely, [(( Me,Ad ArO) 3 N)U V (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. Furthermore, this uranium(V) bis-alkoxide can then be reduced with KC 8 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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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. So 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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗