Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum computing for simulation”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Automating Rabi & Ramsey Measurements via ML

As quantum computers scale up, the manual process of qubit tune-up becomes increasingly impractical due to its time-consuming and repetitive nature. While existing research has explored some automation techniques, many models remain underutilized for this purpose. This research aims to answer the question: is qubit tune-up able to be automated using the Long Short-Term Memory (LSTM) model? For the purposes of this project, only the rabi and ramsey measurement cycle was automated. These measurements are used to fine-tune a rough qubit frequency by repeating them until the optimal qubit frequency is obtained. The LSTM model uses the qubit frequency at one time step to forecast the qubit frequency at the next time step. A rabi-ramsey simulation was made to fabricate a dataset to train and test the LSTM model. As the model was trained, the error of the model decreased. Although there wasn't enough training data to generate perfect predictions, this shows it is possible to utilize forecasting models in automating the tune-up process.

Roberts, Rachel↗
DOE OSTI2019

Integration and Performance of New Technologies in the CMS Simulation [Slides]

CMS upgrade simulation for HL-LHC may be 2–3×slower than Run 2. CMS integration tests of GeantV met all goals: 1. Co-development ensured compatible threading models & interfaces, 2. Similar speedup measured in full experiment software framework, and 3. Efficient path to integration was established. GeantV prototype is a useful demonstrator. The next step is a GeantX project for HPCs w/ GPUs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improved C 1 shape functions for simplex meshes

Using continuously-differentiable (C 1 ) interpolation methods previously developed for computer visualization and other applications, a new class of shape functions is constructed and applied to electrostatic (ES) Particle-In-Cell (PIC) simulations on simplex (triangle or tetrahedral) meshes. The resulting shape functions are non-negative and partition unity,and lead to a reduction of shot noise and aliasing associated with particle interpolation to and from the mesh. For an energy-conserving implementation as chosen, C 1 interpolation of the potential leads to a continuous (C 0 ) electric field, with associated improved fidelity of particle motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantum coherence confined

Quantum coherence refers to the ability of a quantum state to maintain its entanglement and superposition in the face of interactions and the effects of thermalization. It is indispensable for quantum technologies — be they in sensing, transduction or computing — permitting them, in principle, to outperform their classical counterparts. To realize the promise of these technologies, one key challenge is to identify quantum states within systems that are quantum coherent. Now, writing in Nature Physics, Wen Lin Tan and collaborators report that they have used a system of trapped rare-earth-ion arrays to demonstrate a class of disorder-free quantum states that have enhanced coherence and experience anomalous patterns of thermalization1. Here, the pivotal notion behind their observations is the confinement of fundamental excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Role of heating and turbulent heat transfer in HE burning and detonation

Enhanced heat transfer model has been developed due to the turbulence in the cracks and pores in the burning fractured HE. This turbulence then serves in the initialization of the turbulence in the product gas. Enhanced heat transfer is also modeled in the turbulent product gas. Numerical simulations of the onset of HE burning and transition to detonation have been performed with and without the turbulent heat transfer model. Computational results indicate that onset of chemical reactions and transition to detonation is induced by the shockwave formed by the pressure rise caused by the heating. In other words, processes similar to shock-to-detonation transition (SDT) could become the driving mechanism for HE explosion and detonation in the case of heating rather than the temperature alone.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Surface Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer on a Sharp Cone at Angle of Attack

High-fidelity simulations are performed to characterize the turbulence-induced wall pressure fluctuations on a sharp cone at a 5.5° angle-of-attack in a Mach 8 flow. Wall-resolved large-eddy simulation (LES) and wall-modeled large-eddy simulation (WMLES) results are compared to measurements at several locations on the cone body, where comparisons with high-frequency PCB sensors are good, while comparison with low-frequency kulite sensors varies depending on the location. For a given streamwise location, simulation results show significant azimuthal variation in the wall pressure fluctuation statistics. Comparisons between LES and WMLES results indicate that WMLES can accurately reproduce the low-frequency component of the autospectra quite well, while slight sensitivity is apparent in the coherence. A grid sensitivity study for WMLES further demonstrated sensitivity in the wall pressure fluctuation coherence, but not in the autospectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A methodology to generate crystal-based molecular structures for atomistic simulations

Abstract We propose a systematic method to construct crystal-based molecular structures often needed as input for computational chemistry studies. These structures include crystal ‘slabs’ with periodic boundary conditions (PBCs) and non-periodic solids such as Wulff structures. We also introduce a method to build crystal slabs with orthogonal PBC vectors. These methods are integrated into our code, Los Alamos Crystal Cut ( LCC ), which is open source and thus fully available to the community. Examples showing the use of these methods are given throughout the manuscript.

