Showing 1–50 of 139 results for author: Barnes, D

L1 and off Sun-Earth line visible-light imaging of Earth-directed CMEs: An analysis of inconsistent observations

Authors: Richard A. Harrison, Jackie A. Davies, David Barnes, Christian Möstl

Abstract: The efficacy of coronal mass ejection (CME) observations as a key input to space weather forecasting is explored by comparing on and off Sun-Earth line observations from the ESA/NASA SOHO and NASA STEREO spacecraft. A comparison is made of CME catalogues based on L1 coronagraph imagery and off Sun-Earth line coronagraph and heliospheric imager (HI) observations, for the year 2011. Analysis reveals… ▽ More The efficacy of coronal mass ejection (CME) observations as a key input to space weather forecasting is explored by comparing on and off Sun-Earth line observations from the ESA/NASA SOHO and NASA STEREO spacecraft. A comparison is made of CME catalogues based on L1 coronagraph imagery and off Sun-Earth line coronagraph and heliospheric imager (HI) observations, for the year 2011. Analysis reveals inconsistencies in the identification of a number of potentially Earth-directed CMEs. The catalogues reflect our ability to identify and characterise CMEs, so any discrepancies can impact our prediction of Earth-directed CMEs. We show that 15 CMEs, which were observed by STEREO, that had estimated directions compatible with Earth-directed events, had no identified halo/partial halo counterpart listed in the L1 coronagraph CME catalogue. In-situ data confirms that for 9 of these there was a consistent L1 Interplanetary CME (ICME). The number of such "discrepant" events is significant compared to the number of ICMEs recorded at L1 in 2011, stressing the need to address space weather monitoring capabilities, particularly with the inclusion of off Sun-Earth line observation. While the study provides evidence that some halo CMEs are simply not visible in near-Earth coronagraph imagery, there is evidence that some halo CMEs viewed from L1 are compromised by preceding CME remnants or the presence of multiple-CMEs. This underlines (1) the value of multiple vantage point CME observation, and (2) the benefit of off Sun-Earth line platform heliospheric imaging, and coronagraph imaging, for the efficient identification and tracking of Earth-directed events. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 36 pages, 6 figures, in press at AGU Space Weather, 2023

Magnetic Structure and Propagation of Two Interacting CMEs from the Sun to Saturn

Authors: Erika Palmerio, Teresa Nieves-Chinchilla, Emilia K. J. Kilpua, David Barnes, Andrei N. Zhukov, Lan K. Jian, Olivier Witasse, Gabrielle Provan, Chihiro Tao, Laurent Lamy, Thomas J. Bradley, M. Leila Mays, Christian Möstl, Elias Roussos, Yoshifumi Futaana, Adam Masters, Beatriz Sánchez-Cano

Abstract: One of the grand challenges in heliophysics is the characterisation of coronal mass ejection (CME) magnetic structure and evolution from eruption at the Sun through heliospheric propagation. At present, the main difficulties are related to the lack of direct measurements of the coronal magnetic fields and the lack of 3D in-situ measurements of the CME body in interplanetary space. Nevertheless, th… ▽ More One of the grand challenges in heliophysics is the characterisation of coronal mass ejection (CME) magnetic structure and evolution from eruption at the Sun through heliospheric propagation. At present, the main difficulties are related to the lack of direct measurements of the coronal magnetic fields and the lack of 3D in-situ measurements of the CME body in interplanetary space. Nevertheless, the evolution of a CME magnetic structure can be followed using a combination of multi-point remote-sensing observations and multi-spacecraft in-situ measurements as well as modelling. Accordingly, we present in this work the analysis of two CMEs that erupted from the Sun on 28 April 2012. We follow their eruption and early evolution using remote-sensing data, finding indications of CME--CME interaction, and then analyse their interplanetary counterpart(s) using in-situ measurements at Venus, Earth, and Saturn. We observe a seemingly single flux rope at all locations, but find possible signatures of interaction at Earth, where high-cadence plasma data are available. Reconstructions of the in-situ flux ropes provide almost identical results at Venus and Earth but show greater discrepancies at Saturn, suggesting that the CME was highly distorted and/or that further interaction with nearby solar wind structures took place before 10 AU. This work highlights the difficulties in connecting structures from the Sun to the outer heliosphere and demonstrates the importance of multi-spacecraft studies to achieve a deeper understanding of the magnetic configuration of CMEs. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 40 pages, 11 figures, accepted for publication in Journal of Geophysical Research: Space Physics

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Authors: Dhayaa Anbajagane, Han Aung, August E. Evrard, Arya Farahi, Daisuke Nagai, David J. Barnes, Weiguang Cui, Klaus Dolag, Ian G. McCarthy, Elena Rasia, Gustavo Yepes

Abstract: Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the rat… ▽ More Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available. △ Less

Submitted 21 December, 2021; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: [v1]: 8 figures, 14 pages in main text + 5 figures, 3 pages in Appendix. [v2]: Minor text edits. Accepted in MNRAS

Multipoint interplanetary coronal mass ejections observed with Solar Orbiter, BepiColombo, Parker Solar Probe, Wind and STEREO-A

Authors: C. Möstl, A. J. Weiss, M. A. Reiss, T. Amerstorfer, R. L. Bailey, J. Hinterreiter, M. Bauer, D. Barnes, J. A. Davies, R. A. Harrison, J. L. Freiherr von Forstner, E. E. Davies, D. Heyner, T. Horbury, S. D. Bale

Abstract: We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP),… ▽ More We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP), Wind and STEREO-A, in connection with coronagraph and heliospheric imaging observations from STEREO-A/SECCHI and SOHO/LASCO. We identify ICME events that could be unambiguously followed with the STEREO-A heliospheric imagers during their interplanetary propagation to their impact at the aforementioned spacecraft, and look for events where the same ICME is seen in situ by widely separated spacecraft. We highlight two events: (1) a small streamer blowout CME on 2020 June 23 observed with a triple lineup by PSP, BepiColombo and Wind, guided by imaging with STEREO-A, and (2) the first fast CME of solar cycle 25 ($ \approx 1600$ km s$^{-1}$) on 2020 November 29 observed in situ by PSP and STEREO-A. These results are useful for modeling the magnetic structure of ICMEs and the interplanetary evolution and global shape of their flux ropes and shocks, and for studying the propagation of solar energetic particles. The combined data from these missions are already turning out to be a treasure trove for space weather research and are expected to become even more valuable with an increasing number of ICME events expected during the rise and maximum of solar cycle 25. △ Less

Submitted 4 January, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: in press at ApJ Letters (submitted 2021 September 15, revised 2021 November 24, accepted 2021 December 9). 11 pages, 3 figures, 1 table

CME Magnetic Structure and IMF Preconditioning Affecting SEP Transport

Authors: Erika Palmerio, Emilia K. J. Kilpua, Olivier Witasse, David Barnes, Beatriz Sánchez-Cano, Andreas J. Weiss, Teresa Nieves-Chinchilla, Christian Möstl, Lan K. Jian, Marilena Mierla, Andrei N. Zhukov, Jingnan Guo, Luciano Rodriguez, Patrick J. Lowrance, Alexey Isavnin, Lucile Turc, Yoshifumi Futaana, Mats Holmström

Abstract: Coronal mass ejections (CMEs) and solar energetic particles (SEPs) are two phenomena that can cause severe space weather effects throughout the heliosphere. The evolution of CMEs, especially in terms of their magnetic structure, and the configuration of the interplanetary magnetic field (IMF) that influences the transport of SEPs are currently areas of active research. These two aspects are not ne… ▽ More Coronal mass ejections (CMEs) and solar energetic particles (SEPs) are two phenomena that can cause severe space weather effects throughout the heliosphere. The evolution of CMEs, especially in terms of their magnetic structure, and the configuration of the interplanetary magnetic field (IMF) that influences the transport of SEPs are currently areas of active research. These two aspects are not necessarily independent of each other, especially during solar maximum when multiple eruptive events can occur close in time. Accordingly, we present the analysis of a CME that erupted on 2012 May 11 (SOL2012-05-11) and an SEP event following an eruption that took place on 2012 May 17 (SOL2012-05-17). After observing the May 11 CME using remote-sensing data from three viewpoints, we evaluate its propagation through interplanetary space using several models. Then, we analyse in-situ measurements from five predicted impact locations (Venus, Earth, the Spitzer Space Telescope, the Mars Science Laboratory en route to Mars, and Mars) in order to search for CME signatures. We find that all in-situ locations detect signatures of an SEP event, which we trace back to the May 17 eruption. These findings suggest that the May 11 CME provided a direct magnetic connectivity for the efficient transport of SEPs. We discuss the space weather implications of CME evolution, regarding in particular its magnetic structure, and CME-driven IMF preconditioning that facilitates SEP transport. Finally, this work remarks the importance of using data from multiple spacecraft, even those that do not include space weather research as their primary objective. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 50 pages, 14 figures, 2 tables, accepted for publication in Space Weather

In situ multi-spacecraft and remote imaging observations of the first CME detected by Solar Orbiter and BepiColombo

