Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Phase space dynamics”

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 73 records · Page 4

Coherence-Induced Deep Thermalization Transition in Random Permutation Quantum Dynamics

We report a phase transition in the projected ensemble—the collection of postmeasurement wave functions of a local subsystem obtained by measuring its complement. The transition emerges in systems undergoing random permutation dynamics, a type of quantum time evolution wherein computational basis states are shuffled without creating superpositions. It separates a phase exhibiting deep thermalization, where the projected ensemble is distributed over Hilbert space in a maximally entropic fashion (Haar random), from a phase where it is minimally entropic (“classical bit-string ensemble”). Crucially, this deep thermalization transition is invisible to the subsystem’s density matrix, which always exhibits thermalization to infinite temperature across the phase diagram. Through a combination of analytical arguments and numerical simulations, we show that the transition is tuned by the total amount of injected by the input state and the measurement basis, and is exhibited robustly across different microscopic models. Our findings represent a novel form of ergodicity-breaking universality in quantum many-body dynamics, characterized not by a failure of regular thermalization, but rather by a failure of deep thermalization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Learning plasma dynamics and robust rampdown trajectories with predict-first experiments at TCV

The rampdown phase of a tokamak pulse is difficult to simulate and often exacerbates multiple plasma instabilities. To reduce the risk of disrupting operations, we leverage advances in Scientific Machine Learning (SciML) to combine physics with data-driven models, developing a neural state-space model (NSSM) that predicts plasma dynamics during Tokamak à Configuration Variable (TCV) rampdowns. The NSSM efficiently learns dynamics from a modest dataset of 311 pulses with only five pulses in a reactor-relevant high-performance regime. The NSSM is parallelized across uncertainties, and reinforcement learning (RL) is applied to design trajectories that avoid instability limits. High-performance experiments at TCV show statistically significant improvements in relevant metrics. A predict-first experiment, increasing plasma current by 20% from baseline, demonstrates the NSSM’s ability to make small extrapolations. The developed approach paves the way for designing tokamak controls with robustness to considerable uncertainty and demonstrates the relevance of SciML for fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Advanced blade-shaped thermal energy storage device: Development and application

Thermal energy storage (TES) using phase change materials (PCMs) is a promising approach for capturing and reusing excess thermal energy, yet widespread adoption is limited by low thermal conductivity, bulky configurations, and inadequate scalability. Here, this study presents a modular, blade-shaped TES prototype designed to address these challenges. The device integrates a lightweight aluminum shell, an embedded serpentine coil for active or passive heat exchange, and a cost-effective corrugated metal mesh for enhanced PCM thermal conductivity. With thickness-to-length and thickness-to-width ratios of 0.03 and 0.08, respectively, the blade-shaped TES achieves a compact, modular form factor suitable for space-constrained applications. Experimental testing demonstrated the efficient charge and discharge behavior of blade-shaped TES, capturing PCM superheating, phase-change transitions, and subcooling dynamics, with charging and discharging efficiencies of 94.9% and 94.6%, respectively. Also, the system can potentially achieve higher energy density than that of conventional TES designs. When integrated into a household refrigerator during the study, three blade-shaped TES modules successfully shifted 100% of peak-time compressor operation to off-peak hours, reducing energy consumption while maintaining more stable compartment temperatures. The blade-shaped TES's thin geometry, modularity, and enhanced thermal performance support scalable deployment across residential, commercial, and industrial applications, providing a versatile, cost-effective solution for high-efficiency, demand-flexible thermal energy management.

Blade-shaped

Primordial black holes and magnetic fields in conformal neutrino mass models

