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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Experimental determination of the valence band offsets of ZnGeN 2 and (ZnGe) 0.94 Ga 0.12 N 2 with GaN

A predicted type-II staggered band alignment with an approximately 1.4 eV valence band offset at the ZnGeN 2 /GaN heterointerface has inspired novel band-engineered III-N/ZnGeN 2 heterostructure-based device designs for applications in high performance optoelectronics. We report on the determination of the valence band offset between metalorganic chemical vapor deposition grown (ZnGe) 1- x Ga 2 x N 2 , for x = 0 and 0.06, and GaN using x-ray photoemission spectroscopy. The valence band of ZnGeN 2 was found to lie 1.45–1.65 eV above that of GaN. This result agrees well with the value predicted by first-principles density functional theory calculations using the local density approximation for the potential profile and quasiparticle self-consistent GW calculations of the band edge states relative to the potential. For (ZnGe) 0.94 Ga 0.12 N 2 the value was determined to be 1.29 eV, ~10%–20% lower than that of ZnGeN 2 . The experimental determination of the large band offset between ZnGeN 2 and GaN provides promising alternative solutions to address challenges faced with pure III-nitride-based structures and devices.

Physics↗

DOE OSTI2017

Free energy functionals for polarization fluctuations: Pekar factor revisited

The separation of slow nuclear and fast electronic polarization in problems related to electron mobility in polarizable media was considered by Pekar 70 years ago. Within dielectric continuum models, this separation leads to the Pekar factor in the free energy of solvation by the nuclear degrees of freedom. The main qualitative prediction of Pekar’s perspective is a significant, by about a factor of two, drop of the nuclear solvation free energy compared to the total (electronic plus nuclear) free energy of solvation. The Pekar factor enters the solvent reorganization energy of electron transfer reactions and is a significant mechanistic parameter accounting for the solvent effect on electron transfer. Here, we study the separation of the fast and slow polarization modes in polar molecular liquids (polarizable dipolar liquids and polarizable water force fields) without relying on the continuum approximation. We derive the nonlocal free energy functional and use atomistic numerical simulations to obtain nonlocal, reciprocal space electronic and nuclear susceptibilities. A consistent transition to the continuum limit is introduced by extrapolating the results of finite-size numerical simulation to zero wavevector. The continuum nuclear susceptibility extracted from simulations is numerically close to the Pekar factor. However, we derive a new functionality involving the static and high-frequency dielectric constants. The main distinction of our approach from the traditional theories is found for the solvation free energy due to the nuclear polarization: the anticipated significant drop of its magnitude with increasing liquid polarizability does not occur. The reorganization energy of electron transfer is either nearly constant with increasing the solvent polarizability and the corresponding high-frequency dielectric constant (polarizable dipolar liquids) or actually noticeably increases (polarizable force fields of water).

14 SOLAR ENERGY↗

DOE OSTI2022

Analytic gradients for compressed multistate pair-density functional theory

Photochemical reactions often involve states that are closely coupled due to near degeneracies, for example by proximity to conical intersections. Therefore, a multistate method is used to accurately describe these states; for example, ordinary perturbation theory is replaced by quasidegenerate perturbation theory. Multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to approximate the full dynamical correlation energy of strongly correlated systems, and we recently proposed compressed multistate pair-density functional theory (CMS-PDFT) to treat closely coupled states. In the present paper, we report the implementation of analytic gradients for CMS-PDFT in both OpenMolcas and PySCF, and here we illustrate the use of these gradients by applying the method to the excited states of formaldehyde and phenol.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2014

MO-E-17A-10: Evaluation of Body and Head Dimensions of Pediatric Patients as a Function of Age

