Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Methane dissociative chemisorption and detailed balance on Pt(111): Dynamical constraints and the modest influence of tunneling

A dynamically biased (d-) precursor mediated microcanonical trapping (PMMT) model of the activated dissociative chemisorption of methane on Pt(111) is applied to a wide range of dissociative sticking experiments, and, by detailed balance, to the methane product state distributions from the thermal associative desorption of adsorbed hydrogen with coadsorbed methyl radicals. Tunneling pathways were incorporated into the d-PMMT model to better replicate the translational energy distribution of the desorbing methane product from the laser induced thermal reaction of coadsorbed hydrogen and methyl radicals occurring near T{sub s} = 395 K. Although tunneling is predicted to be inconsequential to the thermal dissociative chemisorption of CH{sub 4} on Pt(111) at the high temperatures of catalytic interest, once the temperature drops to 395 K the tunneling fraction of the reactive thermal flux reaches 15%, and as temperatures drop below 275 K the tunneling fraction exceeds 50%. The d-PMMT model parameters of (E{sub 0} = 58.9 kJ/mol, s = 2, η{sub v} = 0.40) describe the apparent threshold energy for CH{sub 4}/Pt(111) dissociative chemisorption, the number of surface oscillators involved in the precursor complex, and the efficacy of molecular vibrational energy to promote reaction, relative to translational energy directed along the surface normal. Molecular translations parallel to the surface and rotations are treated as spectator degrees of freedom. Transition state vibrational frequencies are derived from generalized gradient approximation-density functional theory electronic structure calculations. The d-PMMT model replicates the diverse range of experimental data available with good fidelity, including some new effusive molecular beam and ambient gas dissociative sticking measurements. Nevertheless, there are some indications that closer agreement between theory and experiments could be achieved if a surface efficacy less than one was introduced into the modeling as an additional dynamical constraint.

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

Towards quantifying the role of exact exchange in predictions of transition metal complex properties

We estimate the prediction sensitivity with respect to Hartree-Fock exchange in approximate density functionals for representative Fe(II) and Fe(III) octahedral complexes. Based on the observation that the range of parameters spanned by the most widely employed functionals is relatively narrow, we compute electronic structure property and spin-state orderings across a relatively broad range of Hartree-Fock exchange (0%-50%) ratios. For the entire range considered, we consistently observe linear relationships between spin-state ordering that differ only based on the element of the direct ligand and thus may be broadly employed as measures of functional sensitivity in predictions of organometallic compounds. The role Hartree-Fock exchange in hybrid functionals is often assumed to play is to correct self-interaction error-driven electron delocalization (e.g., from transition metal centers to neighboring ligands). Surprisingly, we instead observe that increasing Hartree-Fock exchange reduces charge on iron centers, corresponding to effective delocalization of charge to ligands, thus challenging notions of the role of Hartree-Fock exchange in shifting predictions of spin-state ordering.

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

The electronic structure, mechanical and thermodynamic properties of Mo{sub 2}XB{sub 2} and MoX{sub 2}B{sub 4} (X = Fe, Co, Ni) ternary borides

The mechanical properties, electronic structure and thermodynamic properties of the Mo{sub 2}XB{sub 2} and MoX{sub 2}B{sub 4} (X = Fe, Co, Ni) ternary borides were calculated by first-principles methods. The elastic constants show that these ternary borides are mechanically stable. Formation enthalpy of Mo{sub 2}XB{sub 2} and MoX{sub 2}B{sub 4} (X = Fe, Co, Ni) ternary borides are at the range of −118.09 kJ/mol to −40.14 kJ/mol. The electronic structures and chemical bonding characteristics are analyzed by the density of states. Mo{sub 2}FeB{sub 2} has the largest shear and Young's modulus because of its strong chemical bonding, and the values are 204.3 GPa and 500.3 GPa, respectively. MoCo{sub 2}B{sub 4} shows the lowest degree of anisotropy due to the lack of strong direction in the bonding. The Debye temperature of MoFe{sub 2}B{sub 4} is the largest among the six phases, which means that MoFe{sub 2}B{sub 4} possesses the best thermal conductivity. Enthalpy shows an approximately linear function of the temperature above 300 K. The entropy of these compounds increase rapidly when the temperature is below 450 K. The Gibbs free energy decreases with the increase in temperature. MoCo{sub 2}B{sub 4} has the lowest Gibbs free energy, which indicates the strongest formation ability in Mo{sub 2}XB{sub 2} and MoX{sub 2}B{sub 4} (X = Fe, Co, Ni) ternary borides.

36 MATERIALS SCIENCE↗
DOE OSTI2016

On the accuracy of the Padé-resummed master equation approach to dissipative quantum dynamics

Well-defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé-resummed master equations with numerically exact results in the context of the spin–boson model demonstrates that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques.

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

Communication: Fitting potential energy surfaces with fundamental invariant neural network

