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 253 records · Page 14

Properties of reactive oxygen species by quantum Monte Carlo

The electronic properties of the oxygen molecule, in its singlet and triplet states, and of many small oxygen-containing radicals and anions have important roles in different fields of chemistry, biology, and atmospheric science. Nevertheless, the electronic structure of such species is a challenge for ab initio computational approaches because of the difficulties to correctly describe the statical and dynamical correlation effects in presence of one or more unpaired electrons. Only the highest-level quantum chemical approaches can yield reliable characterizations of their molecular properties, such as binding energies, equilibrium structures, molecular vibrations, charge distribution, and polarizabilities. In this work we use the variational Monte Carlo (VMC) and the lattice regularized Monte Carlo (LRDMC) methods to investigate the equilibrium geometries and molecular properties of oxygen and oxygen reactive species. Quantum Monte Carlo methods are used in combination with the Jastrow Antisymmetrized Geminal Power (JAGP) wave function ansatz, which has been recently shown to effectively describe the statical and dynamical correlation of different molecular systems. In particular, we have studied the oxygen molecule, the superoxide anion, the nitric oxide radical and anion, the hydroxyl and hydroperoxyl radicals and their corresponding anions, and the hydrotrioxyl radical. Overall, the methodology was able to correctly describe the geometrical and electronic properties of these systems, through compact but fully-optimised basis sets and with a computational cost which scales as N{sup 3} − N{sup 4}, where N is the number of electrons. This work is therefore opening the way to the accurate study of the energetics and of the reactivity of large and complex oxygen species by first principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solving kinetic equations with adaptive mesh in phase space for rarefied gas dynamics and plasma physics (Invited)

The paper describes an Adaptive Mesh in Phase Space (AMPS) technique for solving kinetic equations with deterministic mesh-based methods. The AMPS technique allows automatic generation of adaptive Cartesian mesh in both physical and velocity spaces using a Tree-of-Trees data structure. We illustrate advantages of AMPS for simulations of rarefied gas dynamics and electron kinetics on low temperature plasmas. In particular, we consider formation of the velocity distribution functions in hypersonic flows, particle kinetics near oscillating boundaries, and electron kinetics in a radio-frequency sheath. AMPS provide substantial savings in computational cost and increased efficiency of the mesh-based kinetic solvers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sequential ensemble transform for Bayesian inverse problems

In this work, we present the Sequential Ensemble Transform (SET) method, an approach for generating approximate samples from a Bayesian posterior distribution. The method explores the posterior distribution by solving a sequence of discrete optimal transport problems to produce a series of transport plans which map prior samples to posterior samples. We prove that the sequence of Dirac mixture distributions produced by the SET method converges weakly to the true posterior as the sample size approaches infinity. Furthermore, our numerical results indicate that, when compared to standard Sequential Monte Carlo (SMC) methods, the SET approach is more robust to the choice of Markov mutation kernels and requires less computational efforts to reach a similar accuracy when used to explore complex posterior distributions. Finally, we describe adaptive schemes that allow to completely automate the use of the SET method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Point defects in semiconductors: Friends and foes for quantum technologies

Point defects in semiconductors are both a curse and a blessing in microelectronics: they enable the control of electrical conductivity through doping, yet can also act as trapping and recombination centers that degrade device performance. In quantum information science, defects play a similarly dual role. They can be harnessed as spin–photon interfaces enabling the coupling of electronic and nuclear spins to light and the creation of distributed entanglement for quantum networks or used as atomistic scale sensors for quantum sensing. At the same time, defects are a major source of decoherence for superconducting qubits, one of the leading quantum computing platforms. This article discusses how a deeper materials-level understanding of defects can guide the design of improved quantum devices for communication, sensing, and computation.

