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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Advanced architectures for high-performance quantum networking

As practical quantum networks prepare to serve an ever-expanding number of nodes, there has grown a need for advanced auxiliary classical systems that support the quantum protocols and maintain compatibility with the existing fiber-optic infrastructure. We propose and demonstrate a quantum local area network design that addresses current deployment limitations in timing and security in a scalable fashion using commercial off-the-shelf components. First, we employ White Rabbit switches to synchronize three remote nodes with ultra-low timing jitter, significantly increasing the fidelities of the distributed entangled states over previous work with Global Positioning System clocks. Second, using a parallel quantum key distribution channel, we secure the classical communications needed for instrument control and data management. Therefore, the conventional network that manages our entanglement network is secured using keys generated via an underlying quantum key distribution layer, preserving the integrity of the supporting systems and the relevant data in a future-proof fashion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A Hartree-Fock Application Using UPC++ and the New DArray Library

The Hartree-Fock (HF) method is the fundamental first step for incorporating quantum mechanics into many-electron simulations of atoms and molecules, and it is an important component of computational chemistry toolkits like NWChem. The GTFock code is an HF implementation that, while it does not have all the features in NWChem, represents crucial algorithmic advances that reduce communication and improve load balance by doing an up-front static partitioning of tasks, followed by work stealing whenever necessary. To enable innovations in algorithms and exploit next generation exascale systems, it is crucial to support quantum chemistry codes using expressive and convenient programming models and runtime systems that are also efficient and scalable. Here, this paper presents an HF implementation similar to GTFock using UPC++, a partitioned global address space model that includes flexible communication, asynchronous remote computation, and a powerful multidimensional array library. UPC++ offers runtime features that are useful for HF such as active messages, a rich calculus for array operations, hardware-supported fetch-and-add, and functions for ensuring asynchronous runtime progress. We present a new distributed array abstraction, DArray, that is convenient for the kinds of random-access array updates and linear algebra operations on block-distributed arrays with irregular data ownership. Finally, we analyze the performance of atomic fetch-and-add operations (relevant for load balancing) and runtime attentiveness, then compare various techniques and optimizations for each. Our optimized implementation of HF using UPC++ and the DArrays library shows up to 20% improvement over GTFock with Global Arrays at scales up to 24,000 cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dissolving salt is not equivalent to applying a pressure on water

Abstract Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

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

A new approach to three-dimensional neutron transport solution based on the method of characteristics and linear axial approximation

A new approach based on the method of characteristics (MOC) is proposed to solve the neutron transport equation. A new three-dimensional (3D) spatial discretization is applied to avoid the instability issue of the transverse leakage iteration of the traditional 2D/1D approach. In this new approach, the axial and radial variables are discretized in two different ways: the linear expansion is performed in the axial direction, then, the 3D solution of the angular flux is transformed to be the planar solution of 2D angular expansion moments, which are solved by the planar MOC sweeping. Based on the boundary and interface continuity conditions, the 2D expansion moment solution is equivalently transformed to be the solution of the axially averaged angular flux. Using the piecewise averaged angular flux at the top and bottom surfaces of 3D meshes, the planes are coupled to give the 3D angular flux distribution. The 3D CMFD linear system is established from the surface net current of every 3D pin-mesh to accelerate the convergence of power iteration. The STREAM code is extended to be capable of handling 3D problems based on the new approach. Several benchmarks are tested to verify its feasibility and accuracy, including the 3D homogeneous benchmarks and heterogeneous benchmarks. The computational sensitivity is discussed. The results show good accuracy in all tests. With the CMFD acceleration, the convergence is stable. In addition, a pin-cell problem with void gap is calculated. This shows the advantage compared to the traditional 2D/1D MOC methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Quantization of the Electromagnetic Field in Three-Dimensional Photonic Structures on the Basis of the Scattering Matrix Formalism (S Quantization)

