Showing 1–50 of 166 results for author: Babul, A

The HYENAS project: a prediction for the X-ray undetected galaxy groups

Authors: Weiguang Cui, Fred Jennings, Romeel Dave, Arif Babul, Ghassem Gozaliasl

Abstract: Galaxy groups contain the majority of bound mass with a significant portion of baryons due to the combination of halo mass and abundance (Cui 2024). Hence they serve as a crucial missing piece in the puzzle of galaxy formation and the evolution of large-scale structures in the Universe. In observations, mass-complete group catalogues are normally derived from galaxy redshift surveys detected throu… ▽ More Galaxy groups contain the majority of bound mass with a significant portion of baryons due to the combination of halo mass and abundance (Cui 2024). Hence they serve as a crucial missing piece in the puzzle of galaxy formation and the evolution of large-scale structures in the Universe. In observations, mass-complete group catalogues are normally derived from galaxy redshift surveys detected through various three-dimensional group-finding algorithms. Confirming the reality of such groups, particularly in the X-rays, is critical for ensuring robust studies of galaxy evolution in these environments. Recent works have reported numerous optical groups that are X-ray undetected (see, e.g., Popesso et al. 2024), sparking debates regarding the reasons for the unexpectedly low hot gas fraction in galaxy groups. To address this issue, we utilise zoomed-in simulations of galaxy groups from the novel HYENAS project to explore the range of hot gas fractions within galaxy groups and investigate the intrinsic factors behind the observed variability in X-ray emission. We find that the halo formation time can play a critical role -- we see that groups in halos that formed earlier exhibit up to an order of magnitude brighter X-ray luminosities compared to those formed later. This suggests that undetected X-ray groups are preferentially late-formed halos and highlights the connection between gas fraction and halo formation time in galaxy groups. Accounting for these biases in galaxy group identification is essential for advancing our understanding of galaxy formation and achieving precision in cosmological studies. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 10 pages, 3 figures. Comments are welcome

\texttt{Simba}-\texttt{C}: the evolution of the thermal and chemical properties in the intragroup medium

Authors: Renier T. Hough, Zhiwei Shao, Weiguang Cui, S. Ilani Loubser, Arif Babul, Romeel Davé, Douglas Rennehan, Chiaki Kobayashi

Abstract: The newly updated \texttt{GIZMO} and \texttt{Simba} based simulation, \texttt{Simba-C}, with its new stellar feedback, chemical enrichment, and recalibrated AGN feedback, allows for a detailed study of the intragroup medium X-ray properties. We discuss the impact of various physical mechanisms, e.g. stellar and AGN feedback, and chemical enrichment, on the composition and the global scaling relati… ▽ More The newly updated \texttt{GIZMO} and \texttt{Simba} based simulation, \texttt{Simba-C}, with its new stellar feedback, chemical enrichment, and recalibrated AGN feedback, allows for a detailed study of the intragroup medium X-ray properties. We discuss the impact of various physical mechanisms, e.g. stellar and AGN feedback, and chemical enrichment, on the composition and the global scaling relations of nearby galaxy groups. We also study the evolution ($z=2$ to $0$) of the global properties for the $1\,\mathrm{keV}$ temperature groups. \texttt{Simba-C} shows improved consistent matching with the observations of all X-ray scaling relations compared to \texttt{Simba}. It is well known that AGN feedback has a significant influence on $L_{X,0.5-2.0}-T_{spec,corr}$, $S_{500/2500}-T_{spec,corr}$, and gas mass fractions, with our \texttt{Simba-C} results consistent with it. Our recalibrated AGN feedback strength also showed an additional improvement in gas entropy, which now aligns with CLoGS observations. The updated stellar feedback and chemical enrichment model is shown to play an important role in our understanding of the chemical abundance ratios and their evolution within galaxy groups. In particular, we find that \texttt{Simba-C} produces an increase in the amount of heavier elements (specifically Si and Fe) relative to O, compared to \texttt{Simba}. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 20 pages, 13 figures, 2 tables, accepted by MNRAS on 6 June 2024

A hot core in the group-dominant elliptical galaxy NGC 777

Authors: Ewan O'Sullivan, Kamlesh Rajpurohit, Gerrit Schellenberger, Jan Vrtilek, Laurence P. David, Arif Babul, Valeria Olivares, Francesco Ubertosi, Konstantinos Kolokythas, Iurii Babyk, Ilani Loubser

Abstract: NGC 777 provides an example of a phenomenon observed in some group-central ellipticals, in which the temperature profile shows a central peak, despite the short central cooling time of the intra-group medium. We use deep Chandra X-ray observations of the galaxy, supported by uGMRT 400 MHz radio imaging, to investigate the origin of this hot core. We confirm the centrally-peaked temperature profile… ▽ More NGC 777 provides an example of a phenomenon observed in some group-central ellipticals, in which the temperature profile shows a central peak, despite the short central cooling time of the intra-group medium. We use deep Chandra X-ray observations of the galaxy, supported by uGMRT 400 MHz radio imaging, to investigate the origin of this hot core. We confirm the centrally-peaked temperature profile and find that entropy and cooling time both monotonically decline to low values (2.62 [+0.19, -0.18] keV cm$^2$ and 71.3 [+12.8, -13.1] Myr) in the central ~700 pc. Faint diffuse radio emission surrounds the nuclear point source, with no clear jets or lobes but extending to ~10 kpc on a northwest-southeast axis. This alignment and extent agree well with a previously identified filamentary H$α$+[NII] nebula. While cavities are not firmly detected, we see X-ray surface brightness decrements on the same axis at 10-20 kpc radius which are consistent with the intra-group medium having been pushed aside by expanding radio lobes. Any such outburst must have occurred long enough ago for lobe emission to have faded below detectability. Cavities on this scale would be capable of balancing radiative cooling for at least ~240 Myr. We consider possible causes of the centrally peaked temperature profile, including gravitational heating of gas as the halo relaxes after a period of AGN jet activity, and heating by particles leaking from the remnant relativistic plasma of the old radio jets. △ Less

Submitted 17 June, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: 19 pages, accepted for publication in the Astrophysical Journal, v2 corrects some minor typographical errors identified at the proof stage

Imprints of Supermassive Black Hole Evolution on the Spectral and Spatial Anisotropy of Nano-Hertz Stochastic Gravitational-Wave Background

Authors: Mohit Raj Sah, Suvodip Mukherjee, Vida Saeedzadeh, Arif Babul, Michael Tremmel, Thomas R. Quinn

Abstract: The formation and evolution of supermassive black holes (SMBHs) remains an open question in the field of modern cosmology. The detection of nanohertz (n-Hz) gravitational waves via pulsar timing arrays (PTAs) in the form of individual events and the stochastic gravitational wave background (SGWB) offers a promising avenue for studying SMBH formation across cosmic time, with SGWB signal being the i… ▽ More The formation and evolution of supermassive black holes (SMBHs) remains an open question in the field of modern cosmology. The detection of nanohertz (n-Hz) gravitational waves via pulsar timing arrays (PTAs) in the form of individual events and the stochastic gravitational wave background (SGWB) offers a promising avenue for studying SMBH formation across cosmic time, with SGWB signal being the immediately detectable signal with the currently accessible telescope sensitivities. By connecting the galaxy properties in the large scale (Gpc scale) cosmological simulation such as \texttt{MICECAT} with the small scale ($\sim$ Mpc scale) galaxy simulations from \texttt{ROMULUS}, we show that different scenarios of galaxy--SMBHs evolution with redshift leads to a frequency-dependent spatial anisotropy in the SGWB signal. The presence of slow evolution of the SMBHs in the Universe leads to a pronounced blue anisotropic spectrum of the SGWB. In contrast, if SMBHs grow faster in the Universe in lighter galaxies, the frequency-dependent spatial anisotropy exhibits a more flattened anisotropic spectrum. This additional aspect of the SGWB signal on top of the monopole SGWB signal, can give insight on how the SMBHs form in the high redshift Universe and its interplay with the galaxy formation from future measurements. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 19 page, 27 figures

Dual AGNs: Precursors of Binary Supermassive Black Hole Formation and Mergers

Authors: Vida Saeedzadeh, Arif Babul, Suvodip Mukherjee, Michael Tremmel, Thomas R. Quinn, Lucio Mayer

