Search

The Local Volume Complete Cluster Survey (LoVoCCS) is an on-going program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts $0.0310^{44}$ erg/s) with the Dark Energy Camera (DECam), capturing data in $u,g,r,i,z$ bands with a $5\sigma$ point source depth of approximately 25-26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19 deg within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32 deg within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales., Comment: 40 pages, 16 figures, 5 tables; revised and accepted for publication in ApJ

Galaxy clusters are the products of structure formation through myriad physical processes that affect their growth and evolution throughout cosmic history. As a result, the matter distribution within galaxy clusters, or their shape, is influenced by cosmology and astrophysical processes, in particular the accretion of new material due to gravity. We introduce an analysis method to investigate the 3D triaxial shapes of galaxy clusters from the Cluster HEritage project with XMM-Newton -- Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE). In this work, the first paper of a CHEX-MATE triaxial analysis series, we focus on utilizing X-ray data from XMM and Sunyaev-Zel'dovich (SZ) effect maps from Planck and ACT to obtain a three dimensional triaxial description of the intracluster medium (ICM) gas. We present the forward modeling formalism of our technique, which projects a triaxial ellipsoidal model for the gas density and pressure to compare directly with the observed two dimensional distributions in X-rays and the SZ effect. A Markov chain Monte Carlo is used to estimate the posterior distributions of the model parameters. Using mock X-ray and SZ observations of a smooth model, we demonstrate that the method can reliably recover the true parameter values. In addition, we apply the analysis to reconstruct the gas shape from the observed data of one CHEX-MATE galaxy cluster, Abell 1689, to illustrate the technique. The inferred parameters are in agreement with previous analyses for that cluster, and our results indicate that the geometrical properties, including the axial ratios of the ICM distribution, are constrained to within a few percent. With much better precision than previous studies, we thus further establish that Abell 1689 is significantly elongated along the line of sight, resulting in its exceptional gravitational lensing properties., Comment: accepted for publication in A&A

The HST treasury program BUFFALO provides extended wide-field imaging of the six Hubble Frontier Fields galaxy clusters. Here we present the combined strong and weak-lensing analysis of Abell 370, a massive cluster at z=0.375. From the reconstructed total projected mass distribution in the 6arcmin x 6arcmin BUFFALO field-of-view, we obtain the distribution of massive substructures outside the cluster core and report the presence of a total of seven candidates, each with mass $\sim 5 \times 10^{13}M_{\odot}$. Combining the total mass distribution derived from lensing with multi-wavelength data, we evaluate the physical significance of each candidate substructure, and conclude that 5 out of the 7 substructure candidates seem reliable, and that the mass distribution in Abell 370 is extended along the North-West and South-East directions. While this finding is in general agreement with previous studies, our detailed spatial reconstruction provides new insights into the complex mass distribution at large cluster-centric radius. We explore the impact of the extended mass reconstruction on the model of the cluster core and in particular, we attempt to physically explain the presence of an important external shear component, necessary to obtain a low root-mean-square separation between the model-predicted and observed positions of the multiple images in the cluster core. The substructures can only account for up to half the amplitude of the external shear, suggesting that more effort is needed to fully replace it by more physically motivated mass components. We provide public access to all the lensing data used as well as the different lens models., Comment: 29 pages, 17 figures, 3 tables

We present a statistical study of 180 dust continuum sources identified in 33 massive cluster fields by the ALMA Lensing Cluster Survey (ALCS) over a total of 133 arcmin$^{2}$ area, homogeneously observed at 1.2 mm. ALCS enables us to detect extremely faint mm sources by lensing magnification, including near-infrared (NIR) dark objects showing no counterparts in existing {\it Hubble Space Telescope} and {\it Spitzer} images. The dust continuum sources belong to a blind sample ($N=141$) with S/N $\gtrsim$ 5.0 (a purity of $>$ 0.99) or a secondary sample ($N=39$) with S/N= $4.0-5.0$ screened by priors. With the blind sample, we securely derive 1.2-mm number counts down to $\sim7$ $\mu$Jy, and find that the total integrated 1.2mm flux is 20.7$^{+8.5}_{-6.5}$ Jy deg$^{-2}$, resolving $\simeq$ 80 % of the cosmic infrared background light. The resolved fraction varies by a factor of $0.6-1.1$ due to the completeness correction depending on the spatial size of the mm emission. We also derive infrared (IR) luminosity functions (LFs) at $z=0.6-7.5$ with the $1/V_{\rm max}$ method, finding the redshift evolution of IR LFs characterized by positive luminosity and negative density evolution. The total (=UV+IR) cosmic star-formation rate density (SFRD) at $z>4$ is estimated to be $161^{+25}_{-21}$ % of the established measurements, which were almost exclusively based on optical$-$NIR surveys. Although our general understanding of the cosmic SFRD is unlikely to change beyond a factor of 2, these results add to the weight of evidence for an additional ($\approx 60$ %) SFRD component contributed by the faint-mm population, including NIR dark objects., Comment: ApJS in press, 34 pages, 8 Tables, and 18 Figures (+38 pages, 5 Tables, and 6 Figures in Appendix)

