Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clustering 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 325 records · Page 18

Fuzzy Classification of Ocean Color Satellite Data for Bio-optical Algorithm Constituent Retrievals

The ocean has been traditionally viewed as a 2 class system. Morel and Prieur (1977) classified ocean water according to the dominant absorbent particle suspended in the water column. Case 1 is described as having a high concentration of phytoplankton (and detritus) relative to other particles. Conversely, case 2 is described as having inorganic particles such as suspended sediments in high concentrations. Little work has gone into the problem of mixing bio-optical models for these different water types. An approach is put forth here to blend bio-optical algorithms based on a fuzzy classification scheme. This scheme involves two procedures. First, a clustering procedure identifies classes and builds class statistics from in-situ optical measurements. Next, a classification procedure assigns satellite pixels partial memberships to these classes based on their ocean color reflectance signature. These membership assignments can be used as the basis for a weighting retrievals from class-specific bio-optical algorithms. This technique is demonstrated with in-situ optical measurements and an image from the SeaWiFS ocean color satellite.

Campbell, Janet W.↗

Mg II Spectral Atlas and Flux Catalog for Late-Type Stars in the Hyades Cluster

In the course of a long-running IUE Guest Observer program, UV spectral images were obtained for more than 60 late-type members of the Hyades Cluster in order to investigate their chromospheric emissions. The emission line fluxes extracted from those observations were used to study the dependence of stellar dynamo activity upon age and rotation (IUE Observations of Rapidly Rotating Low-Mass Stars in Young Clusters: The Relation between Chromospheric Activity and Rotation). However, the details of those measurements, including a tabulation of the line fluxes, were never published. The purpose of the investigation summarized here was to extract all of the existing Hyades long-wavelength Mg II spectra in the IUE public archives in order to survey UV chromospheric emission in the cluster, thereby providing a consistent dataset for statistical and correlative studies of the relationship between stellar dynamo activity, rotation, and age over a broad range in mass.

Simon, Theodore↗

A Combined Al-Mg/Pb-Pb Age of the Solar System

