Search

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

Next-generation wide-field spectroscopic surveys will observe the infall regions around large numbers of galaxy clusters with high sampling rates for the first time. Here we assess the feasibility of extracting the large-scale cosmic web around clusters using forthcoming observations, given realistic observational constraints. We use a sample of 324 hydrodynamic zoom-in simulations of massive galaxy clusters from TheThreeHundred project to create a mock-observational catalogue spanning $5R_{200}$ around 160 analogue clusters. These analogues are matched in mass to the 16 clusters targetted by the forthcoming WEAVE Wide-Field Cluster Survey (WWFCS). We consider the effects of the fibre allocation algorithm on our sampling completeness and find that we successfully allocate targets to 81.7 $\% \pm$ 1.3 of the members in the cluster outskirts. We next test the robustness of the filament extraction algorithm by using a metric, $D_{\text{skel}}$, which quantifies the distance to the filament spine. We find that the median positional offset between reference and recovered filament networks is $D_{\text{skel}} = 0.13 \pm 0.02$ Mpc, much smaller than the typical filament radius of $\sim$ 1 Mpc. Cluster connectivity of the recovered network is not substantially affected. Our findings give confidence that the WWFCS will be able to reliably trace cosmic web filaments in the vicinity around massive clusters, forming the basis of environmental studies into the effects of pre-processing on galaxy evolution., 18 pages, 14 figures. Accepted by MNRAS for publication

We investigate how the environment affects the assembly history of massive galaxies. For that purpose, we make use of SHARDS and HST spectro-photometric data, whose depth, spectral resolution, and wavelength coverage allow to perform a detailed analysis of the stellar emission as well as obtaining unprecedentedly accurate photometric redshifts. This expedites a sufficiently accurate estimate of the local environment and a robust derivation of the star formation histories of a complete sample of 332 massive galaxies ($\mathrm{>10^{10}M_{\odot}}$) at redshift $1\leq z \leq 1.5$ in the GOODS-N field. We find that massive galaxies in this redshift range avoid the lowest density environments. Moreover, we observed that the oldest galaxies in our sample with with mass-weighted formation redshift $\mathrm{\overline{z}_{M-w} \geq 2.5}$, avoid the highest density regions, preferring intermediate environments. Younger galaxies, including those with active star formation, tend to live in denser environments ($\Sigma = \mathrm{5.0_{1.1}^{24.8}\times 10^{10}M_{\odot}Mpc^{-2}}$). This behavior could be expected if those massive galaxies starting their formation first would merge with neighbors and sweep their environment earlier. On the other hand, galaxies formed more recently ($\overline{z}_{M-w} < 2.5$) are accreted into large scale structures at later times and we are observing them before sweeping their environment or, alternatively, they are less likely to affect their environment. However, given that both number and mass surface densities of neighbor galaxies is relatively low for the oldest galaxies, our results reveal a very weak correlation between environment and the first formation stages of the earliest massive galaxies., Comment: Accepted for publication in MNRAS

We present angular power spectra and cosmological parameter constraints derived from the Planck PR4 (NPIPE) maps of the Cosmic Microwave Background. NPIPE, released by the Planck Collaboration in 2020, is a new processing pipeline for producing calibrated frequency maps from Planck data. We have created new versions of the CamSpec likelihood using these maps and applied them to constrain LCDM and single-parameter extensions. We find excellent consistency between NPIPE and the Planck 2018 maps at the parameter level, showing that the Planck cosmology is robust to substantial changes in the mapmaking. The lower noise of NPIPE leads to ~10% tighter constraints, and we see both smaller error bars and a shift toward the LCDM values for beyond-LCDM parameters including Omega_K and A_Lens., Comment: 17 pages, 20 figures, updated to match accepted version

