Showing 1–50 of 54 results for author: Harvey, D

A deep-learning algorithm to disentangle self-interacting dark matter and AGN feedback models

Authors: David Harvey

Abstract: Different models of dark matter can alter the distribution of mass in galaxy clusters in a variety of ways. However, so can uncertain astrophysical feedback mechanisms. Here we present a Machine Learning method that ''learns'' how the impact of dark matter self-interactions differs from that of astrophysical feedback in order to break this degeneracy and make inferences on dark matter. We train a… ▽ More Different models of dark matter can alter the distribution of mass in galaxy clusters in a variety of ways. However, so can uncertain astrophysical feedback mechanisms. Here we present a Machine Learning method that ''learns'' how the impact of dark matter self-interactions differs from that of astrophysical feedback in order to break this degeneracy and make inferences on dark matter. We train a Convolutional Neural Network on images of galaxy clusters from hydro-dynamic simulations. In the idealised case our algorithm is 80% accurate at identifying if a galaxy cluster harbours collisionless dark matter, dark matter with $σ_{\rm DM}/m = 0.1$cm$^2/$g or with $σ_{DM}/m = 1$cm$^2$/g. Whilst we find adding X-ray emissivity maps does not improve the performance in differentiating collisional dark matter, it does improve the ability to disentangle different models of astrophysical feedback. We include noise to resemble data expected from Euclid and Chandra and find our model has a statistical error of < 0.01cm$^2$/g and that our algorithm is insensitive to shape measurement bias and photometric redshift errors. This method represents a new way to analyse data from upcoming telescopes that is an order of magnitude more precise and many orders faster, enabling us to explore the dark matter parameter space like never before. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: Accepted Nature Astronomy

Hydrodynamical simulations of merging galaxy clusters: giant dark matter particle colliders, powered by gravity

Authors: Ellen L. Sirks, David Harvey, Richard Massey, Kyle A. Oman, Andrew Robertson, Carlos Frenk, Spencer Everett, Ajay S. Gill, David Lagattuta, Jacqueline McCleary

Abstract: Terrestrial particle accelerators collide charged particles, then watch the trajectory of outgoing debris - but they cannot manipulate dark matter. Fortunately, dark matter is the main component of galaxy clusters, which are continuously pulled together by gravity. We show that galaxy cluster mergers can be exploited as enormous, natural dark matter colliders. We analyse hydrodynamical simulations… ▽ More Terrestrial particle accelerators collide charged particles, then watch the trajectory of outgoing debris - but they cannot manipulate dark matter. Fortunately, dark matter is the main component of galaxy clusters, which are continuously pulled together by gravity. We show that galaxy cluster mergers can be exploited as enormous, natural dark matter colliders. We analyse hydrodynamical simulations of a universe containing self-interacting dark matter (SIDM) in which all particles interact via gravity, and dark matter particles can also scatter off each other via a massive mediator. During cluster collisions, SIDM spreads out and lags behind cluster member galaxies. Individual systems can have quirky dynamics that makes them difficult to interpret. Statistically, however, we find that the mean or median of dark matter's spatial offset in many collisions can be robustly modelled, and is independent of our viewing angle and halo mass even in collisions between unequal-mass systems. If the SIDM cross-section were sigma/m = 0.1cm^2/g = 0.18 barn/GeV, the 'bulleticity' lag would be ~5 percent that of gas due to ram pressure, and could be detected at 95 percent confidence in weak lensing observations of ~100 well-chosen clusters. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 8 pages, 5 figures plus appendices. MNRAS in press

Weak gravitational lensing measurements of Abell 2744 using JWST and shear measurement algorithm pyRRG-JWST

Authors: David R. Harvey, Richard Massey

Abstract: We update the publicly available weak lensing shear measurement algorithm pyRRG for the James Webb Space Telescope, and apply it to UNCOVER DR1 imaging of galaxy cluster Abell 2744. At short wavelengths (< 2.5$μ$m), shear measurements are consistent between independent observations through different JWST bandpasses, and calibrated within 1.5% of those from the Hubble Space Telescope. At longer wav… ▽ More We update the publicly available weak lensing shear measurement algorithm pyRRG for the James Webb Space Telescope, and apply it to UNCOVER DR1 imaging of galaxy cluster Abell 2744. At short wavelengths (< 2.5$μ$m), shear measurements are consistent between independent observations through different JWST bandpasses, and calibrated within 1.5% of those from the Hubble Space Telescope. At longer wavelengths, shear is underestimated by ~5%, probably due to coarser pixellisation. We model the spatially varying Point Spread Function (PSF) using WebbPSF, whose moments are within 0.05 of real stars near the centre of the mosaic, where there are sufficient stars to also generate an empirical model. We measure shear from up to 162 galaxies arcminute$^2$ to derive a map of (dark plus baryonic) mass with 12 arcsecond (55 kpc) spatial resolution. All code, catalogues and maps are available from https://github.com/davidharvey1986/pyRRG. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Accepted MNRAS

Data downloaded via parachute from a NASA super-pressure balloon

Authors: Ellen L. Sirks, Richard Massey, Ajay S. Gill, Jason Anderson, Steven J. Benton, Anthony M. Brown, Paul Clark, Joshua English, Spencer W. Everett, Aurelien A. Fraisse, Hugo Franco, John W. Hartley, David Harvey, Bradley Holder, Andrew Hunter, Eric M. Huff, Andrew Hynous, Mathilde Jauzac, William C. Jones, Nikky Joyce, Duncan Kennedy, David Lagattuta, Jason S. -Y. Leung, Lun Li, Stephen Lishman , et al. (18 additional authors not shown)

Abstract: In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) cap… ▽ More In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) capsules containing 5 TB solid state data storage, plus a GNSS receiver, Iridium transmitter, and parachute. Data from the telescope were copied to these, and two were dropped over Argentina. They drifted 61 km horizontally while they descended 32 km, but we predicted their descent vectors within 2.4 km: in this location, the discrepancy appears irreducible below 2 km because of high speed, gusty winds and local topography. The capsules then reported their own locations to within a few metres. We recovered the capsules and successfully retrieved all of superBIT's data - despite the telescope itself being later destroyed on landing. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 12 pages

Journal ref: Aerospace 2023, 10, 960

Beyond the Ultra-deep Frontier Fields And Legacy Observations (BUFFALO): a high-resolution strong + weak-lensing view of Abell 370

Authors: Anna Niemiec, Mathilde Jauzac, Dominique Eckert, David Lagattuta, Keren Sharon, Anton M. Koekemoer, Keiichi Umetsu, Ana Acebron, Jose M. Diego, David Harvey, Eric Jullo, Vasily Kokorev, Marceau Limousin, Guillaume Mahler, Priyamvada Natarajan, Mario Nonino, Juan D. Remolina, Charles Steinhardt, Sut-Ieng Tam, Adi Zitrin

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

Submitted 7 July, 2023; originally announced July 2023.

Comments: 29 pages, 17 figures, 3 tables

Strong gravitational lensing's 'external shear' is not shear

Authors: Amy Etherington, James W. Nightingale, Richard Massey, Sut-Ieng Tam, XiaoYue Cao, Anna Niemiec, Qiuhan He, Andrew Robertson, Ran Li, Aristeidis Amvrosiadis, Shaun Cole, Jose M. Diego, Carlos S. Frenk, Brenda L. Frye, David Harvey, Mathilde Jauzac, Anton M. Koekemoer, David J. Lagattuta, Marceau Limousin, Guillaume Mahler, Ellen Sirks, Charles L. Steinhardt

Abstract: The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus 'external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Spac… ▽ More The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus 'external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Space Telescope images, or in 50 mock images of lenses with complex distributions of mass. Instead, the best-fit shear is aligned with the major or minor axis of 88% of lens galaxies; and the amplitude of the external shear increases if that galaxy is disky. We conclude that 'external shear' attached to a power law model is not physically meaningful, but a fudge to compensate for lack of model complexity. Since it biases other model parameters that are interpreted as physically meaningful in several science analyses (e.g. measuring galaxy evolution, dark matter physics or cosmological parameters), we recommend that future studies of galaxy-scale strong lensing should employ more flexible mass models. △ Less

Submitted 17 November, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: 13 pages, 11 figures, submitted to MNRAS and Review Response Submitted

AEROS: Oceanographic Hyperspectral Imaging and Argos-Tracking CubeSat

