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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

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

Shock compression of single-crystal austenitic FeCr18Ni12.5 stainless steel to 60 GPa

We measured the austenitic FeCr18Ni12.5 stainless steel Hugoniot as a function of crystallographic direction to approximately 60 GPa. We shock-compressed FeCr18Ni12.5 samples oriented along $\langle$100$\rangle$, $\langle$110$\rangle$, and $\langle$111$\rangle$ to mean stresses ranging 30.5–58.1 GPa via Ta plate impact in a large-bore powder gun and measured the free-surface velocities with laser interferometry. We unambiguously observed the largest post-shock free-surface velocity along $\langle$100$\rangle$ in each experiment, which consequently produced the lowest shock velocity along that orientation. However, the propagation of experimental uncertainties through the impedance matching scheme used to compute the shock velocity produced sufficient uncertainty overlap to preclude definitive conclusion of Hugoniot anisotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects of artificial collisions, filtering, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations

Kinetic simulations of collisionless plasmas are computationally challenging due to phase-space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related to limited velocity resolution in Hermite-based Vlasov–Poisson simulations: artificial collisions, filtering, and nonlocal closure approaches. We evaluate each method's performance in approximating the linear kinetic response function and suppressing recurrence in linear and nonlinear regimes. Numerical simulations of Landau damping demonstrate that artificial collisions, particularly higher orders of the Lenard-Bernstein collisional operator, most effectively recover the correct damping rate across a range of wavenumbers. Moreover, Hou-Li filtering and nonlocal closures underdamp high wavenumber modes in linear simulations, and the Lenard-Bernstein collisional operator overdamps low wavenumber modes in both linear and nonlinear simulations. This study demonstrates that hypercollisions offer a robust approach to kinetic simulations, accurately capturing collisionless dynamics with limited velocity resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Structure and permeation mechanism of a mammalian urea transporter

As an adaptation to infrequent access to water, terrestrial mammals produce urine that is hyperosmotic to plasma. To prevent osmotic diuresis by the large quantity of urea generated by protein catabolism, the kidney epithelia contain facilitative urea transporters (UTs) that allow rapid equilibration between the urinary space and the hyperosmotic interstitium. Here we report the first X-ray crystal structure of a mammalian UT, UT-B, at a resolution of 2.36 {angstrom}. UT-B is a homotrimer and each protomer contains a urea conduction pore with a narrow selectivity filter. Structural analyses and molecular dynamics simulations showed that the selectivity filter has two urea binding sites separated by an approximately 5.0 kcal/mol energy barrier. Functional studies showed that the rate of urea conduction in UT-B is increased by hypoosmotic stress, and that the site of osmoregulation coincides with the location of the energy barrier.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Ab initio molecular dynamics of solvation effects on reactivity at electrified interfaces

Using ab initio molecular dynamics as implemented in periodic, self-consistent (generalized gradient approximation Perdew–Burke–Ernzerhof) density functional theory, we investigated the mechanism of methanol electrooxidation on Pt(111). In this work, we investigated the role of water solvation and electrode potential on the energetics of the first proton transfer step, methanol electrooxidation to methoxy (CH 3 O) or hydroxymethyl (CH 2 OH). The results show that solvation weakens the adsorption of methoxy to uncharged Pt(111), whereas the binding energies of methanol and hydroxymethyl are not significantly affected. The free energies of activation for breaking the C–H and O–H bonds in methanol were calculated through a Blue Moon Ensemble using constrained ab initio molecular dynamics. Calculated barriers for these elementary steps on unsolvated, uncharged Pt(111) are similar to results for climbing-image nudged elastic band calculations from the literature. Water solvation reduces the barriers for both C–H and O–H bond activation steps with respect to their vapor-phase values, although the effect is more pronounced for C–H bond activation, due to less disruption of the hydrogen bond network. The calculated activation energy barriers show that breaking the C–H bond of methanol is more facile than the O–H bond on solvated negatively biased or uncharged Pt(111). However, with positive bias, O–H bond activation is enhanced, becoming slightly more facile than C–H bond activation.

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

