Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Electronic structure, magnetic properties, and pairing tendencies of the copper-based honeycomb lattice Na 2 Cu 2 TeO 6

Spin-1/2 chains with alternating antiferromagnetic (AFM) and ferromagnetic (FM) couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density-matrix renormalization-group (DMRG) methods, we have systematically studied the dimerized chain system Na 2 Cu 2 TeO 6 with a d 9 electronic configuration. Near the Fermi level, in the nonmagnetic phase the dominant states are mainly contributed by the Cu 3d x 2 -y 2 orbitals highly hybridized with the O 2p orbitals, leading to an “effective” single-orbital low-energy model. By calculating the relevant hoping amplitudes, we explain the size and sign of the exchange interactions in Na 2 Cu 2 TeO 6 . In addition, a single-orbital Hubbard model is constructed for this dimerized chain system where the quantum fluctuations are taken into account. Both AFM and FM couplings (leading to an ↑|-|↓-|↓-↑| state) along the chain were found in our DMRG and Lanczos calculations, in agreement with density functional theory and neutron-scattering results. The hole pairing binding energy ΔE is predicted to be negative at Hubbard U~ 11eV, suggesting incipient pairing tendencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2016

Modified Spectral Fatigue Methods for S-N Curves With MIL-HDBK-5J Coefficients

The rainflow method is used for counting fatigue cycles from a stress response time history, where the fatigue cycles are stress-reversals. The rainflow method allows the application of Palmgren-Miner's rule in order to assess the fatigue life of a structure subject to complex loading. The fatigue damage may also be calculated from a stress response power spectral density (PSD) using the semi-empirical Dirlik, Single Moment, Zhao-Baker and other spectral methods. These methods effectively assume that the PSD has a corresponding time history which is stationary with a normal distribution. This paper shows how the probability density function for rainflow stress cycles can be extracted from each of the spectral methods. This extraction allows for the application of the MIL-HDBK-5J fatigue coefficients in the cumulative damage summation. A numerical example is given in this paper for the stress response of a beam undergoing random base excitation, where the excitation is applied separately by a time history and by its corresponding PSD. The fatigue calculation is performed in the time domain, as well as in the frequency domain via the modified spectral methods. The result comparison shows that the modified spectral methods give comparable results to the time domain rainflow counting method.

Irvine, Tom↗
NASA NTRS1983

Measurement of macroscopic plasma parameters with a radio experiment: Interpretation of the quasi-thermal noise spectrum observed in the solar wind

The ISEE-3 SBH radio receiver has provided the first systematic observations of the quasi-thermal (plasma waves) noise in the solar wind plasma. The theoretical interpretation of that noise involves the particle distribution function so that electric noise measurements with long antennas provide a fast and independent method of measuring plasma parameters: densities and temperatures of a two component (core and halo) electron distribution function have been obtained in that way. The polarization of that noise is frequency dependent and sensitive to the drift velocity of the electron population. Below the plasma frequency, there is evidence of a weak noise spectrum with spectral index -1 which is not yet accounted for by the theory. The theoretical treatment of the noise associated with the low energy (thermal) proton population shows that the moving electrical antenna radiates in the surrounding plasma by Carenkov emission which becomes predominant at the low frequencies, below about 0.1 F sub P.

Couturier, P.↗
DOE OSTI2022

A reforecasting-based dynamic reserve estimation for variable renewable generation and demand uncertainty

The installed capacity of renewables-based energy sources has been increasing in traditional power systems. In order to accommodate the increased variability and uncertainty associated with the deeper penetration of renewable sources like solar and wind, adjusted amounts of dynamic reserve are needed. Although probabilistic dynamic reserve estimation methods have been previously developed, most of them consider the uncertainty to be represented by parametric density functions that tend to perform poorly under extreme events and, moreover, neglect uncertainty introduced by the forecasting model itself. Toward addressing these limitations, this work presents, for the first time, a dynamic reserve estimation method for flexibility that incorporates nonparametric density estimation and a machine learning based reforecasting to provide a day-ahead prediction of the mean and spread of uncertainty around the base forecast. The prediction is, in turn, used to estimate the up and down reserve relative to the base forecast. Here, the present method takes various endogenous and exogenous features, including the calendar variables, as input to estimate the day-ahead reserve. Using a combination of reforecasting and dynamic reserve estimation techniques, the method is shown to adjust better to the dynamic nature of reserve requirements providing only what is needed to accommodate the expected deviations. Considering California Independent System Operator (CAISO) solar, wind and load data over an 18 month period, up to 67% reduction in the amount of reserve capacity needed for a one day reserve and reserve penalty for solar uncertainty is demonstrated. Additionally, the risk of reserve insufficiency in meeting the net demand is reduced by 20% with the proposed method.