Purpose: Phantom development in medical physics plays an important role in radiation dose assessment and image quality evaluation, and this is especially true in the pediatric patient population. The purpose of this investigation was to establish the relationship between patient age and patient size, focusing on the abdomen-pelvis and head effective diameters, for patients ranging in age from newborn to 18 years. Methods: A dose reporting tool for computed tomography systems was installed at our institution to achieve compliance with state law commencing on July 1, 2012. The software records a number of patient-specific data, and also reports CT dose metrics (CTDIvol and DLP) into the patients interpretive report. The database generated by the software was mined to determine patient effective diameter as a function of age for pediatric patients aged 0–18 years. CT protocols including abdomen-pelvis and routine head were evaluated, and specific to this study the patients age, gender and equivalent diameter were recorded. Results: Six age ranges were evaluated: A(0–3), B(4–6), C(7–9), D(10–12),E(13–15),F(16–18). For the torso in these groups based upon 694 patients, median effective diameters were 147, 167, 184, 214, 231, 246 mm, respectively. For the head (N=1833), median diameters were 143, 157, 162, 168, 174, and 174, respectively. Conclusion: A solid understanding of the approximate dimensions of pediatric patients as a function of age is useful in the development of age-based imaging protocols and dose assessments. CT dose-reporting tools generate a great deal of data with respect to body dimensions automatically. In this study, median equivalent diameters for the abdomen-pelvis and head of pediatric patients were evaluated. These data may prove useful in the development of both mathematical and physical phantoms for dosimetry and image quality assessment.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TU-H-CAMPUS-IeP1-01: Bias and Computational Efficiency of Variance Reduction Methods for the Monte Carlo Simulation of Imaging Detectors

