Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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

Observing differential spin currents by resonant inelastic X-ray scattering

Controlling spin currents, that is, the flow of spin angular momentum, in small magnetic devices, is the principal objective of spin electronics, a main contender for future energy-efficient information technologies. A pure spin current has never been measured directly because the associated electric stray fields and/or shifts in the non-equilibrium spin-dependent distribution functions are too small for conventional experimental detection methods optimized for charge transport. Here we report that resonant inelastic X-ray scattering (RIXS) can bridge this gap by measuring the spin current carried by magnons—the quanta of the spin wave excitations of the magnetic order—in the presence of temperature gradients across a magnetic insulator. This is possible due to the sensitivity of the momentum- and energy-resolved RIXS intensity to minute changes in the magnon distribution under non-equilibrium conditions. Furthermore, we use the Boltzmann equation in the relaxation time approximation to extract transport parameters, such as the magnon lifetime at finite momentum, essential for the realization of magnon spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multifrequency Analysis of Favored Models for the Messier 87* Accretion Flow

Abstract The polarized images of the supermassive black hole Messier 87* (M87*) produced by the Event Horizon Telescope (EHT) provide a direct view of the near-horizon emission from a black hole accretion and jet system. The EHT theoretical analysis of the polarized M87* images compared thousands of snapshots from numerical models with a variety of spins, magnetization states, viewing inclinations, and electron energy distributions, and found a small subset consistent with the observed image. In this article, we examine two models favored by EHT analyses: a magnetically arrested disk with moderate retrograde spin and a magnetically arrested disk with high prograde spin. Both have electron distribution functions that lead to strong depolarization by cold electrons. We ray trace five snapshots from each model at 22, 43, 86, 230, 345, and 690 GHz to forecast future very long baseline interferometry (VLBI) observations and examine limitations in numerical models. We find that even at low frequencies where optical and Faraday rotation depths are large, approximately rotationally symmetric polarization persists, suggesting that shallow depths dominate the polarization signal. However, morphology and spectra suggest that the assumed thermal electron distribution is not adequate to describe emission from the jet. We find 86 GHz images show a ringlike shape determined by a combination of plasma and spacetime imprints, smaller in diameter than recent results from the Global mm-VLBI Array. We find that the photon ring becomes more apparent with increasing frequency, and is more apparent in the retrograde model, leading to large differences between models in asymmetry and polarization structure.

Astronomy & Astrophysics↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. I. Stochastic Fragmentation with Optical-depth-dependent Rate and Universal Fragment Mass

Fragmentation in a gravitationally unstable accretion disk can be an important pathway for forming stellar/planetary companions. To characterize quantitatively the condition for and outcome of fragmentation under realistic thermodynamics, we perform global 3D simulations of gravitationally unstable disks at various cooling rates and cooling types, including the first global simulations of gravitational instability that employ full radiation transport. We find that fragmentation is a stochastic process, with the fragment generation rate per disk area p frag showing an exponential dependence on the parameter β ≡ Ω K t cool , where Ω K is the Keplerian rotation frequency, and t cool is the average cooling timescale. Compared to a prescribed constant β, radiative cooling in the optically thin/thick regime makes p frag decrease slower/faster in β; the critical β corresponding to ∼1 fragment per orbit is ≈3, 5, and 2 for constant β, optically thin, and optically thick cooling, respectively. The distribution function of the initial fragment mass is remarkably insensitive to disk thermodynamics. Regardless of cooling rate and optical depth, the typical initial fragment mass is m frag ≈ 40M tot h 3 , with M tot being the total (star+disk) mass and h = H/R being the disk aspect ratio. Applying this result to typical Class 0/I protostellar disks, we find m frag ∼ 20M J , suggesting that fragmentation more likely forms brown dwarfs. Given the finite width of the m frag distribution, forming massive planets is also possible.

gravitational instability↗
DOE OSTI2026

Determination of nuclear PDFs using Markov chain Monte Carlo methods

