Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 271 records · Page 15

Gate-Based Quantum Simulation of Gaussian Bosonic Circuits on Exponentially Many Modes

We introduce a framework for simulating, on an ( n + 1 )-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2 n modes. Specifically, we encode the initial bosonic state’s expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼ 8 × 10 9 modes, illustrating the power of our framework. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical calculations of the collapse of nonrotating, magnetic gas clouds

Results of the first self-consistent numerical calculations of the dynamic collapse of a magnetized protostellar gas cloud are presented. Symmetry about an axis parallel to the initial magnetic field direction has been assumed, so that the calculations could be performed on a two-dimensional grid. Also, the cloud was taken to be nonrotating and isothermal, and the magnetic field was assumed to remain frozen in to the gas. As starting models for the calculations, gas spheres with uniform density and magnetic field were used. The time evolution of the clouds has been calculated for roughly two initial free-fall times, at which point the central density has increased by a factor of approximately 10,000 to 1,000,000. Several such calculations have been performed for different values of the cloud's initial thermal, magnetic, and gravitational energies. In virtually all cases it is found that, once a flattened core forms in the cloud, the central magnetic field strength, B, varies with gas density, rho, according to (d log B/d log rho) = 1/2. This behavior is independent of the initial energy ratios mentioned above. It is also found that the magnetic field is able to prevent completely the collapse of part of the outer envelope of the cloud.

Scott, E. H.↗

Multiparameter radar study of rainfall: Potential application to area-time integral studies

Multiparameter radars measure one or more additional parameters in addition to the coventional reflectivity factor. The combination of radar observations from a multiparameter radar is used to study the time evolution of rainstorms. A technique is presented to self-consistently compare the area-time integral (ATI) and rainfall volume estimates from convective storms, using two different measurements from a multiparameter radar. Rainfall volumes for the lifetime of individual storms are computed using the reflectivity at S band (10-cm wavelength) as well as one-way specific attenuation at X band (3-cm wavelength). Area-time integrals are computed by summing all areas in each radar snapshot having reflectivities (S band) in excess of a preselected threshold. The multiparameter radar data used in this study were acquired by the National Center for Atmospheric Research (NCAR) CP-2 radar during the Cooperative Huntsville Meteorological Experiment (COHMEX) and the Convection and Precipitation/Electrification Experiment (CaPE), respectively. ATI studies were accomplished in this work using multiparameter radar data acquired during the lifetime of six convective events that occurred in the COHMEX radar coverage area. A case study from the COMHEX field campaign (20 July 1986) was selected to depict the various stages in the evolution of a storm over which the ATI and rainfall volume computations were performed using multiparameter radar data. Another case study from the CaPE field campaign (12 August 1991) was used to demonstrate the evolution of a convective cell based on differential reflectivity observations.

Raghavan, R.↗

Probing black hole microstate evolution with networks and random walks

We model black hole microstates and quantum tunneling transitions between them with networks and simulate their time evolution using well-established tools in network theory. In particular, we consider two models based on Bena-Warner three-charge multi-centered microstates and one model based on the D1-D5 system; we use network theory methods to determine how many centers (or D1-D5 string strands) we expect to see in a typical late-time state. We find three distinct possible phases in parameter space for the late-time behaviour of these networks, which we call ergodic, trapped, and amplified, depending on the relative importance and connectedness of microstates. We analyze in detail how these different phases of late-time behavior are related to the underlying physics of the black hole microstates. Our results indicate that the expected properties of microstates at late times cannot always be determined simply by entropic arguments; typicality is instead a highly non-trivial, emergent property of the full Hilbert space of microstates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Variable-density buoyancy-driven turbulence with asymmetric initial density distribution

The effects of different initial density distributions on the evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at low (0.05) and high (0.75) Atwood numbers are studied by using high-resolution direct numerical simulations. HVDT aims to mimic the acceleration-driven Rayleigh–Taylor and shock-driven Richtmyer–Meshkov instabilities and reveals new physics that arise from variable-density effects on the turbulent mixing. Here, the initial amounts of pure light and pure heavy flows are altered primarily to mimic the variable-density turbulence at the different locations of the Rayleigh–Taylor and Richtmyer–Meshkov instabilities’ mixing layers where the amounts of the mixing fluids are not equal. It is found that for the low Atwood number cases, the asymmetric initial density distribution has limited effects on both global and local flow evolution for HVDT. However, at high Atwood number, both global flow evolution and the local flow structures are strongly affected by the initial composition ratio. The flow composed of more light fluid reaches higher turbulent levels and the local statistics reach their fully-developed behavior earlier in the time evolution. In contrast, during the late time decay, where most of the flow is well-mixed, all parameters become independent of the initial composition ratio for both low and high Atwood number cases.

42 ENGINEERING↗

Detecting superconductivity out of equilibrium

