Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “classical simulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Three-flavor collective neutrino oscillation simulations on a qubit quantum annealer

Neutrinos are unique among elementary particles in that their flavor-compositions oscillate over time. In extreme environments such as core-collapse supernovae, neutron-star mergers, and the early Universe, neutrinos are dense enough that their self-interactions significantly affect, if not dominate, these oscillations. This has implications for several phenomena within these environments, particularly nucleosynthesis. Simulations of these self-interactions have traditionally approximated neutrinos as having two flavors instead of the physical three. In order to develop techniques for characterizing the resulting quantum entanglement, I present the results of simulations of neutrino-neutrino interactions that include all three physical neutrino flavors and were performed on D-Wave Inc.’s Advantage 5000+ qubit quantum annealer. These results are checked against those from exact classical simulations, which are also used to compare the neutrino-neutrino interactions to neutrino-antineutrino and interactions between Majorana neutrinos, which are their own antiparticles. The D-Wave Advantage annealer is shown to be able to reproduce time evolution with the precision of a classical machine for small numbers of neutrinos and to do so without the Trotter errors present in most simulations of dynamics on quantum devices. Furthermore, it suffers from poor scaling in qubit-count with the number of neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗
DOE OSTI2025

Delocalized Excitation Transfer in Open Quantum Systems with Long-Range Interactions

The interplay between coherence and system-environment interactions is at the basis of a wide range of phenomena, from quantum information processing to charge and energy transfer in molecular systems, biomolecules, and photochemical materials. In this work, we use a Frenkel exciton model with long-range interacting qubits coupled to a damped collective bosonic mode to investigate vibrationally assisted transfer processes in donor-acceptor systems featuring internal substructures analogous to light-harvesting complexes. We find that certain delocalized excitonic states maximize the transfer rate and that the entanglement is preserved during the dissipative transfer over a wide range of parameters. We investigate the reduction in transfer caused by static disorder, white noise, and finite temperature and study how transfer efficiency scales as a function of the number of dimerized monomers and the component number of each monomer, finding which excitonic states lead to optimal transfer. Finally, we provide a realistic experimental setting to realize this model in analog trapped-ion quantum simulators. Analog quantum simulation of systems comprising many and increasingly complex monomers could offer valuable insights into the design of light-harvesting materials, particularly in the nonperturbative intermediate parameter regime examined in this study, where classical simulation methods are resource intensive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Hybrid quantum simulations with qubits and qumodes on trapped-ion platforms

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unex- plored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of free- dom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Furthermore, our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Lower-dimensional field theories↗
DOE OSTI2025

Steps toward quantum simulations of hadronization and energy loss in dense matter

A framework for simulating the real-time dynamics of composite particles in a simple model of dense matter that is amenable to quantum computers is developed. As a demonstration, we perform classical simulations of heavy-hadrons propagating through a dense medium in the Schwinger model. Measurements of the time-dependent energy and charge density are used to identify mechanisms responsible for energy loss and hadron production (hadronization). A study of entanglement dynamics highlights the importance of quantum coherence between the particles that make up the dense medium. Throughout this work, care is taken to isolate, and remove, phenomena that arise solely from a finite lattice spacing. It is found that signatures of entanglement are more sensitive to lattice artifacts than other observables. Toward quantum simulations, we present an efficient method and the corresponding quantum circuits for preparing ground states in the presence of heavy mesons. Finally, these circuits are used to estimate the resources required to simulate in-medium energy loss and hadronization in the Schwinger model using quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Nuclear Quantum Effects on the Electronic Structure of Water and Ice

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

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

State preparation of lattice field theories using quantum optimal control

Here, we explore the application of quantum optimal control (QOC) techniques to state preparation of lattice field theories on quantum computers. As a first example, we focus on the Schwinger model, quantum electrodynamics in 1+1 dimensions. We demonstrate that QOC can significantly speed up the ground state preparation compared to gate-based methods, even for models with long-range interactions. Using classical simulations, we explore the dependence on the interqubit coupling strength and the device connectivity, and we study the optimization in the presence of noise. While our simulations indicate potential speedups, the results strongly depend on the device specifications. In addition, we perform exploratory studies on the preparation of thermal states. Our results motivate further studies of QOC techniques in the context of quantum simulations for fundamental physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Grover-QAOA for 3-SAT: quadratic speedup, fair-sampling, and parameter clustering

