Halo Field Research Articles

Evaluating the variance of individual halo properties in constrained cosmological simulations

ABSTRACT Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation. We introduce an internal method for quantifying the extent to which the variance of individual halo properties is suppressed by the constraints imposed on the initial conditions. We apply it to the Constrained Simulations in BORG (CSiBORG) suite of 101 high-resolution realizations across the posterior probability distribution of initial conditions from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. The method is based on the overlap of the initial Lagrangian patch of a halo in one simulation with those in another, measuring the degree to which the haloes’ particles are initially coincident. This addresses the extent to which the imposed large-scale structure constraints reduce the variance of individual halo properties. We find consistent reconstructions of $M\gtrsim 10^{14}~\mathrm{M}_\odot \, h^{-1}$ haloes, indicating that the constraints from the BORG algorithm are sufficient to pin down the masses, positions, and peculiar velocities of clusters to high precision, though we do not assess how well they reproduce observations of the local Universe. The effect of the constraints tapers off towards lower mass, and the halo spins and concentrations are largely unconstrained at all masses. We document the advantages of evaluating halo consistency in the initial conditions and describe how the method may be used to quantify our knowledge of the halo field given galaxy survey data analysed through the lens of probabilistic inference machines such as BORG.

Discovery of a Metal-poor Red Giant Star with the Highest Ultralithium Enhancement

We present the discovery of 2MASS J05241392−0336543 (hereafter J0524−0336), a very metal-poor ([Fe/H] = −2.43 ± 0.16), highly r-process-enhanced ([Eu/Fe] = +1.34 ± 0.10) Milky Way halo field red giant star, with an ultrahigh Li abundance of A(Li, 3D, NLTE) = 6.15 ± 0.25 and [Li/Fe] = +7.64 ± 0.25, respectively. This makes J0524−0336 the most lithium-enhanced giant star discovered to date. We present a detailed analysis of the star’s atmospheric stellar parameters and chemical abundance determinations. Additionally, we detect indications of infrared excess, as well as observe variable emission in the wings of the Hα absorption line across multiple epochs, indicative of a potential enhanced mass-loss event with possible outflows. Our analysis reveals that J0524−0336 lies either between the bump and the tip of the red giant branch (RGB), or on the early asymptotic giant branch (e-AGB). We investigate the possible sources of lithium enrichment in J0524−0336, including both internal and external sources. Based on current models and on the observational evidence we have collected, our study shows that J0524−0336 may be undergoing the so-called lithium flash that is expected to occur in low-mass stars when they reach the RGB bump and/or the e-AGB.

Enhancing Morphological Measurements of the Cosmic Web with Delaunay Tessellation Field Estimation

The density fields constructed by traditional mass assignment methods are susceptible to irritating discreteness, which hinders morphological measurements of cosmic large-scale structure (LSS) through Minkowski functionals (MFs). To alleviate this issue, fixed-kernel smoothing methods are commonly used in the literature, at the expense of losing substantial structural information. In this work, we propose to measure MFs with the Delaunay tessellation field estimation (DTFE) technique, with the goal of maximizing the extraction of morphological information from sparse tracers. We perform our analyses starting from matter fields and progressively extending to halo fields. At the matter-field level, we elucidate how discreteness affects morphological measurements of LSS. Then, by comparing with the traditional Gaussian smoothing scheme, we preliminarily showcase the advantages of DTFE for enhancing measurements of MFs from sparse tracers. At the halo-field level, we first numerically investigate various systematic effects on MFs of DTFE fields, which are induced by finite voxel sizes, halo number densities, halo weightings, and redshift space distortions (RSDs), respectively. Then, we explore the statistical power of MFs measured with DTFE for extracting the cosmological information encoded in RSDs. We find that MFs measured with DTFE exhibit improvements by ∼2 orders of magnitude in discriminative power for RSD effects and by a factor of ∼3–5 in constraining power on the structure growth rate over the MFs measured with Gaussian smoothing. These findings demonstrate the remarkable enhancements in statistical power of MFs achieved by DTFE, showing enormous application potentials for our method in extracting various key cosmological information from galaxy surveys.