A technique for quantization of the electromagnetic field in photonic nanostructures with three-dimensional modulation of the dielectric constant is developed on the basis of the scattering matrix formalism (S quantization in the three-dimensional case). Quantization is based on equating the eigenvalues of the scattering matrix to unity, which is equivalent to equating each other the sets of Fourier expansions for the fields of the waves incident on the structure and propagating away from the structure. The spatial distribution of electromagnetic fields of the modes in a photonic nanostructure is calculated on the basis of the R and T matrices describing the reflection and transmission of the Fourier components by the structure. To calculate the reflection and transmission coefficients of individual real-space and Fourier-space components, the structure is divided into parallel layers within which the dielectric constant varies as a function of two-dimensional coordinates. Using the Fourier transform, Maxwell’s equations are written in the form of a matrix connecting the Fourier components of the electric field at the boundaries of neighboring layers. Based on the calculated reflection and transmission vectors for all polarizations and Fourier components, the scattering matrix for the entire structure is formed and quantization is carried out by equating the eigenvalues of the scattering matrix to unity. The developed method makes it possible to obtain the spatial profiles of eigenmodes without solving a system of nonlinear integro-differential equations and significantly reduces the computational resources required for calculating the probability of spontaneous emission in three-dimensional systems.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Computational Study of Symmetric Methylation on Histone Arginine Catalyzed by Protein Arginine Methyltransferase PRMT5 through QM/MM MD and Free Energy Simulations

Protein arginine methyltransferases (PRMTs) catalyze the transfer of the methyl group from S-adenosyl-l-methionine (AdoMet) to arginine residues. There are three types of PRMTs (I, II and III) that produce different methylation products, including asymmetric dimethylarginine (ADMA), symmetric dimethylarginine (SDMA) and monomethylarginine (MMA). Since these different methylations can lead to different biological consequences, understanding the origin of product specificity of PRMTs is of considerable interest. In this article, the quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed to study SDMA catalyzed by the Type II PRMT5 on the basis of experimental observation that the dimethylated product is generated through a distributive fashion. The simulations have identified some important interactions and proton transfers during the catalysis. Similar to the cases involving Type I PRMTs, a conserved Glu residue (Glu435) in PRMT5 is suggested to function as general base catalyst based on the result of the simulations. Moreover, our results show that PRMT5 has an energetic preference for the first methylation on N-1 followed by the second methylation on a different -guanidino nitrogen of arginine (N-2).The first and second methyl transfers are estimated to have free energy barriers of 19-20 and 18-19 kcal/mol respectively. The computer simulations suggest a distinctive catalytic mechanism of symmetric dimethylation that seems to be different from asymmetric dimethylation.

ARABIDOPSIS↗
DOE OSTI2015

An excited-state approach within full configuration interaction quantum Monte Carlo

We present a new approach to calculate excited states with the full configuration interaction quantum Monte Carlo (FCIQMC) method. The approach uses a Gram-Schmidt procedure, instantaneously applied to the stochastically evolving distributions of walkers, to orthogonalize higher energy states against lower energy ones. It can thus be used to study several of the lowest-energy states of a system within the same symmetry. This additional step is particularly simple and computationally inexpensive, requiring only a small change to the underlying FCIQMC algorithm. No trial wave functions or partitioning of the space is needed. The approach should allow excited states to be studied for systems similar to those accessible to the ground-state method due to a comparable computational cost. As a first application, we consider the carbon dimer in basis sets up to quadruple-zeta quality and compare to existing results where available.

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

Exact spectral gaps of random one-dimensional quantum circuits

The spectral gap of local random quantum circuits is a fundamental property that determines how close the moments of the circuit's unitaries match those of a Haar random distribution. When studying spectral gaps, it is common to bound these quantities using tools from statistical mechanics or via quantum information-based inequalities. Here, by focusing on the second moment of one-dimensional unitary circuits where nearest-neighboring gates act on sets of qudits (with open and closed boundary conditions), we show that one can exactly compute the associated spectral gaps. Indeed, having access to their functional form allows us to prove several important results, such as the fact that the spectral gap for closed boundary condition is exactly the square of the gap for open boundaries, as well as improve on previously known bounds for approximate design convergence. Finally, we verify our theoretical results by numerically computing the spectral gap for systems of up to 70 qubits, as well as comparing them to gaps of random orthogonal and symplectic circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Adaptive mitigation of time-varying quantum noise

Current quantum computers suffer from non-stationary noise channels with high error rates, which undermines their reliability and reproducibility. We propose a Bayesian inference based adaptive algorithm that can learn and mitigate quantum noise in response to changing channel conditions. Our study emphasizes the need for dynamic inference of critical channel parameters to improve program accuracy. We use the Dirichlet distribution to model the stochasticity of the Pauli channel. This allows us to perform Bayesian inference, which can improve the performance of probabilistic error cancellation (PEC) under time-varying noise. Our work demonstrates the importance of characterizing and mitigating temporal variations in quantum noise, which is crucial for developing more accurate and reliable quantum technologies. Our results demonstrate that Bayesian PEC can outperform non-adaptive approaches by a factor of 4.5x when measured using Hellinger distance from the ideal distribution.