Abstract The SAT problem is a prototypical NP-complete problem of fundamental importance in computational complexity theory with many applications in science and engineering; as such, it has long served as an essential benchmark for classical and quantum algorithms. This study shows numerical evidence for a quadratic speedup of the Grover Quantum Approximate Optimization Algorithm (G-QAOA) over random sampling for finding all solutions to 3-SAT (All-SAT) and Max-SAT problems. G-QAOA is less resource-intensive and more adaptable for these problems than Grover’s algorithm, and it surpasses conventional QAOA in its ability to sample all solutions. We show these benefits by classical simulations of many-round G-QAOA on thousands of random 3-SAT instances. We also observe G-QAOA advantages on the IonQ Aria quantum computer for small instances, finding that current hardware suffices to determine and sample all solutions. Interestingly, a single-angle-pair constraint that uses the same pair of angles at each G-QAOA round greatly reduces the classical computational overhead of optimizing the G-QAOA angles while preserving its quadratic speedup. We also find parameter clustering of the angles. The single-angle-pair protocol and parameter clustering significantly reduce obstacles to classical optimization of the G-QAOA angles.

Zhang, Zewen (ORCID:000000032258613X)↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2026

Quantum utility-scale error mitigation for quantum quench dynamics in Heisenberg spin chains

Here, we implement a quantum error mitigation method termed self-mitigation, which is comparable to zero-noise extrapolation, at large scales to achieve quantum utility on near-term, noisy quantum computers. We investigate the effectiveness of several quantum error mitigation strategies, including self-mitigation, by simulating quantum quench dynamics for Heisenberg spin chains with system sizes up to 104 qubits using IBM quantum processors. In particular, we discuss the limitations of zero-noise extrapolation and the advantages offered by self-mitigation at large scales. The self-mitigation method demonstrates stable accuracy with large systems of 104 qubits comprising more than 3,000 CNOT gates. Also, we combine the discussed quantum error mitigation methods with practical entanglement entropy measuring methods, and it shows a good agreement with the theoretical estimation. Our study illustrates the usefulness of near-term noisy quantum hardware in examining the quantum quench dynamics of many-body systems at large scales and lays the groundwork for surpassing classical simulations with quantum methods prior to the development of fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Analyzing the Quantum Approximate Optimization Algorithm: Ansätze, Symmetries, and Lie Algebras

The quantum approximate optimization algorithm (QAOA) has been proposed as a method to obtain approximate solutions for combinatorial optimization tasks. In this work, we study the underlying algebraic properties of three QAOA ansätze for the maximum-cut problem on connected graphs, while focusing on the generated Lie algebras as well as their invariant subspaces. Specifically, we analyze the standard QAOA ansatz as well as the orbit and multiangle ansätze. We are able to fully characterize the Lie algebras of the multiangle ansatz across arbitrary connected graphs, finding that they only fall into one of just six families. Aside from the cycle and path graphs, the Lie dimensions for every graph are exponentially large in the system size, meaning that multiangle ansätze are extremely prone to exhibiting barren plateaus. Then, a similar quasi-graph-independent Lie-algebraic characterization beyond the multiangle ansatz is impeded as the circuit exhibits additional “hidden” symmetries besides those naturally arising from a certain parity-superselection operator and all automorphisms of the considered graph. Disregarding the “hidden” symmetries, we can upper bound the dimensions of the orbit and the standard Lie algebras, and the dimensions of the associated invariant subspaces are determined via explicit character formulas. To finish, we conjecture that (for most graphs) the standard Lie algebras have only components that are either exponential or that grow, at most, polynomially with the system size. This would imply that the QAOA is either prone to barren plateaus or classically simulable. More generally, our work provides a symmetry framework and tools to analyze any desired variational quantum algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2025

Dissolved CO 2 Modulates the Electrochemical Capacitance on Gold Electrodes

