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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo

We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius–Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. Here, we compare with previous simulations and experimental measurements.

08 HYDROGEN↗
NASA NTRS2003

NASA Tech Briefs, October 2003

Topics covered include: Cryogenic Temperature-Gradient Foam/Substrate Tensile Tester; Flight Test of an Intelligent Flight-Control System; Slat Heater Boxes for Thermal Vacuum Testing; System for Testing Thermal Insulation of Pipes; Electrical-Impedance-Based Ice-Thickness Gauges; Simulation System for Training in Laparoscopic Surgery; Flasher Powered by Photovoltaic Cells and Ultracapacitors; Improved Autoassociative Neural Networks; Toroidal-Core Microinductors Biased by Permanent Magnets; Using Correlated Photons to Suppress Background Noise; Atmospheric-Fade-Tolerant Tracking and Pointing in Wireless Optical Communication; Curved Focal-Plane Arrays Using Back-Illuminated High-Purity Photodetectors; Software for Displaying Data from Planetary Rovers; Software for Refining or Coarsening Computational Grids; Software for Diagnosis of Multiple Coordinated Spacecraft; Software Helps Retrieve Information Relevant to the User; Software for Simulating a Complex Robot; Software for Planning Scientific Activities on Mars; Software for Training in Pre-College Mathematics; Switching and Rectification in Carbon-Nanotube Junctions; Scandia-and-Yttria-Stabilized Zirconia for Thermal Barriers; Environmentally Safer, Less Toxic Fire-Extinguishing Agents; Multiaxial Temperature- and Time-Dependent Failure Model; Cloverleaf Vibratory Microgyroscope with Integrated Post; Single-Vector Calibration of Wind-Tunnel Force Balances; Microgyroscope with Vibrating Post as Rotation Transducer; Continuous Tuning and Calibration of Vibratory Gyroscopes; Compact, Pneumatically Actuated Filter Shuttle; Improved Bearingless Switched-Reluctance Motor; Fluorescent Quantum Dots for Biological Labeling; Growing Three-Dimensional Corneal Tissue in a Bioreactor; Scanning Tunneling Optical Resonance Microscopy; The Micro-Arcsecond Metrology Testbed; Detecting Moving Targets by Use of Soliton Resonances; and Finite-Element Methods for Real-Time Simulation of Surgery.

Source record↗
DOE OSTI2024

Learning a General Model of Single Phase Flow in Complex 3D Porous Media

Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electromagnetic Field Simulation Utilizing VERTEX

The understanding of electromagnetic fields is essential for many core applications at Oak Ridge National Laboratory (ORNL). Electromagnetic fields play an important role in particle accelerators, antennas, electric motors, power grids, and are principal part of nuclear fusion. We develop a numerical solver called able to treat the Maxwell system utilizing ORNL's VERTEX framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Entanglement-based quantum digital signatures over a deployed campus network

The quantum digital signature protocol offers a replacement for most aspects of public-key digital signatures ubiquitous in today’s digital world. A major advantage of a quantum-digital-signatures protocol is that it can have information-theoretic security, whereas public-key cryptography cannot. Here we demonstrate and characterize hardware to implement entanglement-based quantum digital signatures over our campus network. Over 25 hours, we collect measurements on our campus network, where we measure sufficiently low quantum bit error rates (<5% in most cases) which in principle enable quantum digital signatures at over 50 km as shown through rigorous simulation accompanied by a noise model developed specifically for our implementation. These results show quantum digital signatures can be successfully employed over deployed fiber. Moreover, our reported method provides great flexibility in the number of users, but with reduced entanglement rate per user. Finally, while the current implementation of our entanglement-based approach has a low signature rate, feasible upgrades would significantly increase the signature rate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An energetics assessment of benzo[ a ]tetracene and benzo[ a ]pyrene as triplet–triplet annihilation emitters