In general relativity (GR), the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. Recent work has claimed a statistical detection of an `external field effect' (EFE) using galaxy rotation curve data. We show that large uncertainties in rotation curve analyses and inaccuracies in published simulation-based external field estimates compromise the significance of the claimed EFE detection. We further show analytically that a qualitatively similar statistical signal is, in fact, expected in a $\Lambda$-cold dark matter ($\Lambda$CDM) universe without any violation of the strong equivalence principle. Rather, such a signal arises simply because of the inherent correlations between galaxy clustering strength and intrinsic galaxy properties. We explicitly demonstrate the effect in a baryonified mock catalog of a $\Lambda$CDM universe. Although the detection of an EFE-like signal is not, by itself, evidence for physics beyond GR, our work shows that the $\textit{sign}$ of the EFE-like correlation between the external field strength and the shape of the radial acceleration relation can be used to probe new physics: e.g., in MOND, the predicted sign is opposite to that in our $\Lambda$CDM mocks., Comment: 11 pages, 6 figures; v2: minor changes and added discussion, accepted in MNRAS

We present results from nine simulations that compare the standard $\Lambda$ Cold Dark Matter cosmology ($\Lambda$CDM) with counterfactual universes, for approximately $100\,{\rm Gyr}$ using the Enzo simulation code. We vary the value of $\Lambda$ and the fluctuation amplitude to explore the effect on the evolution of the halo mass function (HMF), the intergalactic medium (IGM) and the star formation history (SFH). The distinct peak in star formation rate density (SFRD) and its subsequent decline are both affected by the interplay between gravitational attraction and the accelerating effects of $\Lambda$. The IGM cools down more rapidly in models with a larger $\Lambda$ and also with a lower $\sigma_8$, reflecting the reduced SFRD associated with these changes -- although changing $\sigma_8$ is not degenerate with changing $\Lambda$, either regarding the thermal history of the IGM or the SFH. However, these induced changes to the IGM or ionizing background have little impact on the calculated SFRD. We provide fits for the evolution of the SFRD in these different universes, which we integrate over time to derive an asymptotic star formation efficiency. Together with Weinberg's uniform prior on $\Lambda$, the estimated probability of observers experiencing a value of $\Lambda$ no greater than the observed value is 13%, substantially larger than some alternative estimates. Within the Enzo model framework, then, observer selection within a multiverse is able to account statistically for the small value of the cosmological constant, although $\Lambda$ in our universe does appear to be at the low end of the predicted range., Comment: 18 pages, 10 figures, MNRAS: MN-22-1662-MJ.R1

We study the structural properties like the gravitational mass, radius and tidal deformability of dark matter (DM) admixed strange quark stars (SQSs). For the purpose we consider the vector MIT Bag model to describe the strange quark matter (SQM) and investigate the possible presence of accreted DM in the SQSs consequently forming DM admixed SQSs. We introduce feeble interaction between SQM and the accreted fermionic DM via a vector dark boson mediator. Considering the present literature, in the context of possible presence of DM in SQSs, this work is the first to consider interaction between DM and SQM in the DM admixed SQSs. The mass of the DM fermion ($m_{\chi}$) and the vector mediator ($m_{\xi}$) and the coupling ($y_{\xi}$) between them are determined in accordance with the constraint from Bullet cluster and the present day relic abundance, respectively. We find that the presence of DM reduces both the mass and radius of the star compared to the no-DM case. The massive the DM fermion, the lower the values of maximum mass and radius of the DM admixed SQSs. For the chosen values of $m_{\chi}$ and corresponding values of $m_{\xi}$ and $y_{\xi}$, the computed structural properties of the DM admixed SQSs satisfy all the various present day astrophysical constraints.We obtain massive DM admixed SQSs configurations consistent with the GW190814 observational data. Hence the secondary compact object associated with this event may be a DM admixed SQS., Comment: Accepted for Publication in Monthly Notices of the Royal Astronomical Society

The wavelet scattering transform (WST) has recently gained attention in the context of large-scale structure studies, being a possible generator of summary statistics encapsulating non-Gaussianities beyond the reach of the conventional power spectrum. This work examines the three-dimensional solid harmonic WST in the context of a three-dimensional line-intensity mapping measurement to be undertaken by current and proposed phases of the CO Mapping Array Project (COMAP). The WST coefficients demonstrate interpretable behaviour in the context of noiseless CO line-intensity simulations. The contribution of the cosmological $z\sim3$ signal to these coefficients is also detectable in principle even in the Pathfinder phase of COMAP. Using the peak-patch method to generate large numbers of simulations and incorporating observational noise, we numerically estimate covariance matrices and show that careful choices of WST hyperparameters and rescaled or reduced coefficient sets are both necessary to keep covariances well-conditioned. Fisher forecasts show that even a reduced `shapeless' set of $\ell$-averaged WST coefficients show constraining power that can exceed that of the power spectrum alone even with similar detection significance. The full WST could improve parameter constraints even over the combination of the power spectrum and the voxel intensity distribution, showing that it uniquely encapsulates shape information about the line-intensity field. However, practical applications urgently require further understanding of the WST in key contexts like covariances and cross-correlations., Comment: 15 pages; 6 figures, 5 tables; v2 reflects the version accepted for publication in MNRAS

