Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EQUATIONS OF MOTION”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interacting mesons as degrees of freedom in a chiral model

Here, we study the equation of state of hot and dense hadronic matter using an extended chiral mean field (CMF) model framework where the addition is the inclusion of interactions of thermally excited mesons. This is implemented by calculating the in-medium masses of pseudoscalar and vector mesons, obtained through the explicit chiral symmetry-breaking and vector-interaction terms in the Lagrangian, respectively, prior to applying the mean-field approximation. As a result, the in-medium meson contributions generate a feedback term to the CMF’s equations of motion, which then modifies the equation of state. With this improvement, we quantify the effect on the equation of state of strongly interacting matter through comparisons with state-of-the-art lattice QCD results and other hadronic models like the hadron resonance gas model. We find that the results of the updated hadronic CMF model with an improved meson description (mCMF) provide a better agreement with lattice-QCD data for thermodynamic state variables across a wide range of temperatures and baryon chemical potentials.

phenomenology↗
DOE OSTI2023

Code Demonstration: Approximation of nearly-periodic symplectic maps via structure-preserving neural networks

In this notebook, we use the symplectic gyroceptron architecture from [Duruisseaux et al., 2022] to learn a surrogate map for the nearly-periodic symplectic flow map associated to a nearly-periodic Hamiltonian system composed of two nonlinearly coupled oscillators, where one of them oscillates significantly faster than the other: $\Big\lbrace$ $ {\dot{q}_1 = p_1 \hspace{5mm} \dot{p}_1 = -q_1 - \epsilon \partial_{q1}U(q_1 , q_2 )} \atop {\dot{q}_2 = \epsilon p_2 \hspace{4mm} \dot{p}_2 = -\epsilon q_2 - \epsilon \partial_{q2}U(q_1 , q_2 )}$. These equations of motion are the Hamilton's equations associated to the Hamiltonian $H_{\epsilon}(q_1,q_2,p_1,p_2)=\frac{1}{2}(q^2_1+p^2_1)+\frac{1}{2}\epsilon(q^2_2+p^2_2)+\epsilon$$U(q_1,q_2)$.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Condensation of collective polar vortex modes

The dynamics of extended objects such as domain walls, domain bubbles, vortex structures, etc., can be described by their equations of motion associated with their effective mass and spring constant. Here we analytically derive the equations of motion for the polarization dynamics and elastodynamics for the structural responses of ferroelectric polar vortices, and theoretically extract their effective mass, spring constant, and mode frequencies. We demonstrate two subterahertz phonon modes and predicted their frequencies, both consistent with our recent experimental measurements and phase-field simulations. We show that elastic modulation of the energy function and spring constants leads to a condensation of a collective mode upon a second-order structural transition from symmetric to asymmetric vortices at a critical strain, analogous to the ferroelectric soft phonon mode at a ferroelectric transition. Furthermore, the present work offers a theoretical framework for predicting and manipulating the ultrafast collective dynamics of polar nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Supergravity spectrum of AdS 5 black holes

We embed Kerr-Newman-AdS black holes into N = 8 gauged supergravity and study quadratic fluctuations around the black hole backgrounds of all fields in the larger theory. The equations of motion of the perturbations are partially diagonalized by the group theory of broken symmetry. Nearly all fields in theory have non-minimal couplings, so their equations of motion are not merely massive Klein-Gordon equations with minimal coupling to background gauge fields, and their analogues for fields with spin. In the special case of extremal black holes we identify specific modes of instability, some of which touch supersymmetric locus. For example, we identify scalar fields in supergravity that condense in the near horizon region and transition the black hole into a superconducting phase. We also identify supergravity modes that are susceptible to superradiant instability.

AdS-CFT Correspondence↗
DOE OSTI2020

Horizon radiation reaction forces

Using Effective Field Theory (EFT) methods, we compute the effects of horizon dissipation on the gravitational interactions of relativistic binary black hole systems. We assume that the dynamics is perturbative, i.e it admits an expansion in powers of Newton’s constant (post-Minkowskian, or PM, approximation). As applications, we compute corrections to the scattering angle in a black hole collision due to dissipative effects to leading PM order, as well as the post-Newtonian (PN) corrections to the equations of motion of binary black holes in non-relativistic orbits, which represents the leading order finite size effect in the equations of motion. The methods developed here are also applicable to the case of more general compact objects, eg. neutron stars, where the magnitude of the dissipative effects depends on non-gravitational physics (e.g, the equation of state for nuclear matter).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

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

Conformal wavefunctions for graviton amplitudes

The extended-BMS algebra of asymptotically flat spacetime contains an SO(3, 1) subgroup that acts by conformal transformations on the celestial sphere. It is of interest to study the representations of this subgroup associated with gravitons. To reduce the equation of motion to a Schrodinger-like equation it is necessary to impose a non-covariant gauge condition. Using these solutions, leading-order gauge invariant Weyl scalars are then computed and decomposed into families of unitary principal series representations. An invertible holographic mapping is constructed between these unitary principal series operators and massless spin-2 perturbations of flat spacetime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Predicting Open Quantum Dynamics with Data-Informed Quantum-Classical Dynamics