Astrophysical models of planet formation and protoplanetary disk evolution demand precise and accurate timing of the sequence of events in the solar nebula, relative to a time t=0, usually taken to be during the short epoch of CAI (Ca-rich, Al-rich inclusion) formation. Most CAIs formed withlive26Al (mean-life τ26= 1.034 Myr [1]), with an abundance 26Al/27Al ≈ (26Al/27Al)SS= 5.23 × 10-5[2]. We adopt this as the widespread level of 26Al in the solar nebula at t=0. Assuming spatial homogeneity of 26Al, an inclusion that had less 26Al, (26Al/27Al)0, formed a time Δt26= τ26ln[(26Al/27Al)SS/ (26Al/27Al)0] after t=0.These ages are typical precise to within ±0.1 Myr. Igneous bulk meteorites and inclusions can be relatively dated by the Al-Mg chronometer, but only ifΔt26<6 Myr. The Pb-Pb system is useful as a longer relative chronometer. It yields absolute ages tPb using 207Pb/206Pb, 206Pb/204Pb, and 238U/235U ratios measured indifferent portions of a sample, assuming certain half-lives [4]. These absolute ages are uncertain to within ±9 Myr due to uncertainties in the 235U half-life[3], but times of formation ΔtPb= tCAI–tPb relative to t=0, are more precise(±0.5Myr),iftCAI can be found. Here, tCAI means the Pb-Pb age that would be measured in CAIs using the half-lives the community typically uses, if they achieved isotopic closure at t=0. Unfortunately, direct Pb-Pb dating of CAIs has not definitively determined tCAI. Based on four CAIs with canonical (26Al/27Al)0,[5,6] found tPb= 4567.30 ± 0.16 Myr. No other CAI ages with measured 238U/235U have been reported in the refereed literature, but there are hints of other CAIs with ages tPb= 4568.0 ± 0.2 Myr [7] and tPb= 4568.3 ± 0.2 Myr [8].It is unclear whether anyof these igneous type B CAIs isotopically closed at t=0 or represents tCAI. Instead of measurements, we advocate finding tCAI by minimizing the discrepancies between the Al-Mg and Pb-Pb chronometers. Assuming Δt26=ΔtPb, we find the implied t’CAI= tPb+Δt26, then define t*CAIas the weighted mean of the t’CAI. t*CAIis the best guess for the Pb-Pb age of t=0; the assumption of homogeneity is justified if the t’CAI cluster within errors around t*CAI. This statistical approach is similar to, but improves on, that of[9]. We find t*CAI= 4568.73 ± 0.16 Myr. Below we discuss our methodology and the implications of this age for CAIs, 1.4 Myr older than the reported and typically used age 4567.30±0.16 Myr. Methods: We base our estimate of t*CAIon five achondrites for which published (26Al/27Al)0and Pb-Pb ages exist: the quenched angrites D’Orbigny, Sahara 99555 (SAH 99555), and Northwest Africa (NWA) 1670; the pseudo-eucrite Asuka 881394; and the inner disk achondrite. All are “NC” (non-carbonaceous) achondrites that likely cooled quickly enough that the Al-Mg and Pb-Pb systems achieved isotopic closure simultaneously. We also considered the “CC” (carbonaceous chondrite-like) achondrites NWA 2796 and NWA 6704, butdo not include them in our fit. Al-Mg and Pb-Pb seem not to have closed simultaneously, possibly because formation in the outer disk from volatile-rich composition led to slower cooling. Of the 8 chondrules from NWA 5697 measured by [20,21], we also consider the 4 for which 238U/235U was measured: 2-C1, 5-C2, 3-C5, 11-C1.Depending on their post-formation thermal histories, the Al-Mg and Pb-Pb systems in chondrules may or may not have closed simultaneously. Table 1: (26Al/27Al)0, Pb-Pb ages of selected samplesSample(26Al/27Al)0/ 10-6RefPb-PbRefD’Orbigny3.98±0.15104563.43±0.19♮10-12SAH 995553.64±0.18104563.88±0.2712NWA16705.92±0.59104564.39±0.24*10Asuka 88139413.1±0.5613-154564.98±0.1715NWA 73253.03±0.14164563.4±2.616NWA 27963.94±0.16174562.89±0.5917NWA 67043.03±0.14184562.76±0.26192-C17.56±1.53204567.57±0.56*215-C27.04±1.51204567.54±0.52*213-C58.85±1.83204566.20±0.63*2111-C15.55±1.84204565.84±0.72*21*regression based on one subset of data points ♮weighted mean of two datasets Pb-Pb ages are proportional to the intercept of the line formed by linear regression of 207Pb/206Pb vs. 206Pb/204Pb data from various washes, leachates and residues of acid dissolution of a sample. Because contamination by terrestrial or primordial Pb is pervasive, some fractions must be excluded from regressions to ensure a fit with acceptable mean squares weighted deviation (MSWD). Usually points are excluded based on low [Pb], or low 206Pb/204Pb ratio(low radiogenic component), with single outliers identified [11,12,15,16,17]. In the starred examples (Table 1)and the case of 3 CAI Pb-Pb ages [5], up to half the points were excluded solely because did not fit a pre-determined line. This approach is vulnerable to confirmation bias and produces fits with low MSWD and too-low Pb-Pb age uncertainty. Regressing the same data points as [10], were produce the Pb-Pb age of NWA 1670 of 4564.39±0.24 Myr. But selecting other combinations of data points, other, equally valid, isochrons yield ages from 4563.77±0.21 Myr to 4564.64±0.23 Myr. Similar arguments apply to the Pb-Pb isochrons built by [21] for chondrules 2-C1 (we find 4567.33±0.44 to 4567.85±0.46 Myr), 5-C2 (4566.84±0.53 to 4567.70±0.44 Myr), 3-C5 (4565.84±0.54to 4567.04±0.54)and 11-C1 (4565.36±0.51 to 4565.74±0.45 Myr). Our adopted ages for these and NWA 1670 are listed in Table 2.Table 2. tCAI estimated from various components, using our regressions for the chondrules & NWA 1670.SampleΔt26(Myr)tPb(Myr)t’CAI(Myr)D’Orbigny5.05±0.044563.43±0.194568.48±0.19SAH 995555.14±0.054563.88±0.274569.02±0.27NWA16704.64±0.104564.21±0.634568.85±0.67Asuka 8813943.81±0.044564.98±0.174568.79±0.17NWA 73255.33±0.054563.4±2.64568.7±2.6NWA 27965.06±0.044562.89±0.594567.95±0.59NWA 67045.29±0.134562.76±0.264568.05±0.292-C12.00±0.214567.59±0.704569.59±0.725-C22.07±0.224567.23±0.914569.30±0.933-C51.84±0.214566.44±1.124568.28±1.1411-C12.32±0.344565.52±0.664567.84±0.73achondrite4568.72±0.16chondrules4568.76±0.58combined4568.73±0.16A weighted average of the five NC achondrites(or just D’Orbigny, SAH 99555 and Asuka 881394)yields t*CAI= 4568.72 ± 0.16Myr. All are consistent with this value to within 1.8σ, and MSWD=1.5. Including the 4 U-corrected chondrules, t*CAI= 4568.73± 0.16Myrwith MSWD=1.66, which is statistically significant. All chondrules and NC achondrites are consistent with this to within 1.8σ, (Figure 1).Figure 1. Al-Mg formation times after t=0 vs. Pb-Pb ages. The five NC achondrites and four chondrules are consistent with a Pb-Pb age of t=0 of 4568.7 Myr. Discussion: The data from achondrites and chondrules are consistent with a single Pb-Pb age at t=0, justifying the assumption of 26Al homogeneity. The age, 4568.7 Myr, is ≈1.4 Myr older than the commonly accepted Pb-Pb age of CAIs that formed with canonical 26Al/27Al at t=0 [3]. Others have interpreted the discrepancy to signify 26Al heterogeneity in the CAI-forming region[5,21]. We suggest instead that CAIs were exposed to transient heating events that reset the Pb-Pb system without disturbing the Al-Mg system. Notably, chondrules typically experienced transient heating at these times in the nebula [22]. If so, direct measurements of CAIs will not yield as reliable a Pb-Pb age of t=0 as statistical approaches like this and that of [9].References:[1] Auer et al. 2009. [2] Jacobsen, B et al. 2008, EPSL 272, 353-364. [3] Tissot, Fet al. 2017, GCA 213, 593-617.[4] Villa, I et al. 2016, GCA 172, 387-392.[5] Amelin, Y et al. 2010, EPSL 300, 343-350.[6] Connelly, Jet al. 2012, Science 338, 651.[7] Bouvier, Aet al. 2011, LPICo 1639, 9054. [8] Bouvier, Aand Wadhwa, M2010, Nat Geosci 3, 637-641.[9] Nyquist, Let al. 2009, GCA 73, 5115-5136.[10] Schiller et al. 2015. [11] Wadhwa & Brennecka 2012. [12] Tissot et al. 2017. [13] Nyquist et al. 2003. [14] Wadhwa et al. 2009. [15] Wimpenny et al. 2019, GCA 244, 478-501.[16] Koefoed et al. 2016, GCA 183, 31-45. [17] Bouvier, A et al. 2011, GCA 75, 5310-5323. [18] Sanborn, Met al. 2019, GCA 245, 577-596. [19] Amelin, Yet al. 2019, GCA 245, 628-642.[20] Bollard, Jet al. 2017, Sci Adv 3 ,e1700407. [21] Bollard, Jet al. 2019, GCA 260, 62-83.[22] Villeneuve, J et al. 2009, Science 325, 985

