Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum technologies”

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

On numerical energy conservation for an implicit particle-in-cell method coupled with a binary Monte-Carlo algorithm for Coulomb collisions

We report conventional particle-in-cell (PIC) methods suffer from enhanced numerical heating (explicit PIC) or cooling (semi-implicit PIC) when coupled with a binary Monte-Carlo algorithm for Coulomb collisions. In this work, a fully-implicit θ-PIC scheme (with adjustable time-biasing parameter 1/2 ≤ θ ≤) is considered. The discrete change in energy of a closed system after a time step for this scheme scales with (1/2 - θ)C θ , where C θ is a positive definite quantity that depends on the frequency spectrum of the energy in the fields. Collisions lead to additional energy in the field fluctuations associated with high-frequency light waves produced by a numerical Bremsstrahlung process, which can result in a large increase in the numerical cooling rate for θ > 1/2. However, for θ = 1/2, energy is exactly conserved. The energy in the field fluctuations on long time scales agrees with that calculated using the equipartition theorem for a classical system in thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stability in integrable nonlocal nonlinear equations

Recently a variety of nonlocal integrable systems has been introduced that besides fields located at particular space-time points simultaneously also contain fields that are located at different, but symmetrically related, points. Here we investigate different types of soliton solutions with regard to their stability against linear perturbations obtained for the nonlocal version of the Hirota/nonlinear Schrödinger equation and the so-called Alice and Bob versions of the Korteweg-de Vries and Bousinesq equations. In this work, we encounter different types of scenarios: Solition solutions that are linearly stable or unstable and also solutions that change their stability properties depending on the parameter regime they are in.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Simple, general, realistic, robust, analytic tokamak equilibria. Part 1. Limiter and divertor tokamaks

Tokamak equilibria have been derived that are analytic solutions to the Grad–Shafranov equation. Here in this paper, Part 1, describes a wide range of such equilibria including smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to include edge pedestals and toroidal flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evolution of an arched magnetized laboratory plasma in a sheared magnetic field

Arched magnetized structures are a common occurrence in space and laboratory plasmas. Results from a laboratory experiment on spatio-temporal evolution of an arched magnetized plasma ( $\beta \approx 10^{-3}$ , Lundquist number $\approx 10^{4}$ , plasma radius/ion gyroradius $\approx 20$ ) in a sheared magnetic configuration are presented. The experiment is designed to model conditions relevant to the formation and destabilization of similar structures in the solar atmosphere. The magnitude of a nearly horizontal overlying magnetic field was varied to study its effects on the writhe and twist of the arched plasma. In addition, the direction of the guiding magnetic field along the arch was varied to investigate its role in the formation of either forward- or reverse-S shaped plasma structures. The electrical current in the arched plasma was well below the current required to make it kink unstable. A significant increase in the writhe of the arched plasma was observed with larger magnitudes of overlying magnetic field. A forward-S shaped arched plasma was observed for a guiding magnetic field oriented nearly antiparallel to the initial arched plasma current, while the parallel orientation yielded the reverse-S shaped arched plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Global stochastic optimization of stellarator coil configurations

In the construction of a stellarator, the manufacturing and assembling of the coil system is a dominant cost. These coils need to satisfy strict engineering tolerances, and if those are not met the project could be cancelled as in the case of the National Compact Stellarator Experiment (NCSX) project. Therefore, our goal is to find coil configurations that increase construction tolerances without compromising the performance of the magnetic field. In this paper, we develop a gradient-based stochastic optimization model which seeks robust stellarator coil configurations in high dimensions. In particular, we design a two-step method: first, we perform an approximate global search by a sample efficient trust-region Bayesian optimization; second, we refine the minima found in step one with a stochastic local optimizer. To this end, we introduce two stochastic local optimizers: BFGS applied to the sample average approximation; and Adam, equipped with a control variate for variance reduction. Numerical simulations performed on a W7-X-like coil configuration demonstrate that our global optimization approach finds a variety of promising local solutions at less than 0.1% of the cost of previous work, which considered solely local stochastic optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Partitioning Quantum Chemistry Simulations with Clifford Circuits