Optical upconversion (UC) of low energy photons into high energy photons enables solar cells to harvest photons with energies below the band gap of the absorber, reducing the transmission loss. UC based on triplet–triplet annihilation (TTA) in organic chromophores can upconvert photons from sunlight, albeit with low conversion efficiency. We utilize three energy-based criteria to assess the UC potential of TTA emitters in terms of the quantum yield (QY) and the anti-Stokes shift. The energy loss in the singlet pathway of an emitter encounter complex, where a high energy photon is emitted, determines whether a chromophore may undergo TTA. The energy loss in the triplet pathway, which is the main competing process, impacts the TTA QY. The energy difference between the lowest singlet and triplet excitation states in TTA emitters sets an upper bound for the anti-Stokes shift of TTA-UC. Using the energetic criteria evaluated by time-dependent density functional theory (TDDFT) calculations, we find that benzo[a]tetracene, benzo[a]pyrene, and their derivatives are promising TTA emitters. The energetics assessment and computer simulations could be used to efficiently discover and design more candidate high-performance TTA emitters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Merged-Element Transmon

Transmon qubits are ubiquitous in the pursuit of quantum computing using superconducting circuits. However, they have some drawbacks that still need to be addressed. Most importantly, the scalability of transmons is limited by the large device footprint needed to reduce the participation of the lossy capacitive parts of the circuit. In this work, we investigate and evaluate losses in an alternative device geometry, namely, the merged-element transmon (mergemon). To this end, we replace the large external shunt capacitor of a traditional transmon with the intrinsic capacitance of a Josephson junction and achieve an approximately 100 times reduction in qubit dimensions. We report the implementation of the mergemon using a sputtered Nb–amorphous-Si–Nb trilayer film. In an experiment below 10 mK, the frequency of the readout resonator, capacitively coupled to the mergemon, exhibits a qubit-state-dependent shift in the low-power regime. The device also demonstrates the single-photon and multiphoton transitions that represent a weakly anharmonic system in two-tone spectroscopy. The transition spectra are explained well with master-equation simulations. Here, a participation-ratio analysis identifies the dielectric loss of the amorphous-Si tunnel barrier and its interfaces as the dominant source for qubit relaxation. We expect the mergemon to achieve high coherence for relatively small device dimensions when implemented with a low-loss, epitaxially grown, and lattice-matched trilayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
DOE OSTI2025

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗
DOE OSTI2018

Developing a Computational Chemistry Framework for the Exascale Era

Within computational chemistry the NWChem package has arguably been the de facto standard for running high-accuracy numerical simulations on the most powerful supercomputers. In order to better address the challenges presented by emerging exascale architectures, the decision has been made to rewrite NWChem. Design of the resulting package, NWChemEx, has been driven by exascale computing; however, significant additional design considerations have arisen from the team's involvement with the Molecular Sciences Software Institute (MolSSI). MolSSI is a National Science Foundation initiative focused on establishing coding and data standards for the computational chemistry community. As a result, NWChemEx is built upon a general computational chemistry framework called the Simulation Development Environment (SDE) that is designed with a focus on extensibility and interoperability. The present manuscript describes the modular approach of the SDE and how it has been used to implement the self-consistent field algorithm within NWChemEx.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FY2021 Improvements to the New CTH Code Verification & Validation Test Suite

Over the past few years, the CTH multiphysics hydrocode has overhauled its software quality and testing processes, implementing current best practices in software quality and building a robust V&V test suite comprised of traditional hydrocode verification problems, including ASC Tri-Lab Test Suite and Enhanced Tri-Lab Test Suite problems, as well as validation problems for some of CTH’s most frequently used equations of state, materials models, and other key capabilities. Substantial progress towards building this new test suite was made in FY19 and FY20. In FY21, the test suite has been expanded to include verification and validation tests of the Steinberg-Guinan-Lund (ST) viscoplastic model and the Johnson Cook (JFRAC) fracture model. Additionally, two new verification tests were added, covering hydrodynamics and high explosive (HE) modeling capabilities: the Kidder Gaussian density problem and the Escape of HE Products (EHEP) problem from the Tri-Lab Test Suite. This report discusses each of these test problems in detail. Verification test results are compared to analytic solutions. Validation test results are compared to experimental data. Wherever possible, convergence or mesh refinement studies are included. Additionally, while implementing the Kidder verification problem, a bug was identified that affects the use of tables to initialize pressure or density in 1D or 2D calculations. A brief discussion of the bug and its fix is included. CTH demonstrates good performance overall on the new test suite problems. Simulation results showed good agreement with analytic solutions for the Kidder problem, with convergence rates ranging between 1.8 and sub-linear, and relatively good agreement for the EHEP problem, though convergence rates for pressure and density were nearly 0. The ST and JFRAC strain rate loading verification tests show good agreement with analytic solutions. Likewise, CTH simulation results show good agreement with experimental validation data, including Taylor rod impact testing, for the materials tested. Future V&V work will focus on adding 2D and 3D versions of existing verification tests as well as adding validation tests of other frequently used capabilities such as other fracture models.

