SEARCH · Engineering Papers

Results for “log”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Ubiquity of density slope oscillations in the central regions of galaxy and cluster-sized systems

One usually thinks of a radial density profile as having a monotonically changing logarithmic slope, such as in NFW or Einasto profiles. However, in two different classes of commonly used systems, this is often not the case. These classes exhibit non-monotonic changes in their density profile slopes which we call oscillations for short. We analyze these two unrelated classes separately. Class 1 consists of systems that have density oscillations and that are defined through their distribution function f ( E ), or differential energy distribution N ( E ), such as isothermal spheres, King profiles, or DARKexp, a theoretically derived model for relaxed collisionless systems. Systems defined through f ( E ) or N ( E ) generally have density slope oscillations. Class 1 system oscillations can be found at small, intermediate, or large radii but we focus on a limited set of Class 1 systems that have oscillations in the central regions, usually at log( r / r {sub −2}) ∼< −2, where r {sub −2} is the largest radius where d log(ρ)/ d log( r ) = −2. We show that the shape of their N ( E ) can roughly predict the amplitude of oscillations. Class 2 systems which are a product of dynamical evolution, consist of observed and simulated galaxies and clusters, and pure dark matter halos. Oscillations in the density profile slope seem pervasive in the central regions of Class 2 systems. We argue that in these systems, slope oscillations are an indication that a system is not fully relaxed. We show that these oscillations can be reproduced by small modifications to N ( E ) of DARKexp. These affect a small fraction of systems' mass and are confined to log( r / r {sub −2}) ∼< 0. The size of these modifications serves as a potential diagnostic for quantifying how far a system is from being relaxed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Rotation curves and dynamical masses of MaNGA barred galaxies

In this paper we analyse a sample of 46 barred galaxies of MaNGA. Our goal is to investigate the stellar kinematics of these galaxies and obtain their rotation curves. Additionally, we aim to derive the total stellar and dynamical masses, as well as the maximum rotation velocity, in order to examine their distributions and scaling relations. Using the Pipe3D dataproducts publicly available we obtained the rotation curves, which were fitted considering two components of an axisymmetric Miyamoto–Nagai gravitational potential. We found a wide range of the maximum rotation velocities (117–340 ${\rm km\, s^{-1}}$), with a mean value of 200 ${\rm km\, s^{-1}}$. In addition we found that the total stellar and dynamical masses are in the range of log(M star /M ⊙ ) = 10.1−11.5, with a mean value of log(M star /M ⊙ ) = 10.8, and log(M dyn /M ⊙ ) = 10.4−12.0, with a mean value of log(M dyn /M ⊙ ) = 11.1, respectively. We found a strong correlation between dynamical mass and maximum velocity, between maximum velocity and magnitude, and between stellar mass and maximum velocity. According to these results, barred galaxies exhibit similar behaviour to that of normal spiral galaxies with respect to these relations, as well as in terms of the distribution of their dynamical mass and maximum rotation velocity. However, we found that the distribution of stellar masses of barred galaxies is statistically different from other samples including non-barred galaxies. Finally, analysing the galaxies that show nuclear activity, we find no difference with the rotation curves of normal galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Characterization of sub-micrometre-sized voids in fixed human brain tissue using scanning X-ray microdiffraction

Using a 5 µm-diameter X-ray beam, we collected scanning X-ray microdiffraction in both the small-angle (SAXS) and the wide-angle (WAXS) regimes from thin sections of fixed human brain tissue from Alzheimer's subjects. The intensity of scattering in the SAXS regime of these patterns exhibits essentially no correlation with the observed intensity in the WAXS regime, indicating that the structures responsible for these two portions of the diffraction patterns, which reflect different length scales, are distinct. SAXS scattering exhibits a power-law behavior in which the log of intensity decreases linearly with the log of the scattering angle. The slope of the log–log curve is roughly proportional to the intensity in the SAXS regime and, surprisingly, inversely proportional to the intensity in the WAXS regime. We interpret these observations as being due to the presence of sub-micrometre-sized voids formed during dehydration of the fixed tissue. The SAXS intensity is due largely to scattering from these voids, while the WAXS intensity derives from the secondary structures of macromolecular material surrounding the voids. The ability to detect and map the presence of voids within thin sections of fixed tissue has the potential to provide novel information on the degradation of human brain tissue in neurodegenerative diseases.