S. J. Desch↗

Statistical Mechanics on Lattices (smol) v0.0.1

smol is a Python implementation of the Cluster Expansion Method and extensions of this methodology; which is used to fit applied lattice models from first principle calculations for subsequent use in Monte Carlo simulations for computing thermodynamic and statistical properties of atomic configuration. smol has a several notable advantages over other similar software packages. The first being its Python implementation which makes it easy to use even for users with little software and coding skills. Despite being implemented in Python, essential Monte Carlo routines are implemented in Cython such that performance is not compromised and is competitive even to similar software implemented in C/C++. Furthermore, the package has specific functionality for handling complex ionic materials (such as cathodes and electrolytes) that is not readily available in other packages. Lastly, the software has a flexible and modular design with the intention of making it fast and efficient to develop and extend the methodology.

Barroso-Luque, Luis↗

Low-level processing for real-time image analysis

A system that detects object outlines in television images in real time is described. A high-speed pipeline processor transforms the raw image into an edge map and a microprocessor, which is integrated into the system, clusters the edges, and represents them as chain codes. Image statistics, useful for higher level tasks such as pattern recognition, are computed by the microprocessor. Peak intensity and peak gradient values are extracted within a programmable window and are used for iris and focus control. The algorithms implemented in hardware and the pipeline processor architecture are described. The strategy for partitioning functions in the pipeline was chosen to make the implementation modular. The microprocessor interface allows flexible and adaptive control of the feature extraction process. The software algorithms for clustering edge segments, creating chain codes, and computing image statistics are also discussed. A strategy for real time image analysis that uses this system is given.

