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 793 records · Page 44

Spatially adaptive stochastic methods for fluid–structure interactions subject to thermal fluctuations in domains with complex geometries

We develop stochastic mixed finite element methods for spatially adaptive simulations of fluid–structure interactions when subject to thermal fluctuations. To account for thermal fluctuations, we introduce a discrete fluctuation–dissipation balance condition to develop compatible stochastic driving fields for our discretization. We perform analysis that shows our condition is sufficient to ensure results consistent with statistical mechanics. We show the Gibbs–Boltzmann distribution is invariant under the stochastic dynamics of the semi-discretization. To generate efficiently the required stochastic driving fields, we develop a Gibbs sampler based on iterative methods and multigrid to generate fields with O(N) computational complexity. Our stochastic methods provide an alternative to uniform discretizations on periodic domains that rely on Fast Fourier Transforms. To demonstrate in practice our stochastic computational methods, we investigate within channel geometries having internal obstacles and no-slip walls how the mobility/diffusivity of particles depends on location. Our methods extend the applicability of fluctuating hydrodynamic approaches by allowing for spatially adaptive resolution of the mechanics and for domains that have complex geometries relevant in many applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultracold chemistry with alkali-metal–rare-earth molecules

Here, a first principles study of the dynamics of 6 Li( 2 S) + 6 Li 174 Yb( 2 Σ + ) → 6 Li 2 ( 1 Σ + ) + 174 Yb( 1 S) reaction is presented at cold and ultracold temperatures. The computations involve determination and analytic fitting of a three-dimensional potential energy surface for the Li 2 Yb system and quantum dynamics calculations of varying complexities, ranging from exact quantum dynamics within the close-coupling scheme, to statistical quantum treatment, and universal models. It is demonstrated that the two simplified methods yield zero-temperature limiting reaction rate coefficients in reasonable agreement with the full close-coupling calculations. The effect of the three-body term in the interaction potential is explored by comparing quantum dynamics results from a pairwise potential that neglects the three-body term to that derived from the full interaction potential. Inclusion of the three-body term in the close-coupling calculations was found to reduce the limiting rate coefficients by a factor of two. The reaction exoergicity populates vibrational levels as high as v = 19 of the 6 Li 2 molecule in the limit of zero collision energy. Product vibrational distributions from the close-coupling calculations reveal sensitivity to inclusion of three-body forces in the interaction potential. In conclusion, the results indicate that a simplified model based on the long-range potential is able to yield reliable values of the total reaction rate coefficient in the ultracold limit but a more rigorous approach based on statistical quantum or quantum close-coupling methods is desirable when product rovibrational distribution is required.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗

Rational Chebyshev spectral transform for the dynamics of broad-area laser diodes

This manuscript details the use of the rational Chebyshev transform for describing the transverse dynamics of broad-area laser diodes and amplifiers. This spectral method can be used in combination with the delay algebraic equations approach developed in [1], which substantially reduces the computation time. The theory is presented in such a way that it encompasses the case of the Fourier spectral transform presented in [2] as a particular case. It is also extended to the consideration of index guiding with an arbitrary transverse profile. Because their domain of definition is infinite, the convergence properties of the Chebyshev rational functions allow handling the boundary conditions with higher accuracy than with the previously studied Fourier transform method. As practical examples, we solve the beam propagation problem with and without index guiding: we obtain excellent results and an improvement of the integration time between one and two orders of magnitude as compared with a fully distributed two dimensional model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Explaining an unusual electromigration behavior—A comprehensive experimental and theoretical analysis using finite element method

In metallic interconnects, it is generally assumed that electromigration (EM) failure location is independent of the applied electrical current and always occurs at the highest-current-density area. Our experiments show otherwise. We designed an Al interconnect that alters its failure location by only varying the applied current density. The failure occurs near the high for a current above 2 × 10 7 A/cm 2 , but at a location with 59% of the maximum for lower current densities. Thermoreflectance thermal imaging is employed to gather time-dependent high-resolution spatial temperature distributions of the Al interconnect during EM. More importantly, we propose a computationally inexpensive 2D finite element method that tracks EM evolution in time and matches well with the observations from different experimental conditions. A detailed analysis covering the major driving forces of EM is carried out to understand the complex physics behind EM. The atomic depletion rate contributed by each force is quantitatively studied. By examining the results from every tested experimental condition, the model reveals that the temperature gradient is the key reason causing atomic depletion near the failure location. Graphical illustrations and qualitative analysis are provided to intuitively show the key findings of our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Warm dark matter via ultra-violet freeze-in: reheating temperature and non-thermal distribution for fermionic Higgs portal dark matter