MODELING COLLISIONAL CASCADES IN DEBRIS DISKS: STEEP DUST-SIZE DISTRIBUTIONS

We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi. Our model yields a quasi-steady-state slope of n(m) {approx} m{sup -1.88} [n(a) {approx} a{sup -3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle-mass distribution generates spectral energy distributions that are in agreement with the observed ones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Relativistic redshifts in quasar broad lines

The broad emission lines commonly seen in quasar spectra have velocity widths of a few percent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad Hβ line in the quasar rest frame (determined from the core component of the [O III] line) for over 20,000 quasars from the Sloan Digital Sky Survey Data Release 7 quasar catalog. The mean redshift as a function of line width is approximately consistent with the relativistic redshift that is expected if the line originates in a randomly oriented Keplerian disk that is obscured when the inclination of the disk to the line of sight exceeds ∼30°-45°, consistent with simple active galactic nucleus unification schemes. This result also implies that the net line-of-sight inflow/outflow velocities in the broad-line region are much less than the Keplerian velocity when averaged over a large sample of quasars with a given line width.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Comparative genome analyses suggest a hemibiotrophic lifestyle and virulence differences for the beech bark disease fungal pathogens Neonectria faginata and Neonectria coccinea

Abstract Neonectria faginata and Neonectria coccinea are the causal agents of the insect-fungus disease complex known as beech bark disease (BBD), known to cause mortality in beech forest stands in North America and Europe. These fungal species have been the focus of extensive ecological and disease management studies, yet less progress has been made toward generating genomic resources for both micro- and macro-evolutionary studies. Here, we report a 42.1 and 42.7 mb highly contiguous genome assemblies of N. faginata and N. coccinea, respectively, obtained using Illumina technology. These species share similar gene number counts (12,941 and 12,991) and percentages of predicted genes with assigned functional categories (64 and 65%). Approximately 32% of the predicted proteomes of both species are homologous to proteins involved in pathogenicity, yet N. coccinea shows a higher number of predicted mitogen-activated protein kinase genes, virulence determinants possibly contributing to differences in disease severity between N. faginata and N. coccinea. A wide range of genes encoding for carbohydrate-active enzymes capable of degradation of complex plant polysaccharides and a small number of predicted secretory effector proteins, secondary metabolite biosynthesis clusters and cytochrome oxidase P450 genes were also found. This arsenal of enzymes and effectors correlates with, and reflects, the hemibiotrophic lifestyle of these two fungal pathogens. Phylogenomic analysis and timetree estimations indicated that the N. faginata and N. coccinea species divergence may have occurred at ∼4.1 million years ago. Differences were also observed in the annotated mitochondrial genomes as they were found to be 81.7 kb (N. faginata) and 43.2 kb (N. coccinea) in size. The mitochondrial DNA expansion observed in N. faginata is attributed to the invasion of introns into diverse intra- and intergenic locations. These first draft genomes of N. faginata and N. coccinea serve as valuable tools to increase our understanding of basic genetics, evolutionary mechanisms and molecular physiology of these two nectriaceous plant pathogenic species.

Salgado-Salazar, Catalina↗
DOE OSTI2021

Impacts of random filling on spin squeezing via Rydberg dressing in optical clocks

Here, we analyze spin squeezing via Rydberg dressing in optical lattice clocks with random fractional filling. We compare the achievable clock stability in different lattice geometries, including unity-filled tweezer clock arrays and fractionally filled lattice clocks with varying dimensionality. We provide practical considerations and useful tools in the form of approximate analytical expressions and fitting functions to aid in the experimental implementation of Rydberg-dressed spin squeezing. We demonstrate that spin squeezing via Rydberg dressing in one-, two-, and three-dimensional optical lattices can provide significant improvements in stability in the presence of random fractional filling.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2011

Superfluid-to-Mott-insulator transition in the one-dimensional Bose-Hubbard model for arbitrary integer filling factors

We study the quantum phase transition between the superfluid and the Mott insulator in the one-dimensional (1D) Bose-Hubbard model. Using the time-evolving block-decimation method, we numerically calculate the tunneling splitting of two macroscopically distinct states with different winding numbers. From the scaling of the tunneling splitting with respect to the system size, we determine the critical point of the superfluid-to-Mott-insulator transition for arbitrary integer filling factors. We find that the critical values versus the filling factor in 1D, 2D, and 3D are well approximated by a simple analytical function. We also discuss the condition for determining the transition point from a perspective of the instanton method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Meta-GGA performance in solids at almost GGA cost

A recent modification, r 2 SCAN, of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) exchange-correlation functional mostly eliminates numerical instabilities and attendant integration grid sensitivities exhibited by SCAN. Here, we show that the successful deorbitalization of SCAN to SCAN-L (SCAN with density Laplacian dependence) carries over directly to yield r 2 SCAN-L. A major benefit is that the high iteration counts that hindered the use of SCAN-L are eliminated in r 2 SCAN-L. It therefore is a computationally much faster meta-GGA than its orbital-dependent antecedent. In this work, validation data for molecular heats of formation, bond lengths, and vibration frequencies (G3/99X, T96-R, and T82-F test sets, respectively) and on lattice constants, and cohesive energies (for 55 solids) and bulk moduli (for 40 solids) are provided. In addition, we show that the overmagnetization of bcc Fe, hcp Co, and fcc Ni persists in r 2 SCAN but does not appear in r 2 SCAN-L. Distinct from SCAN, both r 2 SCAN and r 2 SCAN-L give the correct nonmagnetic ground state for bcc V.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Electronic Correlation and Magnetism in the Ferromagnetic Metal Fe 3 GeTe 2

Motivated by the search for design principles of rare-earth-free strong magnets, we present a study of electronic structure and magnetic properties of the ferromagnetic metal Fe3GeTe2 within local density approximation (LDA) of the density functional theory, and its combination with dynamical mean-field theory (DMFT). For comparison to these calculations, we have measured magnetic and thermodynamic properties as well as X-ray magnetic circular dichroism and the photoemission spectrum of single crystal Fe3GeTe2. We find that the experimentally determined Sommerfeld coefficient is enhanced by an order of magnitude with respect to the LDA value. This enhancement can be partially explained by LDA+DMFT. Additionally, the inclusion of dynamical electronic correlation effects provides the experimentally observed magnetic moments, and the spectral density is in better agreement with photoemission data. Lastly, these results establish the importance of electronic correlations in this ferromagnet.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Global fit of electron and neutrino elastic scattering data to determine the strange quark contribution to the vector and axial form factors of the nucleon

We present a global fit of neutral-current elastic (NCE) neutrino-scattering data and parity-violating electron-scattering (PVES) data with the goal of determining the strange quark contribution to the vector and axial form factors of the proton. Previous fits of this form included data from a variety of PVES experiments (PVA4, HAPPEx, G0, SAMPLE) and the NCE neutrino and anti-neutrino data from BNL E734. These fits did not constrain the strangeness contribution to the axial form factor G A s ( Q 2 ) at low Q 2 very well because there was no NCE data for Q 2 < 0.45 GeV 2 . Our new fit includes for the first time MiniBooNE NCE data from both neutrino and antineutrino scattering; this experiment used a hydrocarbon target and so a model of the neutrino interaction with the carbon nucleus was required. Three different nuclear models have been employed: a relativistic Fermi gas model, the superscaling approximation model, and a spectral function model. We find a tremendous improvement in the constraint of G A s ( Q 2 ) at low Q 2 compared to previous work, although more data is needed from NCE measurements that focus on exclusive single-proton final states, for example from MicroBooNE. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Viscous corrections to electromagnetic emissivities in QCD

We provide a general framework for the derivation of the hydrodynamical corrections to the QCD electromagnetic emissivities in a viscous fluid. Here, assuming that the emission times are short in comparison to the fluid evolution time, we show that the leading corrections in the fluid gradients are controlled by the bulk and shear tensors times pertinent response functions involving the energy-momentum tensor. In a hadronic fluid phase, we detail these contributions using spectral functions. Using the vector dominance approximation, we show that the bulk viscosity correction to the photon rate is sizable, while the shear viscosity is negligible for about all frequencies. In the partonic phase near the transition temperature, we provide an assessment of the viscous corrections to the photon and dilepton emissions, using a nonperturbative quark-gluon plasma with soft thermal gluonic corrections in the form of operators of leading mass dimension. Again, the thermal bulk viscosity corrections are found to be larger than the thermal shear viscosity corrections at all energies for both the photon and dilepton in the partonic phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Anisotropy of stress correlation in two-dimensional liquids and a pseudospin model

Liquids are condensed matter in which atoms are strongly correlated in position and momentum. The atomic pair density function (PDF) is used often in describing such correlation. However, elucidation of many properties requires higher degrees of correlation than the pair correlation. For instance, viscosity depends upon the stress correlations in space and time. We examine the cross correlation between the stress correlation at the atomic level and the PDF for two-dimensional liquids. We introduce the concept of the stress-resolved pair distribution function (SRPDF) that uses the sign of atomic-level stress as a selection rule to include particles from density correlations. The connection between SRPDFs and stress correlation function is explained through an approximation in which the shear stress is replaced by a pseudospin. Lastly, we further assess the possibility of interpreting the long-range stress correlation as a consequence of short-range Ising-like pseudospin interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Variational approach to studying solitary waves in the nonlinear Schrödinger equation with complex potentials

Here in this paper, we discuss the behavior of solitary wave solutions of the nonlinear Schrödinger equation (NLSE) as they interact with complex potentials, using a four-parameter variational approximation based on a dissipation functional formulation of the dynamics. We concentrate on spatially periodic potentials with the periods of the real and imaginary part being either the same or different. Our results for the time evolution of the collective coordinates of our variational ansatz are in good agreement with direct numerical simulation of the NLSE. We compare our method with a collective coordinate approach of Kominis and give examples where the two methods give qualitatively different answers. In our variational approach, we are able to give analytic results for the small oscillation frequency of the solitary wave oscillating parameters which agree with the numerical solution of the collective coordinate equations. We also verify that instabilities set in when the slope dp(t)/dv(t) becomes negative when plotted parametrically as a function of time, where p(t) is the momentum of the solitary wave and v(t) the velocity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Unifying first-principles theoretical predictions and experimental measurements of size effects in thermal transport in SiGe alloys

We demonstrate the agreement between first-principles calculations and experimental measurements of size effects in thermal transport in SiGe alloys without fitting parameters. Transient thermal grating (TTG) is used to measure the effect of the grating period on the temperature decay. The virtual crystal approximation under the density-functional-theory framework combined with impurity scattering is used to determine the phonon properties for the exact alloy composition of the measured samples. With these properties, classical size effects are calculated for the experimental geometry of reflection mode TTG using the recently developed variational solution to the phonon Boltzmann transport equation, which is verified against established Monte Carlo simulations. Here, we find agreement between theoretical predictions and experimental measurements in the reduction of thermal conductivity (as much as fourfold of the bulk value) across grating periods spanning one order of magnitude. This paper provides a framework for the study of size effects in thermal transport in opaque materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mirror-centered representation of a focusing hyperbolic mirror for X-ray beamlines

Conic sections are commonly used in reflective X-ray optics. Hyperbolic mirrors can focus a converging light source and are frequently paired with elliptical or parabolic mirrors in Wolter type configurations. This paper derives the closed-form expression for a mirror-centered hyperbolic shape, with zero-slope at the origin. Combined with the slope and curvature, such an expression facilitates metrology, manufacturing and mirror-bending calculations. Previous works consider ellipses, parabolas, magnifying hyperbolas or employ lengthy approximations. Here, the exact shape function is given in terms of the mirror incidence angle and the source and image distances.

47 OTHER INSTRUMENTATION↗