Dasgupta, Samudra↗
DOE OSTI2023

A dynamic parameterization of sulfuric acid–dimethylamine nucleation and its application in three-dimensional modeling

Sulfuric acid (SA) is a governing gaseous precursor for atmospheric new particle formation (NPF), a major source of global ultrafine particles, in environments studied around the world. In polluted urban atmospheres with high condensation sinks (CSs), the formation of stable SA–amine clusters, such as SA–dimethylamine (DMA) clusters, usually initializes intense NPF events. Coagulation scavenging and cluster evaporation are dominant sink processes of SA–amine clusters in urban atmospheres, yet these loss processes are not quantitatively included in the present parameterizations of SA–amine nucleation. We herein report a parameterization of SA–DMA nucleation, based on cluster dynamic simulations and quantum chemistry calculations, with certain simplifications to greatly reduce the computational costs. Compared with previous SA–DMA nucleation parameterizations, this new parameterization was able to reproduce the dependences of particle formation rates on temperature and CSs. We then incorporated it in a three-dimensional (3-D) chemical transport model to simulate the evolution of the particle number size distributions. Simulation results showed good consistency with the observations in the occurrence of NPF events and particle number size distributions in wintertime Beijing and represented a significant improvement compared to that using a parameterization without coagulation scavenging. Quantitative analysis shows that SA–DMA nucleation contributes significantly to nucleation rates and aerosol population during the 3-D simulations in Beijing (>99 % and >60 %, respectively). These results broaden the understanding of NPF in urban atmospheres and stress the necessity of including the effects of coagulation scavenging and cluster stability in simulating SA–DMA nucleation in 3-D simulations. Representing these processes is thus likely to improve model performance in particle source apportionment and quantification of aerosol effects on air quality, human health, and climate.

Li, Yuyang↗
NASA NTRS1988

Fluctuation analysis of relativistic nucleus-nucleus collisions in emulsion chambers

An analytical technique was developed for identifying enhanced fluctuations in the angular distributions of secondary particles produced from relativistic nucleus-nucleus collisions. The method is applied under the assumption that the masses of the produced particles are small compared to their linear momenta. The importance of particles rests in the fact that enhanced fluctuations in the rapidity distributions is considered to be an experimental signal for the creation of the quark-gluon-plasma (QGP), a state of nuclear matter predicted from the quantum chromodynamics theory (QCD). In the approach, Monte Carlo simulations are employed that make use of a portable random member generator that allow the calculations to be performed on a desk-top computer. The method is illustrated with data taken from high altitude emulsion exposures and is immediately applicable to similar data from accelerator-based emulsion exposures.

Mcguire, Stephen C.↗
DOE OSTI2020

Virtualized Logical Qubits: A 2.5D Architecture for Error-Corrected Quantum Computing

Current, near-term quantum devices have shown great progress in the last several years culminating recently with a demonstration of quantum supremacy. In the medium-term, however, quantum machines will need to transition to greater reliability through error correction, likely through promising techniques like surface codes which are well suited for near-term devices with limited qubit connectivity. We discover quantum memory, particularly resonant cavities with transmon qubits arranged in a 2.5D architecture, can efficiently implement surface codes with substantial hardware savings and performance/fidelity gains. Specifically, we virtualize logical qubits by storing them in layers of qubit memories connected to each transmon. Surprisingly, distributing each logical qubit across many memories has a minimal impact on fault tolerance and results in substantially more efficient operations. Our design permits fast transversal application of CNOT operations between logical qubits sharing the same physical address (same set of cavities) which are 6x faster than standard lattice surgery CNOTs. We develop a novel embedding which saves approximately 10x in transmons with another 2x savings from an additional optimization for compactness. Although qubit virtualization pays a 10x penalty in serialization, advantages in the transversal CNOT and in area efficiency result in fault-tolerance and performance comparable to conventional 2D transmon-only architectures. Our simulations show our system can achieve fault tolerance comparable to conventional two-dimensional grids while saving substantial hardware. Furthermore, our architecture can produce magic states at 1.22x the baseline rate given a fixed number of transmon qubits. Here, this is a critical benchmark for future fault-tolerant quantum computers as magic states are essential and machines will spend the majority of their resources continuously producing them. This architecture substantially reduces the hardware requirements for fault-tolerant quantum computing and puts within reach a proof-of-concept experimental demonstration of around 10 logical qubits, requiring only 11 transmons and 9 attached cavities in total.

