Engineering Papers⌕ Search

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

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Momentum space AC Josephson effect and intervalley coherence in multilayer graphene

Electron transport driven by the phase coherence and interference of quantum many-body wavefunctions is a fascinating phenomenon with potential technological significance. Superconductivity, for example, enables dissipationless transport through macroscopic phase twisting, while charge-density waves show broadband noise and AC-DC interference once their phase is depinned. We identify an analogous phase-coherent dynamics in the intervalley-coherent (IVC) state frequently observed in rhombohedral multilayer graphene. We show that a static magnetic field drives an oscillating intervalley current, which produces an oscillating orbital magnetization and thereby a detectable AC Hall response. This mechanism is similar to the Josephson effect observed in superconductors but now happening in momentum space. In this analogy, the static magnetic field acts as the DC voltage, while the oscillating intervalley current is the AC Josephson current. We provide microscopic calculations of all the phase-number free-energy parameters in rhombohedral trilayer graphene, establishing this effect as an experimental signature of IVC order.

electronic properties and materials↗
DOE OSTI2021

Tailored computational approaches to interrogate heavy element chemistry and structure in condensed phase

In this chapter we are presenting a brief review of the challenges encountered in the study of 4f and 5f block elements in the condensed phase. Their recovery, use in molten salt reactors and other interesting applications necessitate the use of molecular dynamics and large-scale models that take into account both the electronic structure and relativistic corrections. Sampling the multitude of electronic and atomic configurational states is at the heart of reliable predictions of structure, reactivity, dynamics and transport of heavy metals in complex environments. We present three examples that combine lanthanide elements with large scale models: i) our recent developments of a versatile adaptive learning method that enables global optimization in high dimensional spaces, ii) results of computed pKa values of lanthanide aqua complexes, and iii) structure and computed EXAFS of heavy elements in molten salts. The latter two employ our recently-optimized lanthanide pseudo-potentials and companion basis sets for condensed phase ab initio molecular dynamics.

Nguyen, Manh Thuong↗
DOE OSTI2023

Phase transition amplification of proton number fluctuations in nuclear collisions from a transport model approach

The time evolution of particle number fluctuations in nuclear collisions at intermediate energies (E lab = 1.23–10⁢A GeV) is studied by means of the UrQMD-3.5 transport model. The transport description incorporates baryonic interactions through a density-dependent potential. This allows for an implementation of a first-order phase transition including a mechanically unstable region at large baryon density. The scaled variance of the baryon and proton number distributions is calculated in the central cubic spatial volume of the collisions at different times. A significant enhancement of fluctuations associated with the unstable region is observed. This enhancement persists to late times, reflecting a memory effect for the fluctuations. The presence of the phase transition has a much smaller influence on the observable event-by-event fluctuations of protons in momentum space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Radiation Effects in III-V Solar Cells Grown by Dynamic Hydride Vapor Phase Epitaxy

The recent development of hydride vapor phase epitaxy (HVPE) is potentially promising as a route to lower the cost of high-efficiency III-V solar cells for space applications. HVPE produces the same materials and device structures as the industry-standard organometallic vapor phase epitaxy (OMVPE) process, and although HVPE has shown promising device efficiencies, it lags the OMVPE growth technique in technical maturity. For example, there are no reports of the performance of HVPE-grown devices in radiation environments. There is an expectation that high-quality (that is, single crystal and low-defect) III-V materials will behave similarly regardless of the growth method. However, it is important to verify assumptions as facts. In this project, we produced materials and devices using both HVPE and OMVPE at NREL and irradiated them using 1 MeV electrons to ascertain the effect of radiation on the materials. We also performed initial work on producing more radiation-hard structures, including devices with the pn junction at the front of the device, and devices with a graded doping profile, both of which help with radiation tolerance. Measurements of the devices post radiation exposure showed that the open-circuit voltage (V OC ) of the HVPE-grown solar cells degraded least on a percent basis, although they also started from a lower baseline than the OMVPE-grown devices. All solar cells reached approximately the same V OC after irradiation, indicating that the dose used was sufficient to degrade all device equally. The short-circuit current in the HVPE devices did degrade more than the other solar cells, and this was attributed to a higher-than-expected doping density in the base layer of that cell. The results of these experiments, while by no means comprehensive, do not show any material difference in the radiation effects in OMVPE- and HVPE-grown materials and devices.

14 SOLAR ENERGY↗
DOE OSTI2021

A Kinetic Monte Carlo Approach for Simulating Cascading Transmission Line Failure

