Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase-space”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Python Script to Read MCNP6.3 Surface-Source Files

This report provides a Python script to read an MCNP ® surface source file created with the SSW card with SYM = 0 (i.e., the default symmetry treatment). For background: the general format of an MCNP surface-source file is described in; however, that document did not provide coding and/or a tool to interrogate such files. The current format will not be given in this document other than through the record-read statements necessary for the script to function. The reader capability in this report is augmented with the ability to directly write a couple demonstrative outputs: 1. A comma-separated value (CSV) file containing particle phase-space state information and 2. A Matplotlib histogram of the energy distribution of the particles. This report also describes accompanying verification work that shows the script performing as required with MCNP6.2, MCNP6.3, and (expected) MCNP6.4 surface-source files. However, users of the enclosed script must still verify that the script is behaving correctly for their own work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Cluster for Research Excellence (CRE) for Accelerator Physics

The purpose of this Cooperative Research and Development Agreement (CRADA) was to establish a Cluster of Research Excellence (CRE) for Accelerator Physics and broaden the accelerator physics collaboration between Northern Illinois University (NIU) and Fermilab. The CRE was based in both Fermilab’s Accelerator Division and in NIU’s Department of Physics. A research agenda and strategy for the CRE for Accelerator Physics included specific research topics, experiments, and research facility upgrades. The scope of the research program covered two thrusts with subtopics in the Accelerator Science Research Program utilizing the infrastructure of FAST, IOTA, the Fermilab Accelerator Complex, and the Illinois Accelerator Research Center (IARC). Specific areas of proposed exploration included: 1.Tests of novel “nonlinear integrable” systems including beam self-fields; 2.Experimental tests of phase-space diffusion and chaos in intense beams; 3.Investigate technically feasible approaches towards next generation intenseneutrino facilities; 4.Understanding optical control of charged particle beams; 5.Exploitation of bright electron beams from FAST for novel applications; 6. Conduct research into accelerator-based ‘precision science’, such as Muon g-2, EDM, Mu2e, etc.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Study of Particle Loss in Synchrotron Phase Space Injection for ESR Using Weak-Strong Beam-Beam Simulation with Nonlinear Lattice

In this report, we use tracking simulations to investigate synchrotron phase space injection for electron accumulation in the electron storage ring of the Electron Ion Collider. Our simulation model accounts for both beam-beam interactions and lattice nonlinearities. Specifically, we examine how particle loss is influenced by various parameters. Additionally, we conduct a theoretical analysis and derive an analytical formula for the rapid evaluation of particle loss. Our results demon strate the feasibility of synchrotron phase-space injection for the electron storage ring and provide insights to guide parameter selection for the design of the injection line.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

LHC Event Generation in the Exascale Era

MCFM is a dedicated Monte-Carlo simulation program for collider phenomenology at highest energies. Designed during the Tevatron era, it has successfully incorporated the latest developments needed for LHC precision calculations and remained on the forefront of collider phenomenology. The Fortran code includes interfaces to modern PDF and loop reduction libraries but has been unchanged structurally compared to the earlier versions. Parallel computing has been enabled using OpenMP and MPI. MCFM provides numerically highly stable one-loop amplitudes and superior phase-space efficiency, leading to excellent performance in NXLO calculations using jettiness or qT subtraction techniques for IR regularization.

Campbell, John [Fermilab]↗
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

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

Reconstructing the Stripping History of the Sagittarius Stream with Neural Networks

The Sagittarius (Sgr) Stream is produced by the ongoing disruption of the Sgr dwarf spheroidal (dSph) galaxy and is thought to contain multiple wraps that were stripped during different pericentric passages. In this study, we introduce a neural-network–based method trained on N-body simulations to infer the stripping time of Sgr Stream stars directly from their phase-space coordinates. We combine spectroscopic data from SEGUE, APOGEE DR17, and LAMOST DR7 low-resolution spectroscopic (LRS) survey with Gaia EDR3 astrometry and distance estimates from the latest StarHorse catalog to identify high-quality Sgr Stream members. Applying our method to these stars, we measure a clear metallicity gradient with stripping time, well described by a linear relation with slope ∼0.3 dex Gyr −1 . We further predict the stripping times of globular clusters previously suggested to originate from the Sgr dSph. M 54, Terzan 7, Terzan 8, and Arp 2 exhibit stripping times consistent with being currently bound to the Sgr remnant. Pal 12, Whiting 1, and NGC 2419 are inferred to have been stripped 0.9 ± 0.1, 1.1 ± 0.2, and 2.1 ± 0.2 Gyr ago, respectively. For NGC 4147 and NGC 5634, whose membership in the Sgr system remains uncertain, our analysis suggests stripping times of 1.1 ± 0.4 and 1.1 ± 0.1 Gyr, respectively, if they are ultimately confirmed as genuine Sgr members. These results demonstrate that data-driven models of dynamical stripping histories offer a promising approach for reconstructing the formation and chemical evolution of the Sgr Stream.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2023

