Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multilayers”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Broadband 920-nm mirror thin film damage competition

The 2023 Laser Damage conference thin-film damage competition was devoted to a survey on the state-of-the-art broadband near-IR multilayer dielectric (MLD) mirrors designed for ultra-short pulsed laser applications. The requirements for the coatings were a minimum reflection of 99.5% at 45-deg incidence angle for S-polarization from 830 nm to 1010 nm and group delay dispersion (GDD) < ± 50fs 2 . The participants were allowed to select the coating materials, coating design, and coating deposition method. Samples were damage tested at a single testing facility to enable direct comparison among the participants using a 25 ± 5 fs optical parametric chirped-pulse amplification (OPCPA) laser system operating at 5 Hz. The testing results from this set of 37 samples showed that dense coatings by ion-beam sputtering (IBS), magnetron sputtering (MS), and electron-beam ion assisted deposition (e-beam IAD) exhibited highest damage initiation onset (laser-induced damage threshold or LIDT) while e-beam coatings were low performers. In addition, multilayer coatings using tantala and/or hafnia as high index materials were top performers. Furthermore, this competition included for the first time the measurement of the damage growth onset (laser-induced damage growth threshold or LDGT). This latter performance metric plays an important role in establishing the safe operational conditions for larger aperture ultrashort pulsed lasers. Information pertaining to the morphology of the damage sites and their evolution under subsequent exposure to different laser fluences leading to damage growth is presented. Finally, not all coating samples in the survey met the GDD requirements stated above and associated measurements are discussed in the context of the present and past thin-film damage competitions focused on similar broadband near-IR MLD coatings.

23 broadband high reflectors↗
DOE OSTI2025

Sieving Hydrogen Isotopes via Machine Learning Assisted Chemical Vapor Deposition (CVD) of High‐Quality Monolayer Hexagonal Boron Nitride (h‐BN) on Iron Foils

Atomically thin two-dimensional (2D) ceramics, such as monolayer hexagonal boron nitride (h-BN), present potential for disruptive advances in separations. However, sub-atomic scale separation of hydrogen isotopes (H + /D + ) require near pristine 2D material membranes, and scalable synthesis of such high-quality h-BN comparable to mechanically exfoliated crystals remains a significant challenge. Here, we report a scalable Fe-catalyzed chemical vapor deposition (CVD) process for bottom-up synthesis of large-area, high-quality monolayer h-BN films, overcoming key limitations of conventional ammonia-based routes. By leveraging mechanistic insights and higher CVD temperatures, we suppress multilayer formation and achieve uniform monolayer h-BN coverage on commercially available Fe foils. Machine learning enables systematic exploration of the complex, multi-dimensional CVD parameter space (growth time, temperature, precursor temperature, multilayer faction, coverage), providing data-driven approaches to visualize and identify process regimes facilitating predominantly monolayer h-BN growth with minimal secondary nuclei/ad-layers. The optimized Fe-catalyzed CVD h-BN membranes show high-quality as observed by proton/deuteron (H + /D + ) selectivity ≈8.45, approaching the highest quality benchmark of mechanically exfoliated h-BN (H + /D + selectivity ≈10) as well as significantly outperforming Cu-catalyzed CVD h-BN membranes (H + /D + selectivity ≈3.62, control selectivity ≈1.7). Our work provides a scalable cost-effective route for high-quality monolayer h-BN synthesis for sub-atomic scale separations (H + /D + ) and demonstrates the broader potential of machine learning-guided optimization of CVD for advancing synthesis of 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular Simulation of Functionalized Covalent Organic Framework Membranes for Inorganic Salt Separation

