Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster finding”

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,405 records · Page 78

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Here, our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of Ru in Improving the Activity of Pd toward Hydrogen Evolution and Oxidation Reactions in Alkaline Solutions

Improving the reaction kinetics of hydrogen evolution and oxidation reactions (HER/HOR) in alkaline media is critical to promote the development of alkaline fuel cells and electrolyzers. Here, we prepared Pd 3 Ru alloy nanocatalysts with Ru segregated on the surfaces, forming adatoms and clusters. This structure dramatically lowered the overpotential of Pd toward HER in 1 M KOH by 104 mV at 10 mA cm –2 . The HER activity was even higher than that of Pt (6 mV improvement at 10 mA cm –2 ). Theoretical simulation results revealed that Ru adatoms/clusters on the surface could weaken the hydrogen-binding energy and promote the OH adsorption, consequently lowering the reaction barrier of the rate-determining step in HER. Here, our findings are of significance for clarifying the role of Ru in bimetallic catalysts and rational design of more active catalysts for HER/HOR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The brown dwarf population in the star-forming region NGC 2264

ABSTRACT The brown dwarf population in the canonical star-forming region NGC 2264 is so far poorly explored. We present a deep, multiwavelength, multiepoch survey of the star-forming cluster NGC 2264, aimed to identify young brown dwarf candidates in this region. Using criteria including optical/near-infrared (IR) colours, variability, Spitzer mid-IR colour excess, extinction, and Gaia parallax and proper motion (in order of relevance), we select 902 faint red sources with indicators of youth. Within this sample, we identify 429 brown dwarf candidates based on their IR colours. The brown dwarf candidates are estimated to span a mass range from 0.01 to 0.08 M⊙. We find rotation periods for 44 sources, 15 of which are brown dwarf candidates, ranging from 3.6 h to 6.5 d. A subset of 38 brown dwarf candidates show high-level irregular variability indicative of ongoing disc accretion, similar to the behaviour of young stars.

Pearson, Samuel↗

The Pristine survey – XVI. The metallicity of 26 stellar streams around the Milky Way detected with the STREAMFINDER in Gaia EDR3

ABSTRACT We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by the application of the STREAMFINDER algorithm to the Gaia Early Data Release 3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, supporting the reliability of STREAMFINDER in finding streams and the power of Pristine to measure precise metallicities. We further present six candidate structures with coherent phase-space and metallicity signals that are very likely streams. The majority of studied streams are very metal-poor (14 structures with [Fe/H] < −2.0) and include three systems with [Fe/H] < −2.9 (C-11, C-19, and C-20). These streams could be the closest debris of low-luminosity dwarf galaxies or may have originated from globular clusters of significantly lower metallicity than any known current Milky Way globular cluster. Our study shows that the promise of the Gaia data for Galactic Archeology studies can be substantially strengthened by quality photometric metallicities, allowing us to peer back into the earliest epochs of the formation of our Galaxy and its stellar halo constituents.

79 ASTRONOMY AND ASTROPHYSICS↗

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

Non-smooth Bayesian optimization in tuning scientific applications

Tuning algorithmic parameters to optimize the performance of large, complicated computational codes is an important problem involving finding the optima and identifying regimes defined by non-smooth boundaries in black-box functions. Within the Bayesian optimization framework, the Gaussian process surrogate model produces smooth mean functions, but functions in the tuning problem are often non-smooth, which is exacerbated by the fact that we usually have limited sequential samples from the black-box function. Here, motivated by these issues encountered in tuning, we propose a novel Gaussian process model called a clustered Gaussian process (cGP), where the components are dynamically updated by clustering. In our studies, the performance of cGP can be better than stationary GPs in nearly 90% of the experiments and better than non-stationary GPs in nearly 70% of the repeated experiments while requiring less computational cost. cGP provides a novel approach for dynamic GP, computes more efficiently than recursive partitioning, and discovers non-smoothness regimes. We provide extensive experiments including high-performance computing (HPC) and industrial simulation functions to show the effectiveness of our methods.

97 MATHEMATICS AND COMPUTING↗

An analytic approximation to the covariance between pre- and post-reconstruction galaxy two-point statistics

