Search

Motivated by the recent release of new results from five different pulsar timing array (PTA) experiments claiming to have found compelling evidence for primordial gravitational waves (GW) at nano-Hz frequencies, we consider the prospects of generating such a signal from inflationary blue-tilted tensor power spectrum in a specific dark matter (DM) scenario dubbed as \textit{Miracle-less WIMP}. While \textit{Miracle-less WIMP}, due to insufficient interaction rate with the standard bath gets thermally overproduced, inflationary blue-tilted gravitational waves (BGW) leads to conflict with cosmological observations if solely responsible for the PTA events. Both these problems are circumvented with late entropy dilution bringing DM abundance within limits while creating a doubly peaked feature of BGW. The blue-tilted part of one of these peaks can fit NANOGrav 15 yr data at $1\sigma$ level. The particle physics setup used here for illustration namely, the gauged $U(1)_{B-L}$ model, naturally leads to \textit{Miracle-less WIMP} and long-lived diluter for entropy dilution while also having GW complementarity due to cosmic strings., Comment: 6 pages, 1 captioned figure

We explore possible signatures of the interaction between dark matter (DM) and massive neutrinos during the post-reionization epoch. Using both Fisher matrix forecast analysis and Markov Chain Monte-Carlo (MCMC) simulation, we conduct a thorough investigation of the constraints and imprints of the scenario on the upcoming post-reionization and galaxy surveys. Our investigation focuses on two key parameters: the strength of the DM-massive neutrino interaction ($u$) and the total neutrino mass ($M_{\rm tot}$), on top of the usual 6 cosmological parameters. We utilize future 21-cm intensity mapping, galaxy clustering as well as cosmic shear observations in order to investigate the possible constraints of these parameters in the future observations: Square Kilometre Array (SKA1 and SKA2) and Euclid, taking both conservative and realistic approaches. All these missions show promise in constraining both the parameters $u$ and $M_{\rm tot}$ by 3-4 orders compared to the current constraints from Planck18 (SKA2 performing the best among them). Although none of the missions help much in addressing the $H_0$ and $\sigma_8$ tensions for the present scenario, SKA2 constrains them better in conservative approach. We further perform a brief investigation of the prospects of some of the next generation Cosmic Microwave Background (CMB) missions in combinations with LSS experiments in improving the constraints. Our analysis reveals that both SKA2 and CMB-S4 + Euclid + SKA1 IM2 combination will put the strongest bounds on the model parameters., Comment: 13 pages, 10 figures

Primordial black holes may arise through ultra slow-roll inflation. In this work we study a toy model of ultra slow-roll inflation with a shallow dip. The ultra slow-roll stage enhances the curvature perturbations and thus the primordial scalar power spectrum. We analyze the features of the power spectrum numerically and analytically, and then give a rough estimate of the lower and upper bound of the enhancement. These large perturbations also produce second order gravitational waves, which are in the scope of future observations., 11pages, 7 figures

Significant evidence for a gravitational-wave background was reported by several pulsar-timing-array collaborations. By assuming that this signal is interpreted by the scalar-induced gravitational waves, we study physical implications of the observed signal for the nature of primordial curvature perturbations and primordial black holes. In particular, we explore the effects of primordial non-Gaussianity on the inferences of model parameters, and obtain the parameter region allowed by the observed signal, i.e., the primordial scalar spectral amplitude $A_S\sim10^{-2}-1$, the primordial non-Gaussian parameter $-10\lesssim f_{\mathrm{NL}} \lesssim 10$, and the mass of primordial black holes $m_{\mathrm{pbh}}\sim10^{-3}-0.1M_{\odot}$. We find that the non-Gaussianity suppressing the abundance of primordial black holes is preferred by the observed signal. We show that the anisotropies of scalar-induced gravitational waves are a powerful probe for measurements of the non-Gaussian parameter $f_{\mathrm{NL}}$, and conduct a complete analysis of the angular power spectrum in the nano-Hertz band. We expect that the Square Kilometre Array project has potentials to measure such anisotropies., 19 pages, 5 figures

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

