Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical systems”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Solving reaction dynamics with quantum computing algorithms

The description of quantum many-body dynamics is extremely challenging on classical computers, as it can involve many degrees of freedom. However, the time evolution of quantum states is a natural application for quantum computers that are designed to efficiently perform unitary transformations. Here, in this paper, we study quantum algorithms for response functions, relevant for describing different reactions governed by linear response. We focus on nuclear-physics applications and consider a qubit-efficient mapping on the lattice, which can efficiently represent the large volumes required for realistic scattering simulations. For the case of a contact interaction, we develop an algorithm for time evolution based on the Trotter approximation that scales logarithmically with the lattice size and is combined with quantum phase estimation. We eventually focus on the nuclear two-body system and a typical response function relevant for electron scattering as an example. We also investigate ground-state preparation and examine the total circuit depth required for a realistic calculation and the hardware noise level required to interpret the signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Magnetized laser–plasma interactions in high-energy-density systems: Parallel propagation

In this study, we investigate parametric processes in magnetized plasmas, driven by a large-amplitude pump light wave. Our focus is on laser–plasma interactions relevant to high-energy-density (HED) systems, such as the National Ignition Facility and the Sandia MagLIF concept. We present a self-contained derivation of a “parametric” dispersion relation for magnetized three-wave interactions, meaning the pump wave is included in the equilibrium, similar to the unmagnetized work of Drake et al., Phys. Fluids 17, 778 (1974). For this, we use a multi-species plasma fluid model and Maxwell's equations. The application of an external B field causes right- and left-polarized light waves to propagate with differing phase velocities. This leads to Faraday rotation of the polarization, which can be significant in HED conditions. Phase-matching and linear wave dispersion relations show that Raman and Brillouin scattering have modified spectra due to the background B field, though this effect is usually small in systems of current practical interest. We study a scattering process we call stimulated whistler scattering, where a light wave decays to an electromagnetic whistler wave (ω≲ω ce ) and a Langmuir wave. This only occurs in the presence of an external B field, which is required for the whistler wave to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

How Atomic Bonding Plays the Hardness Behavior in the Al–Co–Cr–Cu–Fe–Ni High Entropy Family