Warm dark matter (WDM) of order keV mass may be able to resolve the disagreement between structure formation in cold dark matter simulations and observations. The detailed properties of WDM will depend upon its energy distribution, in particular how it deviates from the thermal distribution usually assumed in WDM simulations. Here we focus on WDM production via the Ultra-Violet (UV) freeze-in mechanism, for the case of fermionic Higgs portal dark matter ψ produced via the portal interaction ψ-barψH{sup †}H/Λ. We introduce a new method to simplify the computation of the non-thermal energy distribution of dark matter from freeze-in. We show that the non-thermal energy distribution from UV freeze-in is hotter than the corresponding thermal distribution and has the form of a Bose-Einstein distribution with a non-thermal normalization. The resulting range of dark matter fermion mass consistent with observations is 5–7 keV. The reheating temperature must satisfy T{sub R}≳120 GeV in order to account for the observed dark matter density when m{sub ψ}≈5 keV, where the lower bound on T{sub R} corresponds to the limit where the fermion mass is entirely due to electroweak symmetry breaking via the portal interaction. The corresponding bound on the interaction scale is Λ≳6.0×10{sup 9} GeV.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Aggregation-fragmentation-diffusion model for trail dynamics

We investigate statistical properties of trails formed by a random process incorporating aggregation, fragmentation, and diffusion. In this stochastic process, which takes place in one spatial dimension, two neighboring trails may combine to form a larger one, and also one trail may split into two. In addition, trails move diffusively. The model is defined by two parameters which quantify the fragmentation rate and the fragment size. In the long-time limit, the system reaches a steady state, and our focus is the limiting distribution of trail weights. We find that the density of trail weight has power-law tail P(w)~w –γ for small weight w. We obtain the exponent γ analytically and find that it varies continuously with the two model parameters. In conclusion, the exponent γ can be positive or negative, so that in one range of parameters small-weight trails are abundant and in the complementary range they are rare.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modular Quantum Processor with an All-to-All Reconfigurable Router

Superconducting qubits provide a promising approach to large-scale fault-tolerant quantum computing. However, qubit connectivity on a planar surface is typically restricted to only a few neighboring qubits. Achieving longer-range and more flexible connectivity, which is particularly appealing in light of recent developments in error-correcting codes, however, usually involves complex multilayer packaging and external cabling, which is resource intensive and can impose fidelity limitations. Here, we propose and realize a high-speed on-chip quantum processor that supports reconfigurable all-to-all coupling with a large on-off ratio. We implement the design in a four-node quantum processor, built with a modular design comprising a wiring substrate coupled to two separate qubit-bearing substrates, each including two single-qubit nodes. We use this device to demonstrate reconfigurable controlled-𝑍 gates across all qubit pairs, with a benchmarked average fidelity of 96.00% ± 0.08% and best fidelity of 97.14% ± 0.07%, limited mainly by dephasing in the qubits. We also generate multiqubit entanglement, distributed across the separate modules, demonstrating GHZ-3 and GHZ-4 states with fidelities of 88.15% ± 0.24% and 75.18% ± 0.11%, respectively. This approach promises efficient scaling to larger-scale quantum circuits and offers a pathway for implementing quantum algorithms and error-correction schemes that benefit from enhanced qubit connectivity.

Quantum circuits↗

Freeze-in, glaciation, and UV sensitivity from light mediators

Dark matter (DM) freeze-in through a light mediator is an appealing model with excellent detection prospects at current and future experiments. Light mediator freeze-in is UV-insensitive insofar as most DM is produced at late times, and thus the DM abundance does not depend on the unknown early evolution of our universe. However the final DM yield retains a dependence on the initial DM population, which is usually assumed to be exactly zero. We point out that in models with light mediators, the final DM yield will also depend on the initial conditions assumed for the light mediator population. We describe a class of scenarios we call “glaciation” where DM freezing in from the SM encounters a pre-existing thermal bath of mediators, and study the dependence of the final DM yield on the initial temperature of this dark radiation bath. To compute DM scattering rates in this cosmology, we derive for the first time an exact integral expression for the Boltzmann collision term describing interactions between two species at different temperatures. We quantify the dependence of the DM yield on the initial dark temperature and find that it can be sizeable in regions near the traditional (zero initial abundance) freeze-in curve. We generalize the freeze-in curve to a glaciation band, which can extend as much as an order of magnitude below the traditional freeze-in direct detection target, and point out that the DM phase space distribution as well as the yield can be strongly dependent on initial conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fast equilibration protocol for million atom systems of highly entangled linear polyethylene chains

Equilibrated systems of entangled polymer melts cannot be produced using direct brute force equilibration due to the slow reptation dynamics exhibited by high molecular weight chains. Instead, these dense systems are produced using computational techniques such as Monte Carlo-Molecular Dynamics hybrid algorithms, though the use of soft potentials has also shown promise mainly for coarse-grained polymeric systems. Through the use of soft-potentials, the melt can be equilibrated via molecular dynamics at intermediate and long length scales prior to switching to a Lennard-Jones potential. We will outline two different equilibration protocols, which use various degrees of information to produce the starting configurations. In one protocol, we use only the equilibrium bond angle, bond length, and target density during the construction of the simulation cell, where the information is obtained from available experimental data and extracted from the force field without performing any prior simulation. In the second protocol, we moreover utilize the equilibrium radial distribution function and dihedral angle distribution. This information can be obtained from experimental data or from a simulation of short unentangled chains. Both methods can be used to prepare equilibrated and highly entangled systems, but the second protocol is much more computationally efficient. These systems can be strictly monodisperse or optionally polydisperse depending on the starting chain distribution. Our protocols, which utilize a soft-core harmonic potential, will be applied for the first time to equilibrate a million particle system of polyethylene chains consisting of 1000 united atoms at various temperatures. Calculations of structural and entanglement properties demonstrate that this method can be used as an alternative towards the generation of entangled equilibrium structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transition micromagnetic structures in asymmetric vortexlike domain walls (Static solutions and dynamic reconstructions)