The Coherent Magnetic Field of the Milky Way

We present a suite of models of the coherent magnetic field of the Galaxy based on new divergence-free parametric functions describing the global structure of the field. The model parameters are fit to the latest full-sky Faraday rotation measures (RMs) of extragalactic sources and polarized synchrotron intensity (PI) maps from the Wilkinson Microwave Anisotropy Probe and Planck. We employ multiple models for the density of thermal and cosmic-ray electrons in the Galaxy, needed to predict the sky maps of RMs and PI for a given Galactic magnetic field (GMF) model. The robustness of the inferred properties of the GMF is gauged by studying many combinations of parametric field models and electron density models. We determine the pitch angle of the local magnetic field (11° ± 1°), explore the evidence for a grand-design spiral coherent magnetic field (inconclusive), determine the strength of the toroidal and poloidal magnetic halo fields below and above the disk (magnitudes the same for both hemispheres within ≈10%), set constraints on the half-height of the cosmic-ray diffusion volume (≥2.9 kpc), investigate the compatibility of RM- and PI-derived magnetic field strengths (compatible under certain assumptions), and check if the toroidal halo field could be created by the shear of the poloidal halo field due to the differential rotation of the Galaxy (possibly). A set of eight models is identified to help quantify the present uncertainties in the coherent GMF spanning different functional forms, data products, and auxiliary input. We present the corresponding sky maps of rates for axion–photon conversion in the Galaxy and deflections of ultrahigh-energy cosmic rays.

Cracking the Relation between Mass and 1P Star Fraction of Globular Clusters. I. Present-day Cluster Masses as a First Tool

The phenomenon of multiple stellar populations is exacerbated in massive globular clusters, with the fraction of first-population (1P) stars a decreasing function of the cluster present-day mass. We decipher this relation in far greater detail than has been done so far. We assume (i) a fixed stellar mass threshold for the formation of second-population (2P) stars, (ii) a power-law scaling F1P∝mecl−1 between the mass m ecl of newly formed clusters and their 1P star fraction F 1P, and (iii) a constant F 1P over time. The F 1P(m ecl) relation is then evolved up to an age of 12 Gyr for tidal field strengths representative of the entire Galactic halo. The 12 Gyr old model tracks cover the present-day distribution of Galactic globular clusters in the (mass, F 1P) space extremely well. The distribution is curtailed on its top right side by the scarcity of clusters at large Galactocentric distances and on its bottom left side by the initial scarcity of very high-mass clusters and dynamical friction. Given their distinct dissolution rates, “inner” and “outer” model clusters are offset from each other, as observed. The locus of Magellanic Clouds clusters in the (mass, F 1P) space is as expected for intermediate-age clusters evolving in a gentle tidal field. Given the assumed constancy of F 1P, we conclude that 2P stars do not necessarily form centrally concentrated. We infer a minimum mass of 4 · 105 M ⊙ for multiple-population clusters at secular evolution onset. This high-mass threshold severely limits the number of 2P stars lost from evolving clusters, thereby fitting the low 2P star fraction of the Galactic halo field.

Synchrotron intensity gradient revealing magnetic fields in galaxy clusters

Magnetic fields and their dynamical interplay with matter in galaxy clusters contribute to the physical properties and evolution of the intracluster medium. However, the current understanding of the origin and properties of cluster magnetic fields is still limited by observational challenges. In this article, we map the magnetic fields at hundreds-kpc scales of five clusters RXC J1314.4-2515, Abell 2345, Abell 3376, MCXC J0352.4-7401, and El Gordo using the synchrotron intensity gradient technique in conjunction with high-resolution radio observations from the Jansky Very Large Array (JVLA) and the Karoo Array Telescope (MeerKAT). We demonstrate that the magnetic field orientation of radio relics derived from synchrotron intensity gradient is in agreement with that obtained with synchrotron polarization. Most importantly, the synchrotron intensity gradient is not limited by Faraday depolarization in the cluster central regions and allows us to map magnetic fields in the radio halos of RXC J1314.4-2515 and El Gordo. We find that magnetic fields in radio halos exhibit a preferential direction along the major merger axis and show turbulent structures at higher angular resolution. The results are consistent with expectations from numerical simulations, which predict turbulent magnetic fields in cluster mergers that are stirred and amplified by matter motions.