Sufficiently strong and long-lasting first-order phase transitions can produce primordial black holes (PBHs) that contribute substantially to the dark matter abundance of the Universe, and can produce large-scale primordial magnetic fields. Here, we study these mechanisms in a generic class of conformal U(1)' models that also explain active neutrino oscillation data via the type-I seesaw mechanism. We find that phase transitions that occur at seesaw scales between 10 4 GeV and 10 11 GeV produce gravitational wave signals (from the dynamics of the phase transition and from the decay of cosmic string loops) at LISA/ET that can be correlated with microlensing signals of PBHs at the Roman Space Telescope, while scales near 10 11 GeV can be correlated with Hawking evaporation signals at future gamma-ray telescopes. LISA can probe the entire range of PBH masses between 1 × 10 -16 M ⊙ and 8 × 10 -11 M ⊙ if PBHs fully account for the dark matter abundance. For Z' masses between 40 TeV and 10 4 TeV, and 10 TeV right-handed neutrinos, helical magnetic fields can be produced with magnitudes ≳ 0.5 pG and coherence lengths ≳ 0.008 Mpc, above current blazar lower bounds.

cosmological phase transitions

Quantum fisher information reveals UV-IR mixing in the strange metal

The density-density response in optimally doped Bi 2 Sr 2 CaCu 2 O 8+$x$ has recently been shown to exhibit conformal symmetry. Using, the experimentally inferred conformal dynamic susceptibility, we compute the resultant quantum Fisher information (QFI), a witness to multi-partite entanglement. For a Fermi liquid, we find that the QFI grows quadratically as the temperature increases, consistent then with the phase space available for scattering in the standard theory of metals. By contrast, the QFI in a strange metal increases as a power law at as the temperature decreases, but ultimately extrapolates to a constant at T=0. The constant is of the form, ω$^{2Δ}_g$, where Δ is the conformal dimension and ω g is the UV cutoff which is on the order of the pseudogap. As this constant depends on both UV and IR properties, it illustrates that multipartite entanglement in a strange metal exhibits UV-IR mixing, a benchmark feature of doped Mott insulators as exemplified by dynamical spectral weight transfer. We conclude with a discussion of the implication of our results for low-energy reductions of the Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Pressure-tuned plethora of ferroelectric phases in CuInP 2 S 6

Layered chalcogenides like CuInP 2 S 6 are room temperature ferroelectrics. Modest compression even increases the electric polarization, raising questions about the origin of this unusual trend and other properties under pressure. In this work, we combine synchrotron-based infrared and Raman scattering spectroscopies, x-ray diffraction, and resistivity measurements with first-principles calculations of the lattice dynamics and energy landscape to unravel the influence of pressure on symmetry, polarization, and metallicity. We uncover a remarkable series of phase transitions across a series of polar space groups: monoclinic Cc → trigonal P31c (prismatic sulfur) → P31c (octahedral sulfur). True metallicity develops above 63 GPa, significantly higher than in related MPS 3 materials (M = Mn, Co, Fe, Ni), offering a picture of competing states of matter that is different than previously supposed. Detailed examination of pressure trends within the Cc phase also reveals phonon lifetime changes and streaking of satellite x-ray peaks that correlate with the maximum polarization. We discuss these tendencies in terms of Cu + ion migration, phase formation, and the overall energy landscape. Our findings place the high pressure behavior of CuInP 2 S 6 on a firm foundation and pave the way for the development of structure-property relations in this family of complex chalcogenides.

ferroelectric phases

Probing topological phases in a perturbed Kane-Mele model via RKKY interaction: Application to monolayer jacutingaite Pt 2 ⁢HgSe 3

Quantum spin Hall insulators (QSHIs) leverage strong spin-orbit coupling (SOC) for efficient spin manipulation, making them promising for spintronics. Here, in this study, we investigate the noncollinear Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between two magnetic impurities in a perturbed Kane-Mele model with strong SOC, relevant to monolayer jacutingaite Pt 2 ⁢HgSe 3 as a prominent QSHI. Following the previous studies that mainly focused on the model and its general applications, we provide a systematic examination of the effects of various perturbations and strong spin-orbit hybridizations, which drive phase transitions that have not been extensively explored before. By incorporating these perturbations into the model and accurately accounting for spin-orbit hybridizations through spin-space Green's functions and the RKKY interactions, we uncover distinct, relative (rather than absolute) signatures of different phase transitions. These phase transitions are induced by both static and dynamic perturbations on the magnetic impurities. Notably, we identify additional phases emerging from the interplay with the magnetic substrate. All these influence the switching between ferromagnetic and antiferromagnetic, as well as clockwise and counterclockwise magnetic interactions. Our results provide a practical way to track topological phases through magnetic properties, offering new insights into phase control and spin manipulation in QSHIs.

