Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum programming”

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 361 records · Page 20

Describing adsorption of benzene, thiophene, and xenon on coinage metals by using the Zaremba–Kohn theory-based model

Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, and rVV10) or empirical methods (DFT+D, DFT+vdW, and DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals (Cu, Ag, and Au) is a typical example of systems where vdW interaction is significant. However, it is difficult to find a general method that reasonably describes both adsorption energy and geometry of even the simple prototypes of cyclic and heterocyclic aromatic molecules such as benzene (C 6 H 6 ) and thiophene (C 4 H 4 S), respectively, with reasonable accuracy. Here, we present an alternative scheme based on Zaremba–Kohn theory, called DFT+vdW-dZK. We show that unlike other popular methods, DFT+vdW-dZK and particularly SCAN+vdW-dZK give an accurate description of the physisorption of a rare-gas atom (xenon) and two small albeit diverse prototype organic molecules on the (111) surfaces of the coinage metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

2021 Stewardship Science Academic Programs Annual

I am a staff scientist at Sandia National Laboratories (SNL). I work on multiple fundamental-science projects and lead modeling/data-analysis for the stellar opacity experiments. I also work with theorists to refine plasma material-property calculations to make NNSA simulations more predictable. These are challenging and important problems for national security. In fact, our experiments raised questions about “opacity”, the property of matter that controls energy transport inside stars. This work was published by Nature due to its serious implication over broad applications, and our team received NNSA Defense Program Award of Excellence, both in 2015. Work at SNL is collegial and full of learning through interactions with the world’s finest experimentalists and theorists. When I started my Ph.D. program, working at a national lab was a dream position. Today, I collaborate with top scientists to solve mission-critical problems. The Stewardship Science Academic Programs (SSAP) and National Laser Users’ Facility (NLUF) played a vital role in the efficient learning and seamless transition to SNL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Low energy neutrino detection with a compact water-based liquid scintillator detector

In this study, the conceptual design and physics simulations of a near-field Water-based Liquid Scintillator (WbLS) detector placed 100 m from the Akkuyu Nuclear Power Plant (ANPP), currently under construction and aiming at being Turkey’s first nuclear power plant, is presented. The ANPP is an excellent opportunity for neutrino studies and the development of an R &D program for neutrino detectors in Turkey. The Reactor Neutrino Experiments of Turkey (RNET) program includes a compact detector with a 2.5-ton volume of WbLS and a ~ 30% photo-coverage, and the program is planned to be expanded with a medium-size 30-ton detector that will be an international testbed for WbLS and new detector technologies through low energy neutrino studies. In the following, the focus will be on the smaller ~2.5 ton detector, instrumented with 8-in. high quantum efficiency PMTs and two layers of cosmic veto paddles, covering all sides of the detector, to track and veto cosmic particles. Inverse Beta Decay (IBD) events from electron anti-neutrinos generated in the reactor core are simulated using the RAT-PAC simulation package and several liquids with different percentages of Liquid Scintillator (LS) and Gadolinium (Gd) are investigated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Three individually addressable spin qubits in a single molecule

Here, an asymmetric bis-phenol-β-diketone (H 4 L) has been designed as a ligand programmed to promote the assembly of a molecular arrangement composed of three magnetically exchanged [NiCu] pairs, each exhibiting an S = 1/2 spin. The latter are shown by EPR and magnetometry to be good qubit realizations and non-equivalent within the molecule in the solid state, as required for conditional quantum gates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The ITensor Software Library for Tensor Network Calculations

ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagrams, allowing users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of algorithms. After discussing the philosophy behind the ITensor approach, we show examples of each part of the interface including Index objects, the ITensor product operator, tensor factorizations, tensor storage types, algorithms for matrix product state (MPS) and matrix product operator (MPO) tensor networks, quantum number conserving block sparse tensors, and the NDTensors library. We also review publications that have used ITensor for quantum many-body physics and for other areas where tensor networks are increasingly applied. To conclude we discuss promising features and optimizations to be added in the future.

Fishman, Matthew↗

Codebase release 0.3 for ITensor

ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagrams, allowing users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of algorithms. After discussing the philosophy behind the ITensor approach, we show examples of each part of the interface including Index objects, the ITensor product operator, tensor factorizations, tensor storage types, algorithms for matrix product state (MPS) and matrix product operator (MPO) tensor networks, quantum number conserving block sparse tensors, and the NDTensors library. We also review publications that have used ITensor for quantum many-body physics and for other areas where tensor networks are increasingly applied. To conclude we discuss promising features and optimizations to be added in the future.

Fishman, Matthew↗

qSIEVE: Efficient qLDPC Memory via Systolic Movement in Atom Arrays

As quantum machines have scaled up in their number of qubits, significant research has turned towards increasing their fidelity with quantum error correction codes. Although promising results have been shown with the surface code, which only requires near-neighbor connections between qubits, the high qubit overhead of such local codes promises to be problematic. Consequently, recent work has explored non-local quantum LDPC (qLDPC) codes, which have good asymptotic encoding rates. Despite theoretical progress, hardware implementations of these codes have been a longstanding challenge. At the experimental level, demonstrations of movement based communication on atom arrays suggest this is a powerful new primitive to achieve non-local connectivity. Leveraging this, we present a protocol for implementing non-local qLDPC codes in hardware. Our protocol, qSIEVE, is a co-design of such codes with movement in atom arrays. qSIEVE defines a restricted family of qLDPC codes that can be implemented efficiently with systolic movement. We then quantify the utility of qSIEVE in the context of a complete fault tolerant architecture. We compare the cost of implementing benchmark programs in a standard, surface code only architecture and a mixed architecture where data is stored in qLDPC memory with qSIEVE and loaded to surface codes for computation.

Quantum error correction↗

Axion constraints from quiescent soft gamma-ray emission from magnetars

We report axion-like-particles (ALPs) emitted from the core of a magnetar can convert to photons in its magnetosphere. The resulting photon flux is sensitive to the product of (i) the ALP-nucleon coupling G an which controls the production cross section in the core and (ii) the ALP-photon coupling g aγγ which controls the conversion in the magnetosphere. We study such emissions in the soft-gamma-ray range (300 keV to 1 MeV), where the ALP spectrum peaks and astrophysical backgrounds from resonant Compton up scattering are expected to be suppressed. Using published quiescent soft-gamma-ray flux upper limits in five magnetars obtained with CGRO COMPTEL and INTEGRAL SPI/IBIS/ISGRI, we put limits on the product of the ALP-nucleon and ALP-photon couplings. We also provide a detailed study of the dependence of our results on the magnetar core temperature. We further show projections of our result for future Fermi-GBM observations. Our results motivate a program of studying quiescent soft-gamma-ray emission from magnetars with the Fermi-GBM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two photon exchange with nuclei from $e^+/e^-$ elastic cross section ratios

Hadronic box diagrams pose a significant source of radiative corrections to lepton scattering and β -decay measurements, yet their calculation remains highly model-dependent. Two-photon exchange is a box diagram generally relevant to most lepton scattering experiments, and is directly detectable via multiple observables. Single-spin asymmetries (SSA) are sensitive to the imaginary part of the TPE amplitude, while the ratio of $e^+/e^-$ elastic scattering cross sections is sensitive to the real part of the TPE amplitude. While SSA measurements have been reported for several complex nuclei, measurements of the $e^+/e^-$ ratio exist only for the proton. Proposed here is a measurement of the $e^+/e^-$ ratio on 4 He, a feasible nuclear target for which the ratio’s deviation from unity due to TPE is expected to be on order of 1–2%. Choosing a low Z nucleus minimizes Coulomb distortions, which could greatly overwhelm TPE for heavier nuclei. The proposed 19-day program, including measurements at three kinematic settings between 0.3 < Q 2 < 0.4 GeV 2 , would be made possible by the addition of a positron source at Jefferson Lab.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

INSPIRED: Inelastic neutron scattering prediction for instantaneous results and experimental design

Inelastic neutron scattering (INS) has unique advantages in probing how atoms vibrate and how the vibrations propagate and interact. Such dynamic information is crucial in understanding various material properties, from heat capacity, thermal conductivity, phase transitions, and chemical reactions to more exotic quantum behavior. The analysis and interpretation of the INS spectra often start from a model structure of the sample, followed by a series of calculations to obtain the simulated spectra to compare with experiments. The conventional way to perform such calculations usually requires significant time, computing resources, and specialized expertise. Here, we present a new program named INSPIRED (Inelastic Neutron Scattering Prediction for Instantaneous Results and Experimental Design), which enables users to perform rapid INS simulations in several different ways on their personal computers in just a few clicks, with the crystal structure as the only input file. Specifically, the users can choose a pre-trained symmetry-aware neural network (coupled with an autoencoder) to predict the phonon density of states (DOS), 1D S(E) and 2D S(|Q|,E) spectra for any given structure. One can also choose an existing density functional theory (DFT) calculation from a database (containing over 12,000 crystals), and quickly obtain the simulated INS spectra for single crystals and powders. It is also possible to use pre-trained universal machine learning force fields to relax a given crystal structure, calculate the phonon dispersion and DOS, and, subsequently, the INS spectra. All these functions are implemented with a PyQt graphic user interface. Finally, we expect these new tools will benefit broad user communities and significantly improve the efficiency of experiment design, execution, and data analysis for INS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Building Neutron Scattering Infrastructure in Louisiana for Advanced Materials (Final Report)

The main objective of the Louisiana Consortium for Neutron Scattering (LaCNS) is to build a major neutron scattering infrastructure capable of treating both soft and hard materials. The goal is to create a sustainable effort having the quality, breadth, and depth necessary to produce competitive proposals including collaborative and center type proposals. Our scientific aim is to understand the role of coupling in emergent complex materials and its impact on the structure/property relationship and to explore how to tune the key couplings to guide the design of materials with the desired properties. This naturally includes building a base of users of the Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR) in Louisiana; to train highly talented graduate and post-doctoral students in synthesis, neutron scattering, and simulation and modeling techniques, thereby helping to produce the next generation of scientist who use neutron scattering techniques as a crucial part of their research. The goal of our hard matter program is to understand the interplay between spin, charge, orbital, and lattice degrees of freedom in carefully selected complex materials. We have made considerable progress on number of complex systems including the oxide Sr3(Ru1-xMnx)2O7 where short-range magnetic ordering with anisotropic spin texture is initiated at the metal-insulator transition that clearly indicates a strong spin-charge coupling. Another critical area is in quantum materials. For example, in the topological semimetal BaMnSb2, we found a 3D canted antiferromagnetic Weyl semimetal with a 2D electronic structure and a nontrivial Berry phase. In addition, in the hexagonal chiral Mn1/3NbS2 system, we found soliton/soliton and soliton/antisoliton domain walls where the application of small fields or small currents can be used to control nanoscopic magnetic domains where the control of domain walls is crucially important for information storage. In addition, chemical transformation investigations were performed on the VISION instrument at Spallation Neutron Source (SNS). The soft matter effort was focused on understanding the role of non-covalent interactions on the structure and dynamics of fluid-based soft matter. One key focus was on sequence-defined (SD) amphiphilic peptoid polymers that allow encoding of molecular interactions and thereby systematic investigations of how charge directs the solution self-assembly of amphiphilic polymers in water. Another important area was on the dynamics of lipids self-assembled into membranes for exploring the permeability and mechanical using both neutron spin echo (NSE) spectroscopy and quasielastic neutron scattering (QENS) to distinguished between viscoelasticity and permeability at the molecular scale. Overall, the LaCNS project was quite successful, generating 145 publications and 245 presentations. Our graduate and post-doctoral students were well trained resulting in positions in national laboratories (Oak Ridge National Laboratory (ORNL), Argonne National Laboratory and Los Alamos National Laboratory), major research universities, and industry. We also developed a uniaxial pressure cell along with ORNL for SNS. A critical goal of this project was to establish a foundation for competing nationally in federally funded research programs. To this end, we were quite successful in generating over twenty-two federally and non-federally funded grants including awards from NSF, DOE, and DOD, and two early career awards. Equally important, we were able to secure a key major piece of instrumentation via a large ARO grant for a state-of-the-art electron microscope.