Authors: Sophie Prendergast, Cadence Payne, Miles Lifson, Christian Haughwout, Marcos Tieppo, Paulo Figueiredo, André Guerra, Alexander Costa, Helder Magalhães, Tiago Hormigo, Francisco Câmara, Carlos Mano, Pedro Pinheiro, Alvin D. Harvey, Bruno Macena, Luis F. Azevedo, Miguel Martin, Tiago Miranda, Eduardo Pereira, João Faria, Inês Castelão, Catarina Cecilio, Emanuel Castanho, Kerri Cahoy, Manuel Coutinho , et al. (2 additional authors not shown)

Abstract: AEROS is a 3U CubeSat pathfinder toward a future ocean-observing constellation, targeting the Portuguese Atlantic region. AEROS features a miniaturized, high-resolution Hyperspectral Imager (HSI), a 5MP RGB camera, and a Software Defined Radio (SDR). The sensor generated data will be processed and aggregated for end-users in a new web-based Data Analysis Center (DAC). The HSI has 150 spectrally co… ▽ More AEROS is a 3U CubeSat pathfinder toward a future ocean-observing constellation, targeting the Portuguese Atlantic region. AEROS features a miniaturized, high-resolution Hyperspectral Imager (HSI), a 5MP RGB camera, and a Software Defined Radio (SDR). The sensor generated data will be processed and aggregated for end-users in a new web-based Data Analysis Center (DAC). The HSI has 150 spectrally contiguous bands covering visible to near-infrared with 10 nm bandwidth. The HSI collects ocean color data to support studies of oceanographic characteristics known to influence the spatio-temporal distribution and movement behavior of marine organisms. Usage of an SDR expands AEROS's operational and communication range and allows for remote reconfiguration. The SDR receives, demodulates, and retransmits short duration messages, from sources including tagged marine organisms, autonomous vehicles, subsurface floats, and buoys. The future DAC will collect, store, process, and analyze acquired data, taking advantage of its ability to disseminate data across the stakeholders and the scientific network. Correlation of animal-borne Argos platform locations and oceanographic data will advance fisheries management, ecosystem-based management, monitoring of marine protected areas, and bio-oceanographic research in the face of a rapidly changing environment. For example, correlation of oceanographic data collected by the HSI, geolocated with supplementary images from the RGB camera and fish locations, will provide researchers with near real-time estimates of essential oceanographic variables within areas selected by species of interest. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 13 pages, 16 figures, Manuscript presented at the 73rd International Astronautical Congress, IAC 2022, Paris, France, 18 - 22 September 2022

Constraints on dark matter self-interaction from the internal density profiles of X-COP galaxy clusters

Authors: D. Eckert, S. Ettori, A. Robertson, R. Massey, E. Pointecouteau, D. Harvey, I. G. McCarthy

Abstract: TThe fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally-bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter $α\approx0.18$, in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if dark matt… ▽ More TThe fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally-bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter $α\approx0.18$, in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if dark matter particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measure the structural properties of the 12 massive X-COP galaxy clusters from observations of their hot gaseous atmosphere using the X-ray observatory XMM-Newton, and of the Sunyaev-Zeldovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measure Einasto shape parameters with mean $\langleα\rangle = 0.19 \pm 0.03$ and intrinsic scatter $σ_α=0.06$, in close agreement with the prediction of the cold dark matter paradigm. We use cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We use the fitted relation to turn our measurements of $α$ into constraints on the self-interaction cross section, which imply $σ/m < 0.19$ cm$^2$/g (95% confidence level) at collision velocity $v_\mathrm{DM-DM}\sim1,000$ km/s. This is lower than the interaction cross-section required for dark matter self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: Re-submitted to A&A after revision

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

PyAutoLens: Open-Source Strong Gravitational Lensing

Authors: James. W. Nightingale, Richard G. Hayes, Ashley Kelly, Aristeidis Amvrosiadis, Amy Etherington, Qiuhan He, Nan Li, XiaoYue Cao, Jonathan Frawley, Shaun Cole, Andrea Enia, Carlos S. Frenk, David R. Harvey, Ran Li, Richard J. Massey, Mattia Negrello, Andrew Robertson

