Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function approximation”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Lastly, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

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

A stochastic flow-capturing model to optimize the location of fast-charging stations with uncertain electric vehicle flows

Here, we develop a model to optimize the location of public fast charging stations for electric vehicles (EVs). A difficulty in planning the placement of charging stations is uncertainty in where EV charging demands appear. For this reason, we use a stochastic flow-capturing location model (SFCLM). A sample-average approximation method and an averaged two-replication procedure are used to solve the problem and estimate the solution quality. We demonstrate the use of the SFCLM using a Central-Ohio based case study. We find that most of the stations built are concentrated around the urban core of the region. As the number of stations built increases, some appear on the outskirts of the region to provide an extended charging network. We find that the sets of optimal charging station locations as a function of the number of stations built are approximately nested. We demonstrate the benefits of the charging-station network in terms of how many EVs are able to complete their daily trips by charging midday—six public charging stations allow at least 60% of EVs that would otherwise not be able to complete their daily tours without the stations to do so. We finally compare the SFCLM to a deterministic model, in which EV flows are set equal to their expected values. We show that if a limited number of charging stations are to be built, the SFCLM outperforms the deterministic model. As the number of stations to be built increases, the SFCLM and deterministic model select very similar station locations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

Classification of Energy States of the Exciton in Square Quantum Well

The energy states of the exciton in a single square GaAs-based quantum well are obtained from a numerical solution of the three-dimensional Schrödinger equation for the envelope of the exciton wave function. This equation is based on the exciton effective energy operator with a spherical approximation of the Luttinger Hamiltonian. The calculated states are classified based on the types of one-dimensional functions for the factorized form of the wave function. The upper limit for the energies of the exciton states in a quantum well is confirmed by the complex scaling method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Multicomponent mutual diffusion in the warm, dense matter regime

In this work, we present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Deciphering the small-angle scattering of polydisperse hard spheres using deep learning

We introduce a deep learning approach for analyzing the scattering function of the polydisperse hard sphere system. We use a variational autoencoder-based neural network to learn the bidirectional mapping between the scattering function and the system parameters, including the volume fraction and polydispersity. Such that the trained model serves both as a generator that produces a scattering function from the system parameters and an inferrer that extracts system parameters from the scattering function. We first generate a scattering dataset by carrying out molecular dynamics simulations of the polydisperse hard spheres modeled by the truncated-shifted Lennard-Jones model, then analyze the scattering function dataset using singular value decomposition to confirm the feasibility of dimensional compression. Then, we split the dataset into training and testing sets and train our neural network on the training set only. Our generator model produces a scattering function with significantly higher accuracy compared to the traditional Percus–Yevick approximation and β correction, and the inferrer model can extract the volume fraction and polydispersity with much higher accuracy than traditional model functions.

Ding, Lijie [ORNL] (ORCID:0000000227454606)↗
DOE OSTI2018

Impact of Tensor Interaction on Beta-Delayed Neutron Emission from Neutron-Rich Nickel Isotopes

The Skyrme energy density functional including tensor interaction is used to describe microscopically multineutron delayed-neutron emission accompanying beta decay of even–even neutron-rich nickel isotopes of mass number in the range of A = 74−80. The respective calculations are performed in the quasiparticle random-phase approximation with allowance for the two-phonon components of the wave function for states of the daughter nucleus. The properties of the lowest quadrupole excitation of {sup 74,76,78,80}Ni are also studied. It is shown that a decrease in the strength of neutron–proton tensor interaction leads to a substantial hindrance of beta decay and to an increase in the probability for delayedneutron emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A benchmark of Gutzwiller conjugate gradient minimization method in ground state energy calculations of dimers

Herein we present numerical results of ground-state energies of 9 molecules in the well-established G2 molecule set given by the Gutzwiller conjugate gradient minimization (GCGM) method. The method, beyond the commonly used Gutzwiller approximation, was recently developed based on Gutzwiller variational wave functions. We find that compared to benchmark data given by full configuration interaction, GCGM total energies are reasonably well reproduced with the minimum basis set. To include the dynamical correlation beyond the minimal basis calculations, we adopt the local density approximation for the dynamical correlation energy $E_c$. By comparing the results with benchmark data given by experiments and large-basis configuration interaction, the GCGM total energies with $E_c$ are in general better reproduced, but discrepancies are still observed for some dimers.

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

Integration of transcriptomic and proteomic data identifies biological functions in cell populations from human infant lung