A more flexible neural network (NN) method using the fundamental invariants (FIs) as the input vector is proposed in the construction of potential energy surfaces for molecular systems involving identical atoms. Mathematically, FIs finitely generate the permutation invariant polynomial (PIP) ring. In combination with NN, fundamental invariant neural network (FI-NN) can approximate any function to arbitrary accuracy. Because FI-NN minimizes the size of input permutation invariant polynomials, it can efficiently reduce the evaluation time of potential energy, in particular for polyatomic systems. In this work, we provide the FIs for all possible molecular systems up to five atoms. Potential energy surfaces for OH{sub 3} and CH{sub 4} were constructed with FI-NN, with the accuracy confirmed by full-dimensional quantum dynamic scattering and bound state calculations.

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

A near X-point charge exchange neutral spectroscopy (CENS) system for DIII-D

A 16 channel spectroscopy system has been installed on DIII-D to provide information about the energy distribution of the atomic neutrals using the Doppler shift and broadening of passive Balmer-α emission. Here, the lines of sight are dominantly in the toroidal direction, with tangencies near the lowest point on closed magnetic flux surfaces moving from the lower divertor upward into the confined plasma. This allows the system to provide details of the neutrals as they undergo various atomic physics processes while traveling upward from the lower divertor. The spectrally resolved measurements provide several advantages that complement typical optical filter based measurements of hydrogenic spectral lines. These include more direct measurements of the neutrals that fuel deeper in the plasma and capturing neutral velocity distribution information via the Doppler broadened and shifted line. In addition, the spectral separation of higher energy “thermal” neutrals and bright cold emission from the scrape-off layer allows for a more straightforward calculation of the underlying neutral densities based on the emission because of the dependence of the thermal emission on confined plasma properties that are approximately flux functions and well measured.

Atomic and molecular physics↗
DOE OSTI2026

Mott vs Kondo: Influence of various density functional based methods on the Ce isostructural phase transition mechanism

The cerium iso-structural phase transition (γ–α) is dominated by f-electron localization changes that results in a magnetic ordering change and a volume collapse. Generally, these physics are difficult to capture with ab initio and first principles methods. However, previous works have shown various methods to be successful in predicting at least some of the physics of the γ to α phase transition. Therefore, here, we perform a broad survey of density functional based methods across three levels of theory and types of functions [GGA (generalized gradient approximation), MetaGGA, and hybrid functionals] and compare the results, focusing on hydrostatic compression across the phase boundary at zero kelvin. For the methods that best reproduce experimental results, we directly probe the predicted mechanisms and frame the results in the Mott/Kondo debate, assessing how the underlying methods and assumptions of different functionals can assess the physical drivers in the phase transition, providing insight into the governing processes of this unique phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A near-infrared census of the multicomponent stellar structure of early-type dwarf galaxies in the Virgo cluster

The fraction of star-forming to quiescent dwarf galaxies varies from almost infinity in the field to zero in the centers of rich galaxy clusters. What is causing this pronounced morphology-density relation? What do quiescent dwarf galaxies look like when studied in detail, and what conclusions can be drawn about their formation mechanism? Here we study a nearly magnitude-complete sample (–19 < M{sub r} < –16 mag) of 121 Virgo cluster early types with deep near-infrared images from the SMAKCED project. We fit two-dimensional models with optional inner and outer components, as well as bar and lens components (in ∼15% of the galaxies), to the galaxy images. While a single Sérsic function may approximate the overall galaxy structure, it does not entirely capture the light distribution of two-thirds of our galaxies, for which multicomponent models provide a better fit. This fraction of complex galaxies shows a strong dependence on luminosity, being larger for brighter objects. We analyze the global and component-specific photometric scaling relations of early-type dwarf galaxies and discuss similarities with bright early and late types. The dwarfs' global galaxy parameters show scaling relations that are similar to those of bright disk galaxies. The inner components are mostly fitted with Sérsic n values close to 1. At a given magnitude, they are systematically larger than the bulges of spirals, suggesting that they are not ordinary bulges. We argue that the multicomponent structures in early-type dwarfs are mostly a phenomenon inherent to the disks and may indeed stem from environmental processing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Linear response to long wavelength fluctuations using curvature simulations

We study the local response to long wavelength fluctuations in cosmological N -body simulations, focusing on the matter and halo power spectra, halo abundance and non-linear transformations of the density field. The long wavelength mode is implemented using an effective curved cosmology and a mapping of time and distances. The method provides an alternative, more direct, way to measure the isotropic halo biases. Limiting ourselves to the linear case, we find generally good agreement between the biases obtained from the curvature method and the traditional power spectrum method at the level of a few percent. We also study the response of halo counts to changes in the variance of the field and find that the slope of the relation between the responses to density and variance differs from the naïve derivation assuming a universal mass function by approximately 8–20%. This has implications for measurements of the amplitude of local non-Gaussianity using scale dependent bias. We also analyze the halo power spectrum and halo-dark matter cross-spectrum response to long wavelength fluctuations and derive second order halo bias from it, as well as the super-sample variance contribution to the galaxy power spectrum covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spatiotemporal forecasting of the edge localized modes in tokamak plasmas using neural networks

