Showing 1–50 of 128 results for author: Hudson, M J

Cosmological parameters from the joint analysis of Density Split and Second Order Statistics: an Emulator based on the Halo Occupation Distribution

Authors: Pierre A. Burger, Enrique Paillas, Michael J. Hudson

Abstract: In this work, we develop a simulation-based model to predict the density split (DSS) and second-order shear and clustering statistics. A simulation-based model has the potential to model highly non-linear scales where current DSS models fail. To build this model, we use the \texttt{AbacusSummit} N-body simulation suite from which we measure all necessary statistics and train an emulator based on \… ▽ More In this work, we develop a simulation-based model to predict the density split (DSS) and second-order shear and clustering statistics. A simulation-based model has the potential to model highly non-linear scales where current DSS models fail. To build this model, we use the \texttt{AbacusSummit} N-body simulation suite from which we measure all necessary statistics and train an emulator based on \texttt{CosmoPower}. In that context, we discuss possible improvements for future emulators to make the measurement less noisy and biased, resulting in more accurate and precise model predictions. Regarding the emulator's accuracy, we find that the most important aspect is the average of the summary statistics over multiple-shot noise realizations of the foreground galaxies. However, these results probably depend on the chosen number density of the foreground galaxies. Regarding the parameter forecast based on preliminary LOWZxUNIONS data, we find that DSS has more constraining power to derive cosmological parameters and that the joint analysis with second-order statistics is particularly useful for extracting parameters of the galaxy-halo connection. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 16 pages, 16 Figures, submitted to physical review D

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Black-Hole-to-Halo Mass Relation From UNIONS Weak Lensing

Authors: Qinxun Li, Martin Kilbinger, Wentao Luo, Kai Wang, Huiyuan Wang, Anna Wittje, Hendrik Hildebrandt, Ludovic van Waerbeke, Michael J. Hudson, Samuel Farrens, Tobias I. Liaudat, Huiling Liu, Ziwen Zhang, Qingqing Wang, Elisa Russier, Axel Guinot, Lucie Baumont, Fabian Hervas Peters, Thomas de Boer, Jiaqi Wang

Abstract: This letter presents, for the first time, direct constraints on the black-hole-to-halo-mass relation using weak gravitational lensing measurements. We construct type I and type II Active Galactic Nuclei (AGNs) samples from the Sloan Digital Sky Survey (SDSS), with a mean redshift of 0.4 0.1 for type I (type II) AGNs. This sample is cross-correlated with weak lensing shear from the Ultraviolet Near… ▽ More This letter presents, for the first time, direct constraints on the black-hole-to-halo-mass relation using weak gravitational lensing measurements. We construct type I and type II Active Galactic Nuclei (AGNs) samples from the Sloan Digital Sky Survey (SDSS), with a mean redshift of 0.4 0.1 for type I (type II) AGNs. This sample is cross-correlated with weak lensing shear from the Ultraviolet Near Infrared Northern Survey (UNIONS). We compute the excess surface mass density of the halos associated with $36,181$ AGNs from $94,308,561$ lensed galaxies and fit the halo mass in bins of black-hole mass. We find that more massive AGNs reside in more massive halos. We see no evidence of dependence on AGN type or redshift in the black-hole-to-halo-mass relationship when systematic errors in the measured black-hole masses are included. Our results are consistent with previous measurements for non-AGN galaxies. At a fixed black-hole mass, our weak-lensing halo masses are consistent with galaxy rotation curves, but significantly lower than galaxy clustering measurements. Finally, our results are broadly consistent with state-of-the-art hydro-dynamical cosmological simulations, providing a new constraint for black-hole masses in simulations. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 14 pages, 5 figures. Submitted to ApJL

Cosmological parameters estimated from velocity -- density comparisons: Calibrating 2M++

Authors: Amber M. Hollinger, Michael J. Hudson

Abstract: Cosmological parameters can be measured by comparing peculiar velocities with those predicted from a galaxy density field. Previous work has tested the accuracy of this approach with N-body simulations, but generally on idealised mock galaxy surveys. However, systematic biases may arise solely due to survey selection effects such as flux-limited samples, edge-effects, and complications due to the… ▽ More Cosmological parameters can be measured by comparing peculiar velocities with those predicted from a galaxy density field. Previous work has tested the accuracy of this approach with N-body simulations, but generally on idealised mock galaxy surveys. However, systematic biases may arise solely due to survey selection effects such as flux-limited samples, edge-effects, and complications due to the obscuration of the Galactic plane. In this work, we explore the impact of each of these effects individually well as collectively, using the semi-analytic models from numerical simulations to generate mock catalogues that mimic the 2M++ density field. We find the reconstruction and analysis methods used for our 2M++ mocks produce a value of fsigma_8 that is biased high by a factor 1.04 \pm 0.01 compared to the true value. Moreover, a cosmic volume matching that of 2M++ has a cosmic variance uncertainty in fsigma_8 of ~5%. The systematic bias is a function of distance: it is unbiased close to the origin but is biased slightly high for distances in the range 100-180 Mpc/h. Correcting for this small bias, we find the linear fsigma_8 = 0.362 \pm 0.023. The predicted peculiar velocities from 2M++ have an error of 170 km/s that slowly increases with distance, exceeding 200 km/s only at distances of 180-200 Mpc/h. Finally, the residual bulk flow speeds found in previous work are shown to be not in conflict with those expected in the Lambda cold dark matter model. △ Less

Submitted 25 March, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 18 pages, 10 figures

The discovery of the faintest known Milky Way satellite using UNIONS

Authors: Simon E. T. Smith, William Cerny, Christian R. Hayes, Federico Sestito, Jaclyn Jensen, Alan W. McConnachie, Marla Geha, Julio Navarro, Ting S. Li, Jean-Charles Cuillandre, Raphaël Errani, Ken Chambers, Stephen Gwyn, Francois Hammer, Michael J. Hudson, Eugene Magnier, Nicolas Martin

Abstract: We present the discovery of Ursa Major III/UNIONS 1, the least luminous known satellite of the Milky Way, which is estimated to have an absolute V-band magnitude of $+2.2^{+0.4}_{-0.3}$ mag, equivalent to a total stellar mass of 16$^{+6}_{-5}$ M$_{\odot}$. Ursa Major III/UNIONS 1 was uncovered in the deep, wide-field Ultraviolet Near Infrared Optical Northern Survey (UNIONS) and is consistent with… ▽ More We present the discovery of Ursa Major III/UNIONS 1, the least luminous known satellite of the Milky Way, which is estimated to have an absolute V-band magnitude of $+2.2^{+0.4}_{-0.3}$ mag, equivalent to a total stellar mass of 16$^{+6}_{-5}$ M$_{\odot}$. Ursa Major III/UNIONS 1 was uncovered in the deep, wide-field Ultraviolet Near Infrared Optical Northern Survey (UNIONS) and is consistent with an old ($τ> 11$ Gyr), metal-poor ([Fe/H] $\sim -2.2$) stellar population at a heliocentric distance of $\sim$ 10 kpc. Despite being compact ($r_{\text{h}} = 3\pm1$ pc) and composed of so few stars, we confirm the reality of Ursa Major III/UNIONS 1 with Keck II/DEIMOS follow-up spectroscopy and identify 11 radial velocity members, 8 of which have full astrometric data from $Gaia$ and are co-moving based on their proper motions. Based on these 11 radial velocity members, we derive an intrinsic velocity dispersion of $3.7^{+1.4}_{-1.0}$ km s$^{-1}$ but some caveats preclude this value from being interpreted as a direct indicator of the underlying gravitational potential at this time. Primarily, the exclusion of the largest velocity outlier from the member list drops the velocity dispersion to $1.9^{+1.4}_{-1.1}$ km s$^{-1}$, and the subsequent removal of an additional outlier star produces an unresolved velocity dispersion. While the presence of binary stars may be inflating the measurement, the possibility of a significant velocity dispersion makes Ursa Major III/UNIONS 1 a high priority candidate for multi-epoch spectroscopic follow-ups to deduce to true nature of this incredibly faint satellite. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 21 pages, 9 figures, 3 tables; Accepted for publication in ApJ

How many stars form in galaxy mergers?

Authors: Andrew M. M. Reeves, Michael J. Hudson

Abstract: We forward model the difference in stellar age between post-coalescence mergers and a control sample with the same stellar mass, environmental density, and redshift. In particular, we use a pure sample of 445 post-coalescence mergers from the recent visually-confirmed post-coalescence merger sample identified by Bickley et al. and find that post-coalescence mergers are on average younger than cont… ▽ More We forward model the difference in stellar age between post-coalescence mergers and a control sample with the same stellar mass, environmental density, and redshift. In particular, we use a pure sample of 445 post-coalescence mergers from the recent visually-confirmed post-coalescence merger sample identified by Bickley et al. and find that post-coalescence mergers are on average younger than control galaxies for $10<\log (M_\star/\mathrm{M}_\odot)<11$. The difference in age from matched controls is up to 1.5 Gyr, highest for lower stellar mass galaxies. We forward model this difference using parametric star formation histories, accounting for the pre-coalescence inspiral phase of enhanced star formation using close pair data, and a final additive burst of star formation at coalescence. We find a best-fitting stellar mass burst fraction of $f_\mathrm{burst}=ΔM_\star/M_{\star,\mathrm{merger}}=0.18 \pm 0.02$ for $10<\log (M_\star/\mathrm{M}_\odot)<11$ galaxies, with no evidence of a trend in stellar mass. The modeled burst fraction is robust to choice of parametric star formation history, as well as differences in burst duration. The result appears consistent with some prior observationally-derived values, but is significantly higher than that found in hydrodynamical simulations. Using published Luminous InfraRed Galaxy (LIRG) star formation rates, we find a burst duration increasing with stellar mass, from $120-250$ Myr. A comparison to published cold gas measurements indicates there is enough molecular gas available in very close pairs to fuel the burst. Additionally, given our stellar mass burst estimate, the predicted cold gas fraction remaining after the burst is consistent with observed post-coalescence mergers. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS; 12 pages, 9 figures

Gathering Galaxy Distances in Abundance with Roman Wide-Area Data

Authors: John P. Blakeslee, Michele Cantiello, Michael J. Hudson, Laura Ferrarese, Nandini Hazra, Joseph B. Jensen, Eric W. Peng, Gabriella Raimondo