42 ENGINEERING↗
DOE OSTI2020

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

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

Velocity-Space Hybridization of Direct Simulation Monte Carlo and a Quasi-Particle Boltzmann Solver

This paper presents a new method for modeling rarefied gas flows based on hybridization of direct simulation Monte Carlo (DSMC) and discrete velocity method (DVM)-based quasi-particle representations of the velocity distribution function. It is aimed at improving the resolution of the tails of the distribution function (compared with DSMC) and computational efficiency (compared with DVM). Details of the method, such as the collision algorithm and the particle merging scheme, are discussed. The hybrid approach is applied to the study of noise in a Maxwellian distribution, computation of electron-impact ionization rate coefficient, as well as numerical simulation of a supersonic Couette flow. The hybrid-based solver is compared with pure DSMC and DVM approaches in terms of accuracy, computational speed, and memory use. It is shown that such a hybrid approach can provide a lower computational cost than a pure DVM approach, while being able to retain accuracy in modeling high-velocity tails of the distribution function. For problems where trace species have a significant impact on the flow physics, the proposed method is shown to be capable of providing better computational efficiency and accuracy compared with standard fixed-weight DSMC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

Progress report on analytic and numerical studies of x-ray-induced impulse

In a variety of high-energy-density (HED) systems, x-rays of a given energy are used to generate shockwaves, bulk motion, and impulse in materials. As radiation energy is deposited within the material, the material is heated. The number of photons of a certain wavelength absorbed is determined by the spectral intensity of the radiation and the material’s opacity evaluated at that wavelength. This heating results in a pressure increase dictated by the material’s equation of state. Depending on the intensity of the absorbed radiation, it may also cause the material to change phase into a liquid, gas, or plasma. The increased pressure drives the heated surface layer to blow off, imparting impulse and sending a compressive wave into the bulk of the material. Additionally, the compression wave interacts with the solid boundary of the material, resulting in a tensile wave that may cause the material to spall. The impulse generated by the deposition of x-ray energy within the sample can be modeled using purely analytical methods, e.g. the Bethe, Bade, Averell, and Yost (BBAY) model. However, the blow-off process is rather complicated, and proper modeling efforts must account for material ejected by spallation, vaporization, jetting, and plasma ablation. For this reason, analytical models have an unclosed term describing the final energy of the blown-off material Ef(z). Prior modeling efforts have arbitrarily fixed this at some value or modeled it using a limiting set of thermodynamic assumptions. The work we are currently performing uses validated simulations using sophisticated photon transport, equation of state, and strength models/data to provide a fit for Ef(z) that is useful for predictive calculation of impulse. We will apply our methodology and show results for different materials and x-ray sources. This work is particularly useful for the design of experiments studying x-ray impulse generation. The present report is outlined as follows. Section 1.2 describes a series of HED experiments investigating x-ray-generated impulse in materials. Section 1.3 describes the computational approach we employ in this study, and presents validation results against the aforementioned experiments. Section 1.4 introduces analytical models for impulse generation, as well as our method for utilizing impulse from simulations to close the models. We discuss the concept of impulse-spectrum sensitivity, it’s application to uncertainty quantification, and derive a very useful analytical expression for it in 1.5. We then summarize recent progress in this project and discuss future work in section 1.6.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Haldane topological spin-1 chains in a planar metal-organic framework