Abstract: The presence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. In this study, we use the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs (… ▽ More The presence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. In this study, we use the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs ($L_{bol} > 10^{43} \rm erg/s$) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at $z \leq 2$. This is our underlying population. Subsequently, we applied the luminosity threshold of $L_{bol} > 10^{43} erg/s$ to the neighboring SMBHs thereby identifying dual and multiple AGNs in our simulation. We examined the properties and statistics of dual AGNs in comparison to single AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. All dual AGNs in our sample resulted from major mergers. Compared with the single AGN population, duals are characterized by a lower black hole to halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate (SFR), and gas mass, are generally consistent with those of single AGN halos, albeit tending towards the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective. △ Less

Submitted 25 March, 2024; originally announced March 2024.

The Obsidian model: Three regimes of black hole feedback

Authors: Douglas Rennehan, Arif Babul, Belaid Moa, Romeel Davé

Abstract: In theoretical models of galaxy evolution, black hole feedback is a necessary ingredient in order to explain the observed exponential decline in number density of massive galaxies. Most contemporary black hole feedback models in cosmological simulations rely on a constant radiative efficiency (usually $η\sim 0.1$) at all black hole accretion rates. We present the Obsidian sub-grid model, a synthes… ▽ More In theoretical models of galaxy evolution, black hole feedback is a necessary ingredient in order to explain the observed exponential decline in number density of massive galaxies. Most contemporary black hole feedback models in cosmological simulations rely on a constant radiative efficiency (usually $η\sim 0.1$) at all black hole accretion rates. We present the Obsidian sub-grid model, a synthesis model for the spin-dependent radiative efficiencies of three physical accretion rate regimes, i.e. $η= η(j, \dot{M}_\mathrm{acc})$, for use in large-volume cosmological simulations. The three regimes include: an advection dominated accretion flow ($\dot{M}_\mathrm{acc} < 0.03\,\dot{M}_\mathrm{Edd}$), a quasar-like mode ($0.03 < \dot{M}_\mathrm{acc} / \dot{M}_\mathrm{Edd} < 0.3$), and a slim disc mode ($\dot{M}_\mathrm{acc} > 0.3\,\dot{M}_\mathrm{Edd}$). Additionally, we include a large-scale powerful jet at low accretion rates. The black hole feedback model we present is a kinetic model that prescribes mass loadings but could be used in thermal models directly using the radiative efficiency. We implement the Obsidian model into the Simba galaxy evolution model to determine if it is possible to reproduce galaxy populations successfully, and provide a first calibration for further study. Using a $2\times1024^3$ particle cosmological simulation in a $(150\,\mathrm{cMpc})^3$ volume, we found that the model is successful in reproducing the galaxy stellar mass function, black hole mass-stellar mass relationship, and stellar mass-halo mass relationship. Moving forward, this model opens new avenues for exploration of the impact of black hole feedback on galactic environments. △ Less

Submitted 3 July, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 16 pages, 10 figures, comments more than welcome

Shining Light on the Hosts of the Nano-Hertz Gravitational Wave Sources: A Theoretical Perspective

Authors: Vida Saeedzadeh, Suvodip Mukherjee, Arif Babul, Michael Tremmel, Thomas R. Quinn

Abstract: The formation of supermassive black holes (SMBHs) in the Universe and its role in the properties of the galaxies is one of the open questions in astrophysics and cosmology. Though, traditionally, electromagnetic waves have been instrumental in direct measurements of SMBHs, significantly influencing our comprehension of galaxy formation, gravitational waves (GW) bring an independent avenue to detec… ▽ More The formation of supermassive black holes (SMBHs) in the Universe and its role in the properties of the galaxies is one of the open questions in astrophysics and cosmology. Though, traditionally, electromagnetic waves have been instrumental in direct measurements of SMBHs, significantly influencing our comprehension of galaxy formation, gravitational waves (GW) bring an independent avenue to detect numerous binary SMBHs in the observable Universe in the nano-Hertz range using the pulsar timing array observation. This brings a new way to understand the connection between the formation of binary SMBHs and galaxy formation if we can connect theoretical models with multi-messenger observations namely GW data and galaxy surveys. Along these lines, we present here the first paper on this series based on {\sc Romulus25} cosmological simulation on the properties of the host galaxies of SMBHs and propose on how this can be used to connect with observations of nano-Hertz GW signal and galaxy surveys. We show that the most dominant contribution to the background will arise from sources with high chirp masses which are likely to reside in low redshift early-type galaxies with high stellar mass, largely old stellar population, and low star formation rate, and that reside at centers of galaxy groups and manifest evidence of recent mergers. The masses of the sources show a correlation with the halo mass and stellar mass of the host galaxies. This theoretical study will help in understanding the host properties of the GW sources and can help in establishing a connection with observations. △ Less

Submitted 15 February, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 17 pages, 12 figures. Accepted for publication in MNRAS

A large population of strongly lensed faint submillimetre galaxies in future dark energy surveys inferred from JWST imaging

Authors: James Pearson, Stephen Serjeant, Wei-Hao Wang, Zhen-Kai Gao, Arif Babul, Scott Chapman, Chian-Chou Chen, David L. Clements, Christopher J. Conselice, James Dunlop, Lulu Fan, Luis C. Ho, Ho Seong Hwang, Maciej Koprowski, Michał Michałowski, Hyunjin Shim

Abstract: Bright galaxies at sub-millimetre wavelengths from Herschel are now well known to be predominantly strongly gravitationally lensed. The same models that successfully predicted this strongly lensed population also predict about one percent of faint $450μ$m-selected galaxies from deep James Clerk Maxwell Telescope (JCMT) surveys will also be strongly lensed. Follow-up ALMA campaigns have so far foun… ▽ More Bright galaxies at sub-millimetre wavelengths from Herschel are now well known to be predominantly strongly gravitationally lensed. The same models that successfully predicted this strongly lensed population also predict about one percent of faint $450μ$m-selected galaxies from deep James Clerk Maxwell Telescope (JCMT) surveys will also be strongly lensed. Follow-up ALMA campaigns have so far found one potential lens candidate, but without clear compelling evidence e.g. from lensing arcs. Here we report the discovery of a compelling gravitational lens system confirming the lensing population predictions, with a $z_{s} = 3.4 {\pm} 0.4$ submm source lensed by a $z_{spec} = 0.360$ foreground galaxy within the COSMOS field, identified through public JWST imaging of a $450μ$m source in the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) catalogue. These systems will typically be well within the detectable range of future wide-field surveys such as Euclid and Roman, and since sub-millimetre galaxies are predominantly very red at optical/near-infrared wavelengths, they will tend to appear in near-infrared channels only. Extrapolating to the Euclid-Wide survey, we predict tens of thousands of strongly lensed near-infrared galaxies. This will be transformative for the study of dusty star-forming galaxies at cosmic noon, but will be a contaminant population in searches for strongly lensed ultra-high-redshift galaxies in Euclid and Roman. △ Less

Submitted 9 January, 2024; v1 submitted 2 September, 2023; originally announced September 2023.

Comments: 9 pages, 5 figures, 4 tables; added figure, updated Discussion section, and minor revisions; accepted for publication in MNRAS

SIMBA-C: An updated chemical enrichment model for galactic chemical evolution in the SIMBA simulation

Authors: Renier T. Hough, Douglas Rennehan, Chiaki Kobayashi, S. Ilani Loubser, Romeel Davé, Arif Babul, Weiguang Cui

Abstract: We introduce a new chemical enrichment and stellar feedback model into GIZMO, using the SIMBA sub-grid models as a base. Based on the state-of-the-art chemical evolution model of Kobayashi et al., SIMBA-C tracks 34 elements from H$\rightarrow$Ge and removes SIMBA's instantaneous recycling approximation. Furthermore, we make some minor improvements to SIMBA's base feedback models. SIMBA-C provides… ▽ More We introduce a new chemical enrichment and stellar feedback model into GIZMO, using the SIMBA sub-grid models as a base. Based on the state-of-the-art chemical evolution model of Kobayashi et al., SIMBA-C tracks 34 elements from H$\rightarrow$Ge and removes SIMBA's instantaneous recycling approximation. Furthermore, we make some minor improvements to SIMBA's base feedback models. SIMBA-C provides significant improvements on key diagnostics such as the knee of the $z=0$ galaxy stellar mass function, the faint end of the main sequence, and the ability to track black holes in dwarf galaxies. SIMBA-C also matches better with recent observations of the mass-metallicity relation at $z=0,2$. By not assuming instantaneous recycling, SIMBA-C provides a much better match to galactic abundance ratio measures such as [O/Fe] and [N/O]. SIMBA-C thus opens up new avenues to constrain feedback models using detailed chemical abundance measures across cosmic time. △ Less

Submitted 12 August, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 16 Pages, 9 Figure, 1 Table, Accepted for MNRAS publication: 8 August 2023, doi:10.1093/mnras/stad2394

Cool and gusty, with a chance of rain: Dynamics of multiphase CGM around massive galaxies in the Romulus simulations

Authors: Vida Saeedzadeh, S. Lyla Jung, Douglas Rennehan, Arif Babul, Michael Tremmel, Thomas R. Quinn, Zhiwei Shao, Prateek Sharma, Lucio Mayer, E. OSullivan, S. Ilani Loubser

Abstract: Using high-resolution {\sc Romulus} simulations, we explore the origin and evolution of the circumgalactic medium (CGM) in the region 0.1 $\leq \mathrm{R}/\mathrm{R}_\mathrm{500} \leq$ 1 around massive central galaxies in group-scale halos. We find that the CGM is multiphase and highly dynamic. Investigating the dynamics, we identify seven patterns of evolution. We show that these are robust and d… ▽ More Using high-resolution {\sc Romulus} simulations, we explore the origin and evolution of the circumgalactic medium (CGM) in the region 0.1 $\leq \mathrm{R}/\mathrm{R}_\mathrm{500} \leq$ 1 around massive central galaxies in group-scale halos. We find that the CGM is multiphase and highly dynamic. Investigating the dynamics, we identify seven patterns of evolution. We show that these are robust and detected consistently across various conditions. The gas cools via two pathways: (1) filamentary cooling inflows and (2) condensations forming from rapidly cooling density perturbations. In our cosmological simulations, the perturbations are mainly seeded by orbiting substructures. The condensations can form even when the median $t_\mathrm{cool} / t_\mathrm{ff}$ of the X-ray emitting gas is above 10 or 20. Strong amplitude perturbations can provoke runaway cooling regardless of the state of the background gas. We also find perturbations whose local $t_\mathrm{cool} / t_\mathrm{ff}$ ratios drop below the threshold but which do not condense. Rather, the ratios fall to some minimum value and then bounce. These are weak perturbations that are temporarily swept up in satellite wakes and carried to larger radii. Their $t_\mathrm{cool} / t_\mathrm{ff}$ ratios decrease because $t_\mathrm{ff}$ is increasing, not because $t_\mathrm{cool}$ is decreasing. For structures forming hierarchically, our study highlights the challenge of using a simple threshold argument to infer the CGM's evolution. It also highlights that the median hot gas properties are suboptimal determinants of the CGM's state and dynamics. Realistic CGM models must incorporate the impact of mergers and orbiting satellites, along with the CGM's heating and cooling cycles. △ Less

Submitted 1 September, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

Comments: 25 pages, 12 figures

Brightest Cluster Galaxy Formation in the z=4.3 Protocluster SPT2349-56: Discovery of a Radio-Loud AGN

Authors: Scott C. Chapman, Ryley Hill, Manuel Aravena, Melanie Archipley, Arif Babul, James Burgoyne, Rebecca E. A. Canning, Carlos De Breuck, Anthony H. Gonzalez, Christopher C. Hayward, Seon Woo Kim, Matt Malkan, Dan P. Marrone, Vincent McIntyre, Eric Murphy, Emily Pass, Ryan W. Perry, Kedar A. Phadke, Douglas Rennehan, Cassie Reuter, Kaja M. Rotermund, Douglas Scott, Nick Seymour, Manuel Solimano, Justin Spilker , et al. (7 additional authors not shown)