Measuring the Beam-Spin Asymmetry of Hard Exclusive ¿0 Production off the Neutron with CLAS12

Over the course of the last couple of decades, Generalised Parton Distributions (GPDs) have become a compelling area of focus in both theoretical and experimental research in the ?eld of hadronic structure. GPDs can be indirectly accessed via observables in hard exclusive reactions such as Deeply Virtual Compton Scattering (DVCS) and Deeply Virtual Meson Production (DVMP). In these reactions, an incident lepton scatters o? an individual quark within the target nucleon via the exchange of a virtual photon, and a real photon (DVCS) or a meson (DVMP) are emitted as a result. This thesis presents an analysis of Deeply Virtual ?0 Meson Production o? the neutron in the deuteron (nDV?0P). The data were taken with the CLAS12 detector, which is housed in the experimental Hall B at the Je?erson National Laboratory (JLab) in Virginia, USA. A longitudinally polarised electron beam delivered by JLab?s Continuous Electron Beam accelerator Facility (CEBAF), at an energy of 10.2 and 10.6 GeV, was incident on a liquid deuterium target housed within the CLAS12 detector assembly. The nDV?0P channel is sensitive to the lesser-studied chiral-odd transversity GPDs. GPDs are accessed indirectly via observables related to the cross-sections of hard exclusive processes. One such example is the beam-spin asymmetry (?LU) as a function of Trento-?, the angle between the leptonic plane and the hadronic plane which are de?ned with the beam and scattered electron, and with the recoiling neutron and reconstructed ?0, respectively. Despite the limited statistics, the observable ?LU was successfully extracted for eight kinematic bins (four in Mandelstam ?, and two in ??). This is a ?rst measurement for nDV?0P in the phase-space which is now available with CLAS12 in the 12 GeV era of JLab.

Naidoo, Paul↗
DOE OSTI2024

Computational Modeling, Simulation, and Potential Applications of Optical Stochastic Cooling