Decaying dark matter models generically modify the equation of state around the time of dark matter decay, and this in turn modifies the expansion rate of the Universe through the Friedmann equation. Thus, a priori, these models could solve or alleviate the Hubble tension, and depending on the lifetime of the dark matter, they can be classified as belonging to either the early- or late-time solutions. Moreover, decaying dark matter models can often be realized in particle physics models relatively easily. However, the implementations of these models in Einstein--Boltzmann solver codes are non-trivial, so not all incarnations have been tested. It is well known that models with very late decay of dark matter do not alleviate the Hubble tension, and in fact, cosmological data puts severe constraints on the lifetime of such dark matter scenarios. However, models in which a fraction of the dark matter decays to dark radiation at early times hold the possibility of modifying the effective equation of state around matter-radiation equality without affecting late-time cosmology. This scenario is therefore a simple realization of a possible early-time solution to the Hubble tension, and cosmological parameter estimation with current data in these models yields a value of $H_0 = 68.73^{+0.81}_{-1.3}$ at $68\%$ C.I.. This still leads to a $2.7\sigma$ Gaussian tension with the representative local value of $H_0 = 73.2 \pm 1.3$ km s$^{-1}$ Mpc$^{-1}$. Additional work is, however, required to test more complex decay scenarios, which could potentially prefer higher values of $H_0$ and provide a better solution to the Hubble tension., Comment: Invited chapter for the edited book Hubble Constant Tension (Eds. E. Di Valentino and D. Brout, Springer Singapore, expected in 2024)

