Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling statistics”

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 343 records · Page 19

Structural Studies of Amorphous Materials by Fluctuation Electron Microscopy

Fluctuation Electron Microscopy (FEM) is a technique that examines the fluctuations in electron scattering across a uniformly thin amorphous sample. The statistics of the intensity fluctuations, mean and variance, reveal any underlying medium-range order present in the structure. The goals of this project were: (1) To determine the fundamentals of the scattering physics that gives rise to the variance signal in fluctuation electron microscopy (FEM); (2) To use these discoveries to find ways to quantify FEM; (3) To apply the FEM method to interesting and technologically important families of amorphous materials, particularly those with important applications in energy-related processes. Excellent progress was made in items (1) and (2). In stage (3) we did not examine the metamict zircons, as proposed. Instead, we examined films of polycrystalline and amorphous semi-conducting diamond. Significant accomplishments are: (1) A Reverse Monte Carlo procedure was successfully implemented to invert FEM data into a structural model. This is computer-intensive, but it demonstrated that diffraction and FEM data from amorphous silicon are most consistent with a paracrystallite model. This means that there is more diamond-like topology present in amorphous silicon than is predicted by the continuous random network model. (2) There is significant displacement decoherence arising in diffraction from amorphous silicon and carbon. The samples are being bombarded by the electron beam and atoms do not stay still while being irradiated – much more than was formerly understood. The atom motions cause the destructive and constructive interferences in the diffraction pattern to fluctuate with time, and it is the time-averaged speckle that is being measured. The variance is reduced by a factor m, 4 ≤ m ≤ 1000, relative to that predicted by kinematical scattering theory. (3) Speckle intensity obeys a gamma distribution, where the mean intensity $ \overline{I}\ $ and m are the two parameters governing the shape of the gamma distribution profile. m is determined by the illumination spatial coherence, which is normally very high, and mostly by the displacement decoherence within the sample. (4) Amorphous materials are more affected by the electron beam than are crystalline materials. Different samples exhibit different disruptibility, as measured by the effective values of m that fit the data. (5) Understanding the origin of the displacement decoherence better should lead to efficient methods for computing the observed variance from amorphous materials.

36 MATERIALS SCIENCE↗

A statistical analysis of the Einstein normal galaxy sample. I - Spiral and irregular galaxies

The results of a statistical analysis of 48 spiral and irregular galaxies observed with the Einstein Observatory are reported. It is found that the X-ray luminosity is not directly correlated with the mass of the galaxies, but is strongly correlated with the blue luminosity. This suggests that most X-ray sources, including low-mass binaries, are binary systems belonging to both old disk and young arm Population I components. A real bulge population of X-ray sources exists that can contribute significantly to the X-ray luminosity of early-type spirals. The X-ray emission from a sample of 29 relatively isolated, normal elliptical and S0 galaxies is also studied, and the results are compared to those for the above sample to investigate the origin of the X-ray emission in early-type galaxies and the possible emission mechanisms. The influence of a powerful radio source and the onset of nuclear activity are also examined.

Fabbiano, G.↗

Field‐Induced Ferroelectric Phase Evolution During Polarization “Wake‐Up” in Hf 0.5 Zr 0.5 O 2 Thin Film Capacitors

Abstract As an emerging nonvolatile memory technology, HfO 2 ‐based ferroelectrics exhibit excellent compatibility with silicon CMOS process flows; however, the reliability of polarization switching in these materials remains a major challenge. During repeated field programming and erase of the polarization state of initially pristine HfO 2 ‐based ferroelectric capacitors, the magnitude of the measured polarization increases, a phenomenon known as “wake‐up”. In this study, the authors attempt to understand what causes the wake‐up effect in Hf 0.5 Zr 0.5 O 2 (HZO) capacitors using nondestructive methods that probe statistically significant sample volumes. Synchrotron X‐ray diffraction reveals a concerted shift in HZO Bragg peak position as a function of polarization switching cycle number in films prepared under conditions such that they exhibit extremely large (≈3000%) wake‐up. In contrast, a control sample with insignificant wake‐up shows no such peak shift. Capacitance – voltage measurements show evolution in the capacitance loop with switching cycle number for the wake‐up sample and no change for the control sample. Piezoresponse force microscopy measurements are utilized to visualize the domain switching with wake‐up. The combination of these observations clearly demonstrates that wake‐up is caused by a field‐driven phase transformation of the tetragonal phase to the metastable ferroelectric orthorhombic phase during polarization switching of HZO capacitors.