Abstract: We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.… ▽ More We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.3, consistent with ATCA non-detections at 5.5 and 9GHz, and implying L_1.4(rest)=(2.4pm0.3)x10^{26}W/Hz. This radio luminosity is about 100 times higher than expected from star formation, assuming the usual FIR-radio correlation, which is a clear indication of an AGN driven by a forming brightest cluster galaxy (BCG). None of the SMGs in SPT2349-56 show signs of AGN in any other diagnostics available to us (notably 12CO out to J=16, OH163um, CII/IR, and optical spectra), highlighting the radio continuum as a powerful probe of obscured AGN in high-z protoclusters. No other significant radio detections are found amongst the cluster members, consistent with the FIR-radio correlation. We compare these results to field samples of radio sources and SMGs, along with the 22 SPT-SMG gravitational lenses also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. Our results allow us to better understand the effects of this gas-rich, overdense environment on early supermassive black hole (SMBH) growth and cluster feedback. We estimate that (3.3pm0.7)x10^{38} W of power are injected into the growing ICM by the radio-loud AGN, whose energy over 100Myr is comparable to the binding energy of the gas mass of the central halo. The AGN power is also comparable to the instantaneous energy injection from supernova feedback from the 23 catalogued SMGs in the core region of 120kpc projected radius. The SPT2349-56 radio-loud AGN may be providing strong feedback on a nascent ICM. △ Less

Submitted 4 January, 2023; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: 31 pages, submitted to ApJ, Dec17,2022

Spatially-resolved properties of early-type group-dominant galaxies with MUSE: gas content, ionisation mechanisms and metallicity gradients

Authors: P. Lagos, S. I. Loubser, T. C. Scott, E. O'Sullivan, K. Kolokythas, A. Babul, A. Nigoche-Netro, V. Olivares, C. Sengupta

Abstract: With the goal of a thorough investigation of the ionised gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). This data allowed us to study the spatially-resolved warm gas properties, including the morphology of the ionised gas, EW(H$α$) and kinematics as well as the gas-phase metallicity (12 +… ▽ More With the goal of a thorough investigation of the ionised gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). This data allowed us to study the spatially-resolved warm gas properties, including the morphology of the ionised gas, EW(H$α$) and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionisation mechanisms, we used the emission-line ratios [O III]/H$β$ and [N II]/H$α$ in the BPT diagrams and EW(H$α$). We find that the ionisation sources in our sample have variable impacts at different radii, central regions are more influenced by low-luminosity AGN, while extended regions of LINER-like emission are ionised by other mechanisms with pAGBs photoionisation likely contributing significantly. We classified our sample into three H$α$+[N II] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionisation sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionised gas structures. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS

Merger histories of brightest group galaxies from MUSE stellar kinematics

Authors: S. I. Loubser, P. Lagos, A. Babul, E. O'Sullivan, S. L. Jung, V. Olivares, K. Kolokythas

Abstract: Using Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy, we analyse the stellar kinematics of 18 brightest group early-type (BGEs) galaxies, selected from the Complete Local-Volume Groups Sample (CLoGS). We analyse the kinematic maps for distinct features, and measure specific stellar angular momentum within one effective radius ($λ_{e}$). We classify the BGEs as fast (10/18) or slow (8/18) ro… ▽ More Using Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy, we analyse the stellar kinematics of 18 brightest group early-type (BGEs) galaxies, selected from the Complete Local-Volume Groups Sample (CLoGS). We analyse the kinematic maps for distinct features, and measure specific stellar angular momentum within one effective radius ($λ_{e}$). We classify the BGEs as fast (10/18) or slow (8/18) rotators, suggesting at least two different evolution paths. We quantify the anti-correlation between higher-order kinematic moment $h_{3}$ and V/$σ$ (using the $ξ_{3}$ parameter), and the kinematic misalignment angle between the photometric and kinematic position angles (using the $Ψ$ parameter), and note clear differences between these parameter distributions of the fast and slow rotating BGEs. We find that all 10 of our fast rotators are aligned between the morphological and kinematical axis, consistent with an oblate galaxy shape, whereas the slow rotators are spread over all three classes: oblate (1/8), triaxial (4/8), and prolate (3/8). We place the results into context using known radio properties, X-ray properties, and observations of molecular gas. We find consistent merger histories inferred from observations for the fast-rotating BGEs, indicating that they experienced gas-rich mergers or interactions, and these are very likely the origin of the cold gas. Observational evidence for the slow rotators are consistent with gas-poor mergers. For the slow rotators with cold gas, all evidence point to cold gas cooling from the intragroup medium. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS

Massive central galaxies of galaxy groups in the Romulus simulations: an overview of galaxy properties at z=0

Authors: Seoyoung Lyla Jung, Douglas Rennehan, Vida Saeedzadeh, Arif Babul, Michael Tremmel, Thomas R. Quinn, S. Ilani Loubser, E. O'Sullivan, Sukyoung K. Yi

Abstract: Contrary to many stereotypes about massive galaxies, observed brightest group galaxies (BGGs) are diverse in their star formation rates, kinematic properties, and morphologies. Studying how they evolve into and express such diverse characteristics is an important piece of the galaxy formation puzzle. We use a high-resolution cosmological suite of simulations Romulus and compare simulated central g… ▽ More Contrary to many stereotypes about massive galaxies, observed brightest group galaxies (BGGs) are diverse in their star formation rates, kinematic properties, and morphologies. Studying how they evolve into and express such diverse characteristics is an important piece of the galaxy formation puzzle. We use a high-resolution cosmological suite of simulations Romulus and compare simulated central galaxies in group-scale halos at $z=0$ to observed BGGs. The comparison encompasses the stellar mass-halo mass relation, various kinematic properties and scaling relations, morphologies, and the star formation rates. Generally, we find that Romulus reproduces the full spectrum of diversity in the properties of the BGGs very well, albeit with a tendency toward lower than the observed fraction of quenched BGGs. We find both early-type S0 and elliptical galaxies as well as late-type disk galaxies; we find Romulus galaxies that are fast-rotators as well as slow-rotators; and we observe galaxies transforming from late-type to early-type following strong dynamical interactions with satellites. We also carry out case studies of selected Romulus galaxies to explore the link between their properties, and the recent evolution of the stellar system as well as the surrounding intragroup/circumgalactic medium. In general, mergers/strong interactions quench star-forming activity and disrupt the stellar disk structure. Sometimes, however, such interactions can also trigger star-formation and galaxy rejuvenation. Black hole feedback can also lead to a decline of the star formation rate but by itself, it does not typically lead to complete quenching of the star formation activity in the BGGs. △ Less

Submitted 10 June, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 27 pages, 18 figures, Accepted to MNRAS

Gas condensation in Brightest Group Galaxies unveiled with MUSE

Authors: V. Olivares, P. Salome, S. L. Hamer, F. Combes, M. Gaspari, K. Kolokythas, E. O'Sullivan, R. S. Beckmann, A. Babul, F. L. Polles, M. Lehnert, S. I. Loubser, M. Donahue, M. -L. Gendron-Marsolais, P. Lagos, G. Pineau des Forets, B. Godard, T. Rose, G. Tremblay, G. Ferland, P. Guillard

Abstract: The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extendi… ▽ More The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extending up to 10 kpc to two compact (<3 kpc) and five extended (>5 kpc) disk-dominated structures. Some rotating disks show rings and elongated structures arising from the central disk. The kinematics of the stellar component is mainly rotation-dominated, which is very different from the disturbed kinematics and distribution found in the filamentary sources. The ionized gas is kinematically decoupled from the stellar component for most systems, suggesting an external origin for the gas. We also find that the Halpha luminosity correlates with the cold molecular mass. By exploring the thermodynamical properties of the hot atmospheres, we find that the filamentary sources and compact disks are found in systems with small central entropy values and tcool/teddy ratios. This suggests that, like for Brightest Cluster Galaxies in cool core clusters, the ionized gas are likely formed from hot halo gas condensations, consistently with the Chaotic Cold Accretion simulations (as shown via the C-ratio, Tat, and k-plot). We note that gaseous rotating disks are more frequent than in BCGs. An explanation for the origin of the gas in those objects is a contribution to gas fueling by mergers or group satellites, as qualitatively hinted by some sources of the present sample. Nonetheless, we discuss the possibility that some extended disks could also be a transition stage in an evolutionary sequence including filaments, extended disks and compact disks, as described by hot gas condensation models of cooling flows. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 38 pages, and 16 figures. Resubmitted to A&A after referee report

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

The Complete Local Volume Groups Sample -- IV. Star formation and gas content in group-dominant galaxies

Authors: Konstantinos Kolokythas, Sravani Vaddi, Ewan O'Sullivan, Ilani Loubser, Arif Babul, Somak Raychaudhury, Patricio Lagos, Thomas H. Jarrett

Abstract: Using multi-band data we examine the star formation activity of the nearby group-dominant early-type galaxies of the Complete Local-volume Groups Sample (CLoGS), and the relation between star formation, gas content, and local environment. Only a small fraction of the galaxies (13%; 6/47) are found to be Far-Ultraviolet (FUV) bright, with FUV to near-infrared colours indicative of recent active sta… ▽ More Using multi-band data we examine the star formation activity of the nearby group-dominant early-type galaxies of the Complete Local-volume Groups Sample (CLoGS), and the relation between star formation, gas content, and local environment. Only a small fraction of the galaxies (13%; 6/47) are found to be Far-Ultraviolet (FUV) bright, with FUV to near-infrared colours indicative of recent active star formation (NGC 252, NGC 924, NGC 940, NGC 1106, NGC 7252, and ESO 507-25). These systems are lenticulars presenting the highest FUV specific star-formation rates in the sample (sSFR FUV > 5$\times$10$^{13}$ yr$^{-1}$), significant cold gas reservoirs (M(H2)=0.5 - 61$\times$10$^8$ M$_\odot$), reside in X-ray faint groups, and none of them hosts a powerful radio AGN (P$_{1.4GHz}$<$10^{23}$ W Hz$^{-1}$). The majority of the group-dominant galaxies (87%; 41/47) are FUV faint, with no significant star formation, classified in most cases as spheroids based on their position on the infrared star-forming main sequence (87%; 46/53). Examining the relationships between radio power, SFR FUV, and stellar mass we find a lack of correlation that suggests a combination of origins for the cool gas in these galaxies, including stellar mass loss, cooling from the intra-group medium (IGrM) or galaxy halo, and acquisition through mergers or tidal interactions. X-ray bright systems, in addition to hosting radio powerful AGN, have a range of SFRs but, with the exception of NGC 315, do not rise to the highest rates seen in the FUV bright systems. We suggest that central group galaxy evolution is linked to gas mass availability, with star formation favoured in the absence of a group-scale X-ray halo, but AGN jet launching is more likely in systems with a cooling IGrM. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS, 17 pages with 13 figures and 2 tables

Molecular gas along the old radio jets of the cluster-central type 2 quasar IRAS 09104+4109