Backsplash galaxies are those that traverse and overshoot cluster cores as they fall into these structures. They are affected by environment, and should stand out in contrast to the infalling population. We target galaxies in the vicinity of clusters (R>R200) and select a sample in projected phase space (PPS), from the compilation of Sampaio et al. based on SDSS data. We present a statistical analysis, comparing two regions in PPS, with the same projected distance to the cluster but different velocity. The analysis relies on the presence of variations in the stellar population content of backsplash galaxies. We find a lower limit in the fractional contribution of ~5% with respect to the general sample of infalling galaxies at similar group-centric distance when using single line strength analysis, or ~15-30% when adopting bivariate distributions. The stellar populations show a subtle but significant difference towards older ages, and a higher fraction of quiescent galaxies. We also compare this set with a general field sample, where a substantially larger difference in galaxy properties is found, with the field sample being consistently younger, metal poorer and with a lower fraction of quiescent galaxies. Noting that our "cluster" sample is located outside of the virial radius, we expect this difference to be caused by pre-processing of the infalling galaxies in the overall higher density regions., Comment: 10 pages, 10 figures, accepted for publication in MNRAS

The radial acceleration relation (RAR) represents a tight empirical relation between the inferred total and baryonic centripetal accelerations, $g_{\rm{tot}}=GM_{\rm{tot}}(

The Subaru Hyper Suprime-Cam (HSC) Strategic Survey is the latest-generation multi-band optical imaging survey for galaxy evolution and structure formation. The "Ultra-Deep" component of the HSC survey provides $grizy$ broad-band images over $\sim3.4$ deg$^2$ to detection limits of $\sim26$-28 AB, along with narrow-band images, in the COSMOS and the SXDS fields. These images provide an unprecedented combination of depths and area coverage, for the studies galaxies up to $z\sim7$. However, the lack of coverage at $<4000$ Ang implies incomplete sampling of the rest-frame UV at $z\lesssim 3$, which is critically needed for understanding the buildup of stellar mass in the later cosmic time. We conducted a multi-year CFHT $u^\ast$-band imaging campaign in the two HSC Ultra-Deep fields with CFHT MegaCam. By including shallower archival data, we reach 5-$\sigma$ depths of $u^\ast=28.1$ and 28.4 (AB) at the centers of the COSMOS and SXDS fields, respectively, and $u^\ast=27.7$ and 27.8 in the central 1 deg$^2$ fields. The image quality is $\gtrsim0.90$ arcsec, fairly good for the $u^\ast$ band. Both the photometric and astrometric quality of our data are excellent. We show that the combination of our $u^\ast$-band and HSC data can lead to high-quality photometric redshifts at $z=0$-3, and robust measurements of rest-frame UV on galaxies at $0.4

Extended theories of gravity have been extensively investigated during the last thirty years, aiming at fixing infrared and ultraviolet shortcomings of General Relativity and of the associated $\Lambda$CDM cosmological model. Recently, non-local theories of gravity have drawn increasing attention due to their potential to ameliorate both the ultraviolet and infrared behavior of gravitational interaction. In particular, Integral Kernel theories of Gravity provide a viable mechanism to explain the late time cosmic acceleration so as to avoid the introduction of any form of unknown dark energy. On the other hand, these models represent a natural link towards quantum gravity. Here, we study a scalar-tensor equivalent model of General Relativity corrected with non-local terms, where corrections are selected by the existence of Noether symmetries. After performing the weak field limit and generalizing the results to extended mass distributions, we analyse the non-local model at galaxy cluster scales, by comparing the theoretical predictions with gravitational lensing observations from the CLASH program. We obtain agreement with data at the same level of statistical significance as General Relativity. We also provide constraints for the Navarro--Frenk--White parameters and lower bounds for the non-local length scales. The results are finally compared with those from the literature., Comment: 20 pages, 11 figures, 1 table. New version accepted for publication

We present a detailed weak-lensing and X-ray study of the Frontier Fields galaxy cluster Abell 370, one of the most massive known lenses on the sky, using wide-field BRz Subaru/Sprime-Cam and Chandra X-ray observations. By combining two-dimensional (2D) shear and azimuthally averaged magnification constraints derived from Subaru data, we perform a lensing mass reconstruction in a free-form manner, which allows us to determine both radial structure and 2D morphology of the cluster mass distribution. In a triaxial framework assuming a Navarro-Frenk-White density profile, we constrain the intrinsic structure and geometry of the cluster halo by forward modeling the reconstructed mass map. We obtain a halo mass $M_{200}=(1.54 \pm 0.29)\times 10^{15}h^{-1}M_\odot$, a halo concentration $c_{200}=5.27 \pm 1.28$, and a minor-major axis ratio $q_a=0.62 \pm 0.23$ with uninformative priors. Using a prior on the line-of-sight alignment of the halo major axis derived from binary merger simulations constrained by multi-probe observations, we find that the data favor a more prolate geometry with lower mass and lower concentration. From triaxial lens modeling with the line-of-sight prior, we find a spherically enclosed gas mass fraction of $f_\mathrm{gas}=(8.4 \pm 1.0)\%$ at $0.7h^{-1}$ Mpc $\sim 0.7r_{500}$. When compared to the hydrostatic mass estimate from Chandra observations, our triaxial weak-lensing analysis yields spherically enclosed mass ratios of $1-b \equiv M_\mathrm{HE}/M_\mathrm{WL} = 0.56 \pm 0.09$ and $0.51 \pm 0.09$ at $0.7h^{-1}$ Mpc with and without using the line-of-sight prior, respectively. Since the cluster is in a highly disturbed dynamical state, this represents the likely maximum level of hydrostatic bias in galaxy clusters., Comment: Version accepted by ApJ; 26 pages (main body) + 8 pages (appendices and references), 20 figures, 9 tables

We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxies, are still underdense compared to the mean density from the entire field. In addition, we classify galaxies into star-forming, quiescent, and green valley populations, and find that void galaxies tend to have slightly higher fractions of star-forming galaxies under the mass and redshift control, although the significance of this result is only moderate (2$\sigma$). However, when we focus on the star-forming population, the distribution of the specific star formation rate (sSFR) of void galaxies shows little difference from that of the control galaxies. Similarly, the median sSFR of star-forming void galaxies is also in good agreement with that of the star-forming control galaxies. Moreover, the effective green valley fraction of void galaxies, defined as the number of green valley galaxies over the number of nonquiescent galaxies, is comparable to that of the control ones, supporting the suggestion that void and control galaxies evolve under similar physical processes and quenching frequencies. Our results thus favor a scenario of the galaxy assembly bias., Comment: 17 pages, 10 figures, 3 tables, ApJ accepted

We use high-precision combined strong/weak lensing and kinematics measurements of the total mass profiles of the observed galaxy clusters MACS~J1206.2-0847 and Abell~S1063, to constrain the relativistic sector of the general DHOST dark energy theories, which exhibit a partial breaking of the so called Vainsthein screening mechanism, on the linear level of scalar fluctuations around a cosmological background. In particular, by using the \textsc{MG-MAMMPOSSt} framework developed in Pizzuti et al. 2021, for the kinematics analysis of member galaxies in clusters along with lensing mass profile reconstructions, we provide new constraints on the coupling $Y_2$ which governs the theory's relativistic contribution to the lensing potential. The new bound from the combination of kinematics and lensing measurements of MACS 1206, $Y_2=-0.12^{+0.66}_{-0.67}$ at $2\sigma$, provides about a 2-fold improvement on previous constraints. In the case of Abell~S1063 a $>2\sigma$ tension with the GR expectation arises. We discuss this in some detail, and we investigate the possible sources of systematics which can explain the tension. We further discuss why the combination of kinematics of member galaxies with lensing is capable of providing much tighter bounds compared to kinematics or lensing alone, and we explain how the number density profile of tracers, as well as the choice of the velocity anisotropy profile affects the final results., Comment: 11 pages, 11 figures. Submitted to MNRAS

We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass--total mass ($M_{500}$) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass--total mass relation is shallower than unity, $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of $M_{500}$ agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution., Comment: 34 pages, 16 figures, 8 tables; accepted for the publication in PASJ

Extended Theories of Gravity with additional scalar degrees of freedom have recently acquired increasing interest due to the presence of a screening mechanism that allows suppressing at small scales (e.g., the Solar System scale) every modification restoring General Relativity. In this work, we consider a second-order Extended Theory of Gravity belonging to the family of Degenerate High Order Scalar Tensor theories (DHOST) characterized by a partial breaking of the Vainshtein screening mechanism. We study this model in two different scenarios as a description of dark energy only and as a description of both dark matter and dark energy. Such scenarios have been tested here by analysing a sample of 16 high-mass galaxy clusters targeted by the Cluster Lensing and Supernova survey with Hubble (CLASH) program using two complementary probes, namely X-ray and strong-and-weak gravitational lensing observations. In mass modelling, we adopt a multi-component approach including hot gas and galactic stellar contributions. For the majority of the clusters in our sample, results show mild Bayesian evidence in favour of the DHOST model as a description of dark energy over General Relativity. This model also appears to alleviate the discrepancy present in General Relativity between X-ray hydrostatic and lensing mass estimates. For the second scenario where gravity acts as both dark energy and dark matter due to the partial breaking of the Vainshtein screening mechanism at cluster scales, the model is statistically disfavoured compared to General Relativity., Comment: 20 pages, 2 tables, 6 figures. Comments are welcome

We present the forward cosmological analysis of an $XMM$ selected sample of galaxy clusters out to a redshift of unity. Following our previous 2018 study based on the dn/dz quantity alone, we perform an upgraded cosmological analysis of the same XXL C1 cluster catalogue (178 objects), with a detailed account of the systematic errors. We follow the ASpiX methodology: the distribution of the observed X-ray properties of the cluster population is analysed in a 3D observable space (count rate, hardness ratio, redshift) and modelled as a function of cosmology. Compared to more traditional methods, ASpiX allows the inclusion of clusters down to a few tens of photons. We obtain an improvement by a factor of 2 compared to the previous analysis by letting the normalisation of the M-T relation and the evolution of the L-T relation free. Adding constraints from the XXL cluster 2-point correlation function and the BAO from various surveys decreases the uncertainties by 23 and 53 % respectively, and 62% when adding both. Switching to the scaling relations from the Subaru analysis, and letting free more parameters, our final constraints are $\sigma_8$ = $0.99^{+0.14}_{-0.23}$, $\Omega_m$ = 0.296 $\pm$ 0.034 ($S_8 = 0.98^{+0.11}_{-0.21}$) for the XXL sample alone. Finally, we combine XXL ASpiX, the XXL cluster 2-point correlation function and the BAO, with 11 free parameters, allowing for the cosmological dependence of the scaling relations in the fit. We find $\sigma_8$ = $0.793^{+0.063}_{-0.12}$, $\Omega_m$ = 0.364 $\pm$ 0.015 ($S_8 = 0.872^{+0.068}_{-0.12}$), but still compatible with Planck CMB at 2.2$\sigma$. The results obtained by the ASpiX method are promising; further improvement is expected from the final XXL cosmological analysis involving a cluster sample twice as large. Such a study paves the way for the analysis of the eROSITA and future Athena surveys., Comment: 20 pages, 10 figures, accepted for publication in A&A, A&A version has the unabridged abstract

Likelihood-free inference provides a rigorous approach to preform Bayesian analysis using forward simulations only. The main advantage of likelihood-free methods is its ability to account for complex physical processes and observational effects in forward simulations. Here we explore the potential of likelihood-free forward modeling for Bayesian cosmological inference using the redshift evolution of the cluster abundance combined with weak-lensing mass calibration. We use two complementary likelihood-free methods, namely Approximate Bayesian Computation (ABC) and Density-Estimation Likelihood-Free Inference (DELFI), to develop an analysis procedure for inference of the cosmological parameters $(\Omega_\mathrm{m},\sigma_8)$ and the mass scale of the survey sample. Adopting an eROSITA-like selection function and a 10-percent scatter in the observable-mass relation in a flat $\Lambda$CDM cosmology with $\Omega_\mathrm{m}=0.286$ and $\sigma_8=0.82$, we create a synthetic catalog of observable-selected NFW clusters in a survey area of 50 deg$^2$. The stacked tangential shear profile and the number counts in redshift bins are used as summary statistics for both methods. By performing a series of forward simulations, we obtain convergent solutions for the posterior distribution from both methods. We find that ABC recovers broader posteriors than DELFI, especially for the $\Omega_\mathrm{m}$ parameter. For a weak-lensing survey with a source density of $n_\mathrm{g}=20$ arcmin$^{-2}$, we obtain posterior constraints on $S_8=\sigma_8(\Omega_\mathrm{m}/0.3)^{0.3}$ of $0.836 \pm 0.032$ and $0.810 \pm 0.019$ from ABC and DELFI, respectively. The analysis framework developed in this study will be particularly powerful for cosmological inference with ongoing cluster cosmology programs, such as the XMM-XXL survey and the eROSITA all-sky survey, in combination with wide-field weak-lensing surveys., Comment: Version accepted by ApJ, with textual changes to improve clarity (e.g., more details on the simulation procedure added in Section 3; a general procedure for modeling the mass scaling relations added as Appendix A; toy model simulations moved to Appendix B); results and conclusions remain unchanged; 15 pages, 11 figures

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate., Comment: Published by JOSS. The package is available at https://github.com/skypyproject/skypy. Comments, issues and pull requests are welcome

We perform a systematic study of the recently discovered fundamental plane of galaxy clusters (CFP) using ~250 simulated clusters from The Three Hundred project, focusing on the stability of the plane against different temperature definitions and its dependence on the dynamical relaxation state of clusters. The CFP is characterised by $T\propto M_\mathrm{s}^\alpha\,r_\mathrm{s}^\beta$, defined with the gas temperature ($T$) and the characteristic halo scale radius and mass ($r_\mathrm{s}$ and $M_\mathrm{s}$) assuming a Navarro-Frenk-White halo description. We explore two definitions of weighted temperatures, namely mass-weighted and spectroscopic-like temperatures, in three radial ranges. The Three Hundred project clusters at $z=0$ lie on a thin plane whose parameters ($\alpha,\beta$) and dispersion ($0.015-0.030$ dex) depend on the gas temperature definition. The CFP for mass-weighted temperatures is closer to the virial equilibrium expectation ($\alpha=1,\beta=-1$) with a smaller dispersion. For gas temperatures measured within $500h^{-1}$ kpc, the resulting CFP deviates the most from the virial expectation and shifts towards the similarity solution for a secondary infall model ($\alpha=1.5,\beta=-2$). Independently of the temperature definition, we find that clusters at $z=1$ and relaxed clusters form a CFP similar to the virial expectation, unlike disturbed clusters exhibiting stronger evolution. Only systems formed over the last 4 Gyr present a CFP that is closer to the similarity solution. All these findings are compatible with the CFP obtained for a CLASH subsample excluding the hottest clusters with $T_\mathrm{X}>12$ keV., Comment: 19 pages, 10 figures, 9 tables. Accepted for publication in MNRAS

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) observed by Spectrum Roentgen Gamma/$eROSITA$ over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $\eta$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($\eta$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature $T_{\mathrm{X}}$, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged), Comment: Accepted for publication in A&A. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

