Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum simulations”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Limitations of Hartree–Fock with quantum resources

The Hartree–Fock problem provides the conceptual and mathematical underpinning of a large portion of quantum chemistry. As efforts in quantum technology aim to enhance computational chemistry algorithms, the Hartree–Fock method, central to many other numerical approaches, is a natural target for quantum enhanced algorithms. While quantum computers and quantum simulation offer many prospects for the future of modern chemistry, the non-deterministic polynomial-complete Hartree–Fock problem is not a likely candidate. Here, we highlight this fact from a number of perspectives including computational complexity, practical examples, and the full characterization of energy landscapes for simple systems.

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

Unleashed from constrained optimization: quantum computing for quantum chemistry employing generator coordinate inspired method

Hybrid quantum-classical approaches offer potential solutions to quantum chemistry problems, yet they often manifest as constrained optimization problems. Here, we explore the interconnection between constrained optimization and generalized eigenvalue problems through the Unitary Coupled Cluster (UCC) excitation generators. Inspired by the generator coordinate method, we employ these UCC excitation generators to construct non-orthogonal, overcomplete many-body bases, projecting the system Hamiltonian into an effective Hamiltonian, which bypasses issues such as barren plateaus that heuristic numerical minimizers often encountered in standard variational quantum eigensolver (VQE). Diverging from conventional quantum subspace expansion methods, we introduce an adaptive scheme that robustly constructs the many-body basis sets from a pool of the UCC excitation generators. This scheme supports the development of a hierarchical ADAPT quantum-classical strategy, enabling a balanced interplay between subspace expansion and ansatz optimization to address complex, strongly correlated quantum chemical systems cost-effectively, setting the stage for more advanced quantum simulations in chemistry.

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

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗
DOE OSTI2022

Quantum Phases of Transition Metal Dichalcogenide Moiré Systems

Moiré systems provide a rich platform for studies of strong correlation physics. Recent experiments on heterobilayer transition metal dichalcogenide Moiré systems are exciting in that they manifest a relatively simple model system of an extended Hubbard model on a triangular lattice. Inspired by the prospect of the hetero-transition metal dichalcogenide Moiré system’s potential as a solid-state-based quantum simulator, we explore the extended Hubbard model on the triangular lattice using the density matrix renormalization group. Specifically, we explore the two-dimensional phase space spanned by the key tuning parameters in the extended Hubbard model, namely, the kinetic energy strength and the further-range Coulomb interaction strengths. We find competition between Fermi fluid, chiral spin liquid, spin density wave, and charge order. In particular, our finding of the optimal further-range interaction for the chiral correlation presents a tantalizing possibility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A new “gold standard”: Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the “gold standard” CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.

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

Observation of a prethermal discrete time crystal

Extending the framework of statistical physics to the nonequilibrium setting has led to the discovery of previously unidentified phases of matter, often catalyzed by periodic driving. However, preventing the runaway heating that is associated with driving a strongly interacting quantum system remains a challenge in the investigation of these newly discovered phases. In this work, we utilize a trapped-ion quantum simulator to observe the signatures of a nonequilibrium driven phase without disorder—the prethermal discrete time crystal. Here, the heating problem is circumvented not by disorder-induced many-body localization, but rather by high-frequency driving, which leads to an expansive time window where nonequilibrium phases can emerge. Floquet prethermalization is thus presented as a general strategy for creating, stabilizing, and studying intrinsically out-of-equilibrium phases of matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Characterization of uncertainties in electron-argon collision cross sections