Abstract: The extragalactic distance scale is fundamental to our understanding of astrophysics and cosmology. In recent years, the surface brightness fluctuation (SBF) method, applied in the near-IR, has proven especially powerful for measuring galaxy distances, first with HST and now with a new JWST program to calibrate the method directly from the tip of the red giant branch (TRGB). So far, however, the d… ▽ More The extragalactic distance scale is fundamental to our understanding of astrophysics and cosmology. In recent years, the surface brightness fluctuation (SBF) method, applied in the near-IR, has proven especially powerful for measuring galaxy distances, first with HST and now with a new JWST program to calibrate the method directly from the tip of the red giant branch (TRGB). So far, however, the distances from space have been gathered slowly, one or two at a time. With the Roman Space Telescope, we have the opportunity to measure uniformly high-quality SBF distances to thousands of galaxies out to hundreds of Mpc. The impact of these data on cosmology and galaxy studies depends on the specifics of the survey, including the filter selection, exposure depth, and (especially) the sky coverage. While the baseline HLWAS survey in four filters plus the grism would yield useful data, the impact would be limited by the relatively small area. A more optimal approach would concentrate on the most efficient passband (F146), adopt an exposure time sufficient to measure good quality distances well out into the Hubble flow, and then maximize the sky coverage within the total time constraints. Grism observations over the same area can provide the needed information on redshifts and spectral energy distributions for compact sources, while colors for larger objects can be obtained from lower resolution surveys. The proposed plan will enable accurate determination of the physical properties of thousands of nearby galaxies, an independent measure of the Hubble constant $H_0$ with negligible statistical error, and competitive constraints on $S_8{\,=\,}σ_8(Ω_m/0.3)^{0.5}$. The resulting data set will be a phenomenal resource for a wide range of studies in astrophysics and cosmology. △ Less

Submitted 13 July, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 5 pages; submitted to the call for Roman Core Community Survey white papers (revised version corrects a couple references)

Constraining quenching timescales in galaxy clusters by forward-modelling stellar ages and quiescent fractions in projected phase space

Authors: Andrew M. M. Reeves, Michael J. Hudson, Kyle A. Oman

Abstract: We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in $\log(M_{\mathrm{vir}}/\mathrm{M}_{\odot})>14$ galaxy clusters at $z\sim 0$. We find the average deviation in MWA from the MWA-$M_\star$ relation depends on position in PPS, with a max… ▽ More We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in $\log(M_{\mathrm{vir}}/\mathrm{M}_{\odot})>14$ galaxy clusters at $z\sim 0$. We find the average deviation in MWA from the MWA-$M_\star$ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of $\sim 1$ Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the UniverseMachine model. We find total quenching times of $t_\mathrm{Q}=3.7\pm 0.4$ Gyr and $t_\mathrm{Q}=4.0\pm 0.2$ Gyr after first pericentre, for $9<\log(M_{\star}/\mathrm{M}_{\odot})<10$ and $10<\log(M_{\star}/\mathrm{M}_{\odot})<10.5$ galaxies, respectively. By using MWAs, we break the degeneracy in time of quenching onset and timescale of star formation rate (SFR) decline. We find that time of quenching onset relative to pericentre is $t_{\mathrm{delay}}=3.5^{+0.6}_{-0.9}$ Gyr and $t_{\mathrm{delay}}=-0.3^{+0.8}_{-1.0}$ Gyr for our lower and higher stellar mass bins, respectively, and exponential SFR suppression timescales are $τ_{\mathrm{env}}\leq 1.0$ Gyr and $τ_{\mathrm{env}}\sim 2.3$ Gyr for our lower and higher stellar mass bins, respectively. Stochastic star formation histories remove the need for rapid infall quenching to maintain the bimodality in the SFR of cluster galaxies; the depth of the green valley prefers quenching onsets close to first pericentre and a longer quenching envelope, in slight tension with the MWA-driven results. Taken together these results suggest that quenching begins close to, or just after pericentre, but the timescale for quenching to be fully complete is much longer and therefore ram-pressure stripping is not complete on first pericentric passage. △ Less

Submitted 7 June, 2024; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS (April 1 2023), 21 pages, 13 figures

The shape of dark matter haloes: results from weak lensing in the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS)

Authors: Bailey Robison, Michael J. Hudson, Jean-Charles Cuillandre, Thomas Erben, Sébastien Fabbro, Raphaël Gavazzi, Axel Guinot, Stephen Gwyn, Hendrik Hildebrandt, Martin Kilbinger, Alan McConnachie, Lance Miller, Isaac Spitzer, Ludovic van Waerbeke

Abstract: Cold dark matter haloes are expected to be triaxial, and so appear elliptical in projection. We use weak gravitational lensing from the Canada-France Imaging Survey (CFIS) component of the Ultraviolet-Near Infrared Optical Northern Survey (UNIONS) to measure the ellipticity of the dark matter haloes around Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Data Release 7 (DR7) and from… ▽ More Cold dark matter haloes are expected to be triaxial, and so appear elliptical in projection. We use weak gravitational lensing from the Canada-France Imaging Survey (CFIS) component of the Ultraviolet-Near Infrared Optical Northern Survey (UNIONS) to measure the ellipticity of the dark matter haloes around Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Data Release 7 (DR7) and from the CMASS and LOWZ samples of the Baryon Oscillation Spectroscopic Survey (BOSS), assuming their major axes are aligned with the stellar light. We find that DR7 LRGs with masses $M \sim 2.7\times10^{13} \mathrm{M}_{\odot}/h$ have halo ellipticities $e=0.46\pm0.10$. Expressed as a fraction of the galaxy ellipticity, we find $f_h = 2.2\pm0.6$. For BOSS LRGs, the detection is of marginal significance: $e = 0.20\pm0.10$ and $f_h=0.7\pm0.7$. These results are in agreement with other measurements of halo ellipticity from weak lensing and, taken together with previous results, suggest an increase of halo ellipticity of $0.10\pm0.06$ per decade in halo mass. This trend agrees with the predictions from hydrodynamical simulations, which find that at higher halo masses, not only do dark matter haloes become more elliptical, but that the misalignment between major axis of the stellar light in the central galaxy and that of the dark matter decreases. △ Less

Submitted 19 May, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: MNRAS, accepted

Star formation characteristics of CNN-identified post-mergers in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

Authors: Robert W. Bickley, Sara L. Ellison, David R. Patton, Connor Bottrell, Stephen Gwyn, Michael J. Hudson

Abstract: The importance of the post-merger epoch in galaxy evolution has been well-documented, but post-mergers are notoriously difficult to identify. While the features induced by mergers can sometimes be distinctive, they are frequently missed by visual inspection. In addition, visual classification efforts are highly inefficient because of the inherent rarity of post-mergers (~1% in the low-redshift Uni… ▽ More The importance of the post-merger epoch in galaxy evolution has been well-documented, but post-mergers are notoriously difficult to identify. While the features induced by mergers can sometimes be distinctive, they are frequently missed by visual inspection. In addition, visual classification efforts are highly inefficient because of the inherent rarity of post-mergers (~1% in the low-redshift Universe), and non-parametric statistical merger selection methods do not account for the diversity of post-mergers or the environments in which they appear. To address these issues, we deploy a convolutional neural network (CNN) which has been trained and evaluated on realistic mock observations of simulated galaxies from the IllustrisTNG simulations, to galaxy images from the Canada France Imaging Survey (CFIS), which is part of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS). We present the characteristics of the galaxies with the highest CNN-predicted post-merger certainties, as well as a visually confirmed subset of 699 post-mergers. We find that post-mergers with high CNN merger probabilities (p(x)>0.8) have an average star formation rate that is 0.1 dex higher than a mass- and redshift-matched control sample. The SFR enhancement is even greater in the visually confirmed post-merger sample, a factor of two higher than the control sample. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 14 pages, 9 figures, visually confirmed post-merger catalog available via MNRAS

UNIONS: The impact of systematic errors on weak-lensing peak counts

Authors: Emma Ayçoberry, Virginia Ajani, Axel Guinot, Martin Kilbinger, Valeria Pettorino, Samuel Farrens, Jean-Luc Starck, Raphaël Gavazzi, Michael J. Hudson

Abstract: UNIONS is an ongoing deep photometric multi-band survey of the Northern sky. As part of UNIONS, CFIS provides r-band data which we use to study weak-lensing peak counts for cosmological inference. We assess systematic effects for weak-lensing peak counts and their impact on cosmological parameters for the UNIONS survey. In particular, we present results on local calibration, metacalibration shear… ▽ More UNIONS is an ongoing deep photometric multi-band survey of the Northern sky. As part of UNIONS, CFIS provides r-band data which we use to study weak-lensing peak counts for cosmological inference. We assess systematic effects for weak-lensing peak counts and their impact on cosmological parameters for the UNIONS survey. In particular, we present results on local calibration, metacalibration shear bias, baryonic feedback, the source galaxy redshift estimate, intrinsic alignment, and the cluster member dilution. For each uncertainty and systematic effect, we describe our mitigation scheme and the impact on cosmological parameter constraints. We obtain constraints on cosmological parameters from MCMC using CFIS data and MassiveNuS N-body simulations as a model for peak counts statistics. Depending on the calibration (local versus global, and the inclusion of the residual multiplicative shear bias), the mean matter density parameter $Ω_m$ can shift up to $-0.024$ ($-0.5σ$). We also see that including baryonic corrections can shift $Ω_m$ by $+0.027$ ($+0.5 σ$) with respect to the DM-only simulations. Reducing the impact of the intrinsic alignment and cluster member dilution through signal-to-noise cuts can lead to a shift in $Ω_m$ of $+0.027$ ($+0.5 σ$). Finally, with a mean redshift uncertainty of $Δ\bar{z} = 0.03$, we see that the shift of $Ω_m$ ($+0.001$ which corresponds to $+0.02 σ$) is not significant. This paper investigates for the first time with UNIONS weak-lensing data and peak counts the impact of systematic effects. The value of $Ω_m$ is the most impacted and can shift up to $\sim 0.03$ which corresponds to $0.5σ$ depending on the choices for each systematics. We expect constraints to become more reliable with future (larger) data catalogues, for which the current pipeline will provide a starting point. △ Less

Submitted 22 March, 2023; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: 19 pages, 19 figures

Journal ref: A&A 671, A17 (2023)

ShapePipe: a new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data

Authors: Axel Guinot, Martin Kilbinger, Samuel Farrens, Austin Peel, Arnau Pujol, Morgan Schmitz, Jean-Luc Starck, Thomas Erben, Raphael Gavazzi, Stephen Gwyn, Michael J. Hudson, Hendrik Hiledebrandt, Tobias Liaudat, Lance Miller, Isaac Spitzer, Ludovic Van Waerbeke, Jean-Charles Cuillandre, Sébastien Fabbro, Alan McConnachie

