Search NASA⌕ Search

SEARCH · Search NASA

Results for “Quantum Dynamics”

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 613 records · Page 34

Sparse adaptive basis set methods for solution of the time dependent Schrodinger equation

Scalable numerical solutions to the time dependent Schrodinger equation remain an outstanding goal in theoretical chemistry. Here we present a method which utilises recent breakthroughs in signal processing to consistently adapt a dictionary of basis functions to the dynamics of the system. Further, we show that for two low-dimensional model problems the size of the basis set does not grow quickly with time and appears only weakly dependent on dimensionality. The generality of this finding remains to be seen. The method primarily uses energies and gradients of the potential, opening the possibility for its use in on-the-fly ab initio quantum wavepacket dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Real-Time Time-Dependent Electronic Structure Theory

Real-time electronic structure methods provide an unprecedented view of electron dynamics and ultrafast spectroscopy on the atto- and femtosecond time scale with vast potential to yield new insights into the electronic behavior of molecules and materials. Here, we discuss the fundamental theory underlying various real-time electronic structure methods as well as advantages and disadvantages of each. We give an overview of the numerical techniques that are widely used for real-time propagation of the quantum electron dynamics with an emphasis on Gaussian basis set methods. We also showcase many of the chemical applications and scientific advances made by using real-time electronic structure calculations and provide an outlook of possible new directions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗

XFEL SASE pulses can enhance time-dependent observables

Abstract X-ray free electron lasers (XFELs) have emerged as powerful sources of short and intense x-ray pulses. We propose a simple and robust procedure which takes advantage of the inherent stochasticity of self-amplified stimulated emission (SASE) pulses to enhance the time-resolution and signal strength of the recorded data. Notably, the proposed method is able to enhance the average signal without knowledge of the signal strength of individual shots. Simple metrics for the probe pulses are introduced, such as an effective pulse duration applicable to SASE pulses characterised in the time domain using e.g. an X-band transverse cavity. The approach is evaluated using simulated and real pulse data in the context of ultrafast electron dynamics in a molecule. Utilising H 2 as a model system, we demonstrate the efficacy of the method theoretically, successfully enhancing the predicted nonresonant ultrafast x-ray scattering signal associated with electron dynamics. The method presented is broadly applicable and offers a general strategy for enhancing time-dependent observables at XFELs.

data analysis↗

Three-flavor collective neutrino oscillation simulations on a qubit quantum annealer

Neutrinos are unique among elementary particles in that their flavor-compositions oscillate over time. In extreme environments such as core-collapse supernovae, neutron-star mergers, and the early Universe, neutrinos are dense enough that their self-interactions significantly affect, if not dominate, these oscillations. This has implications for several phenomena within these environments, particularly nucleosynthesis. Simulations of these self-interactions have traditionally approximated neutrinos as having two flavors instead of the physical three. In order to develop techniques for characterizing the resulting quantum entanglement, I present the results of simulations of neutrino-neutrino interactions that include all three physical neutrino flavors and were performed on D-Wave Inc.’s Advantage 5000+ qubit quantum annealer. These results are checked against those from exact classical simulations, which are also used to compare the neutrino-neutrino interactions to neutrino-antineutrino and interactions between Majorana neutrinos, which are their own antiparticles. The D-Wave Advantage annealer is shown to be able to reproduce time evolution with the precision of a classical machine for small numbers of neutrinos and to do so without the Trotter errors present in most simulations of dynamics on quantum devices. Furthermore, it suffers from poor scaling in qubit-count with the number of neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Solver Using Singular Value Decomposition for Computational Fluid Dynamics

Numerical solutions for fluid flow problems are challenging and have been focus of Computational Fluid Dynamics (CFD) research for past several decades. The advent of quantum computing promises exponential speedup in comparison to existing classical methods and alleviate computational constraints posed by CFD problems. Although solutions for most problems of interest in fluid dynamics using quantum computing are distant, recent advances in algorithms, software and hardware provide a path towards realizing this goal. Quantum linear solver algorithms (QLSA) such as Harrow–Hassidim–Lloyd (HHL) and Variational Quantum Linear Solver (VQLS) have been successfully implemented to solve for canonical problems such as Hele-Shaw flow. However, these algorithms still suffer to scale and address problems with ill-conditioned Jacobians. In the current paper, we alleviate these restrictions with a new quantum solver based on Singular Value Decomposition (SVD) and simulate flow past a 2D cylinder. The fidelity of the SVD based quantum solver in predicting the flow past 2D cylinder is computed along with an assessment of errors. Classical and quantum solutions for the flow are compared for different resolutions. Finally, we discuss variation in the solutions based on number of shots used.

