Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space methods”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evaluation protocol for revealing magnonic contrast in TR-STXM measurements

We present a statistically motivated method to extract magnonic contrast from time-resolved scanning transmission x-ray microscopy (TR-STXM) measurements. TR-STXM is an element-specific method for resolving spin-dynamics in space and time. It offers nanometer spatial resolution and picosecond temporal resolution. The presented method makes it possible to obtain phase and amplitude profiles of spin-waves from STXM measurements. Furthermore, it allows for a rigorous transformation to reciprocal magnon k⃗-space, revealing k⃗-dependent magnon properties such as the magnon dispersion in three dimensions and for all directions of the magnetic anisotropy. We demonstrate our method using X-band ferromagnetic resonance on a micrometer-sized permalloy assembly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Outer Radiation Belt Dynamics During the October 2012 Storm Revisited: Rapid Inward Radial Transport From a Dynamic Outer Boundary

Earth's outer radiation belt electron flux is highly variable and can be enhanced by over an order of magnitude over timescales less than one day, as observed during the October 2012 storm. Previous studies of this storm (e.g., Reeves et al., 2013, https://doi.org/10.1126/science.1237743) have invoked local acceleration to explain this. However, here, we argue that the observations can instead be explained by fast inward radial transport. One method often invoked to distinguish between these two acceleration processes is the existence of local peaks in electron phase space density (PSD) as a function of L* at fixed first, M, and second, K, adiabatic invariants. However, this method relies on the assumption that the evolution of the PSD as a function of L* occurs over timescales slower than the satellite orbital period. Here, high spatiotemporal resolution data from the Global Positioning System (GPS) spacecraft constellation is used to show that enhancements in the PSD occur during the October 2012 storm over short timescales not resolvable by the Van Allen Probes. In addition, Geostationary Operational Environmental Satellite spacecraft data also indicate that these enhancements are consistent with relativistic electron injections. A radial diffusion model is shown to reproduce the PSD dynamics observed by the Van Allen Probes, once rapid variations at the simulation outer boundary are included, consistent with GPS data. This verifies that apparently “locally growing” peaks in PSD along high apogee satellite orbits can be produced by fast inward radial transport without requiring the action of any local acceleration processes.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

A Jacobian-free pseudo-arclength continuation method for phase transitions in inhomogeneous thermodynamic systems

Developing phase diagrams for inhomogeneous systems in thermodynamics is difficult, in part, due to the large phase space and the possibility of unstable and metastable solutions arising from first-order phase transitions. Pseudo-arclength continuation (PAC) is a method that allows one to trace out stable and unstable solutions of nonlinear systems. Typically, PAC utilizes the Jacobian in order to implement Newton (or quasi-Newton) steps. In this work, we present a Jacobian-free PAC method that is amenable to the usual workflows in inhomogeneous thermodynamics. We demonstrate our method in systems that have first-order phase transitions, including a novel example of polyelectrolyte complex coacervation in confinement, where multiple surface phase transitions occur and can overlap with one another.

Chemistry↗
DOE OSTI2025

Preliminary study of auto-differentiation algorithm in beam dynamics with stochastic process

Modern particle accelerator optimization requires sophisticated computational methods to address the inherently stochastic nature of beam dynamics. This research develops a framework applying AD to SDEs that specifically addresses beam dynamics challenges in particle accelerators, focusing on accurately modeling and optimizing beam behavior in regimes dominated by stochastic processes. By incorporating key physical phenomena such as synchrotron radiation, wakefield effects, and quantum excitation, the framework aims to provide auto differentiation on the figure of merit of the phase space evolution and beam dynamics. The methodology will enable effective optimization method in a dynamic system with stochastic process.

Accelerator Physics↗
DOE OSTI2022

Fock-Space Schrieffer–Wolff Transformation: Classically-Assisted Rank-Reduced Quantum Phase Estimation Algorithm

We present an extension of many-body downfolding methods to reduce the resources required in the quantum phase estimation (QPE) algorithm. In this paper, we focus on the Schrieffer–Wolff (SW) transformation of the electronic Hamiltonians for molecular systems that provides significant simplifications of quantum circuits for simulations of quantum dynamics. We demonstrate that by employing Fock-space variants of the SW transformation (or rank-reducing similarity transformations (RRST)) one can significantly increase the locality of the qubit-mapped similarity-transformed Hamiltonians. The practical utilization of the SW-RRST formalism is associated with a series of approximations discussed in the manuscript. In particular, amplitudes that define RRST can be evaluated using conventional computers and then encoded on quantum computers. The SW-RRST QPE quantum algorithms can also be viewed as an extension of the standard state-specific coupled-cluster downfolding methods to provide a robust alternative to the traditional QPE algorithms to identify the ground and excited states for systems with various numbers of electrons using the same Fock-space representations of the downfolded Hamiltonian. The RRST formalism serves as a design principle for developing new classes of approximate schemes that reduce the complexity of quantum circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Neutral beam prompt loss in LTX-β

