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 415 records · Page 23

The Nosé–Hoover looped chain thermostat for low temperature thawed Gaussian wave-packet dynamics

We have used a generalised coherent state resolution of the identity to map the quantum canonical statistical average for a general system onto a phase-space average over the centre and width parameters of a thawed Gaussian wave packet. We also propose an artificial phase-space density that has the same behaviour as the canonical phase-space density in the low-temperature limit, and have constructed a novel Nosé–Hoover looped chain thermostat that generates this density in conjunction with variational thawed Gaussian wave-packet dynamics. This forms a new platform for evaluating statistical properties of quantum condensed-phase systems that has an explicit connection to the time-dependent Schrödinger equation, whilst retaining many of the appealing features of path-integral molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A complex systems analysis of stick-slip dynamics of a laboratory fault

We study the stick-slip behavior of a granular bed of photoelastic disks sheared by a rough slider pulled along the surface. Time series of a proxy for granular friction are examined using complex systems methods to characterize the observed stick-slip dynamics of this laboratory fault. Nonlinear surrogate time series methods show that the stick-slip behavior appears more complex than a periodic dynamics description. Phase space embedding methods show that the dynamics can be locally captured within a four to six dimensional subspace. These slider time series also provide an experimental test for recent complex network methods. Phase space networks, constructed by connecting nearby phase space points, proved useful in capturing the key features of the dynamics. In particular, network communities could be associated to slip events and the ranking of small network subgraphs exhibited a heretofore unreported ordering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum Mechanics of a Spherically Symmetric Causal Diamond in Minkowski Spacetime

We construct the phase space of a spherically symmetric causal diamond in ( d + 2 )-dimensional Minkowski spacetime. Utilizing the covariant phase space formalism, we identify the relevant degrees of freedom that localize to the d -dimensional bifurcate horizon and, upon canonical quantization, determine their commutators. On this phase space, we find two Iyer-Wald charges. The first of these charges, proportional to the area of the causal diamond, is responsible for shifting the null time along the horizon and has been well documented in the literature. The second charge is much less understood, being integrable for d ≥ 2 only if we allow for field-dependent diffeomorphisms and is responsible for changing the size of the causal diamond. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Adaptive machine learning for time-varying systems: low dimensional latent space tuning

Machine learning (ML) tools such as encoder-decoder convolutional neural networks (CNN) can represent incredibly complex nonlinear functions which map between combinations of images and scalars. For example, CNNs can be used to map combinations of accelerator parameters and images which are 2D projections of the 6D phase space distributions of charged particle beams as they are transported between various particle accelerator locations. Despite their strengths, applying ML to time-varying systems, or systems with shifting distributions, is an open problem, especially for large systems for which collecting new data for re-training is impractical or interrupts operations. Particle accelerators are one example of large time-varying systems for which collecting detailed training data requires lengthy dedicated beam measurements which may no longer be available during regular operations. We present a novel method of adaptive ML for time-varying systems. Our approach is to map very high (N ≈ 100k) dimensional inputs (a combination of scalar parameters and images) into the low dimensional (N ≈ 2) latent space at the output of the encoder section of an encoder-decoder CNN. We then actively tune the low dimensional latent space-based representation of complex system dynamics by the addition of an adaptively tuned feedback vector directly before the decoder sections builds back up to our image-based high-dimensional phase space density representations. This method allows us to learn correlations within and to quickly tune the characteristics of incredibly large parameter space systems and to track their evolution in real time based on feedback without massive new data sets for re-training. We demonstrate that our method can accurately predict and track the phase space of charged particle beams at various locations in a particle accelerator by adaptively adjusting in real-time while the unknown input beam distribution of the accelerator is changing in shape, charge, and offset and while the RF system of the accelerator itself is also changing in an unpredictable way. For FACET-II we demonstrate that such an approach has the potential to use transverse deflecting cavity and energy spread spectrum beam measurements to accurately predict 2D projections of the 6D phase space of the electron beam at the plasma wakefield acceleration interaction point where such diagnostics are unavailable.

47 OTHER INSTRUMENTATION↗

General guide concepts for compact, high-brilliance neutron moderators

