Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time Evolution of Quantum 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Direct-bonded diamond membranes for heterogeneous quantum and electronic technologies

Diamond has superlative material properties for a broad range of quantum and electronic technologies. However, heteroepitaxial growth of single crystal diamond remains limited, impeding integration and evolution of diamond-based technologies. Here, we directly bond single-crystal diamond membranes to a wide variety of materials including silicon, fused silica, sapphire, thermal oxide, and lithium niobate. Our bonding process combines customized membrane synthesis, transfer, and dry surface functionalization, allowing for minimal contamination while providing pathways for near unity yield and scalability. We generate bonded crystalline membranes with thickness as low as 10 nm, sub-nm interfacial regions, and nanometer-scale thickness variability over 200 by 200 μm 2 areas. We measure spin coherence times T 2 for nitrogen vacancy centers in 150 nm-thick bonded membranes of up to 623 ± 21 μs, suitable for advanced quantum applications. We demonstrate multiple methods for integrating high quality factor nanophotonic cavities with the diamond heterostructures, highlighting the platform versatility in quantum photonic applications. Furthermore, we show that our ultra-thin diamond membranes are compatible with total internal reflection fluorescence (TIRF) microscopy, which enables interfacing coherent diamond quantum sensors with living cells while rejecting unwanted background luminescence. The processes demonstrated herein provide a full toolkit to synthesize heterogeneous diamond-based hybrid systems for quantum and electronic technologies.

color center↗
DOE OSTI2023

Fredkin Staircase: An Integrable System with a Finite-Frequency Drude Peak

We introduce and explore an interacting integrable cellular automaton, the Fredkin staircase, that lies outside the existing classification of such automata, and has a structure that seems to lie beyond that of any existing Bethe-solvable model. The Fredkin staircase has two families of ballistically propagating quasiparticles, each with infinitely many species. Despite the presence of ballistic quasiparticles, charge transport is diffusive in the d.c. limit, albeit with a highly non- gaussian dynamic structure factor. Remarkably, this model exhibits persistent temporal oscillations of the current, leading to a delta-function singularity (Drude peak) in the a.c. conductivity at nonzero frequency. Here, we analytically construct an extensive set of operators that anticommute with the time-evolution operator; the existence of these operators both demonstrates the integrability of the model and allows us to lower-bound the weight of this finite-frequency singularity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗
DOE OSTI2025

Engineering micromotion in Floquet prethermalization via space-time symmetries

We present a systematic framework for Floquet prethermalization under strong resonant driving, emphasizing the pivotal role of dynamical space-time symmetries. Our approach demonstrates how dynamical space-time symmetries map onto the projective static symmetry group of the prethermal Hamiltonian governing the prethermal regime. We introduce techniques for detecting dynamical symmetries through the time evolution of local observables, facilitating a detailed analysis of micromotion within each period and surpassing the limitations of conventional stroboscopic Floquet prethermal dynamics. To implement this framework, we present a prethermal protocol that preserves order-2 dynamical symmetry in a spin-ladder model, confirming the predicted relationships between the expectation values of local observables at distinct temporal points in the Floquet cycle, linked by this symmetry.

floquet systems↗
DOE OSTI2021

Penalized ensemble Kalman filters for high dimensional non-linear systems

The ensemble Kalman filter (EnKF) is a data assimilation technique that uses an ensemble of models, updated with data, to track the time evolution of a usually non-linear system. It does so by using an empirical approximation to the well-known Kalman filter. However, its performance can suffer when the ensemble size is smaller than the state space, as is often necessary for computationally burdensome models. This scenario means that the empirical estimate of the state covariance is not full rank and possibly quite noisy. To solve this problem in this high dimensional regime, we propose a computationally fast and easy to implement algorithm called the penalized ensemble Kalman filter (PEnKF). Under certain conditions, it can be theoretically proven that the PEnKF will be accurate (the estimation error will converge to zero) despite having fewer ensemble members than state dimensions. Further, as contrasted to localization methods, the proposed approach learns the covariance structure associated with the dynamical system. These theoretical results are supported with simulations of several non-linear and high dimensional systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2021

Search for Magnetic Accretion in SW Sextantis Systems

