Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An Evaluation of Bang-Time Measurements from a Multichannel, Triaxial, nTOF Diagnostic for MagLIF Experiments at the Z facility

Neutron bang times for a series of MagLIF (Magnetic Liner Inertial Fusion) experiments with D 2 -filled targets have been measured at the Z facility. The emitted neutrons were detected as current-mode pulses in a multichannel, neutron time-of-flight (nTOF) diagnostic with conventional, scintillator-photomultiplier-tube (PMT) detectors. In these experiments, the detectors were fielded at known, fixed distances L (690-2510 cm) from the target, and on three, non-coplanar (but convergent) lines-of-sight (LOS). The primary goal of this diagnostic was to estimate a fiducial time (bang time) relative to an externally generated time-base for synchronizing all the diagnostics in an experiment. Recorded arrival times (A7) of the pulses were characterized experimentally by three numerical methods: a first-moment estimate (centroid) and two nodal measures — Savitzky-Golay (SG) smoothing and a single point peak estimate of the raw data. These times were corrected for internal detector time delays (transit and impulse-response function) — an adjustment that linked the recorded ATs to the corresponding arrival of uncollided neutrons at each detector. The bang time was then estimated by linearly regressing the arrival times against the associated distances to the source; t bang (on the system timescale) was taken as the temporal intercept of the regression equation at distance L = 0. This article reports the analysis for a representative shot #2584 for which (a) the recorded ATs — even without detector corrections — agreed by method in each channel to within 1-2 ns; (b) internal corrections were each ~3 — 5 ns; and (c) a 95% uncertainty (confidence) interval for t bang in this shot was estimated at ±3 ns with 4 degrees of freedom. A secondary goal for this diagnostic was to check that the bang time measurements corresponded to neutrons emitted by the D(d,n) 3 He reaction in a thermalized DD plasma. According to the theoretical studies by Brysk, such neutrons should be emitted with an isotropic Gaussian distribution of mean kinetic energy $ \overline{E}$ of 2.449 MeV; this energy translates to a mean neutron speed $ \overline{u}$ of 2.160 cm/ns [D. H. Munro, Nuclear Fusion, 56(3) 036001 (2016)]. In the MagLIF series of shots there was no evidence of spatial asymmetry in the time-distance regressions, and it was possible to extract the mean neutron speed from the slope of these fits. In shot 2584 $ \overline{u}$ was estimated at 2.152 cm/ns ± 0.010 cm/ns [95 % confidence, 4 dof] and the mean kinetic energy $ \overline{E}$ (with relativistic corrections) was 2.431 MeV ± 0.022 MeV [95 % confidence, 4 dof] — results supporting the assumption that D-D neutrons were, in fact, measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗
DOE OSTI2015

PROBABILISTIC MODEL OF BEAM–PLASMA INTERACTION IN RANDOMLY INHOMOGENEOUS PLASMA

We propose a new model that describes beam–plasma interaction in the presence of random density fluctuations with a known probability distribution. We use the property that, for the given frequency, the probability distribution of the density fluctuations uniquely determines the probability distribution of the phase velocity of waves. We present the system as discrete and consisting of small, equal spatial intervals with a linear density profile. This approach allows one to estimate variations in wave energy density and particle velocity, depending on the density gradient on any small spatial interval. Because the characteristic time for the evolution of the electron distribution function and the wave energy is much longer than the time required for a single wave–particle resonant interaction over a small interval, we determine the description for the relaxation process in terms of averaged quantities. We derive a system of equations, similar to the quasi-linear approximation, with the conventional velocity diffusion coefficient D and the wave growth rate γ replaced by the average in phase space, by making use of the probability distribution for phase velocities and by assuming that the interaction in each interval is independent of previous interactions. Functions D and γ are completely determined by the distribution function for the amplitudes of the fluctuations. For the Gaussian distribution of the density fluctuations, we show that the relaxation process is determined by the ratio of beam velocity to plasma thermal velocity, the dispersion of the fluctuations, and the width of the beam in the velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗
DOE OSTI2018

Particle image velocimetry (PIV) uncertainty quantification using moment of correlation (MC) plane

