Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAJECTORIES”

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 541 records · Page 30

Solving Newton’s equations of motion with large timesteps using recurrent neural networks based operators

Classical molecular dynamics simulations are based on solving Newton’s equations of motion. Using a small timestep, numerical integrators such as Verlet generate trajectories of particles as solutions to Newton’s equations. We introduce operators derived using recurrent neural networks that accurately solve Newton’s equations utilizing sequences of past trajectory data, and produce energy-conserving dynamics of particles using timesteps up to 4000 times larger compared to the Verlet timestep. We demonstrate significant speedup in many example problems including 3D systems of up to 16 particles.

Newton’s equations↗

Geometric GNNs for charged particle tracking at GlueX

Nuclear physics experiments are aimed at uncovering the fundamental building blocks of matter. The experiments involve high-energy collisions that produce complex events with many particle trajectories. Tracking charged particles resulting from collisions in the presence of a strong magnetic field is critical to enable the reconstruction of particle trajectories and precise determination of interactions. It is traditionally achieved through combinatorial approaches that scale worse than linearly as the number of hits grows. Since particle hit data naturally form a point cloud and can be structured as graphs, graph neural networks (GNNs) emerge as an intuitive and effective choice for this task. In this study, we evaluate the GNN model for track finding on the data from the GlueX experiment at Jefferson Lab. We use simulation data to train the model and test on both simulation and real GlueX measurements. We demonstrate that GNN-based track finding outperforms the currently used traditional method at GlueX in terms of segment-based efficiency at a fixed purity while providing faster inferences. We show that the GNN model can achieve significant speedup by processing multiple events in batches, which exploits the parallel computation capability of graphical processing units (GPUs). Finally, we compare the GNN implementation on GPU and field-programmable gate array and describe the trade-off.

batched GNN pipeline↗

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗

Nonadiabatic strong-field photoionization revisited: Detailed analysis

We develop a theoretical framework to interpret “bump-like” enhancements in the energy distribution of electrons emitted from cesium atoms subjected to intense infrared laser fields, as observed in a recent experiment. These emissions, detected at energies approaching twice the ponderomotive energy (2⁢𝑈 p ), indicate a mechanism where electrons are nonadiabatically ejected with substantial inward velocity following ionization. The bump-like features arise from returning trajectories that scatter off the residual ion. We interpret these experimental results classically using a periodically driven version of the Bohr model and quantum mechanically through Keldysh's adiabaticity parameter (𝛾 K ) and Reiss's continuum-state intensity parameter 𝒵 in a Coulomb-corrected version of the strong-field approximation (SFA) that incorporates “complex” electronic trajectories. Notably, our theoretical analysis shows that these emissions are independent of the specific short-range form of the atomic potential.

74 ATOMIC AND MOLECULAR PHYSICS↗

Reshaping the Quantum Arrow of Time

While the microscopic laws of physics are often symmetric under time reversal, most natural processes that we observe are not. The emergent asymmetry between typical and time-reversed processes is referred to as the arrow of time. In quantum physics, an arrow of time emerges when a sequence of measurements is performed on a system. We introduce quantum control tools that can yield dynamics more consistent with time flowing backward than forward. The control tools are based on the explicit construction of a Hamiltonian that can replicate the stochastic trajectories of a monitored quantum system. Such a Hamiltonian can reverse the effect of monitoring and, via a feedback process, generate trajectories consistent with a reversed arrow of time. It can also be used to simulate the backward-in-time dynamics of an open quantum system. Finally, we design a feedback-driven continuous measurement engine powered by the energy pumped into the system by the monitoring process. We show that the engine can operate under experimentally realizable conditions with feedback delay and finite-efficiency measurements.

Entropy production↗

Anomalous ellipticity dependence of the generation of near-threshold harmonics in noble gases

We measured the ellipticity dependence of the harmonic yield of He, Ne, Ar, Kr, and Xe. It is found that the harmonic yield in general decreases monotonically with increasing laser ellipticity but anomalous dependence does occur for each atom for some harmonic orders. Additionally, we found that the degree of anomaly is weakest in helium but is stronger for the heavier atoms. The anomaly also depends on laser intensity. To explain these features, we developed a quantum trajectory Monte Carlo (QTMC) model based on the quantum path integral theory to study high-order harmonic generation (HHG). The model includes the effect of Coulomb potential from the target ion, thus it is capable of calculating near-threshold harmonics quantitatively. This model reveals that the presence of Coulomb potential would generate orbiting trajectories which are responsible for the anomalous ellipticity dependence near the threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗

Evolution of Efimov states