Abstract The predictive capability of a plasma discharge model depends on accurate representations of electron-impact collision cross sections, which determine the corresponding reaction rates and electron transport properties. The values of cross sections can be known only approximately either through experiments or simulations and are thus subject to uncertainties. Quantifying the uncertainties in plasma simulations allows us to assess the reliability of simulations and to provide a basis for interpreting discrepancies between simulations and experiments. For such uncertainty quantification of plasma simulations, it is essential to quantify the uncertainties of the underlying cross sections. Although much effort has been committed to calibrate the cross section values, their uncertainties are not well investigated. We characterize uncertainties in electron-argon atom collision cross sections using a Bayesian framework. Six collision processes—elastic momentum transfer, ionization, and four excitations—are characterized with semi-empirical models, which effectively capture the features important to the macroscopic properties of the plasma. A probability model for the uncertain parameters of these semi-empirical models is developed. Specifically, a Gaussian-process likelihood model is proposed to capture discrepancies among data sets, as well as the model-form inadequacies of the semi-empirical models. Two other likelihood models are compared with the proposed Gaussian-process model, to illustrate the importance of the choice of the likelihood model. The cross section models are calibrated using the electron-beam experiments and ab-inito quantum simulations. The resulting calibrated uncertainties capture well the scattering among the data sets. The calibrated cross section models are further validated against swarm-parameter experiments and zero-dimensional Boltzmann equation simulations of widely used cross section datasets.

Chung, Seung Whan (ORCID:0000000302501549)↗
DOE OSTI2021

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations

For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.

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

Continuous symmetry breaking in a two-dimensional Rydberg array

Spontaneous symmetry breaking underlies much of our classification of phases of matter and their associated transitions. The nature of the underlying symmetry being broken determines many of the qualitative properties of the phase; this is illustrated by the case of discrete versus continuous symmetry breaking. Indeed, in contrast to the discrete case, the breaking of a continuous symmetry leads to the emergence of gapless Goldstone modes controlling, for instance, the thermodynamic stability of the ordered phase. Here, we realize a two-dimensional dipolar XY model that shows a continuous spin-rotational symmetry using a programmable Rydberg quantum simulator. We demonstrate the adiabatic preparation of correlated low-temperature states of both the XY ferromagnet and the XY antiferromagnet. In the ferromagnetic case, we characterize the presence of a long-range XY order, a feature prohibited in the absence of long-range dipolar interaction. In conclusion, our exploration of the many-body physics of XY interactions complements recent works using the Rydberg-blockade mechanism to realize Ising-type interactions showing discrete spin rotation symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Absence of Off-Diagonal Long-Range Order in hcp 4 He Dislocation Cores

Here, the mass transport properties along dislocation cores in hcp 4 He are revisited by considering two types of edge dislocations as well as a screw dislocation, using a fully correlated quantum simulation approach. Specifically, we employ the zero-temperature path-integral ground state (PIGS) method together with ergodic sampling of the permutation space to investigate the fundamental dislocation core structures and their off-diagonal long-range order properties. It is found that the Bose-Einstein condensate fraction of such defective 4 He systems is practically null (≤ 10 –6 ), just as in the bulk defect-free crystal. These results provide compelling evidence for the absence of intrinsic superfluidity in dislocation cores in hcp 4 He and challenge the superfluid dislocation-network interpretation of the mass-flux-experiment observations, calling for further experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Finite-temperature tensor network study of the Hubbard model on an infinite square lattice

Here, the Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network—an infinite projected entangled pair state—evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With U⁡(1)×U⁡(1) symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators—measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit—provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Atomic-scale control of tunneling in donor-based devices

Atomically precise donor-based quantum devices are a promising candidate for solid-state quantum computing and analog quantum simulations. However, critical challenges in atomically precise fabrication have meant systematic, atomic scale control of the tunneling rates and tunnel coupling has not been demonstrated. Here using a room temperature grown locking layer and precise control over the entire fabrication process, we reduce unintentional dopant movement while achieving high quality epitaxy in scanning tunnelling microscope (STM)-patterned devices. Using the Si(100)2 × 1 surface reconstruction as an atomically-precise ruler to characterize the tunnel gap in precision-patterned single electron transistors, we demonstrate the exponential scaling of the tunneling resistance on the tunnel gap as it is varied from 7 dimer rows to 16 dimer rows. We demonstrate the capability to reproducibly pattern devices with atomic precision and a donor-based fabrication process where atomic scale changes in the patterned tunnel gap result in the expected changes in the tunneling rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Spin transport of a doped Mott insulator in moiré heterostructures