Authors: E. E. Davies, C. Möstl, M. J. Owens, A. J. Weiss, T. Amerstorfer, J. Hinterreiter, M. Bauer, R. L. Bailey, M. A. Reiss, R. J. Forsyth, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini, D. Heyner, I. Richter, H-U. Auster, W. Magnes, W. Baumjohann, D. Fischer, D. Barnes, J. A. Davies, R. A. Harrison

Abstract: On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ i… ▽ More On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ in longitude. The source of the CME, which was launched on April 15th, was an almost entirely isolated streamer blowout. STEREO-A observed the event remotely from -75.1$^\circ$ longitude, which is an exceptionally well suited viewpoint for heliospheric imaging of an Earth directed CME. The configuration of the four spacecraft has provided an exceptionally clean link between remote imaging and in situ observations of the CME. We have used the in situ observations of the CME at Solar Orbiter, Wind, and BepiColombo, and the remote observations of the CME at STEREO-A in combination with flux rope models to determine the global shape of the CME and its evolution as it propagated through the inner heliosphere. A clear flattening of the CME cross-section has been observed by STEREO-A, and further confirmed by comparing profiles of the flux rope models to the in situ data, where the distorted flux rope cross-section qualitatively agrees most with in situ observations of the magnetic field at Solar Orbiter. Comparing in situ observations of the magnetic field between spacecraft, we find that the dependence of the maximum (mean) magnetic field strength decreases with heliocentric distance as $r^{-1.24 \pm 0.50}$ ($r^{-1.12 \pm 0.14}$), in disagreement with previous studies. Further assessment of the axial and poloidal magnetic field strength dependencies suggests that the expansion of the CME is likely neither self-similar nor cylindrically symmetric. △ Less

Submitted 23 February, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Journal ref: A&A 656, A2 (2021)

The splashback boundary of haloes in hydrodynamic simulations

Authors: Stephanie O'Neil, David J. Barnes, Mark Vogelsberger, Benedikt Diemer

Abstract: The splashback radius, $R_{\rm sp}$, is a physically motivated halo boundary that separates infalling and collapsed matter of haloes. We study $R_{\rm sp}$ in the hydrodynamic and dark matter only IllustrisTNG simulations. The most commonly adopted signature of $R_{\rm sp}$ is the radius at which the radial density profiles are steepest. Therefore, we explicitly optimise our density profile fit to… ▽ More The splashback radius, $R_{\rm sp}$, is a physically motivated halo boundary that separates infalling and collapsed matter of haloes. We study $R_{\rm sp}$ in the hydrodynamic and dark matter only IllustrisTNG simulations. The most commonly adopted signature of $R_{\rm sp}$ is the radius at which the radial density profiles are steepest. Therefore, we explicitly optimise our density profile fit to the profile slope and find that this leads to a $\sim5\%$ larger radius compared to other optimisations. We calculate $R_{\rm sp}$ for haloes with masses between $10^{13-15}{\rm M}_{\odot}$ as a function of halo mass, accretion rate and redshift. $R_{\rm sp}$ decreases with mass and with redshift for haloes of similar $M_{\rm200m}$ in agreement with previous work. We also find that $R_{\rm sp}/R_{\rm200m}$ decreases with halo accretion rate. We apply our analysis to dark matter, gas and satellite galaxies associated with haloes to investigate the observational potential of $R_{\rm sp}$. The radius of steepest slope in gas profiles is consistently smaller than the value calculated from dark matter profiles. The steepest slope in galaxy profiles, which are often used in observations, tends to agree with dark matter profiles but is lower for less massive haloes. We compare $R_{\rm sp}$ in hydrodynamic and N-body dark matter only simulations and do not find a significant difference caused by the addition of baryonic physics. Thus, results from dark matter only simulations should be applicable to realistic haloes. △ Less

Submitted 25 May, 2021; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 18 pages, 15 figures, accepted by MNRAS

SHAPing the Gas: Understanding Gas Shapes in Dark Matter Haloes with Interpretable Machine Learning

Authors: Luis Fernando Machado Poletti Valle, Camille Avestruz, David J. Barnes, Arya Farahi, Erwin T. Lau, Daisuke Nagai

Abstract: The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machin… ▽ More The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machine learning approach for modelling the dependence of gas shapes on dark matter and baryonic properties. With data from the IllustrisTNG hydrodynamical cosmological simulations, we develop a machine learning pipeline that applies \pkg{XGBoost}, an implementation of gradient boosted decision trees, to predict radial profiles of gas shapes from halo properties. We show that \pkg{XGBoost} models can accurately predict gas shape profiles in dark matter haloes. We also explore model interpretability with \pkg{SHAP}, a method that identifies the most predictive properties at different halo radii. We find that baryonic properties best predict gas shapes in halo cores, whereas dark matter shapes are the main predictors in the halo outskirts. This work demonstrates the power of interpretable machine learning in modelling observable properties of dark matter haloes in the era of multi-wavelength cosmological surveys. △ Less

Submitted 2 August, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 14 pages, 8 figures, Accepted to MNRAS

The Cluster HEritage project with XMM-Newton: Mass Assembly and Thermodynamics at the Endpoint of structure formation. I. Programme overview

Authors: The CHEX-MATE Collaboration, :, M. Arnaud, S. Ettori, G. W. Pratt, M. Rossetti, D. Eckert, F. Gastaldello, R. Gavazzi, S. T. Kay, L. Lovisari, B. J. Maughan, E. Pointecouteau, M. Sereno, I. Bartalucci, A. Bonafede, H. Bourdin, R. Cassano, R. T. Duffy, A. Iqbal, S. Maurogordato, E. Rasia, J. Sayers, F. Andrade-Santos, H. Aussel , et al. (45 additional authors not shown)

Abstract: The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise limited sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. The programme, described in detail in this paper, aim… ▽ More The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise limited sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. The programme, described in detail in this paper, aims to study the ultimate products of structure formation in time and mass. It is composed of a census of the most recent objects to have formed (Tier-1: 0.05 < z < 0.2; 2 x 10e14 M_sun < M_500 < 9 x 10e14 M_sun), together with a sample of the highest-mass objects in the Universe (Tier-2: z < 0.6; M_500 > 7.25 x 10e14 M_sun). The programme will yield an accurate vision of the statistical properties of the underlying population, measure how the gas properties are shaped by collapse into the dark matter halo, uncover the provenance of non-gravitational heating, and resolve the major uncertainties in mass determination that limit the use of clusters for cosmological parameter estimation. We will acquire X-ray exposures of uniform depth, designed to obtain individual mass measurements accurate to 15-20% under the hydrostatic assumption. We present the project motivations, describe the programme definition, and detail the ongoing multi-wavelength observational (lensing, SZ, radio) and theoretical effort that is being deployed in support of the project. △ Less

Submitted 3 March, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 27 pages, 11 figures; A&A, in press

Journal ref: A&A 650, A104 (2021)

Is there enough star formation in simulated protoclusters?

Authors: Seunghwan Lim, Douglas Scott, Arif Babul, David Barnes, Scott Kay, Ian McCarthy, Douglas Rennehan, Mark Vogelsberger

Abstract: As progenitors of the most massive objects, protoclusters are key to tracing the evolution and star-formation history of the Universe, and are responsible for ${\gtrsim}\,20$ per cent of the cosmic star formation at $z\,{>}\,2$. Using a combination of state-of-the-art hydrodynamical simulations and empirical models, we show that current galaxy-formation models do not produce enough star formation… ▽ More As progenitors of the most massive objects, protoclusters are key to tracing the evolution and star-formation history of the Universe, and are responsible for ${\gtrsim}\,20$ per cent of the cosmic star formation at $z\,{>}\,2$. Using a combination of state-of-the-art hydrodynamical simulations and empirical models, we show that current galaxy-formation models do not produce enough star formation in protoclusters to match observations. We find that the star-formation rates (SFRs) predicted from the models are an order of magnitude lower than what is seen in observations, despite the relatively good agreement found for their mass-accretion histories, specifically that they lie on an evolutionary path to become Coma-like clusters at $z\,{\simeq}\, 0$. Using a well-studied protocluster core at $z\,{=}\,4.3$ as a test case, we find that star-formation efficiency of protocluster galaxies is higher than predicted by the models. We show that a large part of the discrepancy can be attributed to a dependence of SFR on the numerical resolution of the simulations, with a roughly factor of 3 drop in SFR when the spatial resolution decreases by a factor of 4. We also present predictions up to $z\,{\simeq}\,7$. Compared to lower redshifts, we find that centrals (the most massive member galaxies) are more distinct from the other galaxies, while protocluster galaxies are less distinct from field galaxies. All these results suggest that, as a rare and extreme population at high-$z$, protoclusters can help constrain galaxy formation models tuned to match the average population at $z\,{\simeq}\,0$. △ Less

Submitted 24 November, 2020; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: accepted for publication

The high helium abundance and charge states of the interplanetary CME and its material source on the Sun

Authors: Hui Fu, R. A. Harrison, J. A. Davies, LiDong Xia, XiaoShuai Zhu, Bo Li, ZhengHua Huang, D. Barnes

