Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum dynamics 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Limitations for Quantum Algorithms to Solve Turbulent and Chaotic Systems

We investigate the limitations of quantum computers for solving nonlinear dynamical systems. In particular, we tighten the worst-case bounds of the quantum Carleman linearisation (QCL) algorithm answering one of their open questions. We provide a further significant limitation for any quantum algorithm that aims to output a quantum state that approximates the normalized solution vector. Given a natural choice of coordinates for a dynamical system with one or more positive Lyapunov exponents and solutions that grow sub-exponentially, we prove that any such algorithm has complexity scaling at least exponentially in the integration time. As such, an efficient quantum algorithm for simulating chaotic systems or regimes is likely not possible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Editorial: Quantum Information and Quantum Computing for Chemical Systems

Quantum computing has emerged as an exciting inter-disciplinary research topic that cuts across the traditional fields of physics, computer science, and engineering. It is a revolutionary model of computation that has offered new insights into methods for modeling and simulation of chemical systems. Applications of quantum computing to chemistry have demonstrated rapid progress on both theoretical and experimental fronts. Past theoretical efforts have shown how to adapt quantum computation to a variety of problems including electronic structure and molecular dynamics. In addition, the development of quantum algorithms for quantum chemistry has been stimulated by the greater availability of more capable quantum computing devices. Advances in the number and quality of qubits continues to enable remarkable proof-of-concept demonstrations working towards a milestone of quantum computational advantage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Short-Depth QAOA circuits and Quantum Annealing on Higher-Order Ising Models (Rev.2)

The Quantum Alternating Operator Ansatz (QAOA) and Quantum Annealing (QA) are quantum algorithms that are both based on the adiabatic theorem and both have the goal of sampling the optimal solution(s) of combinatorial optimization problems. Quantum annealing has been physically instantiated on D-Wave devices using superconducting flux qubits, and QAOA can be programmed on digital gate-model quantum computers such as the programmable superconducting transmon qubits devices of the IBMQ series, for instance ibm washington. QAOA and QA address the same types of problems, but it is unclear how they will scale to large problem sizes and to larger and higher-fidelity quantum computers. In this article, we present a direct comparison between QAOA, one and two rounds, run on all 127 qubits of ibm washington and QA run on D-Wave Advantage system4.1 and Advantage system6.1. The problems which allow for this comparison are random Ising model problems whose connectivity matches the heavy hexagonal lattice topology of ibm washington and the Pegasus graph connectivity of the two D-Wave devices. We create two classes of problem instances for this comparison: one with higher order terms (ZZZ variable interactions), linear terms, and quadratic terms, and a separate problem type with only linear and quadratic terms. Our QAOA circuits are novel and extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of ibm washington’s heavy hexagonal lattice and can be applied to future heavy-hex chips. We also test the effectiveness of the error suppression technique digital dynamical decoupling on the QAOA circuits. The QAOA circuits compiled to ibm washington are composed of several thousand circuit instructions, approximately 3, 000 depending on the details of the circuit, making these some the largest quantum circuits ever executed on a digital quantum processor. QAOA and QA are compared against the classical heuristic algorithm of simulated annealing and all problem instances are exactly solved using CPLEX in order to evaluate which samplers, if any, correctly found the ground state solution(s) of the problem instances. We find that (i) QA outperforms QAOA on all problem instances, (ii) QAOA samples the problems better than random sampling, and (iii) QAOA angle computation exhibits clear parameter concentration across the ensemble of Ising models.

127 Qubits↗
DOE OSTI2020

Nonadiabatic Electronic Energy Transfer in Chemical Oxygen-Iodine Laser: Powered by Derivative Coupling or Spin-Orbit Coupling?

Derivative couplings near a conical intersection and spin-orbit couplings between different spin states are known to facilitate nonadiabatic transitions in molecular systems. In this report, we investigate a prototypical electronic energy transfer process, I( 2 P 3/2 )+O 2 (a 1 Δ g )→I( 2 P 1/2 )+ O 2 (X 3 Σ g - ), which is of great importance for the chemical oxygen-iodine laser. To understand the nonadiabatic dynamics, this multistate process is investigated in full dimensionality with quantum wave packets using diabatic potential energy surfaces coupled by both derivative and spin-orbit couplings, all determined from first principles. A near quantitative agreement with structural, energetic, and kinetic measurements is achieved. Detailed analyses suggest that the nonadiabatic dynamics is largely controlled by derivative coupling near conical intersections, which leads to a small effective barrier and hence a slightly positive temperature dependence of the rate coefficient. The new findings should extend our understanding of energy transfer, provide a quantitative basis for numerical simulations of the chemical oxygen-iodine laser, and have important implications in other electronic energy transfer processes.

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

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

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

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

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

Machine-Learning Accelerated Studies of Materials with High Performance and Edge Computing

In the studies of materials, experimental measurements often serve as the reference to verify physics theory and modeling; while theory and modeling provide a fundamental understanding of the physics and principles behind. However, the interactions and cross validation between them have long been a challenge even to-date. Not only that inferring a physics model from experimental data is itself a difficult inverse problem, another major challenge is the orders-of-magnitude longer wall-clock time required to carry out high-fidelity computer modeling to match the timescale of experiments. We envisage that by combining high performance computing, data science, and edge computing technology, the current predicament can be alleviated, and a new paradigm of data-driven physics research will open up. For example, we can accelerate computer simulations by first performing the large-scale modeling on high performance computers and train a machine-learned surrogate model. This computationally inexpensive surrogate model can then be transferred to the computing units residing closely to the experimental facilities to perform high-fidelity simulations at a much higher throughout. The model will also be more amenable to analyzing and validating experimental observations in comparable time scales at a much lower computational cost. Further integration of these accelerated computer simulations with an outer machine learning loop can also inform and direct future experiments, while making the inverse problem of physics model inference more tractable. We will demonstrate a proof-of-concept by using a quantum Monte Carlo application, Dynamical Cluster Approximation (DCA++), to machine-learn a surrogate model and accelerate the study of quantum correlated materials.