Abstract: UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking $r$-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS $r$-band data over an area spanni… ▽ More UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking $r$-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS $r$-band data over an area spanning 1700 deg$^2$ of the sky. We create a catalogue with measured shapes for 40 million galaxies, corresponding to an effective density of 6.8 galaxies per square arcminute, and demonstrate a low level of systematic biases. This work serves as the basis for further cosmological studies using the full UNIONS survey of 4800 deg$^2$ when completed. Here we present ShapePipe, a newly developed weak-lensing pipeline. This pipeline makes use of state-of-the-art methods such as Ngmix for accurate galaxy shape measurement. Shear calibration is performed with metacalibration. We carry out extensive validation tests on the Point Spread Function (PSF), and on the galaxy shapes. In addition, we create realistic image simulations to validate the estimated shear. We quantify the PSF model accuracy and show that the level of systematics is low as measured by the PSF residuals. Their effect on the shear two-point correlation function is sub-dominant compared to the cosmological contribution on angular scales <100 arcmin. The additive shear bias is below 5x$10^{-4}$, and the residual multiplicative shear bias is at most $10^{-3}$ as measured on image simulations. Using COSEBIs we show that there are no significant B-modes present in second-order shear statistics. We present convergence maps and see clear correlations of the E-mode with known cluster positions. We measure the stacked tangential shear profile around Planck clusters at a significance higher than $4σ$. △ Less

Submitted 10 April, 2022; originally announced April 2022.

Journal ref: A&A 666, A162 (2022)

A Universal Profile for Stacked Filaments from Cold Dark Matter Simulations

Authors: Tianyi Yang, Michael J. Hudson, Niayesh Afshordi

Abstract: 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), p… ▽ More 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. △ Less

Submitted 7 September, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 15 pages. 11 figures. Accepted by MNRAS. Comments welcome

Reconstructing dark matter distribution with peculiar velocities: Bayesian forward modelling with corrections for inhomogeneous Malmquist bias

Authors: Supranta S. Boruah, Guilhem Lavaux, Michael J. Hudson

Abstract: We present a forward-modelled velocity field reconstruction algorithm that performs the reconstruction of the mass density field using only peculiar velocity data. Our method consistently accounts for the inhomogeneous Malmquist bias using analytic integration along the line-of-sight. By testing our method on a simulation, we show that our method gives an unbiased reconstruction of the velocity fi… ▽ More We present a forward-modelled velocity field reconstruction algorithm that performs the reconstruction of the mass density field using only peculiar velocity data. Our method consistently accounts for the inhomogeneous Malmquist bias using analytic integration along the line-of-sight. By testing our method on a simulation, we show that our method gives an unbiased reconstruction of the velocity field. We show that not accounting for the inhomogeneous Malmquist bias can lead to significant biases in the forward-modelled reconstructions. We applied our method to a peculiar velocity data set consisting of the SFI++ and 2MTF Tully-Fisher catalogues and the A2 supernovae compilation, thus obtaining a novel velocity reconstruction in the local Universe. Our velocity reconstructions have a cosmological power spectrum consistent with the theoretical expectation. Furthermore, we obtain a full description of the uncertainties on reconstruction through samples of the posterior distribution. We validate our velocity reconstruction of the local Universe by comparing it to an independent reconstruction using the 2M++ galaxy catalogue, obtaining good agreement between the two reconstructions. Using Bayesian model comparison, we find that our velocity model performs better than the adaptive kernel smoothed velocity with the same peculiar velocity data. However, our velocity model does not perform as well as the velocity reconstruction from the 2M++ galaxy catalogue, due to the sparse and noisy nature of the peculiar velocity tracer samples. The method presented here provides a way to include peculiar velocity data in initial condition reconstruction frameworks. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 15 pages, 12 figures, To be submitted to MNRAS

Ram Pressure Candidates in UNIONS

Authors: Ian D. Roberts, Laura C. Parker, Stephen Gwyn, Michael J. Hudson, Raymond Carlberg, Alan McConnachie, Jean-Charles Cuillandre, Kenneth C. Chambers, Pierre-Alain Duc, Hisanori Furusawa, Raphael Gavazzi, Vanessa Hill, Mark E. Huber, Rodrigo Ibata, Martin Kilbinger, Simona Mei, Yannick Mellier, Satoshi Miyazaki, Masamune Oguri, Richard J. Wainscoat

Abstract: We present a search for disturbed, candidate ram pressure stripping galaxies across more than 50 spectroscopically selected SDSS groups and clusters. Forty-eight ram pressure candidates are visually identified in these systems using high quality UNIONS imaging from the Canada-France Hawaii Telescope, covering ~6200 and ~2800 square degrees in the u- and r-bands respectively. Ram pressure candidate… ▽ More We present a search for disturbed, candidate ram pressure stripping galaxies across more than 50 spectroscopically selected SDSS groups and clusters. Forty-eight ram pressure candidates are visually identified in these systems using high quality UNIONS imaging from the Canada-France Hawaii Telescope, covering ~6200 and ~2800 square degrees in the u- and r-bands respectively. Ram pressure candidates are found in groups and clusters spanning a wide range in halo mass and include ~30 ram pressure candidates in the group regime ($M_h < 10^{14}$). The observed frequency of ram pressure candidates shows substantial scatter with group/cluster mass, but on average is larger in clusters ($M_h > 10^{14}\,M_\odot$) than groups ($M_h < 10^{14}\,M_\odot$) by a factor of ~2. We find that ram pressure candidates are most commonly low-mass galaxies and have enhanced star formation rates relative to star-forming field galaxies. The enhancement in star formation is largely independent of galaxy mass and strongest for galaxies in clusters. As a result of the large survey footprint and excellent image quality from UNIONS, we are able to identify disturbed galaxies, potentially affected by ram pressure stripping, across a wide range of host environment. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 12 pages, 9 figures, 1 appendix, accepted for publication in MNRAS

Strong lensing in UNIONS: Toward a pipeline from discovery to modeling

Authors: E. Savary, K. Rojas, M. Maus, B. Clément, F. Courbin, R. Gavazzi, J. H. H. Chan, C. Lemon, G. Vernardos, R. Cañameras, S. Schuldt, S. H. Suyu, J. -C. Cuillandre, S. Fabbro, S. Gwyn, M. J. Hudson, M. Kilbinger, D. Scott, C. Stone

Abstract: We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We designed a convolutional neural network (CNN) committee that we applied to a selection of 2 344 002 exquisite-seeing $r$-band images of color-selected luminous red galaxies (LRGs). Our classification uses a realistic training set where the lensing gala… ▽ More We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We designed a convolutional neural network (CNN) committee that we applied to a selection of 2 344 002 exquisite-seeing $r$-band images of color-selected luminous red galaxies (LRGs). Our classification uses a realistic training set where the lensing galaxies and the lensed sources are both taken from real data, namely the CFIS $r$-band images themselves and the Hubble Space Telescope (HST). A total of 9 460 candidates obtain a score above 0.5 with the CNN committee. After a visual inspection of the candidates, we find a total of 133 lens candidates, of which 104 are completely new. The set of false positives mainly contains ring, spiral, and merger galaxies, and to a lesser extent galaxies with nearby companions. We classify 32 of the lens candidates as secure lenses and 101 as maybe lenses. For the 32 highest quality lenses, we also fit a singular isothermal ellipsoid mass profile with external shear along with an elliptical Sersic profile for the lens and source light. This automated modeling step provides distributions of properties for both sources and lenses that have Einstein radii in the range $0.5\arcsec<θ_E<2.5\arcsec$. Finally, we introduce a new lens and/or source single-band deblending algorithm based on auto-encoder representation of our candidates. This is the first time an end-to-end lens-finding and modeling pipeline is assembled together, in view of future lens searches in a single band, as will be possible with Euclid. △ Less

Submitted 24 September, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 29 pages, 21 figures, accepted by A&A, in press

Journal ref: A&A 666, A1 (2022)

Discovery of Strongly Lensed Quasars in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

Authors: J. H. H. Chan, C. Lemon, F. Courbin, R. Gavazzi, B. Clément, M. Millon, E. Paic, K. Rojas, E. Savary, G. Vernardos, J. -C. Cuillandre, S. Fabbro, S. Gwyn, M. J. Hudson, M. Kilbinger, A. McConnachie

Abstract: We report the discovery of five new doubly-imaged lensed quasars from the first 2500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either Gaia pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6'' median-seeing $r$-band imaging of CFIS to confi… ▽ More We report the discovery of five new doubly-imaged lensed quasars from the first 2500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either Gaia pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6'' median-seeing $r$-band imaging of CFIS to confirm the presence of multiple point sources with similar colour of $u-r$, via convolution of the Laplacian of the point spread function. Requiring similar-colour point sources with flux ratios less than 2.5 mag in $r$-band, reduces the number of candidates from 256314 to 7815. After visual inspection we obtain 30 high-grade candidates, and prioritise spectroscopic follow-up for those showing signs of a lensing galaxy upon subtraction of the point sources. We obtain long-slit spectra for 18 candidates with ALFOSC on the 2.56-m Nordic Optical Telescope (NOT), confirming five new doubly lensed quasars with $1.21<z<3.36$ and angular separations from 0.8'' to 2.5''. One additional system is a probable lensed quasar based on the CFIS imaging and existing SDSS spectrum. We further classify six objects as nearly identical quasars -- still possible lenses but without the detection of a lensing galaxy. Given our recovery rate ($83\%$) of existing optically bright lenses within the CFIS footprint, we expect that a similar strategy, coupled with $u-r$ colour-selection from CFIS alone, will provide an efficient and complete discovery of small-separation lensed quasars of source redshifts below $z=2.7$ within the CFIS $r$-band magnitude limit of 24.1 mag. △ Less

Submitted 1 December, 2021; v1 submitted 18 October, 2021; originally announced October 2021.

Comments: 15 pages, 8 figures

Journal ref: A&A 659, A140 (2022)

New Constraints on Anisotropic Expansion from Supernovae Type Ia

Authors: W. Rahman, R. Trotta, S. S. Boruah, M. J. Hudson, D. A. van Dyk

Abstract: We re-examine the contentious question of constraints on anisotropic expansion from Type Ia supernovae (SNIa) in the light of a novel determination of peculiar velocities, which are crucial to test isotropy with supernovae out to distances $\lesssim 200/h$ Mpc. We re-analyze the Joint Light-Curve Analysis (JLA) Supernovae (SNe) data, improving on previous treatments of peculiar velocity correction… ▽ More We re-examine the contentious question of constraints on anisotropic expansion from Type Ia supernovae (SNIa) in the light of a novel determination of peculiar velocities, which are crucial to test isotropy with supernovae out to distances $\lesssim 200/h$ Mpc. We re-analyze the Joint Light-Curve Analysis (JLA) Supernovae (SNe) data, improving on previous treatments of peculiar velocity corrections and their uncertainties (both statistical and systematic) by adopting state-of-the-art flow models constrained independently via the 2M$++$ galaxy redshift compilation. We also introduce a novel procedure to account for colour-based selection effects, and adjust the redshift of low-$z$ SNe self-consistently in the light of our improved peculiar velocity model. We adopt the Bayesian hierarchical model \texttt{BAHAMAS} to constrain a dipole in the distance modulus in the context of the $Λ$CDM model and the deceleration parameter in a phenomenological Cosmographic expansion. We do not find any evidence for anisotropic expansion, and place a tight upper bound on the amplitude of a dipole, $|D_μ| < 5.93 \times 10^{-4}$ (95\% credible interval) in a $Λ$CDM setting, and $|D_{q_0}| < 6.29 \times 10^{-2}$ in the Cosmographic expansion approach. Using Bayesian model comparison, we obtain posterior odds in excess of 900:1 (640:1) against a constant-in-redshift dipole for $Λ$CDM (the Cosmographic expansion). In the isotropic case, an accelerating universe is favoured with odds of $\sim 1100:1$ with respect to a decelerating one. △ Less

