Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum computing for 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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Electronic Structure of InAs and InSb Surfaces: Density Functional Theory and Angle-Resolved Photoemission Spectroscopy

The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, in this work, the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces is studied using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). Large-scale first principles simulations are enabled by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, a “bulk unfolding” scheme is implemented by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, a good agreement is found between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states, which may be detected by ARPES with low photon energies. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. The combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, for example, topological qubits utilizing the Majorana zero modes.

42 ENGINEERING↗
DOE OSTI2015

Structure and dynamics of H{sub 2}O vis-á-vis phenylalanine recognition at a DPPC lipid membrane via interfacial H-bond types: Insights from polarized FT-IRRAS and ADMP simulations

Preferential and enantioselective interactions of L-/D-Phenylalanine (L-Phe and D-Phe) and butoxycarbonyl-protected L-/D-Phenylalanine (LPA and DPA) as guest with 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (L-DPPC) as host were tapped by using real time Fourier transform infrared reflection absorption spectroscopy (FT-IRRAS). Polarization-modulated FT-IRRAS of DPPC monolayers above the phenylalanine modified subphases depicted fine structure/conformation differences under considerations of controlled 2D surface pressure. Selective molecular recognition of D-enantiomer over L-enantiomer driven by the DPPC head group via H-bonding and electrostatic interactions was evident spectroscopically. Accordingly, binding constants (K) of 145, 346, 28, and 56 M{sup −1} for LPA, DPA, L-Phe, and D-Phe, respectively, were estimated. The real time FT-IRRAS water bands were strictly conformation sensitive. The effect of micro-solvation on the structure and stability of the 1:1 diastereomeric L-lipid⋯, LPA/DPA and L-lipid⋯, (L/D)-Phe adducts was investigated with the aid of Atom-centered Density Matrix Propagation (ADMP), a first principle quantum mechanical molecular dynamics approach. The phosphodiester fragment was the primary site of hydration where specific solvent interactions were simulated through single- and triple- “water-phosphate” interactions, as water cluster’s “tetrahedral dice” to a “trimeric motif” transformation as a partial de-clusterization was evident. Under all the hydration patterns considered in both static and dynamic descriptions of density functional theory, L-lipid/D-amino acid enantiomer adducts continued to be stable structures while in dynamic systems, water rearranged without getting “squeezed-out” in the process of recognition. In spite of the challenging computational realm of this multiscale problem, the ADMP simulated molecular interactions complying with polarized vibrational spectroscopy unraveled a novel route to chiral recognition and interfacial water structure.

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

Potential Applications of Quantum Computing at Los Alamos National Laboratory, v0.3.0

Since the scientific revolution in the 16th and 17th centuries, the process of scientific discovery has followed an iterative feedback process of observation, hypothesis development and testing with physical experiments, which is widely referred to as the scientific method. This process remained largely unchanged until the middle of the 20th century, when the emergence of digital computers empowered scientist to build and inspect detailed simulations of physical phenomena. Over the last century, computational tools have transformed modern approaches to scientific discovery by enabling fast and affordable hypothesis testing before physical experiments are conducted, shown in Figure 1-1. Some notable examples include: global climate forecasts to understand how the environment may change over decades [130]; modeling the behavior of plasma to design fusion reactors [59]; and understanding the behavior of molecules in biological processes [161, 223].

36 MATERIALS SCIENCE↗
DOE OSTI2020

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

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

Simulation of injector dynamics during steady inductive helicity injection current drive in the HIT-SI experiment

We present simulations of inductive helicity injection in the Helicity Injected Torus with Steady Inductive helicity injection (HIT-SI) device that treats the entire plasma volume in a single dynamic MHD model. A new fully 3D numerical tool, the PSI-center TETrahedral mesh code, was developed that provides the geometric flexibility required for this investigation. Implementation of a zero-β Hall MHD model using PSI-TET will be presented including formulation of a new self-consistent magnetic boundary condition for the wall of the HIT-SI device. Results from simulations of HIT-SI are presented focusing on injector dynamics that are investigated numerically for the first time. Asymmetries in the plasma loading between the two helicity injectors and progression of field reversal in each injector are observed. Analysis indicates cross-coupling between injectors through confinement volume structures. Injector impedance is found to scale with toroidal current at fixed density, consistent with experimental observation. Comparison to experimental data with an injector drive frequency of 14.5 kHz shows good agreement with magnetic diagnostics. Global mode structures from Bi-Orthogonal decomposition agree well with experimental data for the first four modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Solving independent set problems with photonic quantum circuits

