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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of Nonlinearities in Physics and Demography

Nonlinearities appear in almost all systems. Earlier, we focused on those in plasmas, ionospheric scattering, and the world population. As turned out, the estimate of the population growth made in 1974 is in astonishing agreement with the United Nations estimates and agrees with our present data to within 2%. Here, a particularly important role, both for the population evolution and wave interaction in plasmas, is played by non-Markovian effects (effects depending on the past time). For the population growth, this occurs due to a delay of one generation in the set of population limiting actions, while, for plasmas, it is caused by nonlinear frequency shifts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗
DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Fermionic mean-field dynamics for spin systems beyond free fermions

We introduce the fermionized time-dependent Hartree–Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions, and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

Dutta, Rishab↗
DOE OSTI2024

Effect of causality constraints on Bayesian analyses of heavy-ion collisions

There have long been questions about the limits to the validity of relativistic fluid dynamics and whether it is being used outside its regime of validity in modern simulations of relativistic heavy-ion collisions. An important new tool for answering this question is a causality analysis in the nonlinear regime—if the solutions of the evolution equations do not respect relativistic causality, then they are not a faithful representation of the underlying relativistic theory (in this case, quantum chromodynamics). Using this nonlinear criterion, it has recently been shown that hydrodynamics is indeed being used outside its regime of validity in simulations, at least sometimes. Here we explore the phenomenological implications, particularly the quantitative effects of demanding limits on acausality in modern Bayesian parameter estimation. We find that, while typically only a small fraction of the system's energy is initially in an acausal regime, placing strict limits on the allowed energy fraction significantly changes the preferred properties of the initial condition, which in turn alters the extracted medium properties such as bulk viscosity, where large values are no longer favored. Furthermore, these findings highlight the importance of developing better theoretical descriptions of the early-time, out-of-equilibrium dynamics of relativistic heavy-ion collisions.

Bayesian methods↗
DOE OSTI2024

Local probe of bulk and edge states in a fractional Chern insulator

The fractional quantum Hall effect is a key example of topological quantum many-body phenomena, arising from the interplay between strong electron correlation, topological order and time-reversal symmetry breaking. Recently, a lattice analogue of the fractional quantum Hall effect at zero magnetic field has been observed, confirming the existence of a zero-field fractional Chern insulator (FCI). Despite this, the bulk–edge correspondence—a hallmark of a FCI featuring an insulating bulk with conductive edges—has not been directly observed. In fact, this correspondence has not been visualized in any system for fractional states owing to experimental challenges. Here we report the imaging of FCI edge states in twisted MoTe 2 (t-MoTe 2 ) using microwave impedance microscopy. By tuning the carrier density, we observe the system evolving between metallic and FCI states, the latter of which exhibits insulating bulk and conductive edges, as expected from the bulk–boundary correspondence. Further analysis suggests the composite nature of the FCI edge states. We also observe the evolution of edge states across the topological phase transition as a function of interlayer electric field and reveal exciting prospects of neighbouring domains with different fractional orders. Furthermore, these findings pave the way for research into topologically protected one-dimensional interfaces between various anyonic states at zero magnetic field, such as gapped one-dimensional symmetry-protected phases with non-zero topological entanglement entropy, Halperin–Laughlin interfaces and the creation of non-abelian anyons.

Imaging techniques↗
DOE OSTI2021

Origin and limit of the recovery of damaged information by time reversal

Recently it was found that scrambled information can be partially recovered by a time-reversed evolution, even after being damaged by an intruder. Here we reconsider the origin of the information recovery, and argue that the presence of classical chaos does not preclude it and only leads to a quantitative reduction of the recovery ratio. We also show how decoherence (i.e., entanglement with the intruder) limits the recovery, by proving an upper bound on the recovery ratio in terms of the entangling power of the intruder's action.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Theory and modeling of molecular modes in the NMR relaxation of fluids

Traditional theories of the nuclear magnetic resonance (NMR) autocorrelation function for intra-molecular dipole pairs assume a single-exponential decay, yet the calculated autocorrelation of realistic systems displays a rich, multi-exponential behavior, resulting in anomalous NMR relaxation dispersion (i.e., frequency dependence). We develop an approach to model and interpret the multi-exponential intra-molecular autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational diffusion and translational diffusion in the same framework. Here, we recast the problem of evaluating the autocorrelation in terms of averaging over a diffusion propagator whose evolution is described by a Fokker–Planck equation. The time-independent part admits an eigenfunction expansion, allowing us to write the propagator as a sum over modes. Each mode has a spatial part that depends on the specified eigenfunction and a temporal part that depends on the corresponding eigenvalue (i.e., correlation time) with a simple, exponential decay. The spatial part is a probability distribution of the dipole pair, analogous to the stationary states of a quantum harmonic oscillator. Drawing inspiration from the idea of inherent structures in liquids, we interpret each of the spatial contributions as a specific molecular mode. These modes can be used to model and predict the NMR dipole–dipole relaxation dispersion of fluids by incorporating phenomena on the molecular level. We validate our statistical mechanical description of the distribution in molecular modes with molecular dynamics simulations interpreted without any relaxation models or adjustable parameters: the most important poles in the Padé–Laplace transform of the simulated autocorrelation agree with the eigenvalues predicted by the theory

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

