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 37 records · Page 2

Finite-element time evolution operator for the anharmonic oscillator

The finite-element approach to lattice field theory is both highly accurate (relative errors approximately 1/N(exp 2), where N is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly preserved at the lattice sites). In this talk I construct matrix elements for dynamical variables and for the time evolution operator for the anharmonic oscillator, for which the continuum Hamiltonian is H = p(exp 2)/2 + lambda q(exp 4)/4. Construction of such matrix elements does not require solving the implicit equations of motion. Low order approximations turn out to be extremely accurate. For example, the matrix element of the time evolution operator in the harmonic oscillator ground state gives a results for the anharmonic oscillator ground state energy accurate to better than 1 percent, while a two-state approximation reduces the error to less than 0.1 percent.

Milton, Kimball A.↗

Ground state energy and magnetization curve of a frustrated magnetic system from real-time evolution on a digital quantum processor

Models of interacting many-body quantum systems that may realize new exotic phases of matter, notably quantum spin liquids, are challenging to study using even state-of-the-art classical methods such as tensor network simulations. Quantum computing provides a promising route for overcoming these difficulties to find ground states, dynamics, and more. In this paper, we argue that recently developed hybrid quantum-classical algorithms based on real-time evolution are promising methods for solving a particularly important model in the search for spin liquids, the antiferromagnetic Heisenberg model on the two-dimensional kagome lattice. We show how to construct efficient quantum circuits to implement time evolution for the model and to evaluate key observables on the quantum computer, and we argue that the method has favorable scaling with increasing system size. We then restrict to a 12-spin star plaquette from the kagome lattice and a related 8-spin system, and we give an empirical demonstration on these small systems that the hybrid algorithms can efficiently find the ground state energy and the magnetization curve. For these demonstrations, we use four levels of approximation: exact state vectors, exact state vectors with statistical noise from sampling, noisy classical emulators, and (for the 8-spin system only) real quantum hardware, specifically the Quantinuum H1-1 processor; for the noisy simulations and hardware demonstration, we also employ error mitigation strategies based on the symmetries of the Hamiltonian. Our results strongly suggest that these hybrid algorithms present a promising direction for studying quantum spin liquids and more generally for resolving important unsolved problems in condensed matter theory and beyond.

97 MATHEMATICS AND COMPUTING↗

Inductive Pulsed Plasma Thruster Model with Time-Evolution of Energy and State Properties

A model for pulsed inductive plasma acceleration is presented that consists of a set of circuit equations coupled to both a one-dimensional equation of motion and an equation governing the partitioning of energy. The latter two equations are obtained for the plasma current sheet by treating it as a single element of finite volume and integrating the governing equations over that volume. The integrated terms are replaced where necessary by physically-equivalent quantities that are calculated through the solution of other parts of the governing equation set. The model improves upon previous one-dimensional performance models by permitting the time-evolution of the energy and state properties of the plasma, the latter allowing for the tailoring of the model to different gases that may be chosen as propellants. The time evolution of the various energy modes in the system and the associated plasma properties, calculated for argon propellant, are presented to demonstrate the efficacy of the model. The model produces a result where efficiency is maximized at a given value of the electrodynamic scaling term known as the dynamic impedance parameter. Qualitatively and quantitatively, the model compares favorably with performance measured for two separate inductive pulsed plasma thrusters, with disagreements attributable to simplifying assumptions employed in the generation of the model solution.

Polzin, Kurt A.↗

Improving Spectral Resolution from Real-time Evolution for Correlated Systems

Abstract The quality of numerically simulated spectra using real-time evolution methods for strongly correlated systems is affected by both the length of simulation time and the system size, limiting resolution in both frequency and momentum. In this work, we propose a computationally cheap, linear autoregressive machine learning-based framework to extend short-time and short-distance results over a wider range. We use the proposed method to extend the lesser Green’s function for both the Hubbard model and the much more computationally challenging Hubbard-extended Holstein model. This technique significantly improves both the frequency and momentum resolution of the single-particle removal spectrum $${\mathcal{A}}(k,\omega )$$ A ( k , ω ) , allowing the observation of otherwise obscured spectral features due to electron-phonon coupling.