42 ENGINEERING↗
DOE OSTI2021

Tuning spin–orbit coupling in (6,5) single-walled carbon nanotube doped with sp 3 defects

Single-walled carbon nanotubes (SWCNTs) containing sp 3 defects are a promising class of optoelectronic materials with bright photoluminescence and demonstrated single-photon emission. Using density functional theory simulations, complemented by measurements, we investigate the electronic structure of a series of quantum defects attached to (6,5) SWCNT with the goal of tuning the spin–orbit coupling by introduction of a heavy atom in the defect structure. We characterize the ground state electronic and spin properties of four synthesized and three potential defects on the tube and find that all of the synthesized defects considered introduce a localized midgap defect-centered state containing a single electron, ≈0.2–0.3 eV above the valence band. The spin density is located at the sp 3 defect site with negligible spin–orbit coupling even with the presence of a Pd atom. Three additional functional groups were tested via computation to increase spin localization near the metal, thereby increasing spin–orbit coupling. We predict that only the chlorodiphosphanepalladium(II)– [Cl(PH 3 ) 2 Pd(II)–] defect results in increased spin–orbit splitting of the defect state and the conduction band associated with the pristine-like SWCNT, a measure of the spin–orbit coupling of excited state transitions. This study suggests that for unpassivated sp 3 defects in (6,5) SWCNT, forming a direct bond between a heavy atom and the sp 3 carbon allows for tuning of spin–orbit coupling.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Logic operations with active topological defects

Logic operations performed by semiconductor-based transistors are the basis of modern computing. There is considerable interest in creating autonomous materials systems endowed with the capability to make decisions. In this work, we introduce the concept of using topological defects in active matter to perform logic operations. When an extensile active stress in a nematic liquid crystal is turned on, +1/2 defects can self-propel, in analogy to electron transport under a voltage gradient. By relying on hydrodynamic simulations of active nematics, we demonstrate that patterns of activity, when combined with surfaces imparting certain orientations, can be used to control the formation and transport of +1/2 defects. We further show that asymmetric high- and low-activity patterns can be used to create effective defect gates, tunnels, and amplifiers. The proposed active systems offer the potential to perform computations and transmit information in active soft materials, including actin-, tubulin-, and cell-based systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.

Jaffe, Richard L.↗
NASA NTRSDate not supplied

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗
DOE OSTI2020

Classical Optimizers for Noisy Intermediate-Scale Quantum Devices

We present a collection of optimizers tuned for usage on Noisy Inter-mediate-Scale Quantum (NISQ) devices. Optimizers have a range of applications in quantum computing, including the Variational Quantum Eigensolver (VQE) and Quantum Approximate Optimization (QAOA) algorithms. They have further uses in calibration, hyperparameter tuning, machine learning, etc. We employ the VQE algorithm as a case study. VQE is a hybrid algorithm, with a classical minimizer step driving the next evaluation on the quantum processor. While most results to date concentrated on tuning the quantum VQE circuit, our study indicates that in the presence of quantum noise the classical minimizer step is a weak link and a careful choice combined with tuning is required for correct results. We explore state-of-the-art gradient-free optimizers capable of handling noisy, black-box, cost functions and stress-test them using a quantum circuit simulation environment with noise injection capabilities on individual gates. Our results indicate that specifically tuned optimizers are crucial to obtaining valid science results on NISQ hardware, as well as projecting forward on fault tolerant circuits.

Lavrijsen, Wim↗