Submitted 28 April, 2022; v1 submitted 27 August, 2021; originally announced August 2021.

Comments: 25 Pages, 14 figures, submitted to MNRAS

Peculiar-velocity cosmology with Types Ia and II supernovae

Authors: Benjamin E. Stahl, Thomas de Jaeger, Supranta S. Boruah, WeiKang Zheng, Alexei V. Filippenko, Michael J. Hudson

Abstract: We present the Democratic Samples of Supernovae (DSS), a compilation of 775 low-redshift Type Ia and II supernovae (SNe Ia & II), of which 137 SN Ia distances are derived via the newly developed snapshot distance method. Using the objects in the DSS as tracers of the peculiar-velocity field, we compare against the corresponding reconstruction from the 2M++ galaxy redshift survey. Our analysis -- w… ▽ More We present the Democratic Samples of Supernovae (DSS), a compilation of 775 low-redshift Type Ia and II supernovae (SNe Ia & II), of which 137 SN Ia distances are derived via the newly developed snapshot distance method. Using the objects in the DSS as tracers of the peculiar-velocity field, we compare against the corresponding reconstruction from the 2M++ galaxy redshift survey. Our analysis -- which takes special care to properly weight each DSS subcatalogue and cross-calibrate the relative distance scales between them -- results in a measurement of the cosmological parameter combination $fσ_8 = 0.390_{-0.022}^{+0.022}$ as well as an external bulk flow velocity of $195_{-23}^{+22}$ km s$^{-1}$ in the direction $(\ell, b) = (292_{-7}^{+7}, -6_{-4}^{+5})$ deg, which originates from beyond the 2M++ reconstruction. Similarly, we find a bulk flow of $245_{-31}^{+32}$ km s$^{-1}$ toward $(\ell, b) = (294_{-7}^{+7}, 3_{-5}^{+6})$ deg on a scale of $\sim 30 h^{-1}$ Mpc if we ignore the reconstructed peculiar-velocity field altogether. Our constraint on $fσ_8$ -- the tightest derived from SNe to date (considering only statistical error bars), and the only one to utilise SNe II -- is broadly consistent with other results from the literature. We intend for our data accumulation and treatment techniques to become the prototype for future studies that will exploit the unprecedented data volume from upcoming wide-field surveys. △ Less

Submitted 17 May, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Convolutional neural network identification of galaxy post-mergers in UNIONS using IllustrisTNG

Authors: Robert W. Bickley, Connor Bottrell, Maan H. Hani, Sara L. Ellison, Hossen Teimoorinia, Kwang Moo Yi, Scott Wilkinson, Stephen Gwyn, Michael J. Hudson

Abstract: The Canada-France Imaging Survey (CFIS) will consist of deep, high-resolution r-band imaging over ~5000 square degrees of the sky, representing a first-rate opportunity to identify recently-merged galaxies. Due to the large number of galaxies in CFIS, we investigate the use of a convolutional neural network (CNN) for automated merger classification. Training samples of post-merger and isolated gal… ▽ More The Canada-France Imaging Survey (CFIS) will consist of deep, high-resolution r-band imaging over ~5000 square degrees of the sky, representing a first-rate opportunity to identify recently-merged galaxies. Due to the large number of galaxies in CFIS, we investigate the use of a convolutional neural network (CNN) for automated merger classification. Training samples of post-merger and isolated galaxy images are generated from the IllustrisTNG simulation processed with the observational realism code RealSim. The CNN's overall classification accuracy is 88 percent, remaining stable over a wide range of intrinsic and environmental parameters. We generate a mock galaxy survey from IllustrisTNG in order to explore the expected purity of post-merger samples identified by the CNN. Despite the CNN's good performance in training, the intrinsic rarity of post-mergers leads to a sample that is only ~6 percent pure when the default decision threshold is used. We investigate trade-offs in purity and completeness with a variable decision threshold and find that we recover the statistical distribution of merger-induced star formation rate enhancements. Finally, the performance of the CNN is compared with both traditional automated methods and human classifiers. The CNN is shown to outperform Gini-M20 and asymmetry methods by an order of magnitude in post-merger sample purity on the mock survey data. Although the CNN outperforms the human classifiers on sample completeness, the purity of the post-merger sample identified by humans is frequently higher, indicating that a hybrid approach to classifications may be an effective solution to merger classifications in large surveys. △ Less

Submitted 18 March, 2021; v1 submitted 16 March, 2021; originally announced March 2021.

Comments: 21 pages, 19 figures, 2 tables, Accepted for publication in MNRAS

Evolution of subhalo orbits in a smoothly-growing host halo potential

Authors: Go Ogiya, James E. Taylor, Michael J. Hudson

Abstract: The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other re… ▽ More The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other relevant mechanisms, we study subhalo orbital evolution using numerical calculations in which subhaloes are modelled as massless particles orbiting in a time-varying spherical potential. We find that the radial action of the subhalo orbit decreases over the first few orbits, indicating that the response to the growth of the host halo is not adiabatic during this phase. The subhalo orbits can shrink by a factor of $\sim$1.5 in this phase. Subsequently, the radial action is well conserved and orbital contraction slows down. We propose a model accurately describing the orbital evolution. Given these results, we consider the spatial distribution of the population of subhaloes identified in high-resolution cosmological simulations. We find that it is consistent with this population having been accreted at z < 3, indicating that any subhaloes accreted earlier are unresolved in the simulations. We also discuss tidal stripping as a formation scenario for NGC1052-DF2, an ultra diffuse galaxy significantly lacking DM, and find that its expected DM mass could be consistent with observational constraints if its progenitor was accreted early enough, z > 1.5, although it should still be a relatively rare object. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Comments: 15 pages, 11 figures, 1 table, accepted for publication in MNRAS

Assessing the accuracy of cosmological parameters estimated from velocity -- density comparisons via simulations

Authors: Amber M. Hollinger, Michael J. Hudson

Abstract: A promising method for measuring the cosmological parameter combination fsigma_8 is to compare observed peculiar velocities with peculiar velocities predicted from a galaxy density field using perturbation theory. We use N-body simulations and semi-analytic galaxy formation models to quantify the accuracy and precision of this method. Specifically, we examine a number of technical aspects, includi… ▽ More A promising method for measuring the cosmological parameter combination fsigma_8 is to compare observed peculiar velocities with peculiar velocities predicted from a galaxy density field using perturbation theory. We use N-body simulations and semi-analytic galaxy formation models to quantify the accuracy and precision of this method. Specifically, we examine a number of technical aspects, including the optimal smoothing length applied to the density field, the use of dark matter halos or galaxies as tracers of the density field, the effect of noise in the halo mass estimates or in the stellar-to-halo mass relation, and the effect of finite survey volumes. We find that for a Gaussian smoothing of 4 Mpc/h, the method has only small systematic biases at the level of 5%. Cosmic variance affects current measurements at the 5% level due to the volume of current redshift data sets. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 10 pages, 10 figures. Accepted for publication in MNRAS

Peculiar velocities in the local Universe: comparison of different models and the implications for $H_0$ and dark matter

Authors: Supranta S. Boruah, Michael J. Hudson, Guilhem Lavaux

Abstract: When measuring the value of the Hubble parameter, $H_0$, it is necessary to know the recession velocity free of the effects of peculiar velocities. In this work, we study different models of peculiar velocity in the local Universe. In particular, we compare models based on density reconstruction from galaxy redshift surveys and kernel smoothing of peculiar velocity data. The velocity field from th… ▽ More When measuring the value of the Hubble parameter, $H_0$, it is necessary to know the recession velocity free of the effects of peculiar velocities. In this work, we study different models of peculiar velocity in the local Universe. In particular, we compare models based on density reconstruction from galaxy redshift surveys and kernel smoothing of peculiar velocity data. The velocity field from the density reconstruction is obtained using the 2M++ galaxy redshift compilation, which is compared to two adaptive kernel-smoothed velocity fields: the first obtained from the 6dF Fundamental Plane sample and the other using a Tully-Fisher catalogue obtained by combining SFI++ and 2MTF. We highlight that smoothed velocity fields should be rescaled to obtain unbiased velocity estimates. Comparing the predictions of these models to the observations from a few test sets of peculiar velocity data, obtained from the Second Amendment Supernovae catalogue and the Tully-Fisher catalogues, we find that 2M++ reconstruction provides a better model of the peculiar velocity in the local Universe than the kernel-smoothed peculiar velocity models. We study the impact of peculiar velocities on the measurement of $H_0$ from gravitational waves and megamasers. In doing so, we introduce a probabilistic framework to marginalize over the peculiar velocity corrections along the line-of-sight. For the megamasers, we find $H_0 = 69^{+2.9}_{-2.8}$ km s^{-1} Mpc^{-1} using the 2M++ velocity field. We also study the peculiar velocity of the the galaxy NGC1052-DF2, concluding that a short $\sim$ 13 Mpc distance is not a likely explanation of the anomalously low dark matter fraction of that galaxy. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Submitted to MNRAS

A homogeneous measurement of the delay between the onsets of gas stripping and star formation quenching in satellite galaxies of groups and clusters

Authors: Kyle A. Oman, Yannick M. Bahé, Julia Healy, Kelley M. Hess, Michael J. Hudson, Marc A. W. Verheijen

Abstract: We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster environment on star formation in satellite galaxies. We also consider a model for gas stripping, using the same input supplemented with HI fluxes from the ALFALFA survey. The models are motivated by and tested on the… ▽ More We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster environment on star formation in satellite galaxies. We also consider a model for gas stripping, using the same input supplemented with HI fluxes from the ALFALFA survey. The models are motivated by and tested on the Hydrangea cosmological hydrodynamical simulation suite. We recover the characteristic times when satellite galaxies are stripped and quenched. Stripping in massive ($M_\mathrm{ vir}\sim 10^{14.5}\,\mathrm{M}_\odot$) clusters typically occurs at or just before the first pericentric passage. Lower mass ($\sim10^{13.5}\,\mathrm{M}_\odot$) groups strip their satellites on a significantly longer (by $\sim3\,\mathrm{Gyr}$) timescale. Quenching occurs later: Balmer emission lines typically fade $\sim3.5\,\mathrm{Gyr}$ ($5.5\,\mathrm{Gyr}$) after first pericentre in clusters (groups), followed a few hundred $\mathrm{Myr}$ later by reddenning in $(g-r)$ colour. These `delay timescales' are remarkably constant across the entire satellite stellar mass range probed ($\sim10^{9.5}-10^{11}\,\mathrm{M}_\odot$), a feature closely tied to our treatment of `group pre-processing'. The lowest mass groups in our sample ($\sim10^{12.5}\,\mathrm{M}_\odot$) strip and quench their satellites extremely inefficiently: typical timescales may approach the age of the Universe. Our measurements are qualitatively consistent with the `delayed-then-rapid' quenching scenario advocated for by several other studies, but we find significantly longer delay times. Our combination of a homogeneous analysis and input catalogues yields new insight into the sequence of events leading to quenching across wide intervals in host and satellite mass. △ Less