Recent pump-probe experiments on underdoped cuprates and similar systems suggest the existence of a transient superconducting state above Tc. This poses the question of how to reliably identify the emergence of long-range order, in particular superconductivity, out of equilibrium. We investigate this point by studying a quantum quench in an extended Hubbard model and by computing various observables, which are used to identify (quasi-)long-range order in equilibrium. Our findings imply that, in contrast to current experimental studies, it does not suffice to study the time evolution of the optical conductivity to identify superconductivity. In turn, we suggest to utilize time-resolved angle-resolved photoemission spectroscopy experiments to probe for the formation of a condensate in the two-particle channel.

36 MATERIALS SCIENCE↗

Real-time quantum calculations of phase shifts using wave packet time delays

We present a method to extract the phase shift of a scattering process using the real-time evolution in the early and intermediate stages of the collision in order to estimate the time delay of a wave packet. This procedure is convenient when using noisy quantum computers for which the asymptotic out-state behavior is unreachable. We demonstrate that the challenging Fourier transforms involved in the state preparation and measurements can be implemented in 1+1 dimensions with current trapped ion devices and IBM quantum computers. We compare quantum computations of the time delays obtained in the one-particle quantum mechanics limit and the scalable quantum field theory formulation with accurate numerical results. We discuss the finite volume effects in the Wigner formula connecting time delays to phase shifts. The results reported involve two- and four-qubit calculations, and we discuss the possibility of larger scale computations in the near future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Entanglement entropy production in deep inelastic scattering

Deep inelastic scattering (DIS) samples a part of the wave function of a hadron in the vicinity of the light cone. Lipatov constructed a spin chain which describes the amplitude of DIS in leading logarithmic approximation. Kharzeev and Levin proposed the entanglement entropy as an observable in DIS [Phys. Rev. D 95, 114008 (2017)], and suggested a relation between the entanglement entropy and parton distributions. Here we represent the DIS process as a local quench in Lipatov’s spin chain and study the time evolution of the produced entanglement entropy. We show that the resulting entanglement entropy depends on time logarithmically, $\mathcal{S}(t) = 1/3 ln(t/τ)$ with $τ = 1/m for 1/m ≤ t ≤ (mx)^{–1}$, where m is the proton mass and $\textit{x}$ is the Bjorken $\textit{x}$. The central charge c of Lipatov’s spin chain is determined here to be $\textit{c}$ = 1; using the proposed relation between the entanglement entropy and parton distributions, this corresponds to the gluon structure function growing at small $\textit{x}$ as $xG(x) ~ 1/x^{1/3}$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anelastic response of the earth to a dip slip earthquake

The deformation induced by a vertical dip slip earthquake is examined using a variety of rheologic models. In this way the complications of dipping faults are avoided, and the phenomenon of transient peripheral warping is clearly revealed. A thrust fault dipping at 30 deg is investigated, and the important effects of dip and the existence of a slab on the asymmetry of strain pulses propagated into the overthrust and subducted lithosphere are demonstrated. One of the signal results of the study is the essential similarity of the strain patterns for Newtonian and non-Newtonian flow laws: the two rheologies give nearly identical strain field geometries. The principal difference between the two, which is readily observable, is in their time evolution. Relaxation in non-Newtonian rheologies tends to be initially fast, then slow at times that are late in comparison with relaxation in a Newtonian rheology. The possibility of simply recalling the time dependence of a Newtonian solution to obtain an approximate solution to a non-Newtonian problem is demonstrated.

Melosh, H. J.↗

Time-Variable Gravity Signal Due to Extratropic Pacific Water Mass Redistribution

Cox and Chao [2002] reported the detection of a large anomaly in the form of a positive "jump" in the time series of Earth's lowest-degree gravity harmonic J2, or the dynamic oblateness, during 1998. This prompted us to examine the mass redistribution in the global oceans. We report here a seesaw of the sea-surface height (SSH) in the extratropic north + south Pacific basins -- the leading (nonseasonal) EOF/PC mode in SSH derived from the 10-year TOPEX/Poseidon altimetry data in the extratropic Pacific region. The mode underwent a step-like jump with time evolution that match remarkably well with the observed J2 anomaly. However, the magnitude is several times too small to explain the observed J2, even if assuming the SSH jump was all mass-induced (as opposed to any steric effect which causes no time-variable gravity signal). If one accepts the notion that this extratropic Pacific seesaw is part of the geophysical process that produced the observed 1998 J2 anomaly, then this finding suggests strong geophysical connection of the interannual-to-decadal variation of J2 with the Pacific Decadal Oscillation (PDO), as the time series of the above EOF/PC mode is actually a formally defined PDO Index series.

Chao, Benjamin F.↗

An Analysis of the Orbital Distribution of Solid Rocket Motor Slag