36 MATERIALS SCIENCE↗

GPU Acceleration of the Locally Selfconsistent Multiple Scattering Code for First Principles Calculation of the Ground State and Statistical Physics of Materials

The Locally Self-consistent Multiple Scattering (LSMS) code solves the first principles Density Functional theory Kohn-Sham equation for a wide range of materials with a special focus on metals, alloys and metallic nano-structures. It has traditionally exhibited near perfect scalability on massively parallel high performance computer architectures. We present our efforts to exploit GPUs to accelerate the LSMS code to enable first principles calculations of O(100,000) atoms and statistical physics sampling of finite temperature properties. Using the Cray XK7 system Titan at the Oak Ridge Leadership Computing Facility we achieve a sustained performance of 14.5PFlop/s and a speedup of 8.6 compared to the CPU only code.

Eisenbach, Markus↗

Strangeness enhancement at its extremes: multiple (multi-)strange hadron production in pp collisions at \(\sqrt{s}=5.02\) TeV

The probability to observe a specific number of strange and multi-strange hadrons (nS), denoted as P(nS), is measured by ALICE at midrapidity (|y| < 0.5) in $$\sqrt{s}=5.02$$ TeV proton-proton (pp) collisions, dividing events into several multiplicity-density classes. Exploiting, for the first time, a technique based on counting the number of strange-particle candidates event-by-event, this measurement allows one to extend the study of strangeness production beyond the mean of the distribution. This constitutes a new test bench for production mechanisms, probing events with a large imbalance between strange and non-strange content. The analysis of a large-statistics data sample makes it possible to extract P(nS) up to a maximum nS of 7 for $${\text{K}}_{\text{S}}^{0}$$, 5 for Λ and $$\overline{\Lambda }$$, 4 for Ξ− and $${\overline{\Xi } }^{+}$$, and 2 for Ω− and $${\overline{\Omega } }^{+}$$. From this, the probability of producing strange hadron multiplets per event is calculated, thereby enabling the extension of the study of strangeness enhancement to extreme situations where several strange quarks hadronize in a single event at midrapidity. Moreover, comparing hadron combinations with different u and d quark compositions and equal overall s quark content, the contribution to the enhancement pattern coming from non-strangeness related mechanisms is isolated. The results are compared with state-of-the-art phenomenological models implemented in commonly used Monte Carlo event generators, including PYTHIA 8 Monash 2013, PYTHIA 8 with QCD-based Color Reconnection and Rope Hadronization (QCD-CR + Ropes), and EPOS LHC, which incorporates both partonic interactions and hydrodynamic evolution. These comparisons show that the new approach dramatically enhances the sensitivity to the different underlying physics mechanisms modeled by each generator.

Abualrob, I J↗

The Surface-Topography Challenge: A Multi-Laboratory Benchmark Study to Advance the Characterization of Topography