A systematic study on a face-centered cubic-based compositionally complex alloy system Al–Co–Cr–Cu–Fe–Ni in its single-phase state is carried out, where a mother senary compound Al 8 Co 17 Cr 17 Cu 8 Fe 17 Ni 33 and five of its suballoys, obtained by removing one element at a time, are investigated and exhaustively analyzed determining the contribution of each alloying element in the solid solution. The senary and the quinaries are compared using experimental techniques including X-ray absorption spectroscopy, X-ray diffraction, transmission electron microscopy, and first principles hybrid Monte Carlo/molecular dynamics simulations. Chemical short-range order and bond length distances have been determined both at the experimental and computational level. Electronic structure and local atomic distortions up to 5.2 Å have been correlated to the microhardness values. A linear regression model connecting hardness with local lattice distortions is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Non-adiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (wherethe electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase- space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Furthermore, the F-PSSH algorithm sets the stage for modeling nonadiabatic dynamics under strong externally pumped circular polarization in the future.

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

Resonant instabilities mediated by drag and electrostatic interactions in laboratory and astrophysical dusty plasmas

Dusty plasmas are known to support a diverse range of instabilities, including both generalizations of standard plasma instabilities and ones caused by effects specific to dusty systems. It has been recently demonstrated that a novel broad class of streaming instabilities, termed resonant drag instabilities (RDIs), can be attributed to a particular resonance phenomenon, manifested by defective eigenvalues of the linearized dust/fluid system. In this work, it is demonstrated that this resonance phenomenon is not unique to RDIs and can be used as a framework to understand a wider range of instabilities, termed resonant instabilities. Particular attention is given to the filamentary ionization instability seen in laboratory dusty plasmas and to the two-stream instability. It is shown that, due to the commonalities in underlying physics between the dust-ion-acoustic two-stream instability and the acoustic RDI, these instabilities should be relevant in strongly overlapping regimes in astrophysical dusty plasmas. Further, it is proposed that a similar overlap in the experimental accessibility of these modes (and of the filamentary instability) allows for the possibility of experimental investigation in the laboratory of complex and astrophysically relevant instability dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Diffusion-mediated passing of molecular species in linear nanopores constrained by orientational alignment

For diffusion-mediated catalytic conversion reactions in materials with narrow linear nanopores, e.g., mesoporous silica MCM-41, a key parameter is the propensity for product species to be able to pass reactant species and thus to efficiently exit the pore. For elongated species, this can require orientational alignment with the pore axis. In this work, we perform benchmark analyses for such solution-phase systems where one of these species is elongated in order to quantify the dependence of this passing propensity, P, on pore diameter and on the rotational diffusion coefficient, Dr, of the elongated species. In particular, we consider the passing of a spherical and an elongated spherocylindrical shaped species in a cylindrical pore in an implicit solvent, where these species cannot overlap. Passing is mediated by diffusive Brownian motion of these species as described by strongly damped Langevin dynamics. We quantify scaling of P for pore width just above the threshold where passing is sterically blocked, and also reveal a significant decrease in P for lower Dr. We also consider the dependence of P on the aspect ratio of the elongated species and obtain an exact result in the limiting regime of large aspect ratio.

Brownian motion↗
DOE OSTI2024

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗
DOE OSTI2021

Quasi-one-dimensional uniform spin- 1 2 Heisenberg antiferromagnet KNaCuP 2 O 7 probed by P 31 and Na 23 NMR

Here we present the structural and magnetic properties of KNaCuP 2 O 7 investigated via x-ray diffraction, magnetization, specific heat, and P 31 and Na 23 NMR measurements and complementary electronic structure calculations. The temperature-dependent magnetic susceptibility and P 31 NMR shift could be modeled very well by the uniform spin- 1 2 Heisenberg antiferromagnetic chain model with a nearest-neighbor interaction J / k B ≃ 58.7 K . The corresponding mapping using first-principles electronic structure calculations leads to J DFT / k B ≃ 59 K with negligibly small interchain couplings, further confirming that the system is indeed a one-dimensional uniform spin- 1 2 Heisenberg antiferromagnet. The diverging trend of NMR spin-lattice relaxation rates ( 1 / T 1 31 and 1 / T 1 23 ) implies the onset of a magnetic long-range ordering at around T N ≃ 1 K . From the value of T N , the average interchain coupling is estimated to be J ' / k B ≃ 0.28 K . Moreover, the NMR spin-lattice relaxation rates show the dominant contributions from uniform ( q = 0 ) and staggered ( q = ± π / a ) spin fluctuations in the high- and low-temperature regimes, respectively, mimicking one-dimensionality of the spin lattice. We have also demonstrated that 1 / T 1 31 in high temperatures varies linearly with 1 / H , reflecting the effect of spin diffusion on the dynamic susceptibility. The temperature-dependent unit cell volume could be described well using the Debye approximation with a Debye temperature of Θ D ≃ 294 K , consistent with the heat capacity data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Concept of circular-linear energy recovery accelerator to probe the energy frontier

Energy-frontier accelerators provide powerful tools performing high precision measurements confirming the fundamentals of the physics and broadening new research horizons. Such machines are either driven by circular or linear accelerators. The circular machines, having the centre-of-mass (CM) energy values reaching 200 GeV (for leptons) and above, experience beam energy loss and quality dilution, for example, due to synchrotron radiation, limiting the overall CM energy achievable and requiring a constant energy top-up to compensate the loss and the beam quality dilution. Linear colliders overcome these limitations, while the finite capabilities of generating high average current beams limits the luminosity. This is partially compensated by the quality of the colliding beams. In this work, we suggest a novel design of circular-linear accelerator based on the merging of the “non-emitting”, low-energy storage rings and energy recovery linear accelerators. We suggest using the recently considered dual-axis asymmetric cavities to enable the operation of such a system, and in particular the energy recovery from spent, high-intensity beams. The machine considered, under the scope of the SNOWMASS-2021 initiative, can be potentially used to reach ultimate energy frontiers in high-energy physics as well as to drive next generation light sources. The merging of circular and linear systems, and applications of dual axes cavities, should allow the maintaining of high beam quality, high luminosity, and high energy efficiency, while offering a flexible energy management and opening clear opportunity for reducing the running cost. We note that the numbers shown in the paper are for illustration purpose and can be improved further.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Magnetic field induced phases and spin Hamiltonian in Cs 2 CoBr 4

Magnetic structures and spin excitations are studied across the phase diagram of the geometrically frustrated S = 3/2 quantum antiferromagnet Cs 2 CoBr 4 in magnetic fields applied along the magnetic easy axis, using neutron diffraction, inelastic neutron scattering, and terahertz absorption spectroscopy. Here, the data are analyzed, where appropriate, using extended SU(4) linear spin wave theory. A minimal magnetic Hamiltonian is proposed based on measurements in the high-field paramagnetic state. It deviates considerably from the previously considered models. Additional dilatometry experiments highlight the importance of magnetoelastic coupling in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Physics-informed regularization and structure preservation for learning stable reduced models from data with operator inference

Operator inference learns low-dimensional dynamical-system models with polynomial nonlinear terms from trajectories of high-dimensional physical systems (non-intrusive model reduction). Here, this work focuses on the large class of physical systems that can be well described by models with quadratic and cubic nonlinear terms and proposes a regularizer for operator inference that induces a stability bias onto learned models. The proposed regularizer is physics informed in the sense that it penalizes higher-order terms with large norms and so explicitly leverages the polynomial model form that is given by the underlying physics. This means that the proposed approach judiciously learns from data and physical insights combined, rather than from either data or physics alone. Additionally, a formulation of operator inference is proposed that enforces model constraints for preserving structure such as symmetry and definiteness in linear terms. Numerical results demonstrate that models learned with operator inference and the proposed regularizer and structure preservation are accurate and stable even in cases where using no regularization and Tikhonov regularization leads to models that are unstable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Plasma Wave and Particle Dynamics During Interchange Events in the Jovian Magnetosphere Using Juno Observations

Abstract Interchange instability is known to drive fast radial transport of particles in Jupiter's inner magnetosphere. Magnetic flux tubes associated with the interchange instability often coincide with changes in particle distributions and plasma waves, but further investigations are required to understand their detailed characteristics. We analyze representative interchange events observed by Juno, which exhibit intriguing features of particle distributions and plasma waves, including Z‐mode and whistler‐mode waves. These events occurred at an equatorial radial distance of ∼9 Jovian radii on the nightside, with Z‐mode waves observed at mid‐latitude and whistler‐mode waves near the equator. We calculate the linear growth rate of whistler‐mode and Z‐mode waves based on the observed plasma parameters and electron distributions and find that both waves can be locally generated within the interchanged flux tube. Our findings are important for understanding particle transport and generation of plasma waves in the magnetospheres of Jupiter and other planetary systems.

Geology↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fine-root dynamics vary with soil depth and precipitation in a lown-utrient tropical forest in the Central Amazonia

A common assumption in tropical ecology is that root systems respond rapidly to climatic cues but that most of that response is limited to the uppermost layer of the soil, with relatively limited changes in deeper layers. However, this assumption has not been tested directly, preventing models from accurately predicting the response of tropical forests to environmental change. We measured seasonal dynamics of fine roots in a Central Amazonia mature forest using minirhizotrons to 90 cm depth, which were calibrated with fine roots extracted from soil cores. Root productivity and mortality in superficial soil layers were positively correlated with precipitation, whereas root standing length was greater during the dry periods at the deeper layers. Contrary to historical assumptions, a large fraction of fine-root standing biomass (46%) and productivity (41%) were found in soil layers deeper than 30 cm. Furthermore, root turnover decreased linearly with soil depth. Our findings demonstrate a strong relationship between fine-root dynamics and precipitation regimes in Central Amazonia. Our results also emphasize the importance of deep roots for accurate estimates of primary productivity and the interaction between roots and carbon, water, and nutrients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Theoretical Investigation of the Oxygen Reduction Reaction over Platinum Catalysts Supported by Multi-Edged Vertically Aligned Carbon Nanofiber for Electrocatalyst Preparation

Vertically aligned carbon nanofibers (VACNFs) are promising supports for oxygen reduction reaction (ORR) electrocatalysts in fuel cells. Although experimentally these catalytic systems have shown great potential, there is lack of molecular understanding of the catalytic sites and reaction mechanisms. This work investigated the origin of the ORR reactivities of the platinum catalysts on multi-edged VACNFs (Pt/VACNF) using a multiscale modeling approach combining Density Functional Theory (DFT) and classical Molecular Dynamics (MD) simulations. Based on the ReaxFF potential, here all nanoscale Pt particles (Pt 55 , P 100 , and Pt 147 ) are stabilized by the open edges located axially along the VACNF walls. The calculated first-shell coordination numbers, CN, of surface Pt atoms are 6.63, 7.27, and 7.85, respectively, suggesting that the percentage of low coordination sites increases as the particle size decreases. The adsorption energies of OOH, O, and OH on Pt-55 were systematically probed using DFT calculations. These adsorption energies retain a linear correlation against the generalized coordination numbers ($\overline{CN}$). For Pt nanoparticles supported on VACNF, we found that the OOH and OH bind stronger than on Pt (111) by 0.14 and 0.17 eV, respectively, which can hinder the ORR activity with lower limiting potential than Pt (111). Our theoretical prediction is in good agreement with the linear sweeping voltammetry that revealed a left shift of the half-wave potential.

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

Adaptive Variational Quantum Computing Approaches for Green’s Functions and Nonlinear Susceptibilities

Here, we present and benchmark quantum computing approaches for calculating real-time single-particle Green’s functions and nonlinear susceptibilities of Hamiltonian systems. The approaches leverage adaptive variational quantum algorithms for state preparation and propagation. Using automatically generated compact circuits, the dynamical evolution is performed over sufficiently long times to achieve adequate frequency resolution of the response functions. We showcase accurate Green’s function calculations using a statevector simulator on classical hardware for Fermi-Hubbard chains of 4 and 6 sites, with maximal ansatz circuit depths of 65 and 424 layers, respectively, and for the molecule LiH with a maximal ansatz circuit depth of 81 layers. Additionally, we consider an antiferromagnetic quantum spin-1 model that incorporates the Dzyaloshinskii-Moriya interaction to illustrate calculations of the third-order nonlinear susceptibilities, which can be measured in two-dimensional coherent spectroscopy experiments. These results demonstrate that real-time approaches using adaptive parametrized circuits to evaluate linear and nonlinear response functions can be feasible with near-term quantum processors.

97 MATHEMATICS AND COMPUTING↗