quantum computing↗
DOE OSTI2025

Structural constraint integration in a generative model for the discovery of quantum materials

Billions of organic molecules have been computationally generated, yet functional inorganic materials remain scarce due to limited data and structural complexity. Here, in this work, we introduce Structural Constraint Integration in a GENerative model (SCIGEN), a framework that enforces geometric constraints, such as honeycomb and kagome lattices, within diffusion-based generative models to discover stable quantum materials candidates. SCIGEN enables conditional sampling from the original distribution, preserving output validity while guiding structural motifs. This approach generates ten million inorganic compounds with Archimedean and Lieb lattices, over 10% of which pass multistage stability screening. High-throughput density functional theory calculations on 26,000 candidates shows over 95% convergence and 53% structural stability. A graph neural network classifier detects magnetic ordering in 41% of relaxed structures. Furthermore, we synthesize and characterize two predicted materials, TiPd 0.22 Bi 0.88 and Ti 0.5 Pd 1.5 Sb, which display paramagnetic and diamagnetic behaviour, respectively. Our results indicate that SCIGEN provides a scalable path for generating quantum materials guided by lattice geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theory of resonant x-ray scattering with ultrafast intense pulses

Here, we present a time-dependent Schrödinger equation approach within a nonrelativistic quantum electrodynamics framework to investigate resonant x-ray scattering driven by intense x-ray pulses. This method enables us to explore how coherent x-ray electron dynamics influence scattering signals from Ne + . We account for both resonance fluorescence and elastic scattering channels, while also considering competing photoionization and inner-shell decay processes. By computing the angular distribution and energy spectrum of scattered photons, we uncover interference effects between elastic scattering and resonance fluorescence pathways. Notably, this interference results in a small asymmetry in the energy spectrum. We discuss the experimental potential for detecting signatures of interference. Our findings demonstrate that the x-ray Rabi dynamics can be used to control scattering responses and provide insights into interference mechanisms and scattering efficiency in high-intensity x-ray regimes.

Venkatesh, Akilesh [Argonne National Laboratory (A↗
DOE OSTI2014

Hydration dynamics in water clusters via quantum molecular dynamics simulations

We have investigated the hydration dynamics in size selected water clusters with n = 66, 104, 200, 500, and 1000 water molecules using molecular dynamics simulations. To study the most fundamental aspects of relaxation phenomena in clusters, we choose one of the simplest, still realistic, quantum mechanically treated test solute, an excess electron. The project focuses on the time evolution of the clusters following two processes, electron attachment to neutral equilibrated water clusters and electron detachment from an equilibrated water cluster anion. The relaxation dynamics is significantly different in the two processes, most notably restoring the equilibrium final state is less effective after electron attachment. Nevertheless, in both scenarios only minor cluster size dependence is observed. Significantly different relaxation patterns characterize electron detachment for interior and surface state clusters, interior state clusters relaxing significantly faster. This observation may indicate a potential way to distinguish surface state and interior state water cluster anion isomers experimentally. A comparison of equilibrium and non-equilibrium trajectories suggests that linear response theory breaks down for electron attachment at 200 K, but the results converge to reasonable agreement at higher temperatures. Relaxation following electron detachment clearly belongs to the linear regime. Cluster relaxation was also investigated using two different computational models, one preferring cavity type interior states for the excess electron in bulk water, while the other simulating non-cavity structure. While the cavity model predicts appearance of several different hydrated electron isomers in agreement with experiment, the non-cavity model locates only cluster anions with interior excess electron distribution. The present simulations show that surface isomers computed with the cavity predicting potential show similar dynamical behavior to the interior clusters of the non-cavity type model. Relaxation associated with cavity collapse presents, however, unique dynamical signatures.

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

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗
DOE OSTI2026

Fantômas unconfined: global QCD fits with Bézier parameterizations

Fantômas is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). Fantômas facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the xFitter framework for the global QCD analysis, Fantômas provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the Fantômas toolkit, including its inputs, outputs, and implementation within xFitter. It also provides a practical example of using Fantômas for uncertainty quantification as well as the combination of PDF fits into a single ensemble.

Bézier curves↗