Here, we present a new uncertainty estimation method for particle image velocimetry (PIV), that uses the correlation plane as a model for the probability density function (PDF) of displacements and calculates the second order moment of the correlation (MC). The cross-correlation between particle image patterns is the summation of all particle matches convolved with the apparent particle image diameter. MC uses this property to estimate the PIV uncertainty from the shape of the cross-correlation plane. In this new approach, the generalized cross-correlation (GCC) plane corresponding to a PIV measurement is obtained by removing the particle image diameter contribution. The GCC primary peak represents a discretization of the displacement PDF, from which the standard uncertainty is obtained by convolving the GCC plane with a Gaussian function. Then a Gaussian least-squares-fit is applied to the peak region, accounting for the stretching and rotation of the peak, due to the local velocity gradients and the effect of the convolved Gaussian. The MC method was tested with simulated image sets and the predicted uncertainties show good sensitivity to the error sources and agreement with the expected RMS error. Subsequently, the method was demonstrated in three PIV challenge cases and two experimental datasets and was compared with the published image matching (IM) and correlation statistics (CS) techniques. Results show that the MC method has a better response to spatial variation in RMS error and the predicted uncertainty is in good agreement with the expected standard uncertainty. The uncertainty prediction was also explored as a function of PIV interrogation window size. Overall, the MC method performance establishes itself as a valid uncertainty estimation tool for planar PIV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

A web-based tool for the prediction of rice transcription factor function

Transcription factors (TFs) are an important class of regulatory molecules. Despite their importance, only a small number of genes encoding TFs have been characterized in Oryza sativa (rice), often because gene duplication and functional redundancy complicate their analysis. To address this challenge, we developed a web-based tool called the Rice Transcription Factor Phylogenomics Database (RTFDB) and demonstrate its application for predicting TF function. The RTFDB hosts transcriptome and co-expression analyses. Sources include high-throughput data from oligonucleotide microarray (Affymetrix and Agilent) as well as RNA-Seq-based expression profiles. We used the RTFDB to identify tissue-specific and stress-related gene expression. Subsequently, 273 genes preferentially expressed in specific tissues or organs, 455 genes showing a differential expression pattern in response to 4 abiotic stresses, 179 genes responsive to infection of various pathogens and 512 genes showing differential accumulation in response to various hormone treatments were identified through the meta-expression analysis. Pairwise Pearson correlation coefficient analysis between paralogous genes in a phylogenetic tree was used to assess their expression collinearity and thereby provides a hint on their genetic redundancy. Integrating transcriptome with the gene evolutionary information reveals the possible functional redundancy or dominance played by paralog genes in a highly duplicated genome such as rice. With this method, we estimated a predominant role for 83.3% (65/78) of the TF or transcriptional regulator genes that had been characterized via loss-of-function studies. In this regard, the proposed method is applicable for functional studies of other plant species with annotated genome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Hyperspherical Sparse Approximation Techniques for High-Dimensional Discontinuity Detection

This work proposes a hyperspherical sparse approximation framework for detecting jump discontinuities in functions in high-dimensional spaces. The need for a novel approach results from the theoretical and computational inefficiencies of well-known approaches, such as adaptive sparse grids, for discontinuity detection. Our approach constructs the hyperspherical coordinate representation of the discontinuity surface of a function. Then sparse approximations of the transformed function are built in the hyperspherical coordinate system, with values at each point estimated by solving a one-dimensional discontinuity detection problem. Due to the smoothness of the hypersurface, the new technique can identify jump discontinuities with significantly reduced computational cost, compared to existing methods. Several approaches are used to approximate the transformed discontinuity surface in the hyperspherical system, including adaptive sparse grid and radial basis function interpolation, discrete least squares projection, and compressed sensing approximation. Moreover, hierarchical acceleration techniques are also incorporated to further reduce the overall complexity. In conclusion, rigorous complexity analyses of the new methods are provided, as are several numerical examples that illustrate the effectiveness of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Revisiting the Redshift Distribution of Gamma-Ray Bursts in the Swift Era

Le and Dermer developed a gamma-ray burst (GRB) model to fit the redshift and the jet opening angle distributions measured with pre-Swift and Swift missions and showed that GRBs do not follow the star formation rate. Their fitted results were obtained without the opening angle distribution from Swift with an incomplete Swift sample, and the calculated jet opening angle distribution was obtained by assuming a flat νF{sub ν} spectrum. In this paper, we revisit the work done by Le and Dermer with an assumed broken power law GRB spectrum. Utilizing more than 100 GRBs in the Swift sample that include both the observed estimated redshifts and jet opening angles, we obtain a GRB burst rate functional form that gives acceptable fits to the pre-Swift and Swift redshift and jet opening angle distributions with an indication that an excess of GRBs exists at low redshift below z≈2. The mean redshifts and jet opening angles for pre-Swift (Swift) are 〈z〉∼1.3 (1.7) and 〈θ{sub j}〉∼7{sup ∘} (11{sup ∘}), respectively. Assuming a GRB rate density (SFR9), similar to the Hopkins and Beacom star formation history and as extended by Li, the fraction of high-redshift GRBs is estimated to be below 10% and 5% at z⩾4 and z⩾5, respectively, and below 10% at z⩽1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Approximate Bayesian Model Inversion for PDEs with Heterogeneous and State-Dependent Coefficients