Moiré superlattices of semiconducting transition metal dichalcogenide heterobilayers are model systems for investigating strongly correlated electronic phenomena. Specifically, WSe 2 /WS 2 moiré superlattices have emerged as a quantum simulator for the two-dimensional extended Hubbard model. Experimental studies of charge transport have revealed correlated Mott insulator and generalized Wigner crystal states, but spin transport of the moiré heterostructure has not yet been sufficiently explored. Here, we use spatially and temporally resolved circular dichroism spectroscopy to directly image the spin transport as a function of carrier doping and temperature in WSe 2 /WS 2 moiré heterostructures. We observe diffusive spin transport at all hole concentrations at 11 Kelvin — including the Mott insulator at one hole per moiré unit cell — where charge transport is strongly suppressed. At elevated temperatures the spin diffusion constant remains unchanged in the Mott insulator state, but it increases significantly at finite doping away from the Mott state. The doping- and temperature-dependent spin transport can be qualitatively understood using a t–J model, where spins can move via the hopping of spin-carrying charges and via the exchange interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Two-dimensional spin systems in PECVD-grown diamond with tunable density and long coherence for enhanced quantum sensing and simulation

Systems of spins engineered with tunable density and reduced dimensionality enable a number of advancements in quantum sensing and simulation. Defects in diamond, such as nitrogen-vacancy (NV) centers and substitutional nitrogen (P1 centers), are particularly promising solid-state platforms to explore. However, the ability to controllably create coherent, two-dimensional spin systems and characterize their properties, such as density, depth confinement, and coherence, is an outstanding materials challenge. We present a refined approach to engineer dense (≳1 ppm ∙ nm), 2D nitrogen, and NV layers in diamond using delta-doping during plasma-enhanced chemical vapor deposition epitaxial growth. We employ both traditional materials techniques, e.g., secondary ion mass spectrometry, alongside NV spin decoherence-based measurements to characterize the density and dimensionality of the P1 and NV layers. We find P1 densities of 5–10 ppm ∙ nm, NV densities between 1 and 3.5 ppm ∙ nm tuned via electron irradiation dosage, and depth confinement of the spin layer down to 1.6 nm. We also observe high (up to 0.74) ratios of NV to P1 centers and reproducibly long NV coherence times, dominated by dipolar interactions with the engineered P1 and NV spin baths.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insights into the Mechanism of Chiral‐Induced Spin Selectivity: The Effect of Magnetic Field Direction and Temperature

Chiral oligopeptide monolayers are adsorbed on a ferromagnetic surface and their magnetoresistance is measured as a function of the angle between the magnetization of the ferromagnet and the surface normal. These measurements are conducted as a function of temperature for both enantiomers. The angle dependence is found to follow a changing trend with a period of 360°. Quantum simulations reveal that the angular distribution can be obtained only if the monolayer has significant effective spin orbit coupling (SOC), that includes contribution from the vibrations. The model shows that SOC only in the leads cannot reproduce the observed angular dependence. The simulation can reproduce the experiments if it included electron–phonon interactions and dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Efficient Quantum Gibbs Samplers with Kubo–Martin–Schwinger Detailed Balance Condition

Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Metasurface holographic optical traps for ultracold atoms

We propose metasurface holograms as a novel platform to generate optical trap arrays for cold atoms with high quality, efficiency, and thermal stability. We developed design and fabrication methods to create dielectric, phase-only metasurface holograms based on titanium dioxide. We experimentally demonstrated optical trap arrays of various geometries, including periodic and aperiodic configurations with dimensions ranging from 1D to 3D and up to a few hundred trap sites. We characterized the performance of the holographic metasurfaces in terms of the positioning accuracy, size and intensity uniformity of the generated traps, and power handling capability of the dielectric metasurfaces. Our proposed platform has great potential for enabling fundamental studies of quantum many-body physics, and quantum simulation and computation tasks. The compact form factor, passive nature, good power handling capability, and scalability of generating high-quality, large-scale arrays also make the metasurface platform uniquely suitable for realizing field-deployable devices and systems based on cold atoms.

42 ENGINEERING↗
DOE OSTI2021

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