Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Computational Modeling of Vertical-Cavity Surface-Emitting Lasers

Vertical-Cavity Surface-Emitting Lasers (VCSELs) have become increasingly important for many applications, such as optical interconnects and optical data storages. In many applications, critical requirements are very often high beam quality and/or high output powers. To help the design and optimization of VCSELs for such applications, it is essential to have a simulation model that includes the time evolution of the transverse space dependence of the hewn. Such models allow arbitrary transverse profiles of laser intensity or carrier density to develop in time without any a priori assumption about the number and types of transverse modes. In addition, extensive research in the past has shown that many-body effects are important in influencing the optical gain and the refractive index of semiconductor quantum wells. It is thus important to include the key many-body effects in simulation models for semiconductor lasers. A computational model, called the Maxwell-Effective Bloch Equations (MEBE) was developed recently that simulates the above features and effects. This model combines the space-time resolution with microscopic many-body effects in a way that is both accurate enough and computationally manageable. This model has been used to simulate edge-emitting high power lasers recently. A new finite-difference algorithm has been developed for solving these equations and will be shown in the presentation. That algorithm will be used to calculate the detailed time evolution of the transverse mode profiles of VCSELs. The VCSELs are based on InGaAs/GaAs quantum well structures. Effects of pumping current and refractive index profiles will be investigated in detail. The simulation results will be compared with other results the use the rate-equation model.

Goorjian, Peter M.↗
NASA NTRS1998

Computational Nonlinear of Vertical-Cavity Surface-Emitting Lasers

Vertical-Cavity Surface-Emitting Lasers (VCSELs) have become increasingly important for many applications, such as optical interconnects and optical data storages. In many applications, critical requirements are very often high beam quality and/or high output powers. To help the design and optimization of VCSELs for such applications, it is essential to have a simulation model that includes the time evolution of the transverse space dependence of the beam. Such models allow arbitrary transverse profiles of laser intensity or carrier density to develop in time without any a priori assumption about the number and types of transverse modes. In addition, extensive research in the past has shown that many-body effects are important in influencing the optical gain and the refractive index of semiconductor quantum wells. It is thus important to include the key many-body effects in simulation models for semiconductor lasers. A computational model, called the Maxwell-Effective Bloch Equations (MEBE) was developed recently that simulates the above features and effects. This model combines the space-time resolution with microscopic many-body effects in a way that is both accurate enough and computationally manageable. This model has been used to simulate edge-emitting high power lasers recently. A new finite-difference algorithm has been developed for solving these equations and will be shown in the presentation. That algorithm will be used to calculate the detailed time evolution of the transverse mode profiles of VCSELs. The VCSELs are based on InGaAs/GaAs quantum well structures. Effects of pumping current and refractive index profiles will be investigated in detail. The simulation results will be compared with other results the use the rate-equation model.

Goorjian, Peter M.↗
NASA NTRS1995

GRO: Black hole models for gamma-ray bursts

The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRB's) is isotropic but is bound radially. This finding suggests that the bursts are either cosmological or they originate from an extended Galactic halo. The implied luminosities and the observed variability of the GRB's on time scales as short as one millisecond suggest that they originate from compact objects. We are presently studying black hole models for GRB's. Any such model must produce a non-thermal photon spectrum to agree with the observed properties. For a wide range of burst parameters the assumed bursting source consists of a non-thermal electron-positron-photon plasma of very high density. It seems possible to produce such a plasma in accretion onto black holes. In our on-going work, we are developing the kinetic theory for a non-equilibrium pair plasma. The main new features of our work are as follows: (1) We do not assume the presence of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an arbitrary concentration and energy distribution. (3) There is no soft photon source in our model; initially all the photons in the plasma are either energetic X-rays or gamma-rays. (4) The initial energy distribution of the pairs as well as photons is arbitrary. (5) We collect the analytical expressions for the kinetic kernels for all relevant processes. And (6) we present a different approach to finding the time-evolution of pair and photon spectra, which is a combination of the kinetic-theory and the non-linear Monte-Carlo schemes. We have developed many Monte-Carlo programs to model various process, to take into account the time evolution, and to incorporate various physical effects which are unique to non-thermal plasmas. The hydrodynamics of fireballs in GRB's was studied before. Applying results from kinetic theory will improve our understanding of these systems.