SW Sextantis systems are nova-like cataclysmic variables that have unusual spectroscopic properties, which are thought to be caused by an accretion geometry having part of the mass flux trajectory out of the orbital plane. Accretion onto a magnetic white dwarf is one of the proposed scenarios for these systems. To verify this possibility, we analyzed photometric and polarimetric time-series data for a sample of six SW Sex stars. We report possible modulated circular polarization in BO Cet, SW Sex, and UU Aqr with periods of 11.1, 41.2, and 25.7 minutes, respectively, and less significant periodicities for V380 Oph at 22 minutes and V442 Oph at 19.4 minutes. We confirm previous results that LS Peg shows variable circular polarization. However, we determine a period of 18.8 minutes, which is different from the earlier reported value. We interpret these periods as the spin periods of the white dwarfs. Our polarimetric results indicate that 15% of the SW Sex systems have direct evidence of magnetic accretion. We also discuss SW Sex objects within the perspective of being magnetic systems, considering the latest findings about the demography, formation, and evolution of cataclysmic variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Thermodynamics of continuous spin photons

Special relativity allows massless particles to have states of different integer (or half-integer) helicities that mix under boosts, much like the spin states of a massive particle. Such massless particles are known as “continuous spin” particles (CSPs), a term coined by Wigner, and they are notable for their infinite tower of spin polarizations. The mixing under boosts is controlled by a spin scale 𝜌 with units of momentum. Normally, we assume 𝜌 = 0. The interactions of CSPs are known to satisfy certain simple properties, one of which is that the 𝜌 → 0 limit generically recovers familiar interactions of massless scalars, photons, or gravitons, with all other polarizations decoupling in this limit. Thus, one can ask if the photon of the Standard Model is a CSP at small but nonzero 𝜌. One concern about this possibility—originally raised by Wigner—is that the infinite tower of polarizations could pose problems for thermodynamics. To address this question, we study the thermal evolution of a CSP photon gas coupled to isothermal matter, across CSP helicity modes and phase space. We find that the structure of the interactions dictated by Lorentz symmetry implies well-behaved thermodynamics. When the CSP photon’s interactions with charged matter are turned on, the primary ℎ =±1 helicity modes thermalize quickly, while the other modes require increasingly long time scales to thermalize, set by powers of 𝑇/𝜌. In familiar thermal systems, the CSP photon behaves like the QED photon with small 𝜌- and time-dependent corrections to its effective relativistic degrees of freedom. Sizable departures from familiar thermal behavior arise at energy scales comparable to 𝜌 and could have testable experimental consequences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Non-conformal attractor in boost-invariant plasmas

We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the shear Reynolds number previously observed in massless systems. For generic nonzero particle mass, neither the shear nor the bulk viscous pressure relax quickly to a non-equilibrium attractor; they approach the hydrodynamic limit only late, at small values of the inverse Reynolds numbers. Only the longitudinal pressure, which is a combination of thermal, shear and bulk viscous pressures, continues to show early approach to a far-off-equilibrium attractor, driven by the rapid longitudinal expansion at early times. Second-order dissipative hydrodynamics based on a gradient expansion around locally isotropic thermal equilibrium fails to reproduce this attractor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Pressure-Stabilized MnSb2 with Complex Incommensurate Magnetic Order

Marcasite-type compounds have been proposed as promising hosts of exotic magnetic quantum states, yet experimental realizations in stoichiometric, disorder-free systems remain limited. Here, we report the high-pressure stabilization and magnetic characterization of MnSb2, a marcasite-type compound that is thermodynamically metastable under ambient pressure. Single crystals were synthesized using a cubic multianvil press at 3.3 GPa and 490 °C for 24 h, and powder and single-crystal X-ray diffraction confirm the orthorhombic Pnnm structure. These crystals are stable at ambient pressure for a long time up to between 450 and 500 K. Heat-capacity measurements reveal phase transitions at approximately T0 ∼ 118 K and T1 ∼ 220 K. Neutron diffraction uncovers an unconventional magnetic state below T1 ∼ 220 K. Magnetic powder neutron diffraction refinements reveal possible multiple magnetic configurations that provide comparably acceptable fits to the experimental data. While most solutions are consistent with a spin-density-wave (SDW) description, helical models systematically yield inferior agreement factors. Across a broad range of models, the Mn ordered moment reaches a maximum value of approximately 2 μB and remains predominantly collinear, with minimal canting along the c-axis. At 200 K, the magnetic propagation vector is q = (0, 0.3975, 0.3783); upon cooling, the b component increases toward 0.5, reflecting a temperature-dependent evolution of the modulation. The need for modification of the magnetic model between high and low temperatures further highlights the complex and strongly temperature-dependent nature of the magnetic order in this system. These results establish MnSb2 as a pressure-stabilized marcasite magnet with a tunable, complex magnetic state and a compelling stoichiometric platform for exploring unconventional magnetic behavior, including potential altermagnetism.