Tang, Ta↗

Time Evolution of Modeled Reynolds Stresses in Planar Homogeneous Flows

The analytic expression of the time evolution of the Reynolds stress anisotropy tensor in all planar homogeneous flows is obtained by exact integration of the modeled differential Reynolds stress equations. The procedure is based on results of tensor representation theory, is applicable for general pressure-strain correlation tensors, and can account for any additional turbulence anisotropy effects included in the closure. An explicit solution of the resulting system of scalar ordinary differential equations is obtained for the case of a linear pressure-strain correlation tensor. The properties of this solution are discussed, and the dynamic behavior of the Reynolds stresses is studied, including limit cycles and sensitivity to initial anisotropies.

Jongen, T.↗

Alternative approach to quantum imaginary time evolution

There is increasing interest in quantum algorithms (QAs) that are based on the imaginary time evolution (ITE), a successful classical numerical approach to obtain ground states. However, most of the proposals so far require heavy postprocessing computational steps on a classical computer, such as solving linear equations. Here we provide an alternative approach to implement ITE. A key feature in our approach is the use of an orthogonal basis set: the propagated state is efficiently expressed in terms of orthogonal basis states at every step of the evolution. We argue that the number of basis states needed at those steps to achieve an accurate solution can be kept on the order of n , the number of qubits, by controlling the precision (number of significant digits) and the imaginary time increment. The number of quantum gates per imaginary time step is estimated to be polynomial in n . Additionally, while in many QAs the locality of the Hamiltonian is a key assumption, in our algorithm this restriction is not required. This characteristic of our algorithm renders it useful for studying highly nonlocal systems, such as the occupation-representation nuclear shell model. Here, we illustrate our algorithm through numerical implementation on an IBM quantum simulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time evolution of pulsed far-field patterns of GaAlAs phase-locked laser-diode arrays

The time evolution of pulsed far fields from commercially available gain-guided phase-locked arrays has been studied and compared to experimental results previously obtained by streak camera measurements, and to predictions of coupled-mode theories developed for index- and gain-guided arrays. Although phase locking is evident by 100 ps into the drive pulse, stable operation in a fixed superposition of array modes is not achieved until 1-2 ns after the drive pulse has stabilized.

Forrest, Kathrine A.↗

The Timing Evolution of 4U 1630-47 During Its 1998 Outburst

We report on the evolution of the timing of 4U 1630-47 during its 1998 outburst using data obtained with the Rossi X-ray Timing Explorer (RXTE). The count rate and position in hardness intensity, color-color diagrams, and simple spectral fits are used to track the concurrent spectral changes. The source showed seven distinct types of timing behavior, most of which show differences with the canonical black hole spectral/timing states. In marked contrast to previous outbursts, we find quasi-periodic oscillation (QPO) signals during nearly all stages of the outburst with frequencies between 0.06 and 14 Hz and a remarkable variety of other characteristics. In particular, we find large (up to 23% rms) amplitude QPOs on the early rise. Later, slow 0.1 Hz semiregular short (-5 s), 9%-16% deep dips dominate the light curve. At this time there are two QPOs, one stable near 13.5 Hz and the other whose frequency drops from the 6-8 Hz range to -4.5 Hz during the dips. BeppoSAX observations during the very late declining phase show 4U 1630-47 in a low state.

Dieters, S. W.↗

The application of the phase space time evolution method to electron shielding

A computer technique for treating the motion of charged and neutral particles and called the phase space time evolution method was developed. This technique employs the computer's bookkeeping capacity to keep track of the time development of a phase space distribution of particles. This method was applied to a study of the penetration of electrons. A 1 MeV beam of electrons normally incident on a semi-infinite slab of aluminum was used. Results of the calculation were compared with Monte Carlo calculations and experimental results. Time-dependent PSTE electron penetration results for the same problem are presented.

Cordaro, M. C.↗

Efficient Step-Merged Quantum Imaginary Time Evolution Algorithm for Quantum Chemistry