The presence of CO 2 at an electrified interface between an aqueous electrolyte and a metal electrode is the prerequisite for many electrochemical CO 2 capture technologies. To understand the behavior of dissolved CO 2 at an aqueous electrified interface, we characterized the electrochemical interface of planar gold electrodes with cyclic voltammetry, electrochemical impedance spectroscopy (EIS), electrochemical surface plasmon resonance (EC–SPR), and attenuated total reflectance surface-enhanced infrared spectroscopy (ATR–SEIRAS). Under all investigated conditions, we observed a decrease in the electrochemical capacitance upon saturation of the electrolyte with CO 2 , as compared to an electrolyte saturated with Ar. EIS and EC–SPR showed that this capacitance reduction was also potential dependent: it reached a minimum near the point of zero charge and became more significant as the applied potential moved further away from the point of zero charge. Hybrid quantum–classical simulations of the gold/aqueous electrolyte interface indicate that bicarbonate decreases the capacitance and modifies the composition of the electric double layer. In addition to the binding of bicarbonate under positive bias, we propose that molecular CO 2 can be induced by applied potential to concentrate in the diffuse layer of the electric double layer, leading to a reduction in the electrochemical capacitance under both negative and positive bias. Furthermore, this work advances the understanding of non-Faradaic effects of dissolved CO 2 at aqueous electrified interfaces of relevance for electrochemical CO 2 capture.

25 ENERGY STORAGE↗
DOE OSTI2025

ON-OFF neuromorphic ISING machines using Fowler-Nordheim annealers

We introduce NeuroSA, a neuromorphic architecture specifically designed to ensure asymptotic convergence to the ground state of an Ising problem using a Fowler-Nordheim quantum mechanical tunneling based threshold-annealing process. The core component of NeuroSA consists of a pair of asynchronous ON-OFF neurons, which effectively map classical simulated annealing dynamics onto a network of integrate-and-fire neurons. The threshold of each ON-OFF neuron pair is adaptively adjusted by an FN annealer and the resulting spiking dynamics replicates the optimal escape mechanism and convergence of SA, particularly at low-temperatures. To validate the effectiveness of our neuromorphic Ising machine, we systematically solved benchmark combinatorial optimization problems such as MAX-CUT and Max Independent Set. Across multiple runs, NeuroSA consistently generates distribution of solutions that are concentrated around the state-of-the-art results (within 99%) or surpass the current state-of-the-art solutions for Max Independent Set benchmarks. Furthermore, NeuroSA is able to achieve these superior distributions without any graph-specific hyperparameter tuning. For practical illustration, we present results from an implementation of NeuroSA on the SpiNNaker2 platform, highlighting the feasibility of mapping our proposed architecture onto a standard neuromorphic accelerator platform.

42 ENGINEERING↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗
DOE OSTI2025

Exploring transport mechanisms in atomic precision advanced manufacturing enabled pn junctions

We investigate the different transport mechanisms that can occur in pn junction devices made using atomic precision advanced manufacturing at temperatures ranging from cryogenic to room temperature. We first elucidate the potential cause of the anomalous behavior observed in the forward-bias response of these devices in recent cryogenic temperature measurements, which deviates from the theoretical response of a silicon Esaki diode. These anomalous behaviors include current suppression at low voltages in the forward-bias response and a much lower valley voltage at cryogenic temperatures than theoretically expected for a silicon diode. To investigate the potential causes of these anomalies, we studied the effects of a few possible transport mechanisms, including band-to-band tunneling, bandgap narrowing, potential impact of non-Ohmic contacts, band quantization, impact of leakage, and inelastic trap-assisted tunneling, through semi-classical simulations. We find that a combination of two sets of band-to-band tunneling (BTBT) parameters can qualitatively approximate the shape of the tunneling current at low bias. This can arise from band quantization and realignment due to the strong potential confinement in δ-layers. We also find that the lower-than-theoretically-expected valley voltage can be attributed to modifications in the electronic band structure within the δ-layer regions, leading to a significant bandgap narrowing induced by the high density of dopants. Finally, we extend our analyses to room temperature operation and predict that trap-assisted tunneling (TAT) facilitated by phonon interactions may become significant, leading to a complex superposition of BTBT and TAT transport mechanisms in the electrical measurements.

band gap↗