Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic fields”

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

Elastic cross section is entanglement entropy

We present universal relations between entanglement entropy, which quantifies the quantum correlation between subsystems, and the cross section, which is the primary observable for high-energy particle scattering, by employing a careful formulation of wave packets for the incoming particles. For 2-to-2 elastic scattering with no initial entanglement and subdividing the system along particle labels, we show that both the Rényi and Tsallis entropies in the final states are directly proportional to the elastic cross section in units of the transverse size for the initial wave packets, which is then interpreted as the elastic scattering probability. The relations do not depend on the underlying dynamics of the quantum field theory and are valid to all orders in coupling strengths. Furthermore, computing quantum correlations between momentum and nonkinematic data leads to entanglement entropies expressed as various semi-inclusive elastic cross sections. Our result gives rise to a novel “area law” for entanglement entropy in a two-body system. Published by the American Physical Society 2025

Low, Ian (ORCID:0000000275709597)↗

Entanglement, trace anomaly, and confinement in QCD

We formulate confinement in quantum chromodynamics (QCD) as an entropic surface phenomenon. Quark and gluon quantum information is localized on a transverse, entangling two-sphere of radius 𝑅 𝐸⁢𝐸 ; at this radius the QCD vacuum—partitioned by a hadron into interior and exterior regions—reaches its maximal entanglement entropy. Lattice-QCD determinations of the scalar (trace) gravitational form factors fix both 𝑅 𝐸⁢𝐸 and the transverse trace-anomaly density 𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ), yielding a parameter-free slope 𝑐 ℎ = 8⁢𝜋 2 ⁢𝑅𝑆$^2_{𝐸⁢𝐸}$𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ) and a mechanical entropy 𝑆 𝐸⁢𝐸 ⁡(𝑦) = 𝑐 ℎ ⁢𝑦 that grows linearly with rapidity 𝑦. The entropy gradient ∂ 𝑅 𝑆 𝐸⁢𝐸 changes sign at 𝑅 𝐸⁢𝐸 : it pushes colored degrees of freedom outward for 𝑟 <𝑅 𝐸⁢𝐸 and pulls them inward for 𝑟 >𝑅 𝐸⁢𝐸 , thereby localizing them on the codimension-2 entangling two-sphere Σ ⊥ = 𝑆$^2_{𝑅_{𝐸⁢𝐸}}$ (which, in the infinite-momentum frame (IMF), projects onto the transverse plane)—the “information wall.” This provides a high-energy (large-𝑦) entropic confinement diagnostic that complements—rather than replaces—Wilson’s area-law criterion, which probes long-distance dynamics near the rest frame (𝑦 → 0). Imposing unitarity on an entropic ansatz for the amplitude yields 𝜎⁡(𝑠) ∝ 𝑦 𝛿 . World data favor 𝛿 = 2 for elastic 𝑝⁡𝑝⁡($𝑝\bar{⁡𝑝}$) scattering and heavy-quark photoproduction, whereas 𝜙 photoproduction favors a softer 𝛿 = 0.387. All extracted cross sections remain well below the Froissart-Martin bound. These results provide a confinement criterion quantified directly from nonperturbative QCD inputs, unifying the trace anomaly, entanglement entropy, and high-energy scattering within a single quantitative framework.

Color confinement↗

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗

Regional Subsidence Analysis Through a Multi-Scale Modeling Framework Based on Breakage Mechanics

Although poroelastic models are often used to explain the delay between subsurface fluid depletion and ground subsidence, inelastic compaction involving permanent changes of rock microstructure may exacerbate hydro-mechanical coupling, thus influencing the interpretation of measurements and long-term forecasts. Here, a multi-scale modeling approach is discussed, which accounts for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. A constitutive model built within the framework of breakage mechanics is proposed to link the hydraulic conductivity of granular rocks with inelastic deformations and changes in grading caused by injection-depletion cycles at stress levels far from yielding. The proposed model has been incorporated into large-scale simulation frameworks, thus enabling the spatiotemporal mapping of regional subsidence through a hybrid, semi-analytical approach. Numerical results based on this strategy show that the model allows isolating near-field and far-field effects into the computation of land subsidence and can generate forecasts for different modeling scenarios (e.g., elastic and inelastic compaction, constant permeability, and concurrent change of compressibility and permeability). In particular, examples of simulations for the case of the Groningen gas field are discussed, showing the model capabilities to use both field measurements and laboratory tests for the generation of reasonable subsidence maps, without expensive computational costs. Results indicate that ignoring coupled inelastic effects has major consequences on the predicted timescale of subsidence. Specifically, while all the model scenarios produced similar long-term ground settlements, those ignoring breakage-dependent permeability changes result in a variation of the temporal window of residual subsidence of the order of several decades.

58 GEOSCIENCES↗

Higgs Portal Vector Dark Matter Interpretation: Review of Effective Field Theory Approach and Ultraviolet Complete Models