The contribution made by orbiting solid rocket motors (SRMs) to the orbital debris environment is both potentially significant and insufficiently studied. A combination of rocket motor design and the mechanisms of the combustion process can lead to the emission of sufficiently large and numerous by-products to warrant assessment of their contribution to the orbital debris environment. These particles are formed during SRM tail-off, or the termination of burn, by the rapid expansion, dissemination, and solidification of the molten Al2O3 slag pool accumulated during the main burn phase of SRMs utilizing immersion-type nozzles. Though the usage of SRMs is low compared to the usage of liquid fueled motors, the propensity of SRMs to generate particles in the 100 m and larger size regime has caused concern regarding their contributing to the debris environment. Particle sizes as large as 1 cm have been witnessed in ground tests conducted under vacuum conditions and comparable sizes have been estimated via ground-based telescopic and in-situ observations of sub-orbital SRM tail-off events. Using sub-orbital and post recovery observations, a simplistic number-size-velocity distribution of slag from on-orbit SRM firings was postulated. In this paper we have developed more elaborate distributions and emission scenarios and modeled the resultant orbital population and its time evolution by incorporating a historical database of SRM launches, propellant masses, and likely location and time of particulate deposition. From this analysis a more comprehensive understanding has been obtained of the role of SRM ejecta in the orbital debris environment, indicating that SRM slag is a significant component of the current and future population.

Horstman, Matthew F.↗

Jet-associated deuteron production in pp collisions at $\sqrt{s}$ =13 TeV

Deuteron production in high-energy collisions is sensitive to the space–time evolution of the collision system, and is typically described by a coalescence mechanism. For the first time, we present results on jet-associated deuteron production in pp collisions at $\sqrt{s}$ =13 TeV, providing an opportunity to test the established picture for deuteron production in events with a hard scattering. Using a trigger particle with high transverse-momentum (p T > GeV/c) as a proxy for the presence of a jet at midrapidity, we observe a measurable population of deuterons being produced around the jet proxy. The associated deuteron yield measured in a narrow angular range around the trigger particle differs by 2.4–4.8 standard deviations from the uncorrelated background. The data are described by PYTHIA model calculations featuring baryon coalescence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The turbulence driving parameter of molecular clouds in disc galaxies

ABSTRACT Supersonic turbulence plays a pivotal role during the formation of molecular clouds and stars in galaxies. However, little is known about how the fraction of compressive and solenoidal modes in the velocity field evolves over time and how it depends on properties of the molecular cloud or the galactic environment. In this work, we carry out magnetohydrodynamical simulations of disc galaxies and study the time evolution of the turbulence driving parameter for an ensemble of clouds. We find that the time-averaged turbulence driving parameter is insensitive to the position of the cloud within the galaxy. The ensemble-averaged driving parameter is found to be rather compressive with b ∼ 0.5–0.7, indicating almost time-independent global star formation properties. However, each individual cloud shows a highly fluctuating driving parameter, which would strongly affect the cloud’s star formation rate. We find that the mode of turbulence driving can rapidly change within only a few Myr, both from solenoidal to compressive and vice versa. We attribute these changes to cloud collisions and to tidal interactions with clouds or overdensities in the environment. Last, we find no significant differences in the average driving parameter between hydrodynamic and initially strongly magnetized galaxies. However, the magnetic field tends to reduce the overall fluctuation of the driving parameter. The average driving and its uncertainty are seen to be in agreement with recent constraints on the turbulence driving mode for solar neighbourhood clouds.

Körtgen, Bastian↗

Polarity dependent heating at the phase interface in metal-insulator transitions

We show that in a model of a metal photoexcited by a transient pump pulse resonant with a phonon mode, the absolute dc conductivity may become negative, depending on the interplay between the electronic structure, the phonon frequency, and the pump intensity. The analysis includes the effects of inelastic scattering and thermal relaxation. Results for the time evolution of the negative conductivity state are presented; the associated nonequilibrium physics may persist for long times after the pulse. Furthermore, our findings provide a theoretical justification for previously proposed phenomenology and indicate new routes to the generation and exploration of intrinsically nonequilibrium states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiple energetic injections in a strong spike-like solar burst

An intense and fast spike-like solar burst was built up of short time scale structures superimposed on an underlying gradual emission, the time evolution of which shows remarkable proportionality between hard X-ray and microwave fluxes. The finer time structure were best defined at mm-microwaves. At the peak of the event, the finer structures repeat every 30x60ms. The more slowly varying component with a time scale of about 1 second was identified in microwave hard X-rays throughout the burst duration. It is suggested that X-ray fluxes might also be proportional to the repetition rate of basic units of energy injection (quasi-quantized). The relevant parameters of one primary energy release site are estimated both in the case where hard X-rays are produced primarily by thick-target bremsstrahlung, and when they are purely thermal. The relation of this figure to global energy considerations is discussed.

Kaufmann, P.↗