Abstract: Strong gravitational lensing, which can make a background source galaxy appears multiple times due to its light rays being deflected by the mass of one or more foreground lens galaxies, provides astronomers with a powerful tool to study dark matter, cosmology and the most distant Universe. PyAutoLens is an open-source Python 3.6+ package for strong gravitational lensing, with core features includi… ▽ More Strong gravitational lensing, which can make a background source galaxy appears multiple times due to its light rays being deflected by the mass of one or more foreground lens galaxies, provides astronomers with a powerful tool to study dark matter, cosmology and the most distant Universe. PyAutoLens is an open-source Python 3.6+ package for strong gravitational lensing, with core features including fully automated strong lens modeling of galaxies and galaxy clusters, support for direct imaging and interferometer datasets and comprehensive tools for simulating samples of strong lenses. The API allows users to perform ray-tracing by using analytic light and mass profiles to build strong lens systems. Accompanying PyAutoLens is the autolens workspace (see https://github.com/Jammy2211/autolens_workspace), which includes example scripts, lens datasets and the HowToLens lectures in Jupyter notebook format which introduce non experts to strong lensing using PyAutoLens. Readers can try PyAutoLens right now by going to the introduction Jupyter notebook on Binder (see https://mybinder.org/v2/gh/Jammy2211/autolens_workspace/master) or checkout the readthedocs (see https://pyautolens.readthedocs.io/en/latest/) for a complete overview of PyAutoLens's features. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 2 pages, 1 figure

Journal ref: Journal of Open Source Software, 6(58), 2525 (2021)

The impact of Self-Interacting Dark Matter on the Intrinsic Alignments of Galaxies

Authors: David Harvey, Nora Elisa Chisari, Andrew Robertson, Ian G. McCarthy

Abstract: The formation and evolution of galaxies is known to be sensitive to tidal processes leading to intrinsic correlations between their shapes and orientations. Such correlations can be measured to high significance today, suggesting that cosmological information can be extracted from them. Among the most pressing questions in particle physics and cosmology is the nature of dark matter. If dark matter… ▽ More The formation and evolution of galaxies is known to be sensitive to tidal processes leading to intrinsic correlations between their shapes and orientations. Such correlations can be measured to high significance today, suggesting that cosmological information can be extracted from them. Among the most pressing questions in particle physics and cosmology is the nature of dark matter. If dark matter is self-interacting, it can leave an imprint on galaxy shapes. In this work, we investigate whether self-interactions can produce a long-lasting imprint on intrinsic galaxy shape correlations. We investigate this observable at low redshift ($z<0.4$) using a state-of-the-art suite of cosmological hydro-dynamical simulations where the dark matter model is varied. We find that dark matter self-interactions induce a mass dependent suppression in the intrinsic alignment signal by up to 50\% out to ten's of mega-parsecs, showing that self-interactions can impact structure outside the very core of clusters. We find evidence that self-interactions have a scale-dependent impact on the intrinsic alignment signal that is sufficiently different from signatures introduced by differing baryonic physics prescriptions, suggesting that it is detectable with up-coming all-sky surveys. △ Less

Submitted 25 June, 2021; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: Accepted MNRAS

A 4% measurement of $H_0$ using the cumulative distribution of strong-lensing time delays in doubly-imaged quasars

Authors: David Harvey

Abstract: In the advent of large scale surveys, individually modelling strong-gravitational lenses and their counterpart time-delays in order to precisely measure $H_0$ will become computationally expensive, and highly complex. A complimentary approach is to study the cumulative distribution function (CDF) of time-delays where the global population of lenses is modelled along with $H_0$. In this paper we us… ▽ More In the advent of large scale surveys, individually modelling strong-gravitational lenses and their counterpart time-delays in order to precisely measure $H_0$ will become computationally expensive, and highly complex. A complimentary approach is to study the cumulative distribution function (CDF) of time-delays where the global population of lenses is modelled along with $H_0$. In this paper we use a suite of hydro-dynamical simulations to estimate the CDF of time-delays from doubly-imaged quasars for a realistic distribution of lenses. We find that the CDFs exhibit large amounts of halo-halo variance, regulated by the density profile inner slope and the total mass within $5$kpc. With the objective of fitting to data, we compress the CDFs using Principal Component Analysis and fit a Gaussian Processes Regressor consisting of three physical features: the redshift of the lens, $z_{\rm L}$; the power law index of the halo, $α$, and the mass within $5$kpc, plus four cosmological features. Assuming a flat Universe, we fit our model to 27 doubly-imaged quasars finding $H_0=71^{+2}_{-3}$km/s/Mpc, $z_{\rm L} = 0.36_{-0.09}^{+0.2} $, $α=-1.8_{-0.1}^{+0.1}$, $\log(M(<5$kpc$)/M_\odot)=11.1_{-0.1}^{+0.1} $, $Ω_{\rm M} = 0.3_{-0.04}^{+0.04} $ and $Ω_{\rm Λ}=0.7_{-0.04}^{+0.04}$. We compare our estimates of $z_{\rm L}$ and $\log(M(<5$kpc$)/M_\odot)$ to the data and find that within the sensitivity of the data, they are not systematically biased. We generate mock CDFs and find with that the Vera Rubin Observatory (VRO) could measure $σ/H_0$ to $<3\%$, limited by the precision of the model. If we are to exploit fully VRO, we require simulations that sample a larger proportion of the lens population, with a variety of feedback models, exploring all possible systematics. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: https://www.youtube.com/watch?v=8PFgmXI8g0o ---------------------------- https://www.universiteitleiden.nl/en/news/2020/11/universe-expansion-speed-measured-in-bending-starlight

Reconciling galaxy cluster shapes, measured by theorists vs observers

Authors: David Harvey, Andrew Robertson, Sut-Ieng Tam, Mathilde Jauzac, Richard Massey, Jason Rhodes, Ian G. McCarthy

Abstract: If properly calibrated, the shapes of galaxy clusters can be used to investigate many physical processes: from feedback and quenching of star formation, to the nature of dark matter. Theorists frequently measure shapes using moments of inertia of simulated particles'. We instead create mock (optical, X-ray, strong- and weak-lensing) observations of the twenty-two most massive (… ▽ More If properly calibrated, the shapes of galaxy clusters can be used to investigate many physical processes: from feedback and quenching of star formation, to the nature of dark matter. Theorists frequently measure shapes using moments of inertia of simulated particles'. We instead create mock (optical, X-ray, strong- and weak-lensing) observations of the twenty-two most massive ($\sim10^{14.7}\,M_\odot$) relaxed clusters in the BAHAMAS simulations. We find that observable measures of shape are rounder. Even when moments of inertia are projected into 2D and evaluated at matched radius, they overestimate ellipticity by 56\% (compared to observable strong lensing) and 430\% (compared to observable weak lensing). Therefore, we propose matchable quantities and test them using observations of eight relaxed clusters from the {\emph Hubble Space Telescope} and {\emph Chandra X-Ray Observatory}. We also release our HST data reduction and lensing analysis software to the community. In real clusters, the ellipticity and orientation angle at all radii are strongly correlated. In simulated clusters, the ellipticity of inner ($<r_{\mathrm{vir}}/20$) regions becomes decoupled: for example with greater misalignment of the central cluster galaxy. This may indicate overly efficient implementation of feedback from active galactic nuclei. Future exploitation of cluster shapes as a function of radii will require better understanding of core baryonic processes. Exploitation of shapes on any scale will require calibration on simulations extended all the way to mock observations. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: Accepted to MNRAS

The distribution of dark matter and gas spanning six megaparsecs around the post-merger galaxy cluster MS0451-03

Authors: Sut-Ieng Tam, Mathilde Jauzac, Richard Massey, David Harvey, Dominique Eckert, Harald Ebeling, Richard S. Ellis, Vittorio Ghirardini, Baptiste Klein, Jean-Paul Kneib, David Lagattuta, Priyamvada Natarajan, Andrew Robertson, Graham P. Smith

Abstract: Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures.… ▽ More Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions. △ Less

Submitted 28 July, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 21 pages, 13 figures, updated to match MNRAS version

From dwarf galaxies to galaxy clusters: Self-Interacting Dark Matter over 7 orders of magnitude in halo mass

Authors: Kyrylo Bondarenko, Anastasia Sokolenko, Alexey Boyarsky, Andrew Robertson, David Harvey, Yves Revaz

Abstract: In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from ob… ▽ More In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by CDM simulations that include baryons, while for SIDM with a velocity-independent cross-section of 1 cm$^2/$g, the simulated galaxy clusters have mean maximal surface densities that are below those of observed systems by an amount greater than the standard deviation of the observed maximal surface density at fixed mass. For less massive systems both CDM and SIDM agree with the observation equally well. △ Less

Submitted 22 January, 2021; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: Journal version; minor changes

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

Authors: O. Tihhonova, F. Courbin, D. Harvey, S. Hilbert, A. Peel, C. E. Rusu, C. D. Fassnacht, V. Bonvin, P. J. Marshall, G. Meylan, D. Sluse, S. H. Suyu, T. Treu, K. C. Wong

Abstract: We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse… ▽ More We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $κ_{\rm ext}$ at the position of B1608+656 is $κ_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted MNRAS

The BUFFALO HST Survey

Authors: Charles L. Steinhardt, Mathilde Jauzac, Ana Acebron, Hakim Atek, Peter Capak, Iary Davidzon, Dominique Eckert, David Harvey, Anton M. Koekemoer, Claudia D. P. Lagos, Guillaume Mahler, Mireia Montes, Anna Niemiec, Mario Nonino, P. A. Oesch, Johan Richard, Steven A. Rodney, Matthieu Schaller, Keren Sharon, Louis-Gregory Strolger, Joseph Allingham, Adam Amara, Yannick Bah'e, Celine Boehm, Sownak Bose , et al. (70 additional authors not shown)

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

Submitted 13 February, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: Accepted ApJS; MAST archive will be live concurrent with publication

Setting the scene for BUFFALO: A study of the matter distribution in the HFF galaxy cluster MACS J0416.1-2403 and its parallel field

Authors: Elizabeth J. Gonzalez, Martín Chalela, Mathilde Jauzac, Dominique Eckert, Matthieu Schaller, David Harvey, Anna Niemiec, Anton M. Koekemoer, David Barnes, Doug Clowe, Thomas Connor, José M. Diego, Juan D. Remolina Gonzalez, Charles L. Steinhardt

Abstract: In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter dis… ▽ More In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5σ$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 15 pages, 8 figures, 3 tables, submitted to MNRAS

Exploiting flux ratio anomalies to probe warm dark matter in future large scale surveys

Authors: David Harvey, Wessel Valkenburg, Amelie Tamone, Alexey Boyarsky, Frederic Courbin, Mark Lovell

Abstract: Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic ba… ▽ More Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: Accepted MNRAS

The shape of relaxed galaxy clusters and the public release of a HST shape measurement code, pyRRG

Authors: David Harvey, Sut-Ieng Tam, Mathilde Jauzac, Richard Massey, Jason Rhodes

Abstract: We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupl… ▽ More We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupled with the BCG shape correlated with both the shape predicted by strong gravitational lensing and the X-ray isophote. Conversely we find no such correlation between the shape as predicted by the weak lensing and the other three probes suggesting any coupling between the inner and outer regions is weak. We also present in this paper the public release of the HST weak lensing shape measurement code pyRRG, directly available from PyPi (https://pypi.org/project/pyRRG/). This python3.7 code, based on Rhodes et al. (2000) adopts an automated star-galaxy classifier based on a Random Forest and outputs scientifically useful products such as weak lensing catalogues suitable for the mass mapping algorithm Lenstool. △ Less

Submitted 20 November, 2019; v1 submitted 14 November, 2019; originally announced November 2019.

Galaxy structure with strong gravitational lensing: decomposing the internal mass distribution of massive elliptical galaxies

Authors: James. W. Nightingale, Richard J. Massey, David R. Harvey, Andrew P. Cooper, Amy Etherington, Sut-Ieng Tam, Richard G. Hayes

Abstract: We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central b… ▽ More We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central bulge surrounded by an extended envelope of stellar material. Whilst these two components look photometrically similar, their distinct lensing effects permit a clean decomposition of their mass structure. This allows us to infer two key pieces of information about each lens galaxy: (i) the stellar mass distribution (without invoking stellar populations models) and (ii) the inner dark matter halo mass. We argue that these two measurements are crucial to testing models of ME formation, as the stellar mass profile provides a diagnostic of baryonic accretion and feedback whilst the dark matter mass places each galaxy in the context of LCDM large scale structure formation. We also detect large rotational offsets between the two stellar components and a lopsidedness in their outer mass distributions, which hold further information on the evolution of each ME. Finally, we discuss how this approach can be extended to galaxies of all Hubble types and what implication our results have for studies of strong gravitational lensing. △ Less

Submitted 12 July, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: 19 pages, 4 figures. source code: https://github.com/Jammy2211/PyAutoLens

Observable tests of self-interacting dark matter in galaxy clusters: BCG wobbles in a constant density core

Authors: David Harvey, Andrew Robertson, Richard Massey, Ian G. McCarthy

Abstract: Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed bee… ▽ More Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed been observed - but only interpreted via a simplified, analytic model of oscillations. We compare these observations to the BAHAMAS-SIDM suite of cosmological simulations, which include SIDM and a fully hydrodynamical treatment of star formation and feedback. We predict that the median offset of BCGs increases with the SIDM cross-section, cluster mass and the amount of stellar mass within 10kpc, while CDM exhibits no trend in mass. Interpolating between the simulated cross-sections, we find that the observations (of 10 clusters) is consistent with CDM at the ~1.5$σ$ level, and prefer cross-section $σ$/m < 0.12(0.39)cm$^2$/g at 68% (95%) confidence level. This is on the verge of ruling out velocity-independent dark matter self-interactions as the solution to discrepancies between the predicted and observed behaviour of dwarf galaxies, and will be improved by larger surveys by Euclid or SuperBIT. △ Less

Submitted 13 September, 2019; v1 submitted 17 December, 2018; originally announced December 2018.

The core of the massive cluster merger MACS J0417.5-1154 as seen by VLT/MUSE

Authors: Mathilde Jauzac, Guillaume Mahler, Alastair C. Edge, Keren Sharon, Steven Gillman, Harald Ebeling, David Harvey, Johan Richard, Steven L. Hamer, Michele Fumagalli, A. Mark Swinbank, Jean-Paul Kneib, Richard Massey, Philippe Salome

Abstract: We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partia… ▽ More We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partial and radial image close to the Brightest Cluster Galaxy (BCG) only discernible thanks to its strong [OII] line emission. The best-fit mass model (rms of 0.38\arcsec) yields a two-dimensional enclosed mass of $M({\rm R < 200\,kpc}) = (1.77\pm0.03)\times10^{14}\,\msun$ and almost perfect alignment between the peaks of the BCG light and the dark matter of ($0.5\pm0.5$)\arcsec . Our finding that a significant misalignment results when the radial image of \emph{The Doughnut} is omitted serves as an important caveat for studies of BCG-dark matter offsets in galaxy clusters. Using \emph{Chandra} data to map the intra-cluster gas, we observe an offset between the gas and dark-matter peaks of ($1.7\pm0.5$)\arcsec, and excellent alignment of the X-ray peak with the location of optical emission line associated with the BCG. We interpret all observational evidence in the framework of on-going merger activity, noting specifically that the coincidence between the gas peak and the peak of blue light from the BCG may be evidence of dense, cold gas leading to direct star formation. We use the surface area $σ_μ$ above a given magnification factor $μ$ as a metric to estimate the lensing power of MACS\,J0417. We obtain $σ(μ> 3) = 0.22$\,arcmin$^2$ confirming MACS\,J0417 as an efficient gravitational lens. Finally, we discuss the differences between our mass model and Mahler et al. (2018). △ Less

Submitted 6 December, 2018; v1 submitted 31 October, 2018; originally announced November 2018.

Comments: 15 pages, 6 figures, 8 tables, accepted for publication in MNRAS

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

Authors: Andrew Robertson, David Harvey, Richard Massey, Vincent Eke, Ian G. McCarthy, Mathilde Jauzac, Baojiu Li, Joop Schaye

Abstract: We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations.… ▽ More We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing. △ Less

Submitted 1 June, 2020; v1 submitted 12 October, 2018; originally announced October 2018.

Comments: 19 pages, 13 figures, updated to match MNRAS version

The impact of cored density profiles on the observable quantities of dwarf spheroidal galaxies

Authors: David Harvey, Yves Revaz, Andrew Robertson, Loic Hausammann

Abstract: We modify the chemo-dynamical code GEAR to simulate the impact of self-interacting dark matter on the observable quantities of 19 low mass dwarf galaxies with a variety star forming properties. We employ a relatively high, velocity independent cross-section of $σ/m = 10$cm$^2$/g and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-s… ▽ More We modify the chemo-dynamical code GEAR to simulate the impact of self-interacting dark matter on the observable quantities of 19 low mass dwarf galaxies with a variety star forming properties. We employ a relatively high, velocity independent cross-section of $σ/m = 10$cm$^2$/g and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-sight velocities, the chemical abundance and the star formation history. We find that despite the creation of large cores at the centre of the dark matter haloes, the impact of SIDM on the observable quantities of quenched galaxies is indiscernible, dominated mostly by the stochastic build up of the stellar matter. As such we conclude that it is impossible to make global statements on the density profile of dwarf galaxies from single or small samples. Although based mostly on quenched galaxies, this finding supports other recent work putting into question the reliability of inferred cored density profiles that are derived from observed line-of-sight velocities. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: MNRAS Accepted 2018 August 29. Received 2018 August 13; in original form 2018 May 29

The shape of galaxy dark matter halos in massive galaxy clusters: Insights from strong gravitational lensing

Authors: M. Jauzac, D. Harvey, R. Massey

Abstract: We assess how much unused strong lensing information is available in the deep \emph{Hubble Space Telescope} imaging and VLT/MUSE spectroscopy of the \emph{Frontier Field} clusters. As a pilot study, we analyse galaxy cluster MACS\,J0416.1-2403 ($z$$=$$0.397$, $M(R<200\,{\rm kpc})$$=$$1.6$$\times$$10^{14}\msun$), which has 141 multiple images with spectroscopic redshifts. We find that many addition… ▽ More We assess how much unused strong lensing information is available in the deep \emph{Hubble Space Telescope} imaging and VLT/MUSE spectroscopy of the \emph{Frontier Field} clusters. As a pilot study, we analyse galaxy cluster MACS\,J0416.1-2403 ($z$$=$$0.397$, $M(R<200\,{\rm kpc})$$=$$1.6$$\times$$10^{14}\msun$), which has 141 multiple images with spectroscopic redshifts. We find that many additional parameters in a cluster mass model can be constrained, and that adding even small amounts of extra freedom to a model can dramatically improve its figures of merit. We use this information to constrain the distribution of dark matter around cluster member galaxies, simultaneously with the cluster's large-scale mass distribution. We find tentative evidence that some galaxies' dark matter has surprisingly similar ellipticity to their stars (unlike in the field, where it is more spherical), but that its orientation is often misaligned. When non-coincident dark matter and baryonic halos are allowed, the model improves by 35\%. This technique may provide a new way to investigate the processes and timescales on which dark matter is stripped from galaxies as they fall into a massive cluster. Our preliminary conclusions will be made more robust by analysing the remaining five \emph{Frontier Field} clusters. △ Less

Submitted 3 April, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: 6 pages, 3 figures, 2 tables - Accepted for publication in MNRAS

H0LiCOW VIII. A weak lensing measurement of the external convergence in the field of the lensed quasar HE$\,$0435$-$1223

Authors: O. Tihhonova, F. Courbin, D. Harvey, S. Hilbert, C. E. Rusu, C. D. Fassnacht, V. Bonvin, P. J. Marshall, G. Meylan, D. Sluse, S. H. Suyu, T. Treu, K. C. Wong

Abstract: We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$σ$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and… ▽ More We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$σ$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $κ=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object. △ Less

Submitted 23 November, 2017; originally announced November 2017.

Comments: 14 pages, 8 figures, 2 tables

Growing a `Cosmic Beast': Observations and Simulations of MACS J0717.5+3745

Authors: M. Jauzac, D. Eckert, M. Schaller, J. Schwinn, R. Massey, Y. Bahé, C. Baugh, D. Barnes, C. Dalla Vecchia, H. Ebeling, D. Harvey, E. Jullo, S. T. Kay, J. -P. Kneib, M. Limousin, E. Medezinski, P. Natarajan, M. Nonino, A. Robertson, S. I. Tam, K. Umetsu

Abstract: We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected r… ▽ More We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters. △ Less

Submitted 27 August, 2018; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 18 pages, 6 figures, 8 tables, Accepted for publication in MNRAS

Dark matter dynamics in Abell 3827: new data consistent with standard Cold Dark Matter

Authors: Richard Massey, David Harvey, Jori Liesenborgs, Johan Richard, Stuart Stach, Mark Swinbank, Peter Taylor, Liliya Williams, Douglas Clowe, Frederic Courbin, Alastair Edge, Holger Israel, Mathilde Jauzac, Remy Joseph, Eric Jullo, Thomas D. Kitching, Adrienne Leonard, Julian Merten, Daisuke Nagai, James Nightingale, Andrew Robertson, Luis Javier Romualdez, Prasenjit Saha, Renske Smit, Sut Ieng Tam , et al. (1 additional authors not shown)

Abstract: We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The n… ▽ More We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: 10 pages, 5 figures. Comments welcome

A Potential Recoiling Supermassive Black Hole CXO J101527.2+625911

Authors: D. -C. Kim, Ilsang Yoon, G. C. Privon, A. S. Evans, D. Harvey, S. Stierwalt, Ji Hoon Kim

Abstract: We have carried out a systematic search for recoiling supermassive black holes (rSMBH) using the Chandra Source and SDSS Cross Matched Catalog. From the survey, we have detected a potential rSMBH, 'CXO J101527.2+625911' at z=0.3504. The CXO J101527.2+625911 has a spatially offset (1.26$\pm$0.05 kpc) active SMBH and kinematically offset broad emission lines (175$\pm$25 km s$^{\rm -1}$ relative to s… ▽ More We have carried out a systematic search for recoiling supermassive black holes (rSMBH) using the Chandra Source and SDSS Cross Matched Catalog. From the survey, we have detected a potential rSMBH, 'CXO J101527.2+625911' at z=0.3504. The CXO J101527.2+625911 has a spatially offset (1.26$\pm$0.05 kpc) active SMBH and kinematically offset broad emission lines (175$\pm$25 km s$^{\rm -1}$ relative to systemic velocity). The observed spatial and velocity offsets suggest this galaxy could be a rSMBH, but we also have considered a possibility of dual SMBH scenario. The column density towards the galaxy center was found to be Compton thin, but no X-ray source was detected. The non-detection of the X-ray source in the nucleus suggests either there is no obscured actively accreting SMBH, or there exists an SMBH but has a low accretion rate (i.e. low-luminosity AGN (LLAGN)). The possibility of the LLAGN was investigated and found to be unlikely based on the H$α$ luminosity, radio power, and kinematic arguments. This, along with the null detection of X-ray source in the nucleus supports our hypothesis that the CXO J101527.2+625911 is a rSMBH. Our GALFIT analysis shows the host galaxy to be a bulge-dominated elliptical. The weak morphological disturbance and small spatial and velocity offsets suggest that CXO J101527.2+625911 could be in the final stage of merging process and about to turn into a normal elliptical galaxy. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: ApJ accepted

A detection of wobbling Brightest Cluster Galaxies within massive galaxy clusters

Authors: David Harvey, F. Courbin, J. P. Kneib, Ian G. McCarthy

Abstract: A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a Brightest Cluster Galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomenon is absent with standard cold dark matter where a cuspy den… ▽ More A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a Brightest Cluster Galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomenon is absent with standard cold dark matter where a cuspy density profile keeps a BCG tightly bound at the centre. We test this hypothesis using cosmological simulations and deep observations of 10 galaxy clusters acting as strong gravitational lenses. Modelling the BCG wobble as a simple harmonic oscillator, we measure the wobble amplitude, $A_w$, in the BAHAMAS suite of cosmological hydrodynamical simulations, finding an upper limit for the CDM paradigm of $A_w < 2$kpc at the 95% confidence limit. We carry out the same test on the data finding a non-zero amplitude of $A_w=11.82^{+7.3}_{-3.0}$kpc, with the observations dis-favouring $A_w = 0$ at the 3$σ$ confidence level. This detection of BCG wobbling is evidence for a dark matter core at the heart of galaxy clusters. It also shows that strong lensing models of clusters cannot assume that the BCG is exactly coincident with the large scale halo. While our small sample of galaxy clusters already indicates a non-zero $A_w$, with larger surveys, e.g. Euclid, we will be able to not only to confirm the effect but also to use it to determine whether or not the wobbling finds its origin in new fundamental physics or astrophysical process. △ Less

Submitted 9 August, 2017; v1 submitted 21 March, 2017; originally announced March 2017.

Comments: Accepted MNRAS

A test for skewed distributions of dark matter and a possible detection in galaxy cluster Abell 3827

Authors: Peter Taylor, Richard Massey, Mathilde Jauzac, Frédéric Courbin, David Harvey, Rémy Joseph, Andrew Robertson

Abstract: Simulations of self-interacting dark matter (SIDM) predict that dark matter should lag behind galaxies during a collision. If the interaction is mediated by a high-mass force carrier, the distribution of dark matter can also develop asymmetric dark matter tails. To search for this asymmetry, we compute the gravitational lensing properties of a mass distribution with a free {\em skewness} parameter… ▽ More Simulations of self-interacting dark matter (SIDM) predict that dark matter should lag behind galaxies during a collision. If the interaction is mediated by a high-mass force carrier, the distribution of dark matter can also develop asymmetric dark matter tails. To search for this asymmetry, we compute the gravitational lensing properties of a mass distribution with a free {\em skewness} parameter. We apply this to the dark matter around the four central galaxies in cluster Abell~3827. In the galaxy whose dark matter peak has previously been found to be offset, we tentatively measure a skewness $s=0.23^{+0.05}_{-0.22}$ in the same direction as the peak offset. Our method may be useful in future gravitational lensing analyses of colliding galaxy clusters and merging galaxies. △ Less

Submitted 13 April, 2017; v1 submitted 16 January, 2017; originally announced January 2017.

Abell 2744: Too much substructure for Lambda CDM?

Authors: J. Schwinn, M. Jauzac, C. M. Baugh, M. Bartelmann, D. Eckert, D. Harvey, P. Natarajan, R. Massey

Abstract: The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations i… ▽ More The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations is considered, extending the preliminary evaluation carried out by Jauzac et al. (2016). There are many tens of haloes in the simulation with a total mass comparable with or larger than that of Abell 2744. However, we find no haloes with a number and distribution of massive substructures (> 5 10^13 Msun) that is close to that inferred from the observations of Abell 2744. The application of extreme value statistics suggests that we would need a simulation of at least ten times the volume of the Millennium XXL to find a single dark matter halo with a similar internal structure to Abell 2744. Explaining the distribution of massive substructures in clusters is a new hurdle for hierarchical models to negotiate, which is not weakened by appeals to baryonic physics or uncertainty over the nature of the dark matter particle. △ Less

Submitted 2 February, 2017; v1 submitted 8 November, 2016; originally announced November 2016.

Comments: 12 pages, 5 figures, Accepted for publication in MNRAS, V2: minor changes in the manuscript

Looking for dark matter trails in colliding galaxy clusters

Authors: David Harvey, Andrew Robertson, Richard Massey, Jean-Paul Kneib

Abstract: If dark matter interacts, even weakly, via non-gravitational forces, simulations predict that it will be preferentially scattered towards the trailing edge of the halo during collisions between galaxy clusters. This will temporarily create a non-symmetric mass profile, with a trailing over-density along the direction of motion. To test this hypothesis, we fit (and subtract) symmetric halos to the… ▽ More If dark matter interacts, even weakly, via non-gravitational forces, simulations predict that it will be preferentially scattered towards the trailing edge of the halo during collisions between galaxy clusters. This will temporarily create a non-symmetric mass profile, with a trailing over-density along the direction of motion. To test this hypothesis, we fit (and subtract) symmetric halos to the weak gravitational data of 72 merging galaxy clusters observed with the Hubble Space Telescope. We convert the shear directly into excess κ and project in to a one dimensional profile. We generate numerical simulations and find that the one dimensional profile is well described with simple Gaussian approximations. We detect the weak lensing signal of trailing gas at a 4σ confidence, finding a mean gas fraction of Mgas/Mdm = 0.13 +/- 0.035. We find no evidence for scattered dark matter particles with a estimated scattering fraction of f = 0.03 +/- 0.05. Finally we find that if we can reduce the statistical error on the positional estimate of a single dark matter halo to <2.5", then we will be able to detect a scattering fraction of 10% at the 3σ level with current surveys. This poten- tially interesting new method can provide an important independent test for other complimentary studies of the self-interaction cross-section of dark matter. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: Accepted MNRAS

The Extraordinary Amount of Substructure in the Hubble Frontier Fields Cluster Abell 2744

Authors: M. Jauzac, D. Eckert, J. Schwinn, D. Harvey, C. M. Baugh, A. Robertson, S. Bose, R. Massey, M. Owers, H. Ebeling, H. Y. Shan, E. Jullo, J. -P. Kneib, J. Richard, H. Atek, B. Clément, E. Egami, H. Israel, K. Knowles, M. Limousin, P. Natarajan, M. Rexroth, P. Taylor, C. Tchernin

Abstract: We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to… ▽ More We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $Λ$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015). △ Less

Submitted 14 June, 2016; originally announced June 2016.

Comments: Submitted to MNRAS -- 19 pages, 8 figures, 6 tables

Systematic or Signal? How dark matter misalignments can bias strong lensing models of galaxy clusters

Authors: David Harvey, Jean-Paul Kneib, Mathilde Jauzac

Abstract: We explore how assuming that mass traces light in strong gravitational lensing models can lead to systematic errors in the predicted position of multiple images. Using a model based on the galaxy cluster MACSJ0416 (z = 0.397) from the Hubble Frontier Fields, we split each galactic halo into a baryonic and dark matter component. We then shift the dark matter halo such that it no longer aligns with… ▽ More We explore how assuming that mass traces light in strong gravitational lensing models can lead to systematic errors in the predicted position of multiple images. Using a model based on the galaxy cluster MACSJ0416 (z = 0.397) from the Hubble Frontier Fields, we split each galactic halo into a baryonic and dark matter component. We then shift the dark matter halo such that it no longer aligns with the baryonic halo and investigate how this affects the resulting position of multiple images. We find for physically motivated misalignments in dark halo position, ellipticity, position angle and density profile, that multiple images can move on average by more than 0.2" with individual images moving greater than 1". We finally estimate the full error induced by assuming that light traces mass and find that this assumption leads to an expected RMS error of 0.5", almost the entire error budget observed in the Frontier Fields. Given the large potential contribution from the assumption that light traces mass to the error budget in mass reconstructions, we predict that it should be possible to make a first significant detection and characterisation of dark halo misalignments in the Hubble Frontier Fields with strong lensing. Finally, we find that it may be possible to detect ~1kpc offsets between dark matter and baryons, the smoking gun for self-interacting dark matter, should the correct alignment of multiple images be observed. △ Less

Submitted 25 January, 2016; originally announced January 2016.

Comments: 7 Pages

A weak lensing comparability study of galaxy mergers that host AGNs

Authors: David Harvey, Frederic Courbin

Abstract: We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately lumi… ▽ More We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately luminous quasars at z~0.36 powered by disk instabilities, but with no intense star formation. On an initial test we found no evidence for a connection between the merger state of a galaxy and the profile of the halo, with the PSQ profile comparable to that of the other two samples and consistent with the Leauthaud et al. (2014) study of moderately luminous quasars in COSMOS. Given the compatibility of the two quasar samples, we combined these and found no evidence for any connection between black hole activity and the dark matter halo. All three mass profiles remained compatible with isothermality given the present data. △ Less

Submitted 14 May, 2015; originally announced May 2015.

Comments: 6 pages, 3 figures, 1 table, ACCEPTED MNRAS

The behaviour of dark matter associated with 4 bright cluster galaxies in the 10kpc core of Abell 3827

Authors: Richard Massey, Liliya Williams, Renske Smit, Mark Swinbank, Thomas Kitching, David Harvey, Mathilde Jauzac, Holger Israel, Douglas Clowe, Alastair Edge, Matt Hilton, Eric Jullo, Adrienne Leonard, Jori Liesenborgs, Julian Merten, Irshad Mohammed, Daisuke Nagai, Johan Richard, Andrew Robertson, Prasenjit Saha, Rebecca Santana, John Stott, Eric Tittley

Abstract: Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravit… ▽ More Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration. △ Less

Submitted 13 April, 2015; originally announced April 2015.

Comments: 15 pages, 9 figures

Journal ref: MNRAS 449 (2015) 3393

The non-gravitational interactions of dark matter in colliding galaxy clusters

Authors: David Harvey, Richard Massey, Thomas Kitching, Andy Taylor, Eric Tittley

Abstract: Collisions between galaxy clusters provide a test of the non-gravitational forces acting on dark matter. Dark matter's lack of deceleration in the `bullet cluster collision' constrained its self-interaction cross-section σ_DM/m < 1.25cm2/g (68% confidence limit) for long-ranged forces. Using the Chandra and Hubble Space Telescopes we have now observed 72 collisions, including both `major' and `min… ▽ More Collisions between galaxy clusters provide a test of the non-gravitational forces acting on dark matter. Dark matter's lack of deceleration in the `bullet cluster collision' constrained its self-interaction cross-section σ_DM/m < 1.25cm2/g (68% confidence limit) for long-ranged forces. Using the Chandra and Hubble Space Telescopes we have now observed 72 collisions, including both `major' and `minor' mergers. Combining these measurements statistically, we detect the existence of dark mass at 7.6σsignificance. The position of the dark mass has remained closely aligned within 5.8+/-8.2 kpc of associated stars: implying a self-interaction cross-section σ_DM/m < 0.47 cm2/g (95% CL) and disfavoring some proposed extensions to the standard model. △ Less

Submitted 13 April, 2015; v1 submitted 26 March, 2015; originally announced March 2015.

Comments: 5 Pages, 4 Figures and 18 pages supplementary information

Journal ref: Science, Vol. 347 no. 6229 pp. 1462-1465 (2015)

Hubble Frontier Fields: The Geometry and Dynamics of the Massive Galaxy Cluster Merger MACSJ0416.1-2403

Authors: Mathilde Jauzac, Eric Jullo, Dominique Eckert, Harald Ebeling, Johan Richard, Marceau Limousin, Hakim Atek, Jean-Paul Kneib, Benjamin Clément, Eiichi Egami, David Harvey, Kenda Knowles, Richard Massey, Priyamvada Natarajan, Benoît Neichel, Markus Rexroth

Abstract: We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dyn… ▽ More We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios. △ Less

Submitted 22 October, 2014; v1 submitted 11 June, 2014; originally announced June 2014.

Comments: Submitted to MNRAS - 17 pages, 13 figures, 3 tables

Observing Dark Worlds: A crowdsourcing experiment for dark matter mapping

Authors: David Harvey, Thomas D. Kitching, Joyce Noah-Vanhoucke, Ben Hamner, Tim Salimans

Abstract: We present the results and conclusions from the citizen science competition `Observing Dark Worlds', where we asked participants to calculate the positions of dark matter halos from 120 catalogues of simulated weak lensing galaxy data, using computational methods. In partnership with Kaggle (http://www.kaggle.com), 357 users participated in the competition which saw 2278 downloads of the data and… ▽ More We present the results and conclusions from the citizen science competition `Observing Dark Worlds', where we asked participants to calculate the positions of dark matter halos from 120 catalogues of simulated weak lensing galaxy data, using computational methods. In partnership with Kaggle (http://www.kaggle.com), 357 users participated in the competition which saw 2278 downloads of the data and 3358 submissions. We found that the best algorithms improved on the benchmark code, LENSTOOL by > 30% and could measure the positions of > 3x10^14MSun halos to less than 5'' and < 10^14MSun to within 1'. In this paper, we present a brief overview of the winning algorithms with links to available code. We also discuss the implications of the experiment for future citizen science competitions. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 25 Pages (Large spaced, equiv. 8 pages in MNRAS style), 7 Figures

On the cross-section of Dark Matter using substructure infall into galaxy clusters

Authors: David Harvey, Eric Tittley, Richard Massey, Thomas D. Kitching, Andy Taylor, Simon R. Pike, Scott T. Kay, Erwin T. Lau, Daisuke Nagai

Abstract: We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of… ▽ More We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1. △ Less

Submitted 20 November, 2013; v1 submitted 7 October, 2013; originally announced October 2013.

Comments: 14 pages, 11 figures

Dark Matter Astrometry: Accuracy of sub-halo positions for the measurement of self-interaction cross sections

Authors: David Harvey, Richard Massey, Thomas Kitching, Andy Taylor, Eric Jullo, Jean-Paul Kneib, Eric Tittley, Philip J. Marshall

Abstract: Direct evidence for the existence of dark matter and measurements of its interaction cross-section have been provided by the physical offset between dark matter and intra- cluster gas in merging systems like the Bullet Cluster. Although a smaller signal, this effect is more abundant in minor mergers where infalling substructure dark matter and gas are segregated. In such low-mass systems the gravi… ▽ More Direct evidence for the existence of dark matter and measurements of its interaction cross-section have been provided by the physical offset between dark matter and intra- cluster gas in merging systems like the Bullet Cluster. Although a smaller signal, this effect is more abundant in minor mergers where infalling substructure dark matter and gas are segregated. In such low-mass systems the gravitational lensing signal comes primarily from weak lensing. A fundamental step in determining such an offset in sub- structure is the ability to accurately measure the positions of dark matter sub-peaks. Using simulated Hubble Space Telescope observations, we make a first assessment of the precision and accuracy with which we can measure infalling groups using weak gravitational lensing. We demonstrate that using an existing and well-used mass re- construction algorithm can measure the positions of 1.5x1013 M substructures that have parent halos ten times more massive with a bias of less than 0.3 . In this regime, our analysis suggests the precision is sufficient to detect (at 3 σ statistical significance) the expected mean offset between dark matter and baryonic gas in infalling groups from a sample of 50 massive clusters. △ Less

Submitted 9 May, 2013; originally announced May 2013.

Comments: 12 pages, 10 figures

Journal ref: 2013MNRAS.433.1517H

BK Lyncis: The Oldest Old Nova?... And a Bellwether for Cataclysmic-Variable Evolution

Authors: Joseph Patterson, Helena Uthas, Jonathan Kemp, Enrique de Miguel, Thomas Krajci, Jerry Foote, Franz-Josef Hambsch, Tut Campbell, George Roberts, David Cejudo, Shawn Dvorak, Tonny Vanmunster, Robert Koff, David Skillman, David Harvey, Brian Martin, John Rock, David Boyd, Arto Oksanen, Etienne Morelle, Joseph Ulowetz, Anthony Kroes, Richard Sabo, Lasse Jensen

Abstract: We summarize the results of a 20-year campaign to study the light curves of BK Lyncis, a nova-like star strangely located below the 2-3 hour orbital period gap in the family of cataclysmic variables. Two apparent "superhumps" dominate the nightly light curves - with periods 4.6% longer, and 3.0% shorter, than P_orb. The first appears to be associated with the star's brighter states (V~14), while t… ▽ More We summarize the results of a 20-year campaign to study the light curves of BK Lyncis, a nova-like star strangely located below the 2-3 hour orbital period gap in the family of cataclysmic variables. Two apparent "superhumps" dominate the nightly light curves - with periods 4.6% longer, and 3.0% shorter, than P_orb. The first appears to be associated with the star's brighter states (V~14), while the second appears to be present throughout and becomes very dominant in the low state (V~15.7). Starting in the year 2005, the star's light curve became indistinguishable from that of a dwarf nova - in particular, that of the ER UMa subclass. Reviewing all the star's oddities, we speculate: (a) BK Lyn is the remnant of the probable nova on 30 December 101, and (b) it has been fading ever since, but has taken ~2000 years for the accretion rate to drop sufficiently to permit dwarf-nova eruptions. If such behavior is common, it can explain other puzzles of CV evolution. One: why the ER UMa class even exists (because all members can be remnants of recent novae). Two: why ER UMa stars and short-period novalikes are rare (because their lifetimes, which are essentially cooling times, are short). Three: why short-period novae all decline to luminosity states far above their true quiescence (because they're just getting started in their postnova cooling). Four: why the orbital periods, accretion rates, and white-dwarf temperatures of short-period CVs are somewhat too large to arise purely from the effects of gravitational radiation (because the unexpectedly long interval of enhanced postnova brightness boosts the mean mass-transfer rate). These are substantial rewards in return for one investment of hypothesis: that the second parameter in CV evolution, besides P_orb, is time since the last classical-nova eruption. △ Less

Submitted 23 December, 2012; originally announced December 2012.

Comments: PDF, 46 pages, 4 tables, 10 figures; in preparation; more info at http://cbastro.org/

Origins of weak lensing systematics, and requirements on future instrumentation (or knowledge of instrumentation)

Authors: Richard Massey, Henk Hoekstra, Thomas Kitching, Jason Rhodes, Mark Cropper, Jerome Amiaux, David Harvey, Yannick Mellier, Massimo Meneghetti, Lance Miller, Stephane Paulin-Henriksson, Sandrine Pires, Roberto Scaramella, Tim Schrabback

Abstract: The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicati… ▽ More The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission. △ Less

Submitted 29 October, 2012; originally announced October 2012.

Comments: 19 pages, 5 figures. Submitted to MNRAS. See companion (implementation) paper Cropper et al

Journal ref: MNRAS 429 (2013), 661

Superhumps in Cataclysmic Binaries. XXV. q_crit, epsilon(q), and Mass-Radius

Authors: Joseph Patterson, Jonathan Kemp, David Harvey, Robert Fried, Robert Rea, Berto Monard, Lewis Cook, David Skillman, Tonny Vanmunster, Greg Bolt, Eve Armstrong, Jennie McCormick, Thomas Krajci, Lasse Jensen, Jerry Gunn, Neil Butterworth, Jerry Foote, Marc Bos, Gianluca Masi, Paul Warhurst

Abstract: We report on successes and failures in searching for positive superhumps in cataclysmic variables, and show the superhumping fraction as a function of orbital period. Basically, all short-period systems do, all long-period systems don't, and a 50% success rate is found at P_orb=3.1+-0.2 hr. We can use this to measure the critical mass ratio for the creation of superhumps. With a mass-radius rela… ▽ More We report on successes and failures in searching for positive superhumps in cataclysmic variables, and show the superhumping fraction as a function of orbital period. Basically, all short-period systems do, all long-period systems don't, and a 50% success rate is found at P_orb=3.1+-0.2 hr. We can use this to measure the critical mass ratio for the creation of superhumps. With a mass-radius relation appropriate for cataclysmic variables, and an assumed mean white-dwarf mass of 0.75 M_sol, we find a mass ratio q_crit=0.35+-0.02. We also report superhump studies of several stars of independently known mass ratio: OU Virginis, XZ Eridani, UU Aquarii, and KV UMa (= XTE J1118+480). The latter two are of special interest, because they represent the most extreme mass ratios for which accurate superhump measurements have been made. We use these to improve the epsilon(q) calibration, by which we can infer the elusive q from the easy-to-measure epsilon (the fractional period excess of P_superhump over P_orb). This relation allows mass and radius estimates for the secondary star in any CV showing superhumps. The consequent mass-radius law shows an apparent discontinuity in radius near 0.2 M_sol, as predicted by the disrupted magnetic braking model for the 2.1-2.7 hour period gap. This is effectively the "empirical main sequence" for CV secondaries. △ Less

Submitted 15 July, 2005; originally announced July 2005.

Comments: PDF, 45 pages, 9 tables, 12 figures; accepted, in press, to appear November 2005, PASP; more info at http://cba.phys.columbia.edu/

Superhumps in Cataclysmic Binaries. XXIV. Twenty More Dwarf Novae

Authors: J. Patterson, J. Thorstensen, J. Kemp, D. Skillman, T. Vanmunster, D. Harvey, R. Fried, L. Jensen, L. Cook, R. Rea, B. Monard, J. McCormick, F. Velthuis, S. Walker, B. Martin, G. Bolt, E. Pavlenko, D. O'Donoghue, J. Gunn, R. Novak, G. Masi, G. Garradd, N. Butterworth, T. Krajci, J. Foote , et al. (1 additional authors not shown)

Abstract: We report precise measures of the orbital and superhump period in twenty more dwarf novae. For ten stars, we report new and confirmed spectroscopic periods - signifying the orbital period P_o - as well as the superhump period P_sh. These are GX Cas, HO Del, HS Vir, BC UMa, RZ Leo, KV Dra, KS UMa, TU Crt, QW Ser, and RZ Sge. For the remaining ten, we report a medley of P_o and P_sh measurements f… ▽ More We report precise measures of the orbital and superhump period in twenty more dwarf novae. For ten stars, we report new and confirmed spectroscopic periods - signifying the orbital period P_o - as well as the superhump period P_sh. These are GX Cas, HO Del, HS Vir, BC UMa, RZ Leo, KV Dra, KS UMa, TU Crt, QW Ser, and RZ Sge. For the remaining ten, we report a medley of P_o and P_sh measurements from photometry; most are new, with some confirmations of previous values. These are KV And, LL And, WX Cet, MM Hya, AO Oct, V2051 Oph, NY Ser, KK Tel, HV Vir, and RX J1155.4-5641. Periods, as usual, can be measured to high accuracy, and these are of special interest since they carry dynamical information about the binary. We still have not quite learned how to read the music, but a few things are clear. The fractional superhump excess epsilon [=(P_sh-P_o)/P_o] varies smoothly with P_o. The scatter of the points about that smooth curve is quite low, and can be used to limit the intrinsic scatter in M_1, the white dwarf mass, and the mass-radius relation of the secondary. The dispersion in M_1 does not exceed 24%, and the secondary-star radii scatter by no more than 11% from a fixed mass-radius relation. For the well-behaved part of epsilon(P_o) space, we estimate from superhump theory that the secondaries are 18+-6% larger than theoretical ZAMS stars. This affects some other testable predictions about the secondaries: at a fixed P_o, it suggests that the secondaries are (compared with ZAMS predictions) 40+-14% less massive, 12+-4% smaller, 19+-6% cooler, and less luminous by a factor 2.5(7). The presence of a well-defined mass-radius relation, reflected in a well-defined epsilon(P_o) relation, strongly limits effects of nuclear evolution in the secondaries. △ Less

Submitted 3 September, 2003; originally announced September 2003.

Comments: PDF, 62 pages, 7 tables, 21 figures; accepted, in press, to appear November 2003, PASP; more info at http://cba.phys.columbia.edu/

Envelope Structure of Starless Core L694-2 Derived from a Near-Infrared Extinction Map

Authors: Daniel W. A. Harvey, David J. Wilner, Charles J. Lada, Philip C. Myers, Joao Alves

Abstract: We present a near-infrared extinction study of the dark globule L694-2, a starless core that shows strong evidence for inward motions in molecular line profiles. The J,H, and K band data were taken using the European Southern Observatory New Technology Telescope. The best fit simple spherical power law model has index p=2.6 +/- 0.2, over the 0.036--0.1 pc range in radius sampled in extinction. T… ▽ More We present a near-infrared extinction study of the dark globule L694-2, a starless core that shows strong evidence for inward motions in molecular line profiles. The J,H, and K band data were taken using the European Southern Observatory New Technology Telescope. The best fit simple spherical power law model has index p=2.6 +/- 0.2, over the 0.036--0.1 pc range in radius sampled in extinction. This power law slope is steeper than the value of p=2 for a singular isothermal sphere, the initial condition of the inside-out model for protostellar collapse. Including an additional extinction component along the line of sight further steepens the inferred profile. Fitting a Bonnor-Ebert sphere results in a super-critical value of the dimensionless radius xi_max=25 +/- 3. The unstable configuration of material may be related to the observed inward motions. The Bonnor-Ebert model matches the shape of the observed profile, but significantly underestimates the amount of extinction (by a factor of ~4). This discrepancy in normalization has also been found for the nearby protostellar core B335 (Harvey et al. 2001). A cylindrical density model with scale height H=0.0164+/- 0.002 pc viewed at a small inclination to the cylinder axis provides an equally good radial profile as a power law model, and reproduces the asymmetry of the core remarkably well. In addition, this model provides a basis for understanding the discrepancy in the normalization of the Bonnor-Ebert model, namely that L694-2 has prolate structure, with the full extent (mass) of the core being missed by assuming symmetry between the profiles in the plane of the sky and along the line-of-sight. If the core is sufficiently magnetized then fragmentation may be avoided, and later evolution might produce a protostar similar to B335. △ Less

Submitted 8 August, 2003; originally announced August 2003.

Comments: 38 pages, 7 figures, accepted to Astrophysical Journal

Journal ref: Astrophys.J.598:1112-1126,2003

Inner Structure of Starless Core L694--2 Derived from Millimeter-Wave Interferometry

Authors: Daniel W. A. Harvey, David J. Wilner, Philip C. Myers, Mario Tafalla

Abstract: We study the density structure of candidate contracting starless core L694-2 using 1.3 mm dust continuum observations from the IRAM Plateau de Bure Interferometer and the Berkeley-Illinois-Maryland Array, probing spatial scales from 10000-500 AU. The long baseline PdBI observations detect no emission, and limit the contamination from a compact component F_c < 2.7 mJy. This limit corresponds to a… ▽ More We study the density structure of candidate contracting starless core L694-2 using 1.3 mm dust continuum observations from the IRAM Plateau de Bure Interferometer and the Berkeley-Illinois-Maryland Array, probing spatial scales from 10000-500 AU. The long baseline PdBI observations detect no emission, and limit the contamination from a compact component F_c < 2.7 mJy. This limit corresponds to a very small disk mass, M_disk < 5e-4 M_sun x (60 K / T_disk), and bolsters the ``starless'' interpretation of the core. The shorter baseline BIMA data are compared to density models using a physically motivated temperature distribution with a central minimum. This analysis provides clear evidence for a turn-over from the steep outer density profile observed in dust extinction to much more shallow behavior in the inner regions (<7500 AU). The best fit Bonnor-Ebert, Plummer-like, broken power law, and end-on cylinder models produce very similar flattened profiles and cannot be distinguished. We quantify the sensitivity of the inferred structure to various uncertainties, including the temperature distribution, the central position, and the presence of a weak unresolved central component. The largest uncertainty comes from the temperature assumption; an isothermal model modifies the best fit parameters by \~2 sigma, with the inferred profiles more shallow. Dust emission and extinction profiles are reproduced by an embedded cylinder with scale height H=13.5'' inclined at a small angle to the line of sight. The flat central density profile suggests that pressure forces still support the core, despite the extended inward motions inferred from molecular line observations (Lee, Myers, & Tafalla 2001). In the context of the cylindrical density model, these inward motions may represent the contraction of a prolate core along its major axis. △ Less

Submitted 5 July, 2003; originally announced July 2003.

Comments: 26 pages, 2 figures. Accepted to the Astrophysical Journal

Journal ref: Astrophys.J.597:424-433,2003

Disk Properties and Density Structure of the Star-Forming Dense Core B335

Authors: Daniel W. A. Harvey, David J. Wilner, Philip C. Myers, Mario Tafalla

Abstract: We present subarcsecond resolution dust continuum observations of the protostellar collapse candidate B335 made with the IRAM Plateau de Bure Interferometer at wavelengths of 1.2 and 3.0 mm. These observations probe to < 100 AU size scales and reveal a compact source component that we identify with a circumstellar disk. We analyze these data in concert with previous lower resolution interferomet… ▽ More We present subarcsecond resolution dust continuum observations of the protostellar collapse candidate B335 made with the IRAM Plateau de Bure Interferometer at wavelengths of 1.2 and 3.0 mm. These observations probe to < 100 AU size scales and reveal a compact source component that we identify with a circumstellar disk. We analyze these data in concert with previous lower resolution interferometer observations and find a best fit density structure for B335 that consists of a power law envelope with index p=1.55 +/- 0.04 (r < 5000 AU) together with a disk (r < 100 AU) of flux F_{1.2 mm}=21 +/-2 mJy. We estimate a systematic uncertainty in the power law index delta(p) < 0.15, where the largest error comes from the assumed form of the dust temperature falloff with radius. This determination of the inner density structure of B335 has a precision unique amongst protostellar cores, and it is consistent with the r^{-1.5} profile of gravitational free-fall, in accord with basic expectations for the formation of a star. The flux (and implied mass) of the compact component in B335 is typical of the disks around T Tauri stars. △ Less

Submitted 5 July, 2003; originally announced July 2003.

Comments: 16 pages, 2 figures. Accepted to the Astrophysical Journal, sched v596 (2003 Oct 10)

Journal ref: Astrophys.J.596:383-388,2003

FS Aur - a new class of Cataclysmic Variables or the missing link between Intermediate Polars and SW Sex objects?

Authors: G. Tovmassian, S. Zharikov, R. Michel, ; V. Neustroev, J. Greiner, D. Skillman, D. Harvey, R. Fried, J. Patterson

Abstract: FS Aur is a known dwarf nova with an orbital period of about 85.7 minutes. It has been assumed to be a member of the SU UMa subclass of cataclysmic variables (CVs), but previous searches for superhumps and superoutburst have been unsuccessful. We conducted a series of photometric and spectroscopic observations of FS Aur during quiescence. We confirmed its short orbital period from radial velocit… ▽ More FS Aur is a known dwarf nova with an orbital period of about 85.7 minutes. It has been assumed to be a member of the SU UMa subclass of cataclysmic variables (CVs), but previous searches for superhumps and superoutburst have been unsuccessful. We conducted a series of photometric and spectroscopic observations of FS Aur during quiescence. We confirmed its short orbital period from radial velocity measurements. However, the long-term photometry revealed an unexpected result: the system also shows a distinct 0.24 mag modulation in the BVR photometric bands with a period of 205.5 min, which is 2.4 times longer than the orbital period. We discuss various possible causes for such a peculiar behavior. △ Less

Submitted 26 February, 2003; originally announced February 2003.

Comments: 12 pages, 1 table, 12 figures, accepted for publication in PASP