The trend in neutron sciences is toward integrating compact, high-brightness moderators into new or upgraded facilities. Transporting neutrons from the source to the sample position with a phase-space distribution tailored to specific requirements is crucial to leverage high source brilliance. Here we have investigated four guide concepts using Monte Carlo ray tracing simulations: Montel beamline with nested Kirkpatrick–Baez mirrors, curved-tapered beamline with a bender and straight sections, straight-elliptical beamline, and curved-elliptical beamline. The straight-elliptical (curved-elliptical) beamline features two half-ellipse guides connected by a straight (non-straight) guide section. The neutron transport efficiency and phase space homogeneity have been quantitatively compared. Our results show that the straight-elliptical beamline performs best because of few neutron bounces on the guide surface with small reflection angles, minimizing flux loss. The Montel beamline provides the best spatial confinement of neutrons within the desired region; however, there is a high thermal-neutron loss due to large reflection angles. The curved-tapered beamline suffers from significant flux loss due to high bounces, and it shows a non-uniform angular distribution related to broad ranges of bounces and reflection angles. The non-straight guide section of the curved-elliptical beamline increases the phase space inhomogeneity, leading to a spatially non-uniform beam profile. The results apply to general neutron instruments that require transporting thermal and cold neutrons from a compact, high-brilliance moderator to the sample location with a moderate phase-space volume.

47 OTHER INSTRUMENTATION↗

Topological Langmuir-cyclotron wave

A theory is developed to describe the topological Langmuir-cyclotron wave (TLCW), a topological excitation in magnetized plasmas recently identified by numerical simulations. As a topological wave in a continuous medium, the TLCW propagates unidirectionally without scattering in complex boundaries and can be explored as an effective mechanism to energize particles. We show that, because momentum space in continuous media is contractible in general, the topology of the wave bundles is trivial over momentum space that contains no degeneracy points. This is in stark contrast to condensed matters with periodic lattice structures that impose nontrivial topology on momentum space. In continuous media without lattice structures, nontrivial topology of the eigenmode bundles manifests over phase space, and it is the nontrivial topology over phase space that underpins the topological excitations, such as the TLCW. It is shown that the TLCW can be faithfully modeled by a generic tilted Dirac cone in phase space, whose entire spectrum, including the spectral flow, is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic corrections from analytic non-Maxwellian distribution functions in magnetized plasmas

In magnetized plasma physics, almost all developed analytic theories assume a Maxwellian distribution function (MDF) and in some cases small deviations are described using the perturbation theory. The deviations with respect to the Maxwellian equilibrium, called kinetic effects, are required to be taken into account especially for fusion reactor plasmas. Generally, because the perturbation theory is not consistent with observed steady-state non-Maxwellians, these kinetic effects are numerically evaluated by very central processing unit (CPU)-expensive codes, avoiding the analytic complexity of velocity phase space integrals. We develop here a new method based on analytic non-Maxwellian distribution functions constructed from non-orthogonal basis sets in order to (i) use as few parameters as possible, (ii) increase the efficiency to model numerical and experimental non-Maxwellians, (iii) help to understand unsolved problems such as diagnostics discrepancies from the physical interpretation of the parameters, and (iv) obtain analytic corrections due to kinetic effects given by a small number of terms and removing the numerical error of the evaluation of velocity phase space integrals. This work does not attempt to derive new physical effects even if it could be possible to discover one from the better understandings of some unsolved problems, but here we focus on the analytic prediction of kinetic corrections from analytic non-Maxwellians. As applications, examples of analytic kinetic corrections are shown for the secondary electron emission, the Langmuir probe characteristic curve, and the entropy. This is done by using three analytic representations of the distribution function: the Kappa distribution function, the bi-modal or a new interpreted non-Maxwellian distribution function (INMDF). The existence of INMDFs is proved by new understandings of the experimental discrepancy of the measured electron temperature between two diagnostics in JET. As main results, it is shown that (i) the empirical formula for the secondary electron emission is not consistent with a MDF due to the presence of super-thermal particles, (ii) the super-thermal particles can replace a diffusion parameter in the Langmuir probe current formula, and (iii) the entropy can explicitly decrease in presence of sources only for the introduced INMDF without violating the second law of thermodynamics. Moreover, the first order entropy of an infinite number of super-thermal tails stays the same as the entropy of a MDF. The latter demystifies the Maxwell's demon by statistically describing non-isolated systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Adaptive model tuning studies for non-invasive diagnostics and feedback control of plasma wakefield acceleration at FACET-II