Surface performance is critically influenced by topography in virtually all real-world applications. The current standard practice is to describe topography using one of a few industry-standard parameters. The most commonly reported number is Ra, the average absolute deviation of the height from the mean line (at some, not necessarily known or specified, lateral length scale). However, other parameters, particularly those that are scale-dependent, influence surface and interfacial properties; for example the local surface slope is critical for visual appearance, friction, and wear. The present Surface-Topography Challenge was launched to raise awareness for the need of a multi-scale description, but also to assess the reliability of different metrology techniques. In the resulting international collaborative effort, 153 scientists and engineers from 64 research groups and companies across 20 countries characterized statistically equivalent samples from two different surfaces: a “rough” and a “smooth” surface. The results of the 2088 measurements constitute the most comprehensive surface description ever compiled. We find wide disagreement across measurements and techniques when the lateral scale of the measurement is ignored. Consensus is established through scale-dependent parameters while removing data that violates an established resolution criterion and deviates from the majority measurements at each length scale. Our findings suggest best practices for characterizing and specifying topography. The public release of the accumulated data and presented analyses enables global reuse for further scientific investigation and benchmarking.

42 ENGINEERING↗

GPU acceleration of the Locally Selfconsistent Multiple Scattering code for first principles calculation of the ground state and statistical physics of materials

The Locally Self-consistent Multiple Scattering (LSMS) code solves the first principles Density Functional theory Kohn–Sham equation for a wide range of materials with a special focus on metals, alloys and metallic nano-structures. It has traditionally exhibited near perfect scalability on massively parallel high performance computer architectures. In this paper, we present our efforts to exploit GPUs to accelerate the LSMS code to enable first principles calculations of O(100,000) atoms and statistical physics sampling of finite temperature properties. We reimplement the scattering matrix calculation for GPUs with a block matrix inversion algorithm that only uses accelerator memory. Finally, using the Cray XK7 system Titan at the Oak Ridge Leadership Computing Facility we achieve a sustained performance of 14.5PFlop/s and a speedup of 8.6 compared to the CPU only code.

36 MATERIALS SCIENCE↗

The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations

A window-based steady-state detection algorithm has been developed for application to kinetic Monte Carlo simulation data. The algorithm, termed F-t-Pj-RG sequentially applies an F-test, a t-test, and a projection test on adjacent windows of the data while rolling (or shifting) and growing the windows when any of the tests fail. In aggregate, the algorithm is able to (a) automatically reject the warm-up period as not being at steady-state, as well as (b) determine an appropriate window size for converged statistics when sampling the data, which is necessary for detection of steady-state, and (c) detect steady-state within a particular tolerance. The last step, the projection test, is actually an oscillating-slope projection test, and is performed on j sequential data windows (i.e., more than two adjacent windows). It requires more than simply being within the user defined tolerance: the oscillating-slope projection test includes a condition that the slope must oscillate around zero when 2, which is an additional indication of steady-state. When all three tests are passed, the F-t-Pj test is passed, indicating that the prerequisites of steady-state detection have been met and also that conditions consistent with the definition of steady-state have been realized. This algorithm is applied to a variety of data sets that correspond to the diverse type of data trends that can be produced by kinetic Monte Carlo simulations. The algorithm is shown to be robust in its ability to handle differing functional forms, and is able to detect steady-state with low computational cost. Finally, the low computational cost of this method and its robustness towards varied data trends make it suitable for on-the-fly use in kinetic Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗

Analysis of time-dependent strain heterogeneity in shales using X-ray microscopy and digital volume correlation

In this paper, we report the results of a 4-dimensional (3 spatial dimensions plus time) in-situ creep experiment on a milli-core shale plug from Barnett formation in Texas. An outstanding challenge in shale mechanics is connecting observable characteristics to predictable mechanical and hydraulic behavior. We employ novel imaging and analysis techniques to statistically link sample mineralogy and structure to observed deformation. In particular, we focus on mechanisms of porosity reduction during transient response to step changes in stress, analogous to the aseismic deformation that occurs after an episode of fracturing. Our main goal in this study is to take the first steps to visualize and capture the small-scale deformation mechanisms involved in creep of shales, using imaging techniques. To do so, we conducted a uniaxial creep experiment inside an Xradia Versa 520 micro-CT system for 12 h at room temperature under a constant uniaxial stress of 25 MPa. We employed a Digital Volume Correlation (DVC) method to measure micro-scale deforspmation of the sample by comparing the resulting volumetric images with a high-resolution scan of the sample obtained prior to the creep experiment. The results indicate that the time-dependent deformation in the compliant clay plus kerogen matrix constituents facilitates rigid grain rotation and compaction of intergranular porosity at the boundary between rigid grains and matrix.