In this work, we develop a resource-efficient step-merged quantum imaginary time evolution approach (smQITE) to solve for the ground state of a Hamiltonian on quantum computers. This heuristic method features a fixed shallow quantum circuit depth along the state evolution path. We use this algorithm to determine the binding energy curves of a set of molecules, including H 2 , H 4 , H 6 , LiH, HF, H 2 O, and BeH 2 , and find highly accurate results. The required quantum resources of smQITE calculations can be further reduced by adopting the circuit form of the variational quantum eigensolver (VQE) technique, such as the unitary coupled cluster ansatz. We demonstrate that smQITE achieves a similar computational accuracy as VQE at the same fixed-circuit ansatz, without requiring a generally complicated high-dimensional nonconvex optimization. Finally, smQITE calculations are carried out on Rigetti quantum processing units, demonstrating that the approach is readily applicable on current noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sequency Hierarchy Truncation (SeqHT) for Adiabatic State Preparation and Time Evolution in Quantum Simulations

We introduce the Sequency Hierarchy Truncation (SeqHT) scheme for reducing the resources required for state preparation and time evolution in quantum simulations, based upon a truncation in sequency. For the λϕ 4 interaction in scalar field theory, or any interaction with a polynomial expansion, upper bounds on the contributions of operators of a given sequency are derived. For the systems we have examined, observables computed in sequency-truncated wavefunctions, including quantum correlations as measured by magic, are found to step-wise converge to their exact values with increasing cutoff sequency. The utility of SeqHT is demonstrated in the adiabatic state preparation of the λϕ 4 anharmonic oscillator ground state using IBM's quantum computer ibm_sherbrooke. Using SeqHT, the depth of the required quantum circuits is reduced by ∼ 30 % , leading to significantly improved determinations of observables in the quantum simulations. More generally, SeqHT is expected to lead to a reduction in required resources for quantum simulations of systems with a hierarchy of length scales.

