Search

Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to investigate the sensitivity of the central entropy and the shape of the profiles to changes in the sub-grid model applied to a suite of zoom-in cosmological simulations of a group of mass $M_{500} = 8.8 \times 10^{12}~{\rm M}_\odot$ and a cluster of mass $2.9 \times 10^{14}~{\rm M}_\odot$. Using our reference model, calibrated to match the stellar mass function of field galaxies, we confirm that our simulated groups and clusters contain hot gas with too high entropy in their cores. Additional simulations run without artificial conduction, metal cooling or AGN feedback produce lower entropy levels but still fail to reproduce observed profiles. Conversely, the two objects run without supernova feedback show a significant entropy increase which can be attributed to excessive cooling and star formation. Varying the AGN heating temperature does not greatly affect the profile shape, but only the overall normalisation. Finally, we compared runs with four AGN heating schemes and obtained similar profiles, with the exception of bipolar AGN heating, which produces a higher and more uniform entropy distribution. Our study leaves open the question of whether the entropy core problem in simulations, and particularly the lack of power-law cool-core profiles, arise from incorrect physical assumptions, missing physical processes, or insufficient numerical resolution., Accepted for publication in MNRAS on 2023 February 06. The code used in the analysis is publicly available on the corresponding author's GitHub repository: https://github.com/edoaltamura/entropy-core-problem

We present the deconvolved distribution estimator (DDE), an extension of the voxel intensity distribution (VID), in the context of future observations proposed as part of the CO Mapping Array Project (COMAP). The DDE exploits the fact that the observed VID is a convolution of correlated signal intensity distributions and uncorrelated noise or interloper intensity distributions. By deconvolving the individual VID of two observables away from their joint VID in a Fourier-space operation, the DDE suppresses sensitivity to interloper emission while maintaining sensitivity to correlated components. The DDE thus improves upon the VID by reducing the relative influence of uncorrelated noise and interloper biases, which is useful in the context of COMAP observations that observe different rotational transitions of CO from the same comoving volume in different observing frequency bands. Fisher forecasts suggest that the theoretical sensitivity in the DDE allows significant improvements in constraining power compared to either the cross power spectrum or the individual VID data, and matches the constraining power of the combination of all other one- and two-point summary statistics. Future work should further investigate the covariance and model-dependent behaviour of this novel one-point cross-correlation statistic., Comment: 11 pages + acknowledgements/bibliography/appendix (12 pages total); 10 figures, 2 tables; v3 reflects the version accepted for publication in MNRAS

We present high angular-resolution measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward two galaxy clusters, RCS J2319+0038 at z=0.9 and HSC J0947-0119 at z=1.1, by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. They are supplemented with available Chandra X-ray data, optical data taken by Hyper Suprime-Cam on Subaru, and millimeter-wave SZE data from the Atacama Cosmology Telescope. Taking into account departures from spherical symmetry, we have reconstructed non-parametrically the inner pressure profile of two clusters as well as electron temperature and density profiles for RCS J2319+0038. This is one of the first such measurements for an individual cluster at $z \gtrsim 0.9$. We find that the inner pressure profile of both clusters is much shallower than that of local cool-core clusters. Our results consistently suggest that RCS J2319+0038 hosts a weak cool core, where radiative cooling is less significant than in local cool cores. On the other hand, HSC J0947-0119 exhibits an even shallower pressure profile than RCS J2319+0038 and is more likely a non-cool-core cluster. The SZE centroid position is offset by more than 140 $h_{70}^{-1}$kpc from the peaks of galaxy distribution in HSC J0947-0119, suggesting a stronger influence of mergers in this cluster. We conclude that these distant clusters are at a very early stage of developing the cool cores typically found in clusters at lower redshifts., Comment: 28 pages, 14 figures, 10 tables, revised after referee comments, accepted for publication in PASJ

Primordial black holes (PBHs) form from large-amplitude initial density fluctuations and may comprise some or all of the dark matter. If PBHs have a broadly extended mass spectrum, or in mixed PBH-particle dark matter scenarios, the extreme density fluctuations necessary to produce PBHs also lead to the formation of a much greater abundance of dark matter minihaloes that form during the radiation epoch with internal densities potentially of order $10^{12}$ M$_\odot$ pc$^{-3}$. We develop an analytic description of the formation of these ultradense haloes and use it to quantitatively compare PBH and halo distributions. PBHs that contribute only a per cent level fraction of the dark matter are accompanied by ultradense haloes that nevertheless comprise an order-unity fraction. These haloes would consist of either particle dark matter or much smaller PBHs. This finding significantly alters the predictions of many PBH scenarios, enabling a variety of new observational tests., 6 pages, 3 figures; accepted for publication in MNRAS. Includes expanded discussions of the connection between PBH and halo masses and of the outlook for these results

Where the cosmic baryons lie in and around galactic dark matter halos is only weakly constrained. We develop a method to quickly paint on models for their distribution. Our approach uses the statistical advantages of $N$-body simulations, while painting on the profile of gas around individual halos in ways that can be motivated by semi-analytic models or zoom-in hydrodynamic simulations of galaxies. Possible applications of the algorithm include extragalactic dispersion measures to fast radio bursts (FRBs), the Sunyaev-Zeldovich effect, baryonic effects on weak lensing, and cosmic metal enrichment. As an initial application, we use this tool to investigate how the baryonic profile of foreground galactic-mass halos affects the statistics of the dispersion measure (DM) towards cosmological FRBs. We show that the distribution of DM is sensitive to the distribution of baryons in galactic halos, with viable gas profile models having significantly different probability distributions for DM to a given redshift. We also investigate the requirements to statistically measure the circumgalactic electron profile for FRB analyses that stack DM with impact parameter to foreground galaxies, quantifying the size of the contaminating "two-halo" term from correlated systems and the number of FRBs for a high significance detection. Publicly available Python modules implement our CGMBrush algorithm., Comment: 23 pages, 10 figures, public python code at https://github.com/ianw89/cgm-brush

The early dark energy (EDE) solution to the Hubble tension comes at the cost of an increased clustering amplitude that has been argued to worsen the fit to galaxy clustering data. We explore whether freeing the total neutrino mass $M_{\nu}$, which can suppress small-scale structure growth, improves EDE's fit to galaxy clustering. Using Planck Cosmic Microwave Background and BOSS galaxy clustering data, a Bayesian analysis shows that freeing $M_{\nu}$ does not appreciably increase the inferred EDE fraction $f_{\rm EDE}$: we find the 95% C.L. upper limits $f_{\rm EDE}, Comment: 9 pages, 6 figures. v3: included additional discussions and 3 figures which further strengthen our results and their interpretation, added further references, slight format change. Accepted for publication in MNRAS. and associated Jupyter notebook available at https://github.com/alexreevesy/ede_mnu

As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modeling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the SNANA simulation framework to construct the Simulated Catalogue of Optical Transients and Correlated Hosts (SCOTCH), a publicly-available catalogue of 5 million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range $0, Comment: 30 pages, 21 figures. Updated to reflect version published in MNRAS Vol. 520, Issue 2. This version improves the treatment of transient-host offsets and AGN host correlations. Associated SCOTCH files are available through Zenodo at https://doi.org/10.5281/zenodo.7563623