Xu, Mingyu [Michigan State University , , , ,; Iow↗
DOE OSTI2024

Coexistence of Interacting Charge Density Waves in a Layered Semiconductor

Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe 4 , a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect FS nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that mono- and bi-layers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Finally, our findings not only offer microscopic insights into the correlated ground state of EuTe 4 but also provide a general non-equilibrium approach to understand coexisting, layer-dependent orders in a complex system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Flow-driven spectral chaos (FSC) method for simulating long-time dynamics of arbitrary-order non-linear stochastic dynamical systems

Uncertainty quantification techniques such as the time-dependent generalized polynomial chaos (TD-gPC) use an adaptive orthogonal basis to better represent the stochastic part of the solution space (aka random function space) in time. However, because the random function space is constructed using tensor products, TD-gPC-based methods are known to suffer from the curse of dimensionality. Here, we introduce a new numerical method called the flow-driven spectral chaos (FSC) which overcomes this curse of dimensionality at the random-function-space level. The proposed method is not only computationally more efficient than existing TD-gPC-based methods but is also far more accurate. The FSC method uses the concept of enriched stochastic flow maps to track the evolution of a finite-dimensional random function space efficiently in time. To transfer the probability information from one random function space to another, two approaches are developed and studied herein. In the first approach, the probability information is transferred in the mean-square sense, whereas in the second approach the transfer is done exactly using a new theorem that was developed for this purpose. The FSC method can quantify uncertainties with high fidelity, especially for the long-time response of stochastic dynamical systems governed by ODEs of arbitrary order. Six representative numerical examples, including a nonlinear problem (the Van-der-Pol oscillator), are presented to demonstrate the performance of the FSC method and corroborate the claims of its superior numerical properties. Finally, a parametric, high-dimensional stochastic problem is used to demonstrate that when the FSC method is used in conjunction with Monte Carlo integration, the curse of dimensionality can be overcome altogether.

(nonlinear) stochastic dynamical systems↗
DOE OSTI2020

Distinguishing Electron and Hole Dynamics in Functionalized CdSe/CdS Core/Shell Quantum Dots Using Complementary Ultrafast Spectroscopies and Kinetic Modeling

The evolution of excitation energy and photogenerated charges in semiconductor quantum dots (QDs) functionalized with molecular acceptors can be probed on ultrafast time scales using techniques such as transient absorption (TA) spectroscopy. However, historical interpretations that the 1S(e)-1S 3/2 (h) transition in Cd-chalcogenide QDs is fully attributable to electrons may be misleading, and multiexponential models used to fit TA kinetics do not correspond directly to specific photophysical processes. Here, we present visible-wavelength and mid-IR TA and time-resolved photoluminescence measurements to inform a comprehensive kinetic model of the photoexcited CdSe/CdS core/shell QDs functionalized with passivating oleic acid (OA), hole-accepting ferrocene, or electron-accepting naphthalene bisimide (NBI). We show that ~30% of the 1S signal and 72% of the IR signal can originate from holes in well-passivated core/shell QDs. We also demonstrate evidence of electron trapping in OA-capped core/shell QDs, with additional electron transfer and hole trapping in the QDs functionalized with NBI. Electron (hole) trapping and detrapping occur in 450 ± 100 ps (430 ± 70 ps) and 340 ± 100 ps (1.1 ± 0.4 ns) respectively, while the time constant for electron transfer to NBI is ~1.8 ns. The comprehensive picture of photophysical processes provided by the complementary ultrafast techniques and kinetic modeling can accelerate both the fundamental science and application development of nanostructured and molecular systems.

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

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