The structural and dynamic properties of the dark matter halos, though an important ingredient in understanding large-scale structure formation, require more conservative particle resolution than those required by halo mass alone in a simulation. This reduces the parameter space of the simulations, more severely for high-redshift and large-volume mocks which are required by the next-generation large sky surveys. Here, we incorporate redshift and cosmology dependence into an algorithm that assigns accurate halo properties such as concentration, spin, velocity, and spatial distribution to the sub-resolution haloes in a simulation. By focusing on getting the right correlations with halo mass and local tidal anisotropy $\alpha$ measured at $4 \times$ halo radius, our method will also recover the correlations of these small scale structural properties with the large-scale environment, i.e., the halo assembly bias at all scales greater than $5 \times$ halo radius. We find that the distribution of halo properties is universal with redshift and cosmology. By applying the algorithm to a large volume simulation $(600 h^{-1}{\rm Mpc})$, we can access the $30-500$ particle haloes, thus gaining an order of magnitude in halo mass and two to three orders of magnitude in number density at $z=2-4$. This technique reduces the cost of mocks required for the estimation of covariance matrices, weak lensing studies, or any large-scale clustering analysis with less massive haloes., Comment: 9 figures, 6 tables, accepted by MNRAS, access to code at https://github.com/rsujatha/mockhaloprop

We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs ($M_{\rm vir}\sim 10^{13}M_{\odot}$ and $M_{*}\sim10^{11}M_{\odot}$), simulated in CDM and a SIDM model with constant cross-section of $\sigma_T/m_\chi = 1\ \mathrm{cm}^2 \mathrm{g}^{-1}$. We measure the three-dimensional and projected shapes of the dark matter haloes and their baryonic content using the inertia tensor and compare our measurements to the results of three $HST$ samples of gravitational lenses and $Chandra$ and $XMM-Newton$ X-ray observations. We find that the inclusion of baryons greatly reduces the differences between CDM and a SIDM, together with the ability to draw constraints based on shapes. Lensing measurements reject the predictions of CDM dark-matter-only simulations and prefer one of the hydro scenarios. When we consider the total sample of lenses, observational data prefer the CDM hydro scenario. The shapes of the X-ray emitting gas are compatible with observational results in both hydro runs, with CDM predicting higher elongations only in the very centre. Contrary to previous claims at the scale of elliptical galaxies, we conclude that both CDM and our SIDM model can still explain observed shapes once we include baryons in the simulations. Our results demonstrate that this is essential to derive realistic constraints and that new simulations are needed to confirm and extend our findings., Comment: 17 pages, 13 figures, revised version accepted for publication in MNRAS

Type Ia supernovae (SNe Ia) are powerful tools for measuring the expansion history of the universe, but the impact of dust around SNe Ia remains unknown and is a critical systematic uncertainty. One way to improve our empirical description of dust is to analyse highly reddened SNe Ia ($E(B-V)>0.4$, roughly equivalent to the fitted SALT2 light-curve parameter $c>0.3$). With the recently released Pantheon+ sample, there are 57 SNe Ia that were removed because of their high colour alone (with colours up to $c=1.61$), which can provide enormous leverage on understanding line-of-sight $R_V$. Previous studies have claimed that $R_V$ decreases with redder colour, though it is unclear if this is due to limited statistics, selection effects, or an alternative explanation. To test this claim, we fit two separate colour-luminosity relationships, one for the main cosmological sample ($c0.3$) SNe Ia. We find the change in the colour-luminosity coefficient to be consistent with zero. Additionally, we compare the data to simulations with different colour models, and find that the data prefers a model with a flat dependence of $R_V$ on colour over a declining dependence. Finally, our results strongly support that line-of-sight $R_V$ to SNe Ia is not a single value, but forms a distribution., Accepted to MNRAS, updating to match accepted version. 11 pages, 9 figures, 5 tables

Most cosmic shear analyses to date have relied on summary statistics (e.g. $\xi_+$ and $\xi_-$). These types of analyses are necessarily sub-optimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative to single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps -- and to the extent that the convergence field can be described by the lognormal model -- our map posteriors allow us to reconstruct \it all \rm summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the $\sigma_8$--$\Omega_{\rm m}$ plane by $\approx 30\%$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $S_8=\sigma_8\Omega_{\rm m}^{1/2}$ directions, meaning map-based inference fails to significantly improve constraints on $S_8$., Comment: 12 pages, 11 figures, To be submitted to MNRAS, Comments welcome

The relation between the integrated thermal Sunyaev-Zeldovich (tSZ) $y$-decrement versus halo mass ($Y$--$M$) can potentially constrain galaxy formation models, if theoretical and observational systematics can be properly assessed. We investigate the $Y$--$M$ relation in the SIMBA and IllustrisTNG-100 cosmological hydrodynamic simulations, quantifying the effects of feedback, line-of-sight projection, and beam convolution. We find that SIMBA's AGN jet feedback generates strong deviations from self-similar expectations for the $Y$--$M$ relation, especially at $M_{500}, 14 pages, 7 figures. Accepted by MNRAS, comments welcome