Chemistry↗

DOE OSTI2016

SU-F-T-462: Lessons Learned From a Machine Incident Reporting System

Purpose: Linear accelerators must operate with minimal downtime. Machine incident logs are a crucial tool to meet this requirement. They providing a history of service and demonstrate whether a fix is working. This study investigates the information content of a large department linear accelerator incident log. Methods: Our department uses an electronic reporting system to provide immediate information to both key department staff and the field service department. This study examines reports for five linac logs during 2015. The report attributes for analysis include frequency, level of documentation, who solved the problem, and type of fix used. Results: Of the reports, 36% were documented as resolved. In another 25% the resolution allowed treatment to proceed although the reported problem recurred within days. In 5% only intermediate troubleshooting was documented. The remainder lacked documentation. In 60% of the reports, radiation therapists resolved the problem, often by clearing the appropriate faults or reinitializing a software or hardware service. 22% were resolved by physics and 10% by field service engineers. The remaining 8% were resolved by IT, Facilities, or resolved spontaneously. Typical fixes, in order of scope, included clearing the fault and moving on, closing and re-opening the patient session or software, cycling power to a sub-unit, recalibrating a device (e.g., optical surface imaging), and calling in Field Service (usually resolving the problem through maintenance or component replacement). Conclusion: The reports with undocumented resolution represent a missed opportunity for learning. Frequency of who resolves a problem scales with the proximity of the person’s role (therapist, physicist, or service engineer), which is inversely related to the permanence of the resolution. Review of lessons learned from machine incident logs can form the basis for guidance to radiation therapists and medical physicists to minimize equipment downtime and ensure safe operation.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2018

Bidisperse and polydisperse suspension rheology at large solid fraction

At the same solid volume fraction, bidisperse and polydisperse suspensions display lower viscosities, and weaker normal stress response, compared to monodisperse suspensions. The reduction of viscosity associated with size distribution can be explained by an increase of the maximum flowable, or jamming, solid fraction. In this work, concentrated or "dense" suspensions are simulated under strong shearing, where thermal motion and repulsive forces are negligible, but we allow for particle contact with a mild frictional interaction with interparticle friction coefficient of 0.2. Aspects of bidisperse suspension rheology are first revisited to establish that the approach reproduces established trends; the study of bidisperse suspensions at size ratios of large to small particle radii (2 to 4) shows that a minimum in the viscosity occurs for zeta slightly above 0.5, where zeta=phi_{large}/phi is the fraction of the total solid volume occupied by the large particles. The simple shear flows of polydisperse suspensions with truncated normal and log normal size distributions, and bidisperse suspensions which are statistically equivalent with these polydisperse cases up to third moment of the size distribution, are simulated and the rheologies are extracted. Prior work shows that such distributions with equivalent low-order moments have similar phi_{m}, and the rheological behaviors of normal, log normal and bidisperse cases are shown to be in close agreement for a wide range of standard deviation in particle size, with standard correlations which are functionally dependent on phi/phi_{m} providing excellent agreement with the rheology found in simulation. The close agreement of both viscosity and normal stress response between bi- and polydisperse suspensions demonstrates the controlling in influence of the maximum packing fraction in noncolloidal suspensions. Microstructural investigations and the stress distribution according to particle size are also presented.

Pednekar, Sidhant↗

DOE OSTI2019

Volcanic stratigraphy and age model of the Kimama deep borehole (Project Hotspot): Evidence for 5.8 million years of continuous basalt volcanism, central Snake River Plain, Idaho

