SEARCH · Engineering Papers

Results for “distance scale”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Multi-Functional Distributed Fiber Sensors for Pipeline Monitoring and Methane Detections. Final Report

As an abundant and cheap fossil energy source, natural gas has become a significant energy supply to support the United States’ economy. However, the large-scale extraction and utilization of natural gas also impose significant challenges on methane leakage. This problem is exacerbated by aging gas utility delivery systems, including interstate high-pressure pipelines, storage, and transmission facilities. This project aims to develop a cost-effective fiber optical sensing method that can perform multi-parameter real-time measurements of natural gas pipelines across long interrogation distances up to 100 km with 1-meter spatial resolution. This sensing tool can evaluate overall pipeline efficiency and reduce methane emissions for mid-stream methane infrastructures. To accomplish this objective, research and development efforts funded by this project have resulted in the following accomplishments: This project successfully has developed new functional sensory polymer materials using Metal-Organic Frameworks (MOFs) that can be coated on optical fiber through the reel-to-reel coating process. Functional polymer-coated optical fibers can perform sensitive methane detection through evanescence wave interaction and strain-based measurements to achieve 1% detection sensitivities. The new sensors fibers support both distributed measurements and multiplexed fiber sensors array for multi-point measurements. This project developed and optimized a new multi-core optical fiber that supports simultaneous and distributed measurements of strain and temperatures with 1-meter spatial resolutions across up to 100-km interrogation distance. This new fiber, combined with sensory polymercoated fiber, could perform both distributed temperature and methane detections. This project developed a new artificial intelligence big data algorithm approach that can effectively analyze high-resolution data harnessed by distributed fiber sensors to protect natural gas pipelines against external threats and detect internal defects induced by corrosion. Working with our industry partner, this project developed new optical fibers that support fiber sensor fabrications through polymer coating after the fibers are drawn. These new fibers eliminate the need for direct sensor fabrication when the fiber is fabricated on a fiber draw tower, which drastically expands fiber sensors' applicability. This research project has significantly advanced the distributed fiber sensing technology. It will dramatically increase the applicability and adaptability of distributed fiber sensors for a wide array of applications in energy, sustainability, and environmental science, including structural health monitoring of natural gas pipelines, oil infrastructures, hydrogen facilities, and environmental monitoring of carbon storage sites, water supply systems, and others.

03 NATURAL GAS↗

DOE OSTI2020

The Preservation of Super-Earths and the Emergence of Gas Giants after Their Progenitor Cores Have Entered the Pebble-isolation Phase

The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distributions indicate that some cores have transformed into gas giants through runaway gas accretion at ≳1 au from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically ~ 10M ⊕ ) to enter a pebble-isolation phase. The accumulation of approximately millimeter- to meter-sized pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced submillimeter grains pass through the barrier, elevate the effective dust opacity, and reduce the radiative flux in the core's envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (~0.1 au), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At an intermediate distance (~1 au) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few astronomical unit, the transformation of more-massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗

DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Relaxation Dynamics and Ion Conduction of Poly(Ethylene Carbonate/Ethylene Oxide) Copolymer-Based Electrolytes

Poly(ethylene carbonate/ethylene oxide) P(EC/EO) copolymer was synthesized from EC monomer as a new matrix component of solid polymer electrolytes (SPEs), and the correlation between polymer dynamics and ionic conductivity in the P(EC/EO)-based SPEs at various lithium salt weight fractions (w Li ) was investigated by dielectric spectroscopy, rheology, and molecular dynamics (MD) simulations. With the addition of lithium salt, the molecular motions at various scales (i.e., from local to segmental and global motion scales) were found to change as follows: (i) the local motion of P(EC/EO) was slightly accelerated, and at w Li ≥ 0.22, it was separated into two relaxation modes, one faster and the other slower than that of P(EC/EO), (ii) the segmental motion of P(EC/EO) became significantly slower with the addition of lithium salt, and (iii) the global terminal relaxation became slightly slower with the increase of w Li . Thus, the motion of P(EC/EO) at various scales in SPEs changed with different trends with increasing lithium salt. It was also found that the temperature dependence of the ionic conductivity was described by a Vogel–Fulcher–Tammann-type function, which was correlated with the segmental motion, as is known in other SPE systems. Furthermore, from MD simulations, lithium ions are coordinated with the EC units more likely and with a closer distance than EO units in P(EC/EO), and the coordination number of oxygens around one lithium ion was estimated as 6–7 at the distance of 0.30 nm. Here, these results suggest that the coordination structure of the oxygens around the lithium ion is slightly disrupted due to the presence of EC units compared to the chelating structure of pure PEO electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Cone beam neutron interferometry: From modeling to applications