In this work, cascading transmission line failures are studied through a dynamical model of the power system operating under fixed conditions. The power grid is modeled as a stochastic dynamical system where first-principles electromechanical dynamics are excited by small Gaussian disturbances in demand and generation around a specified operating point. In this context, a single line failure is interpreted in a large deviation context as a first escape event across a surface in phase space defined by line security constraints. The resulting system of stochastic differential equations admits a transverse decomposition of the drift, which leads to considerable simplification in evaluating the quasipotential (rate function) and, consequently, computation of exit rates. Tractable expressions for the rate of transmission line failure in a restricted network are derived from large deviation theory arguments and validated against numerical simulations. Extensions to realistic settings are considered, and individual line failure models are aggregated into a Markov model of cascading failure inspired by chemical kinetics. Cascades are generated by traversing a graph composed of weighted edges representing transitions to degraded network topologies. Numerical results indicate that the Markov model can produce cascades with qualitative power-law properties similar to those observed in empirical cascades.

Kinetic Monte Carlo, rare events, Cascading Failur↗
DOE OSTI2022

Cation reordering instead of phase transitions: Origins and implications of contrasting lithiation mechanisms in 1D ζ- and 2D α-V 2 O 5

Substantial improvements in cycle life, rate performance, accessible voltage, and reversible capacity are required to realize the promise of Li-ion batteries in full measure. Here, we have examined insertion electrodes of the same composition (V 2 O 5 ) prepared according to the same electrode specifications and comprising particles with similar dimensions and geometries that differ only in terms of their atomic connectivity and crystal structure, specifically two-dimensional (2D) layered α-V 2 O 5 that crystallizes in an orthorhombic space group and one-dimensional (1D) tunnel-structured ζ-V 2 O 5 crystallized in a monoclinic space group. By using particles of similar dimensions, we have disentangled the role of specific structural motifs and atomistic diffusion pathways in affecting electrochemical performance by mapping the dynamical evolution of lithiation-induced structural modifications using ex situ scanning transmission X-ray microscopy, operando synchrotron X-ray diffraction measurements, and phase-field modeling. We find the operation of sharply divergent mechanisms to accommodate increasing concentrations of Li-ions: a series of distortive phase transformations that result in puckering and expansion of interlayer spacing in layered α-V 2 O 5 , as compared with cation reordering along interstitial sites in tunnel-structured ζ-V 2 O 5 . By alleviating distortive phase transformations, the ζ-V 2 O 5 cathode shows reduced voltage hysteresis, increased Li-ion diffusivity, alleviation of stress gradients, and improved capacity retention. The findings demonstrate that alternative lithiation mechanisms can be accessed in metastable compounds by dint of their reconfigured atomic connectivity and can unlock substantially improved electrochemical performance not accessible in the thermodynamically stable phase.

25 ENERGY STORAGE↗
DOE OSTI2022

Filming movies of attosecond charge migration in single molecules with high harmonic spectroscopy

Electron migration in molecules is the progenitor of chemical reactions and biological functions after light-matter interaction. Following this ultrafast dynamics, however, has been an enduring endeavor. Here we demonstrate that, by using machine learning algorithm to analyze high-order harmonics generated by two-color laser pulses, we are able to retrieve the complex amplitudes and phases of harmonics of single fixed-in-space molecules. These complex dipoles enable us to construct movies of laser-driven electron migration after tunnel ionization of N 2 and CO 2 molecules at time steps of 50 attoseconds. Moreover, the angular dependence of the migration dynamics is fully resolved. By examining the movies, we observe that electron holes do not just migrate along the laser polarization direction, but may swirl around the atom centers. Our result establishes a general scheme for studying ultrafast electron dynamics in molecules, paving a way for further advance in tracing and controlling photochemical reactions by femtosecond lasers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

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↗
DOE OSTI2024

Protons on the line

Stable regions in four-dimensional phase space were observed by following the motion of accelerated proton beams subject to nonlinear forces. Here, this provides insights into the physics of dynamical systems and may lead to improved accelerator designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

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↗
DOE OSTI2024

Structural Evolution and Photoluminescence Quenching across the FASnI 3– x Br x ( x = 0–3) Perovskites