The Snake River Plain of central Idaho represents the world’s best example of a mantle hotspot track impinging upon continental crust and provides a record of bimodal volcanism extending over 12 Ma to the present. Project Hotspot recovered almost 2 km of continuous drill core from the Kimama borehole, located in central Idaho on the axial volcanic zone. The Kimama drill core represents the most complete record of mafic volcanism along the Yellowstone–Snake River Plain hotspot track. A total of 432 basalt flow units, representing 183 basalt flows, 78 basalt flow groups, and 34 super groups, along with 42 sediment interbeds are recognized using volcanic facies observations, stratigraphic relationships, borehole geophysical logs, and paleosecular variation in magnetostratigraphy. Rhyolite and other non-basaltic volcanic materials were not encountered in the drill core. Ages for six basalt lava flows were determined by 40 Ar/ 39 Ar using incremental heating experiments. Paleomagnetic inclination was measured on over 1200 samples collected at roughly 2-m-depth intervals, yielding mean values of paleosecular variation between ±50° to ±70° in Kimama flow groups, close to the expected 61° axial dipole average for the Kimama borehole location. Twenty-three magnetic reversals were identified and correlated to dated geomagnetic chrons and subchrons and compared with the 40 Ar/ 39 Ar radiometric ages. A linear fit to 40 Ar/ 39 Ar dates, geomagnetic chron and subchron boundaries, and volcanogenic zircon U-Pb ages defines a mean accumulation rate of ~320 m/m.y. and extrapolates to a bottom hole age of 6.3 Ma. Average thicknesses of lithologic units increase from 2.7 m (sediment), 4 m (flow units), 10 m (flows), 23 m (flow groups), to 53 m (super groups). On average, one lava flow inundated the Kimama borehole location every 33 k.y. Intercalated sediments, ranging from 0.06 to 24.5 m thick, make up roughly 6% of the drill core and indicate lulls in local volcanic activity that may have lasted up to 77 k.y. Neutron and gamma-ray logs supplement observations from the drill cores: neutron logs document individual flow units through the contrast between massive flow interiors and more porous flow surfaces, and gamma-ray logs document the depth and thickness of sedimentary interbeds and high–K-Fe basalts. The 5.8 m.y. duration of basaltic volcanism in the Kimama drill core implies a steady rate of volcanism, indicating a relatively stable rate of mantle upflow along the lithosphere-mantle boundary in the wake of Yellowstone–Snake River Plain plume volcanism.

58 GEOSCIENCES↗

DOE OSTI2022

Cretaceous sequence stratigraphy of the northern Baltimore Canyon Trough: Implications for basin evolution and carbon storage

We evaluate the Cretaceous stratigraphy and carbon sequestration potential of the northern Baltimore Canyon Trough (NBCT) using >10,000 km of multi-channel seismic profiles integrated with geophysical logs, biostratigraphy, and lithology from 29 offshore wells. We identify and map six sequences resolved primarily at the stage level. Accommodation was dominated by thermal and non-thermal subsidence, though sequence boundaries correlate with regional and global sea-level changes, and the record is modified by igneous intrusion, active faulting, and changes in sediment supply and sources. Our stratigraphic maps illustrate a primary southern (central Appalachian) Early Cretaceous source that migrated northward during the Aptian and Albian. During the Cenomanian, sedimentation rates in the NBCT increased and depocenters shifted northward and landward. We show that deposition occurred in three phases: (1) earlier Cretaceous paleoenvironments were primarily terrestrial indicated by variable amplitude, chaotic seismic facies, serrated gamma logs, and heterolithic sandstones and mudstones with terrestrial microfossils; (2) the Albian to Cenomanian was dominated by deltaic paleoenvironments indicated by blocky, funnel-shaped, gamma-ray logs and clinoforms characterized by continuous high-amplitude seismic reflections with well-defined terminations; and (3) the Cenomanian and younger was marine shelf, inferred from mudstoneprone lithologies, peak gamma-ray values in well logs, and foraminiferal evidence. Long-term transgression and maximum water depths at the Cenomanian/Turonian boundary correlative with Ocean Anoxic Event 2 were followed by a regression and relative sea-level fall. We show that porous and permeable sandstones of three Aptian to Cenomanian highstand systems tracts are high-volume reservoirs for supercritical CO 2 storage that are confined by overlying deep water mudstones.

58 GEOSCIENCES↗

DOE OSTI2016

Scale-Free Compact Routing Schemes in Networks of Low Doubling Dimension