Chemical Analysis of Tantalum Thin Films compared to Niobium Thin Films using SIMS and XPS

Superconducting qubits are a leading platform for quantum computation. These circuits are typically made from superconducting materials like aluminum or niobium. However, the amorphous niobium oxide and aluminum oxide on the surface of these circuits introduce considerable RF loss due to the presence of two-level systems (TLS), which limits the maximum coherence times T1 to ~100 μs. Capping the niobium qubits with a tantalum layer leads to a 3- 5x improvement [Bal et al., 2024]. But even in this case, tantalum forms an amorphous surface oxide that introduces loss. In an effort to devise strategies to eliminate the presence of this oxide, we present a comprehensive study on the nature of Ta oxide using x-ray photoemission spectroscopy (XPS) and secondary ion mass spectrometry (ToF-SIMS) as a function of heat treatment. The thin films were annealed in ultra-high vacuum conditions and analyzed in situ to characterize the composition and evolution of the native tantalum oxide layer and oxide-metal interface. Our analysis reveals two critical differences between tantalum and niobium oxides: Nb2O5 completely dissolves at 400°C, while Ta2O5 persists even at 800°C. Additionally, tantalum oxide contains only a single suboxide (TaO), in contrast to niobium's two suboxides (NbO and NbO2). The suboxide of tantalum contributes minimally to the total oxide content and shows a relative increase with temperature. Understanding the oxide’s behavior will open new pathways for optimizing coherence times in tantalum qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Out-of-time-order correlators bridge classical transport and quantum dynamics

The out-of-time-order correlator (OTOC) has emerged as a central tool for quantifying decoherence across wide-ranging physical platforms. Here, we demonstrate its direct measurement in a classical ensemble using nuclear magnetic resonance with a modulated gradient spin echo sequence and extend the method into a multidimensional correlation to track exchange phenomena. Position is encoded through magnetic field gradients and momentum through the velocity autocorrelation function, enabling experimental access to OTOCs for proton motion confined within the self-similar lattice of the metal–organic framework MOF-808. Here, water confined to specified geometries within the MOF pores gives rise to spatially distinct diffusive eigenmodes with characteristic relative entropies. We demonstrate that periodic radio frequency driving combined with gradient modulation yields entropy evolution through the selection of distinct diffusion modes. Frequency-resolved diffusion spectra connect these entropy dynamics to classical heat exchange laws, revealing how operational features of quantum systems are mirrored in confined, macroscopic spin ensembles.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2023

Mixed quantum/classical theory for rotational energy exchange in symmetric-top-rotor + linear-rotor collisions and a case study of the ND 3 + D 2 system

The extension of mixed quantum/classical theory (MQCT) to describe collisional energy transfer is developed for a symmetric-top-rotor + linear-rotor system and is applied to ND 3 + D 2 . State-to-state transition cross sections are computed in a broad energy range for all possible processes: when both ND 3 and D 2 molecules are excited or both are quenched, when one is excited while the other is quenched and vice versa, when the ND 3 state changes its parity while D 2 is excited or quenched, and when ND 3 is excited or quenched while D 2 remains in the same state, ground or excited. In all these processes the results of MQCT are found to approximately satisfy the principle of microscopic reversibility. For a set of sixteen state-to-state transitions available from the literature for a collision energy of 800 cm –1 the values of cross sections predicted by MQCT are within 8% of accurate full-quantum results. A useful time-dependent insight is obtained by monitoring the evolution of state populations along MQCT trajectories. It is shown that, if before the collision, D 2 is in its ground state, the excitation of ND 3 rotational states proceeds through a two-step mechanism in which the kinetic energy of molecule–molecule collision is first used to excite D 2 and only then is transferred to the excited rotational states of ND 3 . Finally, it is found that both potential coupling and Coriolis coupling play important roles in ND 3 + D 2 collisions.

inelastic cross-sections↗
DOE OSTI2021

Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators

The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate kx, but the oscillation and damping time scales are longer in local representations, particularly for small kx approaching the system size.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Long-lived period-doubled edge modes of interacting and disorder-free Floquet spin chains