Eskenazi, R.↗

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗

Quantitative assessment of particle dispersion in polymeric composites and its effect on mechanical properties

In this work, an automated image analysis tool is developed to establish quantitative correlations between the particle/cluster size distribution and the mechanical properties of particle reinforced polymer composites (PRPC).This automated image analysis tool is developed within python programming software to process and analyze the microstructural images of the polymerbased composite materials. The spent coffee bean powder (SCBP) reinforced poly-propylene carbonate (PPC) polymer composite with differing wt.% of the filler is selected for the analysis. Detailed statistical analysis of the microstructural images reveals that ‘clustering of clusters’ is also presented in addition to the most commonly reported ‘clustering of particles’, and the distribution of particle/clusters is bimodal. Based on these findings, an effective volume fraction for the filler material is proposed to mainly capture the agglomeration effect. With this effective volume fraction, the standard rule-of-mixture model correctly captures the experimentally measured tensile strength and modulus as a function of filler wt.%. Further, the applicability of this effective volume fraction for other theoretical models is also analyzed. The detailed statistical analysis of the microstructure and the proposed effective volume fraction helps to develop a deeper quantitative understanding of the PRPC than the conventional qualitative correlation of microstructural features with the properties and failure processes.

36 MATERIALS SCIENCE↗

The Second Most Distant Cluster of Galaxies in the Extended Medium Sensitivity Survey

We report on our ASCA, Keck, and ROSAT observations of MS 1137.5+6625, the second most distant cluster of galaxies in the Einstein Extended Medium Sensitivity Survey (EMSS), at redshift 0.78. We now have a full set of X-ray temperatures, optical velocity dispersions, and X-ray images for a complete, high-redshift sample of clusters of galaxies drawn from the EMSS. Our ASCA observations of MS 1137.5 +6625 yield a temperature of 5.7 (+2.1)(-1.1) keV and a metallicity of 0.43 (+40)(-3.7) solar, with 90% confidence limits. Keck II spectroscopy of 22 cluster members reveals a velocity dispersion of 884 (+185)(-124) km 24/s. This cluster is the most distant in the sample with a detected iron line. We also derive a mean abundance at z = 0.8 by simultaneously fitting X-ray data for the two z = 0.8 clusters, and obtain an abundance of Z(sub Fe) = 0.33 (+.26)(-.23). Our ROSAT observations show that MS 1137.5+6625 is regular and highly centrally concentrated. Fitting of a Beta model to the X-ray surface brightness yields a core radius of only 71/h kpc (q(sub o) = 0.1) with Beta = 0.70(+.45)(-.15) The gas mass interior to 0.5/h Mpc is thus 1.2 (+0.2)(-0.3) X 10(exp 13) h(exp - 5/2) Solar Mass (q(sub o) = 0.1). If the cluster's gas is nearly isothermal and in hydrostatic equilibrium with the cluster potential, the total mass of the cluster within this same region is 2.1(+1.5)(-0.8) X 10exp 14)/h Solar Mass, giving a gas fraction of 0.06 +/-0.04 h (exp -3/2). This cluster is the highest redshift EMSS cluster showing evidence for a possible cooling flow (about 20-400 Solar Mass/yr). The velocity dispersion, temperature, gas fraction, and iron abundance of MS 1137.5+6625 are all statistically the same as those properties in lower red- shift clusters of similar luminosity. With this cluster's temperature now in hand, we derive a high-redshift temperature function for EMSS clusters at 0.5 < z < 0.9 and compare it with temperature functions at lower redshifts, showing that the evolution of the temperature function is relatively modest. Supplementing our high-redshift sample with other data from the literature, we demonstrate that neither the cluster luminosity-temperature relation, nor cluster metallicities, nor the cluster gas evolved with redshift. The very modest degree of evolution in the luminosity-temperature relation inferred from these data is inconsistent with the absence of evolution in the X-ray luminosity functions derived from ROSAT cluster surveys if a critical density structure formation model is assumed.

