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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗
DOE OSTI2018

Gradient-based optimization for regression in the functional tensor-train format

Predictive analysis of complex computational models, such as uncertainty quantification (UQ), must often rely on using an existing database of simulation runs. In this paper we consider the task of performing low-multilinear-rank regression on such a database. Specifically we develop and analyze an efficient gradient computation that enables gradient-based optimization procedures, including stochastic gradient descent and quasi-Newton methods, for learning the parameters of a functional tensor-train (FT). We compare our algorithms with 22 other nonparametric and parametric regression methods on 10 real-world data sets and show that for many physical systems, exploiting low-rank structure facilitates efficient construction of surrogate models. Here, we use a number of synthetic functions to build insight into behavior of our algorithms, including the rank adaptation and group-sparsity regularization procedures that we developed to reduce overfitting. Finally we conclude the paper by building a surrogate of a physical model of a propulsion plant on a naval vessel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

FP-LAPW study of the fundamental properties of the cubic spinel CdAl{sub 2}O{sub 4}

Highlights: Black-Right-Pointing-Pointer Fundamental properties of CdAl{sub 2}O{sub 4} are investigated. Black-Right-Pointing-Pointer Results obtained for energy band gaps using EV-GGA are larger than that within PBE-GGA. Black-Right-Pointing-Pointer The decomposition of the dielectric function into individual band-to-band contributions is calculated. Black-Right-Pointing-Pointer The effective charge-carrier masses are estimated from the band structure. Black-Right-Pointing-Pointer Pressure and temperature dependences of some macroscopic parameters are obtained. -- Abstract: We have investigated the structural, elastic, electronic, optical and thermodynamic properties of the cubic spinel CdAl{sub 2}O{sub 4} using accurate ab initio calculations. Computed equilibrium structural parameters are in good agreement with the available experimental data. Single-crystals elastic parameters are calculated for pressure up to 30 GPa using a conserving-volume total energy-strain method. Isotropic elastic parameters for ideal polycrystalline CdAl{sub 2}O{sub 4} aggregates are computed in the framework of the Voigt-Reuss-Hill approximation. Result for band structure using the Engel-Vosko scheme of the GGA shows a significant improvement over the common GGA functionals. Optical spectra have been calculated for the energy range 0-30 eV. The peaks and structures in the optical spectra are assigned to interband transitions. Pressure dependence of the band gaps, static dielectric constant and static refractive index are also investigated. Pressure and thermal effects on some macroscopic properties are predicted using the quasi-harmonic Debye model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2025

Recent Improvements to the NWChem COSMO Module

This work presents recent developments in the COSMO solvation model implementation in NWChem. A new cavity construction approach, based on the solvent-excluding surface (SES) and utilizing the well-established GEPOL algorithm, has been introduced. Additionally, a straightforward procedure to merge surface segments that are too close─often a source of numerical artifacts─has been implemented. The available methods for correcting outlying surface charges have also been reviewed and improved. To validate the new implementation, we computed dielectric solvation energies for a chemically diverse set of approximately 100 molecules, including neutral species, small ions, and common ionic liquid components. Results were compared to those from GAMESS using the double-cavity method as a reference. Although the double-cavity approach can be regarded as more accurate, the simpler correction schemes available in NWChem─based on scaling factors or Lagrange multipliers─can achieve excellent agreement if the potential is also properly corrected, with mean unsigned deviations of around 0.15 kcal/mol. Predictions of typical vapor–liquid and liquid–liquid equilibria using a COSMO-SAC variant based on NWChem also showed very promising results.

Cavities↗
DOE OSTI2022

IMPROVING THE ACCURACY OF COMPOSITE METHODS: A G4MP2 METHOD WITH G4-LIKE ACCURACY AND IMPLICATIONS FOR MACHINE LEARNING