We present an ensemble X-ray analysis of systematic perturbations in the central hot gas properties for a sample of 28 nearby strong cool-core systems selected from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS). We analyze their cool-core features observed with the Chandra X-ray Observatory. All individual systems in our sample exhibit at least a pair of positive and negative excess perturbations in the X-ray residual image after subtracting the global brightness profile. We extract and analyze X-ray spectra of the intracluster medium (ICM) in the detected perturbed regions. To investigate possible origins of the gas perturbations, we characterize thermodynamic properties of the ICM in the perturbed regions and characterize their correlations between positive and negative excess regions. The best-fit relations for temperature and entropy show a clear offset from the one-to-one relation, $T_\mathrm{neg}/T_\mathrm{pos}=1.20^{+0.04}_{-0.03}$ and $K_\mathrm{neg}/K_\mathrm{pos}=1.43\pm 0.07$, whereas the best-fit relation for pressure is found to be remarkably consistent with the one-to-one relation $P_\mathrm{neg}=P_\mathrm{pos}$, indicating that the ICM in the perturbed regions is in pressure equilibrium. These observed features in the HIFLUGCS sample are in agreement with the hypothesis that the gas perturbations in cool cores are generated by gas sloshing. We also analyze synthetic observations of perturbed cluster cores created from binary merger simulations, finding that the observed temperature ratio agrees with the simulations, $T_\mathrm{neg}/T_\mathrm{pos}\sim 1.3$. We conclude that gas sloshing induced by infalling substructures plays a major role in producing the characteristic gas perturbations in cool cores. The ubiquitous presence of gas perturbations in cool cores may suggest a significant contribution of gas sloshing to suppressing runaway cooling of the ICM., Comment: 24 pages, 11 figures, accepted for publication in ApJ