Authors: Ewan O'Sullivan, Françoise Combes, Arif Babul, Scott Chapman, Kedar A. Phadke, Gerrit Schellenberger, Philippe Salomé

Abstract: We present Northern Extended Millimeter Array (NOEMA) CO(2-1) maps of the z=0.4418 cluster-central QSO IRAS 09104+4109, which trace ~4.5x10^10 Msol of molecular gas in and around the galaxy. As in many low redshift cool core clusters, the molecular gas is located in a series of clumps extending along the old radio jets and lobes. It has a relatively low velocity dispersion (336 [+39,-35] km/s FWHM… ▽ More We present Northern Extended Millimeter Array (NOEMA) CO(2-1) maps of the z=0.4418 cluster-central QSO IRAS 09104+4109, which trace ~4.5x10^10 Msol of molecular gas in and around the galaxy. As in many low redshift cool core clusters, the molecular gas is located in a series of clumps extending along the old radio jets and lobes. It has a relatively low velocity dispersion (336 [+39,-35] km/s FWHM) and shows no velocity gradients indicative of outflow or infall. Roughly half the gas is located in a central clump on the northeast side of the galaxy, overlapping a bright ionized gas filament and a spur of excess X-ray emission, suggesting that this is a location of rapid cooling. The molecular gas is unusually extended, out to ~55 kpc radius, comparable to the scale of the filamentary nebula in the Perseus cluster, despite the much higher redshift of this system. The extent falls within the thermal instability radius of the intracluster medium (ICM), with t_cool/t_ff<25 and t_cool}/t_eddy~1 within ~70 kpc. Continuum measurements at 159.9 GHz from NOEMA and 850 micron from the JCMT SCUBA-2 show excess far infrared emission, which we interpret as free-free emission arising from the ongoing starburst. These observations suggest that ICM cooling is not strongly affected by the buried QSO, and that cooling from the ICM can build gas reservoirs sufficient to fuel quasar-mode activity and drive the reorientation of the central AGN. △ Less

Submitted 14 October, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 16 pages, 13 figures, 2 tables, accepted for publication in MNRAS, v2 updates references and corrects minor typographical errors identified at the proof stage

The first nova eruption in a novalike variable: YZ Ret as seen in X-rays and gamma-rays

Authors: Kirill V. Sokolovsky, Kwan-Lok Li, Raimundo Lopes de Oliveira, Jan-Uwe Ness, Koji Mukai, Laura Chomiuk, Elias Aydi, Elad Steinberg, Indrek Vurm, Brian D. Metzger, Aliya-Nur Babul, Adam Kawash, Justin D. Linford, Thomas Nelson, Kim L. Page, Michael P. Rupen, Jennifer L. Sokoloski, Jay Strader, David Kilkenny

Abstract: Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT… ▽ More Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT and NuSTAR observations complemented by XMM-Newton X-ray grating spectroscopy to probe the physical parameters of the shocked ejecta and the nova-hosting white dwarf. The XMM-Newton observations revealed a super-soft X-ray emission which is dominated by emission lines of CV, CVI, NVI, NVII, and OVIII rather than a blackbody-like continuum, suggesting CO-composition of the white dwarf in a high-inclination binary system. Fermi/LAT detected YZ Ret for 15 days with the gamma-ray spectrum best described by a power law with an exponential cut-off at 1.9 +/-0.6 GeV. In stark contrast with theoretical predictions and in keeping with previous NuSTAR observations of Fermi-detected classical novae (V5855 Sgr and V906 Car), the 3.5-78 keV X-ray emission is found to be two orders of magnitude fainter than the GeV emission. The X-ray emission observed by NuSTAR is consistent with a single-temperature thermal plasma. We detect no non-thermal tail of the GeV emission expected to extend down to the NuSTAR band. NuSTAR observations continue to challenge theories of high-energy emission from shocks in novae. △ Less

Submitted 21 May, 2022; v1 submitted 6 August, 2021; originally announced August 2021.

Comments: 22 pages, 5 figures, 4 tables; accepted to MNRAS

Shocks and dust formation in nova V809 Cep

Authors: Aliya-Nur Babul, Jennifer L. Sokoloski, Laura Chomiuk, Justin D. Linford, Jennifer H. S. Weston, Elias Aydi, Kirill V. Sokolovsky, Adam M. Kawash

Abstract: The discovery that many classical novae produce detectable GeV $γ$-ray emission has raised the question of the role of shocks in nova eruptions. Here we use radio observations of nova V809 Cep (Nova Cep 2013) with the Jansky Very Large Array to show that it produced non-thermal emission indicative of particle acceleration in strong shocks for more than a month starting about six weeks into the eru… ▽ More The discovery that many classical novae produce detectable GeV $γ$-ray emission has raised the question of the role of shocks in nova eruptions. Here we use radio observations of nova V809 Cep (Nova Cep 2013) with the Jansky Very Large Array to show that it produced non-thermal emission indicative of particle acceleration in strong shocks for more than a month starting about six weeks into the eruption, quasi-simultaneous with the production of dust. Broadly speaking, the radio emission at late times -- more than a six months or so into the eruption -- is consistent with thermal emission from $10^{-4} M_\odot$ of freely expanding, $10^4$~K ejecta. At 4.6 and 7.4 GHz, however, the radio light-curves display an initial early-time peak 76 days after the discovery of the eruption in the optical ($t_0$). The brightness temperature at 4.6 GHz on day 76 was greater than $10^5 K$, an order of magnitude above what is expected for thermal emission. We argue that the brightness temperature is the result of synchrotron emission due to internal shocks within the ejecta. The evolution of the radio spectrum was consistent with synchrotron emission that peaked at high frequencies before low frequencies, suggesting that the synchrotron from the shock was initially subject to free-free absorption by optically thick ionized material in front of the shock. Dust formation began around day 37, and we suggest that internal shocks in the ejecta were established prior to dust formation and caused the nucleation of dust. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Simulating Groups and the IntraGroup Medium: The Surprisingly Complex and Rich Middle Ground Between Clusters and Galaxies

Authors: Benjamin D. Oppenheimer, Arif Babul, Yannick Bahé, Iryna S. Butsky, Ian G. McCarthy

Abstract: Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories.… ▽ More Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories. Widely used cosmological simulations with ~(100 Mpc)^3 volumes contain statistical samples of groups that provide important tests of galaxy evolution influenced by environmental processes. Larger volumes capable of reproducing LSS while following the redistribution of baryons by cooling and feedback are essential tools necessary to constrain cosmological parameters. Higher resolution simulations can currently model satellite interactions, the processing of cool (T~10^4 K) multi-phase gas, and non-thermal physics including turbulence, magnetic fields, and cosmic ray transport. We review simulation results regarding the gas and stellar contents of groups, cooling flows and the relation to the central galaxy, the formation and processing of multi-phase gas, satellite interactions with the intragroup medium, and the impact of groups for cosmological parameter estimation. Cosmological simulations provide evolutionarily consistent predictions of these observationally difficult-to-define objects, and have untapped potential to accurately model their gaseous, stellar, and dark matter distributions. △ Less

Submitted 24 June, 2021; originally announced June 2021.

Comments: 61 pages, 18 figures, and 2 tables. This review article is part of the special issue "The Physical Properties of the Groups of Galaxies", edited by L. Lovisari and S. Ettori. Published in MDPI - Universe: https://www.mdpi.com/journal/universe/special_issues/PPGG

The Dawn of Disk Formation in a Milky Way-sized Galaxy Halo: Thin Stellar Disks at $z > 4$

Authors: Tomas Tamfal, Lucio Mayer, Thomas R. Quinn, Arif Babul, Piero Madau, Pedro R. Capelo, Sijing Shen, Marius Staub

Abstract: We present results from \textsc{GigaEris}, a cosmological, $N$-body hydrodynamical "zoom-in" simulation of the formation of a Milky Way-sized galaxy halo with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. W… ▽ More We present results from \textsc{GigaEris}, a cosmological, $N$-body hydrodynamical "zoom-in" simulation of the formation of a Milky Way-sized galaxy halo with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We focus on the early assembly of the galaxy, down to redshift $z=4.4$. The simulated galaxy has properties consistent with extrapolations of the main sequence of star-forming galaxies to higher redshifts and levels off to a star formation rate of $\sim$60$\, M_{\odot}$~yr$^{-1}$ at $z=4.4$. A compact, thin rotating stellar disk with properties analogous to those of low-redshift systems arises already at $z \sim 8$. The galaxy rapidly develops a multi-component structure, and the disk, at least at these early stages, does not grow "upside-down" as often reported in the literature. Rather, at any given time, newly born stars contribute to sustain a thin disk. The kinematics reflect the early, ubiquitous presence of a thin disk, as a stellar disk component with $v_φ/σ_R$ larger than unity is already present at $z \sim 9$--10. Our results suggest that high-resolution spectro-photometric observations of very high-redshift galaxies should find thin rotating disks, consistent with the recent discovery of cold rotating gas disks by ALMA. Finally, we present synthetic images for the JWST NIRCam camera, showing how the early disk would be easily detectable already at those early times. △ Less

Submitted 22 February, 2022; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: 19 pages, 14 figures, 1 table, accepted to ApJ

A population of heavily reddened, optically missed novae from Palomar Gattini-IR: Constraints on the Galactic nova rate

Authors: Kishalay De, Mansi M. Kasliwal, Matthew J. Hankins, Jennifer L. Sokoloski, Scott M. Adams, Michael C. B. Ashley, Aliya-Nur Babul, Ashot Bagdasaryan, Alexandre Delacroix, Richard Dekany, Timothee Greffe, David Hale, Jacob E. Jencson, Viraj R. Karambelkar, Ryan M. Lau, Ashish Mahabal, Daniel McKenna, Anna M. Moore, Eran O. Ofek, Manasi Sharma, Roger M. Smith, Jamie Soon, Roberto Soria, Gokul Srinivasaragavan, Samaporn Tinyanont , et al. (3 additional authors not shown)

