Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate computing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,801 records · Page 100

The holar wind

String theory in AdS 3 with purely NS-NS fluxes and vanishing RR moduli has a continuum of winding string excitations in radial plane wave states. BTZ black holes can emit such strings, which then flow out toward the AdS 3 boundary as a stream of massive quanta, and form a black hole analogue of the solar wind. The winding string sector thus provides a decay channel for the black hole to evaporate without having either to couple the system to an external reservoir or to match the AdS 3 throat onto an asymptotically flat region. We compute the emission amplitude of this “holar wind” in the semi-classical approximation, and consider the associated version of the black hole information paradox.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optical lithography

Optical lithography is a photon-based technique comprised of projecting an image into a photosensitive emulsion (photoresist) coated onto a substrate such as a silicon wafer. It is the most widely used lithography process in the high volume manufacturing of nano-electronics by the semiconductor industry. Optical lithography’s ubiquitous use is a direct result of its highly parallel nature allowing vast amounts of information to be transferred very rapidly. For example, a modern leading edge lithography tool produces 150-300-mm patterned wafers per hour with 40-nm two-dimensional pattern resolution, yielding a pixel throughput of approximately 1.8T pixels/s. Continual advances in optical lithography capabilities have enabled the computing revolution over the past 50 years.

Naulleau, P↗

Data preservation at the Fermilab Tevatron

The Fermilab Tevatron collider's data-taking run ended in September 2011, yielding a dataset with rich scientific potential. The CDF and D0 experiments each have approximately 9 PB of collider and simulated data stored on tape. A large computing infrastructure consisting of tape storage, disk cache, and distributed grid computing for physics analysis with the Tevatron data is present at Fermilab. The Fermilab Run II data preservation project intends to keep this analysis capability sustained through the year 2020 and beyond. To achieve this goal, we have implemented a system that utilizes virtualization, automated validation, and migration to new standards in both software and data storage technology and leverages resources available from currently-running experiments at Fermilab. Lastly, these efforts have also provided useful lessons in ensuring long-term data access for numerous experiments, and enable high-quality scientific output for years to come.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Speeding Up Hartree–Fock in JuliaChem with Density Fitting

In this work, the density fitting (DF) approximation is added to the restricted Hartree–Fock (RHF) implementation in the JuliaChem computational chemistry code. Utilizing a DF algorithm that uses symmetry and integral screening, a significant reduction in time to compute the Fock matrix is achieved. The symmetry and screening DF-RHF techniques were adapted to be performed on graphics processing units (GPUs), which are well suited to perform the matrix multiplications that comprise the bulk of the Fock build time in DF-RHF. The JuliaChem DF-RHF GPU algorithm employs a novel approach that automatically switches between two DF-RHF algorithms depending on the number of basis functions in the calculation. The JuliaChem GPU DF-RHF implementation demonstrates up to 2× speedup for Fock build times compared to the existing best-in-class GPU DF-RHF implementation by operating directly on screened intermediate matrices. Due to the high portability of the Julia language code, the JuliaChem CPU and GPU DF-RHF implementations could be benchmarked on a variety of CPU and GPU architectures from multiple hardware vendors.