Phase-grating moiré interferometers (PGMIs) have emerged as promising candidates for the next generation of neutron interferometry, enabling the use of a polychromatic beam and manifesting interference patterns that can be directly imaged by existing neutron cameras. However, the modeling of the various PGMI configurations is limited to cumbersome numerical calculations and backward propagation models which often do not enable one to explore the setup parameters. Here we generalize the Fresnel scaling theorem to introduce a k -space model for PGMI setups illuminated by a cone beam, thus enabling an intuitive forward propagation model for a wide range of parameters and experimental setups. The interference manifested by a PGMI is shown to be a special case of the Talbot effect, and the optimal fringe visibility is shown to occur at the moiré location of the Talbot distances. We derive analytical expressions for the contrast and the propagating intensity profiles in various conditions and provide the first analysis of the PGMI dark-field imaging signal when considering sample characterization. The model's predictions are compared to experimental measurements and good agreement is found between them. Last, we propose and experimentally verify a method to recover contrast at typically inaccessible PGMI autocorrelation lengths. The presented work provides a toolbox for analyzing and understanding existing PGMI setups and their future applications, for example extensions to two-dimensional PGMIs and characterization of samples with nontrivial structures. Published by the American Physical Society 2024

Sarenac, D. (ORCID:0000000185753367)↗

DOE OSTI2021

An Integrated View of Hydraulic Induced Fracture Geometry in Hydraulic Fracture Test Site 2

Until relatively recently, most of our understanding about the geometry of hydraulically induced fractures in unconventional reservoirs has been inferred by patterns observed via microseismic monitoring. Important information, such as fracture-azimuth, fracture-length, height-growth and other geometric data has been determined from microseismic “event clouds”, which in turn has been used to interpret a stimulated reservoir volume (SRV). It is important to note that relatively limited microseismic results have been acquired in conjunction with other diagnostics and validated in a consistent manner. The Hydraulic Fracture Test Site 2 (HFTS2) in the Permian Delaware Basin provides a unique opportunity to compare the frac geometry interpretations derived from multiple frac diagnostic tools and at multiple scales. The integrated view that emerges from all the frac diagnostic tools is that in HFTS2 the geometries of hydraulically induced fractures are not random and with a highly complex branching architecture. In fact, the fracs are mostly vertical, parallel plana domains, and both vertically and horizontally asymmetrical. The hydraulically induced fractures have a very consistent azimuth with a strike that matches SHmax. At the distance corresponding to the spacing between nearby wells in HFTS2, ~ 660 ft, the pattern of strain interceptions observed from Low Frequency Distributed Acoustic Sensing (LF-DAS) shows far-field dimensions that are consistent with the stage length dimensions at the stimulated well. In general, larger stages with higher number of clusters create larger number of interceptions and wider stimulated intervals than that of shorter stages with the same cluster to cluster spacing and fewer clusters. The hydraulically induced fractures also show a clear tendency to preferentially grow upwards as shown by the microseismic data, by the relative lower number of LF-DAS interceptions in nearby deeper wells and finally by the observed LF-DAS strain patterns in the vertical observation well. In the near field, the geometries of the fracture domains are also consistent. As shown by the distributed strain measurements obtained during production via the monitoring of Rayleigh Frequency Shift (RFS), each cluster has a separate, non-overlapping, frac-zone-domain that are generally centered around the locations of each perforation cluster. Overall, the hydraulically induced fractures in HFTS2 are interpreted as occurring in swarms associated with individual clusters, each with its own unique geometry but with similar predictable geometric and propagation tendencies. The comprehensive frac diagnostic program in HFTS2, the well-pad layout that permits the investigation of the Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R003/2477511/urtec-2021-5396-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5396 2 impact that nearby producing wells has on frac geometry, together with an innovative Design of Experiment (DofE) that includes sequential-fracing, consistent cluster-to-cluster distances and few single-perforation cluster stages, has allowed us to gain a unique perspective about the geometry of the induced hydraulic fractures on this pad. Finally, the HFTS2 dataset has also provided us with a better understanding of the applicability and limitations of the different diagnostic tools.