We present a simple analytic approximation for the covariance between pre-reconstruction galaxy power spectrum measurements and post-reconstruction two-point correlation functions. This cross-covariance is essential for joint analyses that combine full-shape clustering information with baryon acoustic oscillation (BAO) measurements, as commonly performed in modern spectroscopic surveys. Our model builds on the disconnected contribution to the covariance and accounts for the damping of correlations due to the BAO reconstruction process. We validate our analytic prescription against numerical simulations from the Dark Energy Spectroscopic Instrument (DESI), testing both idealized cubic geometries and realistic survey configurations including complex footprints and fiber assignment effects. Despite neglecting survey window functions in the analytic calculation, we find excellent agreement with simulation-based covariances and demonstrate that cosmological parameter constraints are virtually unchanged when using our approximation. Our results show that the pre-post cross-covariance is sufficiently small that even approximate treatments are adequate for cosmological inference, opening a pathway toward fully analytic covariance matrices for next-generation galaxy surveys.

baryon acoustic oscillations↗

On the influence of halo mass accretion history on galaxy properties and assembly bias

Halo assembly bias is the secondary dependence of the clustering of dark matter haloes on their assembly histories at fixed halo mass. This established dependence is expected to manifest itself on galaxy clustering, a potential effect commonly known as galaxy assembly bias. Using the IllustrisTNG300 magnetohydrodynamical simulation, in this work we analyse the dependence of the properties and clustering of galaxies on the specific mass accretion history of their hosting haloes (sMAH). We first show that several halo and galaxy properties strongly correlate with the slope of the sMAH (β) at fixed halo mass. Haloes with increasingly steeper β increment their masses faster early on, and their hosted galaxies present larger stellar-to-halo mass ratios, lose their gas faster, reach the peak of their star formation histories at higher redshift, and become quenched earlier. We also demonstrate that β provides a more stable link to these key galaxy formation properties than other broadly employed halo proxies, such as formation time. Finally, we measure the secondary dependence of galaxy clustering on β at fixed halo mass. By tracing back the evolution of individual haloes, we show that the amplitude of the galaxy assembly bias signal for the progenitors of z = 0 galaxies increases with redshift, reaching a factor of 2 at z = 1 for haloes of Mhalo = 10 11.5 –10 12 h –1 M ⊙ . The measurement of the evolution of assembly bias along the merger tree provides a new theoretical perspective to the study of secondary bias. Our findings have also important implications for the generation of mock catalogues for upcoming cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Mass dependence of halo baryon fractions from the kinetic Sunyaev–Zeldovich effect

ABSTRACT We detect the kinetic Sunyaev–Zeldovich imprint of peculiar motions of galaxy groups and clusters, using the photometric Dark Energy Spectroscopic Instrument Legacy Surveys together with cosmic microwave background (CMB) maps from the Atacama Cosmology Telescope (ACT). We develop a comprehensive forward model based on the AbacusSummit cosmological simulations: Mock galaxy group catalogues and synthetic kinetic Sunyaev–Zeldovich (kSZ) maps are generated, together with a reconstructed peculiar velocity field that allows for photo-z errors, redshift-space distortions, and survey masks. We investigate possible contamination from the cosmic infrared background (CIB), finding that CIB effects are subdominant to the kSZ signal in the relevant ACT frequency channel. We then predict the kSZ signal expected when stacking CMB temperature maps around groups, taking account of their estimated radial velocities. Comparing the model with observations, we are able to constrain the total baryon fraction within haloes, as well as their internal gas profiles. We find evidence for mass dependence of the halo baryon fraction within the virial radius. The gas fraction in massive groups is consistent with the universal baryon fraction, but low-mass groups ($10^{12.5} \lesssim M\, /\, h^{-1}\, \mathrm{M}_\odot \lesssim 10^{14}$) are depleted to $0.38 \pm 0.11$ times the universal baryon fraction. We find this low virial baryon fraction to be consistent with an extended gas profile, for which the total baryon content reaches the universal value well beyond the virial radius. This conclusion is consistent with previous analyses using X-ray, kSZ, and weak lensing, and plausibly reflects energetic feedback processes from the galaxies in these haloes.