We present two approximate Bayesian inference methods for parameter estimation in partial differential equation (PDE) models with space-dependent and state-dependent parameters. We demonstrate that these methods provide accurate and cost-effective alternatives to Markov Chain Monte Carlo simulation. We assume a parameterized Gaussian prior on the unknown functions, and approximate the posterior density by a parameterized multivariate Gaussian density. The parameters of the prior and posterior are estimated from sparse observations of the PDE model's states and the unknown functions themselves by maximizing the evidence lower bound (ELBO), a lower bound on the log marginal likelihood of the observations. The first method, Laplace-EM, employs the expectation maximization algorithm to maximize the ELBO, with a Laplace approximation of the posterior on the E-step, and minimization of a Kullback-Leibler divergence on the M-step. The second method, DSVI-EB, employs the doubly stochastic variational inference (DSVI) algorithm, in which the ELBO is maximized via gradient-based stochastic optimization, with nosiy gradients computed via simple Monte Carlo sampling and Gaussian backpropagation. We apply these methods to identifying diffusion coefficients in linear and nonlinear diffusion equations, and we find that both methods provide accurate estimates of posterior densities and the hyperparameters of Gaussian priors. While the Laplace-EM method is more accurate, it requires computing Hessians of the physics model. The DSVI-EB method is found to be less accurate but only requires gradients of the physics model.

approximate Bayesian inference, model inversion, v↗
DOE OSTI2024

Coexistence of low and high spin states in La 18 Co 28 Pb 3

The electronic structure and magnetic properties of a newly predicted stable ternary compound La 18 ⁢Co 28 ⁢Pb 3 are studied using electronic structure analysis. The ground state of this compound is ferromagnetic, with three positions of nonequivalent magnetic Co atoms. A strong dependence of magnetic properties on volume shows that this system is situated near the point of magnetic instability. The coexistence of high- and low-spin ferromagnetic states as a function of volume near equilibrium was discovered. A corresponding spin tunneling splitting was estimated. The stability of the theoretically predicted magnetic ground state was tested by varying the Hubbard parameter. The thermal spin fluctuations were added to estimate the paramagnetic moment and a Curie temperature. Finally, the necessity of experimental verification of the obtained results is emphasized.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A functional global sensitivity measure and efficient reliability sensitivity analysis with respect to statistical parameters

Sensitivity analysis and reliability assessment are two important aspects of structural and system safety. Epistemic uncertainty with respect to probabilistic model of input parameters due to lack of knowledge is present in many scarce-data applications and complicates the characterization of uncertainty in model response. In this article, we present two importance measures to evaluate the impact of distribution parameters on the probability distribution function (PDF) of the output and the failure probability. The epistemic uncertainty associated with the distribution parameters is modeled as random variables. Additionally, a modified extended polynomial chaos expansion (MEPCE) approach is introduced in which aleatory and epistemic random variables are modeled and propagated simultaneously while allowing the separate assessment for any single epistemic variable. A MEPCE-based kernel density estimation (KDE) construction provides a composite map from each epistemic variable to the response PDF. The functional global sensitivity index of the PDF with respect to the distribution parameters is thus derived, as a function of output, which is both more informative and more efficient than standard scalar sensitivity measures. Reliability sensitivity indices can be readily evaluated by integrating the global sensitivity index function over the failure zone. Three illustrative examples are used to demonstrate the proposed methodology.

42 ENGINEERING↗
DOE OSTI2018

Prototype Stilbene Neutron Collar II