In this work, we consider compact routing schemes in networks of low doubling dimension, where the doubling dimension is the least value α such that any ball in the network can be covered by at most 2 α balls of half radius. There are two variants of routing-scheme design: (i) labeled (name-dependent) routing, in which the designer is allowed to rename the nodes so that the names (labels) can contain additional routing information, for example, topological information; and (ii) name-independent routing, which works on top of the arbitrary original node names in the network, that is, the node names are independent of the routing scheme. In this article, given any constant ε ϵ (0, 1) and an n-node edge-weighted network of doubling dimension α ϵ O(loglog n), we present —a (1 + ε)-stretch labeled compact routing scheme with Γlog n-bit routing labels, O(log 2 n/loglog n)-bit packet headers, and ((1/ε) O(α) log 3 n)-bit routing information at each node; —a (9 + ε)-stretch name-independent compact routing scheme with O(log 2 n/loglog n)-bit packet headers, and ((1/ε) O(α) log 3 n)-bit routing information at each node. In addition, we prove a lower bound: any name-independent routing scheme with o(n (ε/60)2 ) bits of storage at each node has stretch no less than 9 - ε for any ε ϵ (0, 8). Therefore, our name-independent routing scheme achieves asymptotically optimal stretch with polylogarithmic storage at each node and packet headers. Note that both schemes are scale-free in the sense that their space requirements do not depend on the normalized diameter Δ of the network. Finally, we also present a simpler nonscale-free (9 + ε)-stretch name-independent compact routing scheme with improved space requirements if Δ is polynomial in n.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Translation-Invariant Quantum Algorithms for Ordered Search are Optimal

Ordered search is the task of finding an item in an ordered list using comparison queries. The best exact classical algorithm for this fundamental problem uses [log 2 n] queries for a list of length n. Quantum computers can achieve a constant-factor speedup, but the best possible coefficient of log 2 n for exact quantum algorithms is only known to lie between (ln2)/π ≈ 0.221 and 4/log 2 605 ≈ 0.4333. We consider a special class of translation-invariant algorithms with no workspace, introduced by Farhi, Goldstone, Gutmann, and Sipser, that has been used to find the best known upper bounds. First, we show that any bounded-error, k-query quantum algorithm for ordered search can be implemented by a k-query algorithm in this special class. Second, we use linear programming to show that the best exact 5-query quantum algorithm can search a list of length 7265, giving an ordered search algorithm that asymptotically uses 5 log 7265 n ≈ 0.390 log 2 n quantum queries.

Translation-invariant quantum algorithms↗

DOE OSTI2023

Impedance Response Influenced by Variability in the Random Distribution of Physical Properties of Coated Materials in Two-Dimensional Space

Heterogeneous physical characteristics of a system featuring a single-layer film on a metallic surface have been explored via its impedance response. The Nyquist plot showed a distorted semicircle, indicative of the system's unique distribution characteristics. Utilizing a copula-based probability method, a two-dimensional deterministic impedance model was successfully integrated, accounting for spatial physical properties such as permittivity and electrical conductivity. This strategy enabled in-depth exploration and mechanistic quantification of a broad spectrum of properties. A quantitative understanding of impedance signal alterations, characterized by normally or log-normally correlated variables, was achieved through the variation in aspect ratio and characteristic frequency of the impedance spectra. Log-normally distributed electrical properties provided a superior representation of the distorted impedance spectra. As coefficient of variation (CV) values fluctuated, the aspect ratio and characteristic frequency showed heightened sensitivity to log-normal permittivity compared to log-normal electrical conductivity. Notably, a marked positive linear correlation between electrical properties resulted in an impedance response that approximated perfect semicircular spectra. The variability in the electrical properties' distribution was demonstrated by considering the correlation coefficient between electrical conductivity and the z-direction position. Furthermore, the highest aspect ratio of the impedance spectra was observed when the electrical conductivity was randomly distributed across the z-direction space.

36 MATERIALS SCIENCE↗

DOE OSTI2019

numericalFunctions

The numericalFunction C library is designed to handle point-wise interpolative data representing a mathematical function. The library currently handles 1-d functions (i.e., functions of the format f(x)). The main object is the ptwXYPoints C struct that holds list of x,y pairs sorted by ascending x-values with an interpolation rule. That is, it contains the 1-d function f(x) = y(x). Functions are included that support common unitary and binary operations (e.g., negation, addition, multiplication). As example, if f(x) and g(x) are two ptwXYPoints instances, then routines to calculate f(x) + g(x) and f(x) * g(x) are included. The library support storage of data with lin-lin, lin-log, log-lin, log-log and flat (or histogram) interpolation. Functions to convert any of the supported interpolation to lin-lin are included. However, many of the binary functions only support lin-lin interpolation. The library also contains a C struct for storing a Legendre series representation of a 1-d function, and functions to convert Legendre series to and from a ptwXYPoints instance. A python interface to the library is included, with the ptwXYPoints C struct treated as a class. This class has many of the standard python math methods (e.g., __add__, __mul__).

