Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initial state”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Efficient quantum circuits based on the quantum natural gradient

Efficient preparation of arbitrary entangled quantum states is crucial for quantum computation. This is particularly important for noisy intermediate-scale quantum simulators relying on variational hybrid quantum-classical algorithms. To that end, we propose symmetry-conserving modified quantum approximate optimization algorithm (SCom-QAOA) circuits. The depths of these circuits depend not only on the desired fidelity to the target state but also on the amount of entanglement the state contains. The parameters of the SCom-QAOA circuits are optimized using the quantum natural gradient method based on the Fubini-Study metric. The SCom-QAOA circuit transforms an unentangled state into a ground state of a gapped one-dimensional Hamiltonian with a circuit depth that depends not on the system size but rather on the finite correlation length. In contrast, the circuit depth grows proportionally to the system size for preparing low-lying states of critical one-dimensional systems. Even in the latter case, SCom-QAOA circuits with depth less than the system size were sufficient to generate states with fidelity in excess of 99%, which is relevant for near-term applications. The proposed scheme enlarges the set of the initial states accessible for variational quantum algorithms and widens the scope of investigation of nonequilibrium phenomena in quantum simulators. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantifying sources of subseasonal prediction skill in CESM2

Abstract Subseasonal prediction fills the gap between weather forecasts and seasonal outlooks. There is evidence that predictability on subseasonal timescales comes from a combination of atmosphere, land, and ocean initial conditions. Predictability from the land is often attributed to slowly varying changes in soil moisture and snowpack, while predictability from the ocean is attributed to sources such as the El Niño Southern Oscillation. Here we use a set of subseasonal reforecast experiments with CESM2 to quantify the respective roles of atmosphere, land, and ocean initial conditions on subseasonal prediction skill over land. These reveal that the majority of prediction skill for global surface temperature in weeks 3–4 comes from the atmosphere, while ocean initial conditions become important after week 4, especially in the Tropics. In the CESM2 subseasonal prediction system, the land initial state does not contribute to surface temperature prediction skill in weeks 3–6 and climatological land conditions lead to higher skill, disagreeing with our current understanding. However, land-atmosphere coupling is important in week 1. Subseasonal precipitation prediction skill also comes primarily from the atmospheric initial condition, except for the Tropics, where after week 4 the ocean state is more important.

Richter, Jadwiga H. (ORCID:0000000170480781)↗
DOE OSTI2026

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Phase transitions and dynamics of the universe

The restoration of symmetry at grand unification in a closed contracting Robertson-Walker universe could slow down and halt the contraction, causing the universe to bounce and avoid the singular state or the big crunch. During subsequent expansion the symmetry is broken and the universe reverts back to its present state. A scenario is presented whereby the present entropy of the substance is created through a series of such bounces along with three possible initial states. The monopole, horizon, and flatness difficulties faced in ordinary Big Bang models are also solved by the model presented.

Petrosian, V.↗
NASA NTRS1981

Modification of average coronal properties in the presence of periodic temperature and density variations near the base

Conductive damping of the transient motion of the solar corona, and the consequences on the average flow field of conductively re-distributed energy are discussed. A particular example is treated. First, a steady solution for solar wind flow is found for a given set of steady boundary conditions, between 1.4 and 28 solar radii. These boundary conditions are that the density is 1.6xEO6 cm-3 and the temperature is 1.6xEO6 degrees. The flow profiles for this initial state are those shown at t=0 in the figures. Then, at t=0, periodic, in phase, sinusoidal variations in temperature and density are initiated and continued indefinitely. The amplitudes of the variations are 5 percent and 7.5 percent for the temperature and density respectively, imposed at 1.4 solar radii. The variations have a period of 1 hour shorter than a coronal transit time, but sufficiently long for the disturbances to propagate for small distances (2 to 3 solar radii). These oscillations are like those that might occur for acoustic oscillations, although no explanation is offered as to how an acoustic oscillation might extend to this radius in the first place.

Suess, S. T.↗
NASA NTRS2003

Reactive Resonances in N+N2 Exchange Reaction

Rich reactive resonances are found in a 3D quantum dynamics study of the N + N2 exchange reaction using a recently developed ab initio potential energy surface. This surface is characterized by a feature in the interaction region called Lake Eyring , that is, two symmetric transition states with a shallow minimum between them. An L2 analysis of the quasibound states associated with the shallow minimum confirms that the quasibound states associated with oscillations in all three degrees of freedom in Lake Eyring are responsible for the reactive resonances in the state-to-state reaction probabilities. The quasibound states, mostly the bending motions, give rise to strong reasonance peaks, whereas other motions contribute to the bumps and shoulders in the resonance structure. The initial state reaction probability further proves that the bending motions are the dominating factors of the reaction probability and have longer life times than the stretching motions. This is the first observation of reactive resonances from a "Lake Eyring" feature in a potential energy surface.

Wang, Dunyou↗
NASA NTRS2008

GEWEX Water and Energy Budget Study