Cosmological parameter constraints from recent galaxy imaging surveys are reaching $2-3\%$-level accuracy. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will produce sub-percent level measurements of cosmological parameters, providing a milestone test of the $\Lambda$CDM model. To supply guidance to the upcoming LSST analysis, it is important to understand thoroughly the results from different recent galaxy imaging surveys and assess their consistencies. In this work we perform a unified catalog-level reanalysis of three cosmic shear datasets: the first year data from the Dark Energy Survey (DES-Y1), the 1,000 deg$^{2}$ dataset from the Kilo-Degree Survey (KiDS-1000), and the first year data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-Y1). We utilize a pipeline developed and rigorously tested by the LSST Dark Energy Science Collaboration to perform the reanalysis and assess the robustness of the results to analysis choices. We find the $S_{8}$ constraint to be robust to two different small-scale modeling approaches, and varying choices of cosmological priors. Our unified analysis allows the consistency of the surveys to be rigorously tested and we find the three surveys to be statistically consistent. Due to the partially overlapping footprint, we model the cross-covariance between KiDS-1000 and HSC-Y1 approximately when combining all three datasets, resulting in a $1.6-1.9\%$ constraint on $S_8$ given different assumptions on the cross-covariance., Comment: 27 pages, 18 figures, to be submitted to MNRAS

We introduce a modification of the Press-Schechter formalism aimed to derive general mass functions for primordial black holes (PBHs). In this case, we start from primordial power spectra (PPS) which include a monochromatic spike, typical of ultra slow-roll inflation models. We consider the PBH formation as being associated to the amplitude of the spike on top of the linear energy density fluctuations, coming from a PPS with a blue index. By modelling the spike with a log-normal function, we study the properties of the resulting mass function spikes, and compare these to the underlying extended mass distributions. When the spike is at PBH masses which are much lower than the exponential cutoff of the extended distribution, very little mass density is held by the PBHs within the spike, and it is not ideal to apply the Press-Schechter formalism in this case as the resulting characteristic overdensity is too different from the threshold for collapse. It is more appropriate to do so when the spike mass is similar to, or larger than the cutoff mass. Additionally, it can hold a similar mass density as the extended part. Such particular mass functions also contain large numbers of small PBHs, especially if stable PBH relics are considered, and they can provide $\sim 1000M_\odot$ seeds for the supermassive black holes at the centres of present-day galaxies. The constraints on the fraction of dark matter in PBHs for monochromatic mass functions are somewhat relaxed when there is an additional underlying extended distribution of masses., 14 pages, 4 figures, 1 table, refereces added, discussion extended and conclusion unchanged. Accepted for publication in MNRAS

We predict the 21-cm global signal and power spectra during the Epoch of Reionisation using the MERAXES semi-analytic galaxy formation and reionisation model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionisation of cosmic hydrogen using semi-analytic models (such as MERAXES) requires N-body simulations within large volumes and high mass resolutions. For this, we use a simulation of side-length $210~h^{-1}$ Mpc with $4320^3$ particles resolving dark matter haloes to masses of $5\times10^8~h^{-1}~M_\odot$. To reach the mass resolution of atomically cooled galaxies, thought to be the dominant population contributing to reionisation, at $z=20$ of $\sim 2\times10^7~h^{-1}~M_\odot$, we augment this simulation using the DARKFOREST Monte-Carlo merger tree algorithm (achieving an effective particle count of $\sim10^{12}$). Using this augmented simulation we explore the impact of mass resolution on the predicted reionisation history as well as the impact of X-ray heating on the 21-cm global signal and the 21-cm power spectra. We also explore the cosmic variance of 21-cm statistics within $70^{3}$ $h^{-3}$ Mpc$^3$ sub-volumes. We find that the midpoint of reionisation varies by $\Delta z\sim0.8$ and that the cosmic variance on the power spectrum is underestimated by a factor of $2-4$ at $k\sim 0.1-0.4$ Mpc$^{-1}$ due to the non-Gaussian nature of the 21-cm signal. To our knowledge, this work represents the first model of both reionisation and galaxy formation which resolves low-mass atomically cooled galaxies while simultaneously sampling sufficiently large scales necessary for exploring the effects of X-rays in the early Universe., Comment: 17 pages, 13 figures, revised and MNRAS accepted version

Inflationary gravitational waves, behaving as additional radiation in the Early Universe, can increase the effective number of relativistic species ($N_{\rm eff}$) by a further correction that depends on the integrated energy-density in gravitational waves over all scales. This effect is typically used to constrain (blue-tilted) models of inflation in light of the bounds resulting from the Big Bang Nucleosynthesis. In this paper, we recompute this contribution, discussing some caveats of the state-of-the-art analyses. Through a parametric investigation, we first demonstrate that the calculation is dominated by the ultraviolet frequencies of the integral and therefore by the behavior of the tensor spectrum on scales corresponding to modes that cross the horizon very close to the end of inflation, when the slow-roll dynamics breaks down and the production of gravitational waves becomes strongly model dependent. Motivated by these results, we realize a theoretical Monte Carlo and, working within the framework of the Effective Field Theory of inflation, we investigate the observable predictions of a very broad class of models. For each model, we solve a system of coupled differential equations whose solution completely specifies the evolution of the spectrum up to the end of inflation. We prove the calculation of $\Delta N_{\rm eff}^{\rm GW}$ to be remarkably model-dependent and therefore conclude that accurate analyses are needed to infer reliable information on the inflationary Universe., Comment: 17 pages, 7 figures, 1 table. Updated to match the MNRAS published version

Modeling the 21cm global signal from the Cosmic Dawn is challenging due to the many poorly constrained physical processes that come into play. We address this problem using the semi-analytical code "Cosmic Archaeology Tool" (CAT). CAT follows the evolution of dark matter halos tracking their merger history and provides an ab initio description of their baryonic evolution, starting from the formation of the first (Pop III) stars and black holes (BHs) in mini-halos at z > 20. The model is anchored to observations of galaxies and AGN at z < 6 and predicts a reionization history consistent with constraints. In this work we compute the evolution of the mean global 21cm signal between $4\leq z \leq 40$ based on the rate of formation and emission properties of stars and accreting black holes. We obtain an absorption profile with a maximum depth $δ{\rm T_b} = -95$ mK at $z \sim 26.5$ (54 MHz). This feature is quickly suppressed turning into an emission signal at $z = 20$ due to the contribution of accreting BHs that efficiently heat the IGM at $z < 27$. The high-$z$ absorption feature is caused by the early coupling between the spin and kinetic temperature of the IGM induced by Pop III star formation episodes in mini-halos. Once we account for an additional radio background from early BHs, we are able to reproduce the timing and the depth of the EDGES signal only if we consider a smaller X-ray background from accreting BHs, but not the shape., 16 pages, 10 figures, Accepted for publication in MNRAS

