SEARCH · Search NASA

Results for “Phase-space distribution”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Minimizing phase-space energies

A primary technical challenge for harnessing fusion energy is to control and extract energy from a nonthermal distribution of charged particles. The fact that phase space evolves by symplectomorphisms fundamentally limits how a distribution may be manipulated. While the constraint of phase-space volume preservation is well understood, other constraints remain to be fully appreciated. To better understand these constraints, we study the problem of extracting energy from a distribution of particles using area-preserving and symplectic linear maps. When a quadratic potential is imposed, we find that the maximal extractable energy can be computed as trace minimization problems. We solve these problems and show that the extractable energy under linear symplectomorphisms may be much smaller than the extractable energy under special linear maps. As a result, the method introduced in the present study enables an energy-based proof of the linear Gromov nonsqueezing theorem.

First-principles calculations in plasma physics↗

DOE OSTI2026

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗

NASA NTRS1982

Plasma electron signature of magnetic connection to the earth's bow shock - ISEE 3

New observations of bowshock modified electron velocity distribution for upstreams are investigated using the Los Alamos ISEE 3 electron analyzer. Examples of two-dimensional electron distributions observed when ISEE 3 magnetically connected to the bowshock are presented, and the velocity moments are determined. The origin of the enhanced backstreaming electron fluxes and their likely effects on the upstream microturbulence spectrum are interpreted. A nearly isotropic component of back-streaming shock-heated electrons is found with energies at least as high as 1 keV. Phase-space density of the low energy parts of distributions measures well within the boundaries of accessibility and is observed to be depressed. Such distributions are likely to result from the direct sampling of electrons originating within the forward edge of the earth's bowshock, and may be unstable to whistler waves.

Feldman, W. C.↗

NASA NTRS1982

A numerical method for solving the Vlasov equation

A numerical procedure is derived for the solution of the Vlasov-Poisson system of equations in two phase-space variables. Derivatives with respect to the phase-space variables are approximated by a weighted sum of the values of the distribution function at property chosen neighboring points. The resulting set of ordinary differential equations is then solved by using an appropriate time intergration scheme. The accuracy of the proposed method is tested with some simple model problems. The results for the free streaming case, linear Landau damping, and nonlinear Landau damping are investigated and compared with those of the splitting scheme. The proposed method is found to be very accurate and efficient.

Satofuka, N.↗

NASA NTRS1976

Lossy radial diffusion of relativistic Jovian electrons

The radial diffusion equation with synchrotron losses is solved by the Laplace-transform method for near equatorially mirroring relativistic electrons. The evolution of a power-law distribution function is found, and the characteristics of synchrotron burnoff are stated in terms of explicit parameters for an arbitrary diffusion coefficient of a specific form. The peaking of the 10.4-cm volume emissivity from Jupiter at an L shell of about 1.8 provides an estimate of the diffusion coefficient in the radiation belts; one value is suggested as the appropriate modification, for an equatorial field strength of 4.2 G, of the Birmingham et al. (1974) result. Nonsynchrotron losses are included phenomenologically; from the phase-space densities reported by McIlwain and Fillius (1975), the particle lifetime is estimated. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.

Barbosa, D. D.↗

DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Characterization and Optimization of the Fitting of Quantum Correlation Functions

This case study presents a characterization and optimization of an application code for extracting parton distribution functions from high energy electron-proton scattering data. Profiling this application code reveals that the phase-space density computation accounts for 93% of the overall execution time for a single iteration on a single core. When executing multiple iterations in parallel on a multicore system, the application spends 78% of its overall execution time idling due to load imbalance. We address these issues by first transforming the application code from Python to C++ and then tackling the application load imbalance via a hybrid scheduling strategy that combines dynamic and static scheduling. These techniques result in a 62% reduction in CPU idle time and a 2.46x speedup in overall execution time per node. In addition, the typically enabled power-management mechanisms in supercomputers (e.g., AMD Turbo Core, Intel Turbo Boost, and RAPL) can significantly impact intra-node scalability when more than 50% of the CPU cores are used. This finding underscores the importance of understanding system interactions with power management, as they can adversely impact application performance, and highlights the necessity of intra-node scaling tests to identify performance degradation that inter-node scaling tests might otherwise overlook.