36 MATERIALS SCIENCE↗

Towards precise and accurate calculations of neutrinoless double-beta decay

We present the results of a National Science Foundation Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a roadmap for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The roadmap includes specific projects in theoretical and computational physics as well as the use of Bayesian methods to quantify both intra- and inter-model uncertainties. Here, the goal of this ambitious program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Programmable topotaxis of magnetic rollers in time-varying fields

We discuss how spatially uniform, time-periodic magnetic fields can be designed to power and direct the migration of ferromagnetic spheres up (or down) local gradients in the topography of a solid substrate. Our results are based on a dynamical model that considers the time-varying magnetic torques on the particle and its motion through the fluid at low Reynolds number. We use both analytical theory and numerical simulation to design magnetic fields that maximize the migration velocity up (or down) an inclined plane. We show how “topotaxis” of spherical particles relies on differences in the hydrodynamic resistance to rotation about axes parallel and perpendicular to the plane. Importantly, the designed fields can drive multiple independent particles to move simultaneously in different directions as determined by gradients in their respective environments. Experiments on ferromagnetic spheres provide evidence for topotactic motions up inclined substrates. The ability to program the autonomous navigation of driven particles within anisotropic environments is relevant to the design of colloidal robots.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superconductor to exciton condensate transition in a model copper-oxide material

Abstract Superconductivity and exciton condensation are fundamental phenomena in condensed matter physics, associated with the condensation of electron–electron and electron–hole pairs, respectively, into coherent quantum states. In this study, we present evidence of a superconductor to exciton condensate transition within the context of the three-band Hubbard model of copper-oxide-like materials. As the electron–electron repulsion increases, the superconducting phase is superseded by exciton condensation. In support of theoretical predictions—not yet realized experimentally—we observe the coexistence of the two condensates in the vicinity of the transition where the quantum states become a superposition of electron–electron and electron–hole condensates. Coexistence is rigorously computed from large eigenvalues and their eigenvectors in both the two-electron reduced density matrix (2-RDM) and the particle-hole RDM, which we obtain from a direct variational ground-state energy minimization with respect to the 2-RDM by semidefinite programming. We further discern that adjacent d orbitals and intervening p orbitals facilitate electron–electron pairing between copper orbitals, thereby supporting the superexchange mechanism for superconductivity. These observations suggest the feasibility of witnessing a superconductor to exciton condensate transition in copper-oxide analogs, bearing significant implications for identifying materials conducive to efficient transport processes.

Schouten, Anna O. (ORCID:0000000264255219)↗

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