Li, Ying Wai↗
DOE OSTI2024

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗
DOE OSTI2024

Two-dimensional-lattice-confined single-molecule-like aggregates

Intermolecular distance largely determines the optoelectronic properties of organic matter. Conventional organic luminescent molecules are commonly used either as aggregates or as single molecules that are diluted in a foreigner matrix. They have garnered great research interest in recent decades for a variety of applications, including light-emitting diodes, lasers and quantum technologies, among others. However, there is still a knowledge gap on how these molecules behave between the aggregation and dilution states. Here we report an unprecedented phase of molecular aggregate that forms in a two-dimensional hybrid perovskite superlattice with a near-equilibrium distance, which we refer to as a single-molecule-like aggregate (SMA). By implementing two-dimensional superlattices, the organic emitters are held in proximity, but, surprisingly, remain electronically isolated, thereby resulting in a near-unity photoluminescence quantum yield, akin to that of single molecules. Moreover, the emitters within the perovskite superlattices demonstrate strong alignment and dense packing resembling aggregates, allowing for the observation of robust directional emission, substantially enhanced radiative recombination and efficient lasing. Molecular dynamics simulations together with single-crystal structure analysis emphasize the critical role of the internal rotational and vibrational degrees of freedom of the molecules in the two-dimensional lattice for creating the exclusive SMA phase. Furthermore, this two-dimensional superlattice unifies the paradoxical properties of single molecules and aggregates, thus offering exciting possibilities for advanced spectroscopic and photonic applications.

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

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

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

Practical guide to replica exchange transition interface sampling and forward flux sampling

Path sampling approaches have become invaluable tools to explore the mechanisms and dynamics of the so-called rare events that are characterized by transitions between metastable states separated by sizable free energy barriers. Their practical application, in particular to ever more complex molecular systems, is, however, not entirely trivial. Here, focusing on replica exchange transition interface sampling (RETIS) and forward flux sampling (FFS), we discuss a range of analysis tools that can be used to assess the quality and convergence of such simulations, which is crucial to obtain reliable results. The basic ideas of a step-wise evaluation are exemplified for the study of nucleation in several systems with different complexities, providing a general guide for the critical assessment of RETIS and FFS simulations.

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

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantum optimal control of superconducting qubits based on machine-learning characterization

Implementing fast and high-fidelity quantum operations using open-loop quantum optimal control relies on having an accurate model of the quantum dynamics. Any deviations between this model and the complete dynamics of the device, such as the presence of spurious modes or pulse distortions, can degrade the performance of optimal controls in practice. Here, we propose an experimentally simple approach to realize optimal quantum controls tailored to the device parameters and environment while specifically characterizing this quantum system. Concretely, we use physics-inspired machine learning to infer an accurate model of the dynamics from experimentally available data and then optimize our experimental controls on this trained model. We show the power and feasibility of this approach by optimizing arbitrary single-qubit operations in detailed numerical simulations of a superconducting transmon qubit. Furthermore, we demonstrate that this framework produces an accurate description of the device dynamics under arbitrary controls, together with the precise pulses achieving arbitrary single-qubit gates with a high fidelity of ∼99.99%.

Artificial neural networks↗
DOE OSTI2023

Circumventing the phonon bottleneck by multiphonon-mediated hot exciton cooling at the nanoscale

Abstract Nonradiative processes govern efficiencies of semiconductor nanocrystal (NC)-based devices. A central process is hot exciton cooling, or the nonradiative relaxation of a highly excited electron/hole pair to form a band-edge exciton. Due to quantum confinement effects, the timescale and mechanism of cooling are not well understood. A mismatch between electronic energy gaps and phonon frequencies has led to the hypothesis of a phonon bottleneck and extremely slow cooling, while enhanced electron-hole interactions have suggested ultrafast cooling. Experimental measurements of the cooling timescale range six orders of magnitude. Here, we develop an atomistic approach to describe phonon-mediated exciton dynamics and simulate cooling in NCs of experimentally relevant sizes. We find that cooling occurs on ~30 fs timescales in CdSe NCs, in agreement with the most recent measurements, and that the phonon bottleneck is circumvented through a cascade of multiphonon-mediated relaxation events. Furthermore, we identify NC handles for tuning the cooling timescale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

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

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Energy and momentum conservation in spin transfer

We utilize simulations of spin-polarized electron scattering by a chain of localized quantum spins to show that energy and linear momentum conservation laws impose strong constraints on the properties of magnetic excitations induced by spin transfer. In turn, orbital and spin dynamics of conduction electrons depends on the dynamical characteristics of the local spins. Furthermore, our results suggest the possibility to achieve precise control of spin-transfer-driven magnetization dynamics by tailoring the spectral characteristics of the magnetic systems and of the driving electrons.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Braiding of Majorana zero modes in vortex cores

Here, we demonstrate the successful simulation of $\sqrt{𝑍-}$, $\sqrt{𝑋-}$, and 𝑋-quantum gates using Majorana zero modes (MZMs) that emerge in magnetic vortices located in topological superconductors. We compute the transition probabilities and geometric phase differences accounting for the full many-body dynamics and show that qubit states can be read out by fusing the vortex core MZMs and measuring the resulting charge density. We visualize the gate processes using the time- and energy-dependent nonequilibrium local density of states. Our results demonstrate the feasibility of employing vortex core MZMs for the realization of fault-tolerant topological quantum computing.

Majorana bound states↗