58 GEOSCIENCES↗

DOE OSTI2021

Guidance for Monitoring Passive Groundwater Remedies Over Extended Time Scales

Passive remediation can be appropriate where natural processes and actions such as institutional controls mitigate exposure to contaminated groundwater, achieving remedial action objectives and protectiveness of human health and the environment. Monitored natural attenuation (MNA) is a prevalent passive remediation strategy supported by a regulatory framework and monitoring design guidance. However, long-term passive remedies are usually selected in combination with at least one active remedy, such as source removal, in situ treatment, or pump-and-treat, functioning as a complementary method for achieving remediation objectives and meeting the applicable statutory and regulatory requirements. However, MNA and existing monitoring guidance primarily target situations where the remedial action objectives are met within a few decades. When time scales for passive remediation extend to many decades (50 years or more), a corresponding change in monitoring strategy is needed to adapt to the extended time scale. This document provides guidance for implementing an extended-scale monitoring (ESM) approach appropriate for long-duration passive remediation. Extended-scale is defined in this document with respect to time (i.e., a long duration of remediation) and a large enough physical scale such that downstream receptors will not be impacted within the remediation timeframe. ESM applies to slow-moving groundwater contaminant plumes and emphasizes monitoring primarily for potential exposure pathways. For this approach, the primary monitoring objective is to demonstrate that the plume diminishes before reaching the receptor zone or point of compliance and/or a receptor does not receive concentrations above the compliance limit. While the overall objectives of protecting human health and the environment are the same as for plumes where remediation can occur over a shorter time period, the time scales between decisions are longer and the dynamics of plume evolution are slower. To this end, a scenario-based strategy is described for different plume and source conditions, defining a containment and receptor zones. The containment zone is the area where the risk of exposure to groundwater contamination can be mitigated (e.g., through institutional controls) during the remediation time period. The receptor zone is defined as the area where exposure to groundwater contamination cannot be mitigated and compliance concentration standards must be met. Within the containment zone, slow plume migration may occur, leading to concentrations that exceed compliance standards. However, where distance to the receptor zone is large relative to plume migration and attenuation rate, this approach can be protective of the receptor zone. Selection of a long-duration passive remedy needs to be based on sufficient understanding of contaminant sources, hydrogeology, and contaminant plumes. A strong technical basis, supported by predictive analysis, is recommended to substantiate that contamination is expected to stay within the containment zone during the active remediation and attainment phase of the remedy, and diminish to meet compliance standards within the extended timeframe prior to reaching the receptor zone (e.g., many decades or even centuries). The ESM approach is based on verification of plume behavior and not on detailed plume dynamics. Monitoring is conducted to confirm expected behavior with an emphasis on exposure pathways to verify that plumes remain contained in areas where the protectiveness objectives can be met. ESM should not be adopted if there is significant risk of the plume extending beyond the containment zone. Given the slow movement within the containment zone, less frequent sampling is required relative to approaches used for conventional-scale remediation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗

DOE OSTI2025

DESI DR2 results. I. Baryon acoustic oscillations from the Lyman alpha forest

We present the baryon acoustic oscillation (BAO) measurements with the Lyman-𝛼 (Ly⁢𝛼) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the autocorrelation of the Ly⁢𝛼 forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of 2 larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped Ly⁢𝛼 absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at 𝑧 eff =2.33. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in Ly⁢𝛼 BAO analysis for the first time. We measure the ratios 𝐷 𝐻 ⁡(𝑧 eff )/𝑟 𝑑 = 8.632 ± 0.098 ± 0.026 and 𝐷 𝑀 ⁡(𝑧 eff )/𝑟 𝑑 = 38.99 ± 0.52 ± 0.12, where 𝐷 𝐻 = 𝑐/𝐻⁡(𝑧) is the Hubble distance, 𝐷 𝑀 is the transverse comoving distance, 𝑟 𝑑 is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.

baryon acoustic oscillations↗

DOE OSTI2025

Nonlocal nucleon matrix elements in the rest frame

