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 91 records · Page 5

Extraction of Drell-Yan Angular Parameters in $pp$ Collisions with a 120 GeV Beam Energy Using a Deep-Learning Unfolding Algorithm

Dilepton production in pp collisions through the Drell-Yan process provides a crucial tool for studying the internal quark-gluon structure of the nucleon. By precisely measuring the $\cos2\phi$ asymmetry, where $\phi$ represents the azimuthal angle of the $l^{+}l^{-}$ pair in the Collins-Soper frame, we can gain valuable insights into the proton’s structure and the transverse momentum ($q_{T}$) dependence of the $\cos2\phi$ asymmetry. SeaQuest, a fixed-target Drell-Yan experiment at Fermilab, involved an unpolarized proton beam colliding with unpolarized LH$_{2}$ and LD$_{2}$ targets. Measurements obtained from experiments typically require corrections for detector inefficiencies, smearing, and acceptance. Traditionally, these corrections involve “unfolding” the detector-level measurements through matrix operations. However, in higher-dimensional phase space, these conventional methods fail to scale effectively. To overcome these limitations, we employ an unbinned unfolding method that utilizes deep neural networks for unfolding higher-dimensional phase space. In this presentation, we will explain the design of the neural network architecture, our training strategies, and outline our plans to achieve conclusive results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗

Shortened Booster Bunch Lengths

Booster typically produces beam bunches that have a phase space profile with a low momentum spread but high time spread, correlating to longer bunch lengths. The goal of this study is to produce beam bunches that have a phase space profile with a high momentum spread but low time spread or narrow bunch lengths. This can be achieved by slowly increasing the RF cavities amplitudes (Adiabatic Excitation) or by oscillating the amplitudes (Quadratic Bunch Rotation, QBR) until the desired bunch lengths are achieved. Results show that when shortening the bunches using QBR, it showed the bunches had a length of 0.6ns (within 1 sigma) but beam losses were being incurred in the 8 GeV line which were a result of beam dispersion.

Abdelhamid, Maan↗

Shortened Booster Bunch Lengths

Booster typically produces beam bunches that have a phase space profile with a low momentum spread but high time spread, correlating to longer bunch lengths. The goal of this study is to produce beam bunches that have a phase space profile with a high momentum spread but low time spread or narrow bunch lengths. This can be achieved by slowly increasing the RF cavities amplitudes (Adiabatic Excitation) or by oscillating the amplitudes (Quadratic Bunch Rotation, QBR) until the desired bunch lengths are achieved. Results show that when shortening the bunches using QBR, it showed the bunches had a length of 0.6ns (within 1 sigma) but beam losses were being incurred in the 8 GeV line which were a result of beam dispersion.

Abdelhamid, Maan↗

Uncertainty-Driven Rapid Thermodynamic Assessment of Nb-Ta-Zr System and Effects of C impurities (L25GF9298S): Annual Progress Report

An integrated computational materials engineering (ICME) method is in development for rapid thermodynamic experimental investigation and high-fidelity computational modeling of refractory multi-principal element alloys (RMPEAs). These ultra-high-temperature (UHT) alloys are of interest for structural applications in extreme environments, but deficiency of reliable data, especially melting temperatures, impedes the prediction of alloys with favorable properties. The method leverages UHT capabilities and computational expertise of LLNL’s Materials Science Division and the McCormack Lab’s UHT conical nozzle levitation (CNL) system to iteratively map the Nb-Ta-Zr phase space, with focus on the liquidus surface, through targeted experiments selected by quantifying uncertainty in the thermodynamic model fitting parameters. This method will reduce the time to map uncharted RMPEA phase space and thereby accelerate discovery and development of advanced materials for applications in extreme environments.

36 MATERIALS SCIENCE↗

Conditional guided generative diffusion for particle accelerator beam diagnostics

Abstract Advanced accelerator-based light sources such as free electron lasers (FEL) accelerate highly relativistic electron beams to generate incredibly short (10s of femtoseconds) coherent flashes of light for dynamic imaging, whose brightness exceeds that of traditional synchrotron-based light sources by orders of magnitude. FEL operation requires precise control of the shape and energy of the extremely short electron bunches whose characteristics directly translate into the properties of the produced light. Control of short intense beams is difficult due to beam characteristics drifting with time and complex collective effects such as space charge and coherent synchrotron radiation. Detailed diagnostics of beam properties are therefore essential for precise beam control. Such measurements typically rely on a destructive approach based on a combination of a transverse deflecting resonant cavity followed by a dipole magnet in order to measure a beam’s 2D time vs energy longitudinal phase-space distribution. In this paper, we develop a non-invasive virtual diagnostic of an electron beam’s longitudinal phase space at megapixel resolution (1024 × 1024) based on a generative conditional diffusion model. We demonstrate the model’s generative ability on experimental data from the European X-ray FEL.