Submitted 7 June, 2024; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: MNRAS accepted version. 25 pages, 11 figures. Appendices included in the source tarball (Download: Other formats: source)

Joint analysis of 6dFGS and SDSS peculiar velocities for the growth rate of cosmic structure and tests of gravity

Authors: Khaled Said, Matthew Colless, Christina Magoulas, John R. Lucey, Michael J. Hudson

Abstract: Measurement of peculiar velocities by combining redshifts and distance indicators is a powerful way to measure the growth rate of cosmic structure and test theories of gravity at low redshift. Here we constrain the growth rate of structure by comparing observed Fundamental Plane peculiar velocities for 15894 galaxies from the 6dF Galaxy Survey (6dFGS) and Sloan Digital Sky Survey (SDSS) with predi… ▽ More Measurement of peculiar velocities by combining redshifts and distance indicators is a powerful way to measure the growth rate of cosmic structure and test theories of gravity at low redshift. Here we constrain the growth rate of structure by comparing observed Fundamental Plane peculiar velocities for 15894 galaxies from the 6dF Galaxy Survey (6dFGS) and Sloan Digital Sky Survey (SDSS) with predicted velocities and densities from the 2M$++$ redshift survey. We measure the velocity scale parameter $β\equiv {Ω_m^γ}/b = 0.372^{+0.034}_{-0.050}$ and $0.314^{+0.031}_{-0.047}$ for 6dFGS and SDSS respectively, where $Ω_m$ is the mass density parameter, $γ$ is the growth index, and $b$ is the bias parameter normalized to the characteristic luminosity of galaxies, $L^*$. Combining 6dFGS and SDSS we obtain $β= 0.341\pm0.024$, implying that the amplitude of the product of the growth rate and the mass fluctuation amplitude is $fσ_8 = 0.338\pm0.027$ at an effective redshift $z=0.035$. Adopting $Ω_m = 0.315\pm0.007$ as favoured by Planck and using $γ=6/11$ for General Relativity and $γ=11/16$ for DGP gravity, we get $S_8(z=0) = σ_8 \sqrt{Ω_m/0.3} =0.637 \pm 0.054$ and $0.741\pm0.062$ for GR and DGP respectively. This measurement agrees with other low-redshift probes of large scale structure but deviates by more than $3σ$ from the latest Planck CMB measurement. Our results favour values of the growth index $γ> 6/11$ or a Hubble constant $H_0 > 70$\,km\,s$^{-1}$\,Mpc$^{-1}$ or a fluctuation amplitude $σ_8 < 0.8$ or some combination of these. Imminent redshift surveys such as Taipan, DESI, WALLABY, and SKA1-MID will help to resolve this tension by measuring the growth rate of cosmic structure to 1\% in the redshift range $0 < z < 1$. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Comments: 19 pages, 15 figures, 2 tables, accepted for publication in MNRAS

How dark are filaments in the cosmic web?

Authors: Tianyi Yang, Michael J. Hudson, Niayesh Afshordi

Abstract: The cold dark matter model predicts that dark matter haloes are connected by filaments. Direct measurements of the masses and structure of these filaments are difficult, but recently several studies have detected these dark-matter-dominated filaments using weak lensing. Here we study the efficiency of galaxy formation within the filaments by measuring their total mass-to-light ratios and stellar m… ▽ More The cold dark matter model predicts that dark matter haloes are connected by filaments. Direct measurements of the masses and structure of these filaments are difficult, but recently several studies have detected these dark-matter-dominated filaments using weak lensing. Here we study the efficiency of galaxy formation within the filaments by measuring their total mass-to-light ratios and stellar mass fractions. Specifically, we stack pairs of Luminous Red Galaxies (LRGs) with a typical separation on the sky of $8 h^{-1}$ Mpc. We stack background galaxy shapes around pairs to obtain mass maps through weak lensing, and we stack galaxies from the Sloan Digital Sky Survey (SDSS) to obtain maps of light and stellar mass. To isolate the signal from the filament, we construct two matched catalogues of physical and non-physical (projected) LRG pairs, with the same distributions of redshift and separation. We then subtract the two stacked maps. Using LRG pair samples from the BOSS survey at two different redshifts, we find that the evolution of the mass in filament is consistent with the predictions from perturbation theory. The filaments are not entirely dark: their mass-to-light ratios ($M/L = 351\pm137$ in solar units in the $r$-band) and stellar mass fractions ($M_{\rm stellar}/M = 0.0073\pm0.0030$) are consistent with the cosmic values (and with their redshift evolutions) △ Less

Submitted 24 August, 2020; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 13 pages, 8 figures. Accepted by MNRAS. Including more discussions regarding the three-point correlation function. Comments are welcomed

Cosmic flows in the nearby Universe: new peculiar velocities from SNe and cosmological constraints

Authors: Supranta S. Boruah, Michael J. Hudson, Guilhem Lavaux

Abstract: The peculiar velocity field offers a unique way to probe dark matter density field on large scales at low redshifts. In this work, we have compiled a new sample of 465 peculiar velocities from low redshift $(z < 0.067)$ Type Ia supernovae. We compare the reconstructed velocity field derived from the 2M++ galaxy redshift compilation to the supernovae, the SFI++ and the 2MTF Tully-Fisher distance ca… ▽ More The peculiar velocity field offers a unique way to probe dark matter density field on large scales at low redshifts. In this work, we have compiled a new sample of 465 peculiar velocities from low redshift $(z < 0.067)$ Type Ia supernovae. We compare the reconstructed velocity field derived from the 2M++ galaxy redshift compilation to the supernovae, the SFI++ and the 2MTF Tully-Fisher distance catalogues. We used a forward method to jointly infer the distances and the velocities of distance indicators by comparing the observations to the reconstruction. Comparison of the reconstructed peculiar velocity fields to observations allows us to infer the cosmological parameter combination $fσ_8$, and the bulk flow velocity arising from outside the survey volume. The residual bulk flow arising from outside the 2M++ volume is inferred to be $171^{+11}_{-11}$ km s$^{-1}$ in the direction $l=301^{\circ} \pm 4^{\circ}$ and $b=0^{\circ} \pm 3^{\circ}$. We obtain $fσ_{8} = 0.400 \pm 0.017$, equivalent to $S_8 \approx σ_8(Ω_m/0.3)^{0.55}=0.776\pm0.033$, which corresponds to an approximately $ 4\%\,$ statistical uncertainty on the value of $fσ_8$. Our inferred value is consistent with other low redshift results in the literature. △ Less

Submitted 2 October, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: MNRAS accepted version

Neural physical engines for inferring the halo mass distribution function

Authors: Tom Charnock, Guilhem Lavaux, Benjamin D. Wandelt, Supranta Sarma Boruah, Jens Jasche, Michael J. Hudson

Abstract: An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter… ▽ More An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: 12 pages, 5 figures

The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019 edition

Authors: The MSE Science Team, Carine Babusiaux, Maria Bergemann, Adam Burgasser, Sara Ellison, Daryl Haggard, Daniel Huber, Manoj Kaplinghat, Ting Li, Jennifer Marshall, Sarah Martell, Alan McConnachie, Will Percival, Aaron Robotham, Yue Shen, Sivarani Thirupathi, Kim-Vy Tran, Christophe Yeche, David Yong, Vardan Adibekyan, Victor Silva Aguirre, George Angelou, Martin Asplund, Michael Balogh, Projjwal Banerjee , et al. (239 additional authors not shown)

Abstract: (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the sc… ▽ More (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

Cosmology with the MaunaKea Spectroscopic Explorer

Authors: Will J. Percival, Christophe Yèche, Maciej Bilicki, Andreu Font-Ribera, Nimish P. Hathi, Cullan Howlett, Michael J. Hudson, Alan W. McConnachie, Faizan Gohar Mohammad, Jeffrey A. Newman, Nathalie Palanque-Delabrouille, Andrei Variu, Yuting Wang, Michael J. Wilson

Abstract: This document summarizes the science cases related to cosmology studies with the MaunaKea Spectroscopic Explorer (MSE), a highly-multiplexed (4332 fibers), wide FOV (1.5 sq deg), large aperture (11.25 m in diameter), optical/NIR (360nm to 1300nm) facility. The MSE High-z Cosmology Survey is designed to probe a large volume of the Universe with a galaxy density sufficient to measure the extremely-l… ▽ More This document summarizes the science cases related to cosmology studies with the MaunaKea Spectroscopic Explorer (MSE), a highly-multiplexed (4332 fibers), wide FOV (1.5 sq deg), large aperture (11.25 m in diameter), optical/NIR (360nm to 1300nm) facility. The MSE High-z Cosmology Survey is designed to probe a large volume of the Universe with a galaxy density sufficient to measure the extremely-large-scale density fluctuations required to explore primordial non-Gaussianity and therefore inflation. We expect a measurement of the local parameter $f_{NL}$ to a precision $σ(f_{NL}) = 1.8$. Combining the MSE High-z Cosmology Survey data with data from a next generation CMB stage 4 experiment and existing DESI data will provide the first $5σ$ confirmation of the neutrino mass hierarchy from astronomical observations. In addition, the Baryonic Acoustic Oscillations (BAO) observed within the sample will provide measurements of the distance-redshift relationship in six different redshift bins between $z=1.6$ and 4.0, each with an accuracy of $\sim0.6\%$. The simultaneous measurements of Redshift Space Distortions (RSD) constrain the amplitude of the fluctuations, at a level ranging from $1.9\%$ to $3.6\%$. The proposed survey covers 10,000 ${\rm deg}^2$, measuring redshifts for three classes of target objects: Emission Line Galaxies (ELGs) with $1.6<z<2.4$, Lyman Break Galaxies (LBGs) with $2.4<z<4.0$, and quasars $2.1<z<3.5$. The ELGs and LBGs will be used as direct tracers of the underlying density field, while the Lyman-$α$ forests in the quasar spectra will be utilized to probe structure. Exposures of duration 1,800sec will guarantee a redshift determination efficiency of $90\%$ for ELGS and at least $50\%$ for LBGs. The survey will represent 100 nights per year for a 5-year MSE program. Finally, three ideas for additional projects of cosmological interest are proposed. △ Less

Submitted 21 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 18 pages, 4 figures, to be appeared as one chapter in "The Detailed Science Case of the Maunakea Spectroscopic Explorer (MSE)"

The SAMI Galaxy Survey: Quenching of star formation in clusters I. Transition galaxies

Authors: Matt S. Owers, Michael J. Hudson, Kyle A. Oman, Joss Bland-Hawthorn, S. Brough, Julia J. Bryant, Luca Cortese, Warrick J. Couch, Scott M. Croom, Jesse van de Sande, Christoph Federrath, Brent Groves, A. M. Hopkins, J. S. Lawrence, Nuria P. F. Lorente, Richard M. McDermid, Anne M. Medling, Samuel N. Richards, Nicholas Scott, Dan S. Taranu, Charlotte Welker, Sukyoung K. Yi

Abstract: We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}δ$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with ste… ▽ More We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}δ$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with stellar mass ${\rm log(}M_*/M_{\odot})>10$. The HDSGs make up a significant fraction of non-passive cluster galaxies (15%; 17/115) and a smaller fraction (2.0%; 8/387) of the non-passive population in low-density environments. The majority (9/17) of cluster HDSGs show evidence for star formation at their centers, with the HDS regions found in the outer parts of the galaxy. Conversely, the $\rm{H}δ$-strong signal is more evenly spread across the galaxy for the majority (6/8) of HDSGs in low-density environments, and is often associated with emission lines that are not due to star formation. We investigate the location of the HDSGs in the clusters, finding that they are exclusively within 0.6$R_{200}$ of the cluster centre, and have a significantly higher velocity dispersion relative to the cluster population. Comparing their distribution in projected-phase-space to those derived from cosmological simulations indicates that the cluster HDSGs are consistent with an infalling population that have entered the central 0.5$r_{200, 3D}$ cluster region within the last $\sim 1\,$Gyr. In the 8/9 cluster HDSGs with central star formation, the extent of star formation is consistent with that expected of outside-in quenching by ram-pressure stripping. Our results indicate that the cluster HDSGs are currently being quenched by ram-pressure stripping on their first passage through the cluster. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: 44 pages (including 8 pages of appendices), 18 figures, 3 tables. Accepted for publicatioin in ApJ