Abstract: The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated in the range of $\approx10-300$ yr$^{-1}$, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galax… ▽ More The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated in the range of $\approx10-300$ yr$^{-1}$, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galaxies; both methods are subject to debatable assumptions. The total discovery rate of optical novae remains much smaller ($\approx5-10$ yr$^{-1}$) than these estimates, even with the advent of all-sky optical time domain surveys. Here, we present a systematic sample of 12 spectroscopically confirmed Galactic novae detected in the first 17 months of Palomar Gattini-IR (PGIR), a wide-field near-infrared time domain survey. Operating in $J$-band ($\approx1.2$ $μ$m) that is relatively immune to dust extinction, the extinction distribution of the PGIR sample is highly skewed to large extinction values ($> 50$% of events obscured by $A_V\gtrsim5$ mag). Using recent estimates for the distribution of mass and dust in the Galaxy, we show that the observed extinction distribution of the PGIR sample is commensurate with that expected from dust models. The PGIR extinction distribution is inconsistent with that reported in previous optical searches (null hypothesis probability $< 0.01$%), suggesting that a large population of highly obscured novae have been systematically missed in previous optical searches. We perform the first quantitative simulation of a $3π$ time domain survey to estimate the Galactic nova rate using PGIR, and derive a rate of $\approx 46.0^{+12.5}_{-12.4}$ yr$^{-1}$. Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population. △ Less

Submitted 24 January, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: 25 pages, 12 figures. Submitted to ApJ

Dynamical masses of brightest cluster galaxies II: constraints on the stellar IMF

Authors: S. I. Loubser, H. Hoekstra, A. Babul, Y. M. Bahé, M. Donahue

Abstract: We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30, to place constraints on their stellar Initial Mass Function (IMF). We measure the spatially-resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($Υ_{\star \rm POP}$). We find… ▽ More We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30, to place constraints on their stellar Initial Mass Function (IMF). We measure the spatially-resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($Υ_{\star \rm POP}$). We find young stellar populations ($<$200 Myr) in the centres of 22 per cent of the sample, and constant $Υ_{\star \rm POP}$ within 15 kpc for 60 per cent of the sample. We further use the stellar mass-to-light ratio from the dynamical mass profiles of the BCGs ($Υ_{\star \rm DYN}$), modelled using a Multi-Gaussian Expansion (MGE) and Jeans Anisotropic Method (JAM), with the dark matter contribution explicitly constrained from weak gravitational lensing measurements. We directly compare the stellar mass-to-light ratios derived from the two independent methods, $Υ_{\star \rm POP}$ (assuming some IMF) to $Υ_{\star \rm DYN}$ for the subsample of BCGs with no young stellar populations and constant $Υ_{\star \rm POP}$. We find that for the majority of these BCGs, a Salpeter (or even more bottom-heavy) IMF is needed to reconcile the stellar population and dynamical modelling results although for a small number of BCGs, a Kroupa (or even lighter) IMF is preferred. For those BCGs better fit with a Salpeter IMF, we find that the mass-excess factor against velocity dispersion falls on an extrapolation (towards higher masses) of known literature correlations. We conclude that there is substantial scatter in the IMF amongst the highest-mass galaxies. △ Less

Submitted 10 November, 2020; originally announced November 2020.

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

LoCuSS: The splashback radius of massive galaxy clusters and its dependence on cluster merger history

Authors: Matteo Bianconi, Riccardo Buscicchio, Graham P. Smith, Sean L. McGee, Chris P. Haines, Alexis Finoguenov, Arif Babul

Abstract: We present the direct detection of the splashback feature using the sample of massive galaxy clusters from the Local Cluster Substructure Survey (LoCuSS). This feature is clearly detected (above $5σ$) in the stacked luminosity density profile obtained using the K-band magnitudes of spectroscopically confirmed cluster members. We obtained the best-fit model by means of Bayesian inference, which ran… ▽ More We present the direct detection of the splashback feature using the sample of massive galaxy clusters from the Local Cluster Substructure Survey (LoCuSS). This feature is clearly detected (above $5σ$) in the stacked luminosity density profile obtained using the K-band magnitudes of spectroscopically confirmed cluster members. We obtained the best-fit model by means of Bayesian inference, which ranked models including the splashback feature as more descriptive of the data with respect to models that do not allow for this transition. In addition, we have assessed the impact of the cluster dynamical state on the occurrence of the splashback feature. We exploited the extensive multi-wavelength LoCuSS dataset to test a wide range of proxies for the cluster formation history, finding the most significant dependence of the splashback feature location and scale according to the presence or absence of X-ray emitting galaxy groups in the cluster infall regions. In particular, we report for the first time that clusters that do not show massive infalling groups present the splashback feature at a smaller clustercentric radius $ r_{\rm{sp}}/r_{\rm{200,m}} = 1.158 \pm 0.071$ than clusters that are actively accreting groups $r_{\rm{sp}}/r_{\rm{200,m}} = 1.291 \pm 0.062$. The difference between these two sub-samples is significant at $4.2σ$, suggesting a correlation between the properties of the cluster potential and its accretion rate and merger history. Similarly, clusters that are classified as old and dynamically inactive present stronger signatures of the splashback feature, with respect to younger, more active clusters. We are directly observing how fundamental dynamical properties of clusters reverberate across vastly different physical scales. △ Less

Submitted 8 March, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 18 pages, 8 figures, 4 tables. Re-submission to ApJ after the first round of comments. Minor text modifications to improve clarity, and updated reference list

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

Revisiting dynamical friction: the role of global modes and local wakes

Authors: Tomas Tamfal, Lucio Mayer, Thomas R. Quinn, Pedro R. Capelo, Stelios Kazantzidis, Arif Babul, Douglas Potter

Abstract: The orbital decay of a perturber within a larger system plays a key role in the dynamics of many astrophysical systems -- from nuclear star clusters or globular clusters in galaxies, to massive black holes in galactic nuclei, to dwarf galaxy satellites within the dark matter halos of more massive galaxies. For many decades, there have been various attempts to determine the underlying physics and t… ▽ More The orbital decay of a perturber within a larger system plays a key role in the dynamics of many astrophysical systems -- from nuclear star clusters or globular clusters in galaxies, to massive black holes in galactic nuclei, to dwarf galaxy satellites within the dark matter halos of more massive galaxies. For many decades, there have been various attempts to determine the underlying physics and time-scales of the drag mechanism, ranging from the local dynamical friction approach to descriptions based on the back-reaction of global modes induced in the background system. We present ultra-high-resolution $N$-body simulations of massive satellites orbiting a Milky Way-like galaxy (with $> 10^8$ particles), that appear to capture both the local "wake" and the global "mode" induced in the primary halo. We address directly the mechanism of orbital decay from the combined action of local and global perturbations and specifically analyze where the bulk of the torque originates. △ Less

Submitted 27 May, 2021; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 18 pages, 11 figures, 3 tables, accepted by ApJ

Optical and near-infrared observations of the SPT2349-56 proto-cluster core at z = 4.3

Authors: K. M. Rotermund, S. C. Chapman, K. A. Phadke, R. Hill, E. Pass, M. Aravena, M. L. N. Ashby, A. Babul, M. Béthermin, R. Canning, C. de Breuck, C. Dong, A. H. Gonzalez, C. C. Hayward, S. Jarugula, D. P. Marrone, D. Narayanan, C. Reuter, D. Scott, J. S. Spilker, J. D. Vieira, G. Wang, A. Weiss

Abstract: We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this pro… ▽ More We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation. △ Less

Submitted 12 January, 2021; v1 submitted 27 June, 2020; originally announced June 2020.

Dynamical masses of brightest cluster galaxies I: stellar velocity anisotropy and mass-to-light ratios

Authors: S. I. Loubser, A. Babul, H. Hoekstra, Y. M. Bahé, E. O'Sullivan, M. Donahue

Abstract: We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for… ▽ More We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for $h_{4}$ for all the BCGs. We derive the stellar mass-to-light ratio ($Υ_{\star \rm DYN}$), and velocity anisotropy ($β$) based on a Multi-Gaussian Expansion (MGE) and axisymmetric Jeans Anisotropic Methods (JAM, cylindrically- and spherically-aligned). We explicitly include a dark matter halo mass component, which is constrained by weak gravitational lensing measurements for these clusters. We find a strong correlation between anisotropy and velocity dispersion profile slope, with rising velocity dispersion profiles corresponding to tangential anisotropy and decreasing velocity dispersion profiles corresponding to radial anisotropy. The rising velocity dispersion profiles can also indicate a significant contribution from the intracluster light (ICL) to the total light (in projection) in the centre of the galaxy. For a small number of BCGs with rising velocity dispersion profiles, a variable stellar mass-to-light ratio can also account for the profile shape, instead of tangential anisotropy or a significant ICL contribution. We note that, for some BCGs, a variable $β_{z}(r)$ (from radial to tangential anisotropy) can improve the model fit to the observed kinematic profiles. The observed diversity in these properties illustrates that BCGs are not the homogeneous class of objects they are often assumed to be. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: 25 pages, accepted for publication in MNRAS

The synchrotron maser emission from relativistic magnetized shocks: Dependence on the pre-shock temperature

Authors: Aliya-Nur Babul, Lorenzo Sironi

Abstract: Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of Fast Radio Bursts. By means of two-dimensional particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron-positron perpendicular shocks as a function of the pre-shock magnet… ▽ More Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of Fast Radio Bursts. By means of two-dimensional particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron-positron perpendicular shocks as a function of the pre-shock magnetization $σ\gtrsim 1$ (i.e., the ratio of incoming Poynting flux to particle energy flux) and thermal spread $Δγ\equiv kT/mc^2=10^{-5}-10^{-1}$. We measure the fraction $f_ξ$ of total incoming energy that is converted into precursor waves, as computed in the post-shock frame. At fixed magnetization, we find that $f_ξ$ is nearly independent of temperature as long as $Δγ\lesssim 10^{-1.5}$ (with only a modest decrease of a factor of three from $Δγ=10^{-5}$ to $Δγ=10^{-1.5}$), but it drops by nearly two orders of magnitude for $Δγ\gtrsim 10^{-1}$. At fixed temperature, the scaling with magnetization $f_ξ\sim 10^{-3}\,σ^{-1}$ is consistent with our earlier one-dimensional results. For our reference $σ=1$, the power spectrum of precursor waves is relatively broad (fractional width $\sim 1-3$) for cold temperatures, whereas it shows pronounced line-like features with fractional width $\sim 0.2$ for $10^{-3} \lesssim Δγ\lesssim 10^{-1.5} $. For $σ\gtrsim 1$, the precursor waves are beamed within an angle $\simeq σ^{-1/2}$ from the shock normal (as measured in the post-shock frame), as required so they can outrun the shock. Our results can provide physically-grounded inputs for FRB emission models based on maser emission from relativistic shocks. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 12 pages, 15 figures, submitted to MNRAS

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) IV: Spatial clustering and halo masses of 450-$μ$m-selected sub-millimeter galaxies