Closing the global water and energy budgets has been an elusive Global Energy and Water-cycle Experiment (GEWEX) goal. It has been difficult to gather many of the needed global water and energy variables and processes, although, because of GEWEX, we now have globally gridded observational estimates for precipitation and radiation and many other relevant variables such as clouds and aerosols. Still, constrained models are required to fill in many of the process and variable gaps. At least there are now several atmospheric reanalyses ranging from the early National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and NCEP/Department of Energy (DOE) reanalyses to the more recent ERA40 and JRA-25 reanalyses. Atmospheric constraints include requirements that the models state variables remain close to in situ observations or observed satellite radiances. This is usually done by making short-term forecasts from an analyzed initial state; these short-term forecasts provide the next guess, which is corrected by comparison to available observations. While this analysis procedure is likely to result in useful global descriptions of atmospheric temperature, wind and humidity, there is no guarantee that relevant hydroclimate processes like precipitation, which we can observe and evaluate, and evaporation over land, which we cannot, have similar verisimilitude. Alternatively, the Global Land Data Assimilation System (GLDAS), drives uncoupled land surface models with precipitation, surface solar radiation, and surface meteorology (from bias-corrected reanalyses during the study period) to simulate terrestrial states and surface fluxes. Further constraints are made when a tuned water balance model is used to characterize the global runoff observational estimates. We use this disparate mix of observational estimates, reanalyses, GLDAS and calibrated water balance simulations to try to characterize and close global and terrestrial atmospheric and surface water and energy budgets to within 10-20% for long term (1986-1995), large-scale global to regional annual means.

Roads, J.↗
DOE OSTI2022

Quantum many-body scars from Einstein-Podolsky-Rosen states in bilayer systems

Quantum many-body scar states are special eigenstates of nonintegrable models with distinctive entanglement features that give rise to infinitely long-lived coherent dynamics under quantum quenches from certain initial states. Here, we elaborate on a construction of quantum many-body scar states in which they emerge from Einstein-Podolsky-Rosen states in systems with two layers, wherein the two layers are maximally entangled. We apply this construction to spin systems as well as systems of itinerant fermions and bosons and demonstrate how symmetries can be harnessed to enhance its versatility. We show that several well-known examples of quantum many-body scars, including the tower of states in the spin-1 XY model and the η-pairing states in the Fermi-Hubbard model, can be understood within this formalism. We also demonstrate how an infinite tower of many-body scar states can emerge in bilayer Bose-Hubbard models with charge conservation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1987

Quantized synchrotron radiation in strong magnetic fields

Spectra from electrons making synchrotron transitions in high magnetic fields are calculated using the exact quantum transition rates for Landau states up to n = 500, and the results are compared with spectra calculated from the classical formula and the asymptotic quantum formula. The behavior of the transition rates are examined at low and high harmonics as a function of the spin state of the electron. The calculations confirm those of Herold, Ruder, and Wunner (1982) for low Landau states, extending them to higher states and individual spin state transitions. The results also confirm the dominance of ground-state transitions at high field strengths noted by White (1976). Single particle emissivities for electrons with both spin-up and spin-down in the initial state are calculated using these transition rates. Spectra for thermal electron distributions at transrelativistic temperatures and for steady state injection of monoenergetic electrons with isotropic pitch angles are also presented.

Harding, Alice K.↗
NASA NTRS1976

Wind laws for shockless initialization

A system of diagnostic equations for the velocity field, or wind laws, was derived for each of a number of models of large-scale atmospheric flow. The derivation in each case is mathematically exact and does not involve any physical assumptions not already present in the prognostic equations, such as nondivergence or vanishing of derivatives of the divergence. Therefore, initial states computed by solving these diagnostic equations should be compatible with the type of motion described by the prognostic equations of the model and should not generate initialization shocks when inserted into the model. Numerical solutions of the diagnostic system corresponding to a barotropic model are exhibited. Some problems concerning the possibility of implementing such a system in operational numerical weather prediction are discussed.

Ghil, M.↗
DOE OSTI2021

Quantum Alternating Operator Ansatz (QAOA) Phase Diagrams and Applications for Quantum Chemistry

Determining Hamiltonian ground states and energies is a challenging task with many possible approaches on quantum computers. While variational quantum eigensolvers are popular approaches for near term hardware, adiabatic state preparation is an alternative that does not require noisy optimization of parameters. Beyond adiabatic schedules, QAOA is an important method for optimization problems. In this work we modify QAOA to apply to finding ground states of molecules and empirically evaluate the modified algorithm on several molecules. This modification applies physical insights used in classical approximations to construct suitable QAOA operators and initial state. We find robust qualitative behavior for QAOA as a function of the number of steps and size of the parameters, and demonstrate this behavior also occurs in standard QAOA applied to combinatorial search. To this end we introduce QAOA phase diagrams that capture its performance and properties in various limits. In particular we show a region in which non-adiabatic schedules perform better than the adiabatic limit while employing lower quantum circuit depth. We further provide evidence our results and insights also apply to QAOA applications beyond chemistry.

Kremenetski, Vladimir↗
NASA NTRS1981