Global QCD analyses of nuclear parton distribution functions (nPDFs) have traditionally relied on the Hessian method for uncertainty estimation. However, the inherent Gaussian approximation and reliance on local curvature often prove insufficient for nPDF fits, which are frequently characterized by limited data constraints and non-Gaussian likelihoods. In this paper, we present the first nPDF determination based on Markov Chain Monte Carlo (MCMC) techniques, implemented within the nCTEQ framework using an adaptive Metropolis-Hastings algorithm. The MCMC approach enables a direct mapping of the posterior distribution and reveals a highly nontrivial parameter-space structure, including multiple modes and pronounced non-Gaussian behavior, particularly for the valence PDFs. We perform the first single-nucleus global analysis of lead PDFs using exclusively lead data and compare it to a multi-nuclei fit employing a standard analytic A dependence. The inclusion of lighter nuclei reduces quark uncertainties and modifies the shape of the lead PDFs, while leaving the gluon distribution largely unaffected. A complementary Hessian analysis exposes systematic limitations of the Gaussian approximation. Our results demonstrate that MCMC methods provide a more reliable framework for uncertainty quantification in nPDF determinations.

Derakhshanian, N. [Institute of Nuclear Physics Po↗
DOE OSTI2023

Evidence of non-Maxwellian ion velocity distributions in spherical shock-driven implosions

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature (< 10 keV), hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. Furthermore, these measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Machine learning-based interatomic potential development and phase transition analysis of ferroelectric hafnium dioxide

The ferroelectric phase (𝑃⁢𝑐⁢𝑎⁢2 1 , which is in orthorhombic symmetry) of hafnium dioxide (HfO 2 ) has gained much attention due to its potential applications in nanoelectronics and advanced memory devices. However, its complex phase behavior under external stimuli, such as pressure and temperature, remains a subject of intense investigation. This study focuses on developing a machine learning-based interatomic potential (MLIP) that is trained with data from density-functional theory (DFT) calculations to simulate phase transitions and mechanical properties of HfO 2 . The developed MLIP predicts lattice parameters, equations of state, bulk and shear moduli, and elastic constants that closely align with DFT predictions for several phases and at various pressures. Once validated, the MLIP is used to investigate the phase transitions of ferroelectric HfO 2 (𝑃⁢𝑐⁢𝑎⁢2 1 ) under both isobaric and constant stress conditions at elevated temperatures ranging from 200 to 2500 K. We used several complementary methods, including local symmetry identification, radial distribution function, and x-ray diffraction characterization, to identify interesting phase transitions among several competitive hafnia phases predicted from our simulations. The suggested methods uniformly reveal that under pure deviatoric condition, the system favors a transition from the orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase to a tetragonal (𝑃⁢4 2 /𝑛⁢𝑚⁢𝑐) phase, whereas a zero stress condition drives the system from the 𝑃⁢𝑐⁢𝑎⁢2 1 phase to another orthorhombic (𝑃⁢𝑏⁢𝑐⁢𝑛) phase. These findings provide crucial insights into stress and temperature-induced phase behavior of hafnia, guiding future experimental and theoretical studies for optimizing hafnia-based ferroelectric devices.

Ferroelectric HfO2↗
DOE OSTI2023

Phosphate‐functionalized Zirconium Metal–Organic Frameworks for Enhancing Lithium–Sulfur Battery Cycling

Abstract Lithium–sulfur batteries are promising candidates for next‐generation energy storage devices due to their outstanding theoretical energy density. However, they suffer from low sulfur utilization and poor cyclability, greatly limiting their practical implementation. Herein, we adopted a phosphate‐functionalized zirconium metal–organic framework (Zr‐MOF) as a sulfur host. With their porous structure, remarkable electrochemical stability, and synthetic versatility, Zr‐MOFs present great potential in preventing soluble polysulfides from leaching. Phosphate groups were introduced to the framework post‐synthetically since they have shown a strong affinity towards lithium polysulfides and an ability to facilitate Li ion transport. The successful incorporation of phosphate in MOF‐808 was demonstrated by a series of techniques including infrared spectroscopy, solid‐state nuclear magnetic resonance spectroscopy, and X‐ray pair distribution function analysis. When employed in batteries, phosphate‐functionalized Zr‐MOF (MOF‐808‐PO4) exhibits significantly enhanced sulfur utilization and ion diffusion compared to the parent framework, leading to higher capacity and rate capability. The improved capacity retention and inhibited self‐discharge rate also demonstrate effective polysulfide encapsulation utilizing MOF‐808‐PO4. Furthermore, we explored their potential towards high‐density batteries by examining the cycling performance at various sulfur loadings. Our approach to correlate structure with function using hybrid inorganic–organic materials offers new chemical design strategies for advancing battery materials.

Liu, Bingqian↗
DOE OSTI2022

Towards understanding the magnetic properties of the breathing pyrochlore compound Ba3Yb2Zn5O11through single-crystal studies

Abstract Ba 3 Yb 2 Zn 5 O 11 is exceptional among breathing pyrochlore compounds for being in the nearly-decoupled limit where inter-tetrahedron interactions are weak, hosting isolated clusters or molecular magnet-like tetrahedra of magnetic ytterbium (Yb 3+ ) ions. In this work, we present the study carried out on single-crystal samples of the breathing pyrochlore Ba 3 Yb 2 Zn 5 O 11 , using a variety of magnetometry and neutron scattering techniques along with theoretical modeling. We employ inelastic neutron scattering to investigate the magnetic dynamics as a function of applied field (with respect to both magnitude and direction) down to a temperature of 70 mK, where inelastic scattering reveals dispersionless bands of excitations as found in earlier powder sample studies, in good agreement with a single-tetrahedron model. However, diffuse neutron scattering at zero field and dc-susceptibility at finite field exhibit features suggesting the presence of excitations at low-energy that are not captured by the single tetrahedron model. Analysis of the local structure down to 2 K via pair distribution function analysis finds no evidence of structural disorder. We conclude that effects beyond the single tetrahedron model are important in describing the low-energy, low-temperature physics of Ba 3 Yb 2 Zn 5 O 11 , but their nature remains undetermined.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic structural views in solar energy materials by femtosecond electron diffraction

Solar materials convert light into other forms of energy through excited state processes occurring on ultrafast time and atomic length scales. Understanding and controlling such nonequilibrium processes is essential for applications ranging from photovoltaics to photocatalysis. These processes are commonly studied by transient optical techniques, which resolve charge carrier dynamics. However, optical spectroscopy does not directly probe the atomic-scale motions that play a central role in the ultimate functionality of these materials. Filling this missing gap, ultrafast electron and x-ray scattering techniques enable monitoring of structural dynamics in materials with femtosecond time resolution. In this work, we focus on the use of ultrafast electron diffraction (UED) techniques to probe photoinduced energy conversion dynamics in solar energy materials. First, we discuss the use of UED as an ultrafast lattice thermometer allowing a direct probe of nonradiative relaxations such as hot carrier cooling and Auger recombination. Second, we present a time-resolved atomic pair distribution function analysis enabled by UED, which uncovers lattice deformations arising from excitation localization in materials. Finally, we provide an outlook looking toward new approaches for resolving in situ energy conversion dynamics and uncovering structure–property relationships in solar energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

From Structure to Function: Zn/Mn-Modified Maghemite as an Advanced Nanoplatform for Magnetic Hyperthermia and Radionuclide Therapy

The development of nanoplatforms capable of efficient heat generation and stable radionuclide delivery is essential for effective bimodal cancer therapy. Here, in this study, binary (Fe–M) and ternary (Fe–M–M′) metal oxide nanoparticles were synthesized via a polyol method optimized to produce flower-like γ-Fe 2 O 3 (maghemite) structures, with M and M′ representing Zn and/or Mn. Comprehensive structural and magnetic characterization was conducted to explain the relationship between composition, defect structure, and hyperthermic performance. The analyses revealed that cation substitution induced an Fe-site vacancy, primarily at octahedral positions, leading to local structural distortions, as confirmed by powder X-ray diffraction and pair distribution function analysis. The optimized composition, with Zn/Mn/Fe = 0.040:0.182:1, exhibited the highest concentration of vacancies and structural disorder. These vacancies altered the bonding environment, enhancing magnetic interactions at tetrahedral sites while weakening those at the octahedral positions. The resulting multicore nanoflowers (20–63 nm; core size 13–18 nm) displayed strong heating performance, with intrinsic loss power ranging from 0.34 to 5.77 nHm 2 kg –1 . The optimized sample achieved a temperature increase of 30 °C within 2 min and a specific absorption rate of 369 W g –1 . This composition was further coated with citrate (CA) and successfully radiolabeled with 177 Lu, achieving a radiolabeling yield of 92.7% and excellent stability, thus forming a robust nanoplatform for combined magnetic hyperthermia and radionuclide therapy. Biological evaluation of the optimized S5 composition revealed selective cytotoxicity toward HeLa and LS174 cells, while toxicity was significantly lower to A549, A375, and normal MRC-5 cells. Citrate coating of S5 nanoparticles (S5@CA) drastically reduced their cytotoxicity across all tested cell lines (IC 50 > 200 μg mL –1 ), confirming their enhanced biocompatibility for therapeutic applications. In HeLa cells subjected to magnetic hyperthermia, the viability decreased to approximately 84% after 30 min and 61% after 60 min of treatment, demonstrating the sustained hyperthermic effect at a controlled working temperature of 48 °C. These results underscore the effectiveness of cation substitution and vacancy engineering in tailoring the functional properties of maghemite-based nanomaterials for advanced multimodal cancer therapies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissolving salt is not equivalent to applying a pressure on water

Abstract Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

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

Generalized parton distributions through universal moment parameterization: non-zero skewness case

We present the first global analysis of generalized parton distributions (GPDs) combing lattice quantum chromodynamics (QCD) calculations and experiment measurements including global parton distribution functions (PDFs), form factors (FFs) and deeply virtual Compton scattering (DVCS) measurements. Following the previous work where we parameterize GPDs in terms of their moments, we extend the framework to allow for the global analysis at non-zero skewness. Together with the constraints at zero skewness, we fit GPDs to global DVCS measurements from both the recent JLab and the earlier Hadron-Electron Ring Accelerator (HERA) experiments with two active quark flavors and leading order QCD evolution. With certain choices of empirical constraints, both sea and valence quark distributions are extracted with the combined inputs, and we present the quark distributions in the proton correspondingly. We also discuss how to extend the framework to accommodate more off-forward constraints beyond the small ξ expansion, especially the lattice calculated GPDs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Local structure study on magnetostrictive material Tb 1–x Dy x Fe 2

Tb 1–x Dy x Fe 2 system has attracted more research interest due to the large magnetostrictive effect. The crystal structures and physical properties have been well studied, but research studies on their local structures are still rare. As such, in this work, the local structure of Tb 1–x Dy x Fe 2 samples was studied using the pair distribution function and x-ray absorption spectroscopy techniques. The results demonstrate that the system owns the same local crystal symmetry with its average structure in the ferromagnetic phase, and the crystal lattice of the system is more ordered with increasing Dy content, indicating that the Dy-rich tetragonal phase is more stable than the Tb-rich rhombohedral phase. The different roles of metallic bonds in affecting the crystal lattice are presented. The weak Fe 1 –Fe 2 bonds influenced by the local environment such as local stress from randomly distributed nanodomains could originate the anomalies in the lattice, resulting in the more ordered and stable Dy-rich phase than the Tb-rich phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influences of δB contribution and parallel inertial term of energetic particles on MHD-kinetic hybrid simulations: a case study of the 1/1 internal kink mode

The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model (Park et al 1992 Phys. Fluids B 4 2033–7) has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling scheme is adopted in this model. However, two noteworthy issues arise in the model application: firstly, the coupled term introduced in the pressure-coupling scheme, (∇•P h ) ⟂ , is often simplified by ∇•P h , which is equivalent to neglecting the parallel inertial term of EPs; secondly, besides the $δf$ contribution caused by changing in the EP distribution function, the magnetic field perturbation (the $δB$ contribution) generated during development of the instabilities should also be considered, but it is often ignored in existing hybrid simulations. In this paper, we derive the analytical formulations under these two coupling schemes and then numerically study the representative case of the linear stability of the $m/n$ = $1/1$ internal kink mode (IKM) (Fu et al 2006 Phys. Plasmas 13 052517) by using the CLT-K code. Further, it is found that the approximated models can still yield reasonable results when EPs are isotopically distributed. But it fails completely in cases with anisotropic EP distributions. In addition, we further investigate the influence of EP's orbit width on the stability of IKM and verify the equivalence between pressure-coupling scheme and the current-coupling scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Detailed study of quark-hadron duality in spin structure functions of the proton and neutron

Background: The response of hadrons, the bound states of the strong force (QCD), to external probes can be described in two different, complementary frameworks: as direct interactions with their fundamental constituents, quarks and gluons, or alternatively as elastic or inelastic coherent scattering that leaves the hadrons in their ground state or in one of their excited (resonance) states. The former picture emerges most clearly in hard processes with high momentum transfer, where the hadron response can be described by the perturbative expansion of QCD, while at lower energy and momentum transfers, the resonant excitations of the hadrons dominate the cross section. The overlap region between these two pictures, where both yield similar predictions, is referred to as quark-hadron duality and has been extensively studied in reactions involving unpolarized hadrons. Some limited information on this phenomenon also exists for polarized protons, deuterons, and 3 He nuclei, but not yet for neutrons. Purpose: In this paper, we present comprehensive and detailed results on the correspondence between the extrapolated deep inelastic structure function g 1 of both the proton and the neutron with the same quantity measured in the nucleon resonance region. Thanks to the fine binning and high precision of our data, and using a well-controlled perturbative QCD (pQCD) fit for the partonic prediction, we can make quantitative statements about the kinematic range of applicability of both local duality and global duality. Method: We use the most updated QCD global analysis results at high x from the Jefferson Lab Angular Momentum Collaboration to extrapolate the spin structure function g 1 into the nucleon resonance region and then integrate over various intervals in the scaling variable x. We compare the results with the large data set collected in the quark-hadron transition region by the CLAS Collaboration, including, for the first time, deconvoluted neutron data, integrated over the same intervals. We present this comparison as a function of the momentum transfer Q 2 . Results: We find that, depending on the integration interval and the minimum momentum transfer chosen, a clear transition to quark-hadron duality can be observed in both nucleon species. Furthermore, we show, for the first time, the approach to scaling behavior for g 1 measured in the resonance region at sufficiently high momentum transfer. Conclusions: Here, our results can be used to quantify the deviations from the applicability of pQCD for data taken at moderate energies and can help with extraction of quark distribution functions from such data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of W±-boson differential cross-sections in proton–proton collisions with low pile-up data at s=5.02 TeV and 13TeV with the ATLAS detector

High precision single-differential W±$$W^\pm $$-boson production cross-sections as a function of electron or muon transverse momentum pT$$p_\textrm{T}$$ or their pseudorapity η$$\eta $$, as well as double-differential cross-sections as functions of these variables, are measured in proton–proton collisions at centre-of-mass energies s=5.02$$\sqrt{s}=5.02$$ TeV and 13 TeV. The W-boson charge asymmetry as a function of lepton η$$\eta $$ is also measured. The data, collected in dedicated runs at reduced instantaneous luminosity with the ATLAS detector at the Large Hadron Collider, correspond to integrated luminosities of 255 pb-1$$^{-1}$$ at 5.02 TeV and 338 pb-1$$^{-1}$$ at 13 TeV. The measurements are in agreement with Standard-Model predictions calculated at next-to-next-to-leading-order in the strong coupling constant αs$$\alpha _s$$ including transverse-momentum resummation at next-to-next-to-leading logarithmic accuracy using several parton distribution functions. The impact of the measured differential cross-sections as a function of lepton η$$\eta $$ on the determination of these functions is studied using a profiling technique.

Aad, G↗