Engineering Papers⌕ Search

Engineering topics

Wang, Gaoxue

Publications and source records attributed to Wang, Gaoxue.

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↗

In Situ High-Temperature Raman Spectroscopy of UCl 3 : A Combined Experimental and Theoretical Study

Uranium trichloride (UCl 3 ) has received growing interest for its use in uranium-fueled molten salt reactors and in the pyrochemical processing of used fuel. In this paper, we report for the first time the experimentally determined Raman spectra of UCl 3 , at both ambient condition and in situ high temperatures up to 871 K. The frequencies of five of the Raman-active vibrational modes (v i ) of UCl 3 exhibit a negative temperature derivative ((∂ν i /∂T) P ) with increasing temperature. Here, this red-shift behavior is likely due to the elongation of U–Cl bonds. The average isobaric mode Grüneisen parameter (γ iP = 0.91 ± 0.02) of UCl 3 was determined through use of the coefficient of thermal expansion published in Vogel et al. (2021) and the (∂ νi /∂T) P values determined in this study. These results are in general agreement with those calculated here by density functional theory (DFT+U). Finally, a comparison of the ambient band positions of UCl 3 to those of isostructural lanthanide (La–Eu) and actinide chlorides (Am–Cf) has been made.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Work function lowering of LaB 6 by monolayer hexagonal boron nitride coating for improved photo- and thermionic-cathodes

We report a lowering of work function for lanthanum hexaboride (LaB 6 ) by monolayer hexagonal boron nitride (hBN) coating. Photoemission electron microcopy (PEEM) and thermionic emission electron microscopy (TEEM) both revealed that the hBN coated region of a LaB 6 (100) single crystal has a lower work function compared to the bare (i.e., non-coated) and graphene coated regions. A broad and uniform brighter image of the hBN coated region in PEEM was quantitatively supported by a 0.4 eV decrease in the work function in photoelectron spectra compared to the bare region. TEEM results were consistent in that the hBN coated region exhibited thermionic emission at 905 °C, whereas the bare and graphene coated regions did not. A larger decrease in the work function for hBN coated LaB 6 (100) compared to graphene coated LaB 6 (100) was qualitatively supported by our density functional theory calculations. Adding an oxide layer in the calculations improved consistency between the calculation and experimental results. We followed up our calculations with synchrotron-radiation x-ray photoelectron spectroscopy and confirmed the presence of an oxide layer on our LaB 6 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rugged bialkali photocathodes encapsulated with graphene and thin metal film

Protection of free-electron sources has been technically challenging due to lack of materials that transmit electrons while preventing corrosive gas molecules. Two-dimensional materials uniquely possess both of required properties. Here, we report three orders of magnitude increase in active pressure and factor of two enhancement in the lifetime of high quantum efficiency (QE) bialkali photocathodes (cesium potassium antimonide (CsK 2 Sb)) by encapsulating them in graphene and thin nickel (Ni) film. The photoelectrons were extracted through the graphene protection layer in a reflection mode, and we achieved QE of ~ 0.17% at ~ 3.4 eV, 1/e lifetime of 188 h with average current of 8.6 nA under continuous illumination, and no decrease of QE at the pressure of as high as ~ 1 × 10 –3 Pa. In comparison, the QE decreased drastically at 10 –6 Pa for bare, non-protected CsK 2 Sb photocathodes and their 1/e lifetime under continuous illumination was ~ 48 h. We attributed the improvements to the gas impermeability and photoelectron transparency of graphene.

36 MATERIALS SCIENCE↗

Near-Threshold Photoemission from Graphene-Coated Cu(110)

In this article, the brightness of electron beams emitted from photocathode sources plays a critical role in determining the performance of x-ray free-electron lasers and ultrafast electron-diffraction applications. In order to achieve the maximum brightness, the electrons need to be emitted from a photocathode with the lowest-possible mean transverse energy (MTE). Recent investigations have shown that capping a Cu(110) photocathode with a monolayer of graphene can protect the quantum efficiency (QE) from long-term exposure to varying vacuum conditions. However, there have been no studies that investigate the effects that a monolayer of graphene has on the MTE. Here, we report on measurements of a graphene-coated Cu(110) single crystal near the photoemission threshold for room and liquid-nitrogen temperatures. At room temperature, a minimum MTE of 25 meV is measured at 295 nm. At liquid-nitrogen temperatures, a minimum MTE of 9 meV is measured at the photoemission threshold of 290 nm.

36 MATERIALS SCIENCE↗

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

Ammonia (NH 3 ) 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 N 2 activation barrier and make NH 3 synthesis more efficient is a long-term goal in the field. Most heterogeneous catalysts for N 2 -to-NH 3 conversion are multicomponent systems with singly dispersed metal clusters on supporting materials to activate N 2 and H 2 molecules. Herein, we report single-component heterogeneous catalysts based on 5 f actinide dioxide surfaces (ThO 2 and UO 2 ) with oxygen vacancies for N 2 -to-NH 3 conversion. The reaction cycle we propose is enabled by a dual-site mechanism, where N 2 and H 2 can be activated at different vacancy sites on the same surface; NH 3 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 NH 3 , making it possible for the catalytic cycle to continue. Our work demonstrates the catalytic activities of oxygen vacancies on 5 f actinide dioxide surfaces for N 2 activation, which may inspire the search for highly efficient, single-component catalysts that are easy to synthesize and control for NH 3 conversion.

36 MATERIALS SCIENCE↗

Ab initio investigation of elastic properties of dilute Cu alloys for high-gradient accelerating structures

The elastic properties and solid solution strengthening of several Cu alloys that may be promising for high gradient applications were studied using ab initio calculations based on density functional theory. Specifically, the temperature and solute concentration dependence of elastic properties including bulk and shear moduli, thermal expansion coefficient, and heat capacity of these alloys were calculated. We found that the bulk and shear moduli are linearly dependent on the solute concentration in the dilute regime and change significantly for different Cu alloys. The thermal expansion coefficient and heat capacity are comparatively less sensitive to solute concentrations. Solid solution strengthening is observed to be linearly correlated with the change of lattice parameter of Cu. This work provides an in-depth understanding of the elastic properties of Cu alloys, which is key for predicting their performance at high gradient in normal conducting accelerating structures.

36 MATERIALS SCIENCE↗

Ab-initio Cu alloy design for high-gradient accelerating structures

Operation of normal conducting accelerator structures at high accelerating gradients is beneficial for many accelerator applications in basic science, industry, medicine, and National Security. RF breakdown is the major factor that limits the achievable accelerating gradients. Previous experiments on copper (Cu) have demonstrated that RF breakdown probability can be significantly decreased by hardening the material and alloying Cu with solutes such as silver (Ag). In this paper, we propose a figure-of-merit (FOM) that characterizes the ability of Cu alloys to withstand high-gradients. The FOM represents a trade-off between hardening through solid solution strengthening and the additional thermal stress induced by incremental RF pulse heating resulting from changes in electronic properties induced by alloying. We performed high-throughput ab initio calculations and computed the FOM for a large number of binary Cu alloys. Several promising candidate alloys for high-gradient accelerating structures were identified, such as CuAg, CuCd, CuHg, CuAu, CuIn, and CuMg. CuAg alloys have previously exhibited low RF breakdown rates in experiments. The results provide guidance for selecting alloys for the future high-gradient normal conducting accelerating structures operating at very high gradients.

36 MATERIALS SCIENCE↗