G4MP2 theory has proven to be a reliable and accurate quantum chemical composite method for the calculation of molecular energies using an approximation based on second-order perturbation theory to lower computational costs compared to G4 theory. However, it has been found to have significantly increased errors when applied to larger organic molecules with 10 or more nonhydrogen atoms. We report here on an investigation of the cause of the failure of G4MP2 theory for such larger molecules. One source of error is found to be the "higher-level correction (HLC)", which is meant to correct for deficiencies in correlation contributions to the calculated energies. This is because the HLC assumes that the contribution is independent of the element and the type of bonding involved, both of which become more important with larger molecules. We address this problem by adding an atom-specific correction, dependent on atom type but not bond type, to the higher-level correction. We find that a G4MP2 method that incorporates this modification of the higher-level correction, referred to as G4MP2A, becomes as accurate as G4 theory (for computing enthalpies of formation) for a test set of molecules with less than 10 nonhydrogen atoms as well as a set with 10-14 such atoms, the set of molecules considered here, with a much lower computational cost. The G4MP2A method is also found to significantly improve ionization potentials and electron affinities. Finally, we implemented the G4MP2A energies in a machine learning method to predict molecular energies.

Dandu, Naveen↗
DOE OSTI2025

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

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

Nuclear quantum effects enter the mainstream

Atomistic simulations of chemical, biological and materials systems have become increasingly precise and predictive owing to the development of accurate and efficient techniques that describe the quantum mechanical behaviour of electrons. Nevertheless, the overwhelming majority of such simulations still assumes that the nuclei behave as classical particles. Historically, this approximation could sometimes be justified owing to the complexity and computational overhead. However, neglecting nuclear quantum effects has become one of the largest sources of error, especially when systems containing light atoms are treated using current state-of-the-art descriptions of chemical interactions. Over the past decade, this realization has spurred a series of methodological advances that have dramatically reduced the cost of including these important physical effects in the structure and dynamics of chemical systems. We discuss how these developments are now allowing nuclear quantum effects to become a mainstream feature of molecular simulations. These advances have led to new insights into phenomena that are relevant to different areas of science — from biochemistry to condensed matter — and open the door to many exciting future opportunities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Time-dependent density functional theory quantum transport simulation in non-orthogonal basis

Basing on the earlier works on the hierarchical equations of motion for quantum transport, we present in this paper a first principles scheme for time-dependent quantum transport by combining time-dependent density functional theory (TDDFT) and Keldysh's non-equilibrium Green's function formalism. This scheme is beyond the wide band limit approximation and is directly applicable to the case of non-orthogonal basis without the need of basis transformation. The overlap between the basis in the lead and the device region is treated properly by including it in the self-energy and it can be shown that this approach is equivalent to a lead-device orthogonalization. This scheme has been implemented at both TDDFT and density functional tight-binding level. Simulation results are presented to demonstrate our method and comparison with wide band limit approximation is made. Finally, the sparsity of the matrices and computational complexity of this method are analyzed.

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

Structural, vibrational, and quasiparticle band structure of 1,1-diamino-2,2-dinitroethelene from ab initio calculations

The effects of pressure on the structural and vibrational properties of the layered molecular crystal 1,1-diamino-2,2-dinitroethelene (FOX-7) are explored by first principles calculations. We observe significant changes in the calculated structural properties with different corrections for treating van der Waals interactions to Density Functional Theory (DFT), as compared with standard DFT functionals. In particular, the calculated ground state lattice parameters, volume and bulk modulus obtained with Grimme's scheme, are found to agree well with experiments. The calculated vibrational frequencies demonstrate the dependence of the intra and inter-molecular interactions on FOX-7 under pressure. In addition, we also found a significant increment in the N–H...O hydrogen bond strength under compression. This is explained by the change in bond lengths between nitrogen, hydrogen, and oxygen atoms, as well as calculated IR spectra under pressure. Finally, the computed band gap is about 2.3 eV with generalized gradient approximation, and is enhanced to 5.1 eV with the GW approximation, which reveals the importance of performing quasiparticle calculations in high energy density materials.

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

Density-functional errors in ionization potential with increasing system size