14 SOLAR ENERGY↗
DOE OSTI2019

Importance of self-interaction-error removal in density functional calculations on water cluster anions

Accurate description of the excess charge in water cluster anions is challenging for standard semi-local and (global) hybrid density functional approximations (DFAs). Using the recent unitary invariant implementation of the Perdew–Zunger self-interaction correction (SIC) method using Fermi-Löwdin orbitals, we assess the effect of self-interaction error on the vertical detachment energies of water cluster anions with the local spin density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation, and the strongly constrained and appropriately normed (SCAN) meta-GGA functionals. Our results show that for the relative energies of isomers with respect to reference CCSD(T) values, the uncorrected SCAN functional has the smallest deviation of 21 meV, better than that for the MP2 method. The performance of SIC-SCAN is comparable to that of MP2 and is better than SIC-LSDA and SIC-PBE, but it reverses the ordering of the two lowest isomers for water hexamer anions. Removing self interaction error (SIE) corrects the tendency of LSDA, PBE, and SCAN to over-bind the extra electron. The vertical detachment energies (VDEs) of water cluster anions, obtained from the total energy differences of corresponding anion and neutral clusters, are significantly improved by removing self-interaction and are better than the hybrid B3LYP functional, but fall short of MP2 accuracy. Removing SIE results in substantial improvement in the position of the eigenvalue of the extra electron. The negative of the highest occupied eigenvalue after SIC provides an excellent approximation to the VDE, especially for SIC-PBE where the mean absolute error with respect to CCSD(T) is only 17 meV, the best among all approximations compared in this work.

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

Calculations of anisotropic magnetic properties using spin-orbit energy variations

We analyze several methods of obtaining the accurate relativistic total energy (TE) variations using traditional perturbation theories (PTs) and proposed coupling constant integration (CCI) methods. For this purpose, we perform benchmark calculations within the density functional theory taking the spin-orbit coupling (SOC) and its derivative as a perturbation. The TE change due to SOC addition obtained from both PTs and CCI is shown to reach the accuracy of fully self-consistent TE calculations. Similar accuracy is also obtained even for the magnetocrystalline anisotropy energy (MAE). The real advantage of the proposed methods is to use PTs and CCI methods in those electronic structure methods where accurate total energies currently cannot be obtained with required accuracy. Correspondingly, we demonstrate the applicability of suggested methods for calculations of MAE in different magnetic materials using a dynamic mean-field method. All suggested PTs and CCI methods also provide convenient site, orbital, and spin decompositions of the TE variation, creating a powerful way to analyze microscopic physics in strongly correlated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Performance of the spin-component-scaled methods for energy bands

The performance of various spin-component-scaled parameterisations is examined for the second-order many-body Green's-function [MBGF(2)] calculations of valence energy bands, taking three of the experimentally well-characterised polymers as examples: polyethylene, polytetrafluoroethylene, and polyacetylene. The parameterisations considered are Grimme's original SCS parameter set, Jung et al.'s original SOS set (retaining the opposite-spin component only), Śmiga et al.'s SCS(IP) set (calibrated specifically for ionization energies), and Śmiga et al.'s SOS(IP) set (calibrated for ionization energies with the opposite-spin component only; implicit in the os-D2 model of Opoku et al.). The SCS(IP) and SOS(IP) parameterisations are found to shift both outer and inner valence bands by up to a few electronvolts away from the experimental data. The original SCS and SOS parameter sets do not improve upon, but largely maintain the accuracy of the unscaled MBGF(2) methods. Given that the SOS-MBGF(2) method can be implemented in a quartic-scaling algorithm (for all roots), it is most promising for solid-state applications. Furthermore this observation is consistent with the success of the quartic-scaling GW methods without the vertex correction based on a density-functional theory reference.