Covalent organic frameworks (COFs) enable molecular-level design of nanochannels for selective separation in pressure-driven membrane processes. Through variation of building blocks and, subsequently, the pore structure and chemistry, membrane performance can be tailored. This study employs nonequilibrium molecular dynamics simulations to theoretically demonstrate the tunable selectivity of COF membranes through a bottom-up functionalization approach. Water and salt transport are evaluated for six β-ketoenamine-linked COFs with varying multilayer thicknesses. For the thinnest multilayer (0.64–0.72 nm), all COFs exhibit low sodium sulfate (Na 2 SO 4 ) rejection (55–66%). However, 20 stacked sheets (6.4–7.2 nm) provide 67–98% rejection, with the sulfonated COF providing the highest Na 2 SO 4 rejection. Analysis of time-resolved ion density profiles reveals that solute rejection is primarily governed by interfacial exclusion arising from pore size and functional group chemistry. Although increasing salt rejection compromises water permeance, the permeance of all COF membranes is at least two orders of magnitude greater than that of a commercially available nanofiltration membrane. Overall, this work guides the rational design of COF membranes for aqueous salt separation.

Nanofiltration↗
DOE OSTI2024

Size Effect of Negative Capacitance State and Subthreshold Swing in Van der Waals Ferrielectric Field–Effect Transistors

Analytical calculations corroborated by the finite element modeling show that thin films of Van der Waals ferrielectrics covered by a 2D-semiconductor are promising candidates for the controllable reduction of the dielectric layer capacitance due to the negative capacitance (NC) effect emerging in the thin films. The NC state is conditioned by energy-degenerated poly-domain states of the ferrielectric polarization induced in the films under incomplete screening conditions in the presence of a dielectric layer. Calculations performed for the FET-type heterostructure “ferrielectric CuInP 2 S 6 film—2D-MoS 2 single-layer—SiO 2 dielectric layer” reveal the pronounced size effect of the multilayer capacitance. Derived analytical expressions for the electric polarization and multilayer capacitance allow to predict the thickness range of the dielectric layer and ferrielectric film for which the NC effect is the most pronounced in various Van der Waals ferrielectrics, and the corresponding subthreshold swing becomes much less than the Boltzmann's limit. Obtained results can be useful for the size and temperature control of the NC effect in the steep-slope ferrielectric FETs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Adsorption of hydroxamic acid ligands for improved extraction of rare earth elements from monazite ores

Efficient separation of rare earth element (REE) ores via froth flotation requires the development of novel ligands with enhanced adsorption capacity and selectivity. To realize these advances, understanding the mechanisms underlying interactions between the ligand and mineral surfaces is essential. This study systematically evaluates the adsorption behavior of alkyl- and aromatic alkyl-substituted hydroxamic acid ligands on monazite surfaces using complementary spectroscopic techniques, including UV–visible (UV–vis) spectroscopy, Raman spectroscopy, infrared spectroscopy, and vibrational sum frequency generation (SFG) spectroscopy, together with the ab initio molecular dynamics (AIMD) simulations. Among the studied ligands, octanohydroxamic acid (OHA) and 4-ethoxy-N,2-dihydroxybenzamide (EDHBA) exhibit high adsorption capacity under basic pH (8–10) by forming multilayers on the surface. OHA has a higher equilibrium adsorption capacity compared to EDHBA, but it forms a less stable multilayer susceptible to disruption in the presence of interfering ions. AIMD results show that OHA adopts a single stable chelating geometry, while EDHBA exhibits multiple binding modes involving distinct interactions with La surface atoms and phosphate-bound oxygens, resulting in more complex adsorption kinetics. The variations in surface binding and intermolecular interactions observed between alkyl and aromatic molecules influence the differences in adsorption kinetics, equilibrium adsorption capacities on the mineral surface, and their flotation performance. This work provides valuable insight into the adsorption mechanism of ligands at mineral interfaces, which is crucial for guiding the design of new ligands with enhanced separation performance.

Zhou, Muchu [ORNL] (ORCID:0000000182650215)↗
DOE OSTI2025

First-Principles-Based Study of the Decomposition of Phenol and Hydroquinone on Pt(111) Combined with Quantitative Information from XPS Spectra to Address the Impact of Coverage and Number of Hydroxyl Functional Groups