First spectroscopic investigation of anomalous Cepheid variables

Context. Anomalous Cepheids (ACEPs) are intermediate-mass metal-poor pulsators that are mostly discovered in dwarf galaxies of the Local Group. However, recent Galactic surveys, including the Gaia Data Release 3, found a few hundred ACEPs in the Milky Way. Their origin is only poorly understood. Aims. We aim to investigate the origin and evolution of Galactic ACEPs by studying the chemical composition of their atmospheres for the first time. Methods. We used UVES at the Very Large Telescope to obtain high-resolution spectra for a sample of nine ACEPs belonging to the Galactic halo. We derived the abundances of 12 elements, C, Na, Mg, Si, Ca, Sc, Ti, Cr, Fe, Ni, Y, and Ba. We complemented these data with literature abundances from high-resolution spectroscopy for an additional three ACEPs that were previously incorrectly classified as type II Cepheids. This increased the sample to a total of 12 stars. Results. All the investigated ACEPs have an iron abundance [Fe/H] < −1.5 dex, as expected from theoretical predictions for these pulsators. The abundance ratios of the different elements to iron show that the chemical composition of ACEPs is generally consistent with that of the Galactic halo field stars, with the exception of sodium, which is found to be overabundant in 9 out of the 11 ACEPs where it was measured. This is very similar to the situation for second-generation stars in Galactic globular clusters. The same comparison with dwarf and ultra-faint satellites of the Milky Way reveals more differences than similarities. It is therefore unlikely that the bulk of Galactic ACEPs originated in a galaxy like this that subsequently dissolved into the Galactic halo. The principal finding of this work is the unexpected overabundance of sodium in ACEPs. We explored several hypotheses to explain this feature, finding that the most promising scenario is the evolution of low-mass stars in a binary system with either mass transfer or merging. Detailed modelling is needed to confirm this hypothesis.

Quijote-PNG: The Information Content of the Halo Mass Function

We study signatures of primordial non-Gaussianity (PNG) in the redshift-space halo field on nonlinear scales using a combination of three summary statistics, namely, the halo mass function (HMF), power spectrum, and bispectrum. The choice of adding the HMF to our previous joint analysis of the power spectrum and bispectrum is driven by a preliminary field-level analysis, in which we train graph neural networks on halo catalogs to infer the PNG f NL parameter. The covariance matrix and the responses of our summaries to changes in model parameters are extracted from a suite of halo catalogs constructed from the Quijote-png N-body simulations. We consider the three main types of PNG: local, equilateral, and orthogonal. Adding the HMF to our previous joint analysis of the power spectrum and bispectrum produces two main effects. First, it reduces the equilateral f NL predicted errors by roughly a factor of 2 while also producing notable, although smaller, improvements for orthogonal PNG. Second, it helps break the degeneracy between the local PNG amplitude, , and assembly bias, b ϕ , without relying on any external prior assumption. Our final forecasts for the PNG parameters are , , , on a cubic volume of , with a halo number density of , at z = 1, and considering scales up to .

CATS: The Hubble Constant from Standardized TRGB and Type Ia Supernova Measurements

The tip of the red giant branch (TRGB) provides a luminous standard candle for constructing distance ladders to measure the Hubble constant. In practice, its measurements via edge-detection response (EDR) are complicated by the apparent fuzziness of the tip and the multipeak landscape of the EDR. Previously, we optimized an unsupervised algorithm, Comparative Analysis of TRGBs, to minimize the variance among multiple halo fields per host without relying on individualized choices, achieving state-of-the-art ∼<0.05 mag distance measures for optimal data. Here we apply this algorithm to an expanded sample of SN Ia hosts to standardize these to multiple fields in the geometric anchor, NGC 4258. In concert with the Pantheon+ SN Ia sample, this analysis produces a (baseline) result of H 0 = 73.22 ± 2.06 km s−1 Mpc−1. The largest difference in H 0 between this and similar studies employing the TRGB derives from corrections for SN survey differences and local flows used in the most recent SN Ia compilations that were absent in earlier studies. The SN-related differences total ∼2.0 km s−1 Mpc−1. A smaller share, ∼1.4 km s−1 Mpc−1, results from the inhomogeneity of the TRGB calibration across the distance ladder. We employ a grid of 108 variants around the optimal TRGB algorithm and find that the median of the variants is 72.94 ± 1.98 km s−1 Mpc−1 with an additional uncertainty due to algorithm choices of 0.83 km s−1 Mpc−1. None of these TRGB variants result in an H 0 of less than 71.6 km s−1 Mpc−1.