Green’s-function theory↗
NASA NTRS2023

Efficient Calibration of Expensive Computational Models

Accounting for uncertainty when calibrating expensive computational models is a common challenge faced by scientists and engineers. Often Bayesian techniques are adopted to estimate a probability density function over the model parameters given noisy empirical data. The methods used to perform this type of probabilistic calibration are computationally prohibitive in that they require a large number of evaluations of the expensive model. In these cases, surrogate modeling -- that is, using a fast-to-evaluate, lower fidelity stand-in for the original computational model -- may be the only option to alleviate this computational burden. However, the upfront cost of generating training data to build a surrogate model can itself be expensive. As such, it is important to be judicious when selecting training points at which the full-fidelity model is evaluated. Here, an active learning approach is proposed that enables efficient selection of training points using approximate samples of the calibrated parameter probability density function. In this way, the training points can be concentrated in regions where the calibration algorithm requires high model accuracy.

active learning↗
DOE OSTI2025

Revisiting the assignment of atomic charges in metal oxides based on core-level x-ray photoelectron spectra: The case of Ti in SrTiO3(001)

We demonstrate that assigning formal charges to transition metal (TM) cations based on core-level (CL) x-ray photoemission binding energies in oxides leads to physically inconsistent pictures of electronic structure. O 2p–TM 3d hybridization is well known to result in significant covalency in TM–O bonds, thereby reducing TM cation charges from their fully ionic values. However, the ionic bonding model remains the working paradigm for assigning TM CL features, and the resulting cation charges are often taken to be representative of the material under study. Here, we show that a more physically meaningful way to assign charges is to extract information about charge distributions utilizing Dirac–Hartree–Fock theory to calculate CL spectra from first principles and then use the resulting wave functions to determine charges based on orbital occupancies. TM cation charges can also be determined using density functional theory and Bader population analysis. We illustrate these two methods using the Ti 2p spectrum for SrTiO3(001) and show that the agreement between them is excellent. Significantly, the resulting Ti charge is considerably lower than the formal charge. The high degree of similarity between the Ti 2p spectrum for SrTiO3 and those for the rutile and anatase polymorphs of TiO2 suggests that the charge densities surrounding Ti in the latter materials are similar to that in SrTiO3. Taking a broader perspective, oxides containing other first-row transition metals also exhibit covalent character, leading to TM cation charges lower than the analogous fully ionic values in these materials as well.

Chambers, Scott A. (ORCID:000000025415043X)↗
NASA NTRS2004

A Kernel-Free Particle-Finite Element Method for Hypervelocity Impact Simulation

An improved hybrid particle-finite element method has been developed for the simulation of hypervelocity impact problems. Unlike alternative methods, the revised formulation computes the density without reference to any kernel or interpolation functions, for either the density or the rate of dilatation. This simplifies the state space model and leads to a significant reduction in computational cost. The improved method introduces internal energy variables as generalized coordinates in a new formulation of the thermomechanical Lagrange equations. Example problems show good agreement with exact solutions in one dimension and good agreement with experimental data in a three dimensional simulation.

Park, Young-Keun↗
DOE OSTI2022

Spin relaxation rate for heavy quarks in weakly coupled QCD plasma

We compute the relaxation rate of the spin density of heavy quarks in a perturbative QCD plasma to leading-log order in the coupling constant g. The spin relaxation rate Γ s in spin hydrodynamics is shown to be Γ s ~ g 4 log(1/g)T(T/M) 2 in the heavy-quark limit T/M $\ll$ 1, which is smaller than the relaxation rate of other non- hydrodynamic modes by additional powers of T/M. We demonstrate three different methods to evaluate the spin relaxation rate: 1) the Green-Kubo formula in the spin hydrodynamic regime, 2) the spin density correlation function in the strict hydrodynamic limit, and 3) quantum kinetic theory of the spin distribution function in momentum space. We highlight the interesting differences between these methods, while they are ultimately connected to each other by the underlying Ward-Takahashi identity for the non-conserved spin density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Mortality correlates with tree functional traits across a wood density gradient in the Central Amazon