Gottiparthi, Kalyan [ORNL] (ORCID:0000000213540255↗

Quantum surrogate models for uncertainty quantification

Surrogate models are a critical ingredient to computation-based design and validation of many DOE mission-relevant physical systems. When first-principles computation of properties of a physical systems becomes pro hibitive, surrogate models are the only path towards achieving tasks such as uncertainty quantification (UQ), exploration of design space, and validation of design choices. In this project we have developed and demonstrated a new surro gate modeling paradigm for complex models that is data-driven, non-intrusive, and has the potential to be versatile and equipped with performance guaran tees. This combination of features is absent in existing surrogate modeling tools. The framework we have developed in this project exploits a quantum-classical correspondence to establish a quantum system that mimics the dynamics of the classical Hamiltonian system from which data in the form of temporal snapshots is provided. Since quantum dynamics propagates distributions over observables, the framework is naturally suited to propagation of epistemic uncertainties in the form of distributions over initial state and parametric uncertainties. In this project, we take the first step in establishing this novel framework by deriving a quantization and de-quantization procedure, demonstrating the accuracy of the quantum surrogate models these define using two model systems, and defining the next steps in maturing the framework towards a tool applicable to Sandia mission-relevant problems.

97 MATHEMATICS AND COMPUTING↗

Applying Quantum Computing to Simulate Power System Dynamics

Power system dynamics are generally modeled by high dimensional nonlinear differential-algebraic equations due to a large number of generators, loads, and transmission lines. Thus, its computational complexity grows exponentially with the system size. This paper demonstrates the potential use of quantum computing algorithms to model the power system dynamics. Leveraging a symbolic programming framework, we equivalently convert the power system dynamics’ differential algebraic equations (DAEs) into ordinary differential equations (ODEs), where the data of the state vector can be encoded into quantum computers via amplitude encoding. The system's nonlinearity is captured by Taylor polynomial expansion, the quantum state tensor, and Hamiltonian simulation, whereas state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can simulate the dynamics of the power system with high accuracy, whereas its complexity is polynomial in the logarithm of the system dimension. Our work also illustrates the use of scientific machine learning tools for implementing scientific computing concepts, e.g., Taylor expansion, DAEs/ODEs transform, and quantum computing solver, in the field of power engineering.

Tran, Huynh↗

Krylov complexity in mixed phase space

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

chaos & nonlinear dynamics↗

IBM Open Science Price - SWAP Gate Challenge

On November 30, 2020, IBM announced an open science price competition aimed at improving the fidelity of a SWAP gate on their 7 qubit Quantum system named ”Casablanca”. The details of the SWAP gate problem was outlined in a Jupyter notebook, which utilizes the IBM developed open source package Qiskit. The participants were allowed to work in teams with up to five people and restricted to only use open source software in solving the problem. Their results were due on April 16, 2021. Our approach to solving the SWAP gate challenge is based on the quantum optimal control techniques implemented in our open source packages Juqbox.jl and Quandary. The results from the optimal control approach is only as good as the accuracy in the description of the quantum system dynamics, characterized by a Hamiltonian model that was provided by IBM. Based on the calibrated control pulses IBM provide for the standard gate set, we developed a reverse engineering approach to calibrate our computational model. Techniques were developed to translate between Qiskit’s pulse representation and the B-spline formulation used in Juqbox and Quandary. The fidelity of the optimized pulse sequences were estimated using Qiskit’s randomized bench marking techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Applied-field magnetic structure and spectroscopy shifts of the effective spin-$\frac{1}{2}$ 𝑋⁢𝑌-like magnet Li 2 ⁢CoCl 4

Insulators containing chains of magnetic transition metal cations provide platforms for probing spin-$\frac{1}{2}$ dynamics and quantum critical behavior. Li 2 ⁢ CoCl 4 contains edge-sharing CoCl 6 octahedra that form chains along the crystallographic 𝑐 axis and orders antiferromagnetically at zero field, but questions remain about its applied-field magnetic structure and the Co 2+ spin state. Here, in this study, we show with neutron diffraction on a polycrystalline sample how the antialigned chains of cobalt moments begin to transition to a ferromagnetic state above 1.6 T. Further, using magnetic resonance absorption measurements and noninteracting spin models, we reveal the strongly anisotropic nature of the Co 2+ ion's 𝑋⁢𝑌-like magnetic behavior (𝑔 ∥ = 2.77 and 𝑔 ⊥ = 5.23) and its 𝐽 = $\frac{1}{2}$ ground state. We therefore supply the magnetic structures and anisotropic description needed to explore the dynamics of the field-driven magnetic phases, laying the foundation for further experimental and theoretical studies.

Riedel, Zachary W. [Univ. of Illinois at Urbana-Ch↗

Redirection and reshaping of intense extreme-ultraviolet radiation

The goal to control short-wavelength radiation for the investigation and manipulation of ultrafast dynamics in quantum systems coevolves with the growing availability of extreme-ultraviolet (XUV) and x-ray sources from high-harmonic generation and free-electron lasers. Here, we present an XUV spatio-spectral phase modulator based on an intense XUV laser beam propagating through an optically thick resonant target, introducing dispersion profile variations around the resonance both perpendicular to and along the laser propagation direction. The resulting dipole radiation gets spectrally reshaped and becomes more divergent as compared to the original beam in the far field. As an experimental demonstration, the intense-XUV–induced double-peak off-axis structure in the far-field spectrum obtained at the Free-Electron Laser in Hamburg (FLASH) shows indications of the underlying XUV-driven Rabi dynamics and resonant pulse propagation effects. The presented work highlights a ubiquitous phenomenon occurring when an intense laser beam passes through a resonant medium.

He, Yu [Max-Planck-Institut für Kernphysik, Heidel↗

Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet–triplet spin–orbit couplings quite well for CH 2 S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

Beam energy dependence of net- Λ fluctuations measured by the STAR experiment at the BNL Relativistic Heavy Ion Collider

The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the quantum chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. Net-protons and net-kaons have been used as proxies for the net-baryon number and net-strangeness, respectively. We report the measurement of efficiency- and centrality-bin width-corrected cumulant ratios (C 2 /C 1 , C 3 /C 2 ) of net- Λ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality, and rapidity. The results are for Au + Au collisions at five beam energies ( √s NN = 19.6 , 27, 39, 62.4, and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to ultrarelativistic quantum molecular dynamics (UrQMD) and hadron resonance gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. UrQMD describes the measured net- Λ C 1 and C 3 at 200 GeV reasonably well but deviates from C 2 , and the deviation increases as a function of collision energy. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon, and net-charge data, indicate Λ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing with temperatures obtained from net-proton fluctuations. The net- Λ cumulants show a weak but finite dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