ATLAS Probe: Breakthrough Science of Galaxy Evolution, Cosmology, Milky Way, and the Solar System

Authors: Yun Wang, Massimo Robberto, Mark Dickinson, Lynne A. Hillenbrand, Wesley Fraser, Peter Behroozi, Jarle Brinchmann, Chia-Hsun Chuang, Andrea Cimatti, Robert Content, Emanuele Daddi, Henry C. Ferguson, Christopher Hirata, Michael J. Hudson, J. Davy Kirkpatrick, Alvaro Orsi, Russell Ryan, Alice Shapley, Mario Ballardini, Robert Barkhouser, James Bartlett, Robert Benjamin, Ranga Chary, Charlie Conroy, Megan Donahue , et al. (15 additional authors not shown)

Abstract: ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, a… ▽ More ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way. △ Less

Submitted 24 January, 2019; v1 submitted 5 February, 2018; originally announced February 2018.

Comments: 46 pages. Accepted by PASA

Journal ref: Publications of the Astronomical Society of Australia (2019), 36, e015, 31 pages

Chemical Mapping of the Milky Way With The Canada-France Imaging Survey: A Non-parametric Metallicity-Distance Decomposition of the Galaxy

Authors: Rodrigo Ibata, Alan McConnachie, Jean-Charles Cuillandre, Nicholas Fantin, Misha Haywood, Nicolas F. Martin, Piere Bergeron, Volker Beckmann, Edouard Bernard, Piercarlo Bonifacio, Elisabetta Caffau, Raymond Carlberg, Patrick Côté, Rémi Cabanac, Scott Chapman, Pierre-Alain Duc, Florence Durret, Benoît Famaey, Sébastien Frabbro, Stephen Gwyn, Francois Hammer, Vanessa Hill, Michael J. Hudson, Ariane Lançon, Geraint Lewis , et al. (19 additional authors not shown)

Abstract: We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicitie… ▽ More We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30°< |b| < 70°$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 21 pages, 17 figures, accepted for publication in ApJ

The Canada-France Imaging Survey: First results from the u-band component

Authors: Rodrigo Ibata, Alan McConnachie, Jean-Charles Cuillandre, Nicholas Fantin, Misha Haywood, Nicolas F. Martin, Piere Bergeron, Volker Beckmann, Edouard Bernard, Piercarlo Bonifacio, Elisabetta Caffau, Raymond Carlberg, Patrick Côté, Rémi Cabanac, Scott Chapman, Pierre-Alain Duc, Florence Durret, Benoît Famaey, Sébastien Frabbro, Stephen Gwyn, Francois Hammer, Vanessa Hill, Michael J. Hudson, Ariane Lançon, Geraint Lewis , et al. (19 additional authors not shown)

Abstract: The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we pre… ▽ More The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 12 pages, 11 figures, accepted for publication in ApJ

The Dependence of Halo Mass on Galaxy Size at Fixed Stellar Mass Using Weak Lensing

Authors: Paul J. L. Charlton, Michael J. Hudson, Michael L. Balogh, Sumeet Khatri

Abstract: Stellar mass has been shown to correlate with halo mass, with non-negligible scatter. The stellar mass-size and luminosity-size relationships of galaxies also show significant scatter in galaxy size at fixed stellar mass. It is possible that, at fixed stellar mass and galaxy colour, the halo mass is correlated with galaxy size. Galaxy-galaxy lensing allows us to measure the mean masses of dark mat… ▽ More Stellar mass has been shown to correlate with halo mass, with non-negligible scatter. The stellar mass-size and luminosity-size relationships of galaxies also show significant scatter in galaxy size at fixed stellar mass. It is possible that, at fixed stellar mass and galaxy colour, the halo mass is correlated with galaxy size. Galaxy-galaxy lensing allows us to measure the mean masses of dark matter haloes for stacked samples of galaxies. We extend the analysis of the galaxies in the CFHTLenS catalogue by fitting single Sérsic surface brightness profiles to the lens galaxies in order to recover half-light radius values, allowing us to determine halo masses for lenses according to their size. Comparing our halo masses and sizes to baselines for that stellar mass yields a differential measurement of the halo mass-galaxy size relationship at fixed stellar mass, defined as $M_{h}(M_{*}) \propto r_{\mathrm{eff}}^η(M_{*})$. We find that on average, our lens galaxies have an $η= 0.42\pm0.12$, i.e. larger galaxies live in more massive dark matter haloes. The $η$ is strongest for high mass luminous red galaxies (LRGs). Investigation of this relationship in hydrodynamical simulations suggests that, at a fixed $M_{*}$, satellite galaxies have a larger $η$ and greater scatter in the $M_{\mathrm{h}}$ and $r_{\mathrm{eff}}$ relationship compared to central galaxies. △ Less

Submitted 16 July, 2017; originally announced July 2017.

Comments: 23 pages, 16 figures

The correlation between the sizes of globular cluster systems and their host dark matter haloes

Authors: Michael J. Hudson, Bailey Robison

Abstract: The sizes of entire systems of globular clusters (GCs) depend on the formation and destruction histories of the GCs themselves, but also on the assembly, merger and accretion history of the dark matter (DM) haloes that they inhabit. Recent work has shown a linear relation between total mass of globular clusters in the globular cluster system and the mass of its host dark matter halo, calibrated fr… ▽ More The sizes of entire systems of globular clusters (GCs) depend on the formation and destruction histories of the GCs themselves, but also on the assembly, merger and accretion history of the dark matter (DM) haloes that they inhabit. Recent work has shown a linear relation between total mass of globular clusters in the globular cluster system and the mass of its host dark matter halo, calibrated from weak lensing. Here we extend this to GC system sizes, by studying the radial density profiles of GCs around galaxies in nearby galaxy groups. We find that radial density profiles of the GC systems are well fit with a de Vaucouleurs profile. Combining our results with those from the literature, we find tight relationship ($\sim 0.2$ dex scatter) between the effective radius of the GC system and the virial radius (or mass) of its host DM halo. The steep non-linear dependence of this relationship ($R_{e, GCS} \propto R_{200}^{2.5 - 3}$) is currently not well understood, but is an important clue regarding the assembly history of DM haloes and of the GC systems that they host. △ Less

Submitted 9 July, 2017; originally announced July 2017.

Comments: Submitted to MNRAS, comments welcome

The Weak Lensing Masses of Filaments between Luminous Red Galaxies

Authors: Seth D. Epps, Michael J. Hudson

Abstract: In the standard model of non-linear structure formation, a cosmic web of dark-matter dominated filaments connects dark matter halos. In this paper, we stack the weak lensing signal of an ensemble of filaments between groups and clusters of galaxies. Specifically, we detect the weak lensing signal, using CFHTLenS galaxy ellipticities, from stacked filaments between SDSS-III/BOSS luminous red galaxi… ▽ More In the standard model of non-linear structure formation, a cosmic web of dark-matter dominated filaments connects dark matter halos. In this paper, we stack the weak lensing signal of an ensemble of filaments between groups and clusters of galaxies. Specifically, we detect the weak lensing signal, using CFHTLenS galaxy ellipticities, from stacked filaments between SDSS-III/BOSS luminous red galaxies (LRGs). As a control, we compare the physical LRG pairs with projected LRG pairs that are more widely separated in redshift space. We detect the excess filament mass density in the projected pairs at the $5σ$ level, finding a mass of $(1.6 \pm 0.3) \times 10^{13} M_{\odot}$ for a stacked filament region 7.1 $h^{-1}$ Mpc long and 2.5 $h^{-1}$ Mpc wide. This filament signal is compared with a model based on the three-point galaxy-galaxy-convergence correlation function, as developed in Clampitt, Jain & Takada (2014), yielding reasonable agreement. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: MNRAS, in press

Journal ref: MNRAS, 468, 2605 (2017)

Dynamical friction in the primordial neutrino sea

Authors: Chiamaka Okoli, Morag I. Scrimgeour, Niayesh Afshordi, Michael J. Hudson

Abstract: Standard big bang cosmology predicts a cosmic neutrino background at $T_ν\simeq 1.95~$K. Given the current neutrino oscillation measurements, we know most neutrinos move at large, but non-relativistic, velocities. Therefore, dark matter haloes moving in the sea of primordial neutrinos form a neutrino wake behind them, which would slow them down, due to the effect of {\it dynamical friction}. In th… ▽ More Standard big bang cosmology predicts a cosmic neutrino background at $T_ν\simeq 1.95~$K. Given the current neutrino oscillation measurements, we know most neutrinos move at large, but non-relativistic, velocities. Therefore, dark matter haloes moving in the sea of primordial neutrinos form a neutrino wake behind them, which would slow them down, due to the effect of {\it dynamical friction}. In this paper, we quantify this effect for realistic haloes, in the context of the halo model of structure formation, and show that it scales as $m_ν^4$ $\times$ relative velocity, and monotonically grows with the halo mass. Galaxy redshift surveys can be sensitive to this effect (at $>3σ$ confidence level, depending on survey properties, neutrino mass and hierarchy) through redshift space distortions (RSD) of distinct galaxy populations. △ Less