Li, Zhiyao [Univ. of Washington, Seattle, WA (Unit↗

Benchmarking Quantum Chemistry Computations with Variational, Imaginary Time Evolution, and Krylov Space Solver Algorithms

Quantum chemistry is a key application area for noisy-intermediate scale quantum (NISQ) devices, and therefore serves as an important benchmark for current and future quantum computer performance. Previous benchmarks in this field have focused on variational methods for computing ground and excited states of various molecules, including a benchmarking suite focused on the performance of computing ground states for alkali-hydrides under an array of error mitigation methods. State-of-the-art methods to reach chemical accuracy in hybrid quantum-classical electronic structure calculations of alkali hydride molecules on NISQ devices from IBM are outlined here. Here it is demonstrated how to extend the reach of variational eigensolvers with symmetry preserving Ansätze. Next, it is outlined how to use quantum imaginary time evolution and Lanczos as a complementary method to variational techniques, highlighting the advantages of each approach. Finally, a new error mitigation method is demonstrated which uses systematic error cancellation via hidden inverse gate constructions, improving the performance of typical variational algorithms. These results show that electronic structure calculations have advanced rapidly, to routine chemical accuracy for simple molecules, from their inception on quantum computers a few short years ago, and they point to further rapid progress to larger molecules as the power of NISQ devices grows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time evolution of a pumped molecular magnet—A time-resolved inelastic neutron scattering study

Introducing an experimental technique of time-resolved inelastic neutron scattering (TRINS), we explore the time-dependent effects of resonant pulsed microwaves on the molecular magnet Cr 8 F 8 Piv 16 . The octagonal rings of magnetic Cr 3+ atoms with antiferromagnetic interactions form a singlet ground state with a weakly split triplet of excitations at 0.8 meV. A 4.6 tesla field was applied to tune the splitting between two members of the triplet excited level |1$\rangle$ ↔ |2$\rangle$ to resonance with 105 GHz (0.434 meV) microwaves. The time-dependent occupations of the ground state |0$\rangle$, lower lying levels |1$\rangle$ and |2$\rangle$, and higher energy states |λ ≥ 3$\rangle$ were extracted during and after 20 s long microwave pulses incident along the (101) direction of a Cr 8 F 8 Piv 16 crystal held at 1.9 K. At significantly elevated spin temperatures, we found underpopulation relative to thermal equilibrium of |2$\rangle$ and spin-lattice thermalization time scales ranging from 1.6(2) s to 5.7(2) s depending on the power level. This contrasts with the relaxation time τ 1 (T → 0) = 27(5) 𝛍 s inferred for |2$\rangle$ from in situ Electron Spin Resonance measurements. By probing a broad range of excited states during intense microwave pumping, TRINS thus provides a first view of long lived excited states in a molecular antiferromagnet.

Science & Technology - Other Topics↗

Efficient method for estimating the time evolution of the proto-neutron star mass and radius from a supernova neutrino signal

In this paper, we present a novel method to estimate the time evolution of the proto-neutron star (PNS) structure from the neutrino signal in a core-collapse supernova (CCSN). Employing recent results from multidimensional CCSN simulations, we delve into a relation between the total emitted neutrino energy (TONE) and PNS mass/radius, and we find that they are strongly correlated with each other. We fit the relation by simple polynomial functions connecting the TONE to the mass and radius of the PNS as a function of time. By combining another fitting function representing the correlation between the TONE and the cumulative number of events at each neutrino observatory, the PNS mass and radius can be retrieved from purely observed neutrino data. We demonstrate retrievals of PNS mass and radius from mock data of the neutrino signal, and we assess the capability of our proposed method. While underlining the limitations of the method, we also discuss the importance of the joint analysis with the gravitational wave signal. This would reduce uncertainties of parameter estimations in our method, and may narrow down the possible neutrino oscillation model. The proposed method is a very easy and inexpensive computation, which will be useful in real data analysis of the CCSN neutrino signal.

79 ASTRONOMY AND ASTROPHYSICS↗

Selected 3D Flash-X Checkpoints for the Long-Time Evolution of a 9.6 Solar-Mass Core-Collapse Supernova Model

This dataset contains selected 3D Flash-X checkpoint files from the long-time evolution of a low-energy core-collapse supernova explosion of a 9.6 Msun zero-metallicity, low-mass iron-core progenitor. The checkpoints span the shock-breakout phase through the young-remnant phase, ending at approximately 3 yr after core bounce. The dataset is intended for follow-up analysis and post-processing, especially radiation-transport calculations using the hydrodynamic and compositional structure of the ejecta. For a full description of the numerical setup, physical assumptions, limitations, and interpretation of the simulation, users should refer to the associated paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Early-stage Testing of Thermal Power Dispatch Simulator for Subjective Mental Workload and Evolution Time

The excess thermal energy produced by nuclear power plants (NPPs) during low electricity demands can be utilized in industrial processes, such as hydrogen production, through a thermal power dispatch (TPD) system. Initial testing of the first iteration of a single-train TPD design with a manual control mode at the Idaho National Lab (INL) revealed a high operator workload and degraded control capability. The current study evaluated the impact of an enhanced dual-train TPD design on operators’ subjective mental workload and evolution task-time while completing two operating scenarios in manual and automatic control modes. The results showed no statistically significant difference between participants’ mental workload using both control modes. Evolution time in automatic control mode took a shorter time than in manual control, with participants completing all evolutions in less than the 10-min set as the design specification limit. The shorter evolution time is discussed within the context of plant safety and operational efficiency.

Gideon, Olugbenga↗

Plio-Pleistocene time evolution of the 100-ky cycle in marine paleoclimate records

To constrain theories for the dynamical evolution of global ice mass through the late Neogene, it is important to determine whether major changes in the record were gradual or rapid. Of particular interest is the evolution of the near 100-ky ice age cycle in the middle Pleistocene. We have applied a new technique based on multiple taper spectrum analysis which allows us to model the time evolution of quasi-periodic signals. This technique uses both phase and amplitude information, and enables us to address the question of abrupt versus gradual onset of the 100-ky periodicity in the middle Pleistocene.

Park, Jeffrey↗