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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Correlation functions from a unified variational principle: Trial Lie groups

Time-dependent expectation values and correlation functions for many-body quantum systems are evaluated by means of a unified variational principle. It optimizes a generating functional depending on sources associated with the observables of interest. It is built by imposing through Lagrange multipliers constraints that account for the initial state (at equilibrium or off equilibrium) and for the backward Heisenberg evolution of the observables. The trial objects are respectively akin to a density operator and to an operator involving the observables of interest and the sources. We work out here the case where trial spaces constitute Lie groups. This choice reduces the original degrees of freedom to those of the underlying Lie algebra, consisting of simple observables; the resulting objects are labeled by the indices of a basis of this algebra. Explicit results are obtained by expanding in powers of the sources. Zeroth and first orders provide thermodynamic quantities and expectation values in the form of mean-field approximations, with dynamical equations having a classical Lie–Poisson structure. At second order, the variational expression for two-time correlation functions separates–as does its exact counterpart–the approximate dynamics of the observables from the approximate correlations in the initial state. Two building blocks are involved: (i) a commutation matrix which stems from the structure constants of the Lie algebra; and (ii) the second-derivative matrix of a free-energy function. The diagonalization of both matrices, required for practical calculations, is worked out, in a way analogous to the standard RPA. The ensuing structure of the variational formulae is the same as for a system of non-interacting bosons (or of harmonic oscillators) plus, at non-zero temperature, classical Gaussian variables. This property is explained by mapping the original Lie algebra onto a simpler Lie algebra. The results, valid for any trial Lie group, fulfill consistency properties and encompass several special cases: linear responses, static and time-dependent fluctuations, zero- and high-temperature limits, static and dynamic stability of small deviations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Linking density functional and mode coupling models for supercooled liquids

We compare predictions from two familiar models of the metastable supercooled liquid, respectively, constructed with thermodynamic and dynamic approaches. In the so called density functional theory the free energy F[ρ] of the liquid is a functional of the inhomogeneous density ρ(r). The metastable state is identified as a local minimum of F[ρ]. The sharp density profile characterizing ρ(r) is identified as a single particle oscillator, whose frequency is obtained from the parameters of the optimum density function. On the other hand, a dynamic approach to supercooled liquids is taken in the mode coupling theory (MCT) which predict a sharp ergodicity-non-ergodicity transition at a critical density. The single particle dynamics in the non-ergodic state, treated approximately, represents a propagating mode whose characteristic frequency is computed from the corresponding memory function of the MCT. The mass localization parameters in the above two models (treated in their simplest forms) are obtained, respectively, in terms of the corresponding natural frequencies depicted and are shown to have comparable magnitudes.

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

Constructing diabatic representations using adiabatic and approximate diabatic data – Coping with diabolical singularities

We have recently introduced a diabatization scheme, which simultaneously fits and diabatizes adiabatic ab initio electronic wave functions, Zhu and Yarkony J. Chem. Phys. 140, 024112 (2014). The algorithm uses derivative couplings in the defining equations for the diabatic Hamiltonian, H{sup d}, and fits all its matrix elements simultaneously to adiabatic state data. This procedure ultimately provides an accurate, quantifiably diabatic, representation of the adiabatic electronic structure data. However, optimizing the large number of nonlinear parameters in the basis functions and adjusting the number and kind of basis functions from which the fit is built, which provide the essential flexibility, has proved challenging. In this work, we introduce a procedure that combines adiabatic state and diabatic state data to efficiently optimize the nonlinear parameters and basis function expansion. Further, we consider using direct properties based diabatizations to initialize the fitting procedure. To address this issue, we introduce a systematic method for eliminating the debilitating (diabolical) singularities in the defining equations of properties based diabatizations. We exploit the observation that if approximate diabatic data are available, the commonly used approach of fitting each matrix element of H{sup d} individually provides a starting point (seed) from which convergence of the full H{sup d} construction algorithm is rapid. The optimization of nonlinear parameters and basis functions and the elimination of debilitating singularities are, respectively, illustrated using the 1,2,3,4{sup 1}A states of phenol and the 1,2{sup 1}A states of NH{sub 3}, states which are coupled by conical intersections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Bootstrapping the 3d Ising model at finite temperature

We estimate thermal one-point functions in the 3d Ising CFT using the operator product expansion (OPE) and the Kubo-Martin-Schwinger (KMS) condition. Several operator dimensions and OPE coefficients of the theory are known from the numerical bootstrap for flat-space four-point functions. Taking this data as input, we use a thermal Lorentzian inversion formula to compute thermal one-point coefficients of the first few Regge trajectories in terms of a small number of unknown parameters. We approximately determine the unknown parameters by imposing the KMS condition on the two-point functions (σσ) and (ϵϵ). As a result, we estimate the one-point functions of the lowest-dimension $\mathbb{Z}_2$-even scalar ϵ and the stress energy tensor T μν . Our result for (σσ) at finite-temperature agrees with Monte-Carlo simulations within a few percent, inside the radius of convergence of the OPE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