We present the first comparison between properties of clusters of galaxies detected by the eROSITA Final Equatorial-Depth Survey (eFEDS) and the XMM Cluster Survey (XCS). We have compared, in an ensemble fashion, properties from the eFEDS X-ray cluster catalogue with those from the Ultimate XMM eXtragaLactic (XXL) survey project (XXL-100-GC). We find the distributions of redshift and X-ray temperature ($T_{\rm X}$) to be broadly similar between the two surveys, with a larger proportion of clusters above 4 keV in the XXL-100-GC sample. We find 62 eFEDS cluster candidates with XMM data (eFEDS-XMM sample); 10 do not have good enough XMM data to confirm or deny, 11 are classed as sample contaminants, and 4 have their X-ray flux contaminated by another source. The majority of eFEDS-XMM sources have a longer exposure in XMM than eFEDS, and the majority of eFEDS positions are within 100 kpc of XCS positions. Our eFEDS-XCS sample of 37 clusters is used to calculate minimum sample contamination fractions of ${\sim}$18% and ${\sim}$9% in the eFEDS X-ray and optically confirmed samples respectively, in general agreement with eFEDS findings. We compare 29 X-ray luminosities ($L_{\rm X}$) measured by eFEDS and XCS, which are in excellent agreement. Eight clusters have a $T_{\rm X}$ measured by XMM and eROSITA, and we find that XMM temperatures are 25$\pm$9% larger than their eROSITA counterparts. Finally, we construct $L_{\rm X}$ - $T_{\rm X}$ scaling relations based on eFEDS and XCS measurements, which are in tension; the tension is decreased when we measure a third scaling relation with calibrated XCS temperatures., 20 pages (13 + 7 appendices), 16 figures, submitted to MNRAS

We present observational constraints on the galaxy-halo connection, focusing particularly on galaxy assembly bias, from a novel combination of counts-in-cylinders statistics, $P(N_{\rm{CIC}})$, with the standard measurements of the projected two-point correlation function, $w_{\rm{p}}(r_{\rm{p}})$, and number density, $n_{\rm{gal}}$, of galaxies. We measure $n_{\rm{gal}}$, $w_{\rm{p}}(r_{\rm{p}})$ and $P(N_{\rm{CIC}})$ for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the $M_r, 16+5 pages, 9+3 figures, 5+2 tables, Fig. 6 shows the main result. To be submitted to MNRAS, comments welcome

The mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) is probing supermassive black holes (SMBHs) in galaxies across the Hubble sequence via molecular gas dynamics. We present the first WISDOM study of a luminous infrared galaxy with an active galactic nuclei (AGN): Fairall 49. We use new ALMA observations of the CO(2-1) line with a spatial resolution of about 80 pc together with ancillary HST imaging. We reach the following results: (1) The CO kinematics are well described by a regularly rotating gas disk with a radial inflow motion, suggesting weak feedback on the cold gas from both AGN and starburst activity; (2) The dynamically inferred SMBH mass is 1.6 +/- 0.4 (rnd) +/- 0.8 (sys) x 10^8 Msun, assuming that we have accurately subtracted the AGN and starburst light contributions, which have a luminosity of about 10^9 Lsun; (3) The SMBH mass agrees with the SMBH-stellar mass relation but is about 50 times higher than previous estimates from X-ray variability; (4) The dynamically inferred molecular gas mass is 30 times smaller than that inferred from adopting the Galactic CO-to-H_2 conversion factor (X_CO) for thermalised gas, suggesting low values of X_CO; (5) the molecular gas inflow rate increases steadily with radius and may be as high as 5 Msun/yr. This work highlights the potential of using high-resolution CO data to estimate, in addition to SMBH masses, the X_CO factor and gas inflow rates in nearby galaxies., 19 pages, 13 figures, accepted for publication in MNRAS

We have applied the method of star counts with Wolf diagrams to determine the interstellar extinction in five Galactic cirri in Sloan Digital Sky Survey (SDSS) Stripe 82. For this purpose, we have used the photometry of stars in the GALEX NUV filter and the photometry of red dwarfs in five SDSS bands and four SkyMapper Southern Sky Survey DR2 bands. We have identified the cirri as sky regions with an enhanced infrared emission from the Schlegel+1998 map. The extinction in them has been calculated relative to the nearby comparison regions with a reduced emission. The results for different filters agree well, giving the range of distances and the extinction law for each cirrus. The distances in the range 140--415 pc found are consistent with the 3D reddening maps. In the range between the $B$ and $V$ filters the extinctions found are consistent with the estimates from Schlegel+1998 for the Cardelli+1989 extinction law with $R_\mathrm{V}=3.1$. However, the extinctions found for all of the filters are best described not by the Cardelli+1989 extinction law with some $R_\mathrm{V}=3.1$, but by the inverse proportionality of the extinction and wavelength with its own coefficient for each cirrus. In one of the cirri our results suggest a very slight decrease in extinction with wavelength, i.e., a large contribution of gray extinction. In the remaining cirri a manifestation of gray extinction is not ruled out either. This is consistent with the previous measurements of the extinction law far from the Galactic midplane., 4 figures, published in Astronomy Letter