Systems biology uses computational approaches to integrate diverse data types to understand cell and organ behavior. Data derived from complementary technologies, for example transcriptomic and proteomic analyses, are providing new insights into development and disease. We compared mRNA and protein profiles from purified endothelial, epithelial, immune, and mesenchymal cells from normal human lung tissue. Signatures for each cell type were identified and compared at both mRNA and protein levels. Cell specific biological processes and pathways were predicted by analysis of concordant and discordant RNA-protein pairs. Cell clustering and gene set enrichment comparisons identified shared versus unique processes associated with transcriptomic and/or proteomic data. Clear cell-cell correlations between mRNA and protein data were obtained from each cell type. Approximately 40% of RNA-protein pairs were coherently expressed. While the correlation between RNA and their protein products was relatively low (approximately 0.4), cell specific signature genes involved in functional processes characteristic of each cell type, were more highly correlated with their protein products. Consistency of cell specific RNA-protein signatures indicated an essential framework for the function of each cell type. Visualization and reutilization of the protein and RNA profiles are supported by a new web application, "LungProteomics," which is incorporated into the LGEA web portal and freely accessible at (https://research.cchmc.org/pbge/lunggens/lungProteomics/human_sorted_profiles.html).

human lung, sorted cells, integrated omics, Transc↗
DOE OSTI2022

How Well Does Kohn–Sham Regularizer Work for Weakly Correlated Systems?

Kohn–Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn–Sham density functional theory that works for strongly correlated systems. Here we test KSR for a weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, and Be 2+ and testing on 1D hydrogen chains, LiH, BeH 2 , and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 mH.

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

Block Krylov Subspace Methods for Functions of Matrices II: Modified Block FOM

We analyze an expansion of the generalized block Krylov subspace framework of [Electron. Trans. Numer. Anal., 47 (2017), pp. 100--126]. This expansion allows the use of low-rank modifications of the matrix projected onto the block Krylov subspace and contains, as special cases, the block GMRES method and the new block Radau--Arnoldi method. Within this general setting, we present results that extend the interpolation property from the nonblock case to a matrix polynomial interpolation property for the block case, and we relate the eigenvalues of the projected matrix to the latent roots of these matrix polynomials. Some error bounds for these modified block FOM methods for solving linear systems are presented. We then show how cospatial residuals can be preserved in the case of families of shifted linear block systems. This result is used to derive computationally practical restarted algorithms for block Krylov approximations that compute the action of a matrix function on a set of several vectors simultaneously. Finally, we prove some error bounds and present numerical results showing that two modifications of FOM, the block harmonic and the block Radau--Arnoldi methods for matrix functions, can significantly improve the convergence behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

First-principles calculation of excited states of diatomic molecules: a benchmark for the Gutzwiller conjugate gradient minimisation method

We recently proposed the Gutzwiller conjugate gradient minimisation (GCGM) method for efficient and accurate calculation of the ground state total energy of molecular and bulk systems. The GCGM method is developed under the framework of Gutzwiller wave function but goes beyond the commonly adopted Gutzwiller approximation to improve the accuracy and flexibility in treating the correlation effects. In this conference proceeding, we benchmark the GCGM method with the calculation of excited state potential energy curves of three diatomic molecules, namely H 2 , N 2 , and O 2 . Our calculations demonstrate the flexibility and reasonable accuracy of the method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
DOE OSTI2020

Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments

The new Bayesian Monte Carlo (MC) evaluation framework described in this abstract has been conceived as an attempt to improve nuclear data evaluations of differential crosssection data by removing the following two approximations conventionally employed for nuclear data evaluations: all probability density functions (PDFs) of all data and model parameters, both prior and posterior, are assumed to be normal (i.e., Gaussian) PDFs, and all uncertainties and covariances are propagated using a linear approximation. With these approximations removed, the Bayesian MC (BMC) framework could be used to account for nonlinear effects and would enable improved evaluations of differential cross sections and IBE data that are presently performed based on the assumptions itemized above. The BMC would also improve upon the uniform sampling of IBE parameters from within ranges defined by their evaluated uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

The quantum Hall effect at a weak magnetic field

Using a weak limit for the hopping integral in one direction in the Hofstadter model, we show that the fermion states in the gaps of the spectrum are determined within the Kitaev chain. The proposed approach allows us to study the behavior of Chern insulators (CI) in different classes of symmetry. We consider the Hofstadter model on the square and honeycomb lattices in the case of rational and irrational magnetic fluxes ϕ, and discuss the behavior of the Hall conductance at a weak magnetic field in a sample of finite size. We show that in the semiclassical limit at the center of the fermion spectrum, the Bloch states of fermions turn into chiral Majorana fermion liquid when the magnetic scale (1/ϕ) is equal to the sample size N. We are talking about the dielectric–metal phase transition, which is determined by the behavior of the Landau levels in 2D fermion systems in a transverse magnetic field. When a magnetic scale, which determines the wave function of fermions, exceeds the size of the sample, a jump in the longitudinal conductance occurs. The wave function describes non-localized states of fermions, the sample becomes a conductor, the system changes from the dielectric state to the metallic one. It is shown, that at 1∕ϕ>N the quantum Hall effect and the Landau levels are not realized, which makes possibility to study the behavior of CI in irrational magnetic fluxes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

On the validity of the guiding-center approximation in the presence of strong magnetic gradients

The motion of a charged particle in a nonuniform straight magnetic field with a constant magnetic-field gradient is solved exactly in terms of elliptic functions. The connection between this problem and the guiding-center approximation is discussed. Here, it is shown that, for this problem, the predictions of higher-order guiding-center theory agree very well with the orbit-averaged particle motion and hold well beyond the standard guiding-center limit $\epsilon$ $\equiv$ p/L << 1, where p is the gyromotion length scale and L is the magnetic-field gradient length scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

The Potential of DAS in Teleseismic Studies: Insights From the Goldstone Experiment

Distributed acoustic sensing (DAS) is a recently developed technique that has demonstrated its utility in the oil and gas industry. Here we demonstrate the potential of DAS in teleseismic studies using the Goldstone OpticaL Fiber Seismic experiment in Goldstone, California. By analyzing teleseismic waveforms from the 10 January 2018 M7.5 Honduras earthquake recorded on ~5,000 DAS channels and the nearby broadband station GSC, we first compute receiver functions for DAS channels using the vertical-component GSC velocity as an approximation for the incident source wavelet. The Moho P -to- s conversions are clearly visible on DAS receiver functions. Here, we then derive meter-scale arrival time measurements along the entire 20-km-long array. We are also able to measure path-averaged Rayleigh wave group velocity and local Rayleigh wave phase velocity. The latter, however, has large uncertainties. Our study suggests that DAS will likely play an important role in many fields of passive seismology in the near future.

58 GEOSCIENCES↗