Purpose: Variance reduction techniques (VRTs) are employed in Monte Carlo simulations to obtain estimates with reduced statistical uncertainty for a given simulation time. In this work, we study the bias and efficiency of a VRT for estimating the response of imaging detectors. Methods: We implemented Directed Sampling (DS), preferentially directing a fraction of emitted optical photons directly towards the detector by altering the isotropic model. The weight of each optical photon is appropriately modified to maintain simulation estimates unbiased. We use a Monte Carlo tool called fastDETECT2 (part of the hybridMANTIS open-source package) for optical transport, modified for VRT. The weight of each photon is calculated as the ratio of original probability (no VRT) and the new probability for a particular direction. For our analysis of bias and efficiency, we use pulse height spectra, point response functions, and Swank factors. We obtain results for a variety of cases including analog (no VRT, isotropic distribution), and DS with 0.2 and 0.8 optical photons directed towards the sensor plane. We used 10,000, 25-keV primaries. Results: The Swank factor for all cases in our simplified model converged fast (within the first 100 primaries) to a stable value of 0.9. The root mean square error per pixel for DS VRT for the point response function between analog and VRT cases was approximately 5e-4. Conclusion: Our preliminary results suggest that DS VRT does not affect the estimate of the mean for the Swank factor. Our findings indicate that it may be possible to design VRTs for imaging detector simulations to increase computational efficiency without introducing bias.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Cluster decomposition of full configuration interaction wave functions: A tool for chemical interpretation of systems with strong correlation

Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size, thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e., those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI algorithm. Furthermore, the cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the π space of polyacenes, and (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with an increasing rank for the first three systems, even connected octuple excitations still appear important in Cr 2 , suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Non-adiabatic ab initio molecular dynamics of supersonic beam epitaxy of silicon carbide at room temperature

In this work, we investigate the processes leading to the room-temperature growth of silicon carbide thin films by supersonic molecular beam epitaxy technique. We present experimental data showing that the collision of fullerene on a silicon surface induces strong chemical-physical perturbations and, for sufficient velocity, disruption of molecular bonds, and cage breaking with formation of nanostructures with different stoichiometric character. We show that in these out-of-equilibrium conditions, it is necessary to go beyond the standard implementations of density functional theory, as ab initio methods based on the Born-Oppenheimer approximation fail to capture the excited-state dynamics. In particular, we analyse the Si-C{sub 60} collision within the non-adiabatic nuclear dynamics framework, where stochastic hops occur between adiabatic surfaces calculated with time-dependent density functional theory. This theoretical description of the C{sub 60} impact on the Si surface is in good agreement with our experimental findings.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Implications of the fractional charge of hydroxide at the electrochemical interface

Rational design of materials that efficiently convert electrical energy into chemical bonds will ultimately depend on a thorough understanding of the electrochemical interface at the atomic level. Towards this goal, the use of density functional theory (DFT) at the generalized gradient approximation (GGA) level has been applied widely in the past 15 years. In the calculation of electrochemical reaction energetics using GGA-DFT, it is frequently implicitly assumed that ions in the Helmholtz plane have unit charge. However, the ion charge is observed to be fractional near the interface through both a capacitor model and through Bader charge partitioning. In this work, we show that this spurious charge transfer can be effectively mitigated by continuum charging of the electrolyte. Finally, we show that, similar to hydronium, the observed fractional charge of hydroxide is not due to a GGA level self-interaction error, as the partial charge is observed even when using hybrid level exchange–correlation functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

A Model for Approximately Stretched-Exponential Relaxation with Continuously Varying Stretching Exponents

Relaxation in glasses is typically approximated by a stretched-exponential form: $f(t) = A \exp [-(t/\tau )^{\beta }]$ . Here, we show that the relaxation in a model of sheared non-Brownian suspensions developed by Corté et al. (Nat Phys 4:420–424, 2008 ) can be well approximated by a stretched exponential with an exponent $\beta $ that depends on the strain amplitude: $0.25< \beta < 1$ . In a one-dimensional version of the model, we show how the relaxation originates from density fluctuations in the initial particle configurations. Our analysis is in good agreement with numerical simulations and reveals a functional form for the relaxation that is distinct from, but well approximated by, a stretched-exponential function.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗

DOE OSTI2014

On the Composition of Bimetallic Near-Surface Alloys in the Presence of Oxygen and Carbon Monoxide

Periodic, self-consistent density functional theory calculations (GGA-PW91) are used to examine surface segregation in close-packed bimetallic Pt-overlayer alloy surfaces (Pt*/M, M = Au, Ag, Cu, Pd, Ir, Rh, Os, Ru, and Re) in different environments. In particular, we find that the thermodynamically stable surface termination in these Pt*/M alloys can be inverted from Pt-terminated in vacuum to M-terminated under exposure to oxygen (for an M that is more oxophillic than Pt). Interestingly, in many of these alloys, Pt is not driven into the bulk, rather it remains in the first subsurface layer where it enhances oxygen binding through a ligand interaction with the surface metal atoms. On the other hand, exposure to CO provides a much milder driving force for the surface composition inversion. To quantify segregation under catalytically relevant conditions, we constructed approximate phase diagrams for the PtRu system as a function of O2 and CO chemical potential (temperature, pressure). The results show that the surface termination inverts with many orders of magnitude higher CO pressure than with O2.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2022

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Optical properties of a nanostructured glass-based film using spectroscopic ellipsometry

Nanostructured glass films, which are fabricated using spinodally phase-separated low-alkali glasses, have several interesting and useful characteristics, including being robust, non-wetting and antireflective. Spectroscopic ellipsometry measurements have been performed on one such film and its optical properties were analyzed using a 5-layer structural model of the near-surface region. Since the glass and the film are transparent over the spectral region of the measurement, the Sellmeier model is used to parameterize the dispersion in the refractive index. To simulate the variation of the optical properties of the film over the spot size of the ellipsometer (~ 3 × 5 mm), the Sellmeier amplitude is convoluted using a Gaussian distribution. The transition layers between the ambient and the film and between the film and the substrate are modeled as graded layers, where the refractive index varies as a function of depth. These layers are modeled using a two-component Bruggeman effective medium approximation where the two components are the layer above and the layer below. Lastly, the fraction is continuous through the transition layer and is modelled using the incomplete beta function.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Improving Estimation of the Koopman Operator with Kolmogorov–Smirnov Indicator Functions

It has become common to perform kinetic analysis using approximate Koopman operators that transform high-dimensional timeseries of observables into ranked dynamical modes. The key to the practical success of the approach is the identification of a set of observables that form a good basis on which to expand the slow relaxation modes. Good observables are, however, difficult to identify a priori and suboptimal choices can lead to significant underestimations of characteristic time scales. Leveraging the representation of slow dynamics in terms of Hidden Markov Models (HMM), we propose a simple and computationally efficient clustering procedure to infer surrogate observables that form a good basis for slow modes. Here, we apply the approach to an analytically solvable model system as well as on three protein systems of different complexities. We consistently demonstrate that the inferred indicator functions can significantly improve the estimation of the leading eigenvalues of Koopman operators and correctly identify key states and transition time scales of stochastic systems, even when good observables are not known a priori.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2011

Strongly correlated states of a small cold-atom cloud from geometric gauge fields

Using exact diagonalization for a small system of cold bosonic atoms, we analyze the emergence of strongly correlated states in the presence of an artificial magnetic field. This gauge field is generated by a laser beam that couples two internal atomic states, and it is related to Berry's geometrical phase that emerges when an atom follows adiabatically one of the two eigenstates of the atom-laser coupling. Our approach allows us to go beyond the adiabatic approximation, and to characterize the generalized Laughlin wave functions that appear in the strong magnetic-field limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2011

Three-point correlators: Examples from Lunin-Maldacena background

Recently there has been progress on the calculation of three-point functions with two 'heavy' operators via semiclassical methods. We extend this analysis to the case of the Lunin-Maldacena background and examine the suggested procedure for several simple string solutions. By making use of AdS/CFT duality, we derive the relevant correlation functions of operators belonging to the dual N=1 superconformal gauge theory, and recover an important relation connecting the dimensions of heavy states and the structure constants.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2015

Debye temperature, thermal expansion, and heat capacity of TcC up to 100 GPa

Highlights: • A number of thermodynamic properties of rocksalt TcC are investigated for the first time. • The quasi-harmonic Debye model is applied to take into account the thermal effect. • The pressure and temperature up to about 100 GPa and 3000 K, respectively. - Abstract: Debye temperature, thermal expansion coefficient, and heat capacity of ideal stoichiometric TcC in the rocksalt structure have been studied systematically by using ab initio plane-wave pseudopotential density functional theory method within the generalized gradient approximation. Through the quasi-harmonic Debye model, in which the phononic effects are considered, the dependences of Debye temperature, thermal expansion coefficient, constant-volume heat capacity, and constant-pressure heat capacity on pressure and temperature are successfully predicted. All the thermodynamic properties of TcC with rocksalt phase have been predicted in the entire temperature range from 300 to 3000 K and pressure up to 100 GPa.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Close-contact melting on an isothermal surface with the inclusion of non-Newtonian effects

The present study deals with a theoretical investigation of a close-contact melting (CCM) process involving a vertical cylinder on a horizontal isothermal surface, where the liquid phase is a non-Newtonian viscoplastic fluid that behaves according to the Bingham model. Accordingly, a new approach is formulated based on the thin layer approximation and different quasi-steady process assumptions. By analytical derivation, an algebraic equation that relates the molten layer thickness and the solid bulk height is developed. The problem is then solved numerically, coupled with another equation for the melting rate. The new model shows that as the yield stress increases the melting rate decreases and the molten layer thickness increases. It is found that under certain conditions, the model can be reduced to a form that allows an analytical solution. The approximate model predicts an exponential dependence of both the melt fraction and the molten layer thickness. Comparison between the numerical and analytical solutions shows that the analytical approximation provides an excellent estimation for sufficiently large values of the yield stress. Dimensional analysis, which is supported by the analytical model results, reveals the dimensionless groups that govern the problem. For the general case, the melt fraction is a function of two dimensionless groups. For the analytical approximation, it is shown that the melt fraction is governed by a single dimensionless group and that the molten layer thickness is governed by two dimensionless groups.

Mechanics↗