Comparative Analysis of TRGBs (CATs) from Unsupervised, Multi-halo-field Measurements: Contrast is Key

The tip of the red giant branch (TRGB) is an apparent discontinuity of the luminosity function (LF) due to the end of the red giant evolutionary phase and is used to measure distances in the local universe. In practice, tip localization via edge detection response (EDR) relies on several methods applied on a case-by-case basis. It is hard to evaluate how individual choices affect a distance estimation using only a single host field while also avoiding confirmation bias. To devise a standardized approach, we compare unsupervised, algorithmic analyses of the TRGB in multiple halo fields per galaxy. We first optimize methods for the lowest field-to-field dispersion, including spatial filtering, smoothing, and weighting of LF, color band selection, and tip selection based on the number of likely RGB stars and the ratio of stars below versus above the tip (R). We find R, which we call the tip contrast, to be the most important indicator of the quality of EDR measurements; higher R selection can decrease field-to-field dispersion. Further, since R is found to correlate with the age or metallicity of the stellar population based on theoretical modeling, it might result in a displacement of the detected tip magnitude. We find a tip-contrast relation with a slope of −0.023 ± 0.0046 mag/ratio, an ∼5σ result that can be used to correct these variations in the detections. When using TRGB to establish a distance ladder, consistent TRGB standardization using tip-contrast relation across rungs is vital to make robust cosmological measurements.

Tidal debris from Omega Centauri discovered with unsupervised machine learning

ABSTRACT The gravitational interactions between the Milky Way and in-falling satellites offer a wealth of information about the formation and evolution of our Galaxy. In this paper, we explore the high-dimensionality of the GALAH DR3 plus Gaia eDR3 data set to identify new tidally stripped candidate stars of the nearby star cluster Omega Centauri ($\omega \, \mathrm{Cen}$). We investigate both the chemical and dynamical parameter space simultaneously, and identify cluster candidates that are spatially separated from the main cluster body, in regions where contamination by halo field stars is high. Most notably, we find candidates for $\omega \, \mathrm{Cen}$ scattered in the halo extending to more than 50° away from the main body of the cluster. Using a grid of simulated streams generated with $\omega \, \mathrm{Cen}$-like orbital properties, we then compare the on sky distribution of these candidates to the models. The results suggest that if $\omega \, \mathrm{Cen}$ had a similar initial mass as its present day mass, then we can place a lower limit on its time of accretion at tacc &amp;gt; 7 Gyr ago. Alternatively, if the initial stellar mass was significantly larger, as would be expected if $\omega \, \mathrm{Cen}$ is the remnant core of a dwarf Galaxy, then we can constrain the accretion time to tacc &amp;gt; 4 Gyr ago. Taken together, these results are consistent with the scenario that $\omega \, \mathrm{Cen}$ is the remnant core of a disrupted dwarf galaxy.

Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field

We study primordial non-Gaussian signatures in the redshift-space halo field on nonlinear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogs constructed from the Quijote-png N-body simulations, which we release to accompany this paper. We verify its unbiasedness and near-optimality for the three main types of primordial non-Gaussianity (PNG): local, equilateral, and orthogonal. We compare the modal bispectrum expansion with a k-binning approach, showing that the former allows for faster convergence of numerical derivatives in the computation of the score function, thus leading to better final constraints. We find, in agreement with previous studies, that the local PNG signal in the halo field is dominated by the scale-dependent bias signature on large scales and saturates at k ∼ 0.2 h Mpc−1, whereas the small-scale bispectrum is the main source of information for equilateral and orthogonal PNG. Combining the power spectrum and bispectrum on nonlinear scales plays an important role in breaking degeneracies between cosmological and PNG parameters; such degeneracies, however, remain strong for equilateral PNG. We forecast that PNG parameters can be constrained with , , and on a cubic volume of at z = 1, considering scales up to .