Floquet spin chains have been a venue for understanding topological states of matter that are qualitatively different from their static counterparts by, for example, hosting π edge modes that show stable period-doubled dynamics. However the stability of these edge modes to interactions has traditionally required the system to be many-body localized in order to suppress heating. In contrast, here we show that even in the absence of disorder, and in the presence of bulk heating, π edge modes are long lived. Their lifetime is extracted from exact diagonalization and is found to be non-perturbative in the interaction strength. A tunneling estimate for the lifetime is obtained by mapping the stroboscopic time-evolution to dynamics of a single particle in Krylov subspace. In this subspace, the π edge mode manifests as the quasi-stable edge mode of an inhomogeneous Su-Schrieffer-Heeger model whose dimerization vanishes in the bulk of the Krylov chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mapping the nonequilibrium order parameter of a quasi-two dimensional charge density wave system

The driving force of a charge density wave (CDW) transition in quasi-two dimensional systems is still debated, while being crucial in understanding electronic correlation in such materials. Here we use femtosecond time- and angle-resolved photoemission spectroscopy combined with computational methods to investigate the coherent lattice dynamics of a prototypical CDW system. The photo-induced temporal evolution of the periodic lattice distortion associated with the amplitude mode reveals the dynamics of the free energy functional governing the order parameter. Our approach establishes that optically-induced screening rather than CDW melting at the electronic level leads to a transiently modified potential which explains the anharmonic behaviour of the amplitude mode and discloses the structural origin of the symmetry-breaking phase transition.

electronic properties and materials↗
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantum simulation of real-space dynamics

Quantum simulation is a prominent application of quantum computers. While there is extensive previous work on simulating finite-dimensional systems, less is known about quantum algorithms for real-space dynamics. We conduct a systematic study of such algorithms. In particular, we show that the dynamics of a d-dimensional Schrödinger equation with η particles can be simulated with gate complexity O ~ (ηdFpoly(log(g'/ϵ))), where ϵ is the discretization error, g' controls the higher-order derivatives of the wave function, and F measures the time-integrated strength of the potential. Compared to the best previous results, this exponentially improves the dependence on ϵ and g' from poly(g'/ϵ) to poly(log(g'/ϵ)) and polynomially improves the dependence on T and d, while maintaining best known performance with respect to η. For the case of Coulomb interactions, we give an algorithm using η 3 (d + η)Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates, and another using η 3 (4d) d/2 Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates and QRAM operations, where T is the evolution time and the parameter Δ regulates the unbounded Coulomb interaction. We give applications to several computational problems, including faster real-space simulation of quantum chemistry, rigorous analysis of discretization error for simulation of a uniform electron gas, and a quadratic improvement to a quantum algorithm for escaping saddle points in nonconvex optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Strong Zero Modes in Integrable Quantum Circuits

It is a classic result that certain interacting integrable spin chains host robust edge modes known as strong zero modes (SZMs). In this Letter, we extend this result to the Floquet setting of local quantum circuits, focusing on a prototypical model providing an integrable Trotterization for the evolution of the XXZ Heisenberg spin chain. By exploiting the algebraic structures of integrability, we show that an exact SZM operator can be constructed for these integrable quantum circuits in certain regions of parameter space. Our construction, which recovers a well-known result by Paul Fendley in the continuous-time limit, relies on a set of commuting transfer matrices known from integrability, and allows us to easily prove important properties of the SZM, including normalizabilty. Our approach is different from previous methods and could be of independent interest even in the Hamiltonian setting. Furthermore, our predictions, which are corroborated by numerical simulations of infinite-temperature autocorrelation functions, are potentially interesting for implementations of the XXZ quantum circuit on available quantum platforms.

1-dimensional spin chains↗
DOE OSTI2022

Impact of Nanoparticle Size and Surface Chemistry on Peptoid Self-Assembly

Self-assembled organic nanomaterials can be generated by bottom-up assembly pathways where the structure is controlled by the organic sequence and altered using pH, temperature, and solvation. In contrast, self-assembled structures based on inorganic nanoparticles typically rely on physical packing and drying effects to achieve uniform superlattices. By combining these two chemistries to access inorganic–organic nanostructures, we aim to understand the key factors that govern the assembly pathway and structural outcomes in hybrid systems. In this work, we outline two assembly regimes between quantum dots (QDs) and reversibly binding peptoids. These regimes can be accessed by changing the solubility and size of the hybrid (peptoid-QD) monomer unit. The hybrid monomers are prepared via ligand exchange and assembled, and the resulting assemblies are studied using ex-situ transmission electron microscopy as a function of assembly time. In aqueous conditions, QDs were found to stabilize certain morphologies of peptoid intermediates and generate a final product consisting of multilayers of small peptoid sheets linked by QDs. The QDs were also seen to facilitate or inhibit assembly in organic solvents based on the relative hydrophobicity of the surface ligands, which ultimately dictated the solubility of the hybrid monomer unit. Increasing the size of the QDs led to large hybrid sheets with regions of highly ordered square-packed QDs. A second, smaller QD species can also be integrated to create binary hybrid lattices. Furthermore, these results create a set of design principles for controlling the structure and structural evolution of hybrid peptoid-QD assemblies and contribute to the predictive synthesis of complex hybrid matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