Submitted 28 April, 2017; v1 submitted 14 November, 2016; originally announced November 2016.

Comments: 12 pages, 15 figures

Journal ref: Mon Not R Astron Soc (2017) 468 (2): 2164-2175

Satellite quenching timescales in clusters from projected phase space measurements matched to simulated orbits

Authors: Kyle A. Oman, Michael J. Hudson

Abstract: We measure the star formation quenching efficiency and timescale in cluster environments. Our method uses N-body simulations to estimate the probability distribution of possible orbits for a sample of observed SDSS galaxies in and around clusters based on their position and velocity offsets from their host cluster. We study the relationship between their star formation rates and their likely orbit… ▽ More We measure the star formation quenching efficiency and timescale in cluster environments. Our method uses N-body simulations to estimate the probability distribution of possible orbits for a sample of observed SDSS galaxies in and around clusters based on their position and velocity offsets from their host cluster. We study the relationship between their star formation rates and their likely orbital histories via a simple model in which star formation is quenched once a delay time after infall has elapsed. Our orbit library method is designed to isolate the environmental effect on the star formation rate due to a galaxy's present-day host cluster from `pre-processing' in previous group hosts. We find that quenching of satellite galaxies of all stellar masses in our sample ($10^{9}-10^{11.5}\,{\rm M}_\odot$) by massive ($> 10^{13}\,{\rm M}_\odot$) clusters is essentially $100$ per cent efficient. Our fits show that all galaxies quench on their first infall, approximately at or within a Gyr of their first pericentric passage. There is little variation in the onset of quenching from galaxy-to-galaxy: the spread in this time is at most $\sim 2$ Gyr at fixed $M_*$. Higher mass satellites quench earlier, with very little dependence on host cluster mass in the range probed by our sample. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 14 pages, 10 figures, 1 table. MNRAS submitted. Comments welcome

Beyond Sérsic + exponential disc morphologies in the Coma Cluster

Authors: Jacob T. C. G. Head, John R. Lucey, Michael J. Hudson