Beck, Bret↗

DOE OSTI2020

eBPF Community ID Calculation

The code titled “eBPF Community ID Calculation” allows easy correlation of monitoring logs with existing network forensics logs on high-performance computing (HPC) systems. This code implements the Community ID Flow Hashing standard to calculate community ID Flow identifiers. These identifiers will be used to correlate host forensics logs with network forensics logs. The correlation provides valuable insight and the ability to easily pivot between disparate data sources. The computer code removes the libpcap dependency from the code base.

Maerz, StefanA↗

DOE OSTI2019

Diventi

Diventi provides a write optimized index of network logs such as netflow or zeek. The program will ingest various network traffic logs (zeek connection logs or NetflowV5, v9 or IPFix logs) and index each IP address and time of its activity. This database can then be queried for timely tracking of network activity. It currently uses the ft-index library from Percona but is adaptable to any key-value store that supports the Berkeley DB api.SAND2019-15037 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Kroeger, Thomas↗

DOE OSTI2019

Application of a convolutional neural network in permeability prediction: A case study in the Jacksonburg-Stringtown oil field, West Virginia, USA

Permeability is a critical parameter for understanding subsurface fluid flow behavior, managing reservoirs, enhancing hydrocarbon recovery, and sequestering carbon dioxide. In general, permeability is measured in the laboratory based on subsurface core samples, calculated from well logs or estimated from well tests. However, laboratory measurements and well tests are expensive, time-consuming, and usually limited to a few core samples or wells in a hydrocarbon field or carbon storage site. Machine-learning techniques are good options for generating a rapid, robust, and cost-effective permeability prediction model because of their strengths to recognize the potential interrelationships between input and output variables. Convolutional neural networks (CNN), as a good pattern recognition algorithm, are widely used in image processing, natural language processing, and speech recognition, but are rarely used with regression problems and even less often in reservoir characterization. We have developed a CNN regression model to estimate the permeability in the Jacksonburg-Stringtown oil field, West Virginia, which is a potential carbon storage site and enhanced oil recovery operations field. We also evaluate the concept of the geologic feature image, which is converted from geophysical well logs. Five variables, including two commonly available conventional well logs (the gamma rays [GRs] and bulk density) and three well-log-derived variables (the slopes of the GR and bulk density curves, and shale content), are used to generate a geologic feature image. The CNN treats the geologic feature image as the input and the permeability as the desired output. In addition, the permeability predicted using traditional backpropagation artificial neural networks, which are optimized by genetic algorithms and particle swarm optimization, is compared with the permeability estimated using our CNN. Our results indicate that the CNN regression model provides more accurate permeability predictions than the traditional neural network.

Geochemistry & Geophysics↗

DOE OSTI2023

Reference data, predictors, and probability grids for forest degradation classes in three sites in the Brazilian Amazon

Forest degradation by fires and selective logging is widespread in the Amazon region. We implemented a gradient boosted classification modeling framework to classify intact, logged, and burned forests at three Amazonian sites: Feliz Natal Municipality and Xingu Indigenous Territory in Mato Grosso State, and Saracá-Taquera National Forest in Pará State. We used forest degradation history from Landsat time-series as reference data and textural metrics derived from PlanetScope images as predictors. Textural metrics were computed using the Gray-Level Co-Occurrence Matrix (GLCM) textural technique. Included in the attached zip file are ten files: - a shapefile containing the reference data (fire and selective logging polygons and year of event) for each site; - a multiband tif file containing the 8 GLCM metrics used as predictors (Mean, Variance, Homogeneity, Contrast, Dissimilarity, Entropy, Angular Second Moment, Correlation) at the original PlanetScope resolution (3.125m) for each site; - a multiband tif file containing the 72 aggregated GLCM metrics used as predictors (Mean, Variance, Homogeneity, Contrast, Dissimilarity, Entropy, Angular Second Moment, and Correlation aggregated using the mean, first quartile, third quartile, maximum, median, minimum, root mean square, standard deviation, and skewness statistics) at 562m resolution for each site; - a multiband tif file containing the 3 probability grids for either intact, logged, or burned forests at the aggregation resolution (562m) for each site.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Field Laboratory for Emerging Stacked Unconventional Plays in Central Appalachia