Hayes, John J. [Ames Laboratory, and Iowa State Un↗

Chemical Imaging of Atmospheric Biomass Burning Particles from North American Wildfires

The effects of biomass burning aerosols (BBA) on radiative forcing and cloud formation depend on chemical composition and the internal structures of individual particles within smoke plumes. To improve our understanding of the chemical and physical properties of BBA emitted at different times of the day and their evolution during atmospheric aging, we conducted a study as a part of the Fire Influence on Regional to Global Environments and Air Quality field campaign. Particle samples were collected onboard a research aircraft from smoke plumes from a wildfire in eastern Oregon during late afternoon and nighttime flights on August 28, 2019. A time-resolved aerosol collector was used to collect samples on substrates for offline spectromicroscopic imaging to investigate the single-particle characteristics of BBA particles. Approximately 20,400 individual particles from 10 selected samples were analyzed using computer-controlled scanning electron microscopy coupled with energy-dispersive X-ray microanalysis, revealing their elemental composition, morphology, and viscosity. Elemental microanalysis indicated that aged potassium is likely found in the form of K 2 SO 4 , KNO 3 , and possible K-organic salts. Further chemical speciation and carbon bonding mapping within individual particles were conducted using synchrotron-based scanning transmission X-ray microscopy (STXM) coupled with near edge X-ray absorption fine structure (NEXAFS) spectroscopy. Real-time, water-soluble light absorption measurements were acquired using a particle-into-liquid sampler instrument coupled to a liquid waveguide capillary cell and total organic analyzer. In the late afternoon samples, 65% of the total particle number population consisted entirely of organic components, compared to 46% in the nighttime particles. These differences were attributed to discrepancies in composition at the time of emission and to the daytime condensation and accumulation of photochemically formed secondary organic material on existing BBA particles, a process that halts at night. Microscopy images indicated that particle viscosity was lower in the nighttime particles (<10 1 Pa·s), likely due to increased relative humidity and a higher contribution from hygroscopic inorganic components. The chemical heterogeneity of individual particles was quantified using STXM-derived mixing state parameters. The nature of carbon bonding within individual particles was inferred from the extent of carbon sp 2 hybridization derived from NEXAFS spectra. Average percentages of sp 2 hybridization range between 40% and 60%, with no noticeable differences between late afternoon and nighttime flights. These findings were compared with the online optical properties of both late afternoon and nighttime smoke plumes, providing valuable insights into the complex relationship between chemical composition and optical properties of BBA particles at different times of the day.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of hydrogenation on the structure and magnetic properties of an iron oxide cluster

The structure and properties of the Fe8O12Hn clusters (n = 0–18) are computed using the all-electron density functional theory with the generalized gradient approximation for the exchange–correlation potential. The ground state of Fe8O12 is found to be a singlet state having a bi-capped octahedron geometry. Upon hydrogenation, the octahedral framework of Fe is retained in Fe8O12Hn up to n < 7, beyond which point the iron octahedron transforms into a cube. Hydrogen atoms are bound to oxygen atoms up to n = 12, but they bind to the faces of the Fe8 cube when n > 12. The total spin magnetic moment of a Fe8O12Hn cluster is larger than 6 μB for 1 ≤ n ≤ 18, except for n = 8 and 10, where the lowest total energy states are antiferromagnetic singlets. The reason for this deviation from the general behavior in the Fe8O12Hn series is attributed to the collective superexchange phenomenon. Surprisingly, the total spin magnetic moment of a Fe8O12Hn cluster is found to be substantially larger than the total spin magnetic moment of the bare Fe8 cluster when n = 12–18. All of the Fe8O12Hn clusters are stable with respect to an abstraction of a single hydrogen atom but are unstable toward the abstraction of an H2 dimer when n =10 and n = 14–18.

Gutsev, G. L.↗

A theoretical and experimental study of pressure broadening of the oxygen A-band by helium

The rotationally resolved magnetic dipole absorption spectrum of the oxygen A-band b{sup 1}Σ{sub g}{sup +}(v=0)←X{sup 3}Σ{sub g}{sup −}(v=0) perturbed by collisions with helium was studied theoretically using the impact approximation. To calculate the relaxation matrix, scattering calculations were performed on a newly computed helium-oxygen (b{sup 1}Σ{sub g}{sup +}) interaction potential as well as on a helium-oxygen (X{sup 3}Σ{sub g}{sup −}) interaction potential from the literature. The calculated integrated line cross sections and broadening coefficients are in good agreement with experimental results from the literature. Additionally, cavity ring-down experiments were performed in the wings of the spectral lines for a quantitative study of line-mixing, i.e., the redistribution of rotational line intensities by helium-oxygen collisions. It is shown that inclusion of line-mixing in the theory is required to reproduce the experimentally determined absolute absorption strengths as a function of the density of the helium gas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient manifolds tracing for planar maps

Here, we present an exact calculation method and an approximation technique for tracing the invariant manifolds of unstable periodic orbits of planar maps. The exact approach relies on an adaptive refinement procedure that prevents redundant calculations occurring in other approaches, and the approximated method relies on a novel interpolation approach based on normal displacement functions. The resulting approximated manifold is precise when compared to the exact one, and its relative computational cost falls like the inverse of the manifold length. To present the tracing method, we obtain the invariant manifolds of the Chirikov-Taylor map, and as an application we illustrate the transition from homoclinic to heteroclinic chaos in the Duffing oscillator that leads from localized chaos to global chaotic motion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Characterization of the 2-methylvinoxy radical + O 2 reaction: A focal point analysis and composite multireference study

Vinoxy radicals are involved in numerous atmospheric and combustion mechanisms. High-level theoretical methods have recently shed new light on the reaction of the unsubstituted vinoxy radical with O 2 . The reactions of 1-methylvinoxy radical and 2-methylvinoxy radical with molecular oxygen have experimental high pressure limiting rate constants, k ∞ , 5–7 times higher than that of the vinoxy plus O 2 reaction. In this work, high-level ab initio quantum chemical computations are applied to the 2-methylvinoxy radical plus O 2 system, namely, the formation and isomerization of the 1-oxo-2-propylperoxy radical, the immediate product of O 2 addition to the 2-methylvinoxy radical. Multireference methods were applied to the entrance channel. No barrier to O 2 addition could be located, and more sophisticated treatment of dynamic electron correlation shows that the principal difference between O 2 addition to the vinoxy and 2-methylvinoxy radicals is a larger steric factor for 2-methylvinoxy + O 2 . This is attributed to the favorable interaction between the incoming O 2 molecule and the methyl group of the 2-methylvinoxy radical. Via the focal point approach, energetics for this reaction were determined, in most cases, to chemical accuracy. The coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] correlation energy and Hartree–Fock energies were independently extrapolated to the complete basis set limit. A correction for the effect of higher excitations was computed at the CCSDT(Q)/6-31G level. Corrections for the frozen-core approximation, the Born–Oppenheimer approximation, the nonrelativistic approximation, and the zero-point vibrational energy were included. From the 1-oxo-2-propylperoxy radical, dissociation to reactants is competitive with the lowest energy isomerization pathway. The lowest energy isomerization pathway ultimately forms acetaldehyde, CO, and ·OH as the final products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗

Approximating large-basis coupled-cluster theory vibrational frequencies using focal-point approximations

The focal-point approximation can be used to estimate a high-accuracy, slow quantum chemistry computation by combining several lower-accuracy, faster computations. We examine the performance of focal-point methods by combining second-order Møller–Plesset perturbation theory (MP2) with coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] for the calculation of harmonic frequencies and that of fundamental frequencies using second-order vibrational perturbation theory (VPT2). In contrast to standard CCSD(T), the focal-point CCSD(T) method approaches the complete basis set (CBS) limit with only triple-ζ basis sets for the coupled-cluster portion of the computation. The predicted harmonic and fundamental frequencies were compared with the experimental values for a set of 20 molecules containing up to six atoms. Here, the focal-point method combining CCSD(T)/aug-cc-pV(T + d)Z with CBS-extrapolated MP2 has mean absolute errors vs experiment of only 7.3 cm −1 for the fundamental frequencies, which are essentially the same as the mean absolute error for CCSD(T) extrapolated to the CBS limit using the aug-cc-pV(Q + d)Z and aug-cc-pV(5 + d)Z basis sets. However, for H 2 O, the focal-point procedure requires only 3% of the computation time as the extrapolated CCSD(T) result, and the cost savings will grow for larger molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange

We report on the use of the adaptively compressed exchange (ACE) operator to accelerate many-body perturbation theory (MBPT) calculations, including G 0 W 0 and the Bethe–Salpeter equation (BSE), for hybrid density functional theory starting points. We show that by approximating the exact exchange operator with the low-rank ACE operator, substantial computational savings can be achieved with systematically controllable errors in the quasiparticle energies computed with full-frequency G 0 W 0 and the optical absorption spectra and vertical excitation energies computed by solving the BSE within density matrix perturbation theory. Our implementation makes use of the ACE-accelerated electronic Hamiltonian to carry out both G 0 W 0 and BSE without explicitly computing empty states. We show the robustness of the approach and present the computational gains obtained on both the central processing unit and graphics processing unit nodes. In conclusion, our work will facilitate the exploration and evaluation of fine-tuned hybrid starting points aimed at enhancing the accuracy of MBPT calculations without involving computationally demanding self-consistency in Hedin’s equations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of reaction schemes using maximum rates of constituent steps

Significance The design of active and selective catalysts is essential for a wide range of industrial applications and societal issues. A fundamental approach is to identify the key elementary steps and to elucidate the predicted reaction kinetics for potential reaction mechanisms. We present a methodology to analyze analytically the performance of catalytic reaction schemes by calculation of the maximum rates of the constituent steps. This proposed methodology can be used to identify the important transition states and adsorbed species, such that more detailed calculations can be carried out for these species, whereas more approximate methods can be used for the remaining species, thereby substantially reducing the computational time required to elucidate how catalyst performance is controlled by the fundamental surface chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast matrix algebra for Bayesian model calibration