Many sources contribute to the diffuse gamma-ray background (DGRB), including star forming galaxies, active galactic nuclei, and cosmic ray interactions in the Milky Way. Exotic sources, such as dark matter annihilation, may also make some contribution. The photon counts-in-pixels distribution is a powerful tool for analyzing the DGRB and determining the relative contributions of different sources. However, including photon energy information in a likelihood analysis of the counts-in-pixels distribution quickly becomes computationally intractable as the number of source types and energy bins increase. Here, we apply the likelihood-free method of Approximate Bayesian Computation (ABC) to the problem. We consider a mock analysis that includes contributions from dark matter annihilation in galactic subhalos as well as astrophysical backgrounds. We show that our results using ABC are consistent with the exact likelihood when energy information is discarded, and that significantly tighter parameter constraints can be obtained with ABC when energy information is included. ABC presents a powerful tool for analyzing the DGRB and understanding its varied origins., 10 pages main text, 2 pages appendices, 7 figures; updated to match version accepted by MNRAS

We study the stacked filaments connecting group-mass halo pairs, using dark-matter-only N-body simulations. We calculate the dark matter over-density profile of these stacked filaments at different redshifts as a function of the distance perpendicular to the filament axis. A four-parameter universal functional form, including three comoving scale radii and one amplitude parameter (core density), provides a good fit out to a radius of 20 cMpc/h for stacked filaments over a range of redshifts, lengths and masses. The scale radii are approximately independent of redshift but increase as power-laws with the comoving filament length. Lastly, we compare the scaling of the filament mass measured directly from the simulations to the predicted scaling from the halo-halo-matter three-point correlation function as a function of redshift and of the mass of the halo pairs. We find that both measured scalings are similar to, but somewhat shallower than the predictions, by 10% and 30%, respectively. These results provide a template to interpret present and upcoming observational results based on stacking, for example, weak lensing, thermal and kinetic Sunyaev-Zel'dovich, or X-ray observations., Comment: 15 pages. 11 figures. Accepted by MNRAS. Comments welcome

The shape of the faint-end of the high-z galaxy luminosity function (LF) informs early star formation and reionization physics during the Cosmic Dawn and Epoch of Reionization. Until recently, based on the strong gravitational lensing cluster deep surveys, the Hubble Frontier Fields (HFF) has found a potential turnover in the ultraviolet (UV) LF at z$\sim$6. In this paper, we analyze the contribution of extremely faint galaxies with the magnitude larger than the turnover magnitude in LF to cosmic reionization. We apply the measurement from HFF to our suppressed star formation efficiency model, including three free parameters: halo mass threshold $M_t$, curvature parameter $\beta$ and a UV conversion factor $l_{\rm UV}$. According to our fit of 68\% confidence level, the high-redshift star formation in haloes smaller than $ M_t=1.82^{+2.86}_{-1.08}\times10^{10} \rm M_{\odot}$ is found to be dampened. The turnover magnitude $\rm \gtrsim -13.99-2.45$, correspondingly the halo mass $\lesssim(4.57+20.03)\times10^{9} \rm M_{\odot}$. We find that the absorption trough in the global 21-cm signal is sensitive to our SFE model parameters. Together with ($\beta$, $l_{\rm UV}$) = ($2.17^{+2.42}_{-1.72}$, $9.33^{+0.43}_{-0.42} \rm ~erg~yr ~s^{-1}M_{\odot}^{-1})$, the trough locates at $\sim$ $134^{+10}_{-17}$ $\rm MHz$ with an amplitude of $\sim$ $-237^{-6}_{+7}$ $\rm mK$, compared to (106\rm MHz, -212\rm mK) in the absence of turnover. Besides, we find that the star formation of faint galaxies has also an impact on the 21-cm power spectra. The best fitting peak power decreases by $\sim4\%$ and shifts towards smaller scales from $0.88 h \rm Mpc^{-1}$ to $0.91 h \rm Mpc^{-1}$. According to our calculation, such impact is distinguishable with the forthcoming Square Kilometre Array., Comment: 11 pages, 5 figures; Accepted by MNRAS