A key unknown of the Milky Way (MW) satellites is their orbital history, and, in particular, the time they were accreted onto the MW system since it marks the point where they experience a multitude of environmental processes. We present a new methodology for determining infall times, namely using a neural network (NN) algorithm. The NN is trained on MW-analogues in the EAGLE hydrodynamical simulation to predict if a dwarf galaxy is at first infall or a backsplash galaxy and to infer its infall time. The resulting NN predicts with 85\% accuracy if a galaxy currently outside the virial radius is a backsplash satellite and determines the infall times with a typical 68\% confidence interval of 4.4 Gyrs. Applying the NN to MW dwarfs with Gaia EDR3 proper motions, we find that all of the dwarfs within 300 kpc had been inside the Galactic halo. The overall MW satellite accretion rate agrees well with the theoretical prediction except for late times when the MW shows a second peak at a lookback time of 1.5 Gyrs corresponding to the infall of the LMC and its satellites. We also find that the quenching times for ultrafaint dwarfs show no significant correlation with infall time and thus supporting the hypothesis that they were quenched during reionisation. In contrast, dwarfs with stellar masses above $10^5~M_\odot$ are found to be consistent with environmental quenching inside the Galactic halo, with star-formation ceasing on average at $0.5^{+0.9}_{-1.2}$ Gyrs after infall., 13 pages, 7 figures, appendix: 4 pages, published version

A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work we calculate the ionizing escape fractions ($f_{\rm esc}$) of reionization-era galaxies from the state-of-the-art THESAN simulations, which combine an accurate radiation-hydrodynamic solver AREPO-RT with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($L_{\rm box} = 95.5\,\rm cMpc$). This allows the formation of numerous massive haloes ($M_{\rm halo} \gtrsim 10^{10}\,{\rm M_{\odot}}$), which are often statistically underrepresented in previous studies but are believed to be important to achieve rapid reionization. We find that low-mass galaxies ($M_{\rm stars} \lesssim 10^7\,{\rm M_{\odot}}$) are the main drivers of reionization above $z \gtrsim 7$, while high-mass galaxies ($M_{\rm stars} \gtrsim 10^8\,{\rm M_{\odot}}$) dominate the escaped ionizing photon budget at lower redshifts. The variation in halo escape fractions decreases for higher-mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. We show that dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global $f_{\rm esc}$ depends on the dust model. Finally, AGN are unimportant for reionization in THESAN and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells., 24 pages, 19 figures, MNRAS. Please visit www.thesan-project.com for more details

General Relativity (GR) has been successfully tested mainly at Solar system scales; however, galaxy-scale tests have become popular in the last few decades. In this work, we investigate the $\eta_\text{PPN}$ parameter, which is commonly defined by the ratio of two scalar potentials that appears in the cosmological linearly perturbed metric. Under the assumption of GR and a vanish anisotropic stress tensor, $\eta_\text{PPN}= 1$. Using ALMA, HST, and VLT/MUSE data, we combine mass measurements, using gravitational lensing and galactic dynamics, for the SDP.81 lens galaxy ($z = 0.299$) to constrain $\eta_\text{PPN}$. By using a flexible and self-consistent mass profile, our fiducial model takes into account the contribution of the stellar mass and a dark matter halo to reconstruct the lensed galaxy and the spatially-resolved stellar kinematics. We infer, after accounting for systematic uncertainties related to the mass model, cosmology and kinematics, $\eta_{\text{PPN}} = 1.13^{+0.03}_{-0.03}\pm0.20\,(\text{sys})$, which is in accordance with GR predictions. Better spectroscopy data are needed to push the systematics down and bring the uncertainty to the percentage level since our analysis shows that the main source of the systematics is related to kinematics, which heavily depends on the signal-to-noise ratio of the spectra., Comment: 18 pages, 8 figures. Accepted to MNRAS

Molecular hydrogen allows cooling in primordial gas, facilitating its collapse into Population III stars within primordial halos. Lyman-Werner (LW) radiation from these stars can escape the halo and delay further star formation by destroying H$_2$ in other halos. As cosmological simulations show that increasing the background LW field strength increases the average halo mass required for star formation, we perform follow-up simulations of selected halos to investigate the knock-on effects this has on the Population III IMF. We follow 5 halos for each of the $J_{21}$ = 0, 0.01 and 0.1 LW field strengths, resolving the pre-stellar core density of $10^{-6}$ g cm$^{-3}$ (10$^{18}$ cm$^{-3}$) before inserting sink particles and following the fragmentation behaviour for hundreds of years further. We find that the mass accreted onto sinks by the end of the simulations is proportional to the mass within the $\sim 10^{-2}$ pc molecular core, which is not correlated to the initial mass of the halo. As such, the IMFs for masses above the brown dwarf limit show little dependence on the LW strength, although they do show variance in the number of low-mass clumps formed. As the range of background LW field strengths tested here covers the most likely values from literature, we conclude that the IMF for so-called Pop III.2 stars is not significantly different from the initial population of Pop III.1 stars. The primordial IMF therefore likely remains unchanged until the formation of the next generation of Population II stars., Comment: MNRAS Accepted 2023 January 16 ref. MN-22-5075-MJ.R2

We consider the unified bulk viscous scenarios and constrain them using the Cosmic Microwave Background observations from Planck 2018 and the Pantheon sample from Type Ia Supernovae. Then we generate the luminosity distance measurements from ${\cal O}(10^3)$ mock Gravitational Wave Standard Sirens (GWSS) events for the proposed Einstein Telescope. We then combine these mock luminosity distance measurements from the GWSS with the current cosmological probes in order to forecast how the mock GWSS data could be effective in constraining these bulk viscous scenarios. Our results show that a non-zero time dependent bulk viscosity in the universe sector is strongly preferred by the current cosmological probes and will possibly be confirmed at many standard deviations by the future GWSS measurements. We further mention that the addition of GWSS data can significantly reduce the uncertainties of the key cosmological parameters obtained from the usual cosmological probes employed in this work., Comment: 11 pages, 3 tables and 4 figures. Version accepted for publication in MNRAS

We study the global agreement between the most recent observations of the Cosmic Microwave Background temperature and polarization anisotropies angular power spectra released by the Atacama Cosmology Telescope and the Planck satellite in various cosmological models that differ by the inclusion of different combinations of additional parameters. By using the Suspiciousness statistic, we show that the global "CMB tension" between the two experiments, quantified at the Gaussian equivalent level of $\sim 2.5\,\sigma$ within the baseline $\Lambda$CDM, is reduced at the level of $1.8\sigma$ when the effective number of relativistic particles ($N_{\rm eff}$) is significantly less than the standard value, while it ranges between $2.3\,\sigma$ and $3.5\,\sigma$ in all the other extended models., Comment: 6 pages, 1 figure, 2 tables. Updated to match the MNRAS published version

The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe. Specifically, we use mock catalogues built from the TNG300 hydrodynamical simulation to show that SHAMe can simultaneously model the multipoles of the redshift-space galaxy correlation function and galaxy-galaxy lensing, without noticeable bias within the statistical sampling uncertainties of a SDSS volume and on scales r = [0.6-30] Mpc/h. Modelling the baryonic processes in galaxy-galaxy lensing with a baryonification scheme allows SHAMe's range of validity to be extended to r = [0.1-30] Mpc/h. Remarkably, our model achieves this level of precision with just five free parameters beyond those describing the baryonification model. At fixed cosmology, we find that galaxy-galaxy lensing provides a general consistency test but little additional information on galaxy modelling parameters beyond that encoded in the redshift-space multipoles. It does, however, improve constraints if only the projected correlation function is available, as in surveys with only photometric redshifts. We expect SHAMe to have a higher fidelity across a wider range of scales than more traditional methods such as Halo Occupation Distribution modelling. Thus it should provide a significantly more powerful and more robust tool for analysing next-generation large-scale surveys., 14 pages, 7 figures. Accepted by MNRAS