Ruderman, Malvin↗
DOE OSTI2024

Nearly optimal state preparation for quantum simulations of lattice gauge theories

Here, we present several improvements to the recently developed ground-state preparation algorithm based on the quantum eigenvalue transformation for unitary matrices (QETU), apply this algorithm to a lattice formulation of U(1) gauge theory in (2+1) dimensions, as well as propose an alternative application of QETU, a highly efficient preparation of Gaussian distributions. The QETU technique was originally proposed as an algorithm for nearly optimal ground-state preparation and ground-state energy estimation on early fault-tolerant devices. It uses the time-evolution input model, which can potentially overcome the large overall prefactor in the asymptotic gate cost arising in similar algorithms based on the Hamiltonian input model. We present modifications to the original QETU algorithm that significantly reduce the cost for the cases of both exact and Trotterized implementation of the time evolution circuit. We use QETU to prepare the ground state of a U(1) lattice gauge theory in two spatial dimensions, explore the dependence of computational resources on the desired precision and system parameters, and discuss the applicability of our results to general lattice gauge theories. We also demonstrate how the QETU technique can be utilized for preparing Gaussian distributions and wave packets in a way which outperforms existing algorithms for as little as n q ≳ 2–5 qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

A numerically exact full wave packet approach to molecule-surface scattering

A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.

Mowrey, R. C.↗
DOE OSTI2026

Recurrent convolutional neural networks for modeling nonadiabatic dynamics of quantum-classical systems

Recurrent neural networks (RNNs) have recently been extensively applied to model the time evolution in fluid dynamics, weather predictions, and even chaotic systems due to their ability to capture temporal dependencies and sequential patterns in data. Here we present an RNN model based on convolutional neural networks for modeling the nonlinear nonadiabatic dynamics of hybrid quantum-classical systems. The dynamical evolution of the hybrid systems is governed by equations of motion for classical degrees of freedom and von Neumann equation for electrons. The Physics-Aware Recurrent Convolution (PARC) neural network structure incorporates a differentiator-integrator architecture that inductively models the spatiotemporal dynamics of generic physical systems. Here, we apply our RNN approach to learn the space-time evolution of a one-dimensional semiclassical Holstein model after an interaction quench. For shallow quenches (small changes in electron-lattice coupling), the deterministic dynamics can be accurately captured using a single-CNN-based recurrent network. In contrast, deep quenches induce chaotic evolution, making long-term trajectory prediction significantly more challenging. Nonetheless, we demonstrate that the PARC-CNN architecture can effectively learn the statistical climate of the Holstein model under deep-quench conditions.

Holstein model↗
DOE OSTI2024

Quantum simulations of hadron dynamics in the Schwinger model using 112 qubits

Hadron wave packets are prepared and time evolved in the Schwinger model using 112 qubits of IBM’s 133-qubit Heron quantum computer ibm_torino. The initialization of the hadron wave packet is performed in two steps. First, the vacuum is prepared across the whole lattice using the recently developed SC-ADAPT-VQE algorithm and workflow. SC-ADAPT-VQE is then extended to the preparation of localized states, and used to establish a hadron wave packet on top of the vacuum. This is done by adaptively constructing low-depth circuits that maximize the overlap with an adiabatically prepared hadron wave packet. Due to the localized nature of the wavepacket, these circuits can be determined on a sequence of small lattices using classical computers, and then robustly scaled to prepare wave packets on large lattices for simulations using quantum computers. Time evolution is implemented with a second-order Trotterization. To reduce both the required qubit connectivity and circuit depth, an approximate quasilocal interaction is introduced. This approximation is made possible by the emergence of confinement at long distances, and converges exponentially with increasing distance of the interactions. Using multiple error-mitigation strategies, up to 14 Trotter steps of time evolution are performed, employing 13,858 two-qubit gates (with a CNOT depth of 370). The propagation of hadrons is clearly identified, with results that compare favorably with Matrix Product State simulations. Finally, prospects for a near-term quantum advantage in simulations of hadron scattering are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2022