Authors: Chen-Fatt Lim, Chian-Chou Chen, Ian Smail, Wei-Hao Wang, Wei-Leong Tee, Yen-Ting Lin, Douglas Scott, Yoshiki Toba, Yu-Yen Chang, YiPing Ao, Arif Babul, Andy Bunker, Scott C. Chapman, David L Clements, Christopher J. Conselice, Yu Gao, Thomas R. Greve, Luis C. Ho, Sungwook E. Hong, Ho Seong Hwang, Maciej Koprowski, Michał J. Michałowski, Hyunjin Shim, Xinwen Shu, James M. Simpson

Abstract: We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS… ▽ More We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$μ$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$μ$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$μ$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$μ$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$μ$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies. △ Less

Submitted 2 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: ApJ accepted

S2COSMOS: Evolution of Gas Mass with Redshift Using Dust Emission

Authors: Jenifer S. Millard, Stephen A. Eales, M. W. L. Smith, H. L. Gomez, K. Małek, J. M. Simpson, Y. Peng, M. Sawicki, R. A. Beeston, Andrew Bunker, Y. Ao, A. Babul, L. C. Ho, Ho Seong Hwang, M. J. Michałowski, N. Scoville, H. Shim, Y. Toba

Abstract: We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$μ$m. We use stacking methodologies on the 850$μ$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature st… ▽ More We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$μ$m. We use stacking methodologies on the 850$μ$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 17 pages (+ 9 pages of appendices), 8 figures (+ 7 figures in appendices), accepted for publication in MNRAS

Fountains and storms: The role of AGN and mergers in disrupting the cool-core in the RomulusC simulation

Authors: Urmila Chadayammuri, Michael Tremmel, Daisuke Nagai, Arif Babul, Thomas Quinn

Abstract: The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star fo… ▽ More The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star formation, especially in the Brightest Cluster Galaxy (BCG), likely because both rely on cold phase gas. A self-consistent model for AGN and star formation in galaxy clusters thus requires cosmological context, higher resolution, and a careful modeling of cooling and heating balance. In this paper, we use the high-resolution hydrodynamical cosmological simulation of the RomulusC galaxy cluster to study in detail the role of AGN and a major, head-on merger in shaping the cluster core. The unprecedented resolution of the RomulusC simulation captures the multiphase structure of the ICM. The realistic large-scale outflows launched by very small-scale thermal injections, the improved modeling of turbulent diffusion and mixing, and the particle nature of the simulation allow us to carefully separate different heating channels. We show that AGN activity, while efficient at regulating star formation, is incapable of destroying a CC. Instead, that process is facilitated by a head-on, 1:8 mass ratio merger. The merger generates bulk and turbulent motions, which in turn mix high entropy gas generated by AGN and merger driven shocks, turbulent dissipation and sloshing of the ICM by infalling substructures. While central cooling times remain shorter than the Hubble time, restoring a CC is made more difficult by the reduced precipitation rates at larger radii, emphasizing that the AGN-ICM connection is truly a multi-scale problem. △ Less

Submitted 7 April, 2021; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 16 pages, 9 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 504, Issue 3, July 2021

SCUBA-2 Ultra Deep Imaging Eao Survey (Studies) III: Multi-wavelength properties, luminosity functions and preliminary source catalog of 450-$μ$m-selected galaxies

Authors: Chen-Fatt Lim, Wei-Hao Wang, Ian Smail, Douglas Scott, Chian-Chou Chen, Yu-Yen Chang, James M. Simpson, Yoshiki Toba, Xinwen Shu, Dave Clements, Josh Greenslade, YiPing Ao, Arif Babul, Jack Birkin, Scott C. Chapman, Tai-An Cheng, Brian S. Cho, Helmut Dannerbauer, Ugnė Dudzevičiūtė, James Dunlop, Yu Gao, Tomotsugu Goto, Luis C. Ho, Li-Ting Hsu, Ho Seong Hwang , et al. (13 additional authors not shown)

Abstract: We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of… ▽ More We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$μ$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$β_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$μ$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$. △ Less

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

Comments: Published in the ApJ

Journal ref: ApJ 889 (2020) 80

Palomar Gattini-IR: Survey overview, data processing system, on-sky performance and first results

Authors: Kishalay De, Matthew J. Hankins, Mansi M. Kasliwal, Anna M. Moore, Eran O. Ofek, Scott M. Adams, Michael C. B. Ashley, Aliya-Nur Babul, Ashot Bagdasaryan, Kevin B. Burdge, Jill Burnham, Richard G. Dekany, Alexander Declacroix, Antony Galla, Tim Greffe, David Hale, Jacob E. Jencson, Ryan M. Lau, Ashish Mahabal, Daniel McKenna, Manasi Sharma, Patrick L. Shopbell, Roger M. Smith, Jamie Soon, Jennifer Sokoloski , et al. (2 additional authors not shown)