The Efimov phenomenon manifests itself as an emergent discrete scaling symmetry in the quantum three-body problem. In the unitarity limit, it leads to an infinite tower of three-body bound states with energies forming a geometric sequence. Here in this work, we study the evolution of these so-called Efimov states using relativistic scattering theory. We identify them as poles of the three-particle S matrix and trace their trajectories in the complex energy plane as they evolve from virtual states through bound states to resonances. We dial the scattering parameters toward the unitarity limit and observe the emergence of the universal scaling of energies and couplings - a behavior known from the non-relativistic case. Interestingly, we find that Efimov resonances follow unusual, cyclic trajectories accumulating at the three-body threshold and then disappear at some values of the two-body scattering length. We propose a partial resolution to this ?missing states? problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nudged elastic band approach to nuclear fission pathways

Here, we study the performance and precision of various methods to determine the minimum-action and minimum-energy fission trajectories in the collective space. To this end, we apply the nudged elastic band method (NEB), grid-based methods, and the Euler-Lagrange approach to the collective action minimization in two- and three-dimensional collective spaces.The NEB method is the tool of choice for finding the least-action and minimum-energy fission trajectories. It will be particularly useful in large-scale static fission calculations of superheavy nuclei and neutron-rich fissioning nuclei contributing to the astrophysical r-process recycling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Far-from-equilibrium search for the QCD critical point

Initial conditions for relativistic heavy-ion collisions may be far from equilibrium (i.e., there are large initial contributions from the shear stress tensor and bulk pressure) but it is expected that on very short time scales the dynamics converge to a universal attractor that defines hydrodynamic behavior. Thus far, studies of this nature have only considered an idealized situation at LHC energies (high temperatures T and vanishing baryon chemical potential µ B = 0) but, in this work, we investigate for the first time how far-from-equilibrium effects may influence experimentally driven searches for the Quantum Chromodynamic critical point at RHIC. We find that the path to the critical point is heavily influenced by far from equilibrium initial conditions where viscous effects lead to dramatically different {T, µ B } trajectories through the QCD phase diagram. We compare hydrodynamic equations of motion with shear and bulk coupled together at finite µB for both DNMR and phenomenological Israel-Stewart equations of motion and discuss their influence on potential attractors at finite µ B and their corresponding {T, µ B } trajectories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Particle production by a relativistic semitransparent mirror in (1+3)D Minkowski spacetime

Production of massless, scalar particles by a relativistic, semitransparent, plane mirror in (1+3)D Minkowski spacetime based on the Barton-Calogeracos (BC) action is investigated. The corresponding Bogoliubov coefficients are derived for a mirror with arbitrary, relativistic trajectories. We apply our derived formula to two specific trajectories. One is commonly used in the (1+1)D literature to mimic gravitational collapse theoretically, and the other is proposed to be realizable experimentally. In addition, we identify the relation between the particle spectrum and the particle production probability, and we demonstrate the equivalence between our approach and the existing approach in the literature, which is restricted to (1+1)D. In short, our treatment extends the study to (1+3)D spacetime for a relativistic, plane mirror. Finally, we offer a third approach for finding the particle spectrum using the S-matrix formalism.

79 ASTRONOMY AND ASTROPHYSICS↗

Superposed metric for spinning black hole binaries approaching merger

Here, we construct an approximate metric that represents the spacetime of spinning binary black holes (BBH) approaching merger. We build the metric as an analytical superposition of two Kerr metrics in harmonic coordinates, where we transform each black hole term with time-dependent boosts describing an inspiral trajectory. The velocities and trajectories of the boost are obtained by solving the post-Newtonian (PN) equations of motion at 3.5 PN order. We analyze the spacetime scalars of the new metric and we show that it is an accurate approximation of Einstein’s field equations in vacuum for a BBH system in the inspiral regime. Furthermore, to prove the effectiveness of our approach, we test the metric in the context of a 3D general relativistic magnetohydrodynamical (GRMHD) simulation of accreting minidisks around the black holes. We compare our results with a previous well-tested spacetime construction based on the asymptotic matching method. We conclude that our new spacetime is well-suited for long-term GRMHD simulations of spinning binary black holes on their way to the merger.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion energy distribution in an electron beam ion trap inferred from simulations of the trapped ion cloud