In this report a review of the Higgs portal-vector dark matter interpretation of the spin-independent dark matter nu-cleon elastic scattering cross section is presented, where the invisible Higgs decay width measured at the LHC is used. Effective Field Theory and ultraviolet complete models are discussed. LHC interpretations show only the scalar and Majorana dark matter scenarios; we propose including interpretation for vec-tor dark matter in the EFT and UV completion theoretical framework. In addition, our studies suggest an extension of the LHC dark matter interpretations to the sub-GeV regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stripe magnetic order and field-induced quantum criticality in the perfect triangular-lattice antiferromagnet CsCeSe 2

The two-dimensional triangular-lattice antiferromagnet (TLAF) is a textbook example of frustrated magnetic systems. Despite its simplicity, the TLAF model exhibits a highly rich and complex magnetic phase diagram, featuring numerous distinct ground states that can be stabilized through frustrated next-nearest-neighbor couplings or anisotropy. In this paper, we report low-temperature magnetic properties of the TLAF material CsCeSe 2 . The inelastic neutron scattering (INS) together with specific heat measurements and density functional theory calculations of crystalline electric field suggest that the ground state of Ce ions is a Kramers doublet with strong easy-plane anisotropy. Elastic neutron scattering measurements demonstrate the presence of stripe-yz magnetic order that develops below T N = 0.35K, with the zero-field ordered moment of m Ce ≈ 0.65μ B . Application of magnetic field first increases the ordering temperature by about 20% at the intermediate field region and eventually suppresses the stripe order in favor of the field-polarized ferromagnetic state via a continuous quantum phase transition (QPT). The field-induced response demonstrates sizable anisotropy for different in-plane directions, B ∥ a and B ⊥ ⁢a, which indicates the presence of bond-dependent coupling in the spin Hamiltonian. As a result, we further show theoretically that the presence of anisotropic bond-dependent interactions can change the universality class of QPT for B ∥ a and B ⊥ ⁢a.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-destructive depth-resolved characterization of residual strain fields in high electron mobility transistors using differential aperture x-ray microscopy

Localized residual stress and elastic strain concentrations in microelectronic devices often affect the electronic performance, resistance to thermomechanical damage, and, likely, radiation tolerance. A primary challenge for the characterization of these concentrations is that they exist over sub-μm length-scales, precluding their characterization by more traditional residual stress measurement techniques. Here, in this study, we demonstrate the use of synchrotron x-ray-based differential aperture x-ray microscopy (DAXM) as a viable, non-destructive means to characterize these stress and strain concentrations in a depth-resolved manner. DAXM is used to map two-dimensional strain fields between the source and the drain in a gallium nitride (GaN) layer within high electron mobility transistors (HEMTs) with sub-μm spatial resolution. Strain fields at various positions in both pristine and irradiated HEMT specimens are presented in addition to a preliminary stress analysis to estimate the distribution of various stress components within the GaN layer. γ-irradiation is found to significantly reduce the lattice plane spacing in the GaN along the sample normal direction, which is attributed to radiation damage in transistor components bonded to the GaN during irradiation.

42 ENGINEERING↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Electron beam irradiation effects on bulk metals: a comparative study of polycrystalline versus single-crystalline structures

This study investigates the effects of electron beam (e-beam) irradiation on the mechanical and structural properties of eight bulk metallic samples, comprising both polycrystalline (PC) and single-crystalline (SC) forms of Ni, Cr, V, and Ti. These metals were evaluated as potential candidates for beam exit windows in high-power (MW-class) particle accelerators. The primary objective is to identify metals capable of withstanding the conditions of high-power/MW-class e-beam accelerators and serve effectively as exit windows. Selection criteria were based on each metal’s intrinsic properties, power dissipation capability, and irradiation-induced changes in mechanical behavior, including hardness, elastic modulus, and defect density. Comprehensive characterization was conducted using field-emission scanning electron microscopy (FESEM), atomic force microscopy (AFM), scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM–EDS), X-ray diffraction (XRD), high-resolution transmission electron microscopy (HRTEM), selected area electron diffraction (SAED), and nanoindentation, performed both before and after exposure to a ~¿66 kGy dose from a 10 MeV e-beam accelerator. Results revealed that e-beam irradiation induced hardening in PC Ni, whereas PC Ti, commonly used in beam exit windows, exhibited softening. The observed softening in PC Ti is attributed to grain coarsening, elongation, and the formation of twins and twin boundaries, in contrast to the smaller, compressed grains in the pristine (Pr) PC Ti samples, consistent with the Hall–Petch relationship. The stresses due to twinning are small and insignificant in influencing the overall hardening of the PC Ti irradiated sample when compared to the stresses due to the dislocation density. Conversely, SC Ti samples exhibited irradiation-induced hardening. The SC Ti irradiated samples developed additional irradiation-induced modifications in crystallographic texture of (100), (101), (110), (200), (112), (004), and (211) as evidenced from the XRD results, which could probably explain the hardening effect that is caused by irradiation.

36 MATERIALS SCIENCE↗

Theory of spin waves in a hexagonal antiferromagnet