43 PARTICLE ACCELERATORS↗

First Order Transport of a Cold Uniform Ellipsoid of Charge

A single particle has a 6-element phase space vector $s = (x, γv)$ where $γv =\frac{p}{m}$ s proportional to the relativistic momentum. The ellipsoid is represented by a 21-element vector $(s, C, V)$ where $s$ is the phase space location of the centroid. $C$ is the 3 × 3 covariance matrix defined by $C_{ij} = ⟨δx_iδx_j⟩$, where $δx$ is the distance of a particle from the centroid $x$. $C$ is symmetric so only counts as 6 elements of the overall vector. $V$ is a general 3 × 3 matrix that determines the velocity distribution in the bunch by $δ(γv) = Vδx$. In total this makes 6 + 6 + 9 = 21 elements.

43 PARTICLE ACCELERATORS↗

Active control of Alfvén eigenmodes by external magnetic perturbations with different spatial spectra

Alfvén eigenmodes have been suppressed and excited in tokamak plasmas by (just) modifying the poloidal spectra of externally applied static magnetic perturbations. This effect is observed experimentally when toroidal spectra of n = 2, n = 4 as well as a mixed spectrum of n = 2 and n = 4 is applied. Under the n = 2 magnetic perturbations, the modes are excited or suppressed by modifying the coil phasing between the upper and the lower set of coils. Regardless of the absolute rotation, an even parity for the n = 4 perturbation is observed to reduce the amplitude of the Alfvénic instabilities, while an odd parity amplifies it. To combine the stabilizing (and destabilizing) effect of n = 2 and n = 4, a mixed spectrum is applied, finding similar reduction (and amplification) trends. However, the impact on the mode amplitude is more subtle, due to the reduced coil current required for a mixed spectrum. The signal level on the fast-ion loss detector is sensitive to the applied poloidal spectrum, which is consistent with Hamiltonian full-orbit modelling of an edge resonant transport layer activated by the 3D perturbative fields. An internal redistribution of the fast-ion population is induced, modifying the phase-space gradients driving the Alfvénic instabilities, and ultimately determining their existence. The calculated edge resonant layers for both n = 2 and n = 4 toroidal spectra are consistent with the observed suppressed and excited phases. Moreover, hybrid kinetic-magnetohydrodynamic (MHD) simulations reveal that this edge resonant transport layer overlaps in phase-space with the population responsible for the fast-ion drive. The results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

Alfvén waves↗

Forecasting generative amplification

Generative networks are perfect tools to enhance the speed and precision of LHC simulations. Especially when generating events beyond the size of the training dataset, it is important to understand their statistical precision. We present two complementary methods to estimate the amplification factor without large holdout datasets. Averaging amplification uses Bayesian networks or ensembling to estimate amplification from the precision of integrals over given phase-space volumes. Differential amplification uses hypothesis testing to quantify amplification without any resolution loss. Applied to state-of-the-art event generators, both methods indicate that amplification is already possible in specific regions of phase space.