A neutron collar using stilbene organic scintillator cells for fast neutron counting is described for the assay of fresh low enriched uranium (LEU) fuel assemblies. The prototype stilbene collar has a form factor similar to standard He-3 based collars and uses an AmLi interrogation neutron source. This report describes the simulation, using a custom LLNL developed MCNPX, of list mode neutron correlation data on various fuel assemblies. Some of these fuel assemblies have neutron absorbers (burnable Gd poisons) and partial defect s (DU, Fe, Pb pins). Calibration curves (doubles vs 235 U linear mass density) are presented for fast (with Cd lining) mode of operation. The modeling incorporates LLNL developed neutron crosstalk correction algorithms and stilbene light output detector response quench function. Pileup issues are addressed using the LLNL code GamGen / RadSrc to estimate the flux of alpha decay gammas from the fuel assemblies on the stilbene cells with and without Pb shielding. One panel of the prototype collar using 10 high-efficiency 4”x2” stilbene cells from Inrad Optics, using LLNL developed technology, has been constructed and calibrated for optimal PSD and data acquisition electronics. Several sources (AmLi, AmBe, Cf, DD, and DT) have been measured with the single panel and the measured data is consistent with the modeling. A novel part of our work was to determine the light output (LO) quench function for the Inrad 4” x 2” stilbene crystals through a LLNL-MCNPX simulation based optimized fitting of the measured LO data.

42 ENGINEERING↗
DOE OSTI2015

SU-E-T-422: Fast Analytical Beamlet Optimization for Volumetric Intensity-Modulated Arc Therapy

Purpose: To implement a fast optimization algorithm on CPU/GPU heterogeneous computing platform and to obtain an optimal fluence for a given target dose distribution from the pre-calculated beamlets in an analytical approach. Methods: The 2D target dose distribution was modeled as an n-dimensional vector and estimated by a linear combination of independent basis vectors. The basis set was composed of the pre-calculated beamlet dose distributions at every 6 degrees of gantry angle and the cost function was set as the magnitude square of the vector difference between the target and the estimated dose distribution. The optimal weighting of the basis, which corresponds to the optimal fluence, was obtained analytically by the least square method. Those basis vectors with a positive weighting were selected for entering into the next level of optimization. Totally, 7 levels of optimization were implemented in the study.Ten head-and-neck and ten prostate carcinoma cases were selected for the study and mapped to a round water phantom with a diameter of 20cm. The Matlab computation was performed in a heterogeneous programming environment with Intel i7 CPU and NVIDIA Geforce 840M GPU. Results: In all selected cases, the estimated dose distribution was in a good agreement with the given target dose distribution and their correlation coefficients were found to be in the range of 0.9992 to 0.9997. Their root-mean-square error was monotonically decreasing and converging after 7 cycles of optimization. The computation took only about 10 seconds and the optimal fluence maps at each gantry angle throughout an arc were quickly obtained. Conclusion: An analytical approach is derived for finding the optimal fluence for a given target dose distribution and a fast optimization algorithm implemented on the CPU/GPU heterogeneous computing environment greatly reduces the optimization time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

In-Home Mobility Frequency and Stability in Older Adults Living Alone With or Without MCI: Introduction of New Metrics

Background: Older adults spend a considerable amount of time inside their residences; however, most research investigates out-of-home mobility and its health correlates. We measured indoor mobility using room-to-room transitions, tested their psychometric properties, and correlated indoor mobility with cognitive and functional status. Materials and Methods: Community-dwelling older adults living alone ( n = 139; age = 78.1 ± 8.6 years) from the Oregon Center for Aging & Technology (ORCATECH) and Minority Aging Research Study (MARS) were included in the study. Two indoor mobility features were developed using non-parametric parameters (frequency; stability): Indoor mobility frequency (room-to-room transitions/day) was detected using passive infrared (PIR) motion sensors fixed on the walls in four geographic locations (bathroom; bedroom; kitchen; living room) and using door contact sensors attached to the egress door in the entrance. Indoor mobility stability was estimated by variances of number of room-to-room transitions over a week. Test-retest reliability (Intra-class coefficient, ICC) and the minimal clinically important difference (MCID) defined as the standard error of measurement (SEM) were generated. Generalized estimating equations models related mobility features with mild cognitive impairment (MCI) and functional status (gait speed). Results: An average of 206 days (±127) of sensor data were analyzed per individual. Indoor mobility frequency and stability showed good to excellent test-retest reliability (ICCs = 0.91[0.88–0.94]; 0.59[0.48–0.70]). The MCIDs of mobility frequency and mobility stability were 18 and 0.09, respectively. On average, a higher indoor mobility frequency was associated with faster gait speed (β = 0.53, p = 0.04), suggesting an increase of 5.3 room-to-room transitions per day was associated with an increase of 10 cm/s gait speed. A decrease in mobility stability was associated with MCI (β = −0.04, p = 0.03). Discussion: Mobility frequency and stability in the home are clinically meaningful and reliable features. Pervasive-sensing systems deployed in homes can objectively reveal cognitive and functional status in older adults who live alone.

Wu, Chao-Yi↗
DOE OSTI2023

MODFLOW6 models used to evaluate potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, Southern Nevada