58 GEOSCIENCES↗

Probing late-stage hadronic interactions at high baryon density via $K^{⁎0}$ production in the RHIC beam energy scan program

A precision measurement of the $K^{⁎0}$ meson yield is reported in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.6, 19.6, and 27 GeV using the high-statistics data sample collected by the STAR experiment during the Beam Energy Scan II (BES-II) program at RHIC. The transeverse momentum (p T )-integrated yield ratios $\large{(}K^{⁎0} + \overline{K^{⁎0}}\large{)}/(K^+ + K^{-})$ in central collisions show a suppression relative to peripheral collisions at the (1.7–3.6) σ level, while a thermal model without final-stage rescattering overpredicts this ratio with a deviation of (6.9–8.2) σ. These results indicate a loss of the measured $K^{⁎0}$ signal in central collisions due to re-scattering of its hadronic decay products in the hadronic phase. The p T -integrated yield of charged kaons exhibits an approximate scaling with charged-particle multiplicity, independent of collision energy and system size. A similar trend is observed for the short-lived $K^{⁎0}$ resonance, although significant deviations emerge at lower energies. At BES energies, the $K^{⁎0}/K$ ratio shows stronger suppression than at the highest RHIC and LHC energies within a given multiplicity bin, particularly in central and mid-central collisions. This behavior is consistent with changes in the effective hadronic interaction cross section and is supported by transport model calculations, which indicate dominant meson–baryon interactions at lower energies and meson–meson interactions at higher energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Water freezing and ice melting

The generalized replica exchange method (gREM) is designed to sample states with coexisting phases and thereby to describe strong first order phase transitions. The isobaric MD version of the gREM is presented and applied to freezing of liquid water, and melting of hexagonal and cubic ice. It is confirmed that coexisting states are well sampled. The statistical temperature as a function of enthalpy, T S (H), is obtained. Hysteresis between freezing and melting is observed and discussed. The entropic analysis of phase transitions is applied and equilibrium transition temperatures, latent heats, and surface tensions are obtained for hexagonal ice↔liquid and cubic ice↔liquid, with excellent agreement with published values. A new method is given to assign water molecules among various symmetry types. As a result, pathways for water freezing, ultimately leading to hexagonal ice, are found to contain intermediate layered structures built from hexagonal and cubic ice.

36 MATERIALS SCIENCE↗

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗

Hidden depths in the local Universe: The Stellar Stream Legacy Survey

Context. Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-cold dark matter cosmology. These events are thought to provide important observational diagnostics of non-linear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistically meaningful sample of tidal streams beyond the Local Group has proven a daunting observational challenge, and the full potential for deepening our understanding of galaxy assembly using stellar streams has yet to be realised. Aims. Here we introduce the Stellar Stream Legacy Survey, a systematic imaging survey of tidal features associated with dwarf galaxy accretion around a sample of ~3100 nearby galaxies within z~0.02, including about 940 Milky Way analogues. Methods. Our survey exploits public deep imaging data from the DESI Legacy Imaging Surveys, which reach surface brightness as faint as ~29 mag arcsec –2 in the r band. As a proof of concept of our survey, we report the detection and broad-band photometry of 24 new stellar streams in the local Universe. Results. We discuss how these observations can yield new constraints on galaxy formation theory through comparison to mock observations from cosmological galaxy simulations. These tests will probe the present-day mass assembly rate of galaxies, the stellar populations and orbits of satellites, the growth of stellar halos, and the resilience of stellar disks to satellite bombardment.

79 ASTRONOMY AND ASTROPHYSICS↗

An X-ray-selected sample of candidate black holes in dwarf galaxies