Due to its proximity, the stellar populations of the Galactic bulge (GB) can be resolved and can be studied in detail. This allows tracing the bulge metallicity distribution function (MDF) for different spatial regions within the bulge, which may give us clues about the bulge formation and evolution scenarios. In this work, we developed a chemical evolution model (CEM), taking into account the mass distribution in the bulge and disc, to derive the radial dependence of this time-scale in the Galaxy. Since the infall rate depends on that time scale in the CEM, the results of the model were used to test a scenario where the bulge is formed inside-out. The obtained results for the $[\alpha/\mbox{Fe}]$ vs. [Fe/H] relationship, the MDF and the [Fe/H] radial gradient in the bulge have been compared to available data in the literature. The model is able to reproduce most of the observational data: the spread in the relation $[\alpha/\mbox{Fe}]$ vs. [Fe/H], the MDF shape in different regions of the bulge, the [Fe/H] radial gradient inside it and the age-metallicity relation, as well as the [$\alpha$/Fe] evolution with age. The results of the model point to a scenario where the bulk of the bulge stars pre-existed the boxy/peanut X-shape bar formation. As a result, the classical origin of the GB is not ruled out and this scenario may be invoked to explain the chemical properties of the Galactic bulge., Comment: 16 pages, 11 figures. Accepted for publication in MNRAS

Gravitational wave (GW) standard sirens may resolve the Hubble tension, provided that standard siren inference of $H_0$ is free from systematic biases. However, standard sirens from binary neutron star (BNS) mergers suffer from two sources of systematic bias, one arising from the anisotropy of GW emission, and the other from the anisotropy of electromagnetic (EM) emission from the kilonova. For an observed sample of BNS mergers, the traditional Bayesian approach to debiasing involves the direct computation of the detection likelihood. This is infeasible for large samples of detected BNS merger due to the high dimensionality of the parameter space governing merger detection. In this study, we bypass this computation by fitting the Hubble constant to forward simulations of the observed GW and EM data under a simulation-based inference (SBI) framework using marginal neural ratio estimation. A key innovation of our method is the inclusion of BNS mergers which were only detected in GW, which allows for estimation of the bias introduced by EM anisotropy. Our method corrects for $\sim$90$\%$ of the bias in the inferred value of $H_0$ when telescope follow-up observations of BNS mergers have extensive tiling of the merger localization region, using known telescope sensitivities and assuming a model of kilonova emission. Our SBI-based method thus enables a debiased inference of the Hubble constant of BNS mergers, including both mergers with detected EM counterparts and those without., Comment: Accepted MNRAS, 14 pages, 10 figures

We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 ${\mu}$K/beam for total intensity and 13-23 ${\mu}$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7${\sigma}$ and W51 at 4.0${\sigma}$ with peak frequencies ${\nu}_{AME}$ = (20.0 $\pm$ 1.4) GHz and ${\nu}_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices ${\alpha}_s$ = -0.67 $\pm$ 0.10 in W49 and ${\alpha}_s$ = -0.51 $\pm$ 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency ${\nu}_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study., Comment: 29 pages, 11 figures. Accepted for publication in MNRAS

The role of baryonic physics, star formation and stellar feedback, in shaping the galaxies and their host halos is an evolving topic. The dark matter aspects are illustrated in this work by showing distribution features in a Milky Way sized halo. We focus on the halo morphology, geometry, and profile as well as the phase space distribution using one dark matter only and five hydrodynamical cosmological high-resolution simulations of the same halo with different subgrid prescriptions for the baryonic physics (Kennicut versus multi-freefall star formation and delayed cooling versus mechanical supernovae feedback). If some general properties like the relative halo-galaxy orientation are similar, the modifications of the gravitational potential due to the presence of baryons are found to induce different dark matter distributions (rounder and more concentrated halo). The mass density profile as well as the velocity distribution are modified distinctively according to the specific resulting baryonic distribution highlighting the variability of those properties (e.g inner power index from 1.3 to 1.8, broader speed distribution). The uncertainties on those features are of paramount importance for dark matter phenomenology, particularly when dealing with dark matter dynamics or direct and indirect detection searches. As a consequence, dark matter properties and prospects using cosmological simulations require improvement on baryonic physics description. Modeling such processes is a key issue not only for galaxy formation but also for dark matter investigations.

[Abridged] We use high-resolution $HST$ imaging and $e$-MERLIN 1.5-GHz observations of galaxy cores from the LeMMINGs survey to investigate the relation between optical structural properties and nuclear radio emission for a large sample of galaxies. We perform accurate, multi-component decompositions of new surface brightness profiles extracted from $HST$ images for 163 LeMMINGs galaxies and fit up to six galaxy components (e.g., bulges, discs, AGN, bars, rings, spiral arms, and nuclear star clusters) simultaneously with Sérsic and/or core-Sérsic models. By adding such decomposition data for 10 LeMMINGs galaxies from our past work, the final sample of 173 nearby galaxies (102 Ss, 42 S0s, 23 Es plus 6 Irr) with bulge stellar mass (typically) M_*, bulge ~ 10^6-10^12.5 M_sun, encompasses all optical spectral classes (LINER, Seyfert, ALG and H II). We show that the bulge mass can be significantly overestimated in many galaxies when components such as bars, rings and spirals are not included in the fits. We additionally implement a Monte Carlo method to determine errors on bulge, disc and other fitted structural parameters. Moving (in the opposite direction) across the Hubble sequence, i.e., from the irregular to elliptical galaxies, we confirm that bulges become larger, more prominent and round. Such bulge dominance is associated with a brighter radio core luminosity. We also find that the radio detection fraction increases with bulge mass. At M_*,bulge > 10^11 M_sun, the radio detection fraction is 77%, declining to 24% for M_bulge < 10^10 M_sun. Furthermore, we observe core-Sérsic bulges tend to be systematically round and to possess high radio core luminosities and boxy-distorted or pure elliptical isophotes., 102 pages, 26 figures, 10 tables, accepted for publication in A&A

Weak gravitational lensing is a powerful probe, which is used to constrain the standard cosmological model and its extensions. With the enhanced statistical precision of current and upcoming surveys, high-accuracy predictions for weak lensing statistics are needed to limit the impact of theoretical uncertainties on cosmological parameter constraints. For this purpose, we present a comparison of the theoretical predictions for the non-linear matter and weak lensing power spectra, based on the widely used fitting functions (mead and rev-halofit), emulators (EuclidEmulator, EuclidEmulator2, BaccoEmulator, and CosmicEmulator), and N-body simulations (PKDGRAV3). We consider the forecasted constraints on the ЛCDM and wCDM models from weak lensing for stage III and stage IV surveys. We study the relative bias on the constraints and their dependence on the assumed prescriptions. Assuming a ЛCDM cosmology, we find that the relative agreement on the S8 parameter is between 0.2 and 0.3σ for a stage III-like survey between the above predictors. For a stage IV-like survey the agreement becomes 1.4-3.0σ. In the wCDM scenario, we find broader S8 constraints, and agreements of 0.18-0.26σ and 0.7-1.7σ for stage III and stage IV surveys, respectively. The accuracies of the above predictors therefore appear adequate for stage III surveys, whereas the fitting functions would need improvements for future stage IV surveys. Furthermore, we find that, of the fitting functions, mead provides the best agreement with the emulators. We discuss the implication of these findings for the preparation of future weak lensing surveys, and the relative impact of theoretical uncertainties to other systematics., Monthly Notices of the Royal Astronomical Society, 522 (3), ISSN:0035-8711, ISSN:1365-2966