Haldane topological materials contain unique antiferromagnetic chains with symmetry-protected energy gaps. Such materials have potential applications in spintronics and future quantum computers. Haldane topological solids typically consist of spin-1 chains embedded in extended three-dimensional (3D) crystal structures. Here, we demonstrate that [Ni(μ-4,4'-bipyridine)(μ-oxalate)] n (NiBO) instead adopts a two-dimensional (2D) metal-organic framework (MOF) structure of Ni 2+ spin-1 chains weakly linked by 4,4'-bipyridine. NiBO exhibits Haldane topological properties with a gap between the singlet ground state and the triplet excited state. The latter is split by weak axial and rhombic anisotropies. Several experimental probes, including single-crystal X-ray diffraction, variable-temperature powder neutron diffraction (VT-PND), VT inelastic neutron scattering (VT-INS), DC susceptibility and specific heat measurements, high-field electron spin resonance, and unbiased quantum Monte Carlo simulations, provide a detailed, comprehensive characterization of NiBO. Vibrational (also known as phonon) properties of NiBO have been probed by INS and density-functional theory (DFT) calculations, indicating the absence of phonons near magnetic excitations in NiBO, suppressing spin-phonon coupling. The work here demonstrates that NiBO is indeed a rare 2D-MOF Haldane topological material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Freedom of the mixer rotation axis improves performance in the quantum approximate optimization algorithm

Variational quantum algorithms such as the quantum approximate optimization algorithm (QAOA) are particularly attractive candidates for implementation on near-term quantum processors. As hardware realities such as error and qubit connectivity will constrain achievable circuit depth in the near future, new ways to achieve high performance at low depth are of great interest. In this work, we present a modification to QAOA that adds additional variational parameters in the form of freedom of the rotation axis in the XY plane of the mixer Hamiltonian. Here, via numerical simulation, we show that this leads to a drastic performance improvement over standard QAOA at finding solutions to the MaxCut problem on graphs of up to seven qubits. Furthermore, we explore the Z-phase error mitigation properties of our modified Ansatz, its performance under a realistic error model for a neutral atom quantum processor, and the class of problems it can solve in a single round.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Atomic structure of an FeCrMoCBY metallic glass revealed by high energy x-ray diffraction

Amorphous bulk metallic glasses with the composition Fe 48 Cr 15 Mo 14 C 15 B 6 Y 2 have been of interest due to their special mechanical and electronic properties, including corrosion resistance, high yield-strength, large elasticity, catalytic performance, and soft ferromagnetism. Here, we apply a reverse Monte Carlo technique to unravel the atomic structure of these glasses. The pair-distribution functions for various atomic pairs are computed based on the high-energy x-ray diffraction data we have taken from an amorphous sample. Monte Carlo cycles are used to move the atomic positions until the model reproduces the experimental pair-distribution function. The resulting fitted model is consistent with our ab initio simulations of the metallic glass. Our study contributes to the understanding of functional properties of Fe-based bulk metallic glasses driven by disorder effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Leveraging Hamiltonian simulation techniques to compile operations on bosonic devices

Circuit quantum electrodynamics enables the combined use of qubits and oscillator modes. Despite a variety of available gate sets, many hybrid qubit-boson (i.e. qubit-oscillator) operations are realizable only through optimal control theory, which is oftentimes intractable and uninterpretable. We introduce an analytic approach with rigorously proven error bounds for realizing specific classes of operations via two matrix product formulas commonly used in Hamiltonian simulation, the Lie–Trotter–Suzuki and Baker–Campbell–Hausdorff product formulas. We show how this technique can be used to realize a number of operations of interest, including polynomials of annihilation and creation operators, namely (a) p (a † ) q for integer p, q. We show examples of this paradigm including obtaining universal control within a subspace of the entire Fock space of an oscillator, state preparation of a fixed photon number in the cavity, simulation of the Jaynes–Cummings Hamiltonian, and simulation of the Hong-Ou-Mandel effect. This work demonstrates how techniques from Hamiltonian simulation can be applied to better control hybrid qubit-boson devices.

bosonic qubits↗
DOE OSTI2023

Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry

Due to intense interest in the potential applications of quantum computing, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics. Through numerical studies of quantum state preparation and empirical complexity analysis (including the error scaling) of classical heuristics, in both ab initio and model Hamiltonian settings, we conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for ground-state quantum chemistry through polynomial speedups, it may be prudent to assume exponential speedups are not generically available for this problem.

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

How Aqueous Solvation Impacts the Frequencies and Intensities of Infrared Absorption Bands in Flavin: The Quest for a Suitable Solvent Model

FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin’s isoalloxazine ring in the diagnostic 1450–1750 cm−1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution.

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

Lossy compression of statistical data using quantum annealer

Abstract We present a new lossy compression algorithm for statistical floating-point data through a representation learning with binary variables. The algorithm finds a set of basis vectors and their binary coefficients that precisely reconstruct the original data. The optimization for the basis vectors is performed classically, while binary coefficients are retrieved through both simulated and quantum annealing for comparison. A bias correction procedure is also presented to estimate and eliminate the error and bias introduced from the inexact reconstruction of the lossy compression for statistical data analyses. The compression algorithm is demonstrated on two different datasets of lattice quantum chromodynamics simulations. The results obtained using simulated annealing show 3–3.5 times better compression performance than the algorithm based on neural-network autoencoder. Calculations using quantum annealing also show promising results, but performance is limited by the integrated control error of the quantum processing unit, which yields large uncertainties in the biases and coupling parameters. Hardware comparison is further studied between the previous generation D-Wave 2000Q and the current D-Wave Advantage system. Our study shows that the Advantage system is more likely to obtain low-energy solutions for the problems than the 2000Q.

97 MATHEMATICS AND COMPUTING↗