We report the non-thermal pressure fraction (Pnt/Ptot) obtained from a three-dimensional triaxial analysis of 16 galaxy clusters in the CLASH sample using gravitational lensing (GL) data primarily from Subaru and HST, X-ray spectroscopic imaging from Chandra, and Sunyaev-Zel'dovich effect (SZE) data from Planck and Bolocam. Our results span the approximate radial range 0.015-0.4R200m (35-1000 kpc). At cluster-centric radii smaller than 0.1R200m the ensemble average Pnt/Ptot is consistent with zero with an upper limit of nine per cent, indicating that heating from active galactic nuclei and other relevant processes does not produce significant deviations from hydrostatic equilibrium (HSE). The ensemble average Pnt/Ptot increases outside of this radius to approximately 20 per cent at 0.4R200m, as expected from simulations, due to newly accreted material thermalizing via a series of shocks. Also in agreement with simulations, we find significant cluster-to-cluster variation in Pnt/Ptot and little difference in the ensemble average Pnt/Ptot based on dynamical state. We conclude that on average, even for diverse samples, HSE-derived masses in the very central regions of galaxy clusters require only modest corrections due to non-thermal motions., Comment: published in MNRAS

We present constraints on the physical properties (including stellar mass, age, and star formation rate) of 207 $6\lesssim z \lesssim8$ galaxy candidates from the Reionization Lensing Cluster Survey (RELICS) and companion Spitzer-RELICS surveys. We measure photometry using T-PHOT and perform spectral energy distribution fitting using EA$z$Y and BAGPIPES. Of the 207 candidates for which we could successfully measure Spitzer fluxes, 23 were demoted to likely low redshift ($z<4$). Among the remaining high redshift candidates, we find intrinsic stellar masses between $1\times10^6\rm{M_{\odot}}$ and $4\times10^9\rm{M_\odot}$, and rest-frame UV absolute magnitudes between $-22.6$ and $-14.5$ mag. While our sample is mostly comprised of $L_{UV}/L^*_{UV}<1$ galaxies, there are a number of brighter objects in the sample, extending to $L_{UV}/L^*_{UV}\sim2$. The galaxies in our sample span approximately four orders of magnitude in stellar mass and star-formation rates, and exhibit ages ranging from maximally young to maximally old. We highlight 11 galaxies which have detections in Spitzer/IRAC imaging and redshift estimates $z\geq6.5$, several of which show evidence for some combination of evolved stellar populations, large contributions of nebular emission lines, and/or dust. Among these is PLCKG287+32-2013, one of the brightest $z\sim7$ candidates known (AB mag 24.9) with a Spitzer 3.6$\mu$m flux excess suggesting strong [OIII] + H-$\beta$ emission ($\sim$1000\AA\ rest-frame equivalent width). We discuss the possible uses and limits of our sample and present a public catalog of Hubble 0.4--1.6$\mu$m + Spitzer 3.6$\mu$m and 4.5$\mu$m photometry along with physical property estimates for all 207 objects in the sample. Because of their apparent brightnesses, high redshifts, and variety of stellar populations, these objects are excellent targets for follow-up with James Webb Space Telescope., Comment: 20 pages, 9 figures

