Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed quantum computing”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Molecular dynamics investigation of hexagonal boron nitride sputtering and sputtered particle characteristics

The sputtering of hexagonal boron nitride (h-BN) by impacts of energetic xenon ions is investigated using a molecular dynamics (MD) model. The model is implemented within an open-source MD framework that utilizes graphics processing units to accelerate its calculations, allowing the sputtering process to be studied in much greater detail than has been feasible in the past. Integrated sputter yields are computed over a range of ion energies from 20 eV to 300 eV, and incidence angles from 0° to 75°. Sputtering of boron is shown to occur at energies as low as 40 eV at normal incidence, and sputtering of nitrogen at as low as 30 eV at normal incidence, suggesting a threshold energy between 20 eV and 40 eV. The sputter yields at 0° incidence are compared to existing experimental data and are shown to agree well over the range of ion energies investigated. The semi-empirical Bohdansky curve and an empirical exponential function are fit to the data at normal incidence, and the threshold energy for sputtering is calculated from the Bohdansky curve fit as 35 ± 2 eV. These results are shown to compare well with experimental observations that the threshold energy lies between 20 eV and 40 eV. It is demonstrated that h-BN sputters predominantly as atomic boron and diatomic nitrogen, and the velocity distribution function (VDF) of sputtered boron atoms is investigated. The calculated VDFs are found to reproduce the Sigmund-Thompson distribution predicted by Sigmund's linear cascade theory of sputtering. The average surface binding energy computed from Sigmund-Thompson curve fits is found to be 4.5 eV for ion energies of 100 eV and greater. This compares well to the value of 4.8 eV determined from independent experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Low-overhead transversal fault tolerance for universal quantum computation

Fast, reliable logical operations are essential for realizing useful quantum computers. By redundantly encoding logical qubits into many physical qubits and using syndrome measurements to detect and correct errors, we can achieve low logical error rates. However, for many practical quantum error correction codes such as the surface code, owing to syndrome measurement errors, standard constructions require multiple extraction rounds—of the order of the code distance d—for fault-tolerant computation, particularly considering fault-tolerant state preparation. Here we show that logical operations can be performed fault-tolerantly with only a constant number of extraction rounds for a broad class of quantum error correction codes, including the surface code with magic state inputs and feedforward, to achieve ‘transversal algorithmic fault tolerance’. Through the combination of transversal operations7 and new strategies for correlated decoding, despite only having access to partial syndrome information, we prove that the deviation from the ideal logical measurement distribution can be made exponentially small in the distance, even if the instantaneous quantum state cannot be made close to a logical codeword because of measurement errors. We supplement this proof with circuit-level simulations in a range of relevant settings, demonstrating the fault tolerance and competitive performance of our approach. Furthermore, our work sheds new light on the theory of quantum fault tolerance and has the potential to reduce the space–time cost of practical fault-tolerant quantum computation by over an order of magnitude.