For this work, we have inferred the energy distribution of trapped ions in an electron beam ion trap (EBIT) from simulations of the spatial distribution of Fe 13+ ions and a comparison with measured visible light images of the ion cloud. We simulated the cloud of Fe 13+ ions by computing ion trajectories in the EBIT for different ion energy distributions used to initialize the trajectories. We then performed a least-squares fit to infer the ion energy distribution that best reproduced the measured ion cloud. These best-fit distributions were typically non-Maxwellian. For electron beam energies of 395–475 eV and electron beam currents of 1–9 mA, we find that the average ion energy is in the range of 10–300 eV. We also find that the average ion energy increases with increasing beam current approximately as $\langle$E$\rangle$ ≈ 25 I e eV, where I e is the electron beam current in mA. We have also compared our results to Maxwell-Boltzmann-distribution ion clouds. We find that our best-fit non-thermal distributions have an $\langle$E$\rangle$ that is less than half that of the T from the best-fit Maxwell-Boltzmann distributions ($\langle$E$\rangle$/q)/T=0.41 ± 0.05.

74 ATOMIC AND MOLECULAR PHYSICS↗

Spin Hall effect of radiofrequency waves in magnetized plasmas

In inhomogeneous media, electromagnetic-wave rays deviate from the trajectories predicted by the leading-order geometrical optics. This effect, called the spin Hall effect of light, is typically neglected in ray-tracing codes used for modeling waves in plasmas. Here, we demonstrate that the spin Hall effect can be significant for radiofrequency waves in toroidal magnetized plasmas whose parameters are in the ballpark of those used in fusion experiments. For example, an electron-cyclotron wave beam can deviate by as large as 10 wavelengths (~ 0.1 m) relative to the lowest-order ray trajectory in the poloidal direction. Finally, we calculate this displacement using gauge-invariant ray equations of extended geometrical optics, and we also compare our theoretical predictions with full-wave simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural units with time-dependent functionality

We show that the time-resolved dynamics of an underdamped harmonic oscillator can be used to do multifunctional computation, performing distinct computations at distinct times within a single dynamical trajectory. We consider the amplitude of an oscillator whose inputs influence its frequency. The activity of the oscillator at fixed times is a nonmonotonic function of its inputs, so it can solve problems such as XOR that are not linearly separable. The activity of the oscillator at fixed input is a nonmonotonic function of time, so it is multifunctional in a temporal sense, and able to carry out distinct nonlinear computations at distinct times within the same dynamical trajectory. We show that a single oscillator, observed at different times, can act as all of the elementary logic gates and perform binary addition, the latter usually implemented in hardware using five logic gates. We show that a set of n oscillators, observed at different times, can perform an arbitrary number of analog-to-n-bit digital conversions. We also show that oscillators can be trained by gradient descent to perform distinct classification tasks at distinct times. Computing with time-dependent functionality can be done in or out of equilibrium, and suggests a way of reducing the number of parameters or devices required to do nonlinear computations.

97 MATHEMATICS AND COMPUTING↗

Statistical properties of homogeneous and isotropic turbulence in He II measured via particle tracking velocimetry

Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hyperuniform Structures Formed by Shearing Colloidal Suspensions

In periodically sheared suspensions there is a dynamical phase transition, characterized by a critical strain amplitude $γ_c$, between an absorbing state where particle trajectories are reversible and an active state where trajectories are chaotic and diffusive. Repulsive nonhydrodynamic interactions between “colliding” particles’ surfaces have been proposed as a source of this broken time reversal symmetry. A simple toy model called random organization qualitatively reproduces the dynamical features of this transition. Random organization and other absorbing state models exhibit hyperuniformity, a strong suppression of density fluctuations on long length scales quantified by a structure factor $S(q → 0) ~q^α$ with $\textit{α}$ > 0, at criticality. In this work, we show experimentally that the particles in periodically sheared suspensions organize into structures with anisotropic short-range order but isotropic, long-range hyperuniform order when oscillatory shear amplitudes approach $γ_c$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU(1,1) Echoes for Breathers in Quantum Gases

Though the celebrated spin echoes have been widely used to reverse quantum dynamics, they are not applicable to systems whose constituents are beyond the control of the su(2) algebra. In this work, we design echoes to reverse quantum dynamics of breathers in three-dimensional unitary fermions and two-dimensional bosons and fermions with contact interactions, which are governed by an underlying su(1, 1) algebra. Geometrically, SU(1, 1) echoes produce closed trajectories on a single or multiple Poincaré disks and thus could recover any initial states without changing the sign of the Hamiltonian. In particular, the initial shape of a breather determines the superposition of trajectories on multiple Poincaré disks and whether the revival time has period multiplication. Our work provides physicists with a recipe to tailor collective excitations of interacting many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities

Fast ion phase-space flow, driven by Alfven eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D tokamak. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E – (ω/n)P ζ surfaces, where E, P ζ are energy and toroidal canonical momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. Here, the measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