Abstract: (Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a p… ▽ More (Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$σ$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 35 pages, 29 figures. Submitted to PASP

Cool-Core Cycles and Phoenix

Authors: Deovrat Prasad, Prateek Sharma, Arif Babul, G. Mark Voit, Brian W. O'Shea

Abstract: Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM… ▽ More Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM) in a state of rough thermal balance. It has been argued that either AGN feedback saturates in very massive clusters or the central supermassive black hole (SMBH) is too small to produce enough kinetic feedback and hence is unable to quench the catastrophic cooling. In this work, we present an alternate scenario wherein intense short-lived cooling and star formation phases followed by strong AGN outbursts are part of the AGN feedback loop. Using results from a 3D hydrodynamic simulation of a standard cool-core cluster ($M_{200}\sim 7\times10^{14}$ M$_\odot$; $z=0$), scaled to account for differences in mass and redshift, we argue that {\it Phoenix} is at the end of a cooling phase in which an AGN outburst has begun but has not yet arrested core cooling. This state of high cooling rate and star formation is expected to last for $\lesssim$ 100 Myr in {\it Phoenix}. △ Less

Submitted 28 April, 2020; v1 submitted 27 September, 2019; originally announced September 2019.

Comments: 6 pages, 5 figures, Accepted for publication in MNRAS

The Complete Local Volume Groups Sample -- III. Characteristics of group central radio galaxies in the Local Universe

Authors: Konstantinos Kolokythas, Ewan O'Sullivan, Huib Intema, Somak Raychaudhury, Arif Babul, Simona Giacintucci, Myriam Gitti

Abstract: Using new 610 MHz and 235 MHz observations from the Giant Metrewave Radio Telescope (GMRT) in combination with archival GMRT and Very Large Array (VLA) survey data we present the radio properties of the dominant early-type galaxies in the low$-$richness sub-sample of the Complete Local-volume Groups Sample (CLoGS; 27 galaxy groups) and provide results for the radio properties of the full CLoGS sam… ▽ More Using new 610 MHz and 235 MHz observations from the Giant Metrewave Radio Telescope (GMRT) in combination with archival GMRT and Very Large Array (VLA) survey data we present the radio properties of the dominant early-type galaxies in the low$-$richness sub-sample of the Complete Local-volume Groups Sample (CLoGS; 27 galaxy groups) and provide results for the radio properties of the full CLoGS sample for the first time. We find a high radio detection rate in the dominant galaxies of the low-richness sub-sample of 82% (22/27); for the full CLoGS sample, the detection rate is 87% (46/53). The group-dominant galaxies exhibit a wide range of radio power, 10$^{20}$ $-$ 10$^{25}$ W Hz$^{-1}$ in the 235 and 610 MHz bands, with the majority (53%) presenting point-like radio emission, 19% hosting currently active radio jets, 6% having remnant jets, 9% being diffuse and 13% having no detected radio emission. The mean spectral index of the detected radio sources in the 235$-$610 MHz frequency range is found to be $α_{235}^{610}\sim$0.68, and $α_{235}^{1400}\sim$0.59 in the 235$-$1400 MHz one. In agreement with earlier studies, we find that the fraction of ultra-steep spectrum sources ($α>$1.3) is $\sim$4%, mostly dependent on the detection limit at 235 MHz. The majority of point-like systems are found to reside in dynamically young groups, whereas jet systems show no preference between spiral-rich and spiral-poor group environments. The mechanical power of the jet sources in the low$-$richness sample groups is estimated to be $\sim$10$^{42}$ $-$ 10$^{44}$ erg s$^{-1}$ with their black hole masses ranging between 2$\times$10$^{8}$ $-$ 5$\times$10$^{9}$ M$_{\odot}$. We confirm previous findings that, while radio jet sources tend to be associated with more massive black holes, black hole mass is not the decisive factor in determining jet activity or power. △ Less

Submitted 26 July, 2019; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS, 48 pages of which 18 Manuscript pages with 10 tables and 8 figures, plus 30 pages of appendices with 25 figures and 1 table

Rapid early coeval star formation and assembly of the most massive galaxies in the universe

Authors: Douglas Rennehan, Arif Babul, Christopher C. Hayward, Connor Bottrell, Maan H. Hani, Scott C. Chapman

Abstract: The current consensus on the formation and evolution of the brightest cluster galaxies is that their stellar mass forms early ($z \gtrsim 4$) in separate galaxies that then eventually assemble the main structure at late times ($z \lesssim 1$). However, advances in observational techniques have led to the discovery of protoclusters out to $z \sim 7$, suggesting that the late-assembly picture may no… ▽ More The current consensus on the formation and evolution of the brightest cluster galaxies is that their stellar mass forms early ($z \gtrsim 4$) in separate galaxies that then eventually assemble the main structure at late times ($z \lesssim 1$). However, advances in observational techniques have led to the discovery of protoclusters out to $z \sim 7$, suggesting that the late-assembly picture may not be fully complete. If these protoclusters assemble rapidly in the early universe, they should form the brightest cluster galaxies much earlier than suspected by the late-assembly picture. Using a combination of observationally constrained hydrodynamical and dark-matter-only simulations, we show that the stellar assembly time of a sub-set of the brightest cluster galaxies occurs at high redshifts ($z > 3$) rather than at low redshifts ($z < 1$), as is commonly thought. We find, using isolated non-cosmological hydrodynamical simulations, that highly overdense protoclusters assemble their stellar mass into brightest cluster galaxies within $\sim 1$ $\mathrm{Gyr}$ of evolution -- producing massive blue elliptical galaxies at high redshifts ($z \gtrsim 1.5$). We argue that there is a downsizing effect on the cluster scale wherein some of the brightest cluster galaxies in the cores of the most-massive clusters assemble earlier than those in lower-mass clusters. In those clusters with $z = 0$ virial mass $\geqslant 5\times10^{14}$ M$_\mathrm{\odot}$, we find that $9.8$% have their cores assembly early, and a higher fraction of $16.4$% in those clusters above $10^{15}$ M$_\mathrm{\odot}$. The James Webb Space Telescope will be able to detect and confirm our prediction in the near future, and we discuss the implications to constraining the value of $σ_\mathrm{8}$. △ Less

Submitted 4 February, 2020; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: 15 pages, 6 figures, revised and submitted to MNRAS

Circumgalactic Gas and the Precipitation Limit

Authors: G. M. Voit, A. Babul, Iu. Babyk, G. L. Bryan, H. -W. Chen, M. Donahue, D. Fielding, M. Gaspari, Y. Li, M. McDonald, B. W. O'Shea, D. Prasad, P. Sharma, M. Sun, G. Tremblay, J. Werk, N. Werner, F. Zahedy

Abstract: During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galax… ▽ More During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution. △ Less

Submitted 29 March, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Comments: Astro2020 White Paper, 8 pages, 2 figures (differences from Astro2020 version: some typos fixed, some references added)

Astro2020 Science White Paper: A Shocking Shift in Paradigm for Classical Novae

Authors: Laura Chomiuk, Elias Aydi, Aliya-Nur Babul, Andrea Derdzinski, Adam Kawash, Kwan-Lok Li, Justin Linford, Brian D. Metzger, Koji Mukai, Michael P. Rupen, Jennifer Sokoloski, Kirill Sokolovsky, Elad Steinberg

Abstract: The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate… ▽ More The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate their bolometric luminosity. Shocks may also explain long-standing mysteries in novae such as dust production, super-Eddington luminosities, and `flares' in optical light curves. Here, we highlight the multi-wavelength facilities of the next decade that will further test our nova shock model and fulfill the promise of novae as powerful astrophysical laboratories. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers, 7 pages, 2 figures

Mass bias evolution in tSZ cluster cosmology

Authors: Laura Salvati, Marian Douspis, Anna Ritz, Nabila Aghanim, Arif Babul

Abstract: Galaxy clusters observed through the thermal Sunyaev-Zeldovich (tSZ) effect are a recent cosmological probe. The precision on the cosmological constraints is affected mainly by the current knowledge of cluster physics, which enters the analysis through the scaling relations. Here we aim to study one of the most important sources of systematic uncertainties, the mass bias, $b$. We have analysed the… ▽ More Galaxy clusters observed through the thermal Sunyaev-Zeldovich (tSZ) effect are a recent cosmological probe. The precision on the cosmological constraints is affected mainly by the current knowledge of cluster physics, which enters the analysis through the scaling relations. Here we aim to study one of the most important sources of systematic uncertainties, the mass bias, $b$. We have analysed the effects of a mass-redshift dependence, adopting a power-law parametrisation. We applied this parametrisation to the combination of tSZ number counts and power spectrum, finding a hint of redshift dependence that leads to a decreasing value of the mass bias for higher redshift. We tested the robustness of our results for different mass bias calibrations and a discrete redshift dependence. We find our results to be dependent on the clusters sample that we are considering, in particular obtaining an inverse (decreasing) redshift dependence when neglecting $z<0.2$ clusters. We analysed the effects of this parametrisation on the combination of cosmic microwave background (CMB) primary anisotropies and tSZ galaxy clusters. We find a preferred constant value of mass bias, having $(1-b) =0.62 \pm 0.05$. The corresponding value of $b$ is too high with respect to weak lensing and numerical simulations estimations. Therefore we conclude that this mass-redshift parametrisation does not help in solving the remaining discrepancy between CMB and tSZ clusters observations. △ Less

Submitted 15 May, 2019; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 14 pages, 9 figures

Journal ref: A&A 626, A27 (2019)

The SCUBA-2 Web Survey: I. Observations of CO(3-2) in hyper-luminous QSO fields

Authors: Ryley Hill, Scott C. Chapman, Douglas Scott, Ian Smail, Charles C. Steidel, Melanie Krips, Arif Babul, Frank Bertoldi, Yu Gao, Kevin Lacaille, Yuichi Matsuda

Abstract: A primary goal of the SCUBA-2 Web survey is to perform tomography of the early inter-galactic medium by studying systems containing some of the brightest quasi-stellar objects (QSOs; 2.5<z<3.0) and nearby submillimetre galaxies. As a first step, this paper aims to characterize the galaxies that host the QSOs. To achieve this, a sample of 13 hyper-luminous (L_AGN>10^14 L_odot) QSOs with previous su… ▽ More A primary goal of the SCUBA-2 Web survey is to perform tomography of the early inter-galactic medium by studying systems containing some of the brightest quasi-stellar objects (QSOs; 2.5<z<3.0) and nearby submillimetre galaxies. As a first step, this paper aims to characterize the galaxies that host the QSOs. To achieve this, a sample of 13 hyper-luminous (L_AGN>10^14 L_odot) QSOs with previous submillimetre continuum detections were followed up with CO(3-2) observations using the NOEMA interferometer. All but two of the QSOs are detected in CO(3-2); for one non-detection, our observations show a tentative 2sigma line at the expected position and redshift, and for the other non-detection we find only continuum flux density an order of magnitude brighter than the other sources. In three of the fields, a companion potentially suitable for tomography is detected in CO line emission within 25 arcsec of the QSO. We derive gas masses, dynamical masses and far-infrared luminosities, and show that the QSOs in our sample have similar properties as compared to less luminous QSOs and SMGs in the literature, despite the fact that their black-hole masses (which are proportional to L_AGN) are 1-2 orders of magnitude larger. We discuss two interpretations of these observations: this is due to selection effects, such as preferential face-on viewing angles and picking out objects in the tail ends of the scatter in host-galaxy mass and black-hole mass relationships; or the black hole masses have been overestimated because the accretion rates are super-Eddington. △ Less

Submitted 18 March, 2019; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 16 pages, 7 figures. Latest version published in MNRAS

Journal ref: MNRAS 485, 753 (2019)

The large-scale structure of the halo of the Andromeda galaxy II. Hierarchical structure in the Pan-Andromeda Archaeological Survey

Authors: Alan W. McConnachie, Rodrigo Ibata, Nicolas Martin, Annette M. N. Ferguson, Michelle Collins, Stephen Gwyn, Mike Irwin, Geraint F. Lewis, A. Dougal Mackey, Tim Davidge, Veronica Arias, Anthony Conn, Patrick Cote, Denija Crnojevic, Avon Huxor, Jorge Penarrubia, Chelsea Spengler, Nial Tanvir, David Valls-Gabaud, Arif Babul, Pauline Barmby, Nicholas F. Bate, Edouard Bernard, Scott Chapman, Aaron Dotter , et al. (7 additional authors not shown)

Abstract: The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections betwee… ▽ More The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products. △ Less

Submitted 25 October, 2018; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: Accepted for publication in the Astrophysical Journal. 51 pages, 24 figures, 5 tables. Some figures have degraded resolution. All PAndAS data products are available via the CADC at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/community/pandas/query.html where you can also find a version of the paper with full resolution figures

Dynamic Localised Turbulent Diffusion and its Impact on the Galactic Ecosystem

Authors: Douglas Rennehan, Arif Babul, Philip F. Hopkins, Romeel Davé, Belaid Moa

Abstract: Modelling the turbulent diffusion of thermal energy, momentum, and metals is required in all galaxy evolution simulations due to the ubiquity of turbulence in galactic environments. The most commonly employed diffusion model, the Smagorinsky model, is known to be over-diffusive due to its strong dependence on the fluid velocity shear. We present a method for dynamically calculating a more accurate… ▽ More Modelling the turbulent diffusion of thermal energy, momentum, and metals is required in all galaxy evolution simulations due to the ubiquity of turbulence in galactic environments. The most commonly employed diffusion model, the Smagorinsky model, is known to be over-diffusive due to its strong dependence on the fluid velocity shear. We present a method for dynamically calculating a more accurate, locally appropriate, turbulent diffusivity: the dynamic localised Smagorinsky model. We investigate a set of standard astrophysically-relevant hydrodynamical tests, and demonstrate that the dynamic model curbs over-diffusion in non-turbulent shear flows and improves the density contrast in our driven turbulence experiments. In galactic discs, we find that the dynamic model maintains the stability of the disc by preventing excessive angular momentum transport, and increases the metal-mixing timescale in the interstellar medium. In both our isolated Milky Way-like galaxies and cosmological simulations, we find that the interstellar and circumgalactic media are particularly sensitive to the treatment of turbulent diffusion. We also examined the global gas enrichment fractions in our cosmological simulations, to gauge the potential effect on the formation sites and population statistics of Population III stars and supermassive black holes, since they are theorised to be sensitive to the metallicity of the gas out of which they form. The dynamic model is, however, not for galaxy evolution studies only. It can be applied to all astrophysical hydrodynamics simulations, including those modelling stellar interiors, planetary formation, and star formation. △ Less

Submitted 10 December, 2018; v1 submitted 30 July, 2018; originally announced July 2018.

Comments: Accepted to MNRAS December 07 2018. 23 pages, 9 figures

The Complete Local Volume Groups Sample - II. A study of the Central Radio Galaxies in the High$-$Richness Sub-sample

Authors: Konstantinos Kolokythas, Ewan O'Sullivan, Somak Raychaudhury, Simona Giacintucci, Myriam Gitti, Arif Babul

Abstract: We present a study of the radio properties of the dominant early-type galaxies in 26 galaxy groups, the high-richness sub-sample of the Complete Local-volume Groups Sample (CLoGS). Combining new 610 MHz and 235 MHz observations of 21 groups from the Giant Metrewave Radio Telescope (GMRT) with archival GMRT and Very Large Array (VLA) survey data, we find a high detection rate, with 92% of the domin… ▽ More We present a study of the radio properties of the dominant early-type galaxies in 26 galaxy groups, the high-richness sub-sample of the Complete Local-volume Groups Sample (CLoGS). Combining new 610 MHz and 235 MHz observations of 21 groups from the Giant Metrewave Radio Telescope (GMRT) with archival GMRT and Very Large Array (VLA) survey data, we find a high detection rate, with 92% of the dominant galaxies hosting radio sources. The sources have a wide range of luminosities, 10$^{20}$ $-$ 10$^{24}$ W/Hz in the 235 and 610 MHz bands. The majority (54%) are point-like, but 23% have radio jets, and another 15% are diffuse radio sources with no clear jet/lobe morphology. The spectral index of the detected radio sources ranges from very flat values of ~0.2 to typical radio synchrotron spectra of ~0.9 with only two presenting steep radio spectra with $α_{235}^{610}$ > 1. We find that jet sources are more common in X-ray bright groups, with radio non-detections found only in X-ray faint systems. Radio point sources appear in all group environments irrespective of their X-ray properties or spiral fraction. We estimate the mechanical power (Pcav) of the jet sources in the X-ray bright groups to be 10$^{41}$ $-$ 10$^{43}$ erg/s, with the two large-scale jet systems (NGC 193 and NGC 4261) showing jet powers two orders of magnitude greater than the radiative losses from the cool cores of their groups. This suggests that central AGN are not always in balance with cooling, but may instead produce powerful periodical bursts of feedback heating. △ Less

Submitted 29 July, 2018; originally announced July 2018.

Comments: Accepted for publication by MNRAS, 19 Manuscript pages with 8 tables and 8 figures, plus 13 pages of appendices

Cold gas in a complete sample of group-dominant early-type galaxies

Authors: E. O'Sullivan, F. Combes, P. Salomé, L. P. David, A. Babul, J. M. Vrtilek, J. Lim, V. Olivares, S. Raychaudhury, G. Schellenberger

Abstract: We present IRAM 30m and APEX telescope observations of CO(1-0) and CO(2-1) lines in 36 group-dominant early-type galaxies, completing our molecular gas survey of dominant galaxies in the Complete Local-volume Groups Sample. We detect CO emission in 12 of the galaxies at >4sigma significance, with molecular gas masses in the range 0.01-6x10^8 Msol, as well as CO in absorption in the non-dominant gr… ▽ More We present IRAM 30m and APEX telescope observations of CO(1-0) and CO(2-1) lines in 36 group-dominant early-type galaxies, completing our molecular gas survey of dominant galaxies in the Complete Local-volume Groups Sample. We detect CO emission in 12 of the galaxies at >4sigma significance, with molecular gas masses in the range 0.01-6x10^8 Msol, as well as CO in absorption in the non-dominant group member galaxy NGC 5354. In total 21 of the 53 CLoGS dominant galaxies are detected in CO and we confirm our previous findings that they have low star formation rates (0.01-1 Msol/yr) but short depletion times (<1Gyr) implying rapid replenishment of their gas reservoirs. Comparing molecular gas mass with radio luminosity, we find that a much higher fraction of our group-dominant galaxies (60+-16%) are AGN-dominated than is the case for the general population of ellipticals, but that there is no clear connection between radio luminosity and the molecular gas mass. Using data from the literature, we find that at least 27 of the 53 CLoGS dominant galaxies contain HI, comparable to the fraction of nearby non-cluster early type galaxies detected in HI and significantly higher that the fraction in the Virgo cluster. We see no correlation between the presence of an X-ray detected intra-group medium and molecular gas in the dominant galaxy, but find that the HI-richest galaxies are located in X-ray faint groups. Morphological data from the literature suggests the cold gas component most commonly takes the form of a disk, but many systems show evidence of galaxy-galaxy interactions, indicating that they may have acquired their gas through stripping or mergers. We provide improved molecular gas mass estimates for two galaxies previously identified as being in the centres of cooling flows, NGC 4636 and NGC 5846, and find that they are relatively molecular gas poor compared to our other detected systems. △ Less

Submitted 14 September, 2018; v1 submitted 24 July, 2018; originally announced July 2018.

Comments: 20 pages, 12 PDF figures, accepted by A&A. v2 adds references and acknowledgements, and corrects minor typographical errors identified at the proof stage

Journal ref: A&A 618, A126 (2018)

Introducing RomulusC: A Cosmological Simulation of a Galaxy Cluster with Unprecedented Resolution

Authors: Michael Tremmel, Thomas R. Quinn, Angelo Ricarte, Arif Babul, Urmila Chadayammuri, Priyamvada Natarajan, Daisuke Nagai, Andrew Pontzen, Marta Volonteri

Abstract: We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to st… ▽ More We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies. △ Less

Submitted 10 December, 2018; v1 submitted 4 June, 2018; originally announced June 2018.

Comments: Submitted to MNRAS, 29 pages, 19 figures, Accepted to MNRAS

A massive core for a cluster of galaxies at a redshift of 4.3

Authors: T. B. Miller, S. C. Chapman, M. Aravena, M. L. N. Ashby, C. C. Hayward, J. D. Vieira, A. Weiß, A. Babul, M. Béthermin, C. M. Bradford, M. Brodwin, J. E. Carlstrom, Chian-Chou Chen, D. J. M. Cunningham, C. De Breuck, A. H. Gonzalez, T. R. Greve, J. Harnett, Y. Hezaveh, K. Lacaille, K. C. Litke, J. Ma, M. Malkan, D. P. Marrone, W. Morningstar , et al. (13 additional authors not shown)

Abstract: Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters wit… ▽ More Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: To appear in April 26 issue of Nature

A limit on the warm dark matter particle mass from the redshifted 21 cm absorption line

Authors: Mohammadtaher Safarzadeh, Evan Scannapieco, Arif Babul

Abstract: The recent EDGES collaboration detection of an absorption signal at a central frequency of $ν= 78 \pm 1$ MHz points to the presence of a significant Lyman-$α$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_χ$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lo… ▽ More The recent EDGES collaboration detection of an absorption signal at a central frequency of $ν= 78 \pm 1$ MHz points to the presence of a significant Lyman-$α$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_χ$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lower limit can be placed on $m_χ,$ based on the presence of a sufficiently strong Ly-$α$ background due to star formation at $z=18$. Our results show that the coupling the spin temperature to the gas through Ly-$α$ pumping requires a minimum mass of $m_χ>3$ keV if atomic cooling halos dominate the star formation rate at $z=18,$ and $m_χ>2$ keV if ${\rm\,H_2}$ cooling halos also form stars efficiently at this redshift. These limits match or exceed the most stringent limits cited to date in the literature, even in the face of the many uncertainties regarding star-formation at high redshift. △ Less

Submitted 16 May, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: Accepted for publication in ApJ Letters

Diversity in the stellar velocity dispersion profiles of a large sample of Brightest Cluster Galaxies $z\leq0.3$

Authors: S. I. Loubser, H. Hoekstra, A. Babul, E. O'Sullivan

Abstract: We analyse spatially-resolved deep optical spectroscopy of Brightest Cluster Galaxies (BCGs) located in 32 massive clusters with redshifts of 0.05 $\leq z \leq$ 0.30, to investigate their velocity dispersion profiles. We compare these measurements to those of other massive early-type galaxies, as well as central group galaxies, where relevant. This unique, large sample extends to the most extreme… ▽ More We analyse spatially-resolved deep optical spectroscopy of Brightest Cluster Galaxies (BCGs) located in 32 massive clusters with redshifts of 0.05 $\leq z \leq$ 0.30, to investigate their velocity dispersion profiles. We compare these measurements to those of other massive early-type galaxies, as well as central group galaxies, where relevant. This unique, large sample extends to the most extreme of massive galaxies, spanning M$_{K}$ between -25.7 to -27.8 mag, and host cluster halo mass M$_{500}$ up to 1.7 $\times$ 10$^{15}$ M$_{\odot}$. To compare the kinematic properties between brightest group and cluster members, we analyse similar spatially-resolved long-slit spectroscopy for 23 nearby Brightest Group Galaxies (BGGs) from the Complete Local-Volume Groups Sample (CLoGS). We find a surprisingly large variety in velocity dispersion slopes for BCGs, with a significantly larger fraction of positive slopes, unique compared to other (non-central) early-type galaxies as well as the majority of the brightest members of the groups. We find that the velocity dispersion slopes of the BCGs and BGGs correlate with the luminosity of the galaxies, and we quantify this correlation. It is not clear whether the full diversity in velocity dispersion slopes that we see is reproduced in simulations. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 25 pages, 17 figures, accepted for publication in MNRAS

Cool-Core Clusters : Role of BCG, Star Formation & AGN-Driven Turbulence

Authors: Deovrat Prasad, Prateek Sharma, Arif Babul

Abstract: Recent analysis shows that it is important to explicitly include the gravitational potential of the central brightest central galaxy (BCG) to infer the acceleration due to gravity ($g$) and the free-fall time ($t_{\rm ff} \equiv [2r/g]^{1/2}$) in cool cluster cores. Accurately measuring $t_{\rm ff}$ is crucial because according to numerical simulations cold gas condensation and strong feedback occ… ▽ More Recent analysis shows that it is important to explicitly include the gravitational potential of the central brightest central galaxy (BCG) to infer the acceleration due to gravity ($g$) and the free-fall time ($t_{\rm ff} \equiv [2r/g]^{1/2}$) in cool cluster cores. Accurately measuring $t_{\rm ff}$ is crucial because according to numerical simulations cold gas condensation and strong feedback occur in cluster cores with min($t_{\rm cool}/t_{\rm ff}$) below a threshold value close to 10. Recent observations which include the BCG gravity show that the observed threshold in min($t_{\rm cool}/t_{\rm ff}$) lies at a somewhat higher value, close to 10-30; there are only a few clusters in which this ratio falls much below 10. In this paper we compare numerical simulations of feedback AGN (Active Galactic Nuclei) jets interacting with the intracluster medium (ICM), with and without a BCG potential. We find that, for a fixed feedback efficiency, the presence of a BCG does not significantly affect the temperature but increases (decreases) the core density (entropy) on average. Most importantly, min($t_{\rm cool}/t_{\rm ff}$) is only affected slightly by the inclusion of the BCG gravity. Also notable is that the lowest value of min($t_{\rm cool}/t_{\rm ff}$) in the NFW+BCG runs are about twice larger than in the NFW runs. We also look at the role of depletion of cold gas due to star formation and show that it only affects the rotationally dominant component (torus), while the radially dominant component (which regulates the feedback cycle) remains largely unaffected. The distribution of metals due to AGN jets in our simulations is predominantly along the jet direction and the radial spread of metals is less. We also show that the turbulence in cool core clusters is weak, consistent with recent Hitomi results on Perseus cluster. △ Less

Submitted 16 March, 2018; v1 submitted 12 January, 2018; originally announced January 2018.

Comments: Submitted to ApJ; Suggestions and comments are welcome; Substantially updated from the previous version; Abstract abridged; 12 pages, 11 figures