We investigate the evolution of gravitational waves through discontinuous evolution (transition) of the Hubble expansion rate $H(z)$ at a sudden cosmological singularity, which may be due to a transition of the value of the gravitational constant. We find the evolution of the scale factor and the gravitational wave waveform through the singularity by imposing the proper boundary conditions. We also use existing cosmological data and mock data of future gravitational wave experiments (the ET) to impose current and anticipated constraints on the magnitude of such a transition. We show that mock data of the Einstein Telescope can reduce the uncertainties by up to a factor of three depending on the cosmological parameter considered., 40 pages, 15 Figures, 6 Tables. Published in Universe (Open Access: https://www.mdpi.com/2218-1997/9/7/317)!

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

The eROSITA Final Equatorial Depth Survey (eFEDS), with a sky area of 140 square degrees with depth equivalent to the equatorial patch of the final eROSITA all-sky survey, represents the largest continuous non-full-sky X-ray fields to-date, making it the premier data set for measuring the angular power spectrum. In this work, we measure the X-ray angular power spectrum of galaxy clusters and groups in the eFEDS field. We show that the measured power spectrum is consistent with past observations, including the ROSAT All Sky Survey, and the Chandra COSMOS and Bootes fields. The predictions of cluster gas halo model that is calibrated from Chandra observations is also consistent with the eFEDS power spectrum. While the eFEDS does not have large enough sky coverage to provide meaningful cosmological constraints, we predict that the X-ray power spectrum from the cycle 4 of the eROSITA All-Sky Survey (eRASS4) will provide constraints on $\Omega_M$ and $\sigma_8$ at the 10% level., Comment: 10 pages, 4 Figures, accepted for publication in MNRAS. Updated arxiv version matched to the published version. Major change: tensions at small-scale power are resolved after point source mask correction

Precision analysis of galaxy-galaxy strong gravitational lensing images provides a unique way of characterizing small-scale dark matter halos, and could allow us to uncover the fundamental properties of dark matter's constituents. Recently, gravitational imaging techniques made it possible to detect a few heavy subhalos. However, gravitational lenses contain numerous subhalos and line-of-sight halos, whose subtle imprint is extremely difficult to detect individually. Existing methods for marginalizing over this large population of sub-threshold perturbers to infer population-level parameters are typically computationally expensive, or require compressing observations into hand-crafted summary statistics, such as a power spectrum of residuals. Here, we present the first analysis pipeline to combine parametric lensing models and a recently-developed neural simulation-based inference technique called truncated marginal neural ratio estimation (TMNRE) to constrain the warm dark matter halo mass function cutoff scale directly from multiple lensing images. Through a proof-of-concept application to simulated data, we show that our approach enables empirically testable inference of the dark matter cutoff mass through marginalization over a large population of realistic perturbers that would be undetectable on their own, and over lens and source parameters uncertainties. To obtain our results, we combine the signal contained in a set of images with Hubble Space Telescope resolution. Our results suggest that TMNRE can be a powerful approach to put tight constraints on the mass of warm dark matter in the multi-keV regime, which will be relevant both for existing lensing data and in the large sample of lenses that will be delivered by near-future telescopes., 15 pages, 11 figures; Updated to match version published in MNRAS

We study evolution of single subhaloes with their masses of $\sim10^9 M_\odot$ in a Milky-Way-sized host halo for self-interacting dark matter (SIDM) models. We perform dark-matter-only N-body simulations of dynamical evolution of individual subhaloes orbiting its host by varying self-scattering cross sections (including a velocity-dependent scenario), subhalo orbits, and internal properties of the subhalo. We calibrate a gravothermal fluid model to predict time evolution in spherical mass density profiles of isolated SIDM haloes with the simulations. We find that tidal effects of SIDM subhaloes can be described with a framework developed for the case of collision-less cold dark matter (CDM), but a shorter typical time scale for the mass loss due to tidal stripping is required to explain our SIDM simulation results. As long as the cross section is less than $\sim10\, \mathrm{cm}^2/\mathrm{g}$ and initial states of subhaloes are set within a $2\sigma$-level scatter at redshifts of $\sim2$ predicted by the standard $\Lambda$CDM cosmology, our simulations do not exhibit a prominent feature of gravothermal collapse in the subhalo central density for 10 Gyr. We develop a semi-analytic model of SIDM subhaloes in a time-evolving density core of the host with tidal stripping and self-scattering ram pressure effects. Our semi-analytic approach provides a simple, efficient and physically-intuitive prediction of SIDM subhaloes, but further improvements are needed to account for baryonic effects in the host and the gravothermal instability accelerated by tidal stripping effects., Comment: 19 pages, 12 figures, 1 table. Accepted for publication in MNRAS. Note that results in the previous manuscript (v1) were affected by a bug in our SIDM implementation. We fixed the bug and updated all results accordingly

This study explores the gravitational lensing effects of supermassive black holes (SMBHs) in galaxy clusters. While the presence of central SMBHs in galaxies is firmly established, recent work from high-resolution simulations predict the existence of an additional population of wandering SMBHs. Though the masses of these SMBHs are a minor perturbation on the larger scale and individual galaxy scale dark matter components in the cluster, they can impact statistical lensing properties and individual lensed image configurations. Probing for these potentially observable signatures, we find that SMBHs imprint detectable signatures in rare, higher-order strong lensing image configurations although they do not manifest any statistically significant detectable evidence in either the magnification distribution or the integrated shear profile. Investigating specific lensed image geometries, we report that a massive, near point-like, potential of an SMBH causes the following detectable effects: (i) image splitting leading to the generation of extra images; (ii) positional and magnification asymmetries in multiply imaged systems; and (iii) the apparent disappearance of a lensed counter-image. Of these, image splitting inside the cluster tangential critical curve, is the most prevalent notable observational signature. We demonstrate these possibilities in two cases of observed giant arcs in $SGAS\,J003341.5+024217$ and $RX\,J1347.5-1145$, wherein specific image configurations seen can be reproduced with SMBHs. Future observations with high-resolution instrumentation (e.g. MAVIS-Very Large Telescope, MICADO-Extremely Large Telescope, and the upgraded ngVLA, along with data from the \textit{Euclid} \& \textit{Nancy Grace Roman} Space Telescopes and the Rubin LSST Observatory are likely to allow us to probe these unique yet rare SMBHs lensing signatures., 12 pages, 12 figures, submitted to MNRAS. For associated online animation, see https://sites.google.com/view/guillaume-mahler-astronomer/paper-animation

Constraints on cosmological parameters are often distilled from sky surveys by fitting templates to summary statistics of the data that are motivated by a fiducial cosmological model. However, recent work has shown how to estimate the distance scale using templates that are more generic: the basis functions used are not explicitly tied to any one cosmological model. We describe a Bayesian framework for (i) determining how many basis functions to use and (ii) comparing one basis set with another. Our formulation provides intuition into how (a) one's degree of belief in different basis sets, (b) the fact that the choice of priors depends on basis set, and (c) the data set itself, together determine the derived constraints. We illustrate our framework using measurements in simulated datasets before applying it to real data., 10 pages, 7 figures; v2: minor clarifications, fixed some typos, accepted in MNRAS

The masses of galaxy clusters can be measured using data obtained exclusively from wide photometric surveys in one of two ways: directly from the amplitude of the weak lensing signal or, indirectly, through the use of scaling relations calibrated using binned lensing measurements. In this paper, we build on a recently proposed idea and implement an alternative method based on the radial profile of the satellite distribution. This technique relies on splashback, a feature associated with the apocentre of recently accreted galaxies that offers a clear window into the phase-space structure of clusters without the use of velocity information. We carry out this dynamical measurement using the stacked satellite distribution around a sample of luminous red galaxies in the fourth data release of the Kilo-Degree Survey and validate our results using abundance-matching and lensing masses. To illustrate the power of this measurement, we combine dynamical and lensing mass estimates to robustly constrain scalar–tensor theories of gravity at cluster scales. Our results exclude departures from General Relativity of the order of unity. We conclude the paper by discussing the implications for future data sets. Because splashback mass measurements scale only with the survey volume, stage-IV photometric surveys are well-positioned to use splashback to provide high-redshift cluster masses.

We introduce the radiative transfer code Sweep for the cosmological simulation suite arepo. Sweep is a discrete ordinates method in which the radiative transfer equation is solved under the infinite speed of light, steady state assumption by a transport sweep across the entire computational grid. Since arepo is based on an adaptive, unstructured grid, the dependency graph induced by the sweep dependencies of the grid cells is non-trivial. In order to solve the topological sorting problem in a distributed manner, we employ a task-based-parallelism approach. The main advantage of the sweep method is that the computational cost scales only with the size of the grid and is independent of the number of sources or the distribution of sources in the computational domain, which is an advantage for radiative transfer in cosmological simulations, where there are large numbers of sparsely distributed sources. We successfully apply the code to a number of physical tests such as the expansion of H ii regions, the formation of shadows behind dense objects, the scattering of light, and its behaviour in the presence of periodic boundary conditions. In addition, we measure its computational performance with a focus on highly parallel, large-scale simulations.

Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project, which exploits the 'free' calibration data from the Atacama Large Millimetre/submillimetre Array (ALMA) to map the lines of sight towards and beyond the ALMA calibrators. ALMACAL has now covered 1,001 calibrators, with a total sky coverage around 0.3 deg2, distributed across the sky accessible from the Atacama desert, and has accumulated more than 1,000h of integration. The depth reached by combining multiple visits to each field makes ALMACAL capable of searching for faint, dusty, star-forming galaxies (DSFGs), with detections at multiple frequencies to constrain the emission mechanism. Based on the most up-to-date ALMACAL database, we report the detection of 186 DSFGs with flux densities down to S870um ~ 0.2mJy, comparable with existing ALMA large surveys but less susceptible to cosmic variance. We report the number counts at five wavelengths between 870um and 3mm, in ALMA bands 3, 4, 5, 6 and 7, providing a benchmark for models of galaxy formation and evolution. By integrating the observed number counts and the best-fitting functions, we also present the resolved fraction of the cosmic infrared background (CIB) and the CIB spectral shape. Combining existing surveys, ALMA has currently resolved about half of the CIB in the submm/mm regime., Comment: 22 pages, 16 figures, accepted for publication in MNRAS

We measure the local correlation between radio emission and Compton-$y$ signal across two galaxy clusters, Abell~399 and Abell~401, using maps from the Low-Frequency Array (LOFAR) and the Atacama Cosmology Telescope (ACT) + \Planck. These datasets allow us to make the first measurement of this kind at $\sim$arcminute resolution. We find that the radio brightness scales as $F_{\mathrm{radio}} \propto y^{1.5}$ for Abell~401 and $F_{\mathrm{radio}} \propto y^{2.8}$ for Abell~399. Furthermore, using \XMM data, we derive a sublinear correlation between radio and X-ray brightness for both the clusters ($F_{\mathrm{radio}} \propto F_{\rm X}^{0.7}$). Finally, we correlate the Compton-$y$ and X-ray data, finding that an isothermal model is consistent with the cluster profiles, $y \propto F_{\rm X}^{0.5}$. By adopting an isothermal--$\beta$ model, we are able, for the first time, to jointly use radio, X-ray, and Compton-$y$ data to estimate the scaling index for the magnetic field profile, $B(r) \propto n_{\mathrm{e}}(r)^{\eta}$ in the injection and re-acceleration scenarios. Applying this model, we find that the combined radio and Compton-$y$ signal exhibits a significantly tighter correlation with the X-ray across the clusters than when the datasets are independently correlated. We find $\eta \sim 0.6{-}0.8$. These results are consistent with the upper limit we derive for the scaling index of the magnetic field using rotation measure values for two radio galaxies in Abell~401. We also measure the radio, Compton-$y$, and X-ray correlations in the filament between the clusters but conclude that deeper data are required for a convincing determination of the correlations in the filament., Comment: 16 pages, 8 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

The recent rapid growth of the black hole (BH) catalog from gravitational waves (GWs), has allowed us to study the substructure of black hole mass function (BHMF) beyond the simplest Power-Law distribution. However, the BH masses inferred from binary BH merger events, may be systematically 'brightened' or 'dimmed' by the gravitational lensing effect. In this work, we investigate the impact of gravitational lensing on the BHMF inference considering the detection of the third-generation GW detector -- the Einstein Telescope (ET). We focus on high redshift, $z=10$, in order to obtain the upper limits of this effect. We use Monte Carlo (MC) method to simulate the data adopting 3 original BHMFs under Un-Lensed and Lensed scenarios, then recover the parameters of BHMFs from the mock data, and compare the difference of results, respectively. We found that all the parameters are well recovered within one standard deviation(std., 1$\sigma$), and all 3 BHMF models are reconstructed within 68\% credible interval, suggesting that lensing would not change the main structure drastically, even at very high redshifts and with high precision of ET. And the modest influence beyond $50M_{\odot}$, depends on the modeling of the high mass tail or substructure of BHMF. We conclude that the impact of lensing on BHMF inference with ET can be safely ignored in the foreseeable future. Careful handling of lensing effects is required only when focusing on an accurate estimation of the high mass end of BHMF at high redshifts., Comment: 13 pages, 11 figures, accepted by MNRAS

M. Dixon et al., We use stacked spectra of the host galaxies of photometrically identified Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilizing full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features (, MD would like to acknowledge support through an Australian Government Research Training Program Scholarship. This research was supported by the Australian Research Council The Centre of Excellence for Dark Matter Particle Physics (CDM; project number CE200100008) and the Australian Research Council Centre of Excellence for Gravitational Wave Discovery (OzGrav; project number CE170100004). This project/publication was made possible through the support of a grant from the John Templeton Foundation. The authors gratefully acknowledge this grant ID 61807, Two Standard Models Meet. The opinions expressed in this publication are those of the author(s) and do not necessarily reflect the views of the John Templeton Foundation. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft, and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, NSF’s NOIRLab, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory at NSF’s NOIRLab (NOIRLab Prop. ID 2012B-0001; PI: J. Frieman), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) do e-Universo (CNPq grant 465376/2014-2). This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. Based on data acquired at the Anglo-Australian Telescope, under program A/2013B/012. We acknowledge the traditional custodians of the land on which the AAT stands, the Gamilaraay people, and pay our respects to elders past and present.

We perform weak lensing tomographic peak studies using the first-year shear data from Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey. The effective area used in our analyses after field selection, mask and boundary exclusions is $\sim 58 \deg^2$. The source galaxies are divided into low- and high-redshift bins with $0.2\le z_p\le0.85$ and $0.85\le z_p\le1.5$, respectively. We utilize our halo-based theoretical peak model including the projection effect of large-scale structures to derive cosmological constraints from the observed tomographic high peak abundances with the signal-to-noise ratio in the range of $\nu_{\rm N}=[3.5,5.5]$. These high peaks are closely associated with the lensing effects of massive clusters of galaxies. Thus the inclusion of their member galaxies in the shear catalog can lead to significant source clustering and dilute their lensing signals. We account for this systematic effect in our theoretical modelling. Additionally, the impacts of baryonic effects, galaxy intrinsic alignments, as well as residual uncertainties in shear and photometric redshift calibrations are also analyzed. Within the flat $\Lambda$CDM model, the derived constraint is $S_8=0.758_{-0.076}^{+0.033}$ and $0.768_{-0.057}^{+0.030}$ with the source clustering information measured from the two cluster catalogs, CAMIRA and WZL, respectively. The asymmetric uncertainties are due to the different degeneracy direction of $(\Omega_{\rm m}, \sigma_8)$ from high peak abundances comparing to that from the cosmic shear two-point correlations which give rise approximately the power index $\alpha=0.5$. Fitting to our constraints, we obtain $\alpha\approx 0.38$ and $\Sigma_8=0.772_{-0.032}^{+0.028}$ (CAMIRA) and $0.781_{-0.033}^{+0.028}$ (WZL). In comparison with the results from non-tomographic peak analyses, the $1\sigma$ uncertainties on $\Sigma_8$ are reduced by a factor of $\sim1.3$., Comment: 19 pages, 11 figures, accepted for publication in MNRAS

The next generation of detectors will detect gravitational waves from binary neutron stars at cosmological distances, for which around a thousand electromagnetic follow-ups may be observed per year. So far, most work devoted to the expected cosmological impact of these standard sirens employed them only as distance indicators. Only recently their use as tracers of clustering, similar to what already proposed for supernovae, has been studied. Focusing on the expected specifications of the Einstein Telescope (ET), we forecast here the performance on cosmological parameters of future standard sirens as both distance and density indicators, with emphasis on the linear perturbation growth index and on spatial curvature. We improve upon previous studies in a number of ways: a more detailed analysis of available telescope time, the inclusion of more cosmological and nuisance parameters, the Alcock-Paczynski correction, the use of sirens also as both velocity and density tracers, and a more accurate estimation of the distance posterior. We find that the analysis of the clustering of sirens improves the constraints on $H_0$ by 30% and on $\Omega_{k0}$ by over an order of magnitude, with respect to their use merely as distance indicators. With 5 years of joint ET and Rubin Observatory follow-ups we could reach precision of 0.1 km/s/Mpc in $H_0$ and 0.02 in $\Omega_{k0}$ using only data in the range $0, Comment: v2: minor clarifications added; matches version to appear in MNRAS

Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ~MeV-TeV, these fluctuations are usually attributed either to 'extrinsic turbulence' (ET) - a cascade from larger scales - or 'self-confinement' (SC) - self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed live numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally-required scattering rates disagree with ET predictions by orders-of-magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvenic streaming are unstable. In live simulations, they revert to either arbitrarily-rapid CR escape with zero secondary production, or to bottleneck solutions with far-too-strong CR confinement and secondary production. Resolving these fundamental issues without discarding basic plasma processes requires invoking different drivers for scattering fluctuations. These must act on a broad range of scales with a power spectrum obeying several specific (but plausible) constraints., 36 pages, 7 figures. Updated to match published version, added section discussing 'meso-scale' phenomenology

The velocity profile of galaxies around voids is a key ingredient for redshift space distortion (RSD) measurements made using the void-galaxy correlation function. In this paper we use simulations to test whether the velocity profile of the tracers used to find the voids matches the velocity profile of the dark matter around these voids. A mismatch is expected and found in the inner part of voids, where tracers are very sparse. We discuss how this difference is caused by a selection effect where the void centre positions are correlated to the particular realization of the sparse tracers and their spatial distribution. In turn, this then affects the RSD void-galaxy correlation analysis. We show this by evaluating the Jacobian of the real to redshift space mapping using the tracer or matter velocity profile. Differences of the order of 20\% in the velocity profile translate into differences of the order of few percent in the Jacobian. This small discrepancy propagates to the monopole and quadrupole of the void-tracer correlation function, producing modifications of comparable magnitude to those from changes in $f\sigma_8$ at the level of the statistical uncertainties from current analyses., Comment: Author's Accepted Manuscript

The discovery of faint FR~I radio jets in the nearby elliptical galaxy NGC 5322 is reported here using 144 MHz LOFAR image. The jets have an angular extent of $\sim40$ arcmin or a projected physical extent of $\sim350$ kpc. The faint jets remain well collimated and disappear in the intergalactic medium, without any visible hotspot or radio lobes. The jets detected up to $\sim20$ kpc extent at higher frequencies are relatively bright within the optical extent of the galaxy but become faint abruptly outside, where detection is made only in the LOFAR image. The total radio luminosity of the galaxy at 144 MHz is estimated to be $3.7(\pm0.4)\times10^{22}$ W Hz$^{-1}$. The 144 MHz radio luminosity of the faint jets outside the optical extent is estimated to be $7.1(\pm2.0)\times10^{21}$ W Hz$^{-1}$. The extent of the jets for its radio luminosity is abnormally large when compared to the general population of radio galaxies. It makes NGC 5322 a member of a rare population of radio galaxies, previously not detected in other radio surveys. A combined effect of stellar core depletion and low-density environment around the jets, inferred from previous studies in other wave-bands, resulting into weak entrainment of surrounding material to the jets could be responsible for its large size despite a low radio luminosity., Comment: MNRAS submitted; 5 pages, 4 figures; comments welcome

In view of the recent measurement of $H_{0}$ from HST and SH0ES team, we explore the possibility of existence of a negative cosmological constant (AdS vacua in the dark energy sector) in the Universe. In this regard, we consider quintessence fields on top of a negative cosmological constant and compare such construction with $\Lambda$CDM model using a different combination of CMB, SnIa, BAO and $H_{0}$ data. Various model comparison estimators show that quintessence models with a negative $\Lambda$ is either preferred over $\Lambda$CDM or performs equally as $\Lambda$CDM model. This suggests that the presence of a negative $\Lambda$ (AdS ground state) in our Universe, which can naturally arise in string theory, is consistent with cosmological observations., Comment: 9 pages, Latex style, 6 figures, 6 tables, Revised version with new discussions and results, Accepted for publication in MNRAS

Primordial gravitational waves with extremely low frequency are expected to origin from inflation in the early Universe. The detection of such kind of gravitational waves is of great significance to verify the inflationary theory and determine the energy scale of inflation. B-mode of polarization of cosmic microwave background shows to be a promising probe of extremely low frequency primordial gravitational wave. In order to find a complementary observational feature induced by primordial gravitational wave with extremely low frequency, we propose an alternative way of detection by investigating the effect of primordial gravitational wave with extremely low frequency on a gravitational lens system with a non-aligned source-deflector-observer configuration. The results show that, with a series of chosen parameters, gravitational lens system with perturbation from extremely low frequency primordial gravitational wave could induce time delay which could deviate from the time delay deduced from theoretical model as much as about one hundred percent, meaning that gravitational lens system with a non-aligned configuration could serve as a potential long-base-line detector of extremely low frequency primordial gravitational wave., 6 pages, 4 figures, accepted for publication in MNRAS