Beam-driven plasma wakefield acceleration (PWFA) achieves the same energy gain in a single meter, for which conventional accelerators require several kilometers, however much work is still required to match the beam quality of conventional accelerators. The PWFA processes to be studied at the FACET-II facility will utilize extremely short (down to a few fs) high energy (10 GeV), high charge (few nC) electron and eventually positron bunches. The PWFA process is extremely sensitive to the detailed longitudinal current profiles of these bunches and it would be of great benefit to have precise measurement and control of these profiles. We present an adaptive model tuning technique for FACET-II to adaptively tune online models based on real time accelerator and beam data to continuously provide a non-invasive diagnostic which predicts the longitudinal phase space (LPS) of extremely short and highly compressed electron beams, which otherwise require destructive X-band transverse deflecting cavity (XTCAV)-based measurements. Based on simulation studies, our method has the potential for: (1). The development of a non-invasive longitudinal phase space diagnostic by adaptively tuning models based on non-invasive measurements such as energy spread spectra which could be recorded at all of the bunch compressors in the facility. (2). Utilize these diagnostic to perform model-independent feedback-based to achieve desired longitudinal phase space distributions. (3). Utilize the model-independent feedback approach to maximize energy gain while minimizing emittance growth and energy gain variance of the PWFA process by tuning accelerator parameters, while monitoring what longitudinal phase space the algorithm has found, information that will be useful for further analytical and simulation studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Verification of 5D continuum gyrokinetic code COGENT: Studies of kinetic drift wave instability

COGENT (Continuum Gyrokinetic Edge New Technology) is a kinetic plasma simulation code that is being developed by the edge simulation laboratory (ESL) collaboration. The original version of the code has been developed in a 4D phase space (2D configuration space and 2D velocity space) to address kinetic plasma phenomena in a complex magnetic field geometry including the core, magnetic separatrix, and scrape-off layer regions. This paper is focused on extending the original 4D phase space to a 5D phase space (3D2V) to address full kinetic turbulences in a tokamak edge region. Here, we report on the current status of 5D COGENT, which presently operates in a shear-less, simple slab geometry. As a verification study, we use the problems of collision-less drift wave instability (universal instability) including its modification in the presence of collisional effects. The simulation model includes the gyrokinetic equations for ion and electron species coupled to the long-wavelength limit of the gyro-Poisson equation. Collisional effects are represented by the Krook collision model. Linear analytical results for the drift mode growth rate and real frequency are recovered, and the non-linear stage is modelled and analysed as well. In addition, extensive 5D runs have been performed to address the effects of the drift wave instability on blob/filamentary structures characteristic of a tokamak edge. Finally, a helical-shaped potential perturbation is observed to grow exponentially in time while spinning around the filament axis with electron drift frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inferring the gravitational potential of the Milky Way with a few precisely measured stars

The dark matter halo of the Milky Way is expected to be triaxial and filled with substructure. It is hoped that streams or shells of stars produced by tidal disruption of stellar systems will provide precise measures of the gravitational potential to test these predictions. We develop a method for inferring the Galactic potential with tidal streams based on the idea that the stream stars were once close in phase space. Our method can flexibly adapt to any form for the Galactic potential: it works in phase-space rather than action-space and hence relies neither on our ability to derive actions nor on the integrability of the potential. Our model is probabilistic, with a likelihood function and priors on the parameters. The method can properly account for finite observational uncertainties and missing data dimensions. We test our method on synthetic data sets generated from N-body simulations of satellite disruption in a static, multi-component Milky Way, including a triaxial dark matter halo with observational uncertainties chosen to mimic current and near-future surveys of various stars. We find that with just eight well-measured stream stars, we can infer properties of a triaxial potential with precisions of the order of 5%-7%. Without proper motions, we obtain 10% constraints on most potential parameters and precisions around 5%-10% for recovering missing phase-space coordinates. These results are encouraging for the goal of using flexible, time-dependent potential models combined with larger data sets to unravel the detailed shape of the dark matter distribution around the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-space Berry phases: Skyrmion soccer (invited)

Berry phases occur when a system adiabatically evolves along a closed curve in parameter space. This tutorial-like article focuses on Berry phases accumulated in real space. In particular, we consider the situation where an electron traverses a smooth magnetic structure, while its magnetic moment adjusts to the local magnetization direction. Mapping the adiabatic physics to an effective problem in terms of emergent fields reveals that certain magnetic textures, skyrmions, are tailormade to study these Berry phase effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity

In this report we compare time sampling and ensemble averaging as two different methods available for phase space sampling. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium molecular dynamics. We introduce two different schemes for the ensemble averaging approach, and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical molecular dynamics, the ensemble generation approaches may find their greatest utility in computationally expensive simulations such as first principles molecular dynamics. For such simulations, where each time step is costly, time sampling can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. On the other hand, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each ensemble is independent. For this reason, particularly when using massively parallel architectures, ensemble sampling can result in much shorter simulation times and exhibits similar overall computational effort.

74 ATOMIC AND MOLECULAR PHYSICS↗

Optimizing 4d Emittance Measurements Using the Pinhole Scan Technique

Accurate measurement of electron beam emittance is essential for optimizing high-brightness electron sources. The Pinhole Scan Technique measures the 4D phase space and hence the emittance by measuring the beam profile after clipping the beam using a pinhole followed by a drift section and then scanning the beam over the pinhole. This technique has been implemented in low energy (< 200 keV) beamlines at both Cornell University and Arizona State University. However, the technique poses several practical challenges. In this work, we analyze and address key issues affecting the 4D phase space and emittance measurements using this technique. We identify and investigate sources of inaccuracies like the pinhole aspect ratio, beam divergence, position-momentum correlations in the phase space, and the point-spread-function of the detector and suggest techniques to minimize them. Our findings offer a pathway to more accurate 4D phase space characterization in advanced electron beam systems.

42 ENGINEERING↗

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗

Crossing the dividing surface of transition state theory. III. Once and only once. Selecting reactive trajectories

The purpose of the present work is to determine initial conditions that generate reacting, recrossing-free trajectories that cross the conventional dividing surface of transition state theory (i.e., the plane in configuration space passing through a saddle point of the potential energy surface and perpendicular to the reaction coordinate) without ever returning to it. Local analytical equations of motion valid in the neighborhood of this planar surface have been derived as an expansion in Poisson brackets. We show that the mere presence of a saddle point implies that reactivity criteria can be quite simply formulated in terms of elements of this series, irrespective of the shape of the potential energy function. Some of these elements are demonstrated to be equal to a sum of squares and thus to be necessarily positive, which has a profound impact on the dynamics. The method is then applied to a three-dimensional model describing an atom-diatom interaction. A particular relation between initial conditions is shown to generate a bundle of reactive trajectories that form reactive cylinders (or conduits) in phase space. This relation considerably reduces the phase space volume of initial conditions that generate recrossing-free trajectories. Loci in phase space of reactive initial conditions are presented. Reactivity is influenced by symmetry, as shown by a comparative study of collinear and bent transition states. Finally, it is argued that the rules that have been derived to generate reactive trajectories in classical mechanics are also useful to build up a reactive wave packet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Null geodesics and thermodynamic phase transition of four-dimensional Gauss–Bonnet AdS black hole

Highlights: • Modified gravity theories. • Null geodesics and the thermodynamic phase transition. • Four-dimensional Gauss–Bonnet AdS black hole. • Correlation between gravity and thermodynamics. Modified gravity theories are of great interest in both observational and theoretical studies. In this article, we study the correlation between the null geodesics and the thermodynamic phase transition of a four-dimensional Gauss–Bonnet AdS black hole. Firstly, we study the phase structure of the black hole, using the coexistence and spinodal curves, to understand the phase transition in an extended phase space. The imprints of these phase transition features can be clearly observed from the functional dependence of photon orbit radius and minimum impact parameter with respect to the Hawking temperature and pressure. The change in these two quantities during the phase transition serve as order parameters which characterise the critical behaviour. The correlation shows that thermodynamic phase transition can be studied by observing the effects of strong gravitational field on the photon orbit and vice versa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Asymptotic structure of higher dimensional Yang-Mills theory

Using the covariant phase space formalism, we construct the phase space for non-Abelian gauge theories in (d+2)-dimensional Minkowski spacetime for any d ≥ 2, including the edge modes that symplectically pair to the low energy degrees of freedom of the gauge field. Despite the fact that the symplectic form in odd and even-dimensional spacetimes appear ostensibly different, we demonstrate that both cases can be treated in a unified manner by utilizing the shadow transform. Upon quantization, we recover the algebra of the vacuum sector of the Hilbert space and derive a Ward identity that implies the leading soft gluon theorem in (d+2)-dimensional spacetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