Artificial intelligence techniques have been increasingly adopted by the plasma and fusion science to address problems like plasma reconstruction, surrogate modeling, and tokamak/stellarator optimization. A key focus in sustained fusion research is the prediction and mitigation of edge-localized-modes (ELMs), instabilities that occur in short, periodic bursts and can cause erosion to the tokamak vessel wall. Recent research has demonstrated the power of neural networks in approximating continuous functions. In this work, we build spatiotemporal forecasting models that can predict the onset of ELMs and their evolution at early stages. We leverage recent advances in generative modeling, sequence-to-sequence modeling, and Fourier neural operators to propose architectures and training strategies that can learn to forecast short to long term dynamics of the noisy signals due to ELMs. We benchmark the developed model against a state-of-the-art foundation model using the beam emission spectroscopy (BES) data that captures the plasma fluctuations due to ELMs over a 8 x 8 spatial grid. Our models demonstrate high accuracy, outperforming the baselines, in predicting the evolution of BES signals during ELM events. Furthermore, the developed models exhibit high accuracy in predicting the rapid rise and relaxation of the signals due to ELMs within 30–80 µs.

edge localized modes↗
DOE OSTI2020

Two-particle correlation functions in cluster perturbation theory: Hubbard spin susceptibilities

Cluster perturbation theory (CPT) is a computationally economic method commonly used to estimate the momentum- and energy-resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. Finally, we benchmark the method with the one-dimensional Fermi-Hubbard model by comparing with known results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of radiation-induced defects on the superfluid density and optical conductivity of overdoped La 2 - x Sr x CuO 4

Using a combination of time-domain THz spectroscopy (TDTS) and mutual inductance measurements, we have investigated the low-energy electrodynamic response of overdoped La$_{2-x}$Sr$_x$CuO$_4$ films that have been exposed to ion irradiation. Films went through three rounds of irradiation (2, 4, and 6 $\times 10^{13}$ ions/cm$^2$) and mututal inductance and TDTS experiments were performed between each step. Together with the as-grown film, this gives four different levels of disorder. The transport scattering rate that is measured directly in the THz experiments is an approximately linear function of the radiation dose at all temperatures. This is consistent with a proportionate increase in elastic scattering from point-like 2D disorder. Here, in the superconducting state we find that the relation between T$_c$, the superfluid density, and the scattering rates are quantitatively at odds with the predictions based on the extant theory of Abrikosov-Gorkov-like pair breaking in a dirty d-wave superconductor. Increasing disorder causes only a small change in T$_c$ for the overdoped films, but the changes to the $\omega \sim 0$ superfluid density are much larger.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Frequency-dependent dielectric function of semiconductors with application to physisorption

Furthermore, the dielectric function is one of the most important quantities that describes the electrical and optical properties of solids. Accurate modeling of the frequency-dependent dielectric function has great significance in the study of the long-range van der Waals (vdW) interaction for solids and adsorption. In this work we calculate the frequency-dependent dielectric functions of semiconductors and insulators using the GW method with and without exciton effects, as well as efficient semilocal density functional theory (DFT), and compare these calculations with a model frequency-dependent dielectric function. We find that for semiconductors with moderate band gaps, the model dielectric functions, GW values, and DFT calculations all agree well with each other. However, for insulators with strong exciton effects, the model dielectric functions have a better agreement with accurate GW values than the DFT calculations, particularly in high-frequency region. To understand this, we repeat the DFT calculations with scissors correction, by shifting the DFT Kohn-Sham energy levels to match the experimental band gap. We find that scissors correction only moderately improves the DFT dielectric function in the low-frequency region. Based on the dielectric functions calculated with different methods, we make a comparative study by applying these dielectric functions to calculate the vdW coefficients (C3 and C5) for adsorption of rare-gas atoms on a variety of surfaces. We find that the vdW coefficients obtained with the nearly free electron gas-based model dielectric function agree quite well with those obtained from the GW dielectric function, in particular for adsorption on semiconductors, leading to an overall error of less than 7% for C 3 and 5% for C 5 . This demonstrates the reliability of the model dielectric function for the study of physisorption.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Resummation of jet veto logarithms at N 3 LL a + NNLO for W + W ? production at the LHC

We compute the resummed on-shell W+W- production cross section under a jet veto at the LHC to partial N3LL order matched to the fixed-order NNLO result. Differential NNLO cross sections are obtained from an implementation of qT subtraction in Sherpa. The two-loop virtual corrections to the qq¯→W+W- amplitude, used in both fixed-order and resummation predictions, are extracted from the public code qqvvamp. We perform resummation using soft collinear effective theory, with approximate beam functions where only the logarithmic terms are included at two-loop. In addition to scale uncertainties from the hard matching scale and the factorization scale, rapidity scale variations are obtained within the analytic regulator approach. Our resummation results show a decrease in the jet veto cross section compared to NNLO fixed-order predictions, with reduced scale uncertainties compared to NNLL+NLO resummed predictions. We include the loop-induced gg contribution with jet veto resummation to NLL+LO. The prediction shows good agreement with recent LHC measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment

36 MATERIALS SCIENCE↗