Bahl, Henning [Heidelberg Univ. (Germany)] (ORCID:↗

Measurement of the inclusive WZ production cross section in pp collisions at $\sqrt{\textrm{s}}$ = 13.6 TeV

The inclusive WZ production cross section is measured in proton-proton collisions at a centre-of-mass energy of 13.6 TeV, using data collected during 2022 with the CMS detector, corresponding to an integrated luminosity of 34.7 fb −1 . The measurement uses multileptonic final states and a simultaneous likelihood fit to the number of events in four different lepton flavour categories: eee, ee μ , μμ e, and μμμ . The selection is optimized to minimize the number of background events, and relies on an efficient prompt lepton discrimination strategy. The WZ production cross section is measured in a phase space defined within a 30 GeV window around the Z boson mass, as σ total (pp → WZ) = 55.2 ± 1.2 (stat) ± 1.2 (syst) ± 0.8 (lumi) ± 0.3 (theo) pb. In addition, the cross section is measured in a fiducial phase space closer to the detector-level requirements. All the measurements presented in this paper are in agreement with standard model predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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 ↗

Giant quantum oscillations in thermal transport in low-density metals via electron absorption of phonons

Oscillations of conductance observed in strong magnetic fields are a striking manifestation of the quantum dynamics of charge carriers in solids. The large charge carrier density in typical metals sets the scale of oscillations in both electrical and thermal conductivity, which characterize the Fermi surface. In semimetals, thermal transport at low-charge carrier density is expected to be phonon dominated, yet several experiments observe giant quantum oscillations in thermal transport. This raises the question of whether there is an overarching mechanism leading to sizable oscillations that survives in phonon-dominated semimetals. In this work, we show that such a mechanism exists. It relies on the peculiar phase-space allowed for phonon scattering by electrons when only a few Landau levels are filled. Our measurements on the Dirac semimetal ZrTe 5 support this counterintuitive mechanism through observation of pronounced thermal quantum oscillations, since they occur in similar magnitude and phase in directions parallel and transverse to the magnetic field. Our phase-space argument applies to all low-density semimetals, topological or not, including graphene and bismuth. Our work illustrates that phonon absorption can be leveraged to reveal degrees of freedom through their imprint on longitudinal thermal transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗

Transverse Kinematic Imbalance in MicroBooNE's New nue CC0pi Measurements

Neutrino-nucleus cross section measurements require accurate modelling of neutrino interactions. Neutrino beams are not monoenergetic, and the energy of each interaction must instead be modelled using nuclear interaction assumptions. This introduces significant systematic uncertainty into cross section measurements. Effects such as Fermi motion, nuclear correlations, and final-state interactions (FSI) smear the underlying quasi-elastic scattering signal, making it difficult to disentangle genuine quasi-elastic kinematics from nuclear effects across the full range of interaction channels (QE, MEC, RES, DIS) probed in these measurements. Transverse Kinematic Imbalance (TKI) variables, such as $\delta p_T$ and $\delta \alpha_T$, probe this same phase space by exploiting the fact that the incoming neutrino has zero transverse momentum ($\vec{p}_T^{\,\nu} = 0$). Any measured transverse imbalance in the final state therefore arises from nuclear effects rather than from uncertainty in the incident neutrino energy, allowing cross section measurements to select a phase space that is rich in quasi-elastic-like events with minimal contamination from FSI and other nuclear effects, independent of energy reconstruction. Recent unfolded MicroBooNE cross section measurements of electron-neutrino charged-current interactions with zero pions and at least one proton ($\nu_e$ CC0$\pi$, 1eNp0$\pi$) show that several leading nuclear interaction generators (including GENIE variants, NuWro, GiBUU, and NEUT) reproduce the differential cross section in electron energy reasonably well, but consistently struggle to describe the differential cross section in the cosine of the leading proton's angle, yielding lower $p$-values across all seven generators tested. This tension points to a more fundamental, kinematics-driven disagreement between data and generators that is not visible in energy-only cross section observables. This is precisely the regime TKI variables are designed to probe. Following previous TKI cross section measurements with muon-neutrino data in MicroBooNE, this poster presents the case for extending the TKI framework to electron-neutrino cross section measurements as a next step to isolate and characterize the source of the observed generator tension in proton kinematics.

Burridge, Jessica [U. Manchester (main)] (ORCID:00↗

Third integer resonant extraction transit time simulation studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab]↗

Third Integer Resonant Extraction Transit Time Simulation Studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

Measurement of C P Violation Observables in D + → K − K + π + Decays

A search for violation of the charge-parity ( C P ) symmetry in the D + → K − K + π + decay is presented, with proton-proton collision data corresponding to an integrated luminosity of 5.4 fb − 1 , collected at a center-of-mass energy of 13 TeV with the LHCb detector. A novel model-independent technique is used to compare the D + and D − phase-space distributions, with instrumental asymmetries subtracted using the D s + → K − K + π + decay as a control channel. The p value for the hypothesis of C P conservation is 8.1%. The C P asymmetry observables A C P | S ϕ π + = ( 0.95 ± 0.4 3 stat ± 0.2 6 syst ) × 10 − 3 and A C P | S K ¯ * 0 K + = ( − 0.26 ± 0.5 6 stat ± 0.1 8 syst ) × 10 − 3 are also measured. These results show no evidence of C P violation and represent the most sensitive search performed through the phase space of a multibody decay. © 2024 CERN, for the LHCb Collaboration 2024 CERN

Aaij, R. (ORCID:0000000305331952)↗

Confined and deconfined chaos in classical spin systems

Weakly perturbed integrable many-body systems are typically chaotic, and thermal at late times. However, there are distinct relationships between the timescales for thermalization and chaos. The typical relationship is confined chaos: when trajectories are still confined to regions in phase space with constant conserved quantities (actions), the conjugate angle variables are already unstable. Chaotic instabilities thus far precede thermalization. In a different relationship, which we term deconfined chaos, chaotic instabilities and thermalization occur on the same timescale. We investigate these two qualitatively distinct scenarios through numerical and analytical studies of two perturbed integrable classical spin models: the Ishimori spin chain (confined chaos), and the central spin model with XX interactions (deconfined chaos). We analytically establish (super)-integrability in the latter model in a microcanonical shell. Deconfined chaos emerges through the separation of phase space into large quasi-integrable regions and a thin chaotic manifold. The latter leads to chaos and thermalization on the fastest possible timescale, which is proportional to the inverse perturbation strength. This behavior is reminiscent of the quantum SYK models and strange metals.

Chaotic Dynamics (nlin.CD)↗