This work investigates the effects of molecular size on the accuracy of density-functional ionization potentials for a set of 28 hydrocarbons, including series of alkanes, alkenes, and oligoacenes. As the system size increases, delocalization error introduces a systematic underestimation of the ionization potential, which is rationalized by considering the fractional-charge behavior of the electronic energies. The computation of the ionization potential with many density-functional approximations is not size-extensive due to excessive delocalization of the incipient positive charge. While inclusion of exact exchange reduces the observed errors, system-specific tuning of long-range corrected functionals does not generally improve accuracy. These results emphasize that good performance of a functional for small molecules is not necessarily transferable to larger systems.

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

Theory of molecule metal nano-particle interaction: Quantum description of plasmonic lasing

The recent quantum description of a few molecules interacting with plasmon excitations of a spherical metal nano-particle (MNP) as presented in the work of Zhang and May [Phys. Rev. B 89, 245441 (2014)] is extended to systems with up to 100 molecules. We demonstrate the possibility of multiple plasmon excitation and describe their conversion into far-field photons. The calculation of the steady-state photon emission spectrum results in an emission line-narrowing with an increasing number of molecules coupled to the MNP. This is considered as an essential criterion for the action of the molecule-MNP system as a nano-laser. To have exact results for systems with up to 20 molecules, we proceed as recently described by Richter et al. [Phys. Rev. B 91, 035306 (2015)] and study a highly symmetric system. It assumes an equatorial and regular position of identical molecules in such a way that their coupling is dominated by that to a single MNP dipole-plasmon excitation. Changing from the exact computation of the system’s complete density matrix to an approximate theory based on the reduced plasmon density matrix, systems with more than 100 molecules can be described. Finally, nonlinear rate equations are proposed which reproduce the mean number of excited plasmons in their dependence of the number of molecules and of the used pump rate. The second order intensity correlation function of emitted photons is related to the respective plasmon correlation function which approaches unity when the system starts lasing.

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

Theoretical study of radiative electron attachment to CN, C{sub 2}H, and C{sub 4}H radicals

A first-principle theoretical approach to study the process of radiative electron attachment is developed and applied to the negative molecular ions CN{sup −}, C{sub 4}H{sup −}, and C{sub 2}H{sup −}. Among these anions, the first two have already been observed in the interstellar space. Cross sections and rate coefficients for formation of these ions by direct radiative electron attachment to the corresponding neutral radicals are calculated. For the CN molecule, we also considered the indirect pathway, in which the electron is initially captured through non-Born-Oppenheimer coupling into a vibrationally resonant excited state of the anion, which then stabilizes by radiative decay. We have shown that the contribution of the indirect pathway to the formation of CN{sup −} is negligible in comparison to the direct mechanism. The obtained rate coefficients for the direct mechanism at 30 K are 7 × 10{sup −16} cm{sup 3}/s for CN{sup −}, 7 × 10{sup −17} cm{sup 3}/s for C{sub 2}H{sup −}, and 2 × 10{sup −16} cm{sup 3}/s for C{sub 4}H{sup −}. These rates weakly depend on temperature between 10 K and 100 K. The validity of our calculations is verified by comparing the present theoretical results with data from recent photodetachment experiments.

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

A one-dimensional study of the evolution of the microwave breakdown in air

The microwave breakdown in air is simulated numerically within a simple 1D model taking into account a perturbation of electromagnetic field by plasma. The simulations were performed using two qualitatively different codes. One of these codes is based on computation of Maxwell equations, whereas the other one utilizes an approximation of quasi-monochromatic electromagnetic field. There is a good agreement between simulation results obtained by using both codes. Calculations have been carried out in a wide range of air pressures and field frequencies; also varied were initial spatial distributions of plasma density. The results reveal strong dependence of the breakdown evolution on the relation between the field frequency and the gas pressure as well as on the presence of extended rarefied background plasma. At relatively low gas pressures (or high field frequencies), the breakdown process is accompanied by the stationary ionization wave propagating towards the incident electromagnetic wave. In the case of a high gas pressure (or a relatively low field frequency), the peculiarities of the breakdown are associated with a formation of plasma filament array. The extended background plasma can suppress formation of the plasma filament array completely even at high pressures (or low frequencies)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear bioheat transfer models and multi-objective numerical optimization of the cryosurgery operations