We present the public data release of the Uchuu-UM galaxy catalogues by applying the UNIVERSEMACHINE algorithm to assign galaxies to the dark matter haloes in the Uchuu N-body cosmological simulation. It includes a variety of baryonic properties for all galaxies down to ∼5 × 108 M⊙ with haloes in a mass range of 1010 < Mhalo/M⊙ < 5 × 1015 up to redshift z = 10. Uchuu-UM includes more than 104 cluster-size haloes in a volume of 8(h−1Gpc)3, reproducing observed stellar mass functions across the redshift range of z = 0−7, galaxy quenched fractions, and clustering statistics at low redshifts. Compared to the previous largest UM catalogue, the Uchuu-UM catalogue includes significantly more massive galaxies hosted by large-mass dark matter haloes. Overall, the number density profile of galaxies in dark matter haloes follows the dark matter profile, with the profile becoming steeper around the splashback radius and flattening at larger radii. The number density profile of galaxies tends to be steeper for larger stellar masses and depends on the colour of galaxies, with red galaxies having steeper slopes at all radii than blue galaxies. The quenched fraction exhibits a strong dependence on the stellar mass and increases towards the inner regions of clusters. The publicly available Uchuu-UM galaxy catalogue presented here can serve to model ongoing and upcoming large galaxy surveys. © 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., DN and HA acknowledge support from Yale University and the Yale Center for Research Computing facilities and staff. HA also acknowledges support from the Zuckerman Postdoctoral Scholars Program while completing this work. T.I. has been supported by IAAR Research Support Program at Chiba University Japan, MEXT/JSPS KAKENHI (grant numbers JP19KK0344, JP21F51024, and JP21H01122), MEXT as ‘Program for Promoting Researches on the Supercomputer Fugaku’ (JPMXP1020200109), and JICFuS. FP, AK, EP, JR thank the support of the Spanish Ministry of Science and Innovation funding grant PGC2018-101931-B-I00., With funding from the Spanish government through the "Severo Ochoa Centre of Excellence" accreditation (CEX2021-001131-S).

Recent attempts to detect [OIII] 88$μ$m emission from super-early ($z>10$) galaxy candidates observed by JWST have been unsuccessful. By using zoom-in simulations, we show that these galaxies are faint, and mostly fall below the local metal-poor $\rm [OIII]-SFR$ relation as a result of their low ionization parameter, $U_{\rm ion}\lesssim 10^{-3}$. Such low $U_{\rm ion}$ values are found in galaxies that are in an early assembly stage, and whose stars are still embedded in high-density natal clouds. However, the most luminous galaxy in our sample ($\rm{log}[L_{\rm{[OIII]}}/L_\odot] = 8.4$, $U_{\rm ion} \approx 0.1$) could be detected by ALMA in only $2.8$ hrs., 6 pages, 3 figures, Submitted to MNRAS (comments are welcome)

The Fix and Pair techniques were designed to generate simulations with reduced variance in the 2-point statistics by modifying the Initial Conditions (ICs). In this paper we show that this technique is also valid when the initial conditions have local Primordial non-Gaussianities (PNG), parametrised by $f_{\rm NL}$, without biasing the 2-point statistics but reducing significantly their variance. We show how to quantitatively use these techniques to test the accuracy of galaxy/halo clustering models down to a much reduced uncertainty and we apply them to test the standard model for halo clustering in the presence of PNG. Additionally, we show that by Matching the stochastic part of the ICs for two different cosmologies (Gaussian and non-Gaussian) we obtain a large correlation between the (2-point) statistics that can explicitly be used to further reduce the uncertainty of the model testing. For our reference analysis ($f_{\rm NL}=100$, $V=1 [h^{-1}{\rm Gpc}]^3$, $n= 2.5\times 10^{-4}[h^{-1}{\rm Mpc}]^{-3}$, $b=2.32$), we obtain an uncertainty of $\sigma(f_{\rm NL})=60$ with a standard simulation, whereas using Fixed [Fixed-Paired] initial conditions it reduces to $\sigma(f_{\rm NL})=12$ [$\sigma(f_{\rm NL})=12$]. When also Matching the ICs we obtain $\sigma (f_{\rm NL})=18$ for the standard case, and $\sigma (f_{\rm NL})=8$ [$\sigma (f_{\rm NL})=7$] for Fixed [Fixed-Paired]. The combination of the Fix, Pair and Match techniques can be used in the context of PNG to create simulations with an effective volume incremented by a factor $\sim 70$ at given computational resources., Comment: 14+6 pages, 8+8 figures. Accepted in MNRAS

We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-alpha forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-alpha forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-alpha forest, such as flat-bottom or double-dip absorption profiles., Version accepted for publication in MNRAS, some discussions added, data added to Fig. 8. 24 pages, 15 figures. Images and movies may be viewed at https://www.nottingham.ac.uk/astronomy/sherwood-relics/

Fast radio bursts (FRBs) are millisecond-duration transients with large dispersion measures. The origin of FRBs is still mysterious. One of the methods to comprehend FRB origin is to probe the physical environments of FRB host galaxies. Mapping molecular-gas kinematics in FRB host galaxies is critical because it results in star formation that is likely connected to the birth of FRB progenitors. However, most previous works of FRB host galaxies have focused on its stellar component. Therefore, we, for the first time, report the molecular gas kinematics in the host galaxy of the non-repeating FRB 180924B at $z= 0.3216$. Two velocity components of the CO (3-2) emission line are detected in its host galaxy with the Atacama Large Millimeter/submillimeter Array (ALMA): the peak of one component ($-155.40$ km s$^{-1}$) is near the centre of the host galaxy, and another ($-7.76$ km s$^{-1}$) is near the FRB position. The CO (3-2) spectrum shows asymmetric profiles with A$_{\rm peak}$ $=2.03\pm 0.39$, where A$_{\rm peak}$ is the peak flux density ratio between the two velocity components. The CO (3-2) velocity map also indicates an asymmetric velocity gradient from $-180$ km s$^{-1}$ to 8 km s$^{-1}$. These results indicate a disturbed kinetic structure of molecular gas in the host galaxy. Such disturbed kinetic structures are reported for repeating FRB host galaxies using HI emission lines in previous works. Our finding indicates that non-repeating and repeating FRBs could commonly appear in disturbed kinetic environments, suggesting a possible link between the gas kinematics and FRB progenitors., Comment: 5 pages, 4 figures, Accepted for publication in MNRAS, https://www.youtube.com/watch?v=CldxLE7Pdwk&t=1s