The goal of the Field Laboratory for Emerging Stacked Unconventional Plays (ESUP) in Central Appalachia project was to investigate and characterize the resource potential for multi-play production of emerging unconventional reservoirs in Central Appalachia. Project activities included drilling, logging and coring of a vertical characterization well drilled to basement to approximately 15,000 feet in depth. The data from drilling and well logs as well recovered core samples were used by the research team to support characterization of the geology and potential pay zones within the Nora Gas Field of southwestern Virginia and the greater Central Appalachian Basin. This project was led by Dr. Nino Ripepi of the Virginia Center for Coal and Energy Research (VCCER) in close collaboration with EnerVest Operating LLC (EnerVest). VCCER is housed in the Department of Mining and Minerals Engineering at Virginia Tech and was created by an Act of the Virginia General Assembly on March 30, 1977, as an interdisciplinary study, research, information, and resource facility for the Commonwealth of Virginia. EnerVest is a top tier, low-cost oil and natural gas company with a long history of traditional operating relationships with institutional investors as well as public and private companies. This report summarizes two major research activities: 1) Log and Core Analysis and 2) Numerical Modeling. Core analysis comprises the characterization efforts of the ESUP Field Laboratory and includes a discussion of the well logs run in the deep characterization wells, an initial analysis of those, the number, size and location of cores recovered, core testing and analysis done to date as well as plans for future analysis. Numerical modeling is comprised of seven major modeling efforts, including (1) Multiphysics shale transport modeling work, (2) the developed in-house compositional simulator studies, (3) reservoir modeling studies using GEM software, (4) fracture modeling work using EFRAC3D software and (5) numerical modeling efforts with ABAQUS software, (6) Numerical modeling efforts using FLAC3D and PFC3D, and (7) Automatic Machine Learning (AutoML) Studies.

02 PETROLEUM↗

DOE OSTI2024

Risk of theft and malicious use of radiation sources in transit

Remote radioactive source applications require frequent transportation of sources from storage locations to remote sites. This introduces risk of theft of a source during the transportation process, with the level of risk proportional to the radioactivity of the source. To that end, theft of smaller sources, such as microcurie-level moisture density gauges, are of minor concern, but larger sources, such as those used for radiography and well logging, present more risk. Radiography sources include 192 Ir, 75 Se, or 60 Co radionuclides with radioactivity amounts at or exceeding IAEA Category 2. Well-logging sources, primarily 241 Am/Be, are used for their neutron-emission properties. 137 Cs is also used in well-logging at lower activities than in radiography but at levels that still present some risk. The vulnerability for malicious use of such sources to cause contamination and associated economic effects is dependent on the elemental chemical and physical properties, especially melting point and bulk modulus. Theft of radiography sources is somewhat common, well-logging sources less so. Theft of sources commonly occurs in concert with theft of the vehicle, with the source subsequently abandoned. There have been some instances where a source appears to have been specifically targeted. There are a variety of security measures and protocols, available and under development, to mitigate the risk of theft and assist in source recovery.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2022

Quantum simulation of real-space dynamics

Quantum simulation is a prominent application of quantum computers. While there is extensive previous work on simulating finite-dimensional systems, less is known about quantum algorithms for real-space dynamics. We conduct a systematic study of such algorithms. In particular, we show that the dynamics of a d-dimensional Schrödinger equation with η particles can be simulated with gate complexity O ~ (ηdFpoly(log(g'/ϵ))), where ϵ is the discretization error, g' controls the higher-order derivatives of the wave function, and F measures the time-integrated strength of the potential. Compared to the best previous results, this exponentially improves the dependence on ϵ and g' from poly(g'/ϵ) to poly(log(g'/ϵ)) and polynomially improves the dependence on T and d, while maintaining best known performance with respect to η. For the case of Coulomb interactions, we give an algorithm using η 3 (d + η)Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates, and another using η 3 (4d) d/2 Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates and QRAM operations, where T is the evolution time and the parameter Δ regulates the unbounded Coulomb interaction. We give applications to several computational problems, including faster real-space simulation of quantum chemistry, rigorous analysis of discretization error for simulation of a uniform electron gas, and a quadratic improvement to a quantum algorithm for escaping saddle points in nonconvex optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