Numerical computation of the three dimensional problem of the freezing interface propagation during the cryosurgery coupled with the multi-objective optimization methods is used in order to improve the efficiency and safety of the cryosurgery operations performing. Prostate cancer treatment and cutaneous cryosurgery are considered. The heat transfer in soft tissue during the thermal exposure to low temperature is described by the Pennes bioheat model and is coupled with an enthalpy method for blurred phase change computations. The finite volume method combined with the control volume approximation of the heat fluxes is applied for the cryosurgery numerical modeling on the tumor tissue of a quite arbitrary shape. The flux relaxation approach is used for the stability improvement of the explicit finite difference schemes. The method of the additional heating elements mounting is studied as an approach to control the cellular necrosis front propagation. Whereas the undestucted tumor tissue and destucted healthy tissue volumes are considered as objective functions, the locations of additional heating elements in cutaneous cryosurgery and cryotips in prostate cancer cryotreatment are considered as objective variables in multi-objective problem. The quasi-gradient method is proposed for the searching of the Pareto front segments as the multi-objective optimization problem solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Nuclear magnetic relaxation by the dipolar EMOR mechanism: Three-spin systems

In aqueous systems with immobilized macromolecules, including biological tissue, the longitudinal spin relaxation of water protons is primarily induced by exchange-mediated orientational randomization (EMOR) of intra- and intermolecular magnetic dipole-dipole couplings. Starting from the stochastic Liouville equation, we have developed a non-perturbative theory that can describe relaxation by the dipolar EMOR mechanism over the full range of exchange rates, dipole couplings, and Larmor frequencies. Here, we implement the general dipolar EMOR theory for a macromolecule-bound three-spin system, where one, two, or all three spins exchange with the bulk solution phase. In contrast to the previously studied two-spin system with a single dipole coupling, there are now three dipole couplings, so relaxation is affected by distinct correlations as well as by self-correlations. Moreover, relaxation can now couple the magnetizations with three-spin modes and, in the presence of a static dipole coupling, with two-spin modes. As a result of this complexity, three secondary dispersion steps with different physical origins can appear in the longitudinal relaxation dispersion profile, in addition to the primary dispersion step at the Larmor frequency matching the exchange rate. Furthermore, and in contrast to the two-spin system, longitudinal relaxation can be significantly affected by chemical shifts and by the odd-valued (“imaginary”) part of the spectral density function. We anticipate that the detailed studies of two-spin and three-spin systems that have now been completed will provide the foundation for developing an approximate multi-spin dipolar EMOR theory sufficiently accurate and computationally efficient to allow quantitative molecular-level interpretation of frequency-dependent water-proton longitudinal relaxation data from biophysical model systems and soft biological tissue.

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

An experimental and theoretical study of the C 1s ionization satellites in CH 3 I

We report the C 1s ionization spectrum of CH 3 I has been studied both experimentally and theoretically. Synchrotron radiation has been employed to record polarization dependent photoelectron spectra at a photon energy of 614 eV. These spectra encompass the main-line due to the C 1s single-hole state and the peaks associated with the shake-up satellites. Vertical ionization energies and relative photoelectron intensities have been computed using the fourth-order algebraic-diagrammatic construction approximation scheme for the one-particle Green’s function and the 6-311++G** basis set. The theoretical spectrum derived from these calculations agrees qualitatively with the experimental results, thereby allowing the principal spectral features to be assigned. According to our calculations, two 2 A 1 shake-up states of the C 1s -1 σ CI → σ$*\atop{CI}$ type with singlet and triplet intermediate coupling of the electron spins (S' = 0, 1) play an important role in the spectrum and contribute significantly to the overall intensity. Both of these states are expected to have dissociative diabatic potential energy surfaces with respect to the C–I separation. Whereas the upper of these states perturbs the manifold of Rydberg states, the lower state forms a band which is characterized by a strongly increased width. Our results indicate that the lowest shake-up peak with significant spectral intensity is due to the pair (S' = 0, 1) of 2 E (C 1s -1 I 5p → σ$*\atop{CI}$) states. We predict that these 2 E states acquire photoelectron intensity due to spin-orbit interaction. Lastly, such interactions play an important role here due to the involvement of the I 5p orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