Optimal control and controller location for distributed parameter elastic systems

A class of systems governed by second order partial differential equations and driven by controllers located at points r sub i, i = 1, ..., k, is considered. A cost functional quadratic in the time derivative of the state and the control is associated with the system. The optimal controller locations are defined as the ones which minimize the maximum of the cost over all possible initial states. An analytical solution to the associated Riccati equation is presented providing a convenient expression for determining the optimal locations.

Hamidi, M.↗
NASA NTRS1997

On the Possibility of Ill-Conditioned Covariance Matrices in the First-Order Two-Step Estimator

The first-order two-step nonlinear estimator, when applied to a problem of orbital navigation, is found to occasionally produce first step covariance matrices with very low eigenvalues at certain trajectory points. This anomaly is the result of the linear approximation to the first step covariance propagation. The study of this anomaly begins with expressing the propagation of the first and second step covariance matrices in terms of a single matrix. This matrix is shown to have a rank equal to the difference between the number of first step states and the number of second step states. Furthermore, under some simplifying assumptions, it is found that the basis of the column space of this matrix remains fixed once the filter has removed the large initial state error. A test matrix containing the basis of this column space and the partial derivative matrix relating first and second step states is derived. This square test matrix, which has dimensions equal to the number of first step states, numerically drops rank at the same locations that the first step covariance does. It is formulated in terms of a set of constant vectors (the basis) and a matrix which can be computed from a reference trajectory (the partial derivative matrix). A simple example problem involving dynamics which are described by two states and a range measurement illustrate the cause of this anomaly and the application of the aforementioned numerical test in more detail.

Garrison, James L.↗
DOE OSTI2024

MHD stability of spherical tokamak equilibria with non-monotonic q-profiles

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n≥3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The “double tearing mode” (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Formation of Supermassive Black Holes in Galactic Bulges: A Rotating Collapse Model Consistent with the M(sub BH-sigma) Relation

Motivated by the observed correlation between black hole masses M(sub BH) and the velocity dispersion sigma of host galaxies, we develop a theoretical model of black hole formation in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an initial state specified by a uniform rotation rate OMEGA and a density distribution of the form rho = a(sup 2)(sub eff)per2piGR(sup 2)(so that a(sub eff)is an effective transport speed). The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole (for Schwarzschild geometry). This model reproduces the observed correlation between the estimated black hole masses and the velocity dispersions of galactic bulges, i.e., M(sub BH) approximately equal to 10(sup 8) solar mass(sigma per 200 kilometers per second)(sup 4) where sigma = the square root of 2a(sub eff). To obtain this normalization, the rotation rate OMEGA approximately equal to 2 x 10(exp -15) rad per second. The model also defines a bulge mass scale M(sub B). If we identify the scale M(sub B) with the bulge mass, the model determines the ratio mu(sub B) of black hole mass to the host mass: mu(sub B) approximately equal to 0.0024(sigma per 200 kilometer per second), again in reasonable agreement with observed values. In this scenario, supermassive black holes form quickly (in approximately 10(exp 5) yr) and are born rapidly rotating (with a per M approximately 0.9). This paper also shown how these results depend on the assumed initial conditions; the most important quantity is the initial distribution of specific angular momentum in the precollapse state.

Adams, Fred C.↗
DOE OSTI2021

Filtering states with total spin on a quantum computer

Starting from a general wave-function described on a set of spins or qubits, we propose several quantum algorithms to extract the components of this state on eigenstates of total spin S 2 and its z component S z . The method plays the role of total spin projection and gives access to the amplitudes of the initial state on a total spin basis. Here, the different algorithms have various degrees of sophistication depending on the requested tasks. They can either solely project onto the subspace with good total spin or completely uplift the degeneracy in this subspace. In the former case, when the projection on the total spin j is made, we show that the number of operations for the projection can be reduced from a quadratic to a linear dependence in j. After each measurement, the state collapses to one of the spin eigenstates that could be used for postprocessing. For this reason, we call the method total quantum spin filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Role of Sea‐Surface Salinity in Simulating Historical Decadal Variations of Atlantic Meridional Overturning Circulation in a Coupled Climate Model

Abstract Extending climate prediction from seasonal to decadal timescales requires realistic initialization of not only upper ocean heat content but also the Atlantic meridional overturning circulation (AMOC). However, it remains a major challenge to realistically initialize AMOC in a coupled system while also maintaining a balanced atmosphere‐ocean initial state. This study demonstrates the feasibility of generating fully coupled historical states with realistic AMOC variability. Employing a forced ocean—sea‐ice (FOSI) model simulation as the “truth,” we show reproducibility of key features of historical AMOC decadal variability in a fully coupled model by restoring sea‐surface salinity, in addition to sea‐surface temperature restoring widely used in seasonal prediction. The atmospheric state of the restored coupled model solution is much closer to that of the free coupled simulation than to the observations used in FOSI, pointing to potential advantages of using this approach for initializing decadal predictions with reduced ocean initialization shock.

54 ENVIRONMENTAL SCIENCES↗