With the rising demand for intense particles beams, much research is being conducted in the area of particle beam cooling. One of these methods, called Stochastic Cooling (SC) (developed at CERN in the 1970's), delivered a feedback method to improve the quality and lifetime of circulating proton beams by reducing their 6D phase-space and has been widely implemented in a number of hadron machines. However, traditional stochastic cooling schemes are limited by the bandwidth of microwave frequency systems. Optical Stochastic Cooling (OSC) is a promising extension of the stochastic cooling beam cooling technique. OSC instead uses optical wavelengths which allows for improved control and increased cooling but creates its own technical challenges. This lays out work conducted toward the experimental demonstration of OSC at Fermilab's Integrable Optics Test Accelerator (IOTA) storage ring. This includes the design and characterization of parts of the optical delay system, the development and validation of a high-fidelity computational model of the OSC process, and the investigations into possible applications of the OSC mechanism to advanced beam manipulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Linac_Gen: integrating machine learning and particle-in-cell methods for enhanced beam dynamics at Fermilab

Here, we introduce Linac_Gen, a tool developed at Fermilab, which combines machine learning algorithms with Particle-in-Cell methods to advance beam dynamics in linacs. Linac_Gen employs techniques such as Random Forest, Genetic Algorithms, Support Vector Machines, and Neural Networks, achieving a tenfold increase in speed for phase-space matching in linacs over traditional methods through the use of genetic algorithms. Crucially, Linac_Gen's adept handling of 3D field maps elevates the precision and realism in simulating beam instabilities and resonances, marking a key advancement in the field. Benchmarked against established codes, Linac_Gen demonstrates not only improved efficiency and precision in beam dynamics studies but also in the design and optimization of linac systems, as evidenced in its application to Fermilab's PIP-II linac project. This work represents a notable advancement in accelerator physics, marrying ML with PIC methods to set new standards for efficiency and accuracy in accelerator design and research. Linac_Gen exemplifies a novel approach in accelerator technology, offering substantial improvements in both theoretical and practical aspects of beam dynamics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

First inclusive triple-differential measurement of the muon-antineutrino charged-current cross section using the NOvA Near Detector

We present the first measurement of the triple-differential muon antineutrino charged-current inclusive cross section, using the NOvA Near Detector and $12.5 \times 10^{20}$ protons on target in the NuMI beam. This sample of muon antineutrino interactions is the largest ever published, with approximately 1 million selected muon antineutrino events. The triple-differential cross section is measured in the final-state kinetic energy, the scattering angle, and the available energy of the interaction. The measurement enables phase-space regions populated by differing neutrino reaction processes to be isolated and the transition regions between them to be defined. The results are compared with the predictions of the main event generators used in the neutrino community and we observe energy- and angle-dependent discrepancies across a broad range of energies and interaction types

Abubakar, S. [Erciyes U.]↗
DOE OSTI2026

Characterizing the GD-1 Stream with DESI DR2 Data: Thin Stream and Hot Cocoon

GD-1 is among the longest, coldest stellar streams in the Milky Way, making it an ideal target for probing dark matter substructure through dynamical heating. We present a catalog of 608 spectroscopically confirmed GD-1 members from the first three years of Dark Energy Spectroscopic Instrument (DESI) observations. This constitutes the largest homogeneous spectroscopic sample of GD-1, doubling the number of members previously available only through heterogeneous compilations combining multiple surveys with different systematics. Using these data, we derive updated stream tracks in sky position, proper motion, and radial velocity that extend over $100^\circ$ of the stream. We apply a Gaussian mixture model to decompose the stream into a dynamically cold thin component ($σ_V = 2.49\pm 0.28$ km s$^{-1}$, width $= 0.23\pm0.01^\circ$) and a kinematically hot cocoon ($σ_V = 6.13\pm0.75$ km s$^{-1}$, width $= 2.18\pm0.17^\circ$). The cocoon contains $\sim30\%$ of members and its velocity dispersion is consistent with $\sim11$ Gyr of heating by cold dark matter subhalos. We also detect a large proper motion dispersion ($41.36\pm4.98$ km s$^{-1}$) along the stream direction in the cocoon component. This feature indicates a significant line-of-sight distance spread in the cocoon, and its origin will be further explored in a forthcoming paper. These measurements demonstrate the power of DESI spectroscopy for characterizing the multi-component phase-space structure of stellar streams and constraining small-scale dark matter substructure.

Jarvis, Emma [Toronto U.] (ORCID:0009000656127336)↗
DOE OSTI2026

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.

Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗
DOE OSTI2026

Double-real corrections to color singlet decay in a parton-shower inspired scheme

We introduce a local infrared subtraction method for next-to-next-to-leading order QCD calculations in color singlet decays, with counterterms based on scalar radiators and pure splitting functions. Overlapping singularities in the multipole radiation pattern are disentangled by partial fractioning, and the kinematics mapping corresponds to iterated next-to-leading order kinematics. We verify that the double-real remainder to $e^+e^-\to\;q\bar{q}$ is rendered finite in the single and double unresolved limits and investigate the numerical convergence of the Monte-Carlo integral. We compute the phase-space integrals of the scalar counterterms in the back-to-back configuration, both analytically and with the help of numerical techniques based on sector decomposition.

Campbell, John M. [Fermilab]↗
DOE OSTI2026

Strong-Drive Limits in Josephson Circuits: From Chaos to an Unbound-Resonance Threshold

Strong microwave drives enable fast measurement and parametric control in superconducting circuits but can induce transitions out of the intended low-energy manifold. We develop a unified description of strong-drive limits in flux- and charge-driven Josephson circuits across drive frequency and dc flux bias. Using classical phase-space analysis and Floquet--Markov simulations, we identify distinct low- and high-frequency mechanisms. At low frequency, we characterize bound-state resonances and separatrix chaos and find that the flux-drive chaos threshold depends strongly on dc flux bias. At high frequency, these mechanisms are suppressed, and the dissipative steady state transfers from the central bound-state sector to outer resonances formed by above-barrier running trajectories. The resulting unbound-resonance threshold is nearly independent of drive frequency and circuit parameters over the regime studied and is controlled primarily by dc flux bias. Coherent simulations show that parametric operation persists beyond this threshold, but at a reduced rate, setting an effective upper bound on the achievable operation speed. We derive analytical criteria for both thresholds, validate them numerically, and experimentally confirm the predicted dc-bias dependence of the low-frequency threshold in a flux-driven SQUID. We also determine the timescales for transfer into the unbound-resonance regime and relaxation back to the bound-state manifold after the drive is removed. Finally, we relate the stability limits to a complementary picture based on the junction critical current and extend the framework to multitone drives and inductively shunted circuits. Together, these results identify the mechanisms limiting strong driving and suggest routes to extend the stable operating range of Josephson circuits.

You, Xinyuan [Fermilab] (ORCID:0000000291789419)↗
DOE OSTI2025

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Dynamics of McMillan mappings I. McMillan multipoles

In this article, we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings, originally proposed by Edwin McMillan. Both represent the simplest symmetric McMillan maps, characterized by a single intrinsic parameter. While these systems find numerous applications across various domains of mathematics and physics, some of their dynamical properties remain unexplored. We aim to bridge this gap by providing a comprehensive description of all stable trajectories, including the parametrization of invariant curves, Poincaré rotation numbers, and canonical action–angle variables. In the second part, we establish connections between these maps and general chaotic maps in standard form. Our investigation reveals that the McMillan sextupole and octupole serve as first-order approximations of the dynamics around the fixed point, akin to the linear map and quadratic invariant (known as the Courant–Snyder invariant in accelerator physics), which represents zeroth-order approximations (referred to as linearization). Furthermore, we propose a novel formalism for nonlinear Twiss parameters, which accounts for the dependence of rotation number on amplitude. This stands in contrast to conventional betatron phase advance used in accelerator physics, which remains independent of amplitude. Notably, in the context of accelerator physics, this new formalism demonstrates its capability in predicting dynamical aperture around low-order resonances for flat beams, a critical aspect in beam injection/extraction scenarios.

43 PARTICLE ACCELERATORS↗