Abstract: [abridged] We explore the diversity of internal galaxy structures in the Coma cluster across a wide range of luminosities ($-17$\,$>$\,$M_g$\,$>$\,$-22$) and cluster-centric radii ($0$\,$<$\,$r_{\rm{cluster}}$\,$<$\,1.3 $r_{200}$) through analysis of deep Canada-France-Hawaii Telescope $i$ band imaging. We present 2D multi-component decomposition via GALFIT, encompassing a wide range of candidate… ▽ More [abridged] We explore the diversity of internal galaxy structures in the Coma cluster across a wide range of luminosities ($-17$\,$>$\,$M_g$\,$>$\,$-22$) and cluster-centric radii ($0$\,$<$\,$r_{\rm{cluster}}$\,$<$\,1.3 $r_{200}$) through analysis of deep Canada-France-Hawaii Telescope $i$ band imaging. We present 2D multi-component decomposition via GALFIT, encompassing a wide range of candidate model morphologies with up to three photometric components. Particular focus is placed on early-type galaxies with outer discs (i.e. S0s), and deviations from simple (`unbroken') exponential discs. Rigorous filtering ensures that each model component provides a statistically significant improvement to the goodness-of-fit. The majority of Coma cluster members in our sample (478 of 631) are reliably fitted by symmetric structural models. Of these, 134 ($28\%$) are single Sérsic objects, 143 ($30\%$) are well-described by 2 component structures, while 201 ($42\%$) require more complex models. Multi-component Sérsic galaxies resemble compact psuedobulges ($n\sim$\,2, $R_e \sim$\, 4 kpc) surrounded by extended Gaussian-like outer structures ($R_e > 10$ kpc). 11\% of galaxies ($N=52$) feature a break in their outer profiles, indicating `truncated' or `anti-truncated' discs. Beyond the break radius, truncated galaxies are structurally consistent with exponential discs, disfavouring physical truncation as their formation mechanism. Bulge luminosity in anti-truncated galaxies correlates strongly with galaxy luminosity, indicating a bulge-enhancing origin for these systems. Both types of broken disc are found overwhelmingly ($>70\%$) in `barred' galaxies, despite a low measured bar fraction for Coma ($20\pm2\%$). Thus, galaxy bars play an important role in formation of broken disc structures. No strong variation in galaxy structure is detected with projected cluster-centric radius. △ Less

Submitted 28 July, 2015; originally announced July 2015.

Comments: 26 pages, 24 figures. Accepted for publication in MNRAS

CFHTLenS: Weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment

Authors: Tim Schrabback, Stefan Hilbert, Henk Hoekstra, Patrick Simon, Edo van Uitert, Thomas Erben, Catherine Heymans, Hendrik Hildebrandt, Thomas D. Kitching, Yannick Mellier, Lance Miller, Ludovic Van Waerbeke, Philip Bett, Jean Coupon, Liping Fu, Michael J. Hudson, Benjamin Joachimi, Martin Kilbinger, Konrad Kuijken

Abstract: We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observ… ▽ More We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models. △ Less

Submitted 20 October, 2015; v1 submitted 15 July, 2015; originally announced July 2015.

Comments: 21 pages, 3 tables, 9 figures. This replacement matches the version accepted for publication in MNRAS

Journal ref: MNRAS 2015, 454, p1432-1452

Cosmological parameters from the comparison of peculiar velocities with predictions from the 2M++ density field

Authors: Jonathan Carrick, Stephen J. Turnbull, Guilhem Lavaux, Michael J. Hudson

Abstract: Peculiar velocity measurements are the only tool available in the low-redshift Universe for mapping the large-scale distribution of matter and can thus be used to constrain cosmology. Using redshifts from the 2M++ redshift compilation, we reconstruct the density of galaxies within 200 Mpc/h, allowing for the first time good sampling of important superclusters such as the Shapley Concentration. We… ▽ More Peculiar velocity measurements are the only tool available in the low-redshift Universe for mapping the large-scale distribution of matter and can thus be used to constrain cosmology. Using redshifts from the 2M++ redshift compilation, we reconstruct the density of galaxies within 200 Mpc/h, allowing for the first time good sampling of important superclusters such as the Shapley Concentration. We compare the predicted peculiar velocities from 2M++ to Tully-Fisher and SNe peculiar velocities. We find a value of $β^* \equiv Ω_{\rm{m}}^{0.55}/b^* = 0.431 \pm 0.021$, suggesting $Ω_{\rm{m}}^{0.55}σ_{\rm{8,lin}} = 0.401 \pm 0.024$, in good agreement with other probes. The predicted peculiar velocity of the Local Group arising from the 2M++ volume alone is $540 \pm 40$ km/s, towards $l = 268 \pm 4, b= 38 \pm 6$, only $10^\circ$ out of alignment with the Cosmic Microwave Background dipole. To account for velocity contributions arising from sources outside the 2M++ volume, we fit simultaneously for $β^*$ and an external bulk flow in our analysis. We find that an external bulk flow is preferred at the 5.1$σ$ level, and the best fit has a velocity of $159\pm23$ km/s towards $l=304 \pm 11, b = 6 \pm 13$. Finally, the predicted bulk flow of a 50 Mpc/h Gaussian-weighted volume centred on the Local Group is $230 \pm 30$ km/s, in the direction $l=293\pm 8, b = 14 \pm 10$, in agreement with predictions from $Λ$CDM. △ Less

Submitted 17 April, 2015; originally announced April 2015.

Comments: 17 pages. MNRAS, in press. Density and peculiar velocity grids available at http://cosmicflows.uwaterloo.ca or http://cosmicflows.iap.fr

The galaxy-halo connection from a joint lensing, clustering and abundance analysis in the CFHTLenS/VIPERS field

Authors: J. Coupon, S. Arnouts, L. van Waerbeke, T. Moutard, O. Ilbert, E. van Uitert, T. Erben, B. Garilli, L. Guzzo, C. Heymans, H. Hildebrandt, H. Hoekstra, M. Kilbinger, T. Kitching, Y. Mellier, L. Miller, M. Scodeggio, C. Bonnett, E. Branchini, I. Davidzon, G. De Lucia, A. Fritz, L. Fu, P. Hudelot, M. J. Hudson , et al. (9 additional authors not shown)

Abstract: We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplem… ▽ More We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by $\sim60\,000$ secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy lensing and the stellar mass function. We interpret our measurements within the halo occupation distribution (HOD) framework, separating the contributions from central and satellite galaxies. We find that the SHMR for the central galaxies peaks at $M_{\rm h, peak} = 1.9^{+0.2}_{-0.1}\times10^{12} M_{\odot}$ with an amplitude of $0.025$, which decreases to $\sim0.001$ for massive halos ($M_{\rm h} > 10^{14} M_{\odot}$). Compared to central galaxies only, the total SHMR (including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size halos), a result consistent with cluster analyses from the literature based on fully independent methods. After properly accounting for differences in modelling, we have compared our results with a large number of results from the literature up to $z=1$: we find good general agreement, independently of the method used, within the typical stellar-mass systematic errors at low to intermediate mass (${M}_{\star} < 10^{11} M_{\odot}$) and the statistical errors above. We have also compared our SHMR results to semi-analytic simulations and found that the SHMR is tilted compared to our measurements in such a way that they over- (under-) predict star formation efficiency in central (satellite) galaxies. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 31 pages, 18 figures, 4 table. Accepted for publication in MNRAS. Online material available at http://www.cfhtlens.org

Kinematic Bias in Cosmological Distance Measurement

Authors: Nick Kaiser, Michael J. Hudson

Abstract: Recent calculations using non-linear relativistic cosmological perturbation theory show biases in the mean luminosity distance and distance modulus at low redshift. We show that these effects may be understood very simply as a non-relativistic, and purely kinematic, Malmquist-like bias, and we describe how the effect changes if one averages over sources that are limited by apparent magnitude. This… ▽ More Recent calculations using non-linear relativistic cosmological perturbation theory show biases in the mean luminosity distance and distance modulus at low redshift. We show that these effects may be understood very simply as a non-relativistic, and purely kinematic, Malmquist-like bias, and we describe how the effect changes if one averages over sources that are limited by apparent magnitude. This effect is essentially identical to the distance bias from small-scale random velocities that has previously been considered by astronomers, though we find that the standard formula overestimates the homogeneous bias by a factor 2. △ Less

Submitted 5 February, 2015; originally announced February 2015.

Comments: 8 pages, 1 figure, submitted to MNRAS

On the perturbation of the luminosity distance by peculiar motions

Authors: Nick Kaiser, Michael J. Hudson

Abstract: We consider some aspects of the perturbation to the luminosity distance $d(z)$ that are of relevance for SN1a cosmology and for future peculiar velocity surveys at non-negligible redshifts. 1) Previous work has shown that the correction to the lowest order perturbation $δd / d = -δv / c z$ has the peculiar characteristic that it appears to depend on the absolute state of motion of sources, rathe… ▽ More We consider some aspects of the perturbation to the luminosity distance $d(z)$ that are of relevance for SN1a cosmology and for future peculiar velocity surveys at non-negligible redshifts. 1) Previous work has shown that the correction to the lowest order perturbation $δd / d = -δv / c z$ has the peculiar characteristic that it appears to depend on the absolute state of motion of sources, rather than on their motion relative to that of the observer. The resolution of this apparent violation of the equivalence principle is that it is necessary to allow for evolution of the velocities with time, and also, when considering perturbations on the scale of the observer-source separation, to include the gravitational redshift effect. We provide an expression for $δd / d$ that provides a physically consistent way to measure peculiar velocities and determine their impact for SN1a cosmology. 2) We then calculate the perturbation to the redshift as a function of source flux density, which has been proposed as an alternative probe of large-scale motions. We show how the inclusion of surface brightness modulation modifies the relation between $δz(m)$ and the peculiar velocity, and that, while the noise properties of this method might appear promising, the velocity signal is swamped by the effect of galaxy clustering for most scales of interest. 3) We show how, in linear theory, peculiar velocity measurements are biased downwards by the effect of smaller scale motions or by measurement errors (such as in photometric redshifts). Our results nicely explain the effects seen in simulations by Koda et al.\ 2013. We critically examine the prospects for extending peculiar velocity studies to larger scales with near-term future surveys. △ Less

Submitted 25 March, 2015; v1 submitted 23 November, 2014; originally announced November 2014.

Comments: 14 pages, 0 figures, revised version for MNRAS

CFHTLenS: Weak lensing calibrated scaling relations for low mass clusters of galaxies

Authors: K. Kettula, S. Giodini, E. van Uitert, H. Hoekstra, A. Finoguenov, M. Lerchster, T. Erben, C. Heymans, H. Hildebrandt, T. D. Kitching, A. Mahdavi, Y. Mellier, L. Miller, M. Mirkazemi, L. Van Waerbeke, J. Coupon, E. Egami, L. Fu, M. J. Hudson, J. P. Kneib, K. Kuijken, H. J. McCracken, M. J. Pereira, B. Rowe, T. Schrabback , et al. (2 additional authors not shown)

Abstract: We present weak lensing and X-ray analysis of 12 low mass clusters from the CFHTLenS and XMM-CFHTLS surveys. We combine these systems with high-mass systems from CCCP and low-mass systems from COSMOS to obtain a sample of 70 systems, spanning over two orders of magnitude in mass. We measure core-excised Lx-Tx, M-Lx and M-Tx scaling relations and include corrections for observational biases. By pro… ▽ More We present weak lensing and X-ray analysis of 12 low mass clusters from the CFHTLenS and XMM-CFHTLS surveys. We combine these systems with high-mass systems from CCCP and low-mass systems from COSMOS to obtain a sample of 70 systems, spanning over two orders of magnitude in mass. We measure core-excised Lx-Tx, M-Lx and M-Tx scaling relations and include corrections for observational biases. By providing fully bias corrected relations, we give the current limitations for Lx and Tx as cluster mass proxies. We demonstrate that Tx benefits from a significantly lower intrinsic scatter at fixed mass than Lx. By studying the residuals of the bias corrected relations, we show for the first time using weak lensing masses that galaxy groups seem more luminous and warmer for their mass than clusters. This implies a steepening of the M-Lx and M-Tx relations at low masses. We verify the inferred steepening using a different high mass sample from the literature and show that variance between samples is the dominant effect leading to discrepant scaling relations. We divide our sample into subsamples of merging and relaxed systems, and find that mergers may have enhanced scatter in lensing measurements, most likely due to stronger triaxiality and more substructure. For the Lx-Tx relation, which is unaffected by lensing measurements, we find the opposite trend in scatter. We also explore the effects of X-ray cross-calibration and find that Chandra calibration leads to flatter Lx-Tx and M-Tx relations than XMM-Newton. △ Less

Submitted 9 April, 2015; v1 submitted 31 October, 2014; originally announced October 2014.

Comments: Analysis modified to include bias correction. Updated CCCP measurements. Comments are welcome, submitted to MNRAS, 26 pages, 20 figures

CFHTLenS: A Weak Lensing Shear Analysis of the 3D-Matched-Filter Galaxy Clusters

Authors: Jes Ford, Ludovic Van Waerbeke, Martha Milkeraitis, Clotilde Laigle, Hendrik Hildebrandt, Thomas Erben, Catherine Heymans, Henk Hoekstra, Thomas Kitching, Yannick Mellier, Lance Miller, Ami Choi, Jean Coupon, Liping Fu, Michael J. Hudson, Konrad Kuijken, Naomi Robertson, Barnaby Rowe, Tim Schrabback, Malin Velander

Abstract: We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2… ▽ More We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$σ$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$σ$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^β$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $β\sim 1.4 \pm 0.1$, both of which are in 1$σ$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org. △ Less

Submitted 9 December, 2014; v1 submitted 11 September, 2014; originally announced September 2014.

Comments: 3D-MF cluster catalog is NOW AVAILABLE at cfhtlens.org. Magnification-shear mass comparison in Figure 10. 19 pages, 10 figures. Accepted to MNRAS

CFHTLenS: Cosmological constraints from a combination of cosmic shear two-point and three-point correlations

Authors: Liping Fu, Martin Kilbinger, Thomas Erben, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Thomas D. Kitching, Yannick Mellier, Lance Miller, Elisabetta Semboloni, Patrick Simon, Ludovic Van Waerbeke, Jean Coupon, Joachim Harnois-Déraps, Michael J. Hudson, Konrad Kuijken, Barnaby Rowe, Tim Schrabback, Sanaz Vafaei, Malin Velander

Abstract: Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2σ$. The combined constr… ▽ More Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2σ$. The combined constraint on $Σ_8 = σ_8 (Ω_{\rm m}/0.27)^α$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $Ω_{\rm m}$ and $σ_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $Ω_{\rm m}$ and $σ_8$ by 26% for flat $Λ$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $Ω_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function. △ Less

Submitted 22 April, 2014; originally announced April 2014.

Comments: Accepted by MNRAS. 21 pages, 14 figues, 8 tables. The data is available at http://www.cfhtlens.org . The software used for the cosmological analysis can be downloaded from http://cosmopmc.info

Dark Matter Halos in Galaxies and Globular Cluster Populations

Authors: Michael J. Hudson, Gretchen L. Harris, William E. Harris

Abstract: We combine a new, comprehensive database for globular cluster populations in all types of galaxies with a new calibration of galaxy halo masses based entirely on weak lensing. Correlating these two sets of data, we find that the mass ratio $η\equiv M_{GCS}/M_{h}$ (total mass in globular clusters, divided by halo mass) is essentially constant at $\langle η\rangle \sim 4 \times 10^{-5}$, strongly co… ▽ More We combine a new, comprehensive database for globular cluster populations in all types of galaxies with a new calibration of galaxy halo masses based entirely on weak lensing. Correlating these two sets of data, we find that the mass ratio $η\equiv M_{GCS}/M_{h}$ (total mass in globular clusters, divided by halo mass) is essentially constant at $\langle η\rangle \sim 4 \times 10^{-5}$, strongly confirming earlier suggestions in the literature. Globular clusters are the only known stellar population that formed in essentially direct proportion to host galaxy halo mass. The intrinsic scatter in $η$ appears to be at most 0.2 dex; we argue that some of this scatter is due to differing degrees of tidal stripping of the globular cluster systems between central and satellite galaxies. We suggest that this correlation can be understood if most globular clusters form at very early stages in galaxy evolution, largely avoiding the feedback processes that inhibited the bulk of field-star formation in their host galaxies. The actual mean value of $η$ also suggests that about $1/4$ of the \emph{initial} gas mass present in protogalaxies collected into GMCs large enough to form massive, dense star clusters. Finally, our calibration of $\langle η\rangle$ indicates that the halo masses of the Milky Way and M31 are $(1.2\pm0.5)\times 10^{12} M_{\odot}$ and $(3.9\pm1.8)\times 10^{12} M_{\odot}$ respectively. △ Less

Submitted 7 April, 2014; originally announced April 2014.

Comments: 6 pages, ApJL in press

The Phase Space and Stellar Populations of Cluster Galaxies at z ~ 1: Simultaneous Constraints on the Location and Timescale of Satellite Quenching

Authors: Adam Muzzin, R. F. J. van der Burg, Sean L. McGee, Michael Balogh, Marijn Franx, Henk Hoekstra, Michael J. Hudson, Allison Noble, Dan Taranu, Tracy Webb, Gillian Wilson, H. K. C. Yee

Abstract: We investigate the velocity vs. position phase space of z ~ 1 cluster galaxies using a set of 424 spectroscopic redshifts in 9 clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories: quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commo… ▽ More We investigate the velocity vs. position phase space of z ~ 1 cluster galaxies using a set of 424 spectroscopic redshifts in 9 clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories: quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commonly found at small clustercentric radii with high clustercentric velocities, and appear to trace a coherent ``ring" in phase space. Using several zoom simulations of clusters we show that the coherent distribution of the poststarbursts can be reasonably well-reproduced using a simple quenching scenario. Specifically, the phase space is best reproduced if satellite quenching occurs on a rapid timescale (0.1 < tau_{Q} < 0.5 Gyr) after galaxies make their first passage of R ~ 0.5R_{200}, a process that takes a total time of ~ 1 Gyr after first infall. We compare this quenching timescale to the timescale implied by the stellar populations of the poststarburst galaxies and find that the poststarburst spectra are well-fit by a rapid quenching (tau_{Q} = 0.4^{+0.3}_{-0.4} Gyr) of a typical star-forming galaxy. The similarity between the quenching timescales derived from these independent indicators is a strong consistency check of the quenching model. Given that the model implies satellite quenching is rapid, and occurs well within R_{200}, this would suggest that ram-pressure stripping of either the hot or cold gas component of galaxies are the most plausible candidates for the physical mechanism. The high cold gas consumption rates at z ~ 1 make it difficult to determine if hot or cold gas stripping is dominant; however, measurements of the redshift evolution of the satellite quenching timescale and location may be capable of distinguishing between the two. △ Less

Submitted 27 February, 2014; originally announced February 2014.

Comments: 10 pages, 4 figures, submitted to the ApJ