High temperature neutron diffraction investigation of U 3 Si 2

The U-Si system is actively undergoing studies due to its promise as a component of an accident tolerant nuclear fuel. In this work, the crystal structure of the U 3 Si 2 compound in the U-Si system was investigated as a function of temperature from room temperature to 1298 K for a sample of approximately 42.02 at.% Si composition and 1373 K for a sample with approximately 42.10 at.% Si using high temperature time-of-flight neutron diffraction on the High-Pressure Preferred Orientation (HIPPO) diffractometer at Los Alamos Neutron Science Center (LANSCE). The simultaneous Rietveld refinement of five histograms from the five HIPPO detector backs (40°, 60°, 90°, 120°, and 145°) provided datasets for the lattice parameters, anisotropic atomic displacement parameters, and atomic positions as a function of temperature. To explore the possibility of a homogeneity range, two sample compositions were analyzed, a 40.02 at.% Si (near stoichiometric) and a potentially hyperstoichiometric (40.10 at.% Si) sample. While minor differences in the anisotropic atomic displacement parameters between the two samples were observed, over the entire investigated temperature range no additional phases were found indicative of a single phase sample and the ability of the U 3 Si 2 compound to accommodate nonstoichiometry, suggesting that the U 3 Si 2 compound is not a line compound confirming previous determinations. However, the differences observed in the average coefficient of linear thermal expansion identified between the two compositions warrant future investigation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Applications of the modified Hulthén-Kohn method for bound and scattering states

We adapt the Hulthén–Kohn method suggested by Efros [Phys. Rev. C 99, 034620 (2019)] for calculating various observables in the continuum and discrete spectrum using two-body interactions in single- and coupled channel systems. We explore the convergence of phase shifts and wave functions as well as the location of S-matrix poles which enables obtaining both resonance and bound state parameters. We find that employing a harmonic oscillator basis, together with an interaction smoothing scheme introduced by Gyarmati et al. [Nucl. Phys. A 326, 119 (1979)], and adopting approximate bound-state solutions for the short-range components of basis wave functions lead to good convergence even with restricted oscillator quanta accessible for modern no-core shell model codes. The adapted Efros method will facilitate ab initio many-body nuclear structure applications.

Nuclear reactions↗
DOE OSTI2016

Microscopic origin of the charge transfer in single crystals based on thiophene derivatives: A combined NEXAFS and density functional theory approach

We have investigated the charge transfer mechanism in single crystals of DTBDT-TCNQ and DTBDT-F{sub 4}TCNQ (where DTBDT is dithieno[2,3-d;2′,3′-d′] benzo[1,2-b;4,5-b′]dithiophene) using a combination of near-edge X-ray absorption spectroscopy (NEXAFS) and density functional theory calculations (DFT) including final state effects beyond the sudden state approximation. In particular, we find that a description that considers the partial screening of the electron-hole Coulomb correlation on a static level as well as the rearrangement of electronic density shows excellent agreement with experiment and allows to uncover the details of the charge transfer mechanism in DTBDT-TCNQ and DTBDT-F{sub 4} TCNQ, as well as a reinterpretation of previous NEXAFS data on pure TCNQ. Finally, we further show that almost the same quality of agreement between theoretical results and experiment is obtained by the much faster Z+1/2 approximation, where the core hole effects are simulated by replacing N or F with atomic number Z with the neighboring atom with atomic number Z+1/2.

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

Principal shapes and squeezed limits in the effective field theory of large scale structure

We apply an orthogonalization procedure on the effective field theory of large scale structure (EFT of LSS) shapes, relevant for the angle-averaged bispectrum and non-Gaussian covariance of the matter power spectrum at one loop. Assuming natural-sized EFT parameters, this identifies a linear combination of EFT shapes—referred to as the principal shape—that gives the dominant contribution for the whole kinematic plane, with subdominant combinations suppressed by a few orders of magnitude. For the covariance, our orthogonal transformation is in excellent agreement with a principal component analysis applied to available data. Additionally we find that, for both observables, the coefficients of the principal shapes are well approximated by the EFT coefficients appearing in the squeezed limit, and are thus measurable from power spectrum response functions. Employing data from N-body simulations for the growth-only response, we measure the single EFT coefficient describing the angle-averaged bispectrum with Ο (10%) precision. These methods of shape orthogonalization and measurement of coefficients from response functions are valuable tools for developing the EFT of LSS framework, and can be applied to more general observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Green’s function methods for excited states and x-ray spectra of functional materials