Primordial black holes (PBHs) formed in the early Universe constitute an attractive candidate for dark matter. Within the gaseous environment of the interstellar medium, PBHs with accretion disks naturally launch outflows such as winds and jets. PBHs with significant spin can sustain powerful relativistic jets and generate associated cocoons. Jets and winds can efficiently deposit their kinetic energies and heat the surrounding gas through shocks. Focusing on the Leo T dwarf galaxy, we demonstrate that these considerations can provide novel tests of PBHs over a significant $\sim 10^{-2} M_{\odot} - 10^6 M_{\odot}$ mass range, including the parameter space associated with gravitational wave observations by the LIGO and VIRGO Collaborations. Observing the morphology of emission could allow to distinguish between jet and wind contributions, and hence indirectly detect spinning PBHs., 10 pages, 4 figures

The merger rate density evolution of binary compact objects and the properties of their host galaxies carry crucial information to understand the sources of gravitational waves. Here, we present galaxyRate, a new code that estimates the merger rate density of binary compact objects and the properties of their host galaxies, based on observational scaling relations. We generate our synthetic galaxies according to the galaxy stellar mass function. We estimate the metallicity according to both the mass-metallicity relation (MZR) and the fundamental metallicity relation (FMR). Also, we take into account galaxy-galaxy mergers and the evolution of the galaxy properties from the formation to the merger of the binary compact object. We find that the merger rate density changes dramatically depending on the choice of the star-forming galaxy main sequence, especially in the case of binary black holes (BBHs) and black hole neutron star systems (BHNSs). The slope of the merger rate density of BBHs and BHNSs is steeper if we assume the MZR with respect to the FMR, because the latter predicts a shallower decrease of metallicity with redshift. In contrast, binary neutron stars (BNSs) are only mildly affected by both the galaxy main sequence and metallicity relation. Overall, BBHs and BHNSs tend to form in low-mass metal-poor galaxies and merge in high-mass metal-rich galaxies, while BNSs form and merge in massive galaxies. We predict that passive galaxies host at least ~5-10%, ~15-25%, and ~15-35% of all BNS, BHNS and BBH mergers in the local Universe., 21 pages, 22 figures (including appendices), 3 tables, published in MNRAS

We revisited the problem of mixing in a gravitational N-body system from the point of view of the ordering of coarse-grained cells in the one-particle energy space, here denoted {\it energy ranking preservation} (ERP). This effect has been noted for some time in simulations, although individual particle energies and their phase-space variables mix considerably. The present investigation aimed to map ERP in terms of parameters involving the collective range in which it is effective, as well as in terms of global and historical characterisations of gravitational systems evolving towards equilibrium. We examined a subset of the IllustrisTNG cosmological magnetohydrodynamical simulations (TNG50-4 and TNG100-3), considering both their full and dark-only versions. For each simulation, we selected the $20$ most massive haloes at redshift $z=0$, tracing their ERP fractions back at selected redshift markers ({\tt z} $= \{1.0, 5.0, 10.0 \}$), and for a coarse-graining set ranging from $5$ to $30$ energy bins. At the redshift marker {\tt z} $= 1$, we found high ERP fractions (above $\sim 80 \%$) in both simulations, regardless of the coarse-graining level. The {\it decline} in ERP fractions with redshift was roughly a function of mass and fractional mass increase in the analysed TNG50-4 haloes, but not in the TNG100-3 ones, indicating a possible relative susceptibility of the ERP effect to mass accretion for haloes less massive than $\sim 10^{14} ~ M_{\odot}$. We confirmed earlier indications in the literature concerning a possible "mesoscopic" constraint operative in a time span of at least several Gyr., 23 pages, 14 figures (+ 6 panels in the Appendix), final accepted version to MNRAS

The cross-correlation between 21-cm intensity mapping experiments and photometric surveys of galaxies (or any other cosmological tracer with a broad radial kernel) is severely degraded by the loss of long-wavelength radial modes due to Galactic foreground contamination. Higher-order correlators are able to restore some of these modes due to the non-linear coupling between them and the local small-scale clustering induced by gravitational collapse. We explore the possibility of recovering information from the bispectrum between a photometric galaxy sample and an intensity mapping experiment, in the context of the clustering-redshifts technique. We demonstrate that the bispectrum is able to calibrate the redshift distribution of the photometric sample to the required accuracy of future experiments such as the Rubin Observatory, using future single-dish and interferometric 21-cm observations, in situations where the two-point function is not able to do so due to foreground contamination. We also show how this calibration is affected by the photometric redshift width $\sigma_{z,0}$ and maximum scale $k_{\mathrm{max}}$. We find that it is important to reach scales $k \gtrsim 0.3\,h\,\mathrm{Mpc}^{-1}$, with the constraints saturating at around $k\sim 1\,h\,\mathrm{Mpc}^{-1}$ for next-generation experiments., Comment: 20 pages, 8 figures