Abstract: Identifying the source of the material within coronal mass ejections (CMEs) and understanding CME onset mechanisms are fundamental issues in solar and space physics. Parameters relating to plasma composition, such as charge states and He abundance (\ahe), may be different for plasmas originating from differing processes or regions on the Sun. Thus, it is crucial to examine the relationship between… ▽ More Identifying the source of the material within coronal mass ejections (CMEs) and understanding CME onset mechanisms are fundamental issues in solar and space physics. Parameters relating to plasma composition, such as charge states and He abundance (\ahe), may be different for plasmas originating from differing processes or regions on the Sun. Thus, it is crucial to examine the relationship between in-situ measurements of CME composition and activity on the Sun. We study the CME that erupted on 2014 September 10, in association with an X1.6 flare, by analyzing AIA imaging and IRIS spectroscopic observations and its in-situ signatures detected by Wind and ACE. We find that during the slow expansion and intensity increase of the sigmoid, plasma temperatures of 9 MK, and higher, first appear at the footpoints of the sigmoid, associated with chromospheric brightening. Then the high-temperature region extends along the sigmoid. IRIS observations confirm that this extension is caused by transportation of hot plasma upflow. Our results show that chromospheric material can be heated to 9 MK, and above, by chromospheric evaporation at the sigmoid footpoints before flare onset. The heated chromospheric material can transport into the sigmoidal structure and supply mass to the CME. The aforementioned CME mass supply scenario provides a reasonable explanation for the detection of high charge states and elevated \ahe\ in the associated ICME. The observations also demonstrate that the quasi-steady evolution in the precursor phase is dominated by magnetic reconnection between the rising flux rope and the overlying magnetic field structure. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 10 pages, 5 figures, accepted for publication in ApJL

First Light And Reionisation Epoch Simulations (FLARES) II: The Photometric Properties of High-Redshift Galaxies

Authors: Aswin P. Vijayan, Christopher C. Lovell, Stephen M. Wilkins, Peter A. Thomas, David J. Barnes, Dimitrios Irodotou, Jussi Kuusisto, Will Roper

Abstract: We present the photometric properties of galaxies in the First Light and Reionisation Epoch Simulations (FLARES). The simulations trace the evolution of galaxies in a range of overdensities through the Epoch of Reionistion (EoR). With a novel weighting scheme we combine these overdensities, extending significantly the dynamic range of observed composite distribution functions compared to periodic… ▽ More We present the photometric properties of galaxies in the First Light and Reionisation Epoch Simulations (FLARES). The simulations trace the evolution of galaxies in a range of overdensities through the Epoch of Reionistion (EoR). With a novel weighting scheme we combine these overdensities, extending significantly the dynamic range of observed composite distribution functions compared to periodic simulation boxes. FLARES predicts a significantly larger number of intrinsically bright galaxies, which can be explained through a simple model linking dust-attenuation to the metal content of the interstellar medium, using a line-of-sight (LOS) extinction model. With this model we present the photometric properties of the FLARES galaxies for $z \in [5,10]$. We show that the ultraviolet (UV) luminosity function (LF) matches the observations at all redshifts. The function is fit by Schechter and double power-law forms, with the latter being favoured at these redshifts by the FLARES composite UV LF. We also present predictions for the UV continuum slope as well as the attenuation in the UV. The impact of environment on the UV LF is also explored, with the brightest galaxies forming in the densest environments. We then present the line luminosity and equivalent widths of some prominent nebular emission lines arising from the galaxies, finding rough agreement with available observations. We also look at the relative contribution of obscured and unobscured star formation, finding comparable contributions at these redshifts. △ Less

Submitted 14 December, 2020; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 21 pages, 22 figures. Data available at https://flaresimulations.github.io/data.html

Properties of the ionized CGM and IGM: a comparison of galaxy formation models via the Sunyaev-Zel'dovich effect

Authors: S. H. Lim, David Barnes, Mark Vogelsberger, H. J. Mo, Dylan Nelson, Annalisa Pillepich, Klaus Dolag, Federico Marinacci

Abstract: We present a comparison of the physical properties of the ionized gas in the circumgalactic (CGM) and intergalactic (IGM) media at $z\sim0$ between observations and four cosmological hydrodynamical simulations: Illustris, TNG300 of the IllustrisTNG project, EAGLE, and one of the Magneticum simulations. For the observational data, we use the gas properties that are inferred from cross-correlating t… ▽ More We present a comparison of the physical properties of the ionized gas in the circumgalactic (CGM) and intergalactic (IGM) media at $z\sim0$ between observations and four cosmological hydrodynamical simulations: Illustris, TNG300 of the IllustrisTNG project, EAGLE, and one of the Magneticum simulations. For the observational data, we use the gas properties that are inferred from cross-correlating the Sunyaev-Zel'dovich effect (SZE) from the {\it Planck} CMB maps with the haloes and the large-scale structure reconstructed from Sloan Digital Sky Survey data. Both the observational and simulation results indicate that the integrated gas pressure in haloes deviates from the self-similar case, showing that feedback impacts haloes with $M_{500}\sim 10^{12-13}\,{\rm M_\odot}$. The simulations predict that more than half the baryons are displaced from haloes, while the gas fraction inferred from our observational data roughly equals the cosmic baryon fraction throughout the $M_{500}\sim 10^{12-14.5}\,{\rm M_\odot}$ halo mass range. All simulations tested here predict that the mean gas temperature in haloes is about the virial temperature, while that inferred from the SZE is up to one order of magnitude lower than that from the simulations (and also from X-ray observations). While a remarkable agreement is found for the average properties of the IGM between the observation and some simulations, we show that their dependence on the large-scale tidal field can break the degeneracy between models that show similar predictions otherwise. Finally, we show that the gas pressure and the electron density profiles from simulations are not well described by a generalized NFW (GNFW) profile. Instead, we present a new model with a mass-dependent shape that fits the profiles accurately. △ Less

Submitted 16 April, 2021; v1 submitted 22 July, 2020; originally announced July 2020.

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

A deep learning view of the census of galaxy clusters in IllustrisTNG

Authors: Y. Su, Y. Zhang, G. Liang, J. A. ZuHone, D. J. Barnes, N. B. Jacobs, M. Ntampaka, W. R. Forman, P. E. J. Nulsen, R. P. Kraft, C. Jones

Abstract: The origin of the diverse population of galaxy clusters remains an unexplained aspect of large-scale structure formation and cluster evolution. We present a novel method of using X-ray images to identify cool core (CC), weak cool core (WCC), and non cool core (NCC) clusters of galaxies, that are defined by their central cooling times. We employ a convolutional neural network, ResNet-18, which is c… ▽ More The origin of the diverse population of galaxy clusters remains an unexplained aspect of large-scale structure formation and cluster evolution. We present a novel method of using X-ray images to identify cool core (CC), weak cool core (WCC), and non cool core (NCC) clusters of galaxies, that are defined by their central cooling times. We employ a convolutional neural network, ResNet-18, which is commonly used for image analysis, to classify clusters. We produce mock Chandra X-ray observations for a sample of 318 massive clusters drawn from the IllustrisTNG simulations. The network is trained and tested with low resolution mock Chandra images covering a central 1 Mpc square for the clusters in our sample. Without any spectral information, the deep learning algorithm is able to identify CC, WCC, and NCC clusters, achieving balanced accuracies (BAcc) of 92%, 81%, and 83%, respectively. The performance is superior to classification by conventional methods using central gas densities, with an average BAcc = 81%, or surface brightness concentrations, giving BAcc = 73%. We use Class Activation Mapping to localize discriminative regions for the classification decision. From this analysis, we observe that the network has utilized regions from cluster centers out to r~300 kpc and r~500 kpc to identify CC and NCC clusters, respectively. It may have recognized features in the intracluster medium that are associated with AGN feedback and disruptive major mergers. △ Less

Submitted 25 August, 2020; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS

CMEs in the Heliosphere: III. A Statistical Analysis of the Kinematic Properties Derived from Stereoscopic Geometrical Modelling Techniques Applied to CMEs Detected in the Heliosphere from 2008 to 2014 by STEREO/HI-1

Authors: D. Barnes, J. A. Davies, R. A. Harrison, J. P. Byrne, C. H. Perry, V. Bothmer, J. P. Eastwood, P. T. Gallagher, E. K. J. Kilpua, C. Möstl, L. Rodriguez, A. P. Rouillard, D. Odstrcil