We introduce a data-informed quantum-classical dynamics (DIQCD) approach for predicting the evolution of an open quantum system. The equation of motion in DIQCD is a Lindblad equation with a flexible, time-dependent Hamiltonian that can be optimized to fit sparse and noisy data from local observations of an extensive open quantum system. We demonstrate the accuracy and efficiency of DIQCD for both experimental and simulated quantum devices. We show that DIQCD can predict entanglement dynamics of ultracold molecules (calcium fluoride) in optical tweezer arrays. DIQCD also successfully predicts carrier mobility in organic semiconductors (rubrene) with accuracy comparable to nearly exact numerical methods.

Lindblad equation↗
DOE OSTI2025

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Perfect spinfluid: A divergence-type approach

We present a new formulation of nondissipative relativistic spin hydrodynamics that incorporates spin degrees of freedom into the divergence-type theory framework. Due to the divergence-type structure, it is straightforward to enforce nonlinear causality and symmetric hyperbolicity of the equations of motion, ensuring local well-posedness of the initial-value problem and stability of the theory. Furthermore, in a specific realization based on spin kinetic theory, we prove that the equations of motion remain nonlinearly causal and symmetric-hyperbolic to all orders in the spin potential, provided a specific thermodynamic constraint is satisfied. Here, this framework can be applied for numerical simulations to study the dynamics of spin-polarized fluids, such as the quark-gluon plasma in heavy-ion collisions.

Chirality↗
DOE OSTI2022

Charged particle motion in spherically symmetric distributions of magnetic monopoles

The classical equations of motion of a charged particle in a spherically symmetric distribution of magnetic monopoles can be transformed into a system of linear equations, thereby providing a type of integrability. In the case of a single monopole, the solution was given long ago by Poincaré. In the case of a uniform distribution of monopoles, the solution can be expressed in terms of parabolic cylinder functions (essentially the eigenfunctions of an inverted harmonic oscillator). Further, this solution is relevant to recent studies of nonassociative star products, symplectic lifts of twisted Poisson structures, and fluids and plasmas of electric and magnetic charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Electron-induced massive dynamics of magnetic domain walls

We study the dynamics of domain walls (DWs) in a metallic, ferromagnetic nanowire, focusing on inertial effects on the DW due to interaction with a conduction electron bath. In this work, we develop a Keldysh collective coordinate technique to describe the effect of conduction electrons on rigid magnetic structures. The effective Lagrangian and Langevin equations of motion for a DW are derived microscopically, including the full response kernel which is nonlocal in time. The DW dynamics is described by two collective degrees of freedom: position and tilt angle. The coupled Langevin equations therefore involve two correlated noise sources, leading to a generalized fluctuation-dissipation theorem (FDT). The DW response kernel due to electrons contains two parts: one related to dissipation via FDT and another reactive part. We prove that the latter term leads to a mass for both degrees of freedom, even though the intrinsic bare mass is zero. The electron-induced mass is present even in a clean system without pinning or specifically engineered potentials. The resulting equations of motion contain rich dynamical solutions and point toward a way to control domain wall motion in metals via the electronic system properties. We discuss two observable consequences of the mass, hysteresis in the DW dynamics, and resonant response to ac current.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗
DOE OSTI2022

Online Machine Learning for Accelerating Molecular Dynamics Modeling of Cells

We developed a biomechanics-informed online learning framework to learn the dynamics with ground truth generated with multiscale modeling simulation. It was built on Summit-like supercomputers, which were also used to benchmark and validate our framework on one physiologically significant modeling of deformable biological cells. We generalized the century-old equation of Jeffery orbits to a new equation of motion with additional parameters to account for the flow conditions and the cell deformability. Using simulation data at particle-based resolutions for flowing cells and the learned parameters from our framework, we validated the new equation by the motions, mostly rotations, of a human platelet in shear blood flow at various shear stresses and platelet deformability. Our online framework, which surrogates redundant computations in the conventional multiscale modeling by solutions of our learned equation, accelerates the conventional modeling by three orders of magnitude without visible loss of accuracy.

multiscale modeling↗
DOE OSTI2021

Symbolic pregression: Discovering physical laws from distorted video

In this work, we present a method for unsupervised learning of equations of motion for objects in raw and optionally distorted unlabeled synthetic video (or, more generally, for discovering and modeling predictable features in time-series data). We first train an autoencoder that maps each video frame into a low-dimensional latent space where the laws of motion are as simple as possible, by minimizing a combination of nonlinearity, acceleration, and prediction error. Differential equations describing the motion are then discovered using Pareto-optimal symbolic regression. We find that our pre-regression (“pregression”) step is able to rediscover Cartesian coordinates of unlabeled moving objects even when the video is distorted by a generalized lens. Using intuition from multidimensional knot theory, we find that the pregression step is facilitated by first adding extra latent space dimensions to avoid topological problems during training and then removing these extra dimensions via principal component analysis. An inertial frame is autodiscovered by minimizing the combined equation complexity for multiple experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