Current quantum computing hardware is restricted by the availability of only few, noisy qubits which limits the investigation of larger, more complex molecules in quantum chemistry calculations on quantum computers in the near term. Here, in this work, we investigate the limits of their classical and near-classical treatment while staying within the framework of quantum circuits and the variational quantum eigensolver. To this end, we consider naive and physically motivated, classically efficient product ansatz for the parametrized wavefunction adapting the separable-pair ansatz form. We combine it with post-treatment to account for interactions between subsystems originating from this ansatz. The classical treatment is given by another quantum circuit that has support between the enforced subsystems and is folded into the Hamiltonian. To avoid an exponential increase in the number of Hamiltonian terms, the entangling operations are constructed from purely Clifford or near-Clifford circuits. While Clifford circuits can be simulated efficiently classically, they are not universal. In order to account for missing expressibility, near-Clifford circuits with only few, selected non-Clifford gates are employed. The exact circuit structure to achieve this objective is molecule-dependent and is constructed using simulated annealing and genetic algorithms. We demonstrate our approach on a set of molecules of interest and investigate the extent of our methodology’s reach.

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

Collection of fluorescence from an ion using trap-integrated photonics

Spontaneously emitted photons are entangled with the electronic and nuclear degrees of freedom of the emitting atom, so interference and measurement of these photons can entangle separate matter-based quantum systems as a resource for quantum information processing. Since confinement in a single-mode facilitates the photon interference needed for generating entanglement, the dipole emission patterns relevant in spontaneous emission present a mode-matching challenge. Current demonstrations rely on bulk photon-collection and manipulation optics that suffer from large component size and system-to-system variability—factors that impede scaling to the large numbers of entangled pairs needed for quantum information processing. To address these limitations, we demonstrate a collection method that enables passive phase stability, straightforward photonic manipulation, and intrinsic reproducibility. Specifically, we engineer a waveguide-integrated grating to couple photons emitted from a trapped ion into a single optical mode within a microfabricated ion-trap chip. Using the integrated collection optic, we characterize the collection efficiency, image the ion, and detect the ion’s quantum state. The integrated optic covers 2.18% of the solid angle and collects 1.97 ± 0.3% of the spontaneously emitted light incident on the grating for a total collection efficiency of 0.043% into a single-mode waveguide. This proof-of-principle demonstration lays the foundation for leveraging the inherent stability and reproducibility of integrated photonics to create, manipulate, and measure multipartite quantum states in arrays of quantum emitters.

Optical physics↗
DOE OSTI2023

Out-of-distribution generalization for learning quantum dynamics

Abstract Generalization bounds are a critical tool to assess the training data requirements of Quantum Machine Learning (QML). Recent work has established guarantees for in-distribution generalization of quantum neural networks (QNNs), where training and testing data are drawn from the same data distribution. However, there are currently no results on out-of-distribution generalization in QML, where we require a trained model to perform well even on data drawn from a different distribution to the training distribution. Here, we prove out-of-distribution generalization for the task of learning an unknown unitary. In particular, we show that one can learn the action of a unitary on entangled states having trained only product states. Since product states can be prepared using only single-qubit gates, this advances the prospects of learning quantum dynamics on near term quantum hardware, and further opens up new methods for both the classical and quantum compilation of quantum circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Practical Hamiltonian learning with unitary dynamics and Gibbs states

We study the problem of learning the parameters for the Hamiltonian of a quantum many-body system, given limited access to the system. In this work, we build upon recent approaches to Hamiltonian learning via derivative estimation. We propose a protocol that improves the scaling dependence of prior works, particularly with respect to parameters relating to the structure of the Hamiltonian (e.g., its locality k). Furthermore, by deriving exact bounds on the performance of our protocol, we are able to provide a precise numerical prescription for theoretically optimal settings of hyperparameters in our learning protocol, such as the maximum evolution time (when learning with unitary dynamics) or minimum temperature (when learning with Gibbs states). Thanks to these improvements, our protocol has practical scaling for large problems: we demonstrate this with a numerical simulation of our protocol on an 80-qubit system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exponential concentration in quantum kernel methods