The hydrograph for well ER-5-3-2 in Frenchman Flat, southern Nevada, has previously unexplained water-level fluctuations. Four, three-dimensional, groundwater models (MODFLOW 6) were developed to evaluate potential stresses and hydrologic conditions affecting the well ER-5-3-2 hydrograph. Four model scenarios were developed that simulated: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) shallow (low transmissivity) and deep (high transmissivity) carbonate rocks, and (4) lateral heterogeneity of carbonate rocks. Input and output files for the four model scenarios are in the model and output directories, respectively. Hydraulic conductivity, specific storage, and wellbore-leakage rates (when simulated) were estimated with parameter estimation (PEST) by minimizing a weighted composite, sum-of-squares objective function. The objective function was informed by measurement and Tikhonov regularization observations. Measurement observations included drawdowns from the constant-rate aquifer test and water-level altitudes measured in well ER-5-3-2 from 2001-2021. Tikhonov regularization informed hydraulic conductivity and specific storage parameters that were insensitive to measurement observations, where homogeneity was the preferred relation. Batch files, executables, and MODFLOW 6, PEST, and post-processing utilities are in the ancillary directory. Supplementary data also are included in the ancillary directory, including site information, high-frequency water-level and aquifer-test data, transmissivity estimates, water-chemistry data, and water-temperature analyses. This USGS data release contains data, analyses, and model files for the simulations and analysis results described in U.S. Geological Survey Scientific Investigations Report (https://doi.org/10.3133/sir20225132).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

On the inward drift of runaway electrons during the plateau phase of runaway current

The well observed inward drift of current carrying runaway electrons during runaway plateau phase after disruption is studied by considering the phase space dynamic of runaways in a large aspect ratio toroidal system. We consider the case where the toroidal field is unperturbed and the toroidal symmetry of the system is preserved. The balance between the change in canonical angular momentum and the input of mechanical angular momentum in such a system requires runaways to drift horizontally in configuration space for any given change in momentum space. The dynamic of this drift can be obtained by integrating the modified Euler-Lagrange equation over one bounce time. It is then found that runaway electrons will always drift inward as long as they are decelerating. This drift motion is essentially non-linear, since the current is carried by runaways themselves, and any runaway drift relative to the magnetic axis will cause further displacement of the axis itself. A simplified analytical model is constructed to describe such inward drift both in the ideal wall case and no wall case, and the runaway current center displacement as a function of parallel momentum variation is obtained. The time scale of such displacement is estimated by considering effective radiation drag, which shows reasonable agreement with the observed displacement time scale. Furthermore, this indicates that the phase space dynamic studied here plays a major role in the horizontal displacement of runaway electrons during plateau phase. (C) 2016 AIP Publishing LLC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

First-principles analysis of anharmonic nuclear motion and thermal transport in thermoelectric materials

We show a first-principles approach for analyzing anharmonic properties of lattice vibrations in solids. We firstly extract harmonic and anharmonic force constants from accurate first-principles calculations based on the density functional theory. Using the many-body perturbation theory of phonons, we then estimate the phonon scattering probability due to anharmonic phonon-phonon interactions. We show the validity of the approach by computing the lattice thermal conductivity of Si, a typical covalent semiconductor, and selected thermoelectric materials PbTe and Bi{sub 2}Te{sub 3} based on the Boltzmann transport equation. We also show that the phonon lifetime and the lattice thermal conductivity of the high-temperature phase of SrTiO{sub 3} can be estimated by employing the perturbation theory on top of the solution of the self-consistent phonon equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Real-space atomic dynamics in metallic liquids investigated by inelastic neutron scattering

Understanding the dynamics of liquids at the atomic level remains a major challenge. Even though viscosity is one of the most fundamental properties of liquids, its atomistic origin is not fully elucidated. Through inelastic neutron scattering experiment on levitated metallic liquid droplets, the time-dependent pair correlation function, the Van Hove function, was determined for Z⁢r 50 ⁢C⁢u 50 and Z⁢r 80 ⁢P⁢t 20 liquids at various temperatures. The time for change in local atomic connectivity, τ LC , which is the timescale of atomic bond cutting and forming, is estimated based on the exponential decay of the nearest neighbor peak of the Van Hove function. At high temperatures above the crossover temperature T A , τ LC is equal to the Maxwell relaxation time, τ M =η/G ∞ , where η is the macroscopic shear viscosity and G ∞ is the high-frequency shear modulus. Below T A the ratio of τ M /τ LC increases with decreasing temperature, indicating increased atomic cooperativity as predicted by molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