An independent set (IS) is a set of vertices in a graph such that no edge connects any two vertices. In adiabatic quantum computation [E. Farhi, et al ., Science 292, 472–475 (2001); A. Das, B. K. Chakrabarti, Rev. Mod. Phys. 80, 1061–1081 (2008)], a given graph G ( V , E ) can be naturally mapped onto a many-body Hamiltonian H IS G ( V , E ) , with edges E being the two-body interactions between adjacent vertices V . Thus, solving the IS problem is equivalent to finding all the computational basis ground states of H IS G ( V , E ) . Very recently, non-Abelian adiabatic mixing (NAAM) has been proposed to address this task, exploiting an emergent non-Abelian gauge symmetry of H IS G ( V , E ) [B. Wu, H. Yu, F. Wilczek, Phys. Rev. A 101, 012318 (2020)]. Here, we solve a representative IS problem G ( 8 , 7 ) by simulating the NAAM digitally using a linear optical quantum network, consisting of three C-Phase gates, four deterministic two-qubit gate arrays (DGA), and ten single rotation gates. The maximum IS has been successfully identified with sufficient Trotterization steps and a carefully chosen evolution path. Remarkably, we find IS with a total probability of 0.875(16), among which the nontrivial ones have a considerable weight of about 31.4%. Our experiment demonstrates the potential advantage of NAAM for solving IS-equivalent problems.

Science & Technology - Other Topics↗
DOE OSTI2015

Spin-flip, tensor equation-of-motion configuration interaction with a density-functional correction: A spin-complete method for exploring excited-state potential energy surfaces

We revisit the formalism of the spin-adapted, spin-flip (SA-SF) configuration-interaction singles (CIS) method based on a tensor equation-of-motion formalism that affords proper spin eigenstates without sacrificing single-reference simplicity. Matrix elements for SA-SF-CIS are then modified in a manner similar to collinear spin-flip time-dependent density functional theory (SF-TDDFT), to include a DFT exchange-correlation correction. The performance of this method, which we call SA-SF-DFT, is evaluated numerically and we find that it systematically improves the energies of electronic states that exhibit significant spin contamination within the conventional SF-TDDFT approach. The new method cures the state assignment problem that plagues geometry optimizations and ab initio molecular dynamics simulations using traditional SF-TDDFT, without sacrificing computational efficiency, and furthermore provides correct topology at conical intersections, including those that involve the ground state, unlike conventional TDDFT. As such, SA-SF-DFT appears to be a promising method for generating excited-state potential energy surfaces at DFT cost.

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

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

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

Parallel hybrid quantum-classical machine learning for kernelized time-series classification

Supervised time-series classification garners widespread interest because of its applicability throughout a broad application domain including finance, astronomy, biosensors, and many others. Here, in this work, we tackle this problem with hybrid quantum-classical machine learning, deducing pairwise temporal relationships between time-series instances using a timeseries Hamiltonian kernel (TSHK). A TSHK is constructed with a sum of inner products generated by quantum states evolved using a parameterized time evolution operator. This sum is then optimally weighted using techniques derived from multiple kernel learning. Because we treat the kernel weighting step as a differentiable convex optimization problem, our method can be regarded as an end-to-end learnable hybrid quantum-classical-convex neural network, or QCC-net, whose output is a data set-generalized kernel function suitable for use in any kernelized machine learning technique such as the support vector machine (SVM). Using our TSHK as input to a SVM, we classify univariate and multivariate time-series using quantum circuit simulators and demonstrate the efficient parallel deployment of the algorithm to 127-qubit superconducting quantum processors using quantum multi-programming.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

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

A quantitative quantum-chemical analysis tool for the distribution of mechanical force in molecules

The promising field of mechanochemistry suffers from a general lack of understanding of the distribution and propagation of force in a stretched molecule, which limits its applicability up to the present day. In this article, we introduce the JEDI (Judgement of Energy DIstribution) analysis, which is the first quantum chemical method that provides a quantitative understanding of the distribution of mechanical stress energy among all degrees of freedom in a molecule. The method is carried out on the basis of static or dynamic calculations under the influence of an external force and makes use of a Hessian matrix in redundant internal coordinates (bond lengths, bond angles, and dihedral angles), so that all relevant degrees of freedom of a molecule are included and mechanochemical processes can be interpreted in a chemically intuitive way. The JEDI method is characterized by its modest computational effort, with the calculation of the Hessian being the rate-determining step, and delivers, except for the harmonic approximation, exact ab initio results. We apply the JEDI analysis to several example molecules in both static quantum chemical calculations and Born-Oppenheimer Molecular Dynamics simulations in which molecules are subject to an external force, thus studying not only the distribution and the propagation of strain in mechanically deformed systems, but also gaining valuable insights into the mechanochemically induced isomerization of trans-3,4-dimethylcyclobutene to trans,trans-2,4-hexadiene. The JEDI analysis can potentially be used in the discussion of sonochemical reactions, molecular motors, mechanophores, and photoswitches as well as in the development of molecular force probes.

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