Weak gravitational lensing of background galaxies provides a direct probe of the projected matter distribution in and around galaxy clusters. Here we present a self-contained pedagogical review of cluster--galaxy weak lensing, covering a range of topics relevant to its cosmological and astrophysical applications. We begin by reviewing the theoretical foundations of gravitational lensing from first principles, with special attention to the basics and advanced techniques of weak gravitational lensing. We summarize and discuss key findings from recent cluster--galaxy weak-lensing studies on both observational and theoretical grounds, with a focus on cluster mass profiles, the concentration--mass relation, the splashback radius, and implications from extensive mass calibration efforts for cluster cosmology., Comment: Invited review article published in The Astronomy and Astrophysics Review (see https://rdcu.be/caSUp); a version that corrects references and errors introduced by the publisher during the production process; 102 pages, 30 figures

We perform a self-calibration of the richness-to-mass ($N$-$M$) relation of CAMIRA galaxy clusters with richness $N\geq15$ at redshift $0.2\leq z<1.1$ by modeling redshift-space two-point correlation functions. These correlation functions are $\xi_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $\xi_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $\xi_{\mathrm{cg}}$ between these two samples. We focus on constraining the normalization $A_{\mathrm{N}}$ of the $N$-$M$ relation in a forward-modeling approach, carefully accounting for the redshift-space distortion, the Finger-of-God effect, and the uncertainty in photometric redshifts of CAMIRA clusters. The modeling also takes into account the projection effect on the halo bias of CAMIRA clusters. The parameter constraints are shown to be unbiased according to validation tests using a large set of mock catalogs constructed from N-body simulations. At the pivotal mass $M_{500}=10^{14}h^{-1}M_{\odot}$ and the pivotal redshift $z_{\mathrm{piv}} = 0.6$, the resulting normalization $A_{\mathrm{N}}$ is constrained as $13.8^{+5.8}_{-4.2}$, $13.2^{+3.4}_{-2.7}$, and $11.9^{+3.0}_{-1.9}$ by modeling $\xi_{\mathrm{cc}}$, $\xi_{\mathrm{cc}}+\xi_{\mathrm{cg}}$, and $\xi_{\mathrm{cc}} + \xi_{\mathrm{cg}} + \xi_{\mathrm{gg}}$, with average uncertainties at levels of $36\%$, $23\%$, and $21\%$, respectively. We find that the resulting $A_{\mathrm{N}}$ is statistically consistent with those independently obtained from weak-lensing magnification and from a joint analysis of shear and cluster abundance, with a preference for a lower value at a level of $\lesssim1.9\sigma$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged], Comment: Accepted for publication in the MNRAS. References to the figures are re-ordered. Add one more grant in the acknowledgement