Kane-Mele model

Dynamic realization of emergent high-dimensional optical vortices

The dimensionality of vortical structures has recently been extended beyond two dimensions, providing additional topological complexity and robustness for high-capacity information processing and turbulence control. The generation of high-dimensional vortical structures has mostly been demonstrated in classical systems through the complex interference of fluidic, acoustic, or electromagnetic waves. However, natural materials rarely support three- or higher-dimensional vortical structures and their physical interactions. Here, we experimentally demonstrate a high-dimensional gradient thickness optical cavity (GTOC) in which the optical coupling of planar metal-dielectric multilayers implements topological interactions across multiple dimensions. At non-trivial topological phases, high-dimensional GTOC induces high-dimensional vortical structures in generalized parameter space in three, four dimensions, and beyond. These emergent high-dimensional vortical structures are observed under electro-optic tomography as optical vortex dynamics in two-dimensional real-space, employing the optical thicknesses of dielectric layers as synthetic dimensions. Our findings hold significant promise for emulating high-dimensional physics and developing active topological photonic devices.

36 MATERIALS SCIENCE

New Muon Campus Simulations for the Muon $g\textrm{-}2$ Experiment at Fermilab

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab conducted high-precision measurements of the muon anomalous magnetic moment $a_{\mu}$ using a storage ring from 2018 to 2023, achieving a remarkable precision of $200\:\mathrm{ppb}$ over Runs 1-3, with analyses for Runs 4-6 ongoing. A comprehensive understanding of the storage ring's beam dynamics and its accurate simulations are crucial for achieving the experiment's ambitious goals. One of the requirements for this effort is a very detailed knowledge of the phase space distribution of the beam. To address this requirement, we performed high-statistics simulations of the Muon $g\textrm{-}2$ Target Station (AP0) and the Muon Campus beamlines: M2 and M3, followed by the Delivery Ring, and then M4 and M5. The resulting muon distribution at the end of the M5 beamline from our previous $3\times10^{12}$ protons-on-target (PoT) simulation serves as an essential input for the storage ring simulations. In 2024, to facilitate the analyses of Runs 4-6, we have updated our Muon Campus models and re-optimised certain parameters to reflect the operational currents and wire chamber measurements of the beam. For these optimisations, we employed the heterogeneous island method, implemented in our evolutionary optimisation tool, glyfada. This key update addresses the need to use the best possible beam for the storage ring simulations and thus supports the experiment's overall precision. The Muon $g\textrm{-}2$ Experiment at Fermilab has successfully achieved its $70\:\mathrm{ppb}$ systematic uncertainty goal and collected 21 times more data than its predecessor at BNL. The updated and improved Muon Campus models and simulations not only facilitate the experiment's efforts to potentially resolve the current tension between experimental measurements and theoretical predictions of $a_{\mu}$, but also provide a basis for future simulations for the Mu2e Experiment (E-973), which utilises shared Muon Campus beamlines.

43 PARTICLE ACCELERATORS

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING

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

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

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

accelerator

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

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Content-aware foveated camera for multi-target tracking

Modern image sensors deliver substantial space-time bandwidth, yet indiscriminate acquisition often overwhelms memory, computation, and downstream perception. We present a content-aware, multi-foveated camera that dynamically reallocates sensing and magnification to multiple regions of interest (ROIs). A phase-only spatial light modulator (SLM) serves as a solid-state, inertia-free beam-steering and lens element, enabling per-frame field-of-view (FOV) reconfiguration and content-aware target tracking. By interleaving frames across foveae, our system preserves a wide-FOV situational context while refreshing each ROI at high rates, thereby reducing data volume without degrading task performance. We constructed a prototype employing a single-SLM, single-sensor architecture and demonstrated its application in real-time multi-object tracking with dynamic ROI maintenance across multiple viewpoints. The approach offers a general pathway to integrate detection, tracking, and segmentation algorithms in the acquisition loop, shifting workload from post hoc processing to intelligent capture.

Zang, Zihan [Univ. of California, Los Angeles, CA