Impact of Meteorological Factors on the Mesoscale Morphology of Cloud Streets during a Cold-Air Outbreak over the Western North Atlantic

Postfrontal clouds (PFC) are ubiquitous in the marine boundary layer, and their morphology is essential to estimating the radiation budget in weather and climate models. Here we examine the roles of sea surface temperature (SST) and meteorological factors in controlling the mesoscale morphology and evolution of shallow clouds associated with a cold-air outbreak that occurred on 1 March 2020 during phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). Our results show that the simulated PFC structure and ambient conditions by the Weather Research and Forecasting (WRF) Model are generally consistent with observations from GOES-16 and dropsonde measurements. We also examine the thermodynamical and dynamical influences in the cloud mesoscale morphology using WRF sensitivity experiments driven by two meteorological forcing datasets with different domain-mean SST and spatial gradients, which lead to dissimilar values of hydrometeor water path and cloud core fraction. The SST from ERA5 leads to weaker stability and higher inversion height than the SST from FNL does. In addition, the use of large-scale meteorological forcings from ERA5 yields a distinctive time evolution of wind direction shear in the inner domain, which favors the formation and persistence of longer cloud rolls. Both factors contribute to a change in the time evolution of domain-mean water path and cloud core fraction of cloud streets. Our study takes advantage of the simulation driven by the differences between two large-scale forcing datasets to illustrate the importance of SST and wind direction shear in the cloud street morphology in a realistic scenario

Jingyi Chen↗
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↗
NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗
NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibrium states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, C. G.↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Semicoherent symmetric quantum processes: Theory and applications

Discovering pragmatic and efficient approaches to construct ε-approximations of quantum operators such as real (imaginary) time-evolution propagators in terms of the basic quantum operations (gates) is challenging. Prior ε-approximations are invaluable, in that they enable the compilation of classical and quantum algorithm modeling of, e.g., dynamical and thermodynamic quantum properties. In parallel, symmetries are powerful tools concisely describing the fundamental laws of nature; the symmetric underpinnings of physical laws have consistently provided profound insights and substantially increased predictive power. In this work, we consider the interplay between the ε-approximate processes and the exact symmetries in a semicoherent context—where measurements occur at each logical clock cycle. Here we draw inspiration from Pascual Jordan's groundbreaking formulation of nonassociative, but commutative, symmetric algebraic form. Our symmetrized formalism is then applied in various domains such as quantum random walks, real-time evolutions, variational algorithm ansatzes, and efficient entanglement verification. Our work paves the way for a deeper understanding and greater appreciation of how symmetries can be used to control quantum dynamics in settings where coherence is a limited resource.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2020

The Late-time Afterglow Evolution of Long Gamma-Ray Bursts GRB 160625B and GRB 160509A

We present post-jet-break Hubble Space Telescope, Very Large Array, and Chandra observations of the afterglow of the long γ-ray bursts GRB 160625B (between 69 and 209 days) and GRB 160509A (between 35 and 80 days). We calculate the post-jet-break decline rates of the light curves and find the afterglow of GRB 160625B is inconsistent with a simple t−3/4 steepening over the break, expected from the geometric effect of the jet edge entering our line of sight. However, the favored optical post-break decline (fv alpha t-1.96+-0.07) is also inconsistent with the fν alpha t−p decline (where p ≈ 2.3 from the pre-break light curve), which is expected from exponential lateral expansion of the jet; perhaps suggesting lateral expansion that only affects a fraction of the jet. The post break decline of GRB 160509A is consistent with both the t−3/4 steepening and with fν alpha t−p. We also use BOXFIT to fit afterglow models to both light curves and find both to be energetically consistent with a millisecond magnetar central engine, but the magnetar parameters need to be extreme (i.e., E ~ 3 × 1052 erg). Finally, the late-time radio light curves of both afterglows are not reproduced well by BOXFIT and are inconsistent with predictions from the standard jet model; instead, both are well represented by a single power-law decline (roughly fν alpha t−1) with no breaks. This requires a highly chromatic jet break (tj, radio > 10 x tj, optical) and possibly a two-component jet for both bursts.