We outline a concept for OWL-Moon, a 50-100m aperture telescope located on the surface of the Moon, to address three major areas in astronomy, namely the detection of biosignatures on habitable exoplanets, the geophysics of exoplanets, and cosmology. Such a large lunar telescope, when coupled with large Earth-based telescopes, would allow Intensity Interferometric measurements, leading to pico-arcsecond angular resolution. This would have applications in many areas of astronomy and is timely in light of the renewed interest of space agencies in returning to the Moon., Comment: Accepted in Experimental Astronomy. arXiv admin note: substantial text overlap with arXiv:1908.02080

Numerical simulations with varying realism indicate an emergent principle -- multiphase condensation and large cavity power occur when the ratio of the cooling time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) falls below a threshold value close to 10. Observations indeed show cool-core signatures when this ratio falls below 20-30, but the prevalence of cores with $t_{\rm cool}/t_{\rm ff}$ ratio below 10 is rare as compared to simulations. In X-ray observations, we obtain projected spectra from which we have to infer radial gas density and temperature profiles. Using idealized models of X-ray cavities and multiphase gas in the core and 3-D hydro jet-ICM simulations, we quantify the biases introduced by deprojection based on the assumption of spherical symmetry in determining $t_{\rm cool}/t_{\rm ff}$. We show that while the used methods are able to recover the $t_{\rm cool}/t_{\rm ff}$ ratio for relaxed clusters, they have an uncertainty of a factor of $2-3$ in systems containing large cavities ($\gtrsim 20$ kpc). We also show that the mass estimates from these methods, in the absence of X-ray spectra close to the virial radius, suffer from a degeneracy between the virial mass ($M_{200}$) and the concentration parameter ($c$) in the form of $M_{200} c^2 \approx$ constant. Additionally, lack of soft-X-ray ($\lesssim 0.5$ keV) coverage and poor spatial resolution make us overestimate min($t_{\rm cool}/t_{\rm ff}$) by a factor of few in clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$. This bias can largely explain the lack of cool-core clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$., Accepted for publication in MNRAS. Comments are welcome

Strong gravitational lensing effect is a powerful tool to probe cosmological models and gravity theories. Recently, the time-delay cosmography from strong lensing and the stellar kinematics of the deflector, which encode the Hubble constant and the post-Newtonian parameter via two distance ratios reflecting the lensing mass and dynamical mass respectively, have been proposed to investigate these two parameters simultaneously. Among strong lensing systems with different sources, strongly lensed fast radio bursts (FRBs) have been proposed as precision probes of the universe since the time delay $\sim$ 10 days between images could be measured extremely precisely because of their short duration of a few milliseconds. In this work, we investigate the ability of strongly lensed FRBs on simultaneously estimating these two parameters via simulations. Take the expected FRB detection rate of upcoming facilities and lensing probability into consideration, it is likely to accumulate 10 lensed FRBs in several years and we find that $H_0$ could be determined to a $\sim1.5\%$ precision and $\gamma_{\rm PPN}$ could be constrained to a $\sim8.7\%$ precision simultaneously from them. These simultaneous estimations will be helpful for properly reflecting the possible correlation between these two fundamental parameters., Comment: 6 pages, 2 figures

We perform a flexion-based weak gravitational analysis of the first two Hubble Frontier Field clusters: Abell 2744 and MACS 0416. A parametric method for using radially projected flexion signals as a probe of cluster member mass is described in detail. The normalization and slope of a L − θE (as a proxy for L − σ) scaling relation in each cluster is determined using measured flexion signals. A parallel field analysis is undertaken concurrently to provide a baseline measure of method effectiveness. We find an agreement in the Faber–Jackson slope ℓ associated with galaxy age and morphology for both clusters, as well as a theoretical distinction in the cluster normalization mass.

Realistic lightcone mocks are important in the clustering analyses of large galaxy surveys. For simulations where only the snapshots are available, it is common to create approximate lightcones by joining together the snapshots in spherical shells. We assess the two-point clustering measurements of central galaxies in approximate lightcones built from the Millennium-XXL simulation, which are constructed using different numbers of snapshots. The monopole and quadrupole of the real-space correlation function is strongly boosted on small scales below 1 Mpc/h, due to some galaxies being duplicated at the boundaries between snapshots in the lightcone. When more snapshots are used, the total number of duplicated galaxies is approximately constant, but they are pushed to smaller separations. The effect of this in redshift space is small, as long as the snapshots are cut into shells in real space. Randomly removing duplicated galaxies is able to reduce the excess clustering signal. Including satellite galaxies will reduce the impact of the duplicates, since many small-scale pairs come from satellites in the same halo. Galaxies that are missing from the lightcone at the boundaries can be added to the lightcone by having a small overlap between each shell. This effect will impact analyses that use very small-scale clustering measurements, and when using mocks to test the impact of fibre collisions., 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