Green valley galaxies represent the population that is likely to transition from the star-forming to the quiescent phases. To investigate the role of the environment in quenching star formation, we use the wide-field data from the Hyper Suprime-Cam Strategic Subaru Proposal survey to quantify the frequency of green valley galaxies in different environments and their redshift evolution. We find that the green valley fraction, in general, is less than 20% in any redshift and environment. The green valley fraction, when normalized to the total population, is higher in the field than that in groups or clusters and decreases with a decreasing redshift and increasing mass. The lower fraction of transitional galaxies in denser environments could be a consequence of the lack of star-forming galaxies, which could be the progenitors of green valley galaxies. To assess the effect of the environment on star formation quenching, we define the effective green valley fraction as the ratio of the number of green valley galaxies to that of nonquiescent galaxies only. The effective green valley fraction for field galaxies is lower than that for group or cluster galaxies, which reveals a strong positive mass dependence and mild redshift evolution. Moreover, the specific star formation rate (sSFR) is reduced by 0.1-0.3 dex in groups or clusters. Our results thus imply that an ongoing slow quenching process has been acting in the dense environment since z~1., Comment: 18 pages, 6 figures, and 5 tables. Accepted by ApJ

We study the dynamics of Abell 370 (A370), a highly massive Hubble Frontier Fields galaxy cluster, using self-consistent three-dimensional N-body/hydrodynamical simulations. Our simulations are constrained by X-ray, optical spectroscopic and gravitational lensing, and Sunyaev-Zel'dovich (SZ) effect observations. Analyzing archival Chandra observations of A370 and comparing the X-ray morphology to the latest gravitational lensing mass reconstruction, we find offsets of ~30 kpc and ~100 kpc between the two X-ray surface brightness peaks and their nearest mass surface density peaks, suggesting that it is a merging system, in agreement with previous studies. Based on our dedicated binary cluster merger simulations, we find that initial conditions of the two progenitors with virial masses of 1.7 x 10^(15) M_sun and 1.6 x 10^(15) M_sun, an infall velocity of 3500 km/s, and an impact parameter of 100 kpc can explain the positions and the offsets between the peaks of the X-ray emission and mass surface density, the amplitude of the integrated SZ signal, and the observed relative line-of-sight velocity. Moreover, our best model reproduces the observed velocity dispersion of cluster member galaxies, which supports the large total mass of A370 derived from weak lensing. Our simulations suggest that A370 is a major merger after the second core passage in the infalling phase, just before the third core passage. In this phase, the gas has not settled down in the gravitational potential well of the cluster, which explains why A370 does not follow closely the galaxy cluster scaling relations., Comment: 9 pages, 6 figures, and 1 table, accepted for publication in the Astrophysical Journal. Comments are welcome

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters., Comment: Accepted ApJS; MAST archive will be live concurrent with publication