Donahue, Megan↗

The DESI N -body simulation project – I. Testing the robustness of simulations for the DESI dark time survey

ABSTRACT Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogues to validate data analysis pipelines and identify potential systematics. We compare three N-body simulation codes, abacus, gadget-2, and swift, to investigate the regimes in which their results agree. We run N-body simulations at three different mass resolutions, 6.25 × 108, 2.11 × 109, and 5.00 × 109 h−1 M⊙, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the Dark Energy Spectroscopic Instrument (DESI) statistical error for $s\ \gt\ 20\ h^{-1}$ Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of 2.1 × 109 h−1 M⊙ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $20\ h^{-1}$ Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1 per cent for k ≤ 10 h Mpc−1. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-N-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy–halo relation will be presented in a follow-up study.

79 ASTRONOMY AND ASTROPHYSICS↗

Very low mass stars and white dwarfs in NGC 6397

Deep Wide Field/Planetary Camera 2 (WFPC2) images in wide bands centered at 606 and 802 nm were taken with the Hubble Space Telescope (HST) 4.6 min from the center of the galactic globular cluster NGC 6397. The images were used to accurately position approximately 2120 stars detected in the field on a color magnitude diagram down to a limiting magnitude m(sub 814) approximately = m(sub I) approximately = 26 determined reliably and solely by counting statistics. A white dwarf sequence and a rich, narrow cluster main sequence are detected for the first time, the latter stretching from m(sub 814) = 18.5 to m(sub 814) = 24.0 where it becomes indistinguishable from the field population. Two changes of slope of the main sequence at m(sub 814) approximately = 20 and m(sub 814) approximately = 22.5 are evident. The corresponding luminosity function increases slowly from M(sub 814) approximately = 6.5 to 8.5 are expected from ground-based observations but then drops sharply from there dwon to the measurement limit. The corresponding mass function obtained bu using the only presently available mass-luminosity function for the cluster's metallicity rises to a plateau between approximately 0.25 and approximately 0.15 solar mass, but drops toward the expected mass limit of the normal hydrogen burning main sequence at approximately 0.1 solar mass. This result is in clear contrast to that obtained from the ground and implies either a substantial modification of the cluster's initial mass function due to dynamical evolution in its lifetime, or that very low mass stars are not produced in any dynamically significant amount by clusters of this type. The white dwarf sequence is in reasonable agreement with a cooling sequence of models of mass 0.5 solar mass at the canonical distance of NGC 6397 with a scatter that is most likely due to photometric errors, but may also reflect real differences in mass or chemical composition. Contamination from unresolved galaxies, which cannot be reliably identified with our filters, makes it difficult to meaningfully compare the observed white dwarf luminosity function with its theoretical counterpart.

Paresce, Francesco↗

Constraining gravity with a new precision 𝐸 𝐺 estimator using Planck + SDSS BOSS data

The 𝐸 𝐺 statistic is a discriminating probe of gravity developed to test the prediction of general relativity (GR) for the relation between gravitational potential and clustering on the largest scales in the observable Universe. We present a novel high-precision estimator for the 𝐸 𝐺 statistic using CMB lensing and galaxy clustering correlations that carefully matches the effective redshifts across the different measurement components to minimize corrections. A suite of detailed tests is performed to characterize the estimator’s accuracy, its sensitivity to assumptions and analysis choices, and the non-Gaussianity of the estimator’s uncertainty is characterized. After finalization of the estimator, it is applied to Planck CMB lensing and SDSS CMASS and LOWZ galaxy data. We report the first harmonic space measurement of 𝐸 𝐺 using the LOWZ sample and CMB lensing and also updated constraints using the final CMASS sample and the latest Planck CMB lensing map. We find $\hat{𝐸}$$^{Planck+CMASS}_{𝐺}$ = 0.3⁢6$^{+0.06}_{−0.05}$⁢(68.27%) and $\hat{𝐸}$$^{Planck+LOWZ}_{𝐺}$ = 0.4⁢0$^{+0.11}_{−0.09}$⁢(68.27%), with additional subdominant systematic error budget estimates of 2% and 3%, respectively. Using Ω m,0 constraints from Planck and SDSS BAO observations, Λ⁢CDM-GR predicts 𝐸$^{GR}_ {𝐺}$⁡(𝑧 =0.555) = 0.401 ± 0.005 and 𝐸$^{GR}_{𝐺}$⁡(𝑧 =0.316) = 0.452 ± 0.005 at the effective redshifts of the CMASS and LOWZ based measurements. We report the measurement to be in good statistical agreement with the Λ⁢CDM-GR prediction and report that the measurement is also consistent with the more general GR prediction of scale independence for 𝐸 𝐺 . Furthermore, this work provides a carefully constructed and calibrated statistic with which 𝐸 𝐺 measurements can be confidently and accurately obtained with upcoming survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