We present a new set of cosmological zoom-in simulations of a MW-like galaxy which for the first time include elastic velocity-dependent self interacting dark matter (SIDM) and IllustrisTNG physics. With these simulations we investigate the interaction between SIDM and baryons and its effects on the galaxy evolution process. We also introduce a novel set of modified DMO simulations which can reasonably replicate the effects of fully realized hydrodynamics on the DM halo while simplifying the analysis and lowering the computational cost. We find that baryons change the thermal structure of the central region of the halo to a greater extent than the SIDM scatterings for MW-like galaxies. Additionally, we find that the new thermal structure of the MW-like halo causes SIDM to create cuspier central densities rather than cores because the SIDM scatterings remove the thermal support by transferring heat away from the center of the galaxy. We find that this effect, caused by baryon contraction, begins to affect galaxies with a stellar mass of $10^8$ M$_\odot$ and increases in strength to the MW-mass scale. This implies that any simulations used to constrain the SIDM cross sections for galaxies with stellar masses between $10^8$ and at least $10^{11}$ M$_\odot$ will require baryons to make accurate predictions., Comment: 14 pages, 11 figures

Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS., Comment: 11 pages, 9 figures. Matching to the published version on MNRAS

Non-thermal emission from relativistic Cosmic Ray (CR) electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of structure formation shocks. Emission caused by CR protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We introduce CRESCENDO: Cosmic Ray Evolution with SpeCtral Electrons aND prOtons, an efficient on-the-fly Fokker-Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR (re-)acceleration at intra-cluster shocks, based on results of recent PIC simulations, adiabatic changes and radiative losses of electrons. We show its performance in test cases as well as idealized galaxy cluster (GC) simulations. We apply the model to an idealized GC merger following best-fit parameters for CIZA J2242.4+5301-1 and study CR injection, radio relic morphology, spectral steepening and synchrotron emission., Comment: 25 pages (21 main paper + 4 appendix), 17 figures. Accepted for publication in MNRAS

We employ several galaxy formation models, in particular, L-GALAXIES, IllustrisTNG, and EAGLE, as well as observational samples from SDSS and DESI, to investigate galactic conformity, the observed large-scale correlation between the star-formation properties of central (primary) galaxies and those of their neighbours. To analyse the models and observations uniformly, we introduce CENSAT, a new algorithm to define whether a galaxy is a central or a satellite system based on an isolation criterion. We find that the conformity signal is present, up to at least 5 Mpc from the centres of low- and intermediate-mass centrals in the latest version of L-GALAXIES (Ayromlou et al. 2021), IllustrisTNG, and EAGLE, as well as in SDSS and DESI observational samples. In comparison, the conformity signal is substantially weaker in an older version of L-GALAXIES (Henriques et al. 2020). One of the main differences between this older model and the other models is its neglect of ram-pressure stripping of the gas reservoirs of galaxies except within the boundaries of massive cluster haloes. Our observational comparisons demonstrate that this difference significantly affects the observed large-scale conformity signal. Furthermore, by examining the contribution of backsplash, fly-by, central, and satellite galaxies to the conformity signal, we show that much, but not all, of it arises from primary galaxies near massive systems. Remaining tensions between the models and observations may be solved by modifying the physical prescriptions for how feedback processes affect the distribution and kinematics of gas and the environment around galaxies out to scales of several Megaparsecs., Comment: Submitted to MNRAS

The next-generation cosmic microwave background (CMB) satellite missions are expected to provide robust constraints on a wide range of cosmological parameters with unprecedented precision. But these constraints on the parameters could weaken if we do not attribute dark energy to a cosmological constant. The cosmological models involving interaction between dark energy and dark matter can give rise to comparable energy densities at the present epoch, thereby alleviating the so-called cosmic coincidence problem. In the present paper, we perform a forecast analysis to test the ability of the future generation high-sensitive CMB, and baryon acoustic oscillation (BAO) experiments to constrain phenomenological interacting dark energy models. We consider cosmic variance limited future CMB experiment Probe of Inflation and Cosmic Origins (PICO) along with BAO information from the Dark Energy Spectroscopic Instrument (DESI), to constrain the parameters of the interacting dark sector. Based on the stability of the cosmological perturbations, we consider two possibilities for the interaction scenario. We investigate the impact of both coupling constant and equation of state parameter of dark energy on CMB temperature power spectrum, matter power spectrum, and fσ8. We have used simulated temperature and polarization data from PICO within the multipole ranges (ℓ = 2–4000), and as expected, we do see PICO alone produces better constraints than Planck on the ΛCDM parameters. With the integration of PICO and DESI missions, we observe a significant improvement in the constraints on several cosmological parameters, especially the equation of state parameter of dark energy. However, we note that additional data are required to constrain a small positive coupling constant.

As one of Stage IV space-based telescopes, China Space Station Telescope (CSST) can perform photometric and spectroscopic surveys simultaneously to efficiently explore the Universe in extreme precision. In this work, we investigate several powerful CSST cosmological probes, including cosmic shear, galaxy-galaxy lensing, photometric and spectroscopic galaxy clustering, and number counts of galaxy clusters, and study the capability of these probes by forecasting the results of joint constraints on the cosmological parameters. By referring to real observational results, we generate mock data and estimate the measured errors based on CSST observational and instrumental designs. To study the systematical effects on the results, we also consider a number of systematics in CSST photometric and spectroscopic surveys, such as the intrinsic alignment, shear calibration uncertainties, photometric redshift uncertainties, galaxy bias, non-linear effects, instrumental effects, etc. The Fisher matrix method is used to derive the constraint results from individual or joint surveys on the cosmological and systematical parameters. We find that the joint constraints by including all these CSST cosmological probes can significantly improve the results from current observations by one order of magnitude at least, which gives $\Omega_m$ and $\sigma_8$ $, Comment: 17 pages, 12 figures, 3 tables. Accepted for publication in MNRAS

Entropy is an advantageous diagnostics to study the thermodynamic history of the intracluster plasma of galaxy clusters. We present the entropy profile of the Abell 2244 galaxy cluster derived both exclusively using X-ray data from the low-background Swift XRT telescope and also using Planck y data. The entropy profile derivation using X-rays only is robust at least to the virial radius because the cluster brightness is large compared to the X-ray background at low energies, temperature is strongly bounded by the lack of cluster X-ray photons at energies kT>3 keV, and the XRT background is low, stable and understood. In the observed solid angle, about one quadrant, the entropy radial profile deviates from a power-law at the virial radius, mainly because of a sharp drop of the cluster temperature. This bending of the entropy profile is confirmed when X-ray spectral information is replaced by the Compton map. Clumping and non-thermal pressure support are insufficient to restore a power law entropy profile because they are bound to be small by: a) the agreement between mass estimates from different tracers (gas and galaxies), b) the agreement between entropy profile determinations based on combinations of observables with different sensitivities and systematics, and c) the low value of clumping as estimated using the azimuthal scatter and the gas fraction. Based on numerical simulations, ion-electron equilibration is also insufficient to restore a linear entropy profile. Therefore, the bending of the entropy profiles seems to be robustly derived and witnesses the teoretically-predicted decrease in the inflow through the virial boundary., MNRAS, in press

