Experimental verification of ion impact angle distribution at divertor surfaces using micro-engineered targets on DiMES at DIII-D
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
SEARCH · Search NASA
Results for “Ion impact angle”
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.
Data-driven surrogate modeling of hPIC ion energy-angle distributions for high-dimensional sensitivity analysis of plasma parameters' uncertainty
In this work, we present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Evaluation of silicon carbide as a divertor armor material in DIII-D H-mode discharges
Silicon carbide (SiC) represents a promising but largely untested plasma-facing material (PFM) for next-step fusion devices. In this work, an analytic mixed-material erosion model is developed by calculating the physical (via SDTrimSP) and chemical (via empirical scalings) sputtering yield from SiC, Si, and C. The Si content in the near-surface SiC layer is predicted to increase during D plasma bombardment due to more efficient physical and chemical sputtering of C relative to Si. Silicon erosion from SiC thereby occurs primarily from sputtering of the enriched Si layer, rather than directly from the SiC itself. SiC coatings on ATJ graphite, manufactured via chemical vapor deposition, were exposed to repeated H-mode plasma discharges in the DIII-D tokamak to test this model. The qualitative trends from analytic modeling are reproduced by the experimental measurements, obtained via spectroscopic inference using the S/XB method. Quantitatively the model slightly under-predicts measured erosion rates, which is attributed to uncertainties in the ion impact angle distribution, as well as the effect of edge-localized modes. After exposure, minimal changes to the macroscopic or microscopic surface morphology of the SiC coatings were observed. Compositional analysis reveals Si enrichment of about 10%, in line with expectations from the erosion model. Extrapolating to a DEMO-type device, an order-of-magnitude decrease in impurity sourcing, and up to a factor of 2 decrease in impurity radiation, is expected with SiC walls, relative to graphite, if low C plasma impurity content can be achieved. Furthermore, these favorable erosion properties motivate further investigations of SiC as a low-Z, non-metallic PFM.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Angle-resolved Auger electron spectra induced by neon ion impact on aluminum
Auger electron emission from aluminum bombarded with 1 to 5 keV neon ions was studied by angle-resolved electron spectroscopy. The position and shape of the spectral features depended on the incident ion energy, angle of ion incidence, and electron take-off angle with respect to the aluminum surface. These spectral dependencies were interpreted in terms of the Doppler shift given to the Auger electron velocity by the excited atom ejected into the vacuum. For oblique ion incidence it is concluded that a flux of high energy atoms are ejected in a direction close to the projection of the ion beam on the target surface. In addition, a new spectral feature was found and identified as due to Auger emission from excited neon in the aluminum matrix.
Pepper, S. V.↗
MPEX AI Digital Twins
All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Non-Equilibrium Processes, Stability, Design and Control of Pulsed Plasmas for Materials Processing
The energy and angle of impact of ions, interspersed with fluxes of neutral radicals, on the surface of wafers during plasma materials processing are critical parameters for the fabrication of microelectronic devices, nanostructure materials and biocompatible materials. This sensitivity extends to nearly all applications of plasmas to materials processing and surface modification. Low temperature plasmas for plasma materials processing have traditionally used continuously excited plasmas. However, pulsed plasma processing is an enabling technology for achieving finer resolution of microelectronics fabrication and for producing unique surface properties in a wide range of applications. This research project is addressing fundamental plasma physics issues in materials processing – the control of pulsed plasmas and the means for optimizing plasma produced reactant fluxes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Resonance in radio frequency sheath admittance and enhanced impurity emission near the ion cyclotron frequency
Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/ R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.
ICRH antenna↗
Energy-angle distribution of the ions in the RF sheath of ICRH antennas
Radio-frequency sheaths forming at the surface of ion cyclotron resonance heating (ICRH) actuators in fusion experiments are associated with enhanced impurity sputtering from ICRH plasma-facing components (PFCs). The minimization of impurity fluxes from the ICRH PFC is a vital task for the usage of ICRH systems. Capturing the ion kinetics is important for an accurate description of plasma–material interactions, because the ion dynamics plays a crucial role in RF sheaths. Here, we present a hybrid particle-in-cell (hPIC) model able to capture the kinetic behavior of the ions. We analyze the kinetic ion energy-angle distributions (IEADs) impacting the RF antenna and its dependence on different plasma and RF sheath parameters. In particular, the IEAD dependence on RF frequency and magnetic field alignment is investigated. Using hPIC, we simulated a case emulating the latest experimental campaign from JET. The simulation showed that under specific plasma and RF parameters, the kinetic motion of ions results in a cusp formation in the IEAD.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors
Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators
We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Electron impact excitation of Hg/+/
A crossed charged beam technique was employed in the investigation reported. The ion beam was crossed at right angles by a magnetically confined electron beam. Cross sections were determined as a function of electron impact energy. Considerable structure in the excitation cross section appears immediately above the threshold. It is pointed out that a primary application of cross sections is related to the calculation of plasma rate coefficients.
Crandall, D. H.↗
Quantitative investigation of surface structure and interatomic potential with impact-collision ion scattering spectroscopy
Helium ion beam interactions with materials have important implications for magnetic confinement fusion, material modification, and helium ion microscopy. These interactions depend on the precise physics of how helium ions channel into the materials, which can vary greatly based on the local crystalline orientation. Here, we performed a dedicated experiment to investigate helium ion channeling in a well-characterized tungsten single crystal. Time-of-flight impact-collision ion scattering spectroscopy was used to obtain multi-angle maps of the backscattering intensity for 3 keV He + → W(111). We found that the backscattering intensity profile arising from helium ion channeling could be well described by a shadow cone analysis. This analysis revealed that subsurface W atoms as deep as the ninth monolayer contributed to the backscattering intensity profile. Binary collision approximation simulations were performed with MARLOWE to model the experimental maps with sufficient accuracy to allow for quantitative comparisons using reliability factors. These quantitative comparisons were applied to investigate how the W lattice structure and He–W interatomic potential affect the multi-angle maps.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine
The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.
Artificial neural networks↗
Surface Energies and Structure of Salt–Brine Interfaces
Permeability of salt formations is controlled by the equilibrium between the salt–brine and salt–salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. In this work, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt–brine interface to provide insight into the stability of salt systems. Mixed NaCl–KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt–brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2–3 times (5.0 M systems). Size differences in Na+ and K+ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt–brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
New experimental approach to understanding the chemical reactivity of oxide surfaces
Metal oxides have been an attractive option for a range of applications, including hydrogen sensors, microelectronics, and catalysis, due to their reactivity and tunability. The properties of metal oxides can vary greatly on their precise surface structure; however, few surface science techniques can achieve atomistic-level determinations of surface structure, and fewer yet can do so for insulator surfaces. Low energy ion beam analysis offers a potential insulator-compatible solution to characterizing the surface structure of metal oxides. As a feasibility study, we apply low energy ion beam analysis to investigate the surface structure of a magnetite single crystal, Fe 3 O 4 (100). We obtain multi-angle maps using both forward-scattering low energy ion scattering (LEIS) and backscattering impact-collision ion scattering spectroscopy (ICISS). Both sets of experimental maps have intensity patterns that reflect the symmetries of the Fe 3 O 4 (100) surface structure. However, analytical interpretation of these intensity patterns to extract details of the surface structure is significantly more complex than previous LEIS and ICISS structural studies of one-component metal crystals, which had far more symmetries to exploit. To gain further insight into the surface structure, we model our experimental measurements with ion-trajectory tracing simulations using molecular dynamics. Our simulations provide a qualitative indication that our experimental measurements agree better with a subsurface cation vacancy model than a distorted bulk model.
36 MATERIALS SCIENCE↗
Proton and oxygen plasmas at Uranus
Despite the presence of several large, ice moons within the Uranian magnetosphere, the Voyager 2 spacecraft did not detect any heavy ion plasma. This paper estimates the heavy ion density that would be consistent with a heavy neutral torus formed by photosputtering, charged particle sputtering, and micrometeoroid impact vaporization of icy surfaces on the moons, taking into account the large 60 deg angle between the satellite orbit plane and the magnetic equator. The expected heavy ion density is unobservably small. The observed proton plasma of the inner magnetosphere can be maintained by ionization of the atomic hydrogen corona and by the ionospheric proton source driven by photoelectron escape, for a plasma residence time of about 30 days. These two proton sources are comparable near L = 5. The same 30-day residence time for energetic protons implies an upper limit of 0.00001 on the fraction of incident solar wind protons that enter the magnetosphere and penetrate to within L less than 6 while conserving their first adiabatic invariants.
Cheng, Andrew F.↗
Direct Write of 3D Nanoscale Mesh Objects with Platinum Precursor via Focused Helium Ion Beam Induced Deposition
The next generation optical, electronic, biological, and sensing devices as well as platforms will inevitably extend their architecture into the 3rd dimension to enhance functionality. In focused ion beam induced deposition (FIBID), a helium gas field ion source can be used with an organometallic precursor gas to fabricate nanoscale structures in 3D with high-precision and smaller critical dimensions than focused electron beam induced deposition (FEBID), traditional liquid metal source FIBID, or other additive manufacturing technology. In this work, we report the effect of beam current, dwell time, and pixel pitch on the resultant segment and angle growth for nanoscale 3D mesh objects. We note subtle beam heating effects, which impact the segment angle and the feature size. Additionally, we investigate the competition of material deposition and sputtering during the 3D FIBID process, with helium ion microscopy experiments and Monte Carlo simulations. Our results show complex 3D mesh structures measuring ~300 nm in the largest dimension, with individual features as small as 16 nm at full width half maximum (FWHM). These assemblies can be completed in minutes, with the underlying fabrication technology compatible with existing lithographic techniques, suggesting a higher-throughput pathway to integrating FIBID with established nanofabrication techniques.
36 MATERIALS SCIENCE↗
Comparison of Single-Ion Conducting Polymer Gel Electrolytes for Sodium, Potassium, and Calcium Batteries: Influence of Polymer Chemistry, Cation Identity, Charge Density, and Solvent on Conductivity
From the standpoint of material diversification and sustainability, the development of so-called “beyond lithium-ion” battery chemistries is important for the future of energy storage. Na, K, and Ca are promising as the basis for battery chemistries in that these elements are highly abundant. Here, a series of single-ion conducting polymer electrolytes (SIPEs) for Na, K, and Ca batteries are synthesized and investigated. The two classes of metal cation neutralized SIPEs compared are crosslinked poly(ethylene glycol) dimethacrylate-x-styrene sulfonate (PEGDMA-SS) and poly(tetrahydrofuran) diacrylate-x-4-styrenesulfonyl (trifluoromethylsulfonyl)imide (PTHFDA-STFSI); three cation types, three charge densities, and four swelling states are examined. The impact on conductivity of all of these parameters is studied, and in conjunction with small angle X-ray scattering (SAXS), it is found that promoting ion dissociation and preventing the formation of dense ionic aggregates facilitates ion transport. These results indicate many of the lessons learned from the Li SIPE literature can be translated to beyond Li chemistries. At 25 °C, the best performing Na/K and Ca exchanged polymers yield active cation conductivity on the order of 10 –4 S/cm and 10 –6 S/cm, respectively, for ethylene carbonate:propylene carbonate gelled SIPEs, and 10 –5 S/cm and 10 –7 S/cm, respectively, for glyme gelled SIPEs.
X-ray scattering↗