Kernel methods in Quantum Machine Learning (QML) have recently gained significant attention as a potential candidate for achieving a quantum advantage in data analysis. Among other attractive properties, when training a kernel-based model one is guaranteed to find the optimal model’s parameters due to the convexity of the training landscape. However, this is based on the assumption that the quantum kernel can be efficiently obtained from quantum hardware. In this work we study the performance of quantum kernel models from the perspective of the resources needed to accurately estimate kernel values. We show that, under certain conditions, values of quantum kernels over different input data can be exponentially concentrated (in the number of qubits) towards some fixed value. Thus on training with a polynomial number of measurements, one ends up with a trivial model where the predictions on unseen inputs are independent of the input data. We identify four sources that can lead to concentration including expressivity of data embedding, global measurements, entanglement and noise. For each source, an associated concentration bound of quantum kernels is analytically derived. Lastly, we show that when dealing with classical data, training a parametrized data embedding with a kernel alignment method is also susceptible to exponential concentration. Our results are verified through numerical simulations for several QML tasks. Altogether, we provide guidelines indicating that certain features should be avoided to ensure the efficient evaluation of quantum kernels and so the performance of quantum kernel methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Large-scale simulations of Floquet physics on near-term quantum computers

Abstract Periodically driven quantum systems exhibit a diverse set of phenomena but are more challenging to simulate than their equilibrium counterparts. Here, we introduce the Quantum High-Frequency Floquet Simulation (QHiFFS) algorithm as a method to simulate fast-driven quantum systems on quantum hardware. Central to QHiFFS is the concept of a kick operator which transforms the system into a basis where the dynamics is governed by a time-independent effective Hamiltonian. This allows prior methods for time-independent simulation to be lifted to simulate Floquet systems. We use the periodically driven biaxial next-nearest neighbor Ising (BNNNI) model, a natural test bed for quantum frustrated magnetism and criticality, as a case study to illustrate our algorithm. We implemented a 20-qubit simulation of the driven two-dimensional BNNNI model on Quantinuum’s trapped ion quantum computer. Our error analysis shows that QHiFFS exhibits not only a cubic advantage in driving frequency ω but also a linear advantage in simulation time t compared to Trotterization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A measurement of the equation of state of carbon envelopes of white dwarfs

White dwarfs represent the final state of evolution for most stars. Certain classes of white dwarfs pulsate, leading to observable brightness variations, and analysis of these variations with theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution. However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. In this paper we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments, and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure–density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Logical quantum processor based on reconfigurable atom arrays

Suppressing errors is the central challenge for useful quantum computing, requiring quantum error correction (QEC) for large-scale processing. However, the overhead in the realization of error-corrected ‘logical’ qubits, in which information is encoded across many physical qubits for redundancy, poses substantial challenges to large-scale logical quantum computing. Here we report the realization of a programmable quantum processor based on encoded logical qubits operating with up to 280 physical qubits. Using logical-level control and a zoned architecture in reconfigurable neutral-atom arrays, our system combines high two-qubit gate fidelities, arbitrary connectivity, as well as fully programmable single-qubit rotations and mid-circuit readout. Operating this logical processor with various types of encoding, we demonstrate improvement of a two-qubit logic gate by scaling surface-code distance from d = 3 to d = 7, preparation of colour-code qubits with break-even fidelities, fault-tolerant creation of logical Greenberger–Horne–Zeilinger (GHZ) states and feedforward entanglement teleportation, as well as operation of 40 colour-code qubits. Finally, using 3D [[8,3,2]] code blocks, we realize computationally complex sampling circuits with up to 48 logical qubits entangled with hypercube connectivity with 228 logical two-qubit gates and 48 logical CCZ gates. We find that this logical encoding substantially improves algorithmic performance with error detection, outperforming physical-qubit fidelities at both cross-entropy benchmarking and quantum simulations of fast scrambling. These results herald the advent of early error-corrected quantum computation and chart a path towards large-scale logical processors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Parallel quantum annealing

Quantum annealers of D-Wave Systems, Inc., offer an efficient way to compute high quality solutions of NP-hard problems. This is done by mapping a problem onto the physical qubits of the quantum chip, from which a solution is obtained after quantum annealing. However, since the connectivity of the physical qubits on the chip is limited, a minor embedding of the problem structure onto the chip is required. In this process, and especially for smaller problems, many qubits will stay unused. We propose a novel method, called parallel quantum annealing, to make better use of available qubits, wherein either the same or several independent problems are solved in the same annealing cycle of a quantum annealer, assuming enough physical qubits are available to embed more than one problem. Although the individual solution quality may be slightly decreased when solving several problems in parallel (as opposed to solving each problem separately), we demonstrate that our method may give dramatic speed-ups in terms of the Time-To-Solution (TTS) metric for solving instances of the Maximum Clique problem when compared to solving each problem sequentially on the quantum annealer. Additionally, we show that solving a single Maximum Clique problem using parallel quantum annealing reduces the TTS significantly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