Prompt loss of beam injected fast ions approaches 100% in lithium tokamak experiment-beta (LTX-β) discharges, though significantly improved confinement is expected for the higher current plasmas made available by a recent upgrade to the Ohmic heating power supply. Modeling of fast ions using TRANSP/NUBEAM finds a maximum coupled beam fraction of 76% at the near-term limits of the LTX-β operating space. The full ion orbit code POET is employed to validate NUBEAM results against possible non-adiabatic effects on fast ion orbits, but corrections to the prompt loss fraction due to collisionless transport are found to be small. The graphical method code CONBEAM is used to investigate the topology of fast ion phase space as it relates to neutral beam deposition, and counter-injected NBI is considered as a way to access a region of high field side beam deposition. A metric is developed within the CONBEAM using a beam filament model to estimate the prompt loss fraction and shown to agree well with both POET and NUBEAM, enabling near real-time analysis and potential feedback to operators between plasma discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Probing structure–property relationship in chemical vapor deposited hybrid perovskites by pressure and temperature

Chemical vapor deposition (CVD), a low-cost and a scalable deposition technique, allows the growth of methylammonium lead iodide (MAPbI 3 ) films without the use of solvents, substantially increasing air stability while also inducing the stable cubic phase at room temperature and at pressures as low as 0.25 GPa. MAPbI 3 thin films were grown by a facile two-step low-pressure vapor deposition process in a single reactor. This method results in films, which are usually in the tetragonal phase (space group: I4/mcm) and occasionally in the cubic phase under ambient conditions. High-pressure synchrotron-based X-ray diffraction studies from CVD-grown MAPbI 3 crystallites show that the sample remains in the cubic phase (space group: Im$\bar{3}$) between 0.25 and 3.0 GPa. Temperature-dependent transport measurements show sharp anomalies, correlating with the structural changes. The transport measurements from the CVD- grown cubic MAPI 3 film is further compared with a film in the tetragonal phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Possibilities and Limitations of Kinematically Identifying Stars from Accreted Ultra-faint Dwarf Galaxies

Abstract The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ( E , L z , J r , J z ) is one method of identifying accreted structure that is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the Caterpillar simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that, even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover 6% of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, z accretion ≲ 0.5. Based on these results, we make recommendations to help guide the search for dynamically linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and J r , providing a way to help identify real versus fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Comparison of heavy-ion transport simulations: Mean-field dynamics in a box

Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared among each other and to exact calculations. For BUU-like codes, employing the test particle method, the results depend on the combination of the number of test particles and the spread of the profile functions that weight integration over space. These parameters can be properly adapted to give a good reproduction of the analytical zero-sound features. QMD-like codes, using molecular dynamics methods, are characterized by large damping effects, attributable to the fluctuations inherent in their phase-space representation. Moreover, for a given nuclear effective interaction, they generally lead to slower density oscillations, as compared to BUU-like codes. The latter problem is mitigated in the more recent lattice formulation of some of the QMD codes. Furthermore, the significance of these results for the description of real heavy-ion collisions is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Computational discovery of ultra-strong, stable, and lightweight refractory multi-principal element alloys. Part II: comprehensive ternary design and validation

Abstract Here the discovery of refractory multi-principal element alloys (MPEAs) with high-temperature strength and stability is pursued within a constrained and application-relevant design space. A comprehensive approach is developed and applied to explore all 165 ternary systems in the Al-Ce-Fe-Hf-Mo-Nb-Ta-Ti-V-W-Zr family. A subset of ternary systems that contain large areas in composition–temperature space with high strength and robust BCC phase stability is found. Twelve sets of high-performing alloys are identified, each set optimized for one combination of phase constraint, optimization target, and temperature range. Preliminary mechanical tests support the viability of the method. This work highlights the importance of considering phase stability, exploring non-equiatomic regions of composition space, and applying application-relevant constraints. Parts I and II provide three down-selection techniques for identifying high-performing BCC refractory MPEAs, design guidelines, and many candidates predicted to have BCC phase stability and strengths 2–3 times higher than any reported to date.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Classification of fragile topology enabled by matrix homotopy