Theoretical modeling of UV-Vis absorption and emission spectra in liquid state systems including vibrational and conformational effects: Explicit treatment of the vibronic transitions

Here, we extend a recently introduced theoretical-computational procedure [M. D’Alessandro, M. Aschi, C. Mazzuca, A. Palleschi, and A. Amadei, J. Chem. Phys. 139, 114102 (2013)] to include quantum vibrational transitions in modelling electronic spectra of atomic molecular systems in condensed phase. The method is based on the combination of Molecular Dynamics simulations and quantum chemical calculations within the Perturbed Matrix Method approach. The main aim of the presented approach is to reproduce as much as possible the spectral line shape which results from a subtle combination of environmental and intrinsic (chromophore) mechanical-dynamical features. As a case study, we were able to model the low energy UV-vis transitions of pyrene in liquid acetonitrile in good agreement with the experimental data.

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

Symmetry Breaking Bias and the Dynamics of a Quantum Phase Transition

The Kibble-Zurek mechanism predicts the formation of topological defects and other excitations that quantify how much a quantum system driven across a quantum critical point fails to be adiabatic. We point out that, thanks to the divergent linear susceptibility at the critical point, even a tiny symmetry breaking bias can restore the adiabaticity. Furthermore, the minimal required bias scales like τ Q - β δ / ( 1 + z ν ) , where β , δ , z , ν are the critical exponents and τ Q is a quench time. We test this prediction by DMRG simulations of the quantum Ising chain. It is directly applicable to the recent emulation of quantum phase transition dynamics in the Ising chain with ultracold Rydberg atoms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Implementing a Laser Stabilization System for Trapping Ca+ Ions: an Internship Reflection

At Lawrence Livermore National Laboratory, I contributed to a project developing 3D printed micro ion traps for quantum computing. I designed, implemented, and assessed a laser stabilization system that locked lasers to the frequencies required for calibrating our High Finesse WS8-10 wavelength meter and for laser cooling and trapping of Ca+ ions. I also programmed a Python interface for hardware communication, data collection, and statistical analysis. Additionally, I optimized and aligned laser beam paths, and I implemented a closed digital feedback loop using Proportional, Integral, and Derivative (PID) control parameters. I analyzed both the long-term and short-term behavior of our locked lasers and adjusted PID parameters to enhance performance. Furthermore, I used COMSOL to simulate the capacitance of a linear Paul trap design and predict our trap’s performance. The procedures I developed for the interface, analysis, and simulations will continue to support the ion trapping experiment after my appointment. I strengthened my skills in data analysis, Python coding, and optical alignment for laser systems. My confidence as a researcher grew, particularly in communicating my research. This experience taught me the importance of careful planning and consideration in research and solidified my desire to continue exploring novel quantum technology as an undergraduate

42 ENGINEERING↗
DOE OSTI2019

Effective matrix model for nuclear physics on a quantum computer

In this project we studied a simplified potential from [1] describing nuclear physics from quark and gluon fields, where spatial degrees of freedom have been integrated over. This potential deals with the phase structure of the strong interaction. We created a matrix representation of the Hamiltonian for this theory using discrete quantum mechanics [2]. We then feed this Hamiltonian to IBM's Qiskit package [3]. We run on IBM's 32-qubit quantum simulator, which approximates the ground state energy using the variational quantum eigensolver (VQE) algorithm [4]. This result is then compared with the exact ground state energy, which we can calculate classically since the matrices we use are still small enough. The bosonic Hamiltonian for a 6-qubit representation took over 1000 Pauli factors to represent fully. To improve our performance, we can look for a more efficient representation of this model's Hamiltonian. We will develop a better variational trial wave function and initial state, as well as include the effect of fermions and heavy quarks on the potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

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

A fast band–Krylov eigensolver for macromolecular functional motion simulation on multicore architectures and graphics processors

We introduce a new iterative Krylov subspace-based eigensolver for the simulation of macromolecular motions on desktop multithreaded platforms equipped with multicore processors and, possibly, a graphics accelerator (GPU). The method consists of two stages, with the original problem first reduced into a simpler band-structured form by means of a high-performance compute-intensive procedure. This is followed by a memory-intensive but low-cost Krylov iteration, which is off-loaded to be computed on the GPU by means of an efficient data-parallel kernel. The experimental results reveal the performance of the new eigensolver. Concretely, when applied to the simulation of macromolecules with a few thousands degrees of freedom and the number of eigenpairs to be computed is small to moderate, the new solver outperforms other methods implemented as part of high-performance numerical linear algebra packages for multithreaded architectures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