Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PPPL”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Zero curvature is a necessary and sufficient condition for a spin-orbital decomposition

There has been an extended debate regarding the existence of a spin-orbital decomposition of the angular momentum of photons and other massless particles. It was recently shown that there are both geometric and topological obstructions preventing any such decomposition. Here we show that any geometric connection on a particle’s state space induces a splitting of the angular momentum into two operators. These operators are well-defined angular momentum operators if and only if the connection has zero curvature. Massive particles have two canonical curved connections corresponding to boosts and rotations, respectively. Furthermore, these can be uniquely combined to produce a flat connection, and this gives a novel derivation of the Newton-Wigner position operator and the corresponding spin and orbital angular momenta for relativistic massive particles. When the mass is taken to zero, transverse boosts and rotations degenerate, leaving only a single connection for massless particles. This connection produces a commonly proposed splitting of the massless angular momentum into two operators. However, the connection is not flat, explaining why these operators do not satisfy the angular momentum commutation relations and are thus not true spin and orbital angular momentum operators.

Angular momentum↗
DOE OSTI2024

Kinetic simulations underestimate the effects of waves during magnetic reconnection

Collisionless plasma systems are often studied using fully kinetic simulations, where protons and electrons are treated as particles. Due to their computational expense, it is necessary to reduce the ion-to-electron mass ratio $m_i/m_e$ or the ratio between plasma and cyclotron frequencies in simulations of large systems. In this Letter we show that when electron-scale waves are present in larger-scale systems, numerical parameters affect their amplitudes and effects on the larger system. Using lower-hybrid drift waves during magnetic reconnection as an example, we find that the ratio between the wave electric field and the reconnection electric field scales as $\sqrt{m_i/m_e}$, while the phase relationship is also affected. The combination of these effects means that the anomalous drag that contributes to momentum balance in the reconnection region can be underestimated by an order of magnitude. The results are relevant to the coupling of electron-scale waves to ion-scale reconnection regions, and other systems such as collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutron Yield of Thermo Scientific P385 D-T Neutron Generator vs . Current and Voltage

The Thermo Scientific P385 Neutron Generator is a compact neutron source, producing 14 MeV neutrons through the deuterium-tritium (DT) fusion reaction. It is important to measure and understand the dependence of the neutron production rate on the accelerator current and voltage. In this study we evaluated neutron production with an absolutely calibrated liquid scintillator neutron spectrometer (BTI N-Probe), an absolutely calibrated He-3 detector surrounded by HDPE shells (Detec Nested Neutron Spectrometer, NNS), and two uncalibrated ZnS fast neutron scintillators (EJ-410), for both A3082 and A3083 sealed tubes. Here we also modeled the neutron yield using the TRIM code, which calculates the trajectory and the energy loss of deuterons and tritons within the target. Experimental results showed an essentially linear dependence on beam current, as expected. A 3.59 ±0.08 power law dependence on the operating voltage was measured, in effective agreement with the modeled value of 3.5. A series of absolute NNS and N-Probe measurements, matched against MCNP calculations, showed that the A3083 and A3082 tubes provide a maximum neutron yield of 8.2 × 10 8 n/s and 4.7 × 10 8 n/s respectively, with estimated uncertainty of ±10%.We showed, through modeling, that tritium decay is not a significant consideration for tubes, such as these, with lifetimes of less than 10 years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Stellarator Design Exploration Using Symbolic-Regression Neutronics Surrogates

Systems codes require fast, simplified models to rapidly evaluate fusion power plant concepts, but neutronics analyses are often a computational bottleneck. Here, to address this, surrogate models for key neutronics responses have been developed using 3-D neutronics-ready models built with the open-source code ParaStell from a database of stellarator equilibria. Neutronics responses such as tritium breeding ratio (TBR), nuclear heating, and neutron-induced radiation damage displacements per atom (dpa) were simulated using OpenMC. Through sensitivity analysis and symbolic regression (SR), simple power-law formulas were derived connecting these neutronics responses to global stellarator parameters, including fusion power, plasma surface area, and plasma elongation. Validation shows these formulas can predict the simulation results with low error, enabling quick and accurate assessment of neutronics requirements in stellarator design exploration activities with systems codes.

Modeling↗
DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

xGFabric: Coupling Sensor Networks and HPC Facilities with Private 5G Wireless Networks for Real-Time Digital Agriculture

Advanced scientific applications require coupling distributed sensor networks with centralized high-performance computing facilities. Citrus Under Protective Screening (CUPS) exemplifies this need in digital agriculture, where citrus research facilities are instrumented with numerous sensors monitoring environmental conditions and detecting protective screening damage. CUPS demands access to computational fluid dynamics codes for modeling environmental conditions and guiding real-time interventions like water application or robotic repairs. These computing domains have contrasting properties: sensor networks provide low-performance, limited-capacity, unreliable data access, while high-performance facilities offer enormous computing power through high-latency batch processing. Private 5G networks present novel capabilities addressing this challenge by providing low latency, high throughput, and reliability necessary for near-real-time coupling of edge sensor networks with HPC simulations. This work presents xGFabric, an end-to-end system coupling sensor networks with HPC facilities through Private 5G networks. The prototype connects remote sensors via 5G network slicing to HPC systems, enabling real-time digital agriculture simulation.