In Bayesian model calibration, evaluation of the likelihood function usually involves finding the inverse and determinant of a covariance matrix. When Markov Chain Monte Carlo (MCMC) methods are used to sample from the posterior, hundreds of thousands of likelihood evaluations may be required. In this paper, we demonstrate that the structure of the covariance matrix can be exploited, leading to substantial time savings in practice. Here, we also derive two simple equations for approximating the inverse of the covariance matrix in this setting, which can be computed in near-quadratic time. The practical implications of these strategies are demonstrated using a simple numerical case study and the "quack" R package. For a covariance matrix with 1000 rows, application of these strategies for a million likelihood evaluations leads to a speedup of roughly 4000 compared to the naive implementation

97 MATHEMATICS AND COMPUTING↗

DYNAMO EFFECTS IN MAGNETOROTATIONAL TURBULENCE WITH FINITE THERMAL DIFFUSIVITY

We investigate the saturation level of hydromagnetic turbulence driven by the magnetorotational instability in the case of vanishing net flux. Motivated by a recent paper of Bodo et al., we here focus on the case of a non-isothermal equation of state with constant thermal diffusivity. The central aim of the paper is to complement the previous result with closure parameters for mean-field dynamo models, and to test the hypothesis that the dynamo is affected by the mode of heat transport. We perform computer simulations of local shearing-box models of stratified accretion disks with approximate treatment of radiative heat transport, which is modeled via thermal conduction. We study the effect of varying the (constant) thermal diffusivity, and apply different vertical boundary conditions (BCs). In the case of impenetrable vertical boundaries, we confirm the transition from mainly conductive to mainly convective vertical heat transport below a critical thermal diffusivity. This transition is however much less dramatic when more natural outflow BCs are applied. Similarly, the enhancement of magnetic activity in this case is less pronounced. Nevertheless, heating via turbulent dissipation determines the thermodynamic structure of accretion disks and clearly affects the properties of the related dynamo. This effect, however, may have been overestimated in previous work, and a careful study of the role played by boundaries will be required.

79 ASTRONOMY AND ASTROPHYSICS↗

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages—conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced, and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber–Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. Here, we discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mixing bispectrum multipoles under geometric distortions

We derive general expressions for how the Alcock–Paczynski distortions affect the power spectrum and the bispectrum of cosmological fields. We compute explicit formulas for the mixing coefficients of bispectrum multipoles in the linear approximation. The leading-order effect for the bispectrum is the uniform dilation of all three wavevectors. The mixing coefficients depend on the shape of the bispectrum triplet. Our results for the bispectrum multipoles are framed in terms of the ‘natural’ basis of the lengths of three wavevectors but can be easily generalized for other bases and reduction schemes. Our validation tests confirm that the linear approximation is extremely accurate for all power spectrum multipoles. The linear approximation is accurate for the bispectrum monopole but results in sub-per cent level inaccuracies for the bispectrum quadrupole and fails for the bispectrum hexadecapole. Our results can be used to simplify the analysis of the bispectrum from galaxy surveys, especially the measurement of the baryon acoustic oscillation peak position. They can be used to replace numeric schemes with exact analytical formulae.

79 ASTRONOMY AND ASTROPHYSICS↗

Black Hole Airy Tail

In Jackiw-Teitelboim (JT) gravity, which is dual to a random matrix ensemble, the annealed entropy differs from the quenched entropy at low temperatures and goes negative. However, computing the quenched entropy in JT gravity requires a replica limit that is poorly understood. To circumvent this, we define an intermediate quantity called the semiquenched entropy, which has the positivity properties of the quenched entropy, while requiring a much simpler replica trick. We compute this in JT gravity in different regimes using (i) a bulk calculation involving wormholes corresponding to the Airy limit of the dual matrix integral and (ii) a boundary calculation involving one-eigenvalue instanton saddles proposed by Hernández-Cuenca, demonstrating consistency between these two calculations in their common regime of validity. We also clarify why similar one-eigenvalue instanton saddles cannot be used to compute the quenched entropy due to a breakdown of the saddle-point approximation for the one-eigenvalue instanton in the replica limit. Our results show how to use the gravitational path integral to prove that black holes in JT gravity have isolated ground states and to study their properties.

Gauge-gravity dualities↗

Refine your search to explore more results.