We construct a field-theoretic description of spin waves in hexagonal antiferromagnets with three magnetic sublattices and coplanar 120° magnetic order. The three Goldstone modes can be separated by point-group symmetry into a singlet α 0 and a doublet β. The α 0 singlet is described by the standard theory of a free relativistic scalar field. The field theory of the β doublet is analogous to the theory of elasticity of a two-dimensional isotropic solid with distinct longitudinal and transverse “speeds of sound.” The well-known Heisenberg models on the triangular and kagome lattices with nearest-neighbor exchange turn out to be special cases with accidental degeneracy of the spin-wave velocities. The speeds of sound can be readily calculated for any lattice model. We apply this approach to the compounds of the Mn 3 X family with stacked kagome layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Two-dimensional quantum universality in the spin-1/2 triangular-lattice quantum antiferromagnet Na 2 BaCo(PO 4 ) 2

An interplay of geometrical frustration and strong quantum fluctuations in a spin-1/2 triangular-lattice antiferromagnet (TAF) can lead to exotic quantum states. Here, we report the neutron-scattering, magnetization, specific heat, and magnetocaloric studies of the recently discovered spin-1/2 TAF Na 2 BaCo(PO 4 ) 2 , which can be described by a spin-1/2 easy axis XXZ model. The zero-field neutron diffraction experiment reveals an incommensurate antiferromagnetic ground state with a significantly reduced ordered moment of about 0.54(2) μ B /Co. Different magnetic phase diagrams with magnetic fields in the a b plane and along the easy c -axis were extracted based on the magnetic susceptibility, specific heat, and elastic neutron-scattering results. In addition, two-dimensional (2D) spin dispersion in the triangular plane was observed in the high-field polarized state, and microscopic exchange parameters of the spin Hamiltonian have been determined through the linear spin wave theory. Consistently, quantum critical behaviors with the universality class of d = 2 and ν z = 1 were established in the vicinity of the saturation field, where a Bose–Einstein condensation (BEC) of diluted magnons occurs. The newly discovered quantum criticality and fractional magnetization phase in this ideal spin-1/2 TAF present exciting opportunities for exploring exotic quantum phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗

Input to the European strategy for particle physics: strong-field quantum electrodynamics

This document sets out the intention of the strong-field QED community to carry out, both experimentally and numerically, high-statistics parametric studies of quantum electrodynamics in the non-perturbative regime, at fields approaching and exceeding the critical or ‘Schwinger’ field of QED (F qed ≈ m 2 c 3 / eh ≈ 1.3 x 10 18 V/m) in the rest frame of a charged particle. In this regime, several exotic and fascinating phenomena are predicted to occur that have never been directly observed in the laboratory. These include Breit–Wheeler pair production, vacuum birefringence, and quantum radiation reaction. This experimental programme will also serve as a stepping stone towards studies of elusive phenomena such as elastic scattering of real photons and the conjectured perturbative breakdown of QED at extreme fields. State-of-the-art high-power laser facilities in Europe and beyond are starting to offer unique opportunities to study this uncharted regime at the intensity frontier, which is highly relevant also for the design of future multi-TeV lepton colliders. A transition from qualitative observational experiments to quantitative and high-statistics measurements can only be performed with large-scale collaborations and with systematic experimental programmes devoted to the optimisation of several aspects of these complex experiments, including detector developments, stability and tolerances studies, and laser technology.

Sarri, G. [Queen’s University Belfast (United King↗

Scattering Observables from One- and Two-body Densities: Formalism and Application to $\pmb \gamma $ ${}^3\hbox {He}$ Scattering

We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Elastic scattering of 3 He + 4 He with the SONIK scattering chamber

Measurements of the elastic scattering cross section of 3 He and 4 He are important in order to improve constraints on theoretical models of 4 He ⁢( 3 He ,𝛾)⁢ 7 Be, a key reaction in Big Bang nucleosynthesis and solar neutrino production. The astrophysical 𝑆 factor for this reaction is a significant source of uncertainty in the standard-solar-model prediction of the 7 Be and 8 B solar neutrino fluxes. The elastic scattering measurements reported in the literature do not extend to low energies and lack proper uncertainty quantification. A new measurement of the 4 He ⁢( 3 He, 3 He)⁢ 4 He reaction has been made at center-of-mass energies 𝐸c.m.=0.38−3.13 MeV using the Scattering of Nuclei in Inverse Kinematics (SONIK) scattering chamber: a windowless, extended gas target surrounded by an array of 30 collimated silicon charged particle detectors situated at TRIUMF. This is the first elastic scattering measurement of 3 He + 4 He made below 500 keV and it has greater angular range and better precision than previous measurements. The elastic scattering data were analyzed using both 𝑅-matrix and halo effective field theory frameworks, and values of the 𝑠-wave scattering length and effective range were extracted. Here, the resulting improvement in knowledge of the 𝑠-wave effective-range function at low energies reduces the overall uncertainty in 𝑆 34 at solar energies.

Elastic scattering reactions↗

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