Flat bands in twisted materials have attracted considerable attention due to the emergence of correlated phases that can be associated with the non-Wannier-representable nature of its single-particle states. Specifically, these bands can exhibit a class of topology that can be nullified by the addition of trivial bands, termed fragile topology, which has required an expansion of prior classification schemes. However, existing approaches for predicting fragile topology rely on momentum-space methods, e.g., Wilson loops, presenting a fundamental challenge for using fragile topology as a predictor of correlated phases in aperiodic systems, such as incommensurate twist angles in moiré materials. Here, we develop a ℤ 2 energy-resolved topological marker for classifying fragile phases using a system’s position-space description, enabling the direct classification of finite, disordered, and aperiodic materials. By translating the physical symmetries protecting the system’s fragile topological phase into matrix symmetries of the system’s Hamiltonian and position operators, we use matrix homotopy to construct our topological marker while simultaneously yielding a quantitative measure of topological robustness. We demonstrate our framework’s effectiveness in both a low-energy tight-binding model and a continuum photonic crystal model of 𝐶 2 ⁢𝒯-symmetric systems, and find that fragile topology can both persist under strong disorder and even exhibit disorder-induced reentrant phase transitions. Our photonic crystal results also demonstrate the robustness of fragile topology, and the applicability of our approach, to heterostructures lacking a bulk spectral gap. Overall, our framework serves as an efficient tool for elucidating fragile topology, offering guidance for the prediction and discovery of correlated phases in both crystalline and aperiodic materials.