Monthly notices of the Royal Astronomical Society 516(2), 1923 - 1940 (2022). doi:10.1093/mnras/stac2207, Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self-scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function, and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall, we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide – for the first time – upper limits on the cross-section for frequent self-interactions., Published by Oxford Univ. Press, Oxford

We introduce AutoEnRichness, a hybrid approach that combines empirical and analytical strategies to determine the richness of galaxy clusters (in the redshift range of $0.1 \leq z \leq 0.35$) using photometry data from the Sloan Digital Sky Survey Data Release 16, where cluster richness can be used as a proxy for cluster mass. In order to reliably estimate cluster richness, it is vital that the background subtraction is as accurate as possible when distinguishing cluster and field galaxies to mitigate severe contamination. AutoEnRichness is comprised of a multi-stage machine learning algorithm that performs background subtraction of interloping field galaxies along the cluster line-of-sight and a conventional luminosity distribution fitting approach that estimates cluster richness based only on the number of galaxies within a magnitude range and search area. In this proof-of-concept study, we obtain a balanced accuracy of $83.20$ per cent when distinguishing between cluster and field galaxies as well as a median absolute percentage error of $33.50$ per cent between our estimated cluster richnesses and known cluster richnesses within $r_{200}$. In the future, we aim for AutoEnRichness to be applied on upcoming large-scale optical surveys, such as the Legacy Survey of Space and Time and $\textit{Euclid}$, to estimate the richness of a large sample of galaxy groups and clusters from across the halo mass function. This would advance our overall understanding of galaxy evolution within overdense environments as well as enable cosmological parameters to be further constrained., Comment: Accepted to MNRAS on 1st August 2022; pp. 1-21. Supplementary material is included; pp. 22-29

The cross-correlation of optical galaxies with the neutral hydrogen (HI) radiation intensity can enhance the signal-to-noise ratio (SNR) of the HI intensity measurement. In this paper, we investigate the cross-correlation of the galaxy samples obtained by the spectroscopic survey of the China Space Station Telescope (CSST) with the HI Intensity mapping (IM) survey of the Five-hundred-meter Aperture Spherical Telescope (FAST). Using the IllusitrisTNG simulation result at redshift $0.2 \sim 0.3$, we generate mock data of the CSST survey and a FAST L-band drift scan survey. The CSST spectroscopic survey can yield a sample of galaxies with a high comoving number density of $10^{-2} (\unit{Mpc}/h)^{-3}$ at $z \sim 0.3$. We cross-correlate the foreground-removed radio intensity with the CSST galaxies, including both the whole sample, and red and blue galaxy sub-samples separately. We find that in all cases the HI and optical galaxies are well correlated. The total HI abundance can be measured with a high precision from this correlation. A relative error of $\sim 0.6\%$ for $\Omega_{\rm HI}$ could be achieved at $z\sim 0.3$ for an overlapping survey area of $10000 \unit{deg}^2$., Comment: 16 pages, 10 figures

Accurate computations of spectral distortions of the cosmic microwave background (CMB) are required for constraining energy release scenarios at redshifts $z\gtrsim 10^3$. The existing literature focuses on distortions that are small perturbations to the background blackbody spectrum. At high redshifts ($z\gtrsim 10^6$), this assumption can be violated, and the CMB spectrum can be significantly distorted at least during part of its cosmic evolution. In this paper, we carry out accurate thermalization computations, evolving the distorted CMB spectrum in a general, fully non-linear way, consistently accounting for the time-dependence of the injection process, modifications to the Hubble expansion rate and relativistic Compton scattering. Specifically, we study single energy injection and decaying particle scenarios, discussing constraints on these cases. We solve the thermalization problem using two independent numerical approaches that are now available in {\tt CosmoTherm} as dedicated setups for computing CMB spectral distortions in the large distortion regime. New non-linear effects at low frequencies are furthermore highlighted, showing that these warrant a more rigorous study. This work eliminates one of the long-standing simplifications in CMB spectral distortion computations, which also opens the way to more rigorous treatments of distortions induced by high-energy particle cascade, soft photon injection and in the vicinity of primordial black holes., Comments welcome