Chuang, Pi-Yueh [Virginia Tech,Dept. of Computer S↗

DOE OSTI2026

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

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 OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗

DOE OSTI2025

First fluctuation measurements using an Imaging Neutral Particle Analyzer on DIII-D

A recent upgrade to the Imaging Neutral Particle Analyzer (INPA) on DIII-D has allowed for the first fluctuation measurements using an INPA to be taken during a neoclassical tearing mode (NTM). The INPA signal tracked the mode over 150 ms as the mode frequency dropped to zero, capturing both the NTM with poloidal and toroidal mode numbers $m/n=2/1$ and a $3/2$ mode at double the frequency. Analysis shows that the signals originate from charge exchange events near the edge of the plasma, and relative fluctuation amplitudes are greater than 25% for the duration of the NTM. Filtered signals show frequency beating patterns that are phase-space dependent. Simulated signal is dominated by prompt transport from the neutral beams to the INPA sightline, while the contribution from the slowing down distribution is significantly lower. Simulated measurements in the range of pitches that the diagnostic is sensitive to (0.5≤|v∥/v|≤ 0.75) show the signal is dominated by trapped orbits that pass through magnetic islands near the edge of the plasma. Calculations of expected fluctuation levels show that only direct interaction with the NTM can provide the strong relative fluctuation levels seen experimentally. The prompt nature of the orbits and thin radial layer found to contribute to synthetic signals suggest INPA passive data may be used to measure the perturbations of confined orbits on a single pass through a plasma instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗

DOE OSTI2025

Unified and optimal frame choice for generalized parton distributions

Reconstructing the internal three-dimensional quark and gluon structures of hadrons through generalized parton distributions (GPDs) from hard exclusive scattering processes is one of the most challenging tasks in nuclear and particle physics. In this paper, we introduce a new optimized reference frame that, for the first time, enables a unified view of all the reactions sensitive to GPDs and facilitates the interpretation of a variety of phase-space patterns that were previously hardly accessible and interpretable. Similarly to how the heliocentric description advanced our understanding of the solar system and gravitation, our new frame centers around a quasireal state, allows for a consistent separation of physical scales, and reveals a novel quantum interference mechanism. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Interaction between photons and leaf canopies

The physics of neutral particle interaction for photons traveling in media consisting of finite-dimensional scattering centers that cross-shade mutually is investigated. A leaf canopy is a typical example of such media. The leaf canopy is idealized as a binary medium consisting of randomly distributed gaps (voids) and regions with phytoelements (turbid phytomedium). In this approach, the leaf canopy is represented by a combination of all possible open oriented spheres. The mathematical approach for characterizing the structure of the host medium is considered. The extinction coefficient at any phase-space location in a leaf canopy is the product of the extinction coefficient in the turbid phytomedium and the probability of absence gaps at that location. Using a similar approach, an expression for the differential scattering coefficient is derived.

Knyazikhin, Yuri V.↗

DOE OSTI2025

Status of QCD precision predictions for Drell–Yan rapidity distributions

We compute differential distributions for Drell–Yan processes at the LHC and the Tevatron colliders at next-to-next-to-leading order in perturbative QCD, including fiducial cuts on the decay leptons in the final state. The comparison of predictions obtained with four different codes shows excellent agreement, once linear power corrections from the fiducial cuts are included in those codes that rely on phase-space slicing subtraction schemes. For Z-boson production we perform a detailed study of the symmetric cuts on the transverse momenta of the decay leptons. Predictions at fixed order in perturbative QCD for those symmetric cuts, typically imposed in experiments, suffer from an instability. We show how this can be remedied by an all-order resummation of the fiducial transverse momentum spectrum, and we comment on the choice of cuts for future experimental analyses.

Alekhin, S. [University of Hamburg (Germany)]↗

DOE OSTI2026

ELG×LRG Distribution through Dark Matter Halo Dynamics

We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at 0.8 < z < 1.1, including their cross-correlation (ELG×LRG), using the A BACUS S UMMIT suite and a new Halo Occupation Model (H OME ) for galaxy multitracers. This integrates intrahalo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semianalytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, Home reproduces the anisotropic clustering down to s = 200 h −1 kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG, and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below 4 h −1 Mpc. The HOD from the best-fit Home has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of M vir ∼ 6.6 × 10 11 (1.2 × 10 13 ) h −1 M ⊙ ; (ii) the ELG×LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on 2–5 h −1 Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. The high sensitivity of Home precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics and the dynamical nature of dark matter, potentially distinguishing gravity models.

Favole, Ginevra [Universidad de La Laguna (Spain);↗

NASA NTRS1978

On the formation and evolution of clumps of galaxies in an expanding universe

Results are derived for the development of phase-space clumps of mass points in a background spectrum of gravitational-potential fluctuations. The Vlasov equation and the pair correlation equation (in the weak coupling limit) are solved exactly in an Einstein-de Sitter cosmology, and the plasma-clumping theory is used to identify terms that yield important collective effects. Various astrophysical implications are discussed, including the formation of large-scale inhomogeneity and the enhanced generation of correlations in the distribution of galaxies.

Norman, C. A.↗

DOE OSTI2025

Measurement of three-dimensional inclusive muon-neutrino charged-current cross sections on argon with the MicroBooNE detector

We report the measurement of the triple-differential cross section d 3 σ/dE vis d cos(θ μ )dP μ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4 x 10 20 protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8 GeV. The mapping from reconstructed kinematics to truth quantities is validated within uncertainties by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after applying a conditional constraint from the muon momentum distribution in data. The success of this validation provides confidence that the energy transfer in the MicroBooNE detector is well-modeled within simulation uncertainties, enabling a reliable unfolding to a triple-differential cross section defined at the nominal neutrino flux over muon momentum, muon scattering angle, and visible neutrino energy. This validation not only supports accurate cross-section extraction, but also establishes a critical foundation for tuning interaction models used in future neutrino oscillation measurements. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future liquid argon time projection chamber experiments measuring neutrino oscillations. Comparisons against a number of commonly used model predictions are included and their performance in different parts of the available phase-space is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