Extracting parton structure from lattice quantum chromodynamics (QCD) calculations requires studying the coordinate scale 𝑧 3 dependence of the matrix elements of bilocal operators. The most significant contribution comes from the 𝑧 3 dependence induced by ultraviolet (UV) renormalization of the Wilson line. We demonstrate that the next-to-leading order perturbative calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice QCD rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. The residual discrepancies can be modeled by a leading effect from the structure of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Design and fabrication of ion traps for low RF power dissipation

Large surface-electrode ion traps with multiple trapping regions and junctions are a natural approach to scaling trapped ion quantum computers, supporting the connectivity and ion counts necessary for complex quantum algorithms. However, a major hurdle in this scaling is on-chip power dissipation from the applied RF voltage, which increases at a rate between linear and cubic relative to trap size, depending on whether the losses are dielectric or Ohmic. Here, we present two versions of a trap with features designed to reduce both types of RF power dissipation. The first variant contains a raised RF electrode that increases the electrode–ground distance to reduce capacitance. Different DC voltage sources are demonstrated on this trap to show that technical noise before the filter remains the dominant source of voltage noise and therefore motional heating. The second variant additionally includes a method for removing dielectric from beneath the RF electrode to further reduce dielectric losses. These traps were demonstrated at room temperature with 40 Ca + ions. In conclusion, the similar heating rates and heating rate axial frequency dependencies between 2.4 and 3.0 MHz illustrate that this dielectric modification is not detrimental to trap performance.