The radial acceleration relation (RAR) in galaxies describes a tight empirical scaling law between the total acceleration $g_\mathrm{tot}(r)=GM_\mathrm{tot}(

We present a systematic study of gas density perturbations in cool cores of high-mass galaxy clusters. We select 12 relaxed clusters from the Cluster Lensing And Supernova survey with Hubble (CLASH) sample and analyze their cool core features observed with the Chandra X-ray Observatory. We focus on the X-ray residual image characteristics after subtracting their global profile of the X-ray surface brightness distribution. We find that all the galaxy clusters in our sample have, at least, both one positive and one negative excess regions in the X-ray residual image, indicating the presence of gas density perturbations. We identify and characterize the locally perturbed regions using our detection algorithm, and extract X-ray spectra of the intracluster medium (ICM). The ICM temperature in the positive excess region is lower than that in the negative excess region, whereas the ICM in both regions is in pressure equilibrium in a systematic manner. These results indicate that gas sloshing in cool cores takes place in more than 80% of relaxed clusters (95% CL). We confirm this physical picture by analyzing synthetic X-ray observations of a cool-core cluster from a hydrodynamic simulation, finding that our detection algorithm can accurately extract both the positive and negative excess regions and can reproduce the temperature difference between them. Our findings support the picture that the gas density perturbations are induced by gas sloshing, and a large fraction of cool-core clusters have undergone gas sloshing, indicating that gas sloshing may be capable of suppressing runaway cooling of the ICM., Comment: 30 pages, 8 figures, accepted for publication in ApJ

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties., Comment: 23 pages, accepted for publication in ApJ

Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($\beta_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($\beta_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($\beta_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $\rho_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $\rho_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$)., Comment: 19 pages; in press on MNRAS. One of two companion papers on HSC-XXL surveys (Keiichi Umetsu et al., arXiv:1909.10524)

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions: a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The main SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that the integrated Compton $y$ parameter and the temperature for the major merger are boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations., Comment: 36 pages, 25 figures, 5 tables; accepted for the publication in MNRAS

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures., Comment: Version matching the one published in ApJ. We recommend to use statistically corrected mass estimates (M200MT, M500MT) of Table 2 for a given individual cluster. One of two companion papers presenting initial HSC-XXL results (Mauro Sereno et al., arXiv:1912.02827)

We present a statistical weak-lensing magnification analysis on an optically selected sample of 3029 \texttt{CAMIRA} galaxy clusters with richness $N>15$ at redshift $0.2\leq z <1.1$ in the Subaru Hyper Suprime-Cam (HSC) survey. We use two distinct populations of color-selected, flux-limited background galaxies, namely the low-$z$ and high-$z$ samples at mean redshifts of $\approx1.1$ and $\approx1.4$, respectively, from which to measure the weak-lensing magnification signal by accounting for cluster contamination as well as masking effects. Our magnification bias measurements are found to be uncontaminated according to validation tests against the "null-test" samples for which the net magnification bias is expected to vanish. The magnification bias for the full \texttt{CAMIRA} sample is detected at a significance level of $9.51\sigma$, which is dominated by the high-$z$ background. We forward-model the observed magnification data to constrain the normalization of the richness-to-mass ($N$--$M$) relation for the \texttt{CAMIRA} sample with informative priors on other parameters. The resulting scaling relation is $N\propto {M_{500}}^{0.92\pm0.13} (1 + z)^{-0.48\pm0.69}$, with a characteristic richness of $N=\left(17.72\pm2.60\right)$ and intrinsic log-normal scatter of $0.15\pm0.07$ at $M_{500} = 10^{14}h^{-1}M_{\odot}$. With the derived $N$--$M$ relation, we provide magnification-calibrated mass estimates of individual \texttt{CAMIRA} clusters, with the typical uncertainty of $\approx39\%$ and $\approx32\%$ at richness$\approx20$ and $\approx40$, respectively. We further compare our magnification-inferred $N$--$M$ relation with those from the shear-based results in the literature, finding good agreement., Comment: 24 pages, 16 figures. Accepted for publication in the MNRAS

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals., Comment: 29 pages, 6 figures, submitted to ApJ. For reduced images, catalogs, lens models, and more, see relics.stsci.edu

According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_\Delta\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses., Comment: Invited review article, to be published in "From Dark Haloes to Visible Galaxies", special issue of Galaxies

We present line-of-sight gas sloshing first found in a cool core in a galaxy cluster. The galaxy cluster Abell 907 is identified as a relaxed cluster owing to its global X-ray surface brightness taken by the Chandra X-ray Observatory. The X-ray residual image after removing the global emission of the intracluster medium (ICM), however, shows an arc-like positive excess and a negative excess surrounding the central positive excess in the cluster core, which in turn indicates a disturbance of the ICM. We analyze the X-ray spectra extracted from both regions and find that (1) the ICM temperature and the metal abundance in the positive excess are lower and higher than those in the negative excess, respectively, and (2) the ICM is nearly in pressure equilibrium. We also find a slight redshift difference between the positive and the negative excesses, which corresponds to the velocity shear of $1680^{+1300}_{-920}$ km s$^{-1}$ ($1\sigma$). The X-ray residual image and the ICM properties are consistent with those expected by line-of-sight gas sloshing. Assuming that the gas is moving toward inverse-parallel to each other along the line-of-sight, the shear velocity is expected to be $\sim 800$ km s$^{-1}$. The velocity field of this level is able to provide non-thermal pressure support by $\sim 34\%$ relative to the thermal one. The total kinetic energy inferred from the shear velocity corresponds to $\sim 30\%$ of the bolometric luminosity of the sloshing ICM. Abell 907 is therefore complementary to galaxy clusters in which gas sloshing takes place in the plane of the sky, and is important for understanding gas dynamics driven by sloshing and its influence on the heating to prevent runaway cooling., Comment: 9 pages, 3 figures, accepted for publication in ApJ

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST)., Comment: 20 pages, 8 figures. Accepted to ApJ

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive., Comment: 15 Pages, 7 Figures, 4 Tables Accepted for publication in ApJ

We investigate the splashback features of dark-matter halos based on cosmic density and velocity fields. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding velocity statistic, alignment momentum correlation function, to take into account the asphericity of halos. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity. On halo scales, $x\sim R_\mathrm{200m} \sim 1Mpc/h$, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. We also find that the splashback radius determined from the density correlation becomes $\sim 3.5\%$ smaller than that from the momentum correlation, with their correlation coefficient being 0.605. Moreover, the orientation-dependent splashback feature due to halo asphericity is measured when the density profile is determined by dark-matter particles, which can be used as a test of collisional CDM since the halo shape is predicted to be rounder in such a model., Comment: 8 pages, 4 figures, published in PRD. This paper is a considerably extended version of the section on Splashback radius of arXiv:1706.08860

To investigate the dynamics of the galaxy cluster A2142, we compile an extended catalog of 2239 spectroscopic redshifts of sources, including newly measured 237 redshifts, within 30 arcmin from the cluster center. With the $\sigma$-plateau algorithm from the caustic method, we identify 868 members and a number of substructures in the galaxy distribution both in the outskirts, out to $\sim$3.5 Mpc from the cluster center, and in the central region. In the outskirts, one substructure overlaps a falling clump of gas previously identified in the X-ray band. These substructures suggests the presence of multiple minor mergers, which are responsible for the complex dynamics of A2142, and the absence of recent or ongoing major mergers. We show that the distribution of the galaxies in the cluster core and in several substructures are consistent with the mass distribution inferred from the weak lensing signal. Moreover, we use spatially-resolved X-ray spectroscopy to measure the redshift of different regions of the intracluster medium within $\sim$3 arcmin from the cluster center. We find a ring of gas near the two X-ray cold fronts identified in previous analyses and measure a velocity of this ring of $810 \pm 330 {\rm km s^{-1}}$ larger than the cluster mean velocity. Our analysis suggests the presence of another ring surrounding the core, whose velocity is $660 \pm 300 {\rm km s^{-1}}$ larger than the cluster velocity. These X-ray features are not associated to any optical substructures, and support the core-sloshing scenario suggested in previous work.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known., Comment: 15 pages, 8 figures 4 tables. ApJ Accepted

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements., Comment: 19 pages, 11 figures, 2 tables, comments are welcome

Observations and numerical simulations have shown that the relation between the mass scaled with the critical density of the universe and the X-ray temperature of galaxy clusters is approximately represented by $M_\Delta \propto T_X^{3/2}$ (e.g. $\Delta=500$). This relation is often interpreted as evidence that clusters are in virial equilibrium. However, the recently discovered fundamental plane (FP) of clusters indicates that the temperature of clusters primarily depends on a combination of the characteristic mass $M_s$ and radius $r_s$ of the Navarro-Frenk-White profile rather than $M_\Delta$. Moreover, the angle of the FP revealed that clusters are not in virial equilibrium because of continuous mass accretion from the surrounding matter. By considering both the FP and the mass dependence of the cluster concentration parameter, we show that this paradox can be solved and the relation $M_\Delta \propto T_X^{3/2}$ actually reflects the central structure of clusters. We also find that the intrinsic scatter in the halo concentration-mass relation can largely account for the spread of clusters on the FP. We also show that X-ray data alone form the FP and the angle and the position are consistent with those of the FP constructed from gravitational lensing data. We demonstrate that a possible shift between the two FPs can be used to calibrate cluster masses obtained via X-ray observations., Comment: Published on ApJ. Matched to published version

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter., Comment: Minor changes to match the version published in ApJ. One of three new companion papers of the CLUMP-3D project (I-Non Chiu et al., arXiv:1804.00676; Mauro Sereno et al., arXiv:1804.00667)

We perform a three-dimensional triaxial analysis of 16 X-ray regular and 4 high-magnification galaxy clusters selected from the CLASH survey by combining two-dimensional weak-lensing and central strong-lensing constraints. In a Bayesian framework, we constrain the intrinsic structure and geometry of each individual cluster assuming a triaxial Navarro-Frenk-White halo with arbitrary orientations, characterized by the mass $M_{200\mathrm{c}}$, halo concentration $C_{200\mathrm{c}}$, and triaxial axis ratios ($q_{\mathrm{a}} \le q_{\mathrm{b}}$), and investigate scaling relations between these halo structural parameters. From triaxial modeling of the X-ray-selected subsample, we find that the halo concentration decreases with increasing cluster mass, with a mean concentration of $C_{200\mathrm{c}} = 4.82\pm0.30$ at the pivot mass $M_{200\mathrm{c}}=10^{15}M_{\odot}h^{-1}$. This is consistent with the result from spherical modeling, $C_{200\mathrm{c}}=4.51\pm 0.14$. Independently of the priors, the minor-to-major axis ratio $q_{\mathrm{a}}$ of our full sample exhibits a clear deviation from the spherical configuration ($q_{\mathrm{a}}=0.52 \pm 0.04$ at $10^{15}M_{\odot}h^{-1}$ with uniform priors), with a weak dependence on the cluster mass. Combining all 20 clusters, we obtain a joint ensemble constraint on the minor-to-major axis ratio of $q_{\mathrm{a}}=0.652^{+0.162}_{-0.078}$ and a lower bound on the intermediate-to-major axis ratio of $q_{\mathrm{b}}>0.63$ at the $2\sigma$ level from an analysis with uniform priors. Assuming priors on the axis ratios derived from numerical simulations, we constrain the degree of triaxiality for the full sample to be $\mathcal{T}=0.79 \pm 0.03$ at $10^{15}M_{\odot}h^{-1}$, indicating a preference for a prolate geometry of cluster halos. We find no statistical evidence for an orientation bias ($f_{\mathrm{geo}}=0.93 \pm 0.07$) (abridged), Comment: Accepted by the ApJ. This is one of the three companion papers, including Umetsu et al. 2018 (arXiv:1804.00664) and Sereno et al. 2018 (arXiv:1804.00667), in the CLUMP-3D project. Go go go, HoHo!

The $\Lambda$CDM model of structure formation makes strong predictions on concentration and shape of DM (dark matter) halos, which are determined by mass accretion processes. Comparison between predicted shapes and observations provides a geometric test of the $\Lambda$CDM model. Accurate and precise measurements needs a full three-dimensional analysis of the cluster mass distribution. We accomplish this with a multi-probe 3D analysis of the X-ray regular CLASH (Cluster Lensing And Supernova survey with Hubble) clusters combining strong and weak lensing, X-ray photometry and spectroscopy, and the Sunyaev-Zel'dovich effect. The cluster shapes and concentrations are consistent with $\Lambda$CDM predictions. The CLASH clusters are randomly oriented, as expected given the sample selection criteria. Shapes agree with numerical results for DM-only halos, which hints at baryonic physics being not so effective in making halos rounder., Comment: v2: 8 pages, in press on ApJL. Extended discussion on regularity. One of three new companion papers of the CLUMP-3D project (Keiichi Umetsu et al., arxiv:1804.00664; I-Non Chiu et al., arXiv:1804.00676)

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s., Comment: 23 pages, 16 figures; V2, accepted for publication in ApJ

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes., Comment: 14 pages, 6 figures, 4 tables, accepted for publication in ApJ

We present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong., Comment: 25 pages, 17 figures Appeared in Publication of Astronomical Society of Japan, Hyper Suprime-Cam Special Issue

In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane., Comment: Replaced with a revised version to match the ApJ accepted version