Digital Agriculture↗
DOE OSTI2025

A Grassroots Network and Community Roadmap for Interconnected Autonomous Science Laboratories for Accelerated Discovery

Scientific discovery is being revolutionized by AI and autonomous systems, yet current autonomous laboratories remain isolated islands unable to collaborate across institutions. We present the Autonomous Interconnected Science Lab Ecosystem (AISLE), a grassroots network transforming fragmented capabilities into a unified system that shorten the path from ideation to innovation to impact and accelerates discovery from decades to months. AISLE addresses five critical dimensions: (1) cross-institutional equipment orchestration, (2) intelligent data management with FAIR compliance, (3) AI-agent driven orchestration grounded in scientific principles, (4) interoperable agent communication interfaces, and (5) AI/ML-integrated scientific education. By connecting autonomous agents across institutional boundaries, autonomous science can unlock research spaces inaccessible to traditional approaches while democratizing cutting-edge technologies. This paradigm shift toward collaborative autonomous science promises breakthroughs in sustainable energy, materials development, and public health.

Ferreira da Silva, Rafael [Oak Ridge National Labo↗
DOE OSTI2025

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗
DOE OSTI2025

Four no-go theorems on the existence of spin and orbital angular momentum of massless bosons

The past decades have seen substantial interest in the so-called orbital angular momentum (OAM) of light, driven largely by its diverse range of applications. However, there are fundamental theoretical issues with decomposing the angular momentum of massless particles, such as photons, into spin (SAM) and orbital angular momentum parts. While the angular momentum of massive particles has a natural splitting into the Wigner SAM and OAM, there are numerous proposed splittings for photons and no consensus about which is correct. Moreover, it has been shown that most of the proposed SAM and OAM operators do not satisfy the defining commutation relations of angular momentum operators and are thus not legitimate splittings. Here, we prove that it is generally impossible to split the total angular momentum operator of massless bosons, such as photons and gravitons, into spin and orbital parts. We prove two further generalizations of this result, showing that there are no SAM-OAM splittings even if (1) the SAM operator generates non-internal symmetries or (2) if one allows the SAM and OAM operators to generate non-SO(3) symmetries.

Chern numbers↗
DOE OSTI2025

Results and lessons learned from accelerating radio frequency modeling using machine learning [slides]

The “advanced tokamak” reactor concept is a leading candidate for a steady state fusion pilot plant. An advanced tokamak (AT) sustains a majority of the required plasma current with effects resulting from maintenance of the peaked pressure at the device center. This current is augmented by auxiliary current drive sources. These auxiliary actuators may consist of neutral particle beams and/or radio frequency (RF) systems such as lower hybrid current drive (LHCD) and high harmonic fast wave (HHFW) current drive using radio and microwaves from antennas. The primary focus of this work is to develop models of RF current profile control suitable for use in integrated modeling frameworks and for real-time control in experiments. Direct physics models of RF current drive can be computationally intensive. In order to achieve predictive times appropriate for the thousands of calls needed in real-time control of experiments and for use in integrated models, we will apply modern machine learning (ML) techniques to accelerate these models and interpolate their results. To generate the fast and accurate models for use in control level algorithms and integrated modeling we need to replace present models with high dimensional interpolation of their results. We will perform additional simulations across a broader parameter range for EAST and other tokamaks in different physics regimes (Alcator C-Mod, DIII-D, WEST, CFETR, ARC, ITER) and combine them into a larger database for training and testing of the ML models. Further testing of the control level models with experimental current profile data from EAST and C-Mod tokamaks will provide additional confirmation of the control level model before integration in a tokamak control system or integrated modeling suite. ML will be used to optimize the selection of training data consisting of RF current driven at different values of density profile, temperature profile, plasma current, and wavenumber. ML will also be used to facilitate classification of current drive from these input data. The output of this effort will be a validated classifier capable of determining the current drive profiles for HHFW CD and LHCD on a mille-second timescale. This will provide a breakthrough capability enabling real-time control of RF driven current profiles in experiments including ITER ICRF and use integrated modeling frameworks requiring thousands of current profile calculations in discharge simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

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↗
DOE OSTI2026

Partially Ionized Plasma Physics and Technological Applications

Partially ionized plasma physics has attracted increased attention recently due to numerous technological applications made possible by the increased sophistication of computer modelling, the depth of the theoretical analysis, and the technological applications to a vast field of manufacturing for computer components. Partially ionized plasma is characterized by a significant presence of neutral particles in contrast to the fully ionized plasma. The theoretical analysis is based upon solutions of the kinetic Boltzmann equation, yielding the non-Maxwellian electron energy distribution function (EEDF), thereby emphasizing the difference with a fully ionized plasma. The impact of the effect on discharges in inert and molecular gases is described in detail, yielding the complex nonlinear phenomena resulting in plasma selforganization. A few examples of such phenomena are given, including the non-monotonic EEDFs in the discharge afterglow in a mixture of argon with the molecular gas NF3; the explosive generation of cold electron populations in capacitive discharges, hysteresis of EEDF in inductively coupled plasmas. Recently, highly advanced computer codes were developed in order to address the outstanding challenges in plasma technology. These developments are briefly described in general terms.

non-Maxwellian electron energy distribution functi↗
DOE OSTI2026

Recent progress in atomic-scale controlled plasma processing

Atomic-scale control in plasma processing is becoming increasingly critical for fabricating of advanced semiconductor devices, particularly as the industry shifts toward three-dimensional (3D) architectures and high-aspect-ratio (HAR) structures. This review presents a comprehensive overview of recent developments in atomic-scale controlled plasma processes, organized along two key directions: the hierarchical structure of plasma–surface interactions and the generational evolution of atomic layer processing (ALP) technologies. We examined the gas phase, where molecular design enables selective generation of ions and radicals; the boundary layer, where transport phenomena govern species delivery into nanoscale features, and the surface, where temperature-dependent reactions and cyclic processing determine etching selectivity and precision. Building on this foundation, we outline five generations of ALP—from thermal atomic layer deposition to transport-aware, temporally and structurally decoupled processes—highlighting the increasing sophistication of process control. The review further explores the transition from empirical recipe development to science-based, data-driven methodologies. By integrating quantum-chemical modeling, advanced diagnostics, and machine learning, we demonstrated how predictive models can link plasma species composition to process outcomes, enabling autonomous and adaptive control strategies. Finally, this review discusses the broader societal implications of plasma process innovation through the E4 quartet: energy and resource efficiency, environmental sustainability, evolutionary advancement, and educational promotion. These principles guide the development of sustainable and intelligent atomic-scale manufacturing technologies that are not only technically advanced but also socially responsible.

Ishikawa, Kenji [Nagoya Univ. (Japan)] (ORCID:0000↗
DOE OSTI2025

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

PaleoSTeHM v1.0: a modern, scalable spatiotemporal hierarchical modeling framework for paleo-environmental data

Abstract. Geological records of past environmental change provide crucial insights into long-term climate variability, trends, non-stationarity, and nonlinear feedback mechanisms. However, reconstructing spatiotemporal fields from these records is statistically challenging due to their sparse, indirect, and noisy nature. Here, we present PaleoSTeHM, a scalable and modern framework for spatiotemporal hierarchical modeling of paleo-environmental data. This framework enables the implementation of flexible statistical models that rigorously quantify spatial and temporal variability from geological data while clearly distinguishing measurement and inferential uncertainty from process variability. We illustrate its application by reconstructing temporal and spatiotemporal paleo-sea-level changes across multiple locations. Using various modeling and analysis choices, PaleoSTeHM demonstrates the impact of different methods on inference results and computational efficiency. Our results highlight the critical role of model selection in addressing specific paleo-environmental questions, showcasing the PaleoSTeHM framework's potential to enhance the robustness and transparency of paleo-environmental reconstructions.

58 GEOSCIENCES↗
DOE OSTI2026

Production of High Yield, Ultra‐Shallow Nitrogen‐Vacancy Centre Ensembles Through High Purity Forming gas Annealing

Dense, ultra‐shallow (<10 nm from the surface) ensembles of nitrogen‐vacancy (NV) centers in diamond are sought‐after for realising practical quantum sensing at the nanoscale. However, these ensembles typically suffer from reduced NV creation efficiency and spin coherence compared to their deeper and bulk counterparts. Consequently, most demonstrations in the literature to date have been restricted to using single NV centers with atypical coherence and charge stability properties not achievable across entire ensembles. Here, in this work, we demonstrate that controlled annealing in a high purity forming gas atmosphere improves near‐surface NV formation efficiency to over 20%. This environment preserves a high quality hydrogen termination throughout the anneal, avoiding material etching over many hours at temperatures up to 1100°C and possibly helping to suppress vacancy out‐diffusion. We show that high temperature annealing is especially vital for improving spin properties in high density, ultra‐shallow implanted samples due to the localised damage profiles produced. By allowing the formation of high yield NV ensembles with improved spin properties, the optimised forming gas anneal can facilitate the reliable production of more sensitive diamond sensors than have previously been accessible.

Healey, Alexander J. [Royal Melbourne Institute of↗
DOE OSTI2026

Vacuum Neutral Transport Model in UEDGE for Tokamak Far Scrape‐Off Layer

A model for neutral transport in the far scrape-off layer (SOL) vacuum region (vacuum neutral model) has been developed and implemented in UEDGE. Free-streaming neutral trajectories between the outermost UEDGE boundary and the vessel wall are preprocessed using DEGAS2 to construct a tele-transport matrix that captures non-local neutral relocation through the vacuum. Here, this matrix is then used in UEDGE as a non-local boundary condition for the neutral equations, preserving the robustness and convergence of the implicit solver without introducing statistical noise. Simulations of a DIII-D lower single-null configuration show that the vacuum neutral model relocates neutrals from the divertor to the upstream region, increasing the outer midplane separatrix density required for detachment onset by about 30%. However, the characteristic target temperature at detachment (T e, osp ~ 3 - 4eV) and the radiation front behavior remain unchanged.

DEGAS2↗