Gamma-ray bursts↗
DOE OSTI2026

Trigonometric continuous-variable gates and hybrid quantum simulations of the sine-Gordon model

Hybrid qubit-qumode quantum computing platforms provide a natural setting for simulating interacting bosonic quantum field theories. However, existing continuous-variable gate constructions rely predominantly on polynomial functions of canonical quadratures. In this work, we introduce a complementary universality paradigm based on trigonometric continuous-variable gates, which enable a Fourier-like representation of bosonic operators and are particularly well suited for periodic and non-perturbative interactions. We present an ancilla-based framework for implementing trigonometric gates with arguments given by arbitrary Hermitian functions of qumode quadratures. The protocol yields unitary gates deterministically, and non-unitary gates through probabilistic post-selection. As a concrete application, we develop a hybrid qubit-qumode quantum simulation of the lattice sine-Gordon model. Using these gates, we prepare ground states via quantum imaginary-time evolution, simulate real-time dynamics, compute time-dependent vertex two-point correlation functions, and extract quantum kink profiles under topological boundary conditions. Our results demonstrate that trigonometric continuous-variable gates provide a physically natural framework for simulating interacting field theories on near-term hybrid quantum hardware, while establishing a parallel route to universality beyond polynomial gate constructions. We expect that the trigonometric gates introduced here to find broader applications, including quantum simulations of condensed matter systems, quantum chemistry, and biological models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Analysis of anion exchange membrane water electrolyzer performance and its evolution over time

Understanding water, evolved gas, and ionic transport in membrane-electrode-assemblies (MEAs) is essential for the development of high performance and durable anion exchange membrane water electrolyzers (AEMWEs). This study evaluates the MEA conditioning process, operating conditions, and short-term stability in a 1 M potassium hydroxide (KOH) electrolyte, focusing on the underlying transport phenomena. We observe a significant initial voltage loss in continuous cell operation, which could be associated with gas bubble accumulation, transport layer or flow field passivation, and changes in the catalyst oxidation state. Further, we investigate the effects of materials and operational configurations, including the membrane type and thickness, and the electrolyte flow rate, including KOH being fed to both electrodes as well as to the anode only. Furthermore, the effect of membrane drying temperature on ex situ as well as in situ electrochemical performance is evaluated. Finally, we discuss 700 h of AEMWE operation at 1 A/cm 2 , highlighting the underlying degradation phenomena.

25 ENERGY STORAGE↗
NASA NTRS1980

Grid evolution in time asymptotic problems

A technique for generating systems of coordinates for solving time asymptotic problems is described which provides a simple way of moving the mesh points in physical space and reduces the error in the solution relative to that obtained using a fixed mesh. First order partial differential equations are formulated for the grid point velocity in transient problems. Local flow information and boundary motion are used to determine the interior grid point motion.

Rai, M. M.↗
NASA NTRS1981

Grid evolution in time asymptotic problems

The selection of the proper coordinate system in solving any fluid flow or heat transfer problem is a very important consideration. A new technique of moving mesh points in physical space is introduced so as to reduce the error in a computed asymptotic solution relative to that obtained using a fixed mesh. The technique has been used to solve the simple viscous Burgers' equation in one and two dimensions. Substantial error reductions are demonstrated. The treatment of boundary points and the effect of using different error criteria in generating grids are discussed.

Rai, M. M.↗