Introduction: Understanding the mechanisms of tree mortality in tropical ecosystems remains challenging, in part due to the high diversity of tree species and the inherently stochastic nature of mortality. Plant functional traits offer a mechanistic link between plant physiology and performance, yet their ability to predict growth and mortality remains poorly understood. Given recent increases in tree mortality rates in the Amazon forest following extreme drought and wind events, we tested if lower wood density and acquisitive plant functional traits were associated with increased growth and mortality for common co-occurring trees in the Central Amazon. Methods: Seventeen trees of different species with similar sizes but a range in wood density (WD) and wood traits were felled, then assessed for 27 different individual functional parameters, including whole tree architecture, stem xylem anatomical and hydraulic traits and leaf traits. Traits of the individual trees were related to stand-level growth and mortality rates collected periodically over 30 years from nearby permanent inventory plots. Results: Higher wood density was associated with smaller leaf size, lower foliar base cations, lower stem water content and sapwood fraction, in agreement with the fast-slow plant economics spectrum. Lower wood density was associated with more acquisitive characteristics with greater hydraulic capacity and foliar nutrient concentrations, correlating with greater growth and mortality rates. Discussion: Our results show that lower wood density is part of a coordinated suite of traits linked to high resource acquisition, fast growth, and increased mortality risk, providing a functional framework for predicting species performance and forest vulnerability under future climate stress.

demographics↗
NASA NTRS1977

Analysis of photographic X-ray images

Some techniques used to extract quantitative data from the information contained in photographic images produced by grazing incidence soft X-ray optical systems are described. The discussion is focussed on the analysis of the data returned by the S-054 X-Ray Spectrographic Telescope Experiment on Skylab. The parameters of the instrument and the procedures used for its calibration are described. The technique used to convert photographic density to focal plane X-ray irradiance is outlined. The deconvolution of the telescope point response function from the image data is discussed. Methods of estimating the temperature, pressure, and number density of coronal plasmas are outlined.

Krieger, A. S.↗
DOE OSTI2019

Molecular Electrostatic Potential and Electron Density of Large Systems in Solution Computed with the Fragment Molecular Orbital Method

A solvent screening model for the molecular electrostatic potential is developed for the fragment molecular orbital method combined with the polarizable continuum model at the Hartree–Fock and density functional theory levels. The accuracy of the generated potentials is established in comparison to calculations without fragmentation. Solvent effects upon the molecular electrostatic potential and electron density of solute are discussed. Finally, the method is applied to two proteins: chignolin (PDB: 1UAO) and ovine prostaglandin H(2) synthase-1 (1EQG).

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

Deuterium Hugoniot: Pitfalls of thermodynamic sampling beyond density functional theory

Outstanding problems in the high-pressure phase diagram of hydrogen have demonstrated the need for more accurate ab initio methods for thermodynamic sampling. One promising method that has been deployed extensively above 100 GPa is coupled electron-ion Monte Carlo (CEIMC), which treats the electronic structure with quantum Monte Carlo (QMC). However, CEIMC predictions of the deuterium principal Hugoniot disagree significantly with experiment, overshooting the experimentally determined peak compression density by 7% and lower temperature gas-gun data by well over 20%. By deriving an equation relating the predicted Hugoniot density to underlying equation of state errors, we show that QMC and many-body methods can easily spoil the error cancellation properties inherent in the Rankine-Hugoniot relation, and very likely suffer from error addition. By cross validating QMC based on systematically improvable trial functions against post-Hartree-Fock many-body methods, we find that these methods introduce errors of the right sign and magnitude to account for much of the observed discrepancy between CEIMC and experiment. In conclusion, we stress that this is not just a CEIMC problem, but that thermodynamic sampling based on other many-body methods is likely to experience similar difficulties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2021

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗
DOE OSTI2020

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