Sterk, J. D. [Sandia National Laboratories (SNL-NM↗

DOE OSTI2025

Preliminary Investigation of Supernova Neutrino Detection Backgrounds in MicroBooNE

MicroBooNE provides a testbed for studying MeV-scale activity relevant to supernova neutrino detection in LArTPCs. This work focuses on classifying isolated mm-scale energy depositions ( blips ), in contrast to longer track depositions with lengths of ~O(10cm), as possible νₑ CC or elastic scatter candidates from a supernova burst. Cosmogenic blip backgrounds are characterized using CRT-tagged muons and spatial correlations with decay products. Selection cuts are defined for rejecting Tl-208 decay radiogenic blip backgrounds and cosmogenic blip backgrounds; they are based on blip energy and distance from nearby cosmic muon tracks. Results affirm strategies for background rejection and signal efficiency in future low-energy searches, including DUNE.

Binau, Amelia [Indiana U.]↗

DOE OSTI2024

Intra- and inter-subtype HIV diversity between 1994 and 2018 in southern Uganda: a longitudinal population-based study

There is limited data on human immunodeficiency virus (HIV) evolutionary trends in African populations. We evaluated changes in HIV viral diversity and genetic divergence in southern Uganda over a 24-year period spanning the introduction and scale-up of HIV prevention and treatment programs using HIV sequence and survey data from the Rakai Community Cohort Study, an open longitudinal population-based HIV surveillance cohort. Gag (p24) and env (gp41) HIV data were generated from people living with HIV (PLHIV) in 31 inland semi-urban trading and agrarian communities (1994–2018) and four hyperendemic Lake Victoria fishing communities (2011–2018) under continuous surveillance. HIV subtype was assigned using the Recombination Identification Program with phylogenetic confirmation. Inter-subtype diversity was evaluated using the Shannon diversity index, and intra-subtype diversity with the nucleotide diversity and pairwise TN93 genetic distance. Genetic divergence was measured using root-to-tip distance and pairwise TN93 genetic distance analyses. Demographic history of HIV was inferred using a coalescent-based Bayesian Skygrid model. Evolutionary dynamics were assessed among demographic and behavioral population subgroups, including by migration status. 9931 HIV sequences were available from 4999 PLHIV, including 3060 and 1939 persons residing in inland and fishing communities, respectively. In inland communities, subtype A1 viruses proportionately increased from 14.3% in 1995 to 25.9% in 2017 (P < .001), while those of subtype D declined from 73.2% in 1995 to 28.2% in 2017 (P < .001). The proportion of viruses classified as recombinants significantly increased by nearly four-fold from 12.2% in 1995 to 44.8% in 2017. Inter-subtype HIV diversity has generally increased. While intra-subtype p24 genetic diversity and divergence leveled off after 2014, intra-subtype gp41 diversity, effective population size, and divergence increased through 2017. Intra- and inter-subtype viral diversity increased across all demographic and behavioral population subgroups, including among individuals with no recent migration history or extra-community sexual partners. This study provides insights into population-level HIV evolutionary dynamics following the scale-up of HIV prevention and treatment programs. Continued molecular surveillance may provide a better understanding of the dynamics driving population HIV evolution and yield important insights for epidemic control and vaccine development.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Gleipnir: toward practical error analysis for Quantum programs

Practical error analysis is essential for the design, optimization, and evaluation of Noisy Intermediate-Scale Quantum(NISQ) computing. However, bounding errors in quantum programs is a grand challenge, because the effects of quantum errors depend on exponentially large quantum states. In this work, we present Gleipnir, a novel methodology toward practically computing verified error bounds in quantum programs. Gleipnir introduces the (ρ,δ)-diamond norm, an error metric constrained by a quantum predicate consisting of the approximate state ρ and its distance δ to the ideal state ρ. This predicate (ρ,δ) can be computed adaptively using tensor networks based on the Matrix Product States. Gleipnir features a lightweight logic for reasoning about error bounds in noisy quantum programs, based on the (ρ,δ)-diamond norm metric. Furthermore, our experimental results show that Gleipnir is able to efficiently generate tight error bounds for real-world quantum programs with 10 to 100 qubits, and can be used to evaluate the error mitigation performance of quantum compiler transformations.

Tao, Runzhou↗

DOE OSTI2021

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

UV-Induced Reaction Pathways in Bromoform Probed with Ultrafast Electron Diffraction

For many chemical reactions, it remains notoriously difficult to predict and experimentally determine the rates and branching ratios between different reaction channels. This is particularly the case for reactions involving short-lived intermediates, whose observation requires ultrafast methods. The UV photochemistry of bromoform (CHBr 3 ) is among the most intensely studied photoreactions. Yet, a detailed understanding of the chemical pathways leading to the production of atomic Br and molecular Br 2 fragments has proven challenging. In particular, the role of isomerization and/or roaming and their competition with direct C–Br bond scission has been a matter of continued debate. Here, in this work, gas-phase ultrafast megaelectronvolt electron diffraction (MeV-UED) is used to directly study structural dynamics in bromoform after single 267 nm photon excitation with femtosecond temporal resolution. The results show unambiguously that isomerization contributes significantly to the early stages of the UV photochemistry of bromoform. In addition to direct C–Br bond breaking within <200 fs, formation of iso-CHBr 3 (Br-CH-Br-Br) is observed on the same time scale and with an isomer lifetime of >1.1 ps. The branching ratio between direct dissociation and isomerization is determined to be 0.4 ± 0.2:0.6 ± 0.2, i.e., approximately 60% of molecules undergo isomerization within the first few hundred femtoseconds after UV excitation. The structure and time of formation of iso-CHBr 3 compare favorably with the results of an ab initio molecular dynamics simulation. The lifetime and interatomic distances of the isomer are consistent with the involvement of a roaming reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

The effects of pre-plasma scale length on the relativistic electron beam directionality

The effects of pre-plasmas on the electron beam directionality was experimentally and numerically investigated. Single material and layered targets made of Ti and/or CH were used to simultaneously measure high-energy (≥3 MeV) electrons along two directions, pre-pulse energy and pre-plasma density. The electron directionality is quantified by using a new parameter, the electron energy ratio of the total kinetic energies along the two directions. Measurements and radiation–hydrodynamic (RH) simulations show that a large (≥3.5 μm) plasma scale length at the critical surface enhances electrons along the laser axis, and such pre-plasma conditions could only be achieved with the CH targets. Particle-in-cell simulations were performed on the RH generated pre-plasmas from Ti and CH targets, and the results show that the CH target provided conditions for higher forward momentum gains by electrons. First, the CH target allowed longer distances for electrons to interact with laser. Second, the intense laser pulse modified the critical surface, but the resulting surface differed. The CH target resulted in a smooth surface where a retro-reflection was observed while the Ti target resulted in a rippled surface that scattered the reflected light. As results, the CH electrons gained higher forward momentum via a direct-laser-acceleration in the counter propagating laser fields. So the results presented in this article show a way of controlling the high-energy electron directionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