New proposed ITER divertor design using carbon insert on tungsten to mitigate ELMs and secondary radiation effects on nearby components

Abstract Building a successful device for the magnetic fusion energy production is a great challenge. ITER is an international project of the tokamak based magnetic fusion design being developed for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. A key obstacle to a successful magnetic fusion energy production is however, the performance during abnormal events including plasma disruptions and edge-localized modes (ELMs). A credible reactor design must tolerate at least a few of these transient events without serious consequences such as melting of the structure. This paper investigates and compares the performance of the current ITER tokamak design during two types of transient events, i.e., ELMs occurring at normal operation and disruptions during abnormal operation. We simulated the divertor components response using our integrated 3D HEIGHTS package. The simulations include self-consistent modeling of the interaction of the released core plasma particles with the initial solid divertor material, energy deposition processes, vaporization of divertor material, secondary plasma formation and MHD evolution, incident core particles collisions and scattering from this dense secondary plasma, photon radiation of secondary plasma, and the resulting heat loads on nearby components. Our simulations showed that using a small carbon insert around the strike point can significantly reduce the overall expected damage on the tungsten dome structure, reflector plates, and prevent tungsten vaporization and its potential core plasma contamination.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Quantum computing and preconditioners for hydrological linear systems

Modeling hydrological fracture networks is a hallmark challenge in computational earth sciences. Accurately predicting critical features of fracture systems, e.g. percolation, can require solving large linear systems far beyond current or future high performance capabilities. Quantum computers can theoretically bypass the memory and speed constraints faced by classical approaches, however several technical issues must first be addressed. Chief amongst these difficulties is that such systems are often ill-conditioned, i.e. small changes in the system can produce large changes in the solution, which can slow down the performance of linear solving algorithms. We test several existing quantum techniques to improve the condition number, but find they are insufficient. We then introduce the inverse Laplacian preconditioner, which improves the scaling of the condition number of the system from O(N) to O($\sqrt{N}$) and admits a quantum implementation. These results are a critical first step in developing a quantum solver for fracture systems, both advancing the state of hydrological modeling and providing a novel real-world application for quantum linear systems algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Plasma and trap-based techniques for science with antimatter

Positrons (i.e., antielectrons) find use in a wide variety of applications, and antiprotons are required for the formation and study of antihydrogen. Available sources of these antiparticles are relatively weak. To optimize their use, most applications require that the antiparticles be accumulated into carefully prepared plasmas. We present an overview of the techniques that have been developed to efficiently accumulate low energy antiparticles and create, in particular, tailored antiparticle plasmas. Techniques are also described to create tailored antiparticle beams. Many of these techniques are based on methods first developed by the nonneutral plasma community using electron plasmas for increased data rate. They have enabled the creation and trapping of antihydrogen, have been critical to studies of positron and positronium interactions with matter, including advanced techniques to characterize materials and material surfaces, and have led to the creation and study of the positronium molecule. Rather than attempting to be comprehensive, we focus on techniques that have proven most useful, applications where there has been significant, recent progress, and areas that hold promise for future advances. Examples of the latter include the ever more precise comparisons of the properties of antihydrogen and hydrogen, tests of gravity using antihydrogen and positronium atoms, and efforts to create and study phases of the many-electron, many-positron system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Flattening of the tokamak current profile by a fast magnetic reconnection with implications for the solar corona

During tokamak disruptions, the profile of the net parallel current is observed to flatten on a time scale that is so fast that it must be due to a fast magnetic reconnection. After a fast magnetic reconnection has broken magnetic surfaces, a single magnetic field line covers an entire volume and not just a magnetic surface. The current profile, given by K ≡ μ 0 j | | / B , relaxes to a constant within that volume by Alfven waves propagating along the chaotic magnetic field lines. The time scale for this relaxation determines the commonly observed disruption phenomena of a current spike and a sudden drop in the plasma internal inductance. An efficient method for studying this relaxation is derived, which allows a better understanding of the information encoded in the current spike and the associated sudden drop in the plasma internal inductance. Implications for coronal heating are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