Abstract: We present an analysis of coronal mass ejections (CMEs) observed by the Heliospheric Imagers (HIs) on board NASA's Solar Terrestrial Relations Observatory (STEREO) spacecraft. Between August 2008 and April 2014 we identify 273 CMEs that are observed simultaneously, by the HIs on both spacecraft. For each CME, we track the observed leading edge, as a function of time, from both vantage points, and… ▽ More We present an analysis of coronal mass ejections (CMEs) observed by the Heliospheric Imagers (HIs) on board NASA's Solar Terrestrial Relations Observatory (STEREO) spacecraft. Between August 2008 and April 2014 we identify 273 CMEs that are observed simultaneously, by the HIs on both spacecraft. For each CME, we track the observed leading edge, as a function of time, from both vantage points, and apply the Stereoscopic Self-Similar Expansion (SSSE) technique to infer their propagation throughout the inner heliosphere. The technique is unable to accurately locate CMEs when their observed leading edge passes between the spacecraft, however, we are able to successfully apply the technique to 151, most of which occur once the spacecraft separation angle exceeds 180 degrees, during solar maximum. We find that using a small half-width to fit the CME can result in observed acceleration to unphysically high velocities and that using a larger half-width can fail to accurately locate the CMEs close to the Sun because the method does not account for CME over-expansion in this region. Observed velocities from SSSE are found to agree well with single-spacecraft (SSEF) analysis techniques applied to the same events. CME propagation directions derived from SSSE and SSEF analysis agree poorly because of known limitations present in the latter. This work was carried out as part of the EU FP7 HELCATS (Heliospheric Cataloguing, Analysis and Techniques Service) project (http://www.helcats-fp7.eu/). △ Less

Submitted 26 June, 2020; originally announced June 2020.

Studying galaxy cluster morphological metrics with Mock-X

Authors: Kaili Cao, David J. Barnes, Mark Vogelsberger

Abstract: Dynamically relaxed galaxy clusters have long played a role in galaxy cluster studies because it is thought their properties can be reconstructed more precisely and with less systematics. As relaxed clusters are desirable, there exist a plethora of criteria for classifying a galaxy cluster as relaxed. In this work, we examine $9$ commonly used observational and theoretical morphological metrics ex… ▽ More Dynamically relaxed galaxy clusters have long played a role in galaxy cluster studies because it is thought their properties can be reconstructed more precisely and with less systematics. As relaxed clusters are desirable, there exist a plethora of criteria for classifying a galaxy cluster as relaxed. In this work, we examine $9$ commonly used observational and theoretical morphological metrics extracted from $54,000$ Mock-X synthetic X-ray images of galaxy clusters taken from the IllustrisTNG, BAHAMAS and MACSIS simulation suites. We find that the simulated criteria distributions are in reasonable agreement with the observed distributions. Many criteria distributions evolve as a function of redshift, cluster mass, numerical resolution and subgrid physics, limiting the effectiveness of a single relaxation threshold value. All criteria are positively correlated with each other, however, the strength of the correlation is sensitive to redshift, mass and numerical choices. Driven by the intrinsic scatter inherent to all morphological metrics and the arbitrary nature of relaxation threshold values, we find the consistency of relaxed subsets defined by the different metrics to be relatively poor. Therefore, the use of relaxed cluster subsets introduces significant selection effects that are non-trivial to resolve. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: 21 pages, 6 figures, submitted to MNRAS. Comments welcome

Redshift evolution of the hot intracluster gas metallicity in the C-EAGLE cluster simulations

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Yannick M. Bahe, Richard G. Bower

Abstract: The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investi… ▽ More The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, submitted to MNRAS

First Light And Reionisation Epoch Simulations (FLARES) I: Environmental Dependence of High-Redshift Galaxy Evolution

Authors: Christopher C. Lovell, Aswin P. Vijayan, Peter A. Thomas, Stephen M. Wilkins, David J. Barnes, Dimitrios Irodotou, Will Roper

Abstract: We introduce the First Light And Reionisation Epoch Simulations (FLARES), a suite of zoom simulations using the EAGLE model. We resimulate a range of overdensities during the Epoch of Reionisation (EoR) in order to build composite distribution functions, as well as explore the environmental dependence of galaxy formation and evolution during this critical period of galaxy assembly. The regions are… ▽ More We introduce the First Light And Reionisation Epoch Simulations (FLARES), a suite of zoom simulations using the EAGLE model. We resimulate a range of overdensities during the Epoch of Reionisation (EoR) in order to build composite distribution functions, as well as explore the environmental dependence of galaxy formation and evolution during this critical period of galaxy assembly. The regions are selected from a large $(3.2 \;\mathrm{cGpc})^{3}$ parent volume, based on their overdensity within a sphere of radius $14\,h^{-1}\;\mathrm{cMpc}$. We then resimulate with full hydrodynamics, and employ a novel weighting scheme that allows the construction of composite distribution functions that are representative of the full parent volume. This significantly extends the dynamic range compared to smaller volume periodic simulations. We present an analysis of the galaxy stellar mass function (GSMF), the star formation rate distribution function (SFRF) and the star forming sequence (SFS) predicted by \flares, and compare to a number of observational and model constraints. We also analyse the environmental dependence over an unprecedented range of overdensity. Both the GSMF and the SFRF exhibit a clear double-Schechter form, up to the highest redshifts ($z = 10$). We also find no environmental dependence of the SFS normalisation. The increased dynamic range probed by FLARES will allow us to make predictions for a number of large area surveys that will probe the EoR in coming years, such as WFIRST and Euclid. △ Less

Submitted 2 September, 2020; v1 submitted 15 April, 2020; originally announced April 2020.

Comments: 21 pages, 19 figures, submitted to MNRAS

The intra-cluster light as a tracer of the total matter density distribution: a view from simulations

Authors: Isaac Alonso Asensio, Claudio Dalla Vecchia, Yannick M. Bahé, David J. Barnes, Scott T. Kay

Abstract: By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-den… ▽ More By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-density contours using the procedure of Montes \& Trujillo (2019), and find that the shape of the stellar mass distribution follows that of the total matter even more closely than observed, although their radial profiles differ substantially. The ratio between stellar and total matter density profiles in circular apertures, shows a slope close to -1, with a small dependence on the cluster's total mass. We propose an indirect method to calculate the halo mass and mass density profile from the radial profile of the intra-cluster stellar mass density. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 7 pages, 5 figures, submitted to MNRAS

Characterizing hydrostatic mass bias with Mock-X

Authors: David J. Barnes, Mark Vogelsberger, Francesca A. Pearce, Ana-Roxana Pop, Rahul Kannan, Kaili Cao, Scott T. Kay, Lars Hernquist

Abstract: Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we co… ▽ More Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the \textsc{Mock-X} analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-$y$ and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias $b\approx0.13$ with cluster mass, independent of hydrodynamical method or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to $b=0.3$ for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 20 pages, 9 figures, submitted to MNRAS

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

Authors: Dhayaa Anbajagane, August E. Evrard, Arya Farahi, David J. Barnes, Klaus Dolag, Ian G. McCarthy, Dylan Nelson, Annalisa Pillepich

Abstract: We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generate… ▽ More We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $γ= -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $γ\simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 17 pages, 15 figures

Relativistic SZ temperature scaling relations of groups and clusters derived from the BAHAMAS and MACSIS simulations

Authors: Elizabeth Lee, Jens Chluba, Scott T. Kay, David J. Barnes

Abstract: The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (… ▽ More The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (tSZ) effect. The third measure determines the lowest order relativistic correction to the tSZ signal, which is seeing increased observational relevance. Our procedure allows us to accurately model the relativistic SZ (rSZ) contribution and we show that a $\gtrsim 10\%-40\%$ underestimation of this rSZ cluster temperature is expected when applying standard X-ray relations. The correction also exhibits significant mass and redshift evolution, as we demonstrate here. We present the mass dependence of each temperature measure alongside their profiles and a short analysis of the temperature dispersion as derived from the aforementioned simulations. We also discuss a new relation connecting the temperature and Compton-$y$ parameter, which can be directly used for rSZ modelling. Simple fits to the obtained scaling relations and profiles are provided. These should be useful for future studies of the rSZ effect and its relevance to cluster cosmology. △ Less

Submitted 2 March, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: Accepted by MNRAS

CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

Authors: Ricardo Herbonnet, Cristóbal Sifón, Henk Hoekstra, Yannick Bahé, Remco F. J. van der Burg, Jean-Baptiste Melin, Anja von der Linden, David Sand, Scott Kay, David Barnes

Abstract: Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and r… ▽ More Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$σ$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 20 pages, 9 figures, 4 tables, submitted to MNRAS

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

Hydrostatic mass estimates of massive galaxy clusters: a study with varying hydrodynamics flavours and non-thermal pressure support

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Richard G. Bower, Matthieu Schaller

Abstract: We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but s… ▽ More We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 21 pages, 14 figures. Accepted for publication by MNRAS

Early-type galaxy density profiles from IllustrisTNG: II. Evolutionary trend of the total density profile

Authors: Yunchong Wang, Mark Vogelsberger, Dandan Xu, Xuejian Shen, Shude Mao, David Barnes, Hui Li, Federico Marinacci, Paul Torrey, Volker Springel, Lars Hernquist

Abstract: We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $γ^{\prime}$ of the total density profile for their main progenitors. We find that their $γ^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost const… ▽ More We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $γ^{\prime}$ of the total density profile for their main progenitors. We find that their $γ^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost constant, aside from a residual trend of becoming shallower towards $z=0$. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since $z=3$, and after $z\approx 1$, they remain nearly invariant with a mild increase towards $z=0$. These evolutionary trends are mass-dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analyzing the impact of mergers and AGN feedback on the progenitors' evolution, we conjecture a multi-phase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen $γ^{\prime}$ from $z=4$ to $z=2$, whereas subsequent AGN feedback (especially in the kinetic mode) makes $γ^{\prime}$ shallower again from $z=2$ to $z=1$. Afterwards, passive evolution from $z=1$ to $z=0$, mainly through gas-poor mergers, mildly decreases $γ^{\prime}$ and maintains the overall mass distribution close to isothermal. △ Less

Submitted 27 November, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 19 pages, 8 figures, 4 tables. Accepted for publication in MNRAS, updated to match published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 490, Issue 4, December 2019, Pages 5722-5738

Constraining the inner density slope of massive galaxy clusters

Authors: Qiuhan He, Hongyu Li, Ran Li, Carlos S. Frenk, Matthieu Schaller, David Barnes, Yannick Bahé, Scott T. Kay, Liang Gao, Claudio Dalla Vecchia

Abstract: We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from… ▽ More We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, $γ_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of $γ_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of $γ_{\rm gNFW}$. △ Less

Submitted 21 June, 2020; v1 submitted 2 July, 2019; originally announced July 2019.

Comments: 19 pages, 16 figures, accepted by MNRAS

Multipoint study of successive coronal mass ejections driving moderate disturbances at 1 AU

Authors: Erika Palmerio, Camilla Scolini, David Barnes, Jasmina Magdalenić, Matthew J. West, Andrei N. Zhukov, Luciano Rodriguez, Marilena Mierla, Simon W. Good, Diana E. Morosan, Emilia K. J. Kilpua, Jens Pomoell, Stefaan Poedts

Abstract: We analyse in this work the propagation and geoeffectiveness of four successive coronal mass ejections (CMEs) that erupted from the Sun during 21--23 May 2013 and that were detected in interplanetary space by the Wind and/or STEREO-A spacecraft. All these CMEs featured critical aspects for understanding so-called "problem space weather storms" at Earth. In the first three events a limb CMEs result… ▽ More We analyse in this work the propagation and geoeffectiveness of four successive coronal mass ejections (CMEs) that erupted from the Sun during 21--23 May 2013 and that were detected in interplanetary space by the Wind and/or STEREO-A spacecraft. All these CMEs featured critical aspects for understanding so-called "problem space weather storms" at Earth. In the first three events a limb CMEs resulted in moderately geoeffective in-situ structures at their target location in terms of the disturbance storm time (Dst) index (either measured or estimated). The fourth CME, which also caused a moderate geomagnetic response, erupted from close to the disc centre as seen from Earth, but it was not visible in coronagraph images from the spacecraft along the Sun--Earth line and appeared narrow and faint from off-angle viewpoints. Making the correct connection between CMEs at the Sun and their in-situ counterparts is often difficult for problem storms. We investigate these four CMEs using multiwavelength and multipoint remote-sensing observations (extreme ultraviolet, white light, and radio), aided by 3D heliospheric modelling, in order to follow their propagation in the corona and in interplanetary space and to assess their impact at 1 AU. Finally, we emphasise the difficulties in forecasting moderate space weather effects provoked by problematic and ambiguous events and the importance of multispacecraft data for observing and modelling problem storms. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: 20 pages, 10 figures, 3 tables, accepted for publication in The Astrophysical Journal

Simulating the dark matter decay signal from the Perseus galaxy cluster

Authors: Mark R. Lovell, Dmytro Iakubovskyi, David Barnes, Sownak Bose, Carlos S. Frenk, Tom Theuns, Wojciech A. Hellwing

Abstract: The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show tha… ▽ More The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show that the line width in C-EAGLE clusters similar to Perseus is typically [600-800]$\mathrm{kms^{-1}}$, and therefore narrower than the amplitude of the velocity dispersion of galaxies in the cluster. Haloes that are significantly disturbed can, however, exhibit galaxy velocity dispersions higher than $1000\mathrm{kms^{-1}}$, and in this case will show a large difference between the line profiles of on- and off-center observations. We show that the line profile is likely to be slightly asymmetric, but still well approximated by a Gaussian at the 10% level, and that the halo asymmetry can lead to fluxes that vary by a factor of two. In summary, we predict that, if the previously reported 3.55keV line detections do originate from dark matter decay, the XRISM mission will detect a line with a roughly Gaussian profile at a rest frame energy of 3.55keV, with a width $>600\mathrm{kms^{-1}}$ and flux approximately in the range $[4-9]\times10^{-8}\mathrm{counts/sec/cm^{2}}$. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: 6 pages and 3 figures, plus 5 pages and 6 figures in supplemental material. Accepted for publication in ApJL

Covariant polarized radiative transfer on cosmological scales for investigating large-scale magnetic field structures

Authors: Jennifer Y. H. Chan, Kinwah Wu, Alvina Y. L. On, David J. Barnes, Jason D. McEwen, Thomas D. Kitching

Abstract: Polarization of radiation is a powerful tool to study cosmic magnetism and analysis of polarization can be used as a diagnostic tool for large-scale structures. In this paper, we present a solid theoretical foundation for using polarized light to investigate large-scale magnetic field structures: the cosmological polarized radiative transfer (CPRT) formulation. The CPRT formulation is fully covari… ▽ More Polarization of radiation is a powerful tool to study cosmic magnetism and analysis of polarization can be used as a diagnostic tool for large-scale structures. In this paper, we present a solid theoretical foundation for using polarized light to investigate large-scale magnetic field structures: the cosmological polarized radiative transfer (CPRT) formulation. The CPRT formulation is fully covariant. It accounts for cosmological and relativistic effects in a self-consistent manner and explicitly treats Faraday rotation, as well as Faraday conversion, emission, and absorption processes. The formulation is derived from the first principles of conservation of phase-space volume and photon number. Without loss of generality, we consider a flat Friedmann-Robertson-Walker (FRW) space-time metric and construct the corresponding polarized radiative transfer equations. We propose an all-sky CPRT calculation algorithm, based on a ray-tracing method, which allows cosmological simulation results to be incorporated and, thereby, model templates of polarization maps to be constructed. Such maps will be crucial in our interpretation of polarized data, such as those to be collected by the Square Kilometer Array (SKA). We describe several tests which are used for verifying the code and demonstrate applications in the study of the polarization signatures in different distributions of electron number density and magnetic fields. We present a pencil-beam CPRT calculation and an all-sky calculation, using a simulated galaxy cluster or a model magnetized universe obtained from GCMHD+ simulations as the respective input structures. The implications on large-scale magnetic field studies are discussed; remarks on the standard methods using rotation measure are highlighted. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 32 pages, 14 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 484, Issue 2, 1 April 2019, Pages 1427-1455

Disruption of satellite galaxies in simulated groups and clusters: the roles of accretion time, baryons, and pre-processing

Authors: Yannick M. Bahé, Joop Schaye, David J. Barnes, Claudio Dalla Vecchia, Scott T. Kay, Richard G. Bower, Henk Hoekstra, Sean L. McGee, Tom Theuns

Abstract: We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters t… ▽ More We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes. △ Less

Submitted 1 February, 2019; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 26 pages, 24 figures (16 pages, 14 figures without appendices); accepted by MNRAS (a few extra references added)

Galaxies with monstrous black holes in galaxy cluster environments

Authors: Lieke A. C. van Son, Christopher Barber, Yannick M. Bahe, Joop Schaye, David J. Barnes, Robert A. Crain, Scott T. Kay, Tom Theuns, Claudio Dalla Vecchia

Abstract: Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations… ▽ More Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations, designed to enable the study of high-mass galaxy formation and evolution in cluster environments. We find 69 $M_{\rm BH}(M_{\star})$ outliers at $z=0$, defined as those with $ \rm M_{BH} >10^{7} M_{\odot}$ and $\rm M_{BH}/\rm M_{\star}> 0.01$. This paper focusses on a sample of 5 extreme outliers, that have been selected based on their $\rm M_{BH}$ and $\rm M_{\star}$ values, which are comparable to the most recent estimates of observed positive outliers. This sample of 5 outliers, classified as `Black hole monster galaxies' (BMGs), was traced back in time to study their origin and evolution. In agreement with the results of previous simulations for lower-mass $\rm M_{BH}(\rm M_{\star})$ outliers, we find that these galaxies became outliers due to a combination of their early formation times and tidal stripping. For BMGs with $\rm M_{BH} > 10^9 M_{\odot}$, major mergers (with a stellar mass ratio of $μ> 0.25$) at early times ($z>2$) precede the rapid growth of their supermassive BHs. Furthermore, the scatter in the relation between $\rm M_{BH}$ and stellar velocity dispersion, $σ$, correlates positively with the scatter in [Mg/Fe]($σ$). This indicates that the alpha enhancement of these galaxies, which is closely related to their star formation history, is related to the growth of their central BHs. △ Less

Submitted 10 January, 2019; originally announced January 2019.

Journal ref: Monthly Notices of the Royal Astronomical Society 485 (2019) 396 - 407

Radiative AGN feedback on a moving mesh: the impact of the galactic disc and dust physics on outflow properties

Authors: David J. Barnes, Rahul Kannan, Mark Vogelsberger, Federico Marinacci

Abstract: Feedback from accreting supermassive black holes, active galactic nuclei (AGN), is now a cornerstone of galaxy formation models. In this work, we present radiation-hydrodynamic simulations of radiative AGN feedback using the novel Arepo-RT. A central black hole emits radiation at a constant luminosity and drives an outflow via radiation pressure on dust grains. Utilising an isolated NFW halo we va… ▽ More Feedback from accreting supermassive black holes, active galactic nuclei (AGN), is now a cornerstone of galaxy formation models. In this work, we present radiation-hydrodynamic simulations of radiative AGN feedback using the novel Arepo-RT. A central black hole emits radiation at a constant luminosity and drives an outflow via radiation pressure on dust grains. Utilising an isolated NFW halo we validate our setup in the single and multi-scattering regimes, with the simulated shock front propagation in excellent agreement with the expected analytic result. For a spherically symmetric NFW halo, an examination of the simulated outflow properties generated by radiative feedback demonstrates that they are lower than typically observed at a fixed AGN luminosity, regardless of the collimation of the radiation. We then explore the impact of a central disc galaxy and the assumed dust model on the outflow properties. The contraction of the halo during the galaxy's formation and modelling the production of dust grains results in a factor $100$ increase in the halo's optical depth. Radiation is then able to couple momentum more efficiently to the gas, driving a stronger shock and producing a mass-loaded $\sim10^{3}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$ outflow with a velocity of $\sim2000\,\mathrm{km}\,\mathrm{s}^{-1}$, in agreement with observations. However, the inclusion of dust destruction mechanisms, like thermal sputtering, leads to the rapid destruction of dust grains within the outflow, reducing its properties below typically observed values. We conclude that radiative AGN feedback can drive outflows, but a thorough numerical and physical treatment is required to assess its true impact. △ Less

Submitted 4 December, 2018; originally announced December 2018.

Comments: 21 pages, 19 figures, submitted to MNRAS

SPMHD simulations of Structure Formation

Authors: David J. Barnes, Alvina Y. L. On, Kinwah Wu, Daisuke Kawata

Abstract: The intracluster medium of galaxy clusters is permeated by μG magnetic fields. Observations with current and future facilities have the potential to illuminate the role of these magnetic fields play in the astrophysical processes of galaxy clusters. To obtain a greater understanding of how the initial seed fields evolve to the magnetic fields in the intracluster medium requires magnetohydrodynamic… ▽ More The intracluster medium of galaxy clusters is permeated by μG magnetic fields. Observations with current and future facilities have the potential to illuminate the role of these magnetic fields play in the astrophysical processes of galaxy clusters. To obtain a greater understanding of how the initial seed fields evolve to the magnetic fields in the intracluster medium requires magnetohydrodynamic simulations. We critically assess the current Smoothed Particle Magneto-Hydrodynamics (SPMHD) schemes, especially highlighting the impact of a hyperbolic divergence cleaning scheme and artificial resistivity switch on the magnetic field evolution in cosmological simulations of the formation of a galaxy cluster using the N-body/SPMHD code gcmhd++. The impact and performance of the cleaning scheme and two different schemes for the artificial resistivity switch is demonstrated via idealized test cases and cosmological simulations. We demonstrate that the hyperbolic divergence cleaning scheme is effective at suppressing the growth of the numerical divergence error of the magnetic field and should be applied to any SPMHD simulation. Although the artificial resistivity is important in the strong field regime, it can suppress the growth of the magnetic field in the weak field regime, such as galaxy clusters. With sufficient resolution, simulations with divergence cleaning can reproduce observed magnetic fields. We conclude that the cleaning scheme alone is sufficient for galaxy cluster simulations, but our results indicate that the SPMHD scheme must be carefully chosen depending on the regime of the magnetic field. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: 15 pages, 11 figures, published (MNRAS 476 2890)

Early-type galaxy density profiles from IllustrisTNG: I. Galaxy correlations and the impact of baryons

Authors: Yunchong Wang, Mark Vogelsberger, Dandan Xu, Shude Mao, Volker Springel, Hui Li, David Barnes, Lars Hernquist, Annalisa Pillepich, Federico Marinacci, Rüediger Pakmor, Rainer Weinberger, Paul Torrey

Abstract: We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langleγ^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $σ_{γ^{\prime}} = 0.171$ over the radial ran… ▽ More We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langleγ^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $σ_{γ^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $γ^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $γ^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $γ^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $γ^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons. △ Less

Submitted 9 January, 2020; v1 submitted 15 November, 2018; originally announced November 2018.

Comments: 31 pages, 20 figures, 9 tables. Accepted for publication in MNRAS. Major revision, added the effects of feedback on the total density profile in Section 4. Pay attention to changes in Figures 7, 10, 12, and 16

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 491, Issue 4, February 2020, Pages 5188-5215

An application of machine learning techniques to galaxy cluster mass estimation using the MACSIS simulations

Authors: Thomas J. Armitage, Scott T. Kay, David J. Barnes

Abstract: Machine learning (ML) techniques, in particular supervised regression algorithms, are a promising new way to use multiple observables to predict a cluster's mass or other key features. To investigate this approach we use the \textsc{MACSIS} sample of simulated hydrodynamical galaxy clusters to train a variety of ML models, mimicking different datasets. We find that compared to predicting the clust… ▽ More Machine learning (ML) techniques, in particular supervised regression algorithms, are a promising new way to use multiple observables to predict a cluster's mass or other key features. To investigate this approach we use the \textsc{MACSIS} sample of simulated hydrodynamical galaxy clusters to train a variety of ML models, mimicking different datasets. We find that compared to predicting the cluster mass from the $σ-M$ relation, the scatter in the predicted-to-true mass ratio is reduced by a factor of 4, from $0.130\pm0.004$ dex (${\simeq} 35$ per cent) to $0.031 \pm 0.001$ dex (${\simeq} 7$ per cent) when using the same, interloper contaminated, spectroscopic galaxy sample. Interestingly, omitting line-of-sight galaxy velocities from the training set has no effect on the scatter when the galaxies are taken from within $r_{200c}$. We also train ML models to reproduce estimated masses derived from mock X-ray and weak lensing analyses. While the weak lensing masses can be recovered with a similar scatter to that when training on the true mass, the hydrostatic mass suffers from significantly higher scatter of ${\simeq} 0.13$ dex (${\simeq} 35$ per cent). Training models using dark matter only simulations does not significantly increase the scatter in predicted cluster mass compared to training on simulated clusters with hydrodynamics. In summary, we find ML techniques to offer a powerful method to predict masses for large samples of clusters, a vital requirement for cosmological analysis with future surveys. △ Less

Submitted 19 October, 2018; originally announced October 2018.

Comments: Submitted to MNRAS, 12 pages, 7 figures, example pipeline can be found here: https://github.com/TomArmitage/ML_Template

The signal of decaying dark matter with hydrodynamical simulations

Authors: Mark R. Lovell, David Barnes, Yannick Bahé, Joop Schaye, Matthieu Schaller, Tom Theuns, Sownak Bose, Robert A. Crain, Claudio dalla Vecchia, Carlos S. Frenk, Wojciech Hellwing, Scott T. Kay, Aaron D. Ludlow, Richard G. Bower

Abstract: Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate tha… ▽ More Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 20 pages, 16 Figures. Highlights: Figs 4, 13, 15. To be submitted to MNRAS, comments welcome. Contact: lovell@hi.is

The Cluster-EAGLE project: a comparison of dynamical mass estimators using simulated clusters

Authors: Thomas J. Armitage, Scott T. Kay, David J. Barnes, Yannick M. Bahé, Claudio Dalla Vecchia

Abstract: Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass,… ▽ More Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass, $M_{200c} = 10^{14.7} \, \mathrm{M_\odot}$. We quantify the bias and scatter of the Jeans, virial, and caustic mass estimators using all galaxies with a stellar mass $M_*> 10^9 \, \mathrm{M_\odot}$, both in the ideal 3D case and in the more realistic projected case. On average we find our mass estimates are unbiased, though relative to the true mass within $r_{200c}$ the scatter is large with a range of $0.09$ - $0.15$ dex. We see a slight increase in the scatter when projecting the clusters. Selecting galaxies using the same criteria, we find no significant difference in the mass bias or scatter when comparing results from hydrodynamical and dark matter only simulations. However, selecting galaxies by stellar mass reduces the bias compared to selecting by total mass. Comparing X-ray derived hydrostatic and dynamical masses, the former are ${\sim} 30$ per cent lower. We find a slight dependence between substructure, measured using two different metrics, and mass bias. In conclusion, we find that dynamical mass estimators, when averaged together, are unbiased with a scatter of $0.11 \pm 0.02$ dex when including interloper galaxies and with no prior knowledge of $r_{200c}$. △ Less

Submitted 25 October, 2018; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: Accepted in MNRAS, 14 figures in main text

The Deformable Mirror Demonstration Mission (DeMi) CubeSat: optomechanical design validation and laboratory calibration

Authors: Gregory Allan, Ewan S. Douglas, Derek Barnes, Mark Egan, Gabor Furesz, Warren Grunwald, Jennifer Gubner, Christian Haughwout, Bobby G. Holden, Paula do Vale Pereira, Abigail J. Stein, Kerri L. Cahoy

Abstract: Coronagraphs on future space telescopes will require precise wavefront correction to detect Earth-like exoplanets near their host stars. High-actuator count microelectromechanical system (MEMS) deformable mirrors provide wavefront control with low size, weight, and power. The Deformable Mirror Demonstration Mission (DeMi) payload will demonstrate a 140 actuator MEMS deformable mirror (DM) with \SI… ▽ More Coronagraphs on future space telescopes will require precise wavefront correction to detect Earth-like exoplanets near their host stars. High-actuator count microelectromechanical system (MEMS) deformable mirrors provide wavefront control with low size, weight, and power. The Deformable Mirror Demonstration Mission (DeMi) payload will demonstrate a 140 actuator MEMS deformable mirror (DM) with \SI{5.5}{\micro\meter} maximum stroke. We present the flight optomechanical design, lab tests of the flight wavefront sensor and wavefront reconstructor, and simulations of closed-loop control of wavefront aberrations. We also present the compact flight DM controller, capable of driving up to 192 actuator channels at 0-250V with 14-bit resolution. Two embedded Raspberry Pi 3 compute modules are used for task management and wavefront reconstruction. The spacecraft is a 6U CubeSat (30 cm x 20 cm x 10 cm) and launch is planned for 2019. △ Less

Submitted 10 July, 2018; v1 submitted 7 July, 2018; originally announced July 2018.

Comments: 15 pages, 10 figues. Presented at SPIE Astronomical Telescopes + Instrumentation, Austin, Texas, USA

Enhancing AGN efficiency and cool-core formation with anisotropic thermal conduction

Authors: David J. Barnes, Rahul Kannan, Mark Vogelsberger, Christoph Pfrommer, Ewald Puchwein, Rainer Weinberger, Volker Springel, Rüdiger Pakmor, Dylan Nelson, Federico Marinacci, Annalisa Pillepich, Paul Torrey, Lars Hernquist

Abstract: Understanding how baryonic processes shape the intracluster medium (ICM) is of critical importance to the next generation of galaxy cluster surveys. However, most models of structure formation neglect potentially important physical processes, like anisotropic thermal conduction (ATC). In this letter, we explore the impact of ATC on the prevalence of cool-cores (CCs) using 12 pairs of magnetohydrod… ▽ More Understanding how baryonic processes shape the intracluster medium (ICM) is of critical importance to the next generation of galaxy cluster surveys. However, most models of structure formation neglect potentially important physical processes, like anisotropic thermal conduction (ATC). In this letter, we explore the impact of ATC on the prevalence of cool-cores (CCs) using 12 pairs of magnetohydrodynamical galaxy cluster simulations, simulated using the IllustrisTNG model with and without ATC. Although the impact of ATC varies from cluster to cluster and with CC criterion, its inclusion produces a systematic shift to larger CC fractions at z = 0 for all CC criteria considered. Additionally, the inclusion of ATC yields a flatter CC fraction redshift evolution, easing the tension with the observed evolution. With ATC included, the energy required for the central black hole to achieve self-regulation is reduced and the gas fraction in the cluster core increases, resulting in larger CC fractions. ATC makes the ICM unstable to perturbations and the increased efficiency of AGN feedback suggests that its inclusion results in a greater level of mixing in the ICM. Therefore, ATC is potentially an important physical process in reproducing the thermal structure of the ICM. △ Less

Submitted 14 May, 2018; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 5 pages, 4 figures, submitted to MNRAS, comments welcome

Immersive Virtual Reality Experiences for All-Sky Data

Authors: C. J. Fluke, D. G. Barnes

Abstract: Spherical coordinate systems, which are ubiquitous in astronomy, cannot be shown without distortion on flat, two-dimensional surfaces. This poses challenges for the two complementary phases of visual exploration -- making discoveries in data by looking for relationships, patterns or anomalies -- and publication -- where the results of an exploration are made available for scientific scrutiny or co… ▽ More Spherical coordinate systems, which are ubiquitous in astronomy, cannot be shown without distortion on flat, two-dimensional surfaces. This poses challenges for the two complementary phases of visual exploration -- making discoveries in data by looking for relationships, patterns or anomalies -- and publication -- where the results of an exploration are made available for scientific scrutiny or communication. This is a long-standing problem, and many practical solutions have been developed. Our allskyVR approach provides a workflow for experimentation with commodity virtual reality head-mounted displays. Using the free, open source S2PLOT programming library, and the A-Frame WebVR browser-based framework, we provide a straightforward way to visualise all-sky catalogues on a user-centred, virtual celestial sphere. The allskyVR distribution contains both a quickstart option, complete with a gaze-based menu system, and a fully customisable mode for those who need more control of the immersive experience. The software is available for download from: https://github.com/cfluke/allskyVR △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 11 pages, 3 figures, accepted for publication in Publications of the Astronomical Society of Australia

CMEs in the Heliosphere: I. A Statistical Analysis of the Observational Properties of CMEs Detected in the Heliosphere from 2007 to 2017 by STEREO/HI-1

Authors: R. A. Harrison, J. A. Davies, D. Barnes, J. P. Byrne, C. H. Perry, V. Bothmer, J. P. Eastwood, P. T. Gallagher, E. K. J. Kilpua, C. Möstl, L. Rodriguez, A. P. Rouillard, D. Odstrcil

Abstract: We present a statistical analysis of coronal mass ejections (CMEs) imaged by the Heliospheric Imager (HI) instruments aboard NASAs twin-spacecraft STEREO mission between April 2007 and August 2017 for STEREO-A and between April 2007 and September 2014 for STEREO-B. The analysis exploits a catalogue that was generated within the FP7 HELCATS project. Here, we focus on the observational characteristi… ▽ More We present a statistical analysis of coronal mass ejections (CMEs) imaged by the Heliospheric Imager (HI) instruments aboard NASAs twin-spacecraft STEREO mission between April 2007 and August 2017 for STEREO-A and between April 2007 and September 2014 for STEREO-B. The analysis exploits a catalogue that was generated within the FP7 HELCATS project. Here, we focus on the observational characteristics of CMEs imaged in the heliosphere by the inner (HI-1) cameras. More specifically, in this paper we present distributions of the basic observational parameters - namely occurrence frequency, central position angle (PA) and PA span - derived from nearly 2000 detections of CMEs in the heliosphere by HI-1 on STEREO-A or STEREO-B from the minimum between Solar Cycles 23 and 24 to the maximum of Cycle 24; STEREO-A analysis includes a further 158 CME detections from the descending phase of Cycle 24, by which time communication with STEREO-B had been lost. We compare heliospheric CME characteristics with properties of CMEs observed at coronal altitudes, and with sunspot number. As expected, heliospheric CME rates correlate with sunspot number, and are not inconsistent with coronal rates once instrumental factors/differences in cataloguing philosophy are considered. As well as being more abundant, heliospheric CMEs, like their coronal counterparts, tend to be wider during solar maximum. Our results confirm previous coronagraph analyses suggesting that CME launch sites don't simply migrate to higher latitudes with increasing solar activity. At solar minimum, CMEs tend to be launched from equatorial latitudes while, at maximum, CMEs appear to be launched over a much wider latitude range; this has implications for understanding the CME/solar source association. Our analysis provides some supporting evidence for the systematic dragging of CMEs to lower latitude as they propagate outwards. △ Less

Submitted 9 April, 2018; v1 submitted 6 April, 2018; originally announced April 2018.

Comments: 26 pages, 8 figures, in press at Solar Physics. The abstract has been shortened to comply with the allowed length at arxiv.org

The connection between mass, environment and slow rotation in simulated galaxies

Authors: Claudia del P. Lagos, Joop Schaye, Yannick Bahe, Jesse van de Sande, Scott Kay, David Barnes, Timothy Davis, Claudio Dalla Vecchia

Abstract: Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of… ▽ More Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of volumes up to $(100\,\rm Mpc)^3$, while Hydrangea consists of $24$ cosmological simulations of galaxy clusters and their environment. Together they provide a statistically significant sample in the stellar mass range $10^{9.5}\,\rm M_{\odot}-10^{12.3}\,\rm M_{\odot}$, of $16,358$ galaxies. We construct IFS-like cubes and measure stellar spin parameters, $λ_{\rm R}$, and ellipticities, allowing us to classify galaxies into slow/fast rotators as in observations. The simulations display a primary dependence of $F_{\rm SR}$ on stellar mass, with a weak dependence on environment. At fixed stellar mass, satellite galaxies are more likely to be SRs than centrals. $F_{\rm SR}$ shows a dependence on halo mass at fixed stellar mass for central galaxies, while no such trend is seen for satellites. We find that $\approx 70$% of SRs at $z=0$ have experienced at least one merger with mass ratio $\ge 0.1$, with dry mergers being at least twice more common than wet mergers. Individual dry mergers tend to decrease $λ_{\rm R}$, while wet mergers mostly increase it. However, $30$% of SRs at $z=0$ have not experienced mergers, and those inhabit halos with median spins twice smaller than the halos hosting the rest of the SRs. Thus, although the formation paths of SRs can be varied, dry mergers and/or halos with small spins dominate. △ Less

Submitted 15 February, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: Accepted for publications in MNRAS (20 pages, 17 main body, 1.5 appendix). Changes include analysis of the orbital angular momentum effect on lambdaR and slightly different ellipticity calculation

The diverse density profiles of galaxy clusters with self-interacting dark matter plus baryons

Authors: Andrew Robertson, Richard Massey, Vincent Eke, Sean Tulin, Hai-Bo Yu, Yannick Bahé, David J. Barnes, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Scott T. Kay, Matthieu Schaller, Joop Schaye

Abstract: We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. D… ▽ More We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions. △ Less

Submitted 22 February, 2018; v1 submitted 24 November, 2017; originally announced November 2017.

Comments: 6 pages, 1 figure, v2 matches accepted MNRAS: Letters version

Localized massive halo properties in Bahamas and Macsis simulations: scalings, log-normality, and covariance

Authors: Arya Farahi, August E. Evrard, Ian McCarthy, David J. Barnes, Scott T. Kay

Abstract: Using tens of thousands of halos realized in the BAHAMAS and MACSIS simulations produced with a consistent astrophysics treatment that includes AGN feedback, we validate a multi-property statistical model for the stellar and hot gas mass behavior in halos hosting groups and clusters of galaxies. The large sample size allows us to extract fine-scale mass--property relations (MPRs) by performing loc… ▽ More Using tens of thousands of halos realized in the BAHAMAS and MACSIS simulations produced with a consistent astrophysics treatment that includes AGN feedback, we validate a multi-property statistical model for the stellar and hot gas mass behavior in halos hosting groups and clusters of galaxies. The large sample size allows us to extract fine-scale mass--property relations (MPRs) by performing local linear regression (LLR) on individual halo stellar mass (${\rm M}_{\rm star}$) and hot gas mass (${\rm M}_{\rm gas}$) as a function of total halo mass (${\rm M}_{\rm halo}$). We find that: 1) both the local slope and variance of the MPRs run with mass (primarily) and redshift (secondarily); 2) the conditional likelihood, $p({\rm M}_{\rm star},\ {\rm M}_{\rm gas} | \ {\rm M}_{\rm halo}, z)$ is accurately described by a multivariate, log-normal distribution, and; 3) the covariance of ${\rm M}_{\rm star}$ and ${\rm M}_{\rm gas}$ at fixed ${\rm M}_{\rm halo}$ is generally negative, reflecting a partially closed baryon box model for high mass halos. We validate the analytical population model of Evrard et al. (2014), finding sub-percent accuracy in the log-mean halo mass selected at fixed property, $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm gas} \rangle$ or $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm star} \rangle$, when scale-dependent MPR parameters are employed. This work highlights the potential importance of allowing for running in the slope and scatter of MPRs when modeling cluster counts for cosmological studies. We tabulate LLR fit parameters as a function of halo mass at $z=0$, $0.5$ and 1 for two popular mass conventions. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Comments: 17 pages, 7 figures, submitted to MNRAS, comments and suggestions are welcome

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

A census of cool core galaxy clusters in IllustrisTNG

Authors: David J. Barnes, Mark Vogelsberger, Rahul Kannan, Federico Marinacci, Rainer Weinberger, Volker Springel, Paul Torrey, Annalisa Pillepich, Dylan Nelson, Rüdiger Pakmor, Jill Naiman, Lars Hernquist, Michael McDonald

Abstract: The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are class… ▽ More The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are classified as CC and non-cool-core (NCC), and the criteria are increasingly correlated for more massive clusters. At $z=0$, the CC fractions for $2$ criteria are in reasonable agreement with the observed fractions but the other $4$ CC fractions are lower than observed. This result is partly driven by systematic differences between the simulated and observed gas fraction profiles. The simulated CC fractions with redshift show tentative agreement with the observed fractions, but linear fits demonstrate that the simulated evolution is steeper than observed. The conversion of CCs to NCCs appears to begin later and act more rapidly in the simulations. Examining the fraction of CCs and NCCs defined as relaxed we find no evidence that CCs are more relaxed, suggesting that mergers are not solely responsible for disrupting CCs. A comparison of the median thermodynamic profiles defined by different CC criteria shows that the extent to which they evolve in the cluster core is dependent on the CC criteria. We conclude that the thermodynamic structure of galaxy clusters in IllustrisTNG shares many similarities with observations, but achieving better agreement most likely requires modifications of the underlying galaxy formation model. △ Less

Submitted 21 November, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: 24 pages, 17 figures, MNRAS published version. The IllustrisTNG project website can be found at http://www.tng-project.org/

The Cluster-EAGLE project: velocity bias and the velocity dispersion - mass relation of cluster galaxies

Authors: Thomas Joshua Armitage, David. J. Barnes, Scott. T. Kay, Yannick M. Bahé, Claudio Dalla Vecchia, Robert A. Crain, Tom Theuns

Abstract: We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated s… ▽ More We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated sub-grid physics modelling and high numerical resolution (sub-kpc gravitational softening) making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5-10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased ($<5$ per cent) estimator of the velocity dispersion. This result holds out to $z=1.5$ and is relatively insensitive to the choice of cluster aperture, varying by less than 5 per cent between $r_{\rm 500c}$ and $r_{\rm 200m}$. We show that the velocity bias is a function of the time spent by a galaxy inside the cluster environment. Selecting galaxies by their total mass results in a larger bias because a larger fraction of objects have only recently entered the cluster and these have a velocity bias above unity. Galaxies that entered more than $4 \, \mathrm{Gyr}$ ago become progressively colder with time, as expected from dynamical friction. We conclude that velocity bias should not be a major issue when estimating cluster masses from kinematic methods. △ Less

Submitted 21 November, 2017; v1 submitted 1 August, 2017; originally announced August 2017.

Comments: 15 pages, 10 figures, accepted for publication in MNRAS

Real-time colouring and filtering with graphics shaders

Authors: Dany Vohl, Christopher J. Fluke, David G. Barnes, Amr H. Hassan

Abstract: Despite the popularity of the Graphics Processing Unit (GPU) for general purpose computing, one should not forget about the practicality of the GPU for fast scientific visualisation. As astronomers have increasing access to three dimensional (3D) data from instruments and facilities like integral field units and radio interferometers, visualisation techniques such as volume rendering offer means t… ▽ More Despite the popularity of the Graphics Processing Unit (GPU) for general purpose computing, one should not forget about the practicality of the GPU for fast scientific visualisation. As astronomers have increasing access to three dimensional (3D) data from instruments and facilities like integral field units and radio interferometers, visualisation techniques such as volume rendering offer means to quickly explore spectral cubes as a whole. As most 3D visualisation techniques have been developed in fields of research like medical imaging and fluid dynamics, many transfer functions are not optimal for astronomical data. We demonstrate how transfer functions and graphics shaders can be exploited to provide new astronomy-specific explorative colouring methods. We present 12 shaders, including four novel transfer functions specifically designed to produce intuitive and informative 3D visualisations of spectral cube data. We compare their utility to classic colour mapping. The remaining shaders highlight how common computation like filtering, smoothing and line ratio algorithms can be integrated as part of the graphics pipeline. We discuss how this can be achieved by utilising the parallelism of modern GPUs along with a shading language, letting astronomers apply these new techniques at interactive frame rates. All shaders investigated in this work are included in the open source software shwirl (Vohl, 2017). △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: 24 pages, 19 figures, 14 algorithms, 1 table. Accepted for publication by Monthly Notices of the Royal Astronomical Society

Reducing biases on $H_0$ measurements using strong lensing and galaxy dynamics: results from the EAGLE simulation

Authors: Amitpal S. Tagore, David J. Barnes, Neal Jackson, Scott T. Kay, Matthieu Schaller, Joop Schaye, Tom Theuns

Abstract: Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the r… ▽ More Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the radial density profiles of lenses, and by selecting lenses near isothermality, we find that the bias on H0 can be reduced to 5% with an intrinsic scatter of 10%, confirming previous results performed on a different simulation data set. We further investigate whether combining lensing observables with kinematic constraints helps to minimize this bias. We do not detect any significant dependence of the bias on lens model parameters or observational properties of the galaxy, but depending on the source--lens configuration, a bias may still exist. Cross lenses provide an accurate estimate of the Hubble constant, while fold (double) lenses tend to be biased low (high). With kinematic constraints, double lenses show bias and intrinsic scatter of 6% and 10%, respectively, while quad lenses show bias and intrinsic scatter of 0.5% and 10%, respectively. For lenses with a reduced $χ^2 > 1$, a power-law dependence of the $χ^2$ on the lens environment (number of nearby galaxies) is seen. Lastly, we model, in greater detail, the cases of two double lenses that are significantly biased. We are able to remove the bias, suggesting that the remaining biases could also be reduced by carefully taking into account additional sources of systematic uncertainty. △ Less

Submitted 23 June, 2017; originally announced June 2017.

Comments: 21 pages, 13 figures, 9 tables; revised version submitted to MNRAS