We present a sample of hard X-ray-selected candidate black holes (BHs) in 19 dwarf galaxies. BH candidates are identified by cross-matching a parent sample of ∼44,000 local dwarf galaxies (M{sub ⋆}⩽3×10{sup 9} M{sub ⊙}, z<0.055) with the Chandra Source Catalog and subsequently analyzing the original X-ray data products for matched sources. Of the 19 dwarf galaxies in our sample, eight have X-ray detections reported here for the first time. We find a total of 43 point-like hard X-ray sources with individual luminosities L{sub 2−10 keV}∼10{sup 37}10{sup 40} erg s{sup −1}. Hard X-ray luminosities in this range can be attained by stellar-mass X-ray binaries (XRBs) and by massive BHs accreting at low Eddington ratio. We place an upper limit of 53% (10/19) on the fraction of galaxies in our sample hosting a detectable hard X-ray source consistent with the optical nucleus, although the galaxy center is poorly defined in many of our objects. We also find that 42% (8/19) of the galaxies in our sample exhibit statistically significant enhanced hard X-ray emission relative to the expected galaxy-wide contribution from low-mass and high-mass XRBs, based on the L{sub 2−10 keV}{sup XRB}−M{sub ⋆}−star formation rate relation defined by more massive and luminous systems. For the majority of these X-ray-enhanced dwarf galaxies, the excess emission is consistent with (but not necessarily due to) a nuclear X-ray source. Follow-up observations are necessary to distinguish between stellar-mass XRBs and active galactic nuclei powered by more massive BHs. In any case, our results support the notion that X-ray-emitting BHs in low-mass dwarf galaxies may have had an appreciable impact on reionization in the early universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy cluster profiles: a Gaussian mixture model approach to halo miscentering

Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center defined as the minimum in the potential well. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys. Our Python package for the Gaussian mixture model is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model.

Bayesian reasoning↗

First astrometric constraints on parity-violation in the gravitational wave background

Astrometry, the precise measurement of stellar positions and velocities, offers a promising approach to probing the low-frequency stochastic gravitational wave background (SGWB). Notably, astrometric vector sky maps are sensitive to parity-violating SGWB signals, which cannot be distinguished using pulsar timing array observations in an isotropic SGWB. We present the first astrometric constraints on parity-violating SGWB using quasar catalogs from Gaia DR3 and VLBA data. By analyzing the EB correlation in the two-point correlation function of the proper motions of the quasars, we find 2σ constraints on the parity-violating SGWB amplitude h 70 2 Ω V = -0.020 ± 0.025 from Gaia DR3 and h 70 2 Ω V = -0.004 ± 0.010 from VLBA. These constraints are valid in the frequency range 4.2 × 10 -18 Hz < f < 1.1 × 10 -8 Hz. Although not currently a tight constraint on theoretical models, this first attempt lays the groundwork for future investigations using more precise astrometric data.

Gravitational waves in GR and beyond: theory↗

Quantifying bias due to non-Gaussian foregrounds in an optimal reconstruction of CMB lensing and temperature power spectra

We estimate the magnitude of the bias due to non-Gaussian extragalactic foregrounds on the optimal reconstruction of the cosmic microwave background (CMB) lensing potential and temperature power spectra. The reconstruction is performed using a Bayesian inference method known as the marginal unbiased score expansion (MUSE). We apply MUSE to a minimum variance combination of multifrequency maps drawn from the Agora publicly available simulations of the lensed CMB and correlated extragalactic foreground emission. Taking noise levels appropriate to the SPT-3G D1 release, we find non-Gaussian foregrounds may bias the MUSE reconstruction of the lensing potential amplitude at the level of (0.7 ± 0.3)σ when using modes up to ℓ max = 3500. We do not detect a statistically significant bias, finding a value of (-0.4 ± 0.3)σ, when restricted to lower angular multipoles, ℓ max = 3000. This work is a first step toward understanding the impact of extragalactic foregrounds on optimal reconstructions of CMB temperature and lensing potential power spectra.

Statistical sampling techniques↗