Zhou, Hengyun [QuEra Computing, Boston, MA (United↗
DOE OSTI2016

Hierarchical optimization for neutron scattering problems

We present a scalable optimization method for neutron scattering problems that determines confidence regions of simulation parameters in lattice dynamics models used to fit neutron scattering data for crystalline solids. The method uses physics-based hierarchical dimension reduction in both the computational simulation domain and the parameter space. We demonstrate for silicon that after a few iterations the method converges to parameters values (interatomic force-constants) computed with density functional theory simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Validating quantum-classical programming models with tensor network simulations

The exploration of hybrid quantum-classical algorithms and programming models on noisy near-term quantum hardware has begun. As hybrid programs scale towards classical intractability, validation and benchmarking are critical to understanding the utility of the hybrid computational model. In this paper, we demonstrate a newly developed quantum circuit simulator based on tensor network theory that enables intermediate-scale verification and validation of hybrid quantum-classical computing frameworks and programming models. We present our tensor-network quantum virtual machine (TNQVM) simulator which stores a multi-qubit wavefunction in a compressed (factorized) form as a matrix product state, thus enabling single-node simulations of larger qubit registers, as compared to brute-force statevector simulators. Our simulator is designed to be extensible in both the tensor network form and the classical hardware used to run the simulation (multicore, GPU, distributed). The extensibility of the TNQVM simulator with respect to the simulation hardware type is achieved via a pluggable interface for different numerical backends (e.g., ITensor and ExaTENSOR numerical libraries). We demonstrate the utility of our TNQVM quantum circuit simulator through the verification of randomized quantum circuits and the variational quantum eigensolver algorithm, both expressed within the eXtreme-scale ACCelerator (XACC) programming model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Binding complexity and multiparty entanglement

We introduce “binding complexity”, a new notion of circuit complexity which quantifies the difficulty of distributing entanglement among multiple parties, each consisting of many local degrees of freedom. We define binding complexity of a given state as the minimal number of quantum gates that must act between parties to prepare it. To illustrate the new notion we compute it in a toy model for a scalar field theory, using certain multiparty entangled states which are analogous to configurations that are known in AdS/CFT to correspond to multiboundary wormholes. Pursuing this analogy, we show that our states can be prepared by the Euclidean path integral in (0 + 1)-dimensional quantum mechanics on graphs with wormhole-like structure. We compute the binding complexity of our states by adapting the Euler-Arnold approach to Nielsen’s geometrization of gate counting, and find a scaling with entropy that resembles a result for the interior volume of holographic multiboundary wormholes. We also compute the binding complexity of general coherent states in perturbation theory, and show that for “double-trace deformations” of the Hamiltonian the effects resemble expansion of a wormhole interior in holographic theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

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

Precision test of statistical dynamics with state-to-state ultracold chemistry

Chemical reactions represent a class of quantum problems that challenge both the current theoretical understanding and computational capabilities1. Reactions that occur at ultralow temperatures provide an ideal testing ground for quantum chemistry and scattering theories, because they can be experimentally studied with unprecedented control, yet display dynamics that are highly complex. Here we report the full product state distribution for the reaction 2KRb → K 2 + Rb 2 . Ultracold preparation of the reactants allows us complete control over their initial quantum degrees of freedom, whereas state-resolved, coincident detection of both products enables the probability of scattering into each of the 57 allowed rotational state-pairs to be measured. Our results show an overall agreement with a state-counting model based on statistical theory, but also reveal several deviating state-pairs. In particular, we observe a strong suppression of population in the state-pair closest to the exoergicity limit as a result of the long-range potential inhibiting the escape of products. Finally, the completeness of our measurements provides a benchmark for quantum dynamics calculations beyond the current state of the art.

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

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Parallel Fokker–Planck-DSMC algorithm for rarefied gas flow simulation in complex domains at all Knudsen numbers

A major challenge for the conventional Direct Simulation Monte Carlo (DSMC) technique lies in the fact that its computational cost becomes prohibitive in the near continuum regime, where the Knudsen number (Kn)—characterizing the degree of rarefaction—becomes small. In contrast, the Fokker–Planck (FP) based particle Monte Carlo scheme allows for computationally efficient simulations of rarefied gas flows in the low and intermediate Kn regime. The Fokker–Planck collision operator—instead of performing binary collisions employed by the DSMC method—integrates continuous stochastic processes for the phase space evolution in time. This allows for time step and grid cell sizes larger than the respective collisional scales required by DSMC. Dynamically switching between the FP and the DSMC collision operators in each computational cell is the basis of the combined FP-DSMC method, which has been proven successful in simulating flows covering the whole Kn range. Until recently, this algorithm had only been applied to two-dimensional test cases. In this contribution, we present the first general purpose implementation of the combined FP-DSMC method. Utilizing both shared- and distributed-memory parallelization, this implementation provides the capability for simulations involving many particles and complex geometries by exploiting state of the art computer cluster technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Successive phase transitions and kink solutions in Φ⁸, Φ¹⁰, and Φ¹² field theories

We obtain exact solutions for kinks in Φ⁸, Φ¹⁰, and Φ¹² field theories with degenerate minima, which can describe a second-order phase transition followed by a first-order one, a succession of two first-order phase transitions and a second-order phase transition followed by two first-order phase transitions, respectively. Such phase transitions are known to occur in ferroelastic and ferroelectric crystals and in meson physics. In particular, we find that the higher-order field theories have kink solutions with algebraically-decaying tails and also asymmetric cases with mixed exponential-algebraic tail decay, unlike the lower-order Φ⁴ and Φ⁶ theories. Additionally, we construct distinct kinks with equal energies in all three field theories considered, and we show the co-existence of up to three distinct kinks (for a Φ¹² potential with six degenerate minima). We also summarize phonon dispersion relations for these systems, showing that the higher-order field theories have specific cases in which only nonlinear phonons are allowed. For the Φ¹⁰ field theory, which is a quasi-exactly solvable (QES) model akin to Φ⁶, we are also able to obtain three analytical solutions for the classical free energy as well as the probability distribution function in the thermodynamic limit.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Dimension-independent likelihood-informed MCMC

Many Bayesian inference problems require exploring the posterior distribution of high-dimensional parameters that represent the discretization of an underlying function. This work introduces a family of Markov chain Monte Carlo (MCMC) samplers that can adapt to the particular structure of a posterior distribution over functions. Two distinct lines of research intersect in the methods developed here. First, we introduce a general class of operator-weighted proposal distributions that are well defined on function space, such that the performance of the resulting MCMC samplers is independent of the discretization of the function. Second, by exploiting local Hessian information and any associated low-dimensional structure in the change from prior to posterior distributions, we develop an inhomogeneous discretization scheme for the Langevin stochastic differential equation that yields operator-weighted proposals adapted to the non-Gaussian structure of the posterior. The resulting dimension-independent and likelihood-informed (DILI) MCMC samplers may be useful for a large class of high-dimensional problems where the target probability measure has a density with respect to a Gaussian reference measure. Two nonlinear inverse problems are used to demonstrate the efficiency of these DILI samplers: an elliptic PDE coefficient inverse problem and path reconstruction in a conditioned diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

Cyclotron Resonant Scattering Feature Simulations I. Thermally Averaged Cyclotron Scattering Cross Sections, Mean Free Photon-Path Tables, and Electron Momentum Sampling

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01 B/B(sub crit) is less than or equal to 0.12, where B(sub crit) = 4.413 x 1013 G, and 3 keV kBT 15 keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗
DOE OSTI2022

Error rate reduction of single-qubit gates via noise-aware decomposition into native gates

In the current era of Noisy Intermediate-Scale Quantum (NISQ) technology, the practical use of quantum computers remains inhibited by our inability to aptly decouple qubits from their environment to mitigate computational errors. In this paper, we introduce an approach by which knowledge of a qubit’s initial quantum state and the standard parameters describing its decoherence can be leveraged to mitigate the noise present during the execution of a single-qubit gate. We benchmark our protocol using cloud-based access to IBM quantum processors. On ibmq_rome, we demonstrate a reduction of the single-qubit error rate by 38%, from $1.6 \times 10 ^{-3}$ to $1.0 \times 10 ^{-3}$ , provided the initial state of the input qubit is known. On ibmq_bogota, we prove that our protocol will never decrease gate fidelity, provided the system’s T 1 and T 2 times have not drifted above 100 times their assumed values. The protocol can be used to reduce quantum state preparation errors, as well as to improve the fidelity of quantum circuits for which some knowledge of the qubits’ intermediate states can be inferred. This paper presents a pathway to using information about noise levels and quantum state distributions to significantly reduce error rates associated with quantum gates via optimized decomposition into native hardware gates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Numerically exact full counting statistics of the nonequilibrium Anderson impurity model

The time-dependent full counting statistics of charge transport through an interacting quantum junction is evaluated from its generating function, controllably computed with the inchworm Monte Carlo method. Exact noninteracting results are reproduced; then, we continue to explore the effect of electron-electron interactions on the time-dependent charge cumulants, first-passage time distributions, and n-electron transfer distributions. We observe a crossover in the noise from Coulomb blockade to Kondo-dominated physics as the temperature is decreased. In addition, we uncover long-tailed spin distributions in the Kondo regime and analyze queuing behavior caused by correlations between single-electron transfer events

36 MATERIALS SCIENCE↗
DOE OSTI2015

First principles-based multiparadigm, multiscale strategy for simulating complex materials processes with applications to amorphous SiC films

Progress has recently been made in developing reactive force fields to describe chemical reactions in systems too large for quantum mechanical (QM) methods. In particular, ReaxFF, a force field with parameters that are obtained solely from fitting QM reaction data, has been used to predict structures and properties of many materials. Important applications require, however, determination of the final structures produced by such complex processes as chemical vapor deposition, atomic layer deposition, and formation of ceramic films by pyrolysis of polymers. This requires the force field to properly describe the formation of other products of the process, in addition to yielding the final structure of the material. We describe a strategy for accomplishing this and present an example of its use for forming amorphous SiC films that have a wide variety of applications. Extensive reactive molecular dynamics (MD) simulations have been carried out to simulate the pyrolysis of hydridopolycarbosilane. The reaction products all agree with the experimental data. After removing the reaction products, the system is cooled down to room temperature at which it produces amorphous SiC film, for which the computed radial distribution function, x-ray diffraction pattern, and the equation of state describing the three main SiC polytypes agree with the data and with the QM calculations. Extensive MD simulations have also been carried out to compute other structural properties, as well the effective diffusivities of light gases in the amorphous SiC film.

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

Realization of polarization evolution on higher-order Poincaré sphere with metasurface

We present a simple and convenient method to yield cylindrical vector (CV) beams and realize its polarization evolution on higher-order Poincaré sphere based on inhomogeneous birefringent metasurface. By means of local polarization transformation of the metasurface, it is possible to convert a light beam with homogeneous elliptical polarization into a vector beam with any desired polarization distribution. The Stokes parameters of the output light are measured to verify our scheme, which show well agreement with the theoretical prediction. Our method may provide a convenient way to generate CV beams, which is expected to have potential applications in encoding information and quantum computation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