We analyze the physical properties of gas in the circumgalactic medium (CGM) of 132 Milky Way (MW)-like central galaxies at $z=0$ from the cosmological magneto-hydrodynamical simulation TNG50, part of the IllustrisTNG project. The properties and abundance of CGM gas across the sample are diverse, and the fractional budgets of different phases (cold, warm, and hot), as well as neutral HI mass and metal mass, vary considerably. Over our stellar mass range of $10^{10.5} < M_\star / \rm{M}_\odot < 10^{10.9}$, radial profiles of gas physical properties from $0.15 < R\rm{ / R_{\rm 200c}} < 1.0$ reveal great CGM structural complexity, with significant variations both at fixed distance around individual galaxies, and across different galaxies. CGM gas is multi-phase: the distributions of density, temperature and entropy are all multimodal, while metallicity and thermal pressure distributions are unimodal; all are broad. We present predictions for magnetic fields in MW-like halos: a median field strength of $|B|\sim\,1μ$G in the inner halo decreases rapidly at larger distance, while magnetic pressure dominates over thermal pressure only within $\sim0.2 \times \rm{R_{200c}}$. Virial temperature gas at $\sim 10^6\,$K coexists with a sub-dominant cool, $< 10^5\,$K component in approximate pressure equilibrium. Finally, the physical properties of the CGM are tightly connected to the galactic star formation rate, in turn dependent on feedback from supermassive black holes (SMBHs). In TNG50, we find that energy from SMBH-driven kinetic winds generates high-velocity outflows ($\gtrsim 500-2000$ km/s), heats gas to super-virial temperatures ($> 10^{6.5-7}$ K), and regulates the net balance of inflows versus outflows in otherwise quasi-static gaseous halos., Minor typos fixed. Additional visuals at www.tng-project.org/ramesh22!

Measurements of the one-point probability distribution function and higher-order moments (variance, skewness, and kurtosis) of the high-redshift 21 cm fluctuations are among the most direct statistical probes of the non-Gaussian nature of structure formation and evolution during reionization. However, contamination from astrophysical foregrounds and instrument systematics pose significant challenges in measuring these statistics in real observations. In this work, we use forward modelling to investigate the feasibility of measuring 21 cm one-point statistics through a foreground avoidance strategy. Leveraging the characteristic wedge-shape of the foregrounds in k-space, we apply a wedge-cut filter that removes the foreground contaminated modes from a mock data set based on the Hydrogen Epoch of Reionization Array (HERA) instrument, and measure the one-point statistics from the image-space representation of the remaining non-contaminated modes. We experiment with varying degrees of wedge-cutting over different frequency bandwidths and find that the centre of the band is the least susceptible to bias from wedge-cutting. Based on this finding, we introduce a rolling filter method that allows reconstruction of an optimal wedge-cut 21~cm intensity map over the full bandwidth using outputs from wedge-cutting over multiple sub-bands. We perform Monte Carlo simulations to show that HERA should be able to measure the rise in skewness and kurtosis near the end of reionization with the rolling wedge-cut method if foreground leakage from the Fourier transform window function can be controlled., 13 pages, 10 figures, Accepted to MNRAS

In this paper, we report an interesting kinematic phenomenon around the halos' edge related to the splashback radius. After the shell-crossing, cosmic flow exhibits various rotational morphologies via stream-mixing. Vorticity is generated in a particular way that coincides with the large-scale structure. Notably, one specific flow morphology, which is spiraling inward and compressing in the third direction, concentrates around halos. A detailed examination that reveals a sharp change in the logarithmic derivative of its volume fraction, coincides with the location of the splashback radius defined as the outermost caustic structure. Such a feature encodes valuable phase space information and provides a new perspective on understanding the dynamical evolution of halos. As a volume-weighted quantity, the profile of flow morphology is purely kinematic. And unlike other related studies, the rotational flow morphologies capture the anisotropic phase structure in the multi-stream region., Comment: Accepted by MNRAS

Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptions. The complexity of these models makes it difficult to assess the validity of any approximation and to quantify uncertainties on the model results. This impedes scientific rigour, in particular, when comparing models to observations, or when using their results as input for other models. We present a probabilistic numerical approach to address these issues by treating radiative transfer as a Bayesian linear regression problem. This allows us to model uncertainties on the input and output of the model with the variances of the associated probability distributions. Furthermore, this approach naturally allows us to create reduced-order radiative transfer models with a quantifiable accuracy. These are approximate solutions to exact radiative transfer models, in contrast to the exact solutions to approximate models that are often used. As a first demonstration, we derive a probabilistic version of the method of characteristics, a commonly-used technique to solve radiative transfer problems., 16 pages, 1 figure, accepted to Monthly Notices of the Royal Astronomical Society (MNRAS)

This paper reports our discovery of the most massive supercluster, termed the King Ghidorah Supercluster (KGSc), at $z=0.50-0.64$ in the Third Public Data Release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP PDR3) over 690 deg$^2$, as well as an initial result for a galaxy and dark matter mapping. The primary structure of the KGSc comprises triple broad weak-lensing (WL) peaks over 70 comoving Mpc. Such extensive WL detection at $z>0.5$ can only currently be achieved using the wide-field high-quality images produced by the HSC-SSP. The structure is also contiguous with multiple large-scale structures across a $\sim400$ comoving Mpc scale. The entire field has a notable overdensity ($\delta=14.7\pm4.5$) of red-sequence clusters. Additionally, large-scale underdensities can be found in the foreground along the line of sight. We confirmed the overdensities in stellar mass and dark matter distributions to be tightly coupled and estimated the total mass of the main structure to be $1\times10^{16}$ solar masses, according to the mock data analyses based on large-volume cosmological simulations. Further, upcoming wide-field multi-object spectrographs such as the Subaru Prime Focus Spectrograph may aid in providing additional insights into distant superclusters beyond the 100 Mpc scale., Comment: 5 pages, 5 figures, accepted for publication in MNRAS Letter

We present a Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. As our approach features a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the final density and velocity fields as natural byproducts. We implement this field-based physical inference approach by adapting the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm, which explores the high-dimensional posterior through the use of Hamiltonian Monte Carlo sampling. We test the self-consistency of the method using random sets of mock tracers, and assess its accuracy in a more complex scenario where peculiar velocity tracers are non-linearly evolved mock haloes. We find that our framework self-consistently infers the initial conditions, density and velocity fields, and shows some robustness to model mis-specification. As compared to the state-of-the-art approach of constrained Gaussian random fields/Wiener filtering, our method produces more accurate final density and velocity field reconstructions. It also allows us to constrain the initial conditions by peculiar velocity observations, complementing in this aspect previous field-based approaches based on other cosmological observables., Comment: 23 pages, 15 figures. Accepted for publication in MNRAS