Many interesting properties of functional materials, such as dynamic response and thermodynamic behavior, depend on their excited state properties. These functional properties are often related to excitations in the system, such as phonons and plasmons, which lead to inelastic losses, lifetime, and other dynamic effects. The excitations are pure many-body correlation effects that are missing from independent particle theories. They are revealed in x-ray spectra such as photoemission and absorption, where they show up as satellites beyond the quasi-particle approximation. Our main focus in this work is the use of Green's function methods to describe these effects. In particular, we discuss how the cumulant Green's function provides a unified treatment of such dynamic correlation effects in many contexts. Besides a robust theoretical framework, these methods also yield widely applicable tools for practical calculations of many functional properties of materials. Furthermore, this methodology is illustrated with a number of applications ranging from optical and x-ray spectra to thermodynamic properties, and dynamic response. Some recent extensions for more correlated systems are also briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

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

The Acceleration of Charged Particles at a Spherical Shock Moving through an Irregular Magnetic Field

We investigate the physics of charged-particle acceleration at spherical shocks moving into a uniform plasma containing a turbulent magnetic field with a uniform mean. This has applications to particle acceleration at astrophysical shocks, most notably, to supernovae blast waves. We numerically integrate the equations of motion of a large number of test protons moving under the influence of electric and magnetic fields determined from a kinematically defined plasma flow associated with a radially propagating blast wave. Distribution functions are determined from the positions and velocities of the protons. The unshocked plasma contains a magnetic field with a uniform mean and an irregular component having a Kolmogorov-like power spectrum. The field inside the blast wave is determined from Maxwell’s equations. The angle between the average magnetic field and unit normal to the shock varies with position along its surface. It is quasi-perpendicular to the unit normal near the sphere’s equator, and quasi-parallel to it near the poles. We find that the highest intensities of particles, accelerated by the shock, are at the poles of the blast wave. The particles “collect” at the poles as they approximately adhere to magnetic field lines that move poleward from their initial encounter with the shock at the equator, as the shock expands. The field lines at the poles have been connected to the shock the longest. We also find that the highest-energy protons are initially accelerated near the equator or near the quasi-perpendicular portion of the shock, where the acceleration is more rapid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-J-238: First-Order Approximation of Time-Resolved 4DMRI From Cine 2DMRI and Respiratory-Correlated 4DMRI

Purpose: Cine 2DMRI is useful in MR-guided radiotherapy but it lacks volumetric information. We explore the feasibility of estimating timeresolved (TR) 4DMRI based on cine 2DMRI and respiratory-correlated (RC) 4DMRI though simulation. Methods: We hypothesize that a volumetric image during free breathing can be approximated by interpolation among 3DMRI image sets generated from a RC-4DMRI. Two patients’ RC-4DMRI with 4 or 5 phases were used to generate additional 3DMRI by interpolation. For each patient, six libraries were created to have total 5-to-35 3DMRI images by 0–6 equi-spaced tri-linear interpolation between adjacent and full-inhalation/full-exhalation phases. Sagittal cine 2DMRI were generated from reference 3DMRIs created from separate, unique interpolations from the original RC-4DMRI. To test if accurate 3DMRI could be generated through rigid registration of the cine 2DMRI to the 3DMRI libraries, each sagittal 2DMRI was registered to sagittal cuts in the same location in the 3DMRI within each library to identify the two best matches: one with greater lung volume and one with smaller. A final interpolation between the corresponding 3DMRI was then performed to produce the first-order-approximation (FOA) 3DMRI. The quality and performance of the FOA as a function of library size was assessed using both the difference in lung volume and average voxel intensity between the FOA and the reference 3DMRI. Results: The discrepancy between the FOA and reference 3DMRI decreases as the library size increases. The 3D lung volume difference decreases from 5–15% to 1–2% as the library size increases from 5 to 35 image sets. The average difference in lung voxel intensity decreases from 7–8 to 5–6 with the lung intensity being 0–135. Conclusion: This study indicates that the quality of FOA 3DMRI improves with increasing 3DMRI library size. On-going investigations will test this approach using actual cine 2DMRI and introduce a higher order approximation for improvements. This study is in part supported by NIH (U54CA137788 and U54CA132378)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Continuation of probability density functions using a generalized Lyapunov approach

Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the efficient solution of a generalized Lyapunov equation using an iterative method involving low-rank approximations. We apply and illustrate the capabilities of the method using a problem in physical oceanography, i.e. the occurrence of multiple steady states of the Atlantic Ocean circulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

What is a parton shower?