New output improvements for CLASSY

Additional output data and formats for the CLASSY clustering algorithm were developed. Four such aids to the CLASSY user are described. These are: (1) statistical measures; (2) special map types; (3) formats for standard output; and (4) special cluster display method.

Rassbach, M. E.↗

Constraints on the Energy Content of the Universe from a Combination of Galaxy Cluster Observables

We demonstrate that constraints on cosmological parameters from the distribution of clusters as a function of redshift (dN/dz) are complementary to accurate angular diameter distance (D(sub A)) measurements to clusters, and their combination significantly tightens constraints on the energy density content of the Universe. The number counts can be obtained from X-ray and/or SZ (Sunyaev-Ze'dovich effect) surveys, and the angular diameter distances can be determined from deep observations of the intra-cluster gas using their thermal bremsstrahlung X-ray emission and the SZ effect. We combine constraints from simulated cluster number counts expected from a 12 deg(sup 2) SZ cluster survey and constraints from simulated angular diameter distance measurements based on the X-ray/SZ method assuming a statistical accuracy of 10% in the angular diameter distance determination of 100 clusters with redshifts less than 1.5. We find that Omega(sub m), can be determined within about 25%, Omega(sub lambda) within 20% and w within 16%. We show that combined dN/dz+(sub lambda) constraints can be used to constrain the different energy densities in the Universe even in the presence of a few percent redshift dependent systematic error in D(sub lambda). We also address the question of how best to select clusters of galaxies for accurate diameter distance determinations. We show that the joint dN/dz+ D(lambda) constraints on cosmological parameters for a fixed target accuracy in the energy density parameters are optimized by selecting clusters with redshift upper cut-offs in the range 0.55 approx. less than 1. Subject headings: cosmological parameters - cosmology: theory - galaxies:clusters: general

Molnar, Sandor M.↗

Clustering and voids in the Lyman-alpha forest

In a statistical study of quasar absorption lines, no evidence is found for large voids in the Ly-alpha forest, contrary to some recent suggestions. The 28 Mpc gap found by Crotts (1987) in the spectrum of Q0420-388 is not statistically significant, and it is possible, from a sample of 18 quasar spectra, to reject at the 2 sigma level the hypothesis that comoving voids of 30 Mpc fill as much as 15 percent of the universe. However, a significant (at the 99 percent level) excess of small line intervals of velocity width 200-600 km/s is found. This provides evidence either for the gravitationally induced correlations proposed in some models or for the fragmented shells proposed in others.

Ostriker, Jeremiah P.↗

On the Determination of Poisson Statistics for Haystack Radar Observations of Orbital Debris

A convenient and powerful method is used to determine if radar detections of orbital debris are observed according to Poisson statistics. This is done by analyzing the time interval between detection events. For Poisson statistics, the probability distribution of the time interval between events is shown to be an exponential distribution. This distribution is a special case of the Erlang distribution that is used in estimating traffic loads on telecommunication networks. Poisson statistics form the basis of many orbital debris models but the statistical basis of these models has not been clearly demonstrated empirically until now. Interestingly, during the fiscal year 2003 observations with the Haystack radar in a fixed staring mode, there are no statistically significant deviations observed from that expected with Poisson statistics, either independent or dependent of altitude or inclination. One would potentially expect some significant clustering of events in time as a result of satellite breakups, but the presence of Poisson statistics indicates that such debris disperse rapidly with respect to Haystack's very narrow radar beam. An exception to Poisson statistics is observed in the months following the intentional breakup of the Fengyun satellite in January 2007.

Stokely, Christopher L.↗

Bayesian prior construction for uncertainty quantification in first-principles statistical mechanics

First-principles statistical mechanics enables the prediction of thermodynamic and kinetic properties of materials, but is computationally expensive. Many approaches require surrogate models to calculate energies within Monte Carlo or molecular dynamics simulations. Inexpensive surrogates such as cluster expansions enable otherwise intractable calculations by interpolating data from higher accuracy methods, such as Density Functional Theory (DFT). Surrogate models introduce uncertainty into downstream calculations, in addition to any uncertainty inherent to DFT calculations. Bayesian frameworks address this by quantifying uncertainty and incorporating expert knowledge through priors. However, constructing effective priors remains challenging. This work introduces and describes practical strategies for building Bayesian cluster expansions, focusing on basis truncation, hyperparameter selection, and ground state replication. We analyze multiple basis truncation schemes, compare cross-validation to the evidence-approximation for hyperparameter optimization, and provide methods to find and enforce ground-state-preserving models through priors. Additionally, we compare the uncertainties between different approximations to DFT (LDA, PBE, SCAN) against the uncertainty introduced with the use of cluster expansion surrogate models. These approaches are demonstrated on the BCC Li x Mg 1-x and Li x Al 1-x alloys, which are both of interest for solid-state Li batteries. Our results provide guidelines for constructing and utilizing Bayesian cluster expansions, thereby improving the transparency of materials modeling. Furthermore, the approaches and insights developed in this work can be transferred to a wide range of cluster expansion surrogate models, including the atomic cluster expansion and related machine-learned interatomic potential architectures.

Alloy theory↗

Using factor analysis to identify neuromuscular synergies during treadmill walking

Neuroscientists are often interested in grouping variables to facilitate understanding of a particular phenomenon. Factor analysis is a powerful statistical technique that groups variables into conceptually meaningful clusters, but remains underutilized by neuroscience researchers presumably due to its complicated concepts and procedures. This paper illustrates an application of factor analysis to identify coordinated patterns of whole-body muscle activation during treadmill walking. Ten male subjects walked on a treadmill (6.4 km/h) for 20 s during which surface electromyographic (EMG) activity was obtained from the left side sternocleidomastoid, neck extensors, erector spinae, and right side biceps femoris, rectus femoris, tibialis anterior, and medial gastrocnemius. Factor analysis revealed 65% of the variance of seven muscles sampled aligned with two orthogonal factors, labeled 'transition control' and 'loading'. These two factors describe coordinated patterns of muscular activity across body segments that would not be evident by evaluating individual muscle patterns. The results show that factor analysis can be effectively used to explore relationships among muscle patterns across all body segments to increase understanding of the complex coordination necessary for smooth and efficient locomotion. We encourage neuroscientists to consider using factor analysis to identify coordinated patterns of neuromuscular activation that would be obscured using more traditional EMG analyses.

NASA Center JSC↗

H I Depletion Begins Well Beyond the Virial Radius: A FAST Stacking Study of 36 Galaxy Clusters to 5 × R 200

Abstract We present a stacking study of the neutral atomic hydrogen (H i ) content in and around 36 local galaxy clusters at z < 0.07, using a combination of the FAST All Sky H i survey (FASHI) and the extensive spectroscopic catalog mainly from the Dark Energy Spectroscopic Instrument (DESI). We employ spectral stacking techniques to probe the average H i mass and HI-to-stellar mass ratio ( M HI / M * ) for member galaxies down to stellar masses of M * ∼ 10 9 M ⊙ , spanning a projected cluster-centric distance of up to 5 R 200 . Our analysis reveals a pronounced environmental effect; both M HI and M HI / M * decrease steadily toward the cluster center, dropping by ∼0.5 dex on average from the outskirts to the core. Crucially, we find that M HI / M * of galaxies remain lower than the field galaxies even at the 5 R 200 . This provides direct, statistical evidence for substantial gas stripping and preprocessing in the cluster outskirts, likely occurring in infalling groups and large-scale filaments. By further splitting the sample by g − r color, we show that the H i deficiency persists at fixed galaxy color; even the bluest cluster members exhibit ∼0.5 dex lower M HI / M * than field galaxies of similar color, reflecting environmental effects on the cold gas reservoir prior to full optical transformation. The total H i mass within clusters and their outskirts agrees broadly with predictions from cosmological simulations. Our results underscore the critical role of the extended cluster environment in quenching galaxies by depleting their cold gas reservoirs well before they enter the dense cluster core.

Cheng, Cheng [Chinese Academy of Sciences South Am↗