Zhu, Yizhi [Rice University, Houston, TX (United S↗

Characterization and Optimization of the Fitting of Quantum Correlation Functions

This case study presents a characterization and optimization of an application code for extracting parton distribution functions from high energy electron-proton scattering data. Profiling this application code reveals that the phase-space density computation accounts for 93% of the overall execution time for a single iteration on a single core. When executing multiple iterations in parallel on a multicore system, the application spends 78% of its overall execution time idling due to load imbalance. We address these issues by first transforming the application code from Python to C++ and then tackling the application load imbalance via a hybrid scheduling strategy that combines dynamic and static scheduling. These techniques result in a 62% reduction in CPU idle time and a 2.46x speedup in overall execution time per node. In addition, the typically enabled power-management mechanisms in supercomputers (e.g., AMD Turbo Core, Intel Turbo Boost, and RAPL) can significantly impact intra-node scalability when more than 50% of the CPU cores are used. This finding underscores the importance of understanding system interactions with power management, as they can adversely impact application performance, and highlights the necessity of intra-node scaling tests to identify performance degradation that inter-node scaling tests might otherwise overlook.

Chuang, Pi-Yueh [Virginia Tech,Dept. of Computer S↗

Continuous Variable Cluster State Generation over the Optical Spatial Mode Comb

One way quantum computing uses single qubit projective measurements performed on a cluster state (a highly entangled state of multiple qubits) in order to enact quantum gates. The model is promising due to its potential scalability; the cluster state may be produced at the beginning of the computation and operated on over time. Continuous variables (CV) offer another potential benefit in the form of deterministic entanglement generation. This determinism can lead to robust cluster states and scalable quantum computation. Recent demonstrations of CV cluster states have made great strides on the path to scalability utilizing either time or frequency multiplexing in optical parametric oscillators (OPO) both above and below threshold. The techniques relied on a combination of entangling operators and beam splitter transformations. Here we show that an analogous transformation exists for amplifiers with Gaussian inputs states operating on multiple spatial modes. By judicious selection of local oscillators (LOs), the spatial mode distribution is analogous to the optical frequency comb consisting of axial modes in an OPO cavity. We outline an experimental system that generates cluster states across the spatial frequency comb which can also scale the amount of quantum noise reduction to potentially larger than in other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neural Networks and Graph Algorithms with Next-Generation Processors

The use of graphical processors for distributed computation revolutionized the field of high performance scientific computing. As the Moore's Law era of computing draws to a close, the development of non-Von Neumann systems: neuromorphic processing units, and quantum annealers; again are redefining new territory for computational methods. While these technologies are still in their nascent stages, we discuss their potential to advance computing in two domains: machine learning, and solving constraint satisfaction problems. Each of these processors utilize fundamentally different theoretical models of computation. This raises questions about how to best use them in the design and implementation of applications. While many processors are being developed with a specific domain target, the ubiquity of spin-glass models and neural networks provides an avenue for multi-functional applications. This provides hints at the future infrastructure needed to integrate many next-generation processing units into conventional high-performance computing systems.

Hamilton, Kathleen↗

Optimization of SiV Defect Yield in Diamond Substrates

Color centers (defect complexes such as SiV) in diamond have shown potential in fields ranging from metrology, cybersecurity to quantum computation. Demonstrations in these fields have pushed the envelope of state-of-the-art operations – for example, single photon sources (SPS) making use of SiV centers in diamond for quantum key distribution have demonstrated all the requirements for SPS operation including: (1) stable operation with second correlation function <<1, (2) electrically driven single photon emission and (3) compatibility with frequenc y conversion to telecommunication frequencies. To-date, however, all these demonstrations have been on lab-scale one-off devices. The key question behind how to deterministically fabricate these devices, namely activation yield has been overlooked. For context, Si based semiconductor devices are hugely successful because we have a high activation yield for implanted dopants. This is not yet true for diamond color centers. As currently understood, the color center yield is dominated by a lack of vacancies in the immediate area of the implantation. We propose to optimize the activation yield of color center using a combination of (1) focused single ion implantation with in-situ detection to count the number of implanted Si ions and (2) localized point defect (vacancy) creation using a focused Li ion beam to improve the yield. These experiments build on the unique capabilities of the SNL nanoImplanter (nI) to produce focused ion beam with spatial resolution of < 10 nm of both Si and Li ions. This work will also leverage our world-leading single ion implantation and detection capabilities.

36 MATERIALS SCIENCE↗

Electron-impact vibrational relaxation in high-temperature nitrogen

Vibrational relaxation process of N2 molecules by electron-impact is examined for the future planetary entry environments. Multiple-quantum transitions from excited states to higher/lower states are considered for the electronic ground state of the nitrogen molecule N2 (X 1Sigma-g(+)). Vibrational excitation and deexcitation rate coefficients obtained by computational quantum chemistry are incorporated into the 'diffusion model' to evaluate the time variations of vibrational number densities of each energy state and total vibrational energy. Results show a non-Boltzmann distribution of number densities at the earlier stage of relaxation, which in turn suppresses the equilibrium process but affects little the time variation of total vibrational energy. An approximate rate equation and a corresponding relaxation time from the excited states, compatible with the system of flow conservation equations, are derived. The relaxation time from the excited states indicates the weak dependency of the initial vibrational temperature. The empirical curve-fit formula for the improved e-V relaxation time is obtained.

Lee, Jong-Hun↗

Quantitative single-particle digital autoradiography with α-particle emitters for targeted radionuclide therapy using the iQID camera

Purpose: Alpha-emitting radionuclides exhibit a potential advantage for cancer treatments because they release large amounts of ionizing energy over a few cell diameters (50–80 μm), causing localized, irreparable double-strand DNA breaks that lead to cell death. Radioimmunotherapy (RIT) approaches using monoclonal antibodies labeled with α emitters may thus inactivate targeted cells with minimal radiation damage to surrounding tissues. Tools are needed to visualize and quantify the radioactivity distribution and absorbed doses to targeted and nontargeted cells for accurate dosimetry of all treatment regimens utilizing α particles, including RIT and others (e.g., Ra-223), especially for organs and tumors with heterogeneous radionuclide distributions. The aim of this study was to evaluate and characterize a novel single-particle digital autoradiography imager, the ionizing-radiation quantum imaging detector (iQID) camera, for use in α-RIT experiments. Methods: The iQID camera is a scintillator-based radiation detection system that images and identifies charged-particle and gamma-ray/x-ray emissions spatially and temporally on an event-by-event basis. It employs CCD-CMOS cameras and high-performance computing hardware for real-time imaging and activity quantification of tissue sections, approaching cellular resolutions. In this work, the authors evaluated its characteristics for α-particle imaging, including measurements of intrinsic detector spatial resolutions and background count rates at various detector configurations and quantification of activity distributions. The technique was assessed for quantitative imaging of astatine-211 ({sup 211}At) activity distributions in cryosections of murine and canine tissue samples. Results: The highest spatial resolution was measured at ∼20 μm full width at half maximum and the α-particle background was measured at a rate as low as (2.6 ± 0.5) × 10{sup −4} cpm/cm{sup 2} (40 mm diameter detector area). Simultaneous imaging of multiple tissue sections was performed using a large-area iQID configuration (ø 11.5 cm). Estimation of the {sup 211}At activity distribution was demonstrated at mBq/μg-levels. Conclusions: Single-particle digital autoradiography of α emitters has advantages over traditional film-based autoradiographic techniques that use phosphor screens, in terms of spatial resolution, sensitivity, and activity quantification capability. The system features and characterization results presented in this study show that the iQID is a promising technology for microdosimetry, because it provides necessary information for interpreting alpha-RIT outcomes and for predicting the therapeutic efficacy of cell-targeted approaches using α emitters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗

Application of Gauss's law space-charge limited emission model in iterative particle tracking method

The particle tracking method with a so-called gun iteration for modeling the space charge is discussed in the following paper. We suggest to apply the emission model based on the Gauss's law for the calculation of the space charge limited current density distribution using considered method. Based on the presented emission model we have developed a numerical algorithm for this calculations. This approach allows us to perform accurate and low time consumpting numerical simulations for different vacuum sources with the curved emitting surfaces and also in the presence of additional physical effects such as bipolar flows and backscattered electrons. The results of the simulations of the cylindrical diode and diode with elliptical emitter with the use of axysimmetric coordinates are presented. The high efficiency and accuracy of the suggested approach are confirmed by the obtained results and comparisons with the analytical solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Defect Structure in Quantum-Cutting Yb3+-Doped CsPbCl3 Perovskites Probed by X-Ray Absorption and Atomic Pair Distribution Function Analysis

Ytterbium-doping in all-inorganic lead-halide perovskites (CsPb(Cl1-xBrx)3) generates novel properties including quantum cutting and narrow line emission, making these materials attractive spectral down-converters for solar photovoltaics. The relationship between this optical efficiency and the defect structure(s) associated with Yb3+ dopants within perovskites is not well understood. Various charge-neutral doping motifs have previously been proposed and studied computationally, including clusters involving two substitutional Yb3+ ions charge-compensated by a single local Pb2+ vacancy. Near-band-edge defect states associated with such motifs are believed to play an important mechanistic role in quantum cutting itself. Here, we report the results of X-ray absorption and X-ray total-scattering measurements on ytterbium-doped CsPbCl3. XANES shows that the dopant oxidation state is exclusively Yb3+, and a combination of Yb L3 and Pb L3 EXAFS shows that this Yb3+ substitutes exclusively at Pb2+ sites, where it adopts a pseudo-octahedral [YbCl6]3- coordination environment. Shell-by-shell fits to the data show a short Yb-Cl bond distance of 2.58 Å compared to the Pb-Cl bond distance of 2.83 Å. We confirm this finding by X-ray pair distribution function analysis, which also shows evidence of additional Pb2+ vacancy formation induced by Yb3+ doping. We evaluate whether this is the primary mechanism of charge compensation using simulated EXAFS and pair distribution function data for several computed defect structures. Together, these results resolve the local dopant structures and charge-compensation mechanisms in lanthanide-doped all-inorganic lead-halide perovskites, and, as such, significantly advance the understanding of structure-function relationships in this important class of materials.

Kluherz, Kyle (ORCID:0000000279865167)↗

Single T gate in a Clifford circuit drives transition to universal entanglement spectrum statistics

Clifford circuits are insufficient for universal quantum computation or creating t t -designs with t\ge 4 t ≥ 4 . While the entanglement entropy is not a telltale of this insufficiency, the entanglement spectrum of a time evolved random product state is: the entanglement levels are Poisson-distributed for circuits restricted to the Clifford gate-set, while the levels follow Wigner-Dyson statistics when universal gates are used. In this paper we show, using finite-size scaling analysis of different measures of level spacing statistics, that in the thermodynamic limit, inserting a single T (\pi/8) ( π / 8 ) gate in the middle of a random Clifford circuit is sufficient to alter the entanglement spectrum from a Poisson to a Wigner-Dyson distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Training a Quantum Annealing Based Restricted Boltzmann Machine on Cybersecurity Data

A restricted Boltzmann machine (RBM) is a generative model that could be used in effectively balancing a cybersecurity dataset because the synthetic data a RBM generates follows the probability distribution of the training data. RBM training can be performed using contrastive divergence (CD) and quantum annealing (QA). QA-based RBM training is fundamentally different from CD and requires samples from a quantum computer. We present a real-world application that uses a quantum computer. Specifically, we train a RBM using QA for cybersecurity applications. The D-Wave 2000Q has been used to implement QA. RBMs are trained on the ISCX data, which is a benchmark dataset for cybersecurity. For comparison, RBMs are also trained using CD. CD is a commonly used method for RBM training. Our analysis of the ISCX data shows that the dataset is imbalanced. We present two different schemes to balance the training dataset before feeding it to a classifier. The first scheme is based on the undersampling of benign instances. The imbalanced training dataset is divided into five sub-datasets that are trained separately. A majority voting is then performed to get the result. Our results show the majority vote increases the classification accuracy up from 90.24% to 95.68%, in the case of CD. For the case of QA, the classification accuracy increases from 74.14% to 80.04%. In the second scheme, a RBM is used to generate synthetic data to balance the training dataset. We show that both QA and CD-trained RBM can be used to generate useful synthetic data. Balanced training data is used to evaluate several classifiers. Among the classifiers investigated, K-Nearest Neighbor (KNN) and Neural Network (NN) perform better than other classifiers. They both show an accuracy of 93%. Our results show a proof-of-concept that a QA-based RBM can be trained on a 64-bit binary dataset. The illustrative example suggests the possibility to migrate many practical classification problems to QA-based techniques. Further, we show that synthetic data generated from a RBM can be used to balance the original dataset.

97 MATHEMATICS AND COMPUTING↗

Photodissociation of HBr/LiF(001): A quantum mechanical model

The photodissociation dynamics of HBr adsorbed on a LiF(001) surface is studied using time-independent quantum mechanics. The photodissociation lineshape and the Br(P(sub 1/2)-2)/Br(P(sub 3/2)-2) yield ratio are computed and compared with the corresponding quantities for gas phase photodissociation. The angular distribution of the hydrogen photofragments following excitation of adsorbed HBr is computed and found to agree qualitatively with experimental data. The effect of polarization of the photon is illustrated and discussed. The field polarization is found to affect significantly the magnitude of the photodissociation signal but not the angular dependence of the photofragment distribution, in agreement with experiment and in accord with expectations for a strongly aligned adsorbed phase.

Seideman, Tamar↗

O(N) Quantum Molecular Dynamics (QMD) Algorithms, Code Requirements, Release Notes, and Library Deployment

The Exascale Computing Project (ECP) Co-design Center for Particle Applications (CoPA) is addressing the challenges for particle applications to run on upcoming exascale computing architectures. This includes short- and long-range molecular dynamics (MD), particle-in-cell (PIC), and quantum molecular dynamics (QMD). This document addresses $\mathcal{O}(N)$ QMD. The goal of this CoPA team is to develop $\mathcal{O}(N)$ linear scaling methods for QMD optimized for new and next generation architectures. This will provide the capability to simulate larger systems and run for longer timeframes than is possible today. Focus is initially on the density matrix construction. This requires efficient matrix formats and linear algebra operations for distributed architectures. A significant reduction in the wallclock time for an MD timestep is necessary to achieve practical large scale QMD. Equally challenging is to develop scalable distributed implementations. This requires the exploration of efficient data formats (dense, sparse, etc.), data decomposition and distribution across compute nodes using block- and graph-based methods. Traditional programming models such as MPI/OpenMP, hybrid MPI/OpenMP/Accelerator, as well as asynchronous task-based models will be explored to achieve maximum parallelism. Parallelism at the level of a matrix operation, solver, and/or application needs to be explored to gain the most benefit on new and emerging architectures. This report serves as a “living document” of the additions and enhancements throughout the ECP CoPA project for $\mathcal{O}(N)$ QMD. The following sections address the applications we are working with, the quantum chemistry libraries (BML and PROGRESS), and our co-design exploration through the use of a proxy application.

36 MATERIALS SCIENCE↗

The Quest for Accurate Liquid Water Properties from First Principles

Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies, as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of Density Functional Theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, displays accuracy comparable to a computationally expensive dispersion corrected hybrid functional, revPBE0-D 3 . We show that the meta-GGA DFT reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared (IR) spectra with comparable accuracy of a much more expensive hybrid functional. In conclusion, this work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei, and that high temperature simulations are a poor proxy for properly treating nuclear quantum effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