Non-parametric Lagrangian biasing from the insights of neural nets

We present a Lagrangian model of galaxy clustering bias in which we train a neural net using the local properties of the smoothed initial density field to predict the late-time mass-weighted halo field.By fitting the mass-weighted halo field in the AbacusSummit simulations at z = 0.5, we find that including three coarsely spaced smoothing scales gives the best recovery of the halo power spectrum. Adding more smoothing scales may lead to 2–5% underestimation of the large-scale power and can cause the neural net to overfit.We find that the fitted halo-to-mass ratio can be well described by two directions in the original high-dimension feature space.Projecting the original features into these two principal components and re-training the neural net either reproduces the original training result, or outperforms it with a better match of the halo power spectrum. The elements of the principal components are unlikely to be assigned physical meanings, partly owing to the features being highly correlated between different smoothing scales.Our work illustrates a potential need to include multiple smoothing scales when studying galaxy bias, and this can be done easily with machine-learning methods that can take in high dimensional input feature space.

RR Lyrae Mid-infrared Period–Luminosity–Metallicity and Period–Wesenheit–Metallicity Relations Based on Gaia DR3 Parallaxes

We present new empirical infrared period–luminosity–metallicity (PLZ) and period–Wesenheit–metallicity (PWZ) relations for RR Lyae based on the latest Gaia Early Data Release 3 (EDR3) parallaxes. The relations are provided in the Wide-field Infrared Survey Explorer (WISE) W1 and W2 bands, as well as in the W(W1, V − W1) and W(W2, V − W2) Wesenheit magnitudes. The relations are calibrated using a very large sample of Galactic halo field RR Lyrae stars with homogeneous spectroscopic [Fe/H] abundances (over 1000 stars in the W1 band), covering a broad range of metallicities (−2.5 ≲ [Fe/H] ≲ 0.0). We test the performance of our PLZ and PWZ relations by determining the distance moduli of both galactic and extragalactic stellar associations: the Sculptor dwarf spheroidal galaxy in the Local Group (finding ), the Galactic globular clusters M4 (), and the Reticulum globular cluster in the Large Magellanic Cloud (). The distance moduli determined through all our relations are internally self-consistent (within ≲0.05 mag) but are systematically smaller (by ∼2–3σ) than previous literature measurements taken from a variety of methods/anchors. However, a comparison with similar recent RR Lyrae empirical relations anchored with EDR3 likewise shows, to varying extents, a systematically smaller distance modulus for PLZ/PWZ RR Lyrae relations.

Electromagnetic load evaluation of K-DEMO divertor for MD and VDE plasma disruption scenarios

Electromagnetic (EM) force caused by the instability of plasma operation is one of the most significant external loads in the DEMO, as well as the extremely high heat flux (more than 10 MW/m2). Abrupt disruption of plasma current of several MA causes thermal quench and current quench sequentially. This paper presents the EM force on the divertor calculated by considering the Poloidal field variation (PFV), halo current (HC), and toroidal field variation (TFV) generated during the quenching process for the K-DEMO. The K-DEMO with the fusion power of 2,200 MW has a major radius of 6.8 m and a plasma current of 12 MA in steady-state operation. The K-DEMO divertor basically employs a water-cooled tungsten monoblock concept with a reduced activation ferritic martensitic RAFM steel heat sink. EM force was calculated using ANSYS Mechanical/Emag solver. For the convenience of calculation, the existing inconvenient GUI environment based on MAPDL was changed to ANSYS Workbench and improved. To evaluate the EM forces generated in the K-DEMO divertor when various disruption scenarios occur, three scenarios were selected for major disruption (MD) and vertical displacement event (VDE), respectively. The variation of PFV, HC, and TFV in the disruption are significant input data for EM analysis. The major disruption data were converted from ITER disruption scenarios by using the in-house code since the plasma data for the disruption scenarios has not been developed for the K-DEMO. This study presents the preliminary evaluation of EM force for the K-DEMO and improved calculation platform although the disruption data derived from ITER are used. The EM forces generated by the selected six scenarios were calculated and compared. In the MD_DW_slow_fast and VDE_DW_slow_fast scenarios, the maximum EM force was 4.8 MN on the K-DEMO divertor module along the vertical direction, and the HC influence was dominant.