The possible types of transition structures with a three-dimensional magnetization distribution over regions in the vortex asymmetric domain walls that exist in magnetically uniaxial soft magnetic films with in-plane anisotropy are studied by computer simulation in terms of a micromagnetic approach. It is shown that the possible structure types include both the type of vertical Bloch lines that was discussed earlier in other works and new types, namely, singular (Bloch) points and clusters consisting of vertical Bloch lines and Bloch points. The spatial configurations of the transition structures are calculated and their topological properties are found. The numerical simulation of the dynamics of closely spaced substructure regions reveals various scenarios of their interaction, including annihilation accompanied by energy release and the excitation of nonlinear waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Symmetric Blind Information Reconciliation for Quantum Key Distribution

Quantum key distribution (QKD) is a quantum-proof key-exchange scheme which is fast approaching the communication industry. An essential component in QKD is the information reconciliation step, which is used for correcting the quantum-channel noise errors. The recently suggested blind-reconciliation technique, based on low-density parity-check codes, offers remarkable prospectives for efficient information reconciliation without an a priori quantum bit error rate estimation. We suggest an improvement of the blind-information-reconciliation protocol promoting a significant increase in the efficiency of the procedure and reducing its interactivity. Finally, the proposed technique is based on introducing symmetry in operations of parties, and the consideration of results of unsuccessful belief-propagation decodings.

97 MATHEMATICS AND COMPUTING↗

Π4U: A high performance computing framework for Bayesian uncertainty quantification of complex models

We present Π4U,{sup 1} an extensible framework, for non-intrusive Bayesian Uncertainty Quantification and Propagation (UQ+P) of complex and computationally demanding physical models, that can exploit massively parallel computer architectures. The framework incorporates Laplace asymptotic approximations as well as stochastic algorithms, along with distributed numerical differentiation and task-based parallelism for heterogeneous clusters. Sampling is based on the Transitional Markov Chain Monte Carlo (TMCMC) algorithm and its variants. The optimization tasks associated with the asymptotic approximations are treated via the Covariance Matrix Adaptation Evolution Strategy (CMA-ES). A modified subset simulation method is used for posterior reliability measurements of rare events. The framework accommodates scheduling of multiple physical model evaluations based on an adaptive load balancing library and shows excellent scalability. In addition to the software framework, we also provide guidelines as to the applicability and efficiency of Bayesian tools when applied to computationally demanding physical models. Theoretical and computational developments are demonstrated with applications drawn from molecular dynamics, structural dynamics and granular flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Characterization of the Reflected Quasipotential

Recent interest in the reflected quasipotential comes from the queueing theory literature, specifically the analysis of so-called (b,A,D) reflected Brownian motion where it is the large deviation rate function for the stationary distribution. Our purpose here is to characterize the reflected quasipotential in terms of a first-order Hamilton–Jacobi equation. Using conventional dynamic programming ideas, along with a complementarity problem formulation of the effect of the Skorokhod map on absolutely continuous paths, we will derive necessary conditions in the form of viscosity-sense boundary conditions. It turns out that even with these boundary conditions solutions are not unique. Thus a unique characterization needs to refer to some additional property of V(⋅). We establish such a characterization in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Superior model for fault tolerance computation in designing nano-sized circuit systems

As CMOS technology scales nano-metrically, reliability turns out to be a decisive subject in the design methodology of nano-sized circuit systems. As a result, several computational approaches have been developed to compute and evaluate reliability of desired nano-electronic circuits. The process of computing reliability becomes very troublesome and time consuming as the computational complexity build ups with the desired circuit size. Therefore, being able to measure reliability instantly and superiorly is fast becoming necessary in designing modern logic integrated circuits. For this purpose, the paper firstly looks into the development of an automated reliability evaluation tool based on the generalization of Probabilistic Gate Model (PGM) and Boolean Difference-based Error Calculator (BDEC) models. The Matlab-based tool allows users to significantly speed-up the task of reliability analysis for very large number of nano-electronic circuits. Secondly, by using the developed automated tool, the paper explores into a comparative study involving reliability computation and evaluation by PGM and, BDEC models for different implementations of same functionality circuits. Based on the reliability analysis, BDEC gives exact and transparent reliability measures, but as the complexity of the same functionality circuits with respect to gate error increases, reliability measure by BDEC tends to be lower than the reliability measure by PGM. The lesser reliability measure by BDEC is well explained in this paper using distribution of different signal input patterns overtime for same functionality circuits. Simulation results conclude that the reliability measure by BDEC depends not only on faulty gates but it also depends on circuit topology, probability of input signals being one or zero and also probability of error on signal lines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