Roper, Finn A [Institute for Astronomy, University↗

Application of sulfur SAD to small crystals with a large asymmetric unit and anomalous substructure

The application of sulfur single-wavelength anomalous dispersion (S-SAD) to determine the crystal structures of macromolecules can be challenging if the asymmetric unit is large, the crystals are small, the size of the anomalously scattering sulfur structure is large and the resolution at which the anomalous signals can be accurately measured is modest. Here, as a study of such a case, approaches to the SAD phasing of orthorhombic Ric-8A crystals are described. The structure of Ric-8A was published with only a brief description of the phasing process [Zeng et al. (2019), Structure , 27 , 1137–1141]. Here, alternative approaches to determining the 40-atom sulfur substructure of the 103 kDa Ric-8A dimer that composes the asymmetric unit are explored. At the data-collection wavelength of 1.77 Å measured at the Frontier micro-focusing Macromolecular Crystallography (FMX) beamline at National Synchrotron Light Source II, the sulfur anomalous signal strength, |Δ ano |/σΔ ano ( d ′′/sig), approaches 1.4 at 3.4 Å resolution. The highly redundant, 11 000 000-reflection data set measured from 18 crystals was segmented into isomorphous clusters using BLEND in the CCP 4 program suite. Data sets within clusters or sets of clusters were scaled and merged using AIMLESS from CCP 4 or, alternatively, the phenix.scale_and_merge tool from the Phenix suite. The latter proved to be the more effective in extracting anomalous signals. The HySS tool in Phenix , SHELXC / D and PRASA as implemented in the CRANK 2 program suite were each employed to determine the sulfur substructure. All of these approaches were effective, although HySS , as a component of the phenix.autosol tool, required data from all crystals to find the positions of the sulfur atoms. Critical contributors in this case study to successful phase determination by SAD included (i) the high-flux FMX beamline, featuring helical-mode data collection and a helium-filled beam path, (ii) as recognized by many authors, a very highly redundant, multiple-crystal data set and (iii) the inclusion within that data set of data from crystals that were scanned over large ω ranges, yielding highly isomorphous and highly redundant intensity measurements.

59 BASIC BIOLOGICAL SCIENCES↗

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗

The impact of the fiducial cosmology assumption on BAO distance scale measurements

ABSTRACT Standard analysis pipelines for measurements of Baryon Acoustic Oscillations (BAO) in galaxy surveys make use of a fiducial cosmological model to guide the data compression required to transform from observed redshifts and angles to the measured angular and radial BAO peak positions. In order to remove any dependence on the fiducial cosmology from the results, all models compared to the data should mimic the compression and its dependence on the fiducial model. In practice, approximations are made when testing models: (1) There is assumed to be no residual dependence on the fiducial cosmology after reconstruction, (2) differences in the distance–redshift relationship are assumed to match a linear scaling, and (3) differences in clustering between true and fiducial models are assumed to be removed by the free parameters used to null the non-BAO signal. We test these approximations using the current standard measurement procedure with a set of halo catalogues from the aemulus suite of N-body simulations, which span a range of wCDM cosmological models. We focus on reconstruction of the primordial BAO and locating the BAO. For the range of wCDM cosmologies covered by the aemulus suite, we find no evidence for systematic errors in the measured BAO shift parameters α∥ and α⊥ to $\lt 0.1\%$. However, the measured errors $\sigma _{\alpha _{\parallel }}$ and $\sigma _{\alpha _{\bot }}$ show a notable absolute increase by up to +0.001 and +0.002, respectively, in the case that the fiducial cosmology does not match the truth. These effects on the inferred BAO scale will be important, given the precision of measurements expected from future surveys including DESI, Euclid, and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Energy Survey Year 1 results: constraints on intrinsic alignments and their colour dependence from galaxy clustering and weak lensing

Abstract We perform a joint analysis of intrinsic alignments and cosmology using tomographic weak lensing, galaxy clustering, and galaxy–galaxy lensing measurements from Year 1 (Y1) of the Dark Energy Survey. We define early- and late-type subsamples, which are found to pass a series of systematics tests, including for spurious photometric redshift error and point spread function correlations. We analyse these split data alongside the fiducial mixed Y1 sample using a range of intrinsic alignment models. In a fiducial non-linear alignment model analysis, assuming a flat Λ cold dark matter cosmology, we find a significant difference in intrinsic alignment amplitude, with early-type galaxies favouring $A_\mathrm{IA} = 2.38^{+0.32}_{-0.31}$ and late-type galaxies consistent with no intrinsic alignments at $0.05^{+0.10}_{-0.09}$. The analysis is repeated using a number of extended model spaces, including a physically motivated model that includes both tidal torquing and tidal alignment mechanisms. In multiprobe likelihood chains in which cosmology, intrinsic alignments in both galaxy samples and all other relevant systematics are varied simultaneously, we find the tidal alignment and tidal torquing parts of the intrinsic alignment signal have amplitudes $A_1 = 2.66 ^{+0.67}_{-0.66}$, $A_2=-2.94^{+1.94}_{-1.83}$, respectively, for early-type galaxies and $A_1 = 0.62 ^{+0.41}_{-0.41}$, $A_2 = -2.26^{+1.30}_{-1.16}$ for late-type galaxies. In the full (mixed) Y1 sample the best constraints are $A_1 = 0.70 ^{+0.41}_{-0.38}$, $A_2 = -1.36 ^{+1.08}_{-1.41}$. For all galaxy splits and IA models considered, we report cosmological parameter constraints consistent with the results of the main DES Y1 cosmic shear and multiprobe cosmology papers.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmic distance inference from purely geometric BAO methods: Linear point standard ruler and correlation function model fitting

Leveraging the baryon acoustic oscillations (BAO) feature present in clustering 2-point statistics, we aim to measure cosmological distances independently of the underlying background cosmological model. However this inference is complicated by late-time nonlinearities that introduce model and tracer dependencies in the clustering correlation function and power spectrum, which must be properly accounted for. With this in mind, we introduce the “purely geometric-BAO,” which provides a rigorous tool to measure cosmological distances without assuming a specific background cosmology. We focus on the 2-point clustering correlation function monopole, and show how to implement such an inference scheme employing two different methodologies: the linear point standard ruler (LP) and correlation-function model-fitting (CF-MF). For the first time we demonstrate how, by means of the CF-MF, we can measure very precisely the sound-horizon/isotropic-volume-distance ratio, r d /DV(z¯), while correctly propagating all the uncertainties. Using synthetic data, we compare the outcomes of the two methodologies, and find that the LP provides up to 50% more precise measurements than the CF-MF. Finally, we test a procedure widely employed in BAO analyses: fitting the 2-point function while fixing the cosmological and the non-linear-damping parameters at fiducial values. We find that this underestimates the distance errors by nearly a factor of 2. Furthermore, we recommend that this practice be reconsidered, whether for parameter determination or model selection.

79 ASTRONOMY AND ASTROPHYSICS↗

VAST: the Void Analysis Software Toolkit

Voids are expansive regions in the universe containing significantly fewer galaxies than surrounding galaxy clusters and filaments. They are a fundamental feature of the cosmic web and provide important information about galaxy physics and cosmology. For example, correlations between voids and luminous tracers of large-scale structure improve constraints on the expansion of the universe as compared to using tracers alone, and numerous studies have shown that the void environment influences the evolution of galaxies. However, what constitutes a void is vague and formulating a concrete definition to use in a void-finding algorithm is not trivial. As a result, several different algorithms exist to identify these cosmic underdensities. Our Void Analysis Software Toolkit, or VAST, provides Python 3 implementations of two such algorithms: VoidFinder and V 2 . This consolidation of two popular void-finding algorithms allows the user to, for example, easily compare the results of their analysis using different void definitions.

97 MATHEMATICS AND COMPUTING↗

Constraining the galaxy-halo connection of infrared-selected unWISE galaxies with galaxy clustering and galaxy-CMB lensing power spectra

We present the first detailed analysis of the connection between galaxies and their dark matter halos for the unWISE galaxy catalog—a full-sky, infrared-selected sample built from WISE data, containing over 500 million galaxies. Using unWISE galaxy-galaxy autocorrelation and Planck CMB lensing-galaxy cross-correlation measurements down to 10 arcmin angular scales, we constrain the halo occupation distribution (HOD), a model describing how central and satellite galaxies are distributed within dark matter halos, for three unWISE galaxy samples at mean redshifts z ¯ ≈ 0.6 , 1.1, and 1.5, assuming a fixed cosmology at the best-fit Planck Λ CDM values. We constrain the characteristic minimum halo mass to host a central galaxy, M min HOD = 1.8 3 - 1.63 + 0.41 × 10 12 M ⊙ / h , 5.22 - 4.80 + 0.34 × 10 12 M ⊙ / h , 6.60 - 1.11 + 0.30 × 10 13 M ⊙ / h and the mass scale at which one satellite galaxy per halo is found, M 1 ' = 1.1 3 - 0.70 + 0.32 × 10 13 M ⊙ / h , 1.18 - 1.11 + 0.30 × 10 13 M ⊙ / h , 1.23 - 1.17 + 0.14 × 10 14 M ⊙ / h for the unWISE samples at z ¯ ≈ 0.6 , 1.1, and 1.5, respectively. We find that all three samples are dominated by central galaxies, rather than satellites. Using our constrained HOD models, we infer the effective linear galaxy bias for each unWISE sample, and find that it does not evolve as steeply with redshift as found in previous perturbation-theory-based analyses of these galaxies. We discuss possible sources of systematic uncertainty in our results, the most significant of which is the uncertainty on the galaxy redshift distribution. Our HOD constraints provide a detailed, quantitative understanding of how the unWISE galaxies populate the underlying dark matter halo distribution. These constraints will have a direct impact on future studies employing the unWISE galaxies as a cosmological and astrophysical probe, including measurements of ionized gas thermodynamics and dark matter profiles via Sunyaev-Zel’dovich and lensing cross-correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Ensemble Effects on Hydroxide Bond Dissociation Free Energies in Polyoxovanadate Clusters

Understanding structure-property relationships is foundational to numerous modern chemistries, such as proton-coupled electron transfer (PCET). However, an experimentally measured property is the result of the behavior from an ensemble of molecules. Neglecting ensemble effects, especially under complex chemical environments, may obfuscate these relationships and lead to discrepancies between theory and experiment. In this work, we demonstrate the impact of configurational entropy and local chemical environments on hydroxide bond dissociation free energies [BDFE- (O−H)] for a set of polyoxovanadate nanoclusters, at ambient conditions. The O−H bond strengths are investigated via density functional theory (DFT) coupled with statistical thermodynamic analysis and bilinear modeling, and compared with previous experimental results on the same systems, namely electrochemical solutions of: [V 6 O 13−x (OH) x (TRIOL R ) 2 ] −2 (x = 2, 4, 6; R = NO 2 , Me) and [V 6 O 11−x (OMe) 2 (OH) x (TRIOL NO 2 ) 2 ] −2 (x = 2, 4). Interestingly, we find that ensemble effects, even at room temperature, can account for a significant portion of the BDFE(O−H) trend with the degree of reduction via H atom binding, which cannot be fully captured by single-structure, static DFT calculations. Moreover, we find that the ensemble effects may be replicated statistically, requiring only enumeration of energetically accessible H-binding sites. With the ensemble effects resolved, we present a simple bilinear model to reconcile remaining biases between experiment and ensemble-informed theory, which corelate with clusterspecific electronic environment differences. The bilinear model achieves outstanding accuracy vs experiments with a root-mean squared error of 0.4 kcal/mol. Finally, based on the physicochemical characteristics of hydrogen interaction with polyoxometalates, we present a simple methodology that captures the BDFE(O−H) trend while dramatically reducing required DFT calculations by 98% and achieving accuracy within 1 kcal/mol. Overall, this work elucidates the roles and structural origins of configurational entropy and chemical effects on polyoxometalate hydroxide bond energies, with potential applicability to various atomically precise metal oxide systems. Importantly, it introduces models for rapid and highly accurate property calculations in connection with experiments.

Cluster chemistry↗

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

ABSTRACT We implement euclidemulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the Markov chain Monte Carlo forecasting code montepython. We compare the performance of halofit, hmcode, and euclidemulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wavenumber cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wavenumber cut-off of $0.2\, h\, {\rm Mpc}^{-1}$. We find that euclidemulator1 is on average up to 17 per cent more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to 42 per cent and the equation of state of dark energy of up to 26 per cent, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5σ and 2σ (from below 1σ up to 6σ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of ≤1.4 in all cases.

Astronomy & Astrophysics↗