The theoretical basis of dark energy remains unknown and could signify a need to modify the laws of gravity on cosmological scales. In this study we investigate how the clustering and motions of galaxies can be used as probes of modified gravity theories, using galaxy and direct peculiar velocity auto- and cross-correlation functions. We measure and fit these correlation functions in simulations of $\Lambda$CDM, DGP, and $f(R)$ cosmologies and, by extracting the characteristic parameters of each model, we show that these theories can be distinguished from General Relativity using these measurements. We present forecasts showing that with sufficiently large data samples, this analysis technique is a competitive probe that can help place limits on allowed deviations from GR. For example, a peculiar velocity survey reaching to $z=0.5$ with $20\%$ distance accuracy would constrain model parameters to 3-$\sigma$ confidence limits $\log_{10}|f_{R0}| < -6.45$ for $f(R)$ gravity and $r_c > 2.88 \, c/H_0$ for nDGP, assuming a fiducial GR model., Comment: 14 pages, 12 figures, 3 tables, accepted for publication in MNRAS

Ultralight axion-like particles $m_a \sim 10^{-22}$ eV, or Fuzzy Dark Matter (FDM), behave comparably to cold dark matter (CDM) on cosmological scales and exhibit a kpc-size de Broglie wavelength capable of alleviating established (sub-)galactic-scale problems of CDM. Substructures inside an FDM halo incur gravitational potential perturbations, resulting in stellar heating sufficient to account for the Galactic disc thickening over a Hubble time, as first demonstrated by Church et al. We present a more sophisticated treatment that incorporates the full baryon and dark matter distributions of the Milky Way and adopts stellar disc kinematics inferred from recent Gaia, APOGEE, and LAMOST surveys. Ubiquitous density granulation and subhalo passages respectively drive inner disc thickening and flaring of the outer disc, resulting in an observationally consistent `U-shaped' disc vertical velocity dispersion profile with the global minimum located near the solar radius. The observed age-velocity dispersion relation in the solar vicinity can be explained by the FDM-substructure-induced heating and places an exclusion bound $m_a \gtrsim 0.4\times10^{-22}$ eV. We assess non-trivial uncertainties in the empirical core-halo relation, FDM subhalo mass function and tidal stripping, and stellar heating estimate. The mass range $m_a\simeq 0.5-0.7\times10^{-22}$ eV favoured by the observed thick disc kinematics is in tension with several exclusion bounds inferred from dwarf density profiles, stellar streams, and Milky Way satellite populations, which could be significantly relaxed due to the aforesaid uncertainties. Additionally, strongly anisotropic heating could help explain the formation of ultra-thin disc galaxies., Comment: 19 pages, 9 figures; Accepted for publication in MNRAS. Matches published version

In this paper we present an extension to the $\texttt{matryoshka}$ suite of neural-network-based emulators. The new editions have been developed to accelerate EFTofLSS analyses of galaxy power spectrum multipoles in redshift space. They are collectively referred to as the $\texttt{EFTEMU}$. We test the $\texttt{EFTEMU}$ at the power spectrum level and achieve a prediction accuracy of better than 1\% with BOSS-like bias parameters and counterterms on scales $0.001\ h\ \mathrm{Mpc}^{-1} \leq k \leq 0.19\ h\ \mathrm{Mpc}^{-1}$. We also run a series of mock full shape analyses to test the performance of the $\texttt{EFTEMU}$ when carrying out parameter inference. Through these mock analyses we verify that the $\texttt{EFTEMU}$ recovers the true cosmology within $1\sigma$ at several redshifts ($z=[0.38,0.51,0.61]$), and with several noise levels (the most stringent of which is Gaussian covariance associated with a volume of $5000^3 \ \mathrm{Mpc}^3 \ h^{-3}$). We compare the mock inference results from the $\texttt{EFTEMU}$ to those obtained with a fully analytic EFTofLSS model and again find no significant bias, whilst speeding up the inference by three orders of magnitude. The $\texttt{EFTEMU}$ is publicly available as part of the $\texttt{matryoshka}$ $\texttt{Python}$ package., Comment: MNRAS accepted version. 11 pages, 8 figures, 3 tables. Code available at https://github.com/JDonaldM/Matryoshka

We present constraints on primordial B modes from large angular scale cosmic microwave background polarisation anisotropies measured with the Planck satellite. To remove Galactic polarised foregrounds, we use a Bayesian parametric component separation method, modelling synchrotron radiation as a power law and thermal dust emission as a modified blackbody. This method propagates uncertainties from the foreground cleaning into the noise covariance matrices of the maps. We construct two likelihoods: (i) a semi-analytical cross-spectrum-based likelihood-approximation scheme (momento) and (ii) an exact polarisation-only pixel-based likelihood (pixlike). Since momento is based on cross-spectra it is statistically less powerful than pixlike, but is less sensitive to systematic errors correlated across frequencies. Both likelihoods give a tensor-to-scalar ratio, r, that is consistent with zero from low multipole (2, 10 pages, 7 figures

We study systematics associated with estimating simple stellar population (SSP) parameters -- age, metallicity [M/H], $\alpha$-enhancement [$\alpha$/Fe] and IMF shape -- and associated $M_*/L$ gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs) and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H$_\beta$, $\langle$Fe$\rangle$, TiO$_{\rm 2SDSS}$ and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius $R_e$: (a) ages are younger and [$\alpha$/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about $R_e/2$; (c) this makes $M_*/L$ about $2\times$ larger in the central regions than beyond $R_e/2$. While the models of Conroy et al. (2018) return similar [M/H] and [$\alpha$/Fe] profiles, the age and (hence) $M_*/L$ profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs and S0s: younger ages and less enhanced [$\alpha$/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred $M_*/L$ gradients significantly., Comment: 22 pages, 23 figures, accepted for publication in MNRAS

We conduct the first case study towards developing optimal foreground mitigation strategies for neutral hydrogen (HI) intensity mapping using radio interferometers at low redshifts. A pipeline for simulation, foreground mitigation and power spectrum estimation is built, which can be used for ongoing and future surveys using MeerKAT and Square Kilometre Array Observatory (SKAO). It simulates realistic sky signals to generate visibility data given instrument and observation specifications, which is subsequently used to perform foreground mitigation and power spectrum estimation. A quadratic estimator formalism is developed to estimate the temperature power spectrum in visibility space. Using MeerKAT telescope specifications for observations in the redshift range z~0.25-0.30 corresponding to the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey, we present a case study where we compare different approaches of foreground mitigation. We find that component separation in visibility space provides a more accurate estimation of HI clustering comparing to foreground avoidance, with the uncertainties being 30 per cent smaller. Power spectrum estimation from image is found to be less robust with larger bias and more information loss when compared to estimation in visibility. We conclude that for z~0.25-0.30, the MIGHTEE survey will be capable of measuring the HI power spectrum from k~0.5 Mpc$^{-1}$ to k~10 Mpc$^{-1}$ with high accuracy. We are the first to show that, at low redshift, component separation in visibility space suppresses foreground contamination at large line-of-sight scales, allowing measurement of HI power spectrum closer to the foreground wedge, crucial for data analysis towards future detections., 20 pages, 19 figures. Updated to match the published version in MNRAS

We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a $\sim 2 \sigma$ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c > 10^{13.2} h^{-1} M_{\odot}$ (95\% C.L.). The best-fit baryonic suppression is $\sim 5\%$ at $k=1.0 {\rm Mpc}\ h^{-1}$ and $\sim 15\%$ at $k=5.0 {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments., Comment: 20 pages, 10 figures. DES Collaboration, Year-3 analysis