We consider idealized parton shower event generators that treat parton spin and color exactly, leaving aside the choice of practical approximations for spin and color. We investigate how the structure of such a parton shower generator is related to the structure of QCD. We argue that a parton shower with splitting functions proportional to ${{\alpha}}_{\mathrm{s}}$ can be viewed not just as a model, but as the lowest order approximation to a shower that is defined at any perturbative order. To support this argument, we present a formulation for a parton shower at order ${{\alpha}}_{\mathrm{s}}^{k}$ for any $k$. Since some of the input functions needed are specified by their properties but not calculated, this formulation does not provide a useful recipe for an order ${{\alpha}}_{\mathrm{s}}^{k}$ parton shower algorithm. However, in this formulation we see how the operators that generate the shower are related to operators that specify the infrared singularities of QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Evolution of Chain Dynamics and Oxidation States with Increasing Chain Length for a Donor–Acceptor-Conjugated Oligomer Series

While it is known that the chain length strongly affects the properties of π-conjugated polymers, the effects of chain length on the molecular structure, chain conformation, and oxidation state properties in donor–acceptor-type conjugated structures remain unclear. This limits our understanding of how the polymer molecular weight impacts material properties. Here, a discrete and monodisperse oligomer series (n = 3, 5, 7, 9, 15, and 21) and polymers (nPB), composed of the donor 3,4-propylenedioxythiophene (ProDOT), the acceptor benzothiadiazole (BTD), and methylthio end-capping groups, are synthesized by C–H-activated cross-coupling. The molecular structure, molecular weight, and dispersity of each oligomer/polymer are thoroughly characterized by nuclear magnetic resonance spectroscopy, mass spectrometry, and gel permeation chromatography. This series reveals a rod-to-coil transition at n = 15 and coil formation at polymer length scales of ~28 units via solution small-angle neutron scattering characterization. The oxidation states are deciphered via cyclic voltammetry, differential pulse voltammetry, spectroelectrochemistry, and density functional theory calculations. Oligomers 3–9 undergo successive one-electron oxidation steps, while 15 and higher undergo multielectron oxidations per step in CH 2 Cl 2 –TBAPF6 at a Pt ultramicroelectrode. The electronic transition of each oxidation state (1+, 2+, 3+, etc.) is tracked by absorption spectroscopy, revealing a “bipolaron to di-polaron” transition at n = 7 at which shorter oligomers prefer bipolaron formation and longer oligomers prefer di-polaron formation in their dication states. Furthermore, oxidized 15 has a lower electronic transition energy compared to its polymer homologue, as shown by spectroelectrochemistry, revealing a synergy between the chain length and the oxidation state properties. This study proves that the convergence limit between small molecule and polymer behavior occurs at approximately 15 units and highlights the property transitions that occur as a function of chain length for a donor–acceptor class of conjugated organic materials.

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

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Proposal for Universality in the Viscosity of Metallic Liquids

The range of magnitude of the liquid viscosity, η, as a function of temperature is one of the most impressive of any physical property, changing by approximately 17 orders of magnitude from its extrapolated value at infinite temperature (η o ) to that at the glass transition temperature, T g . We present experimental measurements of containerlessly processed metallic liquids that suggest that log(η/η o ) as a function of T A /T is a potentially universal scaled curve. In stark contrast to previous approaches, the scaling requires only two fitting parameters, which are on average predictable. The temperature T A corresponds to the onset of cooperative motion and is strongly correlated with T g , suggesting that the processes underlying the glass transition first appear in the high temperature liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Assessment of TD-DFT and LF-DFT for study of d − d transitions in first row transition metal hexaaqua complexes

Herein, we present the systematic, comparative computational study of the d − d transitions in a series of first row transition metal hexaaqua complexes, [M(H{sub 2}O){sub 6}]{sup n+} (M{sup 2+/3+} = V {sup 2+/3+}, Cr{sup 2+/3+}, Mn{sup 2+/3+}, Fe{sup 2+/3+}, Co{sup 2+/3+}, Ni{sup 2+}) by the means of Time-dependent Density Functional Theory (TD-DFT) and Ligand Field Density Functional Theory (LF-DFT). Influence of various exchange-correlation (XC) approximations have been studied, and results have been compared to the experimental transition energies, as well as, to the previous high-level ab initio calculations. TD-DFT gives satisfactory results in the cases of d{sup 2}, d{sup 4}, and low-spin d{sup 6} complexes, but fails in the cases when transitions depend only on the ligand field splitting, and for states with strong character of double excitation. LF-DFT, as a non-empirical approach to the ligand field theory, takes into account in a balanced way both dynamic and non-dynamic correlation effects and hence accurately describes the multiplets of transition metal complexes, even in difficult cases such as sextet-quartet splitting in d{sup 5} complexes. Use of the XC functionals designed for the accurate description of the spin-state splitting, e.g., OPBE, OPBE0, or SSB-D, is found to be crucial for proper prediction of the spin-forbidden excitations by LF-DFT. It is shown that LF-DFT is a valuable alternative to both TD-DFT and ab initio methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