Quijote-PNG: The Information Content of the Halo Power Spectrum and Bispectrum

We investigate how much can be learnt about four types of primordial non-Gaussianity (PNG) from small-scale measurements of the halo field. Using the quijote-png simulations, we quantify the information content accessible with measurements of the halo power spectrum monopole and quadrupole, the matter power spectrum, the halo–matter cross spectrum, and the halo bispectrum monopole. This analysis is the first to include small, nonlinear scales, up to , and to explore whether these scales can break degeneracies with cosmological and nuisance parameters making use of thousands of N-body simulations. We perform all the halo measurements in redshift space with a single sample comprised of all halos with mass >3.2 × 1013 h −1 M ⊙. For local PNG, measurements of the scale-dependent bias effect from the power spectrum using sample variance cancellation provide significantly tighter constraints than measurements of the halo bispectrum. In this case measurements of the small scales add minimal additional constraining power. In contrast, the information on equilateral and orthogonal PNG is primarily accessible through the bispectrum. For these shapes, small-scale measurements increase the constraining power of the halo bispectrum by up to 4×, though the addition of scales beyond k ≈ 0.3 h Mpc−1 improves constraints largely through reducing degeneracies between PNG and the other parameters. These degeneracies are even more powerfully mitigated through combining power spectrum and bispectrum measurements. However, even with combined measurements and small-scale information, equilateral non-Gaussianity remains highly degenerate with σ 8 and our bias model.

Anisotropy and characteristic scales in halo density gradient profiles

We use a large N-body simulation to study the characteristic scales in the density gradient profiles in and around halos with masses ranging from 1012 to 1015 M⊙. We investigate the profiles separately along the major (T1) and minor (T3) axes of the local tidal tensor and how the characteristic scales depend on halo mass, formation time, and environment. We find two prominent features in the gradient profiles: a deep “valley” and a prominent “peak.” We use the Gaussian process regression to fit the gradient profiles and identify the local extrema in order to determine the scales associated with these features. Around the valley, we identify three types of distinct local minima, corresponding to caustics of particles orbiting around halos. The appearance and depth of the three caustics depend on the direction defined by the local tidal field, formation time, and environment of halos. The first caustic is located at r &gt; 0.8R200, corresponding to the splashback feature, and is dominated by particles at their first apocenter after infall. The second and third caustics, around 0.6R200 and 0.4R200, respectively, can be determined reliably only for old halos. The three caustics are consistent with the prediction of self-similar gravitational collapse. The first caustic is always the most prominent feature along T3, but may not be true along T1 or in azimuthally averaged profiles, suggesting that caution must be taken when using averaged profiles to investigate the splashback radius. We find that the splashback feature is approximately isotropic when proper separations are made between the first and the other caustics. We also identify a peak feature located at ∼2.5R200 in the density gradient profile. This feature is the most prominent along T1 and is produced by mass accumulations from the structure outside halos. We also discuss the origins of these features and their observational implications.

The effect of local Universe constraints on halo abundance and clustering