One of the primary methods for band gap tuning in metal halide perovskites has been halide (I/Br) mixing. Despite widespread usage of this type of chemical substitution in perovskite photovoltaics, there is still little understanding of the structural impacts of halide alloying, with the assumption being the formation of ideal solid solutions. The FASnI 3–x Br x (x = 0–3) family of compounds provides the first example where the assumption breaks down, as the composition space is broken into two unique regimes (x = 0–2.9; x = 2.9–3) based on their average structure with the former having a 3D and the latter having an extended 3D (pseudo 0D) structure. Pair distribution function (PDF) analyses further suggest a dynamic 5s 2 lone pair expression resulting in increasing levels of off-centering of the central Sn as the Br concentration is increased. These antiferroelectric distortions indicate that even the x = 0–2.9 phase space behaves as a nonideal solid-solution on a more local scale. Solid-state NMR confirms the difference in local structure yielding greater insight into the chemical nature and local distributions of the FA + cation. In contrast to the FAPbI 3–x Br x series, a drastic photoluminescence (PL) quenching is observed with x ≥ 1.9 compounds having no observable PL. In conclusion, our detailed studies attribute this quenching to structural transitions induced by the distortions of the [SnBr 6 ] octahedra in response to stereochemically expressed lone pairs of electrons. This is confirmed through density functional theory, having a direct impact on the electronic structure.

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

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radiative classical gravitational observables at $ \mathcal{O} $(G 3 ) from scattering amplitudes

We compute classical gravitational observables for the scattering of two spinless black holes in general relativity and N =8 supergravity in the formalism of Kosower, Maybee, and O’Connell (KMOC). We focus on the gravitational impulse with radiation reaction and the radiated momentum in black hole scattering at $ \mathcal{O} $(G 3 ) to all orders in the velocity. These classical observables require the construction and evaluation of certain loop-level quantities which are greatly simplified by harnessing recent advances from scattering amplitudes and collider physics. In particular, we make use of generalized unitarity to construct the relevant loop integrands, employ reverse unitarity, the method of regions, integration-by-parts (IBP), and (canonical) differential equations to simplify and evaluate all loop and phase-space integrals to obtain the classical gravitational observables of interest to two-loop order. The KMOC formalism naturally incorporates radiation effects which enables us to explore these classical quantities beyond the conservative two-body dynamics. From the impulse and the radiated momentum, we extract the scattering angle and the radiated energy. Finally, we discuss universality of the impulse in the high-energy limit and the relation to the eikonal phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Magnetic charge propagation upon a 3D artificial spin-ice

Magnetic charge propagation in spin-ice materials has yielded a paradigm-shift in science, allowing the symmetry between electricity and magnetism to be studied. Recent work is now suggesting the spin-ice surface may be important in mediating the ordering and associated phase space in such materials. Here, we detail a 3D artificial spin-ice, which captures the exact geometry of bulk systems, allowing magnetic charge dynamics to be directly visualized upon the surface. Using magnetic force microscopy, we observe vastly different magnetic charge dynamics along two principal directions. For a field applied along the surface termination, local energetics force magnetic charges to nucleate over a larger characteristic distance, reducing their magnetic Coulomb interaction and producing uncorrelated monopoles. In contrast, applying a field transverse to the surface termination yields highly correlated monopole-antimonopole pairs. Detailed simulations suggest it is the difference in effective chemical potential as well as the energy landscape experienced during dynamics that yields the striking differences in monopole transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Computational Exploration of High Entropy Alloys as Promising Materials for Future Beam Windows

With the ever-increasing demand for high beam power, the currently used beam-intercepting devices (BIDs) such as targets, and beam windows may not be able to handle the high power required for future accelerator complexes or the lifetime may be reduced drastically. As beam power increases, the damage incurred by BIDs, including thermal shock, fatigue, and irradiation damage, also rises. Therefore, it is imperative to design materials that can withstand high beam power for longer lifetimes. High entropy alloys (HEAs) have emerged as potential alternative materials for designing next-generation BIDs. In this study, we primarily focus on materials for developing beam windows for next-generation accelerator complexes. We propose an integrated approach that combines various computational techniques to study and design new materials. Specifically, we use CALPHAD, density functional theory (DFT), and molecular dynamics (MD) to comprehensively investigate the defect properties of suitable HEAs, offering potential alternatives for future beam windows. We begin by scanning the extensive phase space provided by Cr-Mn-V-Ti-Al-Co HEAs, selecting 8 compositions after evaluating approximately 120,000 unique compositions using CALPHAD. We, then employ DFT-informed machine learning techniques to develop force-field parameters. Finally, MD simulations using these developed force-field parameters will be used to study the effects of radiation damage on the defect and mechanical properties of the selected alloys. This research explains the use of the CALPHAD approach and shows how critical modeling (DFT and MD) is in developing novel material such as HEAs. It also highlights the promising role of machine learning in this field. The results from this study will greatly improve the novel materials development to be used in next-generation accelerator components, leading to higher beam power and longer operational times of BIDs.

43 PARTICLE ACCELERATORS↗