Lee, Ki Young [Sandia National Laboratories (SNL-N↗
DOE OSTI2021

Measurements of $\overline{\nu}_{\mu}$ and $\overline{\nu}_{\mu} + \nu_{\mu}$ charged-current cross-sections without detected pions or protons on water and hydrocarbon at a mean anti-neutrino energy of 0.86 GeV

We report measurements of the flux-integrated $\overline{\nu}_\mu$ and $\overline{\nu}_\mu+\nu_\mu$ charged-current cross-sections on water and hydrocarbon targets using the T2K anti-neutrino beam with a mean beam energy of 0.86 GeV. The signal is defined as the (anti-)neutrino charged-current interaction with one induced $\mu^\pm$ and no detected charged pion or proton. These measurements are performed using a new WAGASCI module recently added to the T2K setup in combination with the INGRID Proton Module. The phase space of muons is restricted to the high-detection efficiency region, $p_{\mu}>400~{\rm MeV}/c$ and $\theta_{\mu}<30^{\circ}$, in the laboratory frame. An absence of pions and protons in the detectable phase spaces of $p_{\pi}>200~{\rm MeV}/c$, $\theta_{\pi}<70^{\circ}$ and $p_{\rm p}>600~{\rm MeV}/c$, $\theta_{\rm p}<70^{\circ}$ is required. In this paper, both the $\overline{\nu}_\mu$ cross-sections and $\overline{\nu}_\mu+\nu_\mu$ cross-sections on water and hydrocarbon targets and their ratios are provided by using the D’Agostini unfolding method. The results of the integrated $\overline{\nu}_\mu$ cross-section measurements over this phase space are $\sigma_{\rm H_{2}O}=(1.082\pm0.068(\rm stat.)^{+0.145}_{-0.128}(\rm syst.)) \times 10^{-39}\,{\rm cm^{2} / nucleon}$, $\sigma_{\rm CH}=(1.096\pm0.054(\rm stat.)^{+0.132}_{-0.117}(\rm syst.)) \times 10^{-39}\,{\rm cm^{2} / nucleon}$, and $\sigma_{\rm H_{2}O}/\sigma_{\rm CH} = 0.987\pm0.078(\rm stat.)^{+0.093}_{-0.090}(\rm syst.)$. The $\overline{\nu}_\mu+\nu_\mu$ cross-section is $\sigma_{\rm H_{2}O} = (1.155\pm0.064(\rm stat.)^{+0.148}_{-0.129}(\rm syst.)) \times 10^{-39}\,{\rm cm^{2} / nucleon}$, $\sigma_{\rm CH}=(1.159\pm0.049(\rm stat.)^{+0.129}_{-0.115}(\rm syst.)) \times 10^{-39}\,{\rm cm^{2} / nucleon}$, and $\sigma_{\rm H_{2}O}/\sigma_{\rm CH}=0.996\pm0.069(\rm stat.)^{+0.083}_{-0.078}(\rm syst.)$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Topological terms with qubit regularization and relativistic quantum circuits

Qubit regularization provides a rich framework to explore quantum field theories. The freedom to choose how the important symmetries of the theory are embedded in the qubit regularization scheme allows us to construct new lattice models with rich phase diagrams. Some of the phases can contain topological terms which lead to critical phases. In this work we introduce and study the SU(3)-F qubit regularization scheme to embed the SO(3) spin-symmetry. We argue that qubit models in this regularization scheme contain several phases including a critical phase which describes the k = 1 Wess-Zumino-Witten (WZW) conformal field theory (CFT) at long distances, and two massive phases one of which is trvially gapped and the other which breaks the lattice translation symmetry. We construct a simple space-time Euclidean lattice model with a single coupling U and study it using the Monte Carlo method. We show the model has a critical phase at small U and a trivially massive phase at large U with a first order transition separating the two. Another feature of our model is that it is symmetric under space-time rotations, which means the temporal and spatial lattice spacing are connected to each other. The unitary time evolution operator obtained by a Wick rotation of the transfer matrix of our model can help us compute the physics of the k = 1 WZW CFT in real time without the need for tuning the temporal lattice spacing to zero. We use this idea to introduce the concept of a relativistic quantum circuit on a discrete space-time lattice.

Bhattacharya, Tanmoy↗
DOE OSTI2024

Efficient six-dimensional phase space reconstructions from experimental measurements using generative machine learning

Next-generation accelerator concepts, which hinge on the precise shaping of beam distributions, demand equally precise diagnostic methods capable of reconstructing beam distributions within six-dimensional position-momentum spaces. However, the characterization of intricate features within six-dimensional beam distributions using current diagnostic techniques necessitates a substantial number of measurements, using many hours of valuable beam time. Novel phase space reconstruction techniques are needed to reduce the number of measurements required to reconstruct detailed, high-dimensional beam features in order to resolve complex beam phenomena and as a feedback in precision beam shaping applications. In this study, we present a novel approach to reconstructing detailed six-dimensional phase space distributions from experimental measurements using generative machine learning and differentiable beam dynamics simulations. We demonstrate that this approach can be used to resolve six-dimensional phase space distributions from scratch, using basic beam manipulations and as few as 20 two-dimensional measurements of the beam profile. We also demonstrate an application of the reconstruction method in an experimental setting at the Argonne Wakefield Accelerator, where it is able to reconstruct the beam distribution and accurately predict previously unseen measurements 75× faster than previous methods.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2023

Adaptive autoencoder latent space tuning for more robust machine learning beyond the training set for six-dimensional phase space diagnostics of a time-varying ultrafast electron-diffraction compact accelerator

In this work, we present a general adaptive latent space tuning approach for improving the robustness of machine learning tools with respect to time variation and distribution shift. We demonstrate our approach by developing an encoder-decoder convolutional neural network-based virtual 6D phase space diagnostic of charged particle beams in the HiRES ultrafast electron diffraction (UED) compact particle accelerator with uncertainty quantification. Our method utilizes model-independent adaptive feedback to tune a low dimensional 2D latent space representation of ~1 million dimensional objects which are the 15 unique 2D projections (x, y),...,(z, p z ) of the 6D phase space (x, y, z, p x , p y , p z ) of the charged particle beams. We demonstrate our method with numerical studies of short electron bunches utilizing experimentally measured UED input beam distributions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

A machine learning study on spinodal clumping in heavy ion collisions

Possible observables of baryon number clustering due to the instabilities occurring at a first order QCD phase transition are discussed. The dynamical formation of baryon clusters at a QCD phase transition can be described by numerical fluid dynamics, augmented with a gradient term and an equation of state with a mechanically unstable region. It is shown that the dynamical description of this phase transition, in nuclear collisions, will lead to the formation of dense baryon clusters at the phase boundary. State-of-the-art machine learning methods find that the coordinate space clumping leaves characteristic imprints on the spatial net density distribution in almost every event. On the other hand the momentum distributions do not show any clear event-by-event features. Lastly, it is shown that the 'third order' cumulant, the skewness, shows a peak at the beam energy where the system, created in the heavy ion collision, reaches the deconfinement phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