ABSTRACT Cosmological N-body simulations of the dark matter component of the universe typically use initial conditions with a fixed power spectrum and random phases of the density field, leading to structure consistent with the local distribution of galaxies only in a statistical sense. It is, however, possible to infer the initial phases which lead to the configuration of galaxies and clusters that we see around us. We analyse the CSiBORG suite of 101 simulations, formed by constraining the density field within 155 Mpc h−1 with dark matter particle mass 4.38 × 109 M⊙, to quantify the degree to which constraints imposed on 2.65 Mpc h−1 scales reduce variance in the halo mass function and halo–halo cross-correlation function on a range of scales. This is achieved by contrasting CSiBORG with a subset of the unconstrained Quijote simulations and expectations for the ΛCDM average. Using the FOF, PHEW, and HOP halofinders, we show that the CSiBORG suite beats cosmic variance at large mass scales (≳1014 M⊙ h−1), which are most strongly constrained by the initial conditions, and exhibits a significant halo–halo cross-correlation out to ∼30 Mpc h−1. Moreover, the effect of the constraints percolates down to lower mass objects and to scales below those on which they are imposed. Finally, we develop an algorithm to ‘twin’ haloes between realizations and show that approximately 50 per cent of haloes with mass greater than 1015 M⊙ h−1 can be identified in all realizations of the CSiBORG suite. We make the CSiBORG halo catalogues publicly available for future applications requiring knowledge of the local halo field.

Beware of fakeν’s: The effect of massive neutrinos on the nonlinear evolution of cosmic structure

Massive neutrinos suppress the growth of cosmic structure on small, nonlinear scales. It is thus often proposed that using statistics beyond the power spectrum can tighten constraints on the neutrino mass by extracting additional information from these nonlinear scales. We study the information content regarding neutrino mass at the field level, quantifying how much of this information arises from the difference in nonlinear evolution between a cosmology with one fluid [cold dark matter (CDM)] and two fluids ($\mathrm{CDM}+\text{neutrinos}$). We do so by running two $N$-body simulations, one with and one without massive neutrinos, both with the same phases, and matching their linear power spectrum at a given low redshift. This effectively isolates the information encoded in the linear initial conditions from the nonlinear cosmic evolution. We demonstrate that, for $k\ensuremath{\lesssim}1\text{ }\text{ }h/\mathrm{Mpc}$, and for a single redshift, there is negligible difference in the real-space CDM field between the two simulations. This suggests that all the information regarding neutrino mass is in the linear power spectrum set by the initial conditions. Thus, any probe based on the CDM field alone will have negligible constraining power beyond that which exists at the linear level over the same range of scales. Consequently, any probe based on the halo field will contain little information beyond the linear power. We find similar results for the matter field responsible for weak lensing. We also demonstrate that there may be much information beyond the power spectrum in the 3D matter field; however, this is not observable in modern surveys via dark matter halos or weak lensing. Finally, we show that there is additional information to be found in redshift space.

Metallicity of Galactic RR Lyrae from Optical and Infrared Light Curves. II. Period–Fourier–Metallicity Relations for First Overtone RR Lyrae

We present new period-ϕ 31-[Fe/H] relations for first-overtone RRL stars (RRc), calibrated over a broad range of metallicities (−2.5 ≲ [Fe/H] ≲ 0.0) using the largest currently available set of Galactic halo field RRL with homogeneous spectroscopic metallicities. Our relations are defined in the optical (ASAS-SN V band) and, inaugurally, in the infrared (WISE W1 and W2 bands). Our V-band relation can reproduce individual RRc spectroscopic metallicities with a dispersion of 0.30 dex over the entire metallicity range of our calibrator sample (an rms smaller than what we found for other relations in literature including nonlinear terms). Our infrared relation has a similar dispersion in the low- and intermediate-metallicity range ([Fe/H] ≲ −0.5), but tends to underestimate the [Fe/H] abundance around solar metallicity. We tested our relations by measuring both the metallicity of the Sculptor dSph and a sample of Galactic globular clusters, rich in both RRc and RRab stars. The average metallicity we obtain for the combined RRL sample in each cluster is within ±0.08 dex of their spectroscopic metallicities. The infrared and optical relations presented in this work will enable deriving reliable photometric RRL metallicities in conditions where spectroscopic measurements are not feasible; e.g., in distant galaxies or reddened regions (observed with upcoming Extremely Large Telescopes and the James Webb Space Telescope), or in the large sample of new RRL that will be discovered in large-area time-domain photometric surveys (such as the LSST and the Roman space telescope).