A combined first-principles-based and experimental X-ray photoelectron spectroscopy approach was used to investigate the thermal decomposition of two model biofuel compounds, phenol and hydroquinone, on Pt(111) at both low and high coverages. The DFT-based approach yields adsorption geometries and energies, activation barriers and core-level binding energy shifts for C 1s and O 1s. Increasing the coverage in the theoretical model leads to slight shifts in core-level binding energies─toward higher values for C 1s and lower values for O 1s. It also alters the energy profiles of the decomposition reaction pathway, resulting in weaker adsorption energies and changes in both reaction and activation barriers. At low temperatures, we observe a multilayer for phenol and hydroquinone upon adsorption, with desorption occurring at 200 and 270 K, respectively. Following desorption of the multilayer, decomposition proceeds via initial O–H bond scission, followed by two parallel pathways involving either C–H or C–C bond scission, whereby in the case of phenol C–H bond scission occurs first. Here, we further provide characteristic core level binding energies by theoretical calculations that are subsequently used in experimental analyses, establishing a reference database for key spectra of phenolic functionalities applicable to a range of catalytic reactions.

09 BIOMASS FUELS↗
DOE OSTI2026

Syngas Production at Si Hybrid Photoelectrodes Modified with Re(I) and Mn(I) Tricarbonyl Phenanthroline Complexes Containing Reactive Aryl Azide Groups

We installed molecular CO 2 reduction (CO 2 R) catalysts directly onto Si (photo)electrodes. The highly reactive M(5-azido-1,10-phenanthroline)(CO) 3 X (where M = Mn or Re, X = Br or Cl) complexes readily bubbled when dissolved in polar organic solvents, in both the presence and absence of an ultraviolet light source. When placed on hydrogen-terminated Si (H-Si) and native silicon oxide (SiOx), similar amounts of the complex were attached to the surface under illumination (367 nm, 50–200 mW/cm 2 ) or in the dark. Surprisingly, these films revealed submonolayer coverages instead of the multilayered structures we expected. DFT analyses support monolayer formation, showing that the triplet-state nitrene of the complex is more energetically favorable than the singlet state. Using controlled-potential electrolysis experiments, we showed that Re- and Mn-containing films on pSi photoelectrodes generated small amounts of CO when exposed to 1 atm of CO 2 and 1 sun illumination. These amounts of CO were an order of magnitude greater than control surfaces, producing 5.59 × 10 –7 mol CO/h for Re(az-phen) and 7.83 × 10 –7 mol CO/h for Mn(az-phen) films. Much of the charge passed at the pSi electrodes was consumed by the competing hydrogen evolution reaction, which we attribute to the low molecular coverage and the presence of native oxide on the electrode surface after attachment. Finally, this work demonstrates the feasibility of reacting azide-containing ligands with Si surfaces. Still, it highlights the need for alternative ligand structures and reaction conditions to form multilayer films.

azides↗
DOE OSTI2025

Synthesis and Self-Assembly of Monodisperse Graphene Nanoribbons: Access to Submicron Architectures with Long-Range Order and Uniform Orientation

Fabricating organic semiconducting materials into large-scale well-organized architectures is critical for building high performance molecular electronics. While graphene nanoribbons (GNRs) hold enormous promise for various device applications, their assembly into a well-structured monolayer or multilayer architecture poses a substantial challenge. Here we report the preparation of length-defined monodisperse GNRs via the integrated iterative binomial synthesis (IIBS) strategy and their self-assembly into submicron-architectures with long-range order, uniform orientation as well as regular layers. Further, the use of short alkyl side chains benefits forming stable multi-layers through interlocking structures. By changing the length and backbone shapes of these monodisperse GNRs, various three-dimensional assemblies including multilayer stripes, monolayer stripes, and nanowires, can be achieved, leading to different photophysical properties and band gaps. The discovery of these intriguing self-assembly behaviors of length-defined GNRs is expected to open the door for various future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Curvature Induced Modifications of Chirality and Magnetic Configuration in Perpendicular Films

Designing curvature in three-dimensional (3D) magnetic nanostructures enables controlled manipulation of local energy landscapes, allowing for the modification of noncollinear spin textures relevant for next-generation spintronic devices. In this study, we experimentally investigate 3D magnetization textures in a Co/Pd multilayer film, exhibiting strong perpendicular magnetic anisotropy (PMA), deposited onto curved Cu nanowire meshes with diameters as small as 50 nm and lengths of several microns. Utilizing magnetic soft X-ray nanotomography, we achieve reconstructions of 3D magnetic domain patterns at approximately 30 nm spatial resolution. This approach provides detailed information on both the orientation and magnitude of magnetization within the film. Our results reveal that interfacial anisotropy in the Co/Pd multilayers drives the magnetization toward the local surface normal. In contrast to typical labyrinth domains observed in planar films, the presence of curved nanowires significantly alters the domain structure, with domains preferentially aligning along the nanowire axis in close proximity, while adopting random orientations farther away. We report direct experimental observation of a curvature-induced Dzyaloshinskii-Moriya interaction (DMI), which is quantified to be approximately one-third of the intrinsic DMI in Co/Pd stacks. The curvature induced DMI enhances stability of Néel-type domain walls. These experimental observations are further supported by micromagnetic simulations. Altogether, our findings demonstrate that introducing curvature into magnetic nanostructures provides a powerful strategy for tailoring complex magnetic behaviors, paving the way for the design of advanced 3D racetrack memory and neuromorphic computing devices.

Raftrey, David↗
DOE OSTI2025

Mo Atom Rearrangement Drives Layer-Dependent Reactivity in Two-Dimensional MoS 2

Two-dimensional (2D) materials offer a valuable platform for manipulating and studying chemical reactions at the atomic level, owing to the ease of controlling their microscopic structure at the nanometer scale. While extensive research has been conducted on the structure-dependent chemical activity of 2D materials, the influence of structural transformation during the reaction has remained largely unexplored. In this work, we report the layer-dependent chemical reactivity of MoS 2 during a nitridation atomic substitution reaction and attribute it to the rearrangement of Mo atoms. Our results show that the chemical reactivity of MoS 2 decreases as the number of layers is reduced in the few-layer regime. In particular, monolayer MoS 2 exhibits significantly lower reactivity compared with its few-layer and multilayer counterparts. Atomic-resolution transmission electron microscopy (TEM) reveals that MoN nanonetworks form as reaction products from monolayer and bilayer MoS 2 , with the continuity of the MoN crystals increasing with layer number, consistent with the local conductivity mapping data. The layer-dependent reactivity is attributed to the relative stability of the hypothetically formed MoN phase, which retains the number of Mo atomic layers present in the precursor. Specifically, the low chemical reactivity of monolayer MoS 2 is attributed to the high energy cost associated with Mo atom diffusion and migration necessary to form multilayer Mo lattices in the thermodynamically stable MoN phase. In conclusion, this study underscores the critical role of lattice rearrangement in governing chemical reactivity and highlights the potential of 2D materials as versatile platforms for advancing the understanding of materials chemistry at the atomic scale.

Chemical reactivity↗
DOE OSTI2024

THz generation by exchange-coupled spintronic emitters

Abstract The mechanism of THz generation in ferromagnet/metal (F/M) bilayers has been typically ascribed to the inverse spin Hall effect (ISHE). Here, we fabricated Pt/Fe/Cr/Fe/Pt multilayers containing two back-to-back spintronic THz emitters separated by a thin ( t Cr ≤ 3nm) wedge-shaped Cr spacer. In such an arrangement, magnetization alignment of the two Fe films can be controlled by the interplay between Cr-mediated interlayer exchange coupling (IEC) and an external magnetic field. This in turn results in a strong variation of the THz amplitude A , with A ↑ ↓ reaching up to 14 times A ↑ ↑ (arrows indicate the relative alignment of the magnetization of the two magnetic layers). This observed functionality is ascribed to the interference of THz transients generated by two closely spaced THz emitters. Moreover, the magnetic field dependence A ( H ) shows a strong asymmetry that points to an additional performance modulation of the THz emitter via IEC and multilayer design.

Adam, Roman↗
DOE OSTI2025

Physical patterning of high-Q superconducting niobium resonators via ion beam etching

The development of superconducting quantum circuits increasingly involves the exploration of chemically distinct materials and complex multilayered structures. Accelerating this trend may benefit from low-damage, materials-agnostic patterning techniques that are compatible with a broad range of materials. Here, in this work, we investigate the utility of low-energy ion beam etching (IBE), a physical patterning technique, as an alternative to reactive ion etching for fabricating low-loss superconducting resonators. We use niobium (Nb) resonators as a test platform, leveraging their well-characterized performance metrics for benchmarking. To address IBE-induced surface redeposition, we introduce an in situ aluminum capping layer combined with targeted post-fabrication chemical treatment. This strategy yields resonators with internal quality factors as high as 6 × 10 5 in the single-photon regime at 50 mK. These results establish low-energy IBE as a promising patterning technique for superconducting devices, with the potential to accelerate development across chemically diverse and multilayered material platforms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Photo-induced non-collinear interlayer RKKY coupling in bulk Rashba semiconductors

Abstract The interplay between light-matter, spin-orbit, and magnetic interactions allows the investigation of light-induced magnetic phenomena that are otherwise absent without irradiation. We present our analysis of light-driven effects on the interlayer exchange coupling mediated by a bulk Rashba semiconductor in a magnetic multilayer. The collinear magnetic exchange coupling mediated by the photon-dressed spin-orbit coupled electrons of BiTeI develops light-induced oscillation periods and displays new decay power laws, both of which are enhanced with an increasing light-matter coupling. For magnetic layers with non-collinear magnetization, we find a non-collinear magnetic exchange coupling uniquely generated by light-driving of the multilayer. As the non-collinear magnetic exchange coupling mediated by the photon-dressed electrons of BiTeI is unique to the irradiated system and it is enhanced with increasing light-matter coupling, this effect offers a promising platform of investigation of light-driven effects on magnetic phenomena in spin-orbit coupled systems. In this platform, light properties, such as its intensity, can serve as external knobs for inducing non-collinear couplings of the interlayer exchange and for modulating the collinear couplings. Both of these effects signify the photo-generated modification in the spin textures of spin-orbit coupled electrons in BiTeI.

Physics↗
DOE OSTI2026

Chiral superconductivity from a parent Chern band and its non-Abelian generalization

Here, we propose a minimal model starting from a parent Chern band with quartic dispersion that can describe the spin-valley polarized electrons in rhombohedral tetralayer graphene. The interplay between repulsive and attractive interactions on top of that parent Chern band is studied. We conduct standard self-consistent mean-field calculations, and find a rich phase diagram that consists of metal, quantum anomalous Hall crystal, chiral topological superconductor, as well as trivial gapped Bose-Einstein condensate. In particular, there exists a topological phase transition from the chiral superconductor to the Bose-Einstein condensate at zero temperature. Motivated by the recent experimental and theoretical studies of composite Fermi liquid in rhombohedral stacked multilayer graphene, we further generalize the physical electron model to its composite fermion counterpart based on a field theory analysis. The chiral superconductor phase of the composite fermion becomes the non-abelian Moore-Read quantum Hall phase. We argue that a chiral (pseudo-)spin liquid phase can emerge in the vicinity of this Moore-Read quantum Hall phase. Our work suggests rhombohedral multilayer graphene as a potential platform for rich correlated topological phases.

Wang, Yan-Qi [University of Maryland, College Park↗
DOE OSTI2024

Bilayer Crystals of Trapped Ions for Quantum Information Processing

Trapped-ion systems are a leading platform for quantum information processing, but they are currently limited to 1D and 2D arrays, which imposes restrictions on both their scalability and their range of applications. Here, we propose a path to overcome this limitation by demonstrating that Penning traps can be used to realize remarkably clean bilayer crystals, wherein hundreds of ions self-organize into two well-defined layers. These bilayer crystals are made possible by the inclusion of an anharmonic trapping potential, which is readily implementable with current technology. We study the normal modes of this system and discover salient differences compared to the modes of single-plane crystals. The bilayer geometry and the unique properties of the normal modes open new opportunities—in particular, in quantum sensing and quantum simulation—that are not straightforward in single-plane crystals. Furthermore, we illustrate that it may be possible to extend the ideas presented here to realize multilayer crystals with more than two layers. Our work increases the dimensionality of trapped-ion systems by efficiently utilizing all three spatial dimensions, and it lays the foundation for a new generation of quantum information processing experiments with multilayer 3D crystals of trapped ions. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Impact of high-temperature annealing on hafnia-silica composite coatings deposited via ion beam sputtering for high-peak power 1064 nm lasers

The maximum power handling fluence of high-peak and average power laser systems is often limited by the laser damage of the coatings on optical components. Furthermore, these multilayer dielectric coatings are limited in their maximum power handling due to laser-damage-prone defects in the lower optical bandgap, higher optical index material. Some of these defects can be mitigated by thermal annealing to high temperatures, which can greatly reduce the linear absorbative precursors. Typically, hafnia and silica are the materials of choice for high-peak and average power laser systems in the ultraviolet through infrared spectral range; however, hafnia crystalizes readily when annealed at high temperatures. In this study, we prepare composite HfO 2 -SiO 2 coatings by co-sputtering hafnia and silica in an ion beam sputtering system and compare them to pure hafnia-based coatings. We demonstrate that crystallinity in hafnia can be completely suppressed when it is mixed with silica, such as the composite coatings in this study. High reflectors were fabricated and annealed, demonstrating that the multilayer dielectric stacks can survive high-temperature annealing and exhibit an excellent linear absorption of 0.2 +/- 1 ppm at 1064 nm. Short- and long-pulse laser damage was explored, demonstrating the complex relationship between linear absorption and the non-linear absorption which drives pulsed laser damage. These results provide an excellent route to the creation of very low linear absorption optical coatings, which also utilize low scattering materials that are best suited for high-peak and average power applications.

Harthcock, Colin [Lawrence Livermore National Labo↗
DOE OSTI2026

Influence of Oxygen Flow and Stoichiometry on Optical Properties and Damage Resistance of Hafnium Oxide Thin Films

Hafnium oxide (HfO 2 ) is predominantly used as a high-index material in multi-layer dielectric coatings for high-peak- and high-average-power lasers, but laser damage often initiates within the HfO 2 layers despite their wide bandgap. Oxygen deficiency during deposition can introduce vacancy-related sub-bandgap states and absorptive defects, lowering damage resistance. This study investigates how oxygen flow during HfO 2 deposition with ion beam sputtering (IBS) affects its stoichiometry, defect formation, and nanosecond laser-induced damage threshold (LIDT) and whether single-layer trends predict multilayer performance. Single layers were deposited at varying oxygen flows, characterized for optical and structural properties, and tested for the LIDT at 1064 nm and 355 nm. Increasing oxygen flow drove the layer toward near-stoichiometric HfO 2 , reduced the refractive index, and altered the density of surface pinhole-like features. The single-layer LIDT at 355 nm increased with oxygen, whereas the 1064 nm LIDT was comparatively less sensitive to oxygen flow, consistent with the wavelength-dependent roles of absorptive precursors and microstructural defects. In contrast, a HfO 2 -based high-reflector (HR) showed a higher LIDT at lower oxygen flow, indicating that the family of damage precursors changes between single layers and multilayers; in stacks, structural properties such as stress, gas entrapment and thermal dissipation may outweigh the isolated absorptive defects found in single layers. These results demonstrate that the optimal oxygen flow condition depends on both LIDT wavelength and film architecture. We identified, for single layers, a 15–35 sccm window for maximizing the 1064 nm LIDT and a high-flow optimum (45 sccm) for the 355 nm LIDT and, for 355 nm HR stacks, a distinct lower-flow regime (~10 sccm).

Optics and optical instruments↗
DOE OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