Search

We present millimeter observations of the host galaxy of the most distant blazar known, VLASSJ041009.05-013919.88 (hereafter J0410-0139) at z=7, using ALMA and NOEMA observations. The ALMA data reveal a 2e42 erg/s [CII] 158um emission line at z=6.9964 with a [CII]-inferred star-formation rate of 58 Msun/yr. We estimate a dynamical mass of 4.6e9 Msun, implying a black hole mass to host a dynamical mass ratio of 0.15. The 238 GHz continuum (rest-frame IR) decreased by ~33% from the NOEMA to the ALMA observations taken ~10 months apart. The VLA 3-10 GHz radio flux densities showed a ~37% decrease in a similar time frame, suggesting a causal connection. At face value, J0410-0139 would have the lowest [CII]-to-IR luminosity ratio of a z>5.7 quasar reported to date (~1e-4). However, if only <20% of the measured IR luminosity were due to thermal emission from dust, the [CII]-to-IR luminosity ratio would be typical of (U)LIRGS, and the star formation rates derived from [CII] and IR luminosities would be consistent. These results provide further evidence that synchrotron emission significantly contributes to the observed rest-frame IR emission of J0410-0139, similar to what has been reported in some radio-loud AGN at z<1., Comment: Version after addressing referee report. See also companion paper: arXiv:2407.07236

Relativistic jets are thought to play a crucial role in the formation of massive galaxies and supermassive black holes. Here we report multi-wavelength and multi-epoch observations of the quasar VLASSJ0410-0139 at redshift z=7, powered by a 7e8 solar-mass black hole. Its radio variability, X-ray properties, and compact radio emission on parsec scales reveal that J0410-0139 is a blazar with a relativistic jet aligned with our line of sight. This blazar's existence implies that many more similar (unaligned) jetted sources must exist at z=7. One scenario is that we observe an intrinsically low-power radio jet, but we see it at high luminosity due to relativistic beaming effects. In this case, a large fraction (>80%) of the UV bright quasars must have a similar jet to match the number density expected from the UV quasar luminosity function. These jets can enhance the growth of supermassive black holes and substantially affect their host galaxies. However, the implications would be even more severe if the quasar belongs to the top 10% radio luminous quasars, as measured if the beaming enhancement is less than a factor of 10-15. In this scenario, there should be hundreds to thousands of radio-quiet quasars at z=7 with intrinsic properties similar to J0410-0139 -- in strong tension with the number density of bright quasars derived from their UV luminosity function. To reconcile these results, most black hole growth at z=7 must happen in an obscured phase, as some models predict. The existence of supermassive black holes in the epoch of reionization is facilitated by significant jet-enhanced or obscured super-Eddington accretion., Comment: Submitted

We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$\alpha$, H$\beta$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe., Comment: 15 pages, 16 figures. Accepted for publication in A&A

The evolution of galaxies depends on their masses and local environments; understanding when and how environmental quenching starts to operate remains a challenge. Furthermore, studies of the high-redshift regime have been limited to massive cluster members, owing to sensitivity limits or small fields of views when the sensitivity is sufficient, intrinsically biasing the picture of cluster evolution. In this work, we use stacking to investigate the average star formation history of more than 10,000 groups and clusters drawn from the Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2). Our analysis covers near ultraviolet to far infrared wavelengths, for galaxy overdensities at $0.5 \lesssim z \lesssim 2.54$. We employ SED fitting to measure the specific star formation rates (sSFR) in four annular apertures with radii between 0 and 1000 kpc. At $z \gtrsim 1.6$, the average sSFR evolves similarly to the field in both the core and the cluster outskirts. Between $\overline{z} = 1.60$ and $\overline{z} = 1.35$, the sSFR in the core drops sharply, and continues to fall relative to the field sSFR at lower redshifts. We interpret this change as evidence that the impact of environmental quenching dramatically increases at $z \sim 1.5$, with the short time span of the transition suggesting that the environmental quenching mechanism dominant at this redshift operates on a rapid timescale. We find indications that the sSFR may decrease with increasing host halo mass, but lower-scatter mass tracers than the signal-to-noise ratio (S/N) are needed to confirm this relationship., Comment: 22 pages, 14 figures, accepted for publication in ApJ (corrected Popescu et al., 2023 sSFR on Figure 14)

We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80-month (NSS80) catalog has an unprecedented $\sim$ 62 Ms of effective exposure time over 894 unique fields (a factor of three increase over the 40-month catalog), with an areal coverage of $\sim $36 deg$^2$, larger than all NuSTAR extragalactic surveys. NSS80 provides 1274 hard X-ray sources in the $3-24$ keV band (822 new detections compared to the previous 40-month catalog). Approximately 76% of the NuSTAR sources have lower-energy ($<10$ keV) X-ray counterparts from Chandra, XMM-Newton, and Swift-XRT. We have undertaken an extensive campaign of ground-based spectroscopic follow-up to obtain new source redshifts and classifications for 427 sources. Combining these with existing archival spectroscopy provides redshifts for 550 NSS80 sources, of which 547 are classified. The sample is primarily composed of active galactic nuclei (AGN), detected over a large range in redshift ($z$ = 0.012-3.43), but also includes 58 spectroscopically confirmed Galactic sources. In addition, five AGN/galaxy pairs, one dual AGN system, one BL Lac candidate, and a hotspot of 4C 74.26 (radio quasar) have been identified. The median rest-frame $10-40$ keV luminosity and redshift of the NSS80 are $\langle{L_\mathrm{10-40 keV}}\rangle$ = 1.2 $\times$ 10$^{44}$ erg s$^{-1}$ and $\langle z \rangle = 0.56$. We investigate the optical properties and construct composite optical spectra to search for subtle signatures not present in the individual spectra, finding an excess of redder BL AGN compared to optical quasar surveys predominantly due to the presence of the host-galaxy and, at least in part, due to dust obscuration., Comment: Accepted for publication in ApJ:S. 57 pages, 32 figures

The Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2) is a new survey designed as the successor of the original MaDCoWS survey. MaDCoWS2 improves upon its predecessor by using deeper optical and infrared data and a more powerful detection algorithm (PZWav). As input to the search, we use grz photometry from DECaLS in combination with W1 and W2 photometry from the CatWISE2020 catalog to derive the photometric redshifts with full redshift probability distribution functions for WISE-selected galaxies. Cluster candidates are then detected using the PZWav algorithm to find three-dimensional galaxy overdensities from the sky positions and photometric redshifts. This paper provides the first MaDCoWS2 data release, covering 1461 (1838 without masking) deg^2 centered on the Hyper-SuprimeCam Subaru Strategic Program equatorial fields. Within this region, we derive a catalog of 22,970 galaxy cluster candidates detected at S/N>5. These clusters span the redshift range 0.11.5. We compare MaDCoWS2 to six existing catalogs in the area. We rediscover 60%-92% of the clusters in these surveys at S/N>5. The medians of the absolute redshift offset are <0.02 relative to these surveys, while the standard deviations are less than 0.06. The median offsets between the detection position from MaDCoWS2 and other surveys are less than 0.25 Mpc. We quantify the relation between S/N and gas mass, total mass, luminosity, and richness from other surveys using a redshift-dependent power law relation. We find that the S/N-richness relation exhibits the lowest scatter., Comment: 27 pages, 7 figures. Typo corrected. Accepted for publication in ApJ

We investigate the potential of large language models (LLMs) to disentangle text variables--to remove the textual traces of an undesired forbidden variable in a task sometimes known as text distillation and closely related to the fairness in AI and causal inference literature. We employ a range of various LLM approaches in an attempt to disentangle text by identifying and removing information about a target variable while preserving other relevant signals. We show that in the strong test of removing sentiment, the statistical association between the processed text and sentiment is still detectable to machine learning classifiers post-LLM-disentanglement. Furthermore, we find that human annotators also struggle to disentangle sentiment while preserving other semantic content. This suggests there may be limited separability between concept variables in some text contexts, highlighting limitations of methods relying on text-level transformations and also raising questions about the robustness of disentanglement methods that achieve statistical independence in representation space., Comment: To appear as: Nicolas Audinet de Pieuchon, Adel Daoud, Connor T. Jerzak, Moa Johansson, Richard Johansson. Can Large Language Models (or Humans) Disentangle Text? In: Sixth Workshop on NLP and Computational Social Science at NAACL, 2024

We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm \mu m$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$\beta$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,H\beta}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$\alpha$-based black hole mass from the same spectrum ($M_{\rm BH, H\alpha}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $\lambda 2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $\lambda_{\rm Edd,H\beta}\sim 0.67$ ($\lambda_{\rm Edd,H\alpha}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$\beta$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $\lambda 5007$ emission with a velocity offset $\Delta v_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $\Delta v_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$\beta$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail., Comment: 17 pages, 9 figures, accepted for publication in A&A

We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII). They are combined with the published absorbing systems in Davies et al. (2023a) at the same redshift range to form of a sample of nine metal absorbers at z = 6.03 to 6.49. We identify eight galaxies within 1000 km s$^{-1}$ and 350 kpc around the absorbing gas from the ASPIRE spectroscopic data, with their redshifts secured by [OIII]($\lambda\lambda$4959, 5007) doublets and H$\beta$ emission lines. Our spectral energy distribution fitting indicates that the absorbing galaxies have stellar mass ranging from 10$^{7.2}$ to 10$^{8.8}M_{\odot}$ and metallicity between 0.02 and 0.4 solar. Notably, the z = 6.2713 system in the J0305-3150 field resides in a galaxy overdensity region, which contains two (tentatively) merging galaxies within 350 kpc and seven galaxies within 1 Mpc. We measure the relative abundances of $\alpha$ elements to iron ([$\alpha$/Fe]) and find that the CGM gas in the most overdense region exhibits a lower [$\alpha$/Fe] ratio. Our modeling of the galaxy's chemical abundance favors a top-heavy stellar initial mass function, and hints that we may be witnessing the contribution of the first generation Population III stars to the CGM at the end of reionization epoch., Comment: 21 pages, 4 figures in the main text. Accepted for publication in ApJL

With two central galaxies engaged in a major merger and a remarkable chain of 19 young stellar superclusters wound around them in projection, the galaxy cluster SDSS J1531+3414 ($z=0.335$) offers an excellent laboratory to study the interplay between mergers, AGN feedback, and star formation. New Chandra X-ray imaging reveals rapidly cooling hot ($T\sim 10^6$ K) intracluster gas, with two "wings" forming a concave density discontinuity near the edge of the cool core. LOFAR $144$ MHz observations uncover diffuse radio emission strikingly aligned with the "wings," suggesting that the "wings" are actually the opening to a giant X-ray supercavity. The steep radio emission is likely an ancient relic of one of the most energetic AGN outbursts observed, with $4pV > 10^{61}$ erg. To the north of the supercavity, GMOS detects warm ($T\sim 10^4$ K) ionized gas that enshrouds the stellar superclusters but is redshifted up to $+ 800$ km s$^{-1}$ with respect to the southern central galaxy. ALMA detects a similarly redshifted $\sim 10^{10}$ M$_\odot$ reservoir of cold ($T\sim 10^2$ K) molecular gas, but it is offset from the young stars by $\sim 1{-}3$ kpc. We propose that the multiphase gas originated from low-entropy gas entrained by the X-ray supercavity, attribute the offset between the young stars and the molecular gas to turbulent intracluster gas motions, and suggest that tidal interactions stimulated the "beads on a string" star formation morphology., Comment: Accepted by ApJ, 36 pages, 23 figures

A complete accounting of nearby objects -- from the highest-mass white dwarf progenitors down to low-mass brown dwarfs -- is now possible, thanks to an almost complete set of trigonometric parallax determinations from Gaia, ground-based surveys, and Spitzer follow-up. We create a census of objects within a Sun-centered sphere of 20-pc radius and check published literature to decompose each binary or higher-order system into its separate components. The result is a volume-limited census of $\sim$3,600 individual star formation products useful in measuring the initial mass function across the stellar ($<8 M_\odot$) and substellar ($\gtrsim 5 M_{Jup}$) regimes. Comparing our resulting initial mass function to previous measurements shows good agreement above 0.8$M_\odot$ and a divergence at lower masses. Our 20-pc space densities are best fit with a quadripartite power law, $\xi(M) = dN/dM \propto M^{-\alpha}$ with long-established values of $\alpha = 2.3$ at high masses ($0.55 < M < 8.00 M_\odot$) and $\alpha = 1.3$ at intermediate masses ($0.22 < M < 0.55 M_\odot$), but at lower masses we find $\alpha = 0.25$ for $0.05 < M <0.22 M_\odot$ and $\alpha = 0.6$ for $0.01 < M < 0.05 M_\odot$. This implies that the rate of production as a function of decreasing mass diminishes in the low-mass star/high-mass brown dwarf regime before increasing again in the low-mass brown dwarf regime. Correcting for completeness, we find a star to brown dwarf number ratio of, currently, 4:1, and an average mass per object of 0.41 $M_\odot$., Comment: Accepted for publication in The Astrophysical Journal Supplement Series. 123 pages with four ancillary files

Jets powered by AGN in the early Universe ($z \gtrsim 6$) have the potential to not only define the trajectories of the first-forming massive galaxies but to enable the accelerated growth of their associated SMBHs. Under typical assumptions, jets could even rectify observed quasars with light seed formation scenarios; however, not only are constraints on the parameters of the first jets lacking, observations of these objects are scarce. Owing to the significant energy density of the CMB at these epochs capable of quenching radio emission, observations will require powerful, high angular resolution X-ray imaging to map and characterize these jets. As such, \textit{AXIS} will be necessary to understand early SMBH growth and feedback., Comment: 7 Pages, 2 Figures. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website (http://axis.astro.umd.edu/) with a mission overview here: arXiv:2311.00780

We present radio observations of 24 confirmed and candidate strongly lensed quasars identified by the Gaia Gravitational Lenses (GraL) working group. We detect radio emission from 8 systems in 5.5 and 9 GHz observations with the Australia Telescope Compact Array (ATCA), and 12 systems in 6 GHz observations with the Karl G. Jansky Very Large Array (VLA). The resolution of our ATCA observations is insufficient to resolve the radio emission into multiple lensed images, but we do detect multiple images from 11 VLA targets. We have analysed these systems using our observations in conjunction with existing optical measurements, including measuring offsets between the radio and optical positions, for each image and building updated lens models. These observations significantly expand the existing sample of lensed radio quasars, suggest that most lensed systems are detectable at radio wavelengths with targeted observations, and demonstrate the feasibility of population studies with high resolution radio imaging.

We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $\mu$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $\delta_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3

Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 \r{A}. The profiles of these quasars' broad H$\beta$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$\beta$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$\beta$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $\lambda\lambda$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,\mu$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties., Comment: 13 pages, 4 figures, accepted for publication in ApJL

We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low energy X-ray properties and the larger AGN population. However some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs.

We present multi-wavelength high-spatial resolution (~0.1'', 70 pc) observations of UGC 4211 at z=0.03474, a late-stage major galaxy merger at the closest nuclear separation yet found in near-IR imaging (0.32'', ~230 pc projected separation). Using Hubble Space Telescope/STIS, VLT/MUSE+AO, Keck/OSIRIS+AO spectroscopy, and ALMA observations, we show that the spatial distribution, optical and NIR emission lines, and millimeter continuum emission are all consistent with both nuclei being powered by accreting supermassive black holes (SMBHs). Our data, combined with common black hole mass prescriptions, suggests that both SMBHs have similar masses, log MBH~8.1 (south) and log MBH~8.3 (north), respectively. The projected separation of 230 pc (~6X the black hole sphere of influence) represents the closest-separation dual AGN studied to date with multi-wavelength resolved spectroscopy and shows the potential of nuclear (<50 pc) continuum observations with ALMA to discover hidden growing SMBH pairs. While the exact occurrence rate of close-separation dual AGN is not yet known, it may be surprisingly high, given that UGC 4211 was found within a small, volume-limited sample of nearby hard X-ray detected AGN. Observations of dual SMBH binaries in the sub-kpc regime at the final stages of dynamical friction provide important constraints for future gravitational wave observatories., Comment: 19 pages, 7 figures, 2 tables, published in ApJL

The identification of bright quasars at z>6 enables detailed studies of supermassive black holes, massive galaxies, structure formation, and the state of the intergalactic medium within the first billion years after the Big Bang. We present the spectroscopic confirmation of 55 quasars at redshifts 5.6

We present the $z\!\approx\!6$ type-1 quasar luminosity function (QLF) based on the Pan-STARRS1 (PS1) quasar survey. The PS1 sample includes 125 quasars at $z\approx5.7-6.2$ with $-28\lesssim M_{1450}\lesssim-25$. Complemented by 48 fainter quasars from the SHELLQs survey, we evaluate the $z\approx6$ QLF over $-28\lesssim M_{1450}\lesssim-22$. Adopting a double power law with an exponential evolution of the quasar density ($\Phi(z)\propto10^{k(z-6)}$; $k=-0.7$), we use a maximum likelihood method to model our data. We find a break magnitude of $M^*=-26.38_{-0.60}^{+0.79}\,\text{mag}$, a faint end slope of $\alpha=-1.70_{-0.19}^{+0.29}$, and a steep bright end slope with $\beta=-3.84_{-1.21}^{+0.63}$. Based on our new QLF model we determine the quasar comoving spatial density at $z\!\approx\!6$ to be $n( M_{1450}<-26)=1.16_{-0.12}^{+0.13}\,\text{cGpc}^{-3}$. In comparison with the literature, we find the quasar density to evolve with a constant value of $k\approx-0.7$ from $z\approx7$ to $z\approx4$. Additionally, we derive an ionizing emissivity of $\epsilon_{912}(z=6) =7.23_{-1.02}^{+1.65}\times 10^{22}\,\text{erg}\,\text{s}^{-1}\text{Hz}^{-1}\text{cMpc}^{-3}$ based on the QLF measurement. Given standard assumptions and the recent measurement of the mean free path of Becker et al. (2021) at $z\approx6$ we calculate an HI photoionizing rate of $\Gamma_{\text{HI}}(z{=}6) \approx 6 \times10^{-16}\,\text{s}^{-1}$, strongly disfavoring a dominant role of quasars in hydrogen reionization., Comment: accepted by the Astrophysical Journal; see companion paper by Ba\~nados et al. (arXiv:2212.04452)

We investigate the strong-lensing cluster Abell 370 (A370) using a wide Integral Field Unit (IFU) spectroscopic mosaic from the Multi-Unit Spectroscopic Explorer (MUSE). IFU spectroscopy provides significant insight into the structure and mass content of galaxy clusters, yet IFU-based cluster studies focus almost exclusively on the central Einstein-radius region. Covering over 14 arcmin$^2$, the new MUSE mosaic extends significantly beyond the A370 Einstein radius, providing, for the first time, a detailed look at the cluster outskirts. Combining these data with wide-field, multi-band Hubble Space Telescope (HST) imaging from the BUFFALO project, we analyse the distribution of objects within the cluster and along the line of sight. Identifying 416 cluster galaxies, we use kinematics to trace the radial mass profile of the halo, providing a mass estimate independent from the lens model. We also measure radially-averaged properties of the cluster members, tracking their evolution as a function of infall. Thanks to the high spatial resolution of our data, we identify six cluster members acting as galaxy-galaxy lenses, which constrain localized mass distributions beyond the Einstein radius. Finally, taking advantage of MUSE's 3D capabilities, we detect and analyse multiple spatially extended overdensities outside of the cluster that influence lensing-derived halo mass estimates. We stress that much of this work is only possible thanks to the robust, extended IFU coverage, highlighting its importance even in less optically dense cluster regions. Overall, this work showcases the power of combining HST+MUSE, and serves as the initial step towards a larger and wider program targeting several clusters., Comment: 20 pages, 12 figures, 3 tables (including 1 in an appendix). Accepted in MNRAS. Data Release products available at https://astro.dur.ac.uk/~hbpn39/pilot-wings.html

We present stellar mass fractions and composite luminosity functions (LFs) for a sample of \Ncl\ clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) at a redshift range of $0.951 \leq z \leq 1.43$. Using SED fitting of optical and deep mid-infrared photometry, we establish the membership of objects along the lines-of-sight to these clusters and calculate the stellar masses of member galaxies. We find stellar mass fractions for these clusters largely consistent with previous works, including appearing to display a negative correlation with total cluster mass. We measure a composite $3.6~\mathrm{\mu m}$ LF down to $m^*+2.5$ for all 12 clusters. Fitting a Schechter function to the LF, we find a characteristic $3.6~\mathrm{\mu m}$ magnitude of $m^*=19.83\pm0.12$ and faint-end slope of $\alpha=-0.81\pm0.10$ for the full sample at a mean redshift of $\bar{z} = 1.18$. We also divide the clusters into high- and low-redshift bins at $\bar{z}=1.29$ and $\bar{z}=1.06$ respectively and measure a composite LF for each bin. We see a small, but statistically significant evolution in $m^*$ and $\alpha$ -- consistent with passive evolution -- when we study the joint fit to the two parameters, which is probing the evolution of faint cluster galaxies at $z\sim1$. This highlights the importance of deep IR data in studying the evolution of cluster galaxy populations at high-redshift., Comment: 13 pages, 11 figures, Submitted to ApJ

The $z=6.327$ quasar SDSS J010013.02+280225.8 (hereafter J0100+2802) is believed to be powered by a black hole more massive than $10^{10}\ {\rm M}_\odot$, making it the most massive black hole known in the first billion years of the Universe. However, recent high-resolution ALMA imaging shows four structures at the location of this quasar, potentially implying that it is lensed with a magnification of $\mu\sim450$ and thus its black hole is significantly less massive. Furthermore, for the underlying distribution of magnifications of $z\gtrsim6$ quasars to produce such an extreme value, theoretical models predict that a larger number of quasars in this epoch should be lensed, implying further overestimates of early black hole masses. To provide an independent constraint on the possibility that J0100+2802 is lensed, we re-analyzed archival XMM-Newton observations of the quasar and compared the expected ratios of X-ray luminosity to rest-frame UV and IR luminosities. For both cases, J0100+2802's X-ray flux is consistent with the no-lensing scenario; while this could be explained by J0100+2802 being X-ray faint, we find it does not have the X-ray or optical spectral features expected for an X-ray faint quasar. Finally, we compare the overall distribution of X-ray fluxes for known, typical $z\gtrsim6$ quasars. We find a $3\sigma$ tension between the observed and predicted X-ray-to-UV flux ratios when adopting the magnification probability distribution required to produce a $\mu=450$ quasar., Comment: 8 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

We present XMM-Newton X-ray observations of nine confirmed lensed quasars at $1 \lesssim z \lesssim 3$ identified by the Gaia Gravitational Lens program. Eight systems are strongly detected, with 0.3--8.0 keV fluxes $F_{0.3-8.0} \gtrsim 5 \times 10^{-14}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm s}^{-1}$. Modeling the X-ray spectra with an absorbed power law, we derive power law photon indices and 2--10 keV luminosities for the eight detected quasars. In addition to presenting sample properties for larger quasar population studies and for use in planning for future caustic crossing events, we also identify three quasars of interest: a quasar that shows evidence of flux variability from previous ROSAT observations, the most closely-separated individual lensed sources resolved by XMM-Newton, and one of the X-ray brightest quasars known at $z>3$. These sources represent the tip of discovery that will be enabled by SRG/eROSITA., Comment: 11 pages with 11 pages of appendices. Submitted to the Astrophysical Journal. Comments welcome

The interactions between radio jets and the interstellar medium play a defining role for the co-evolution of central supermassive black holes and their host galaxies, but observational constraints on these feedback processes are still very limited at redshifts $z > 2$. We investigate the radio-loud quasar PSO J352.4034-15.3373 at $z \sim 6$ at the edge of the Epoch of Reionization. This quasar is among the most powerful radio emitters and the first one with direct evidence of extended radio jets ($\sim$1.6 kpc) at these high redshifts. We analyze NOEMA and ALMA millimeter data targeting the CO (6-5) and [CII] far-infrared emission lines, respectively, and the underlying continuum. The broad $440\pm 80$ km s$^{-1}$ and marginally resolved [CII] emission line yields a systemic redshift of $z\!=\!5.832 \pm 0.001$. Additionally, we report a strong 215 MHz radio continuum detection, $88\pm 7$ mJy, using the GMRT. This measurement significantly improves the constraints at the low-frequency end of the spectral energy distribution of this quasar. In contrast to what is typically observed in high-redshift radio-quiet quasars, we show that cold dust emission alone cannot reproduce the millimeter continuum measurements. This is evidence that the strong synchrotron emission from the quasar contributes substantially to the emission even at millimeter (far-infrared in the rest-frame) wavelengths. This quasar is an ideal system to probe the effects of radio jets during the formation of a massive galaxy within the first Gyr of the Universe., Comment: 14 pages, 5 figures

We present an analysis of the intracluster light in the Frontier Field Cluster MACS J1149.5+2223 (z=0.544), which combines new and archival Hubble WFC3/IR imaging to provide continuous radial coverage out to 2.8 Mpc from the brightest cluster galaxy. Employing careful treatment of potential systematic biases and using data at the largest radii to determine the background sky level, we reconstruct the surface brightness profile out to a radius of 2 Mpc. This radius is the largest to which the intracluster light (ICL) has been measured for an individual cluster. Within this radius, we measure a total luminosity of 1.5e13 Lsun for the brightest cluster galaxy plus ICL light. From the profile and its logarithmic slope, we identify the transition from the brightest cluster galaxy to ICL at r~70 kpc. Remarkably, we also detect an inflection in the profile centered in the 1.2-1.7 Mpc (0.37-0.52 r200m) radial bin, a signature of an infall caustic in the stellar distribution. Based upon the shape and strength of the feature, we interpret it as potentially being at the splashback radius, although the radius is smaller than theoretical predictions. If this is the splashback radius, then it is the first such detection in the ICL and the first detection of the splashback radius for an individual cluster. Similar analyses should be possible with the other Frontier Field clusters, and eventually with clusters from the Euclid and Roman missions., Comment: 8 pages, 5 figures, submission to MNRAS

The discovery of hundreds of QSOs in the first Gyr of the Universe powered by already grown SMBHs challenges our knowledge of SMBH formation. In particular, investigations of $z>6$ QSOs presenting notable properties can provide unique information on the physics of fast SMBH growth in the early universe. We present the results of follow-up observations of the $z=6.515$ radio-quiet QSO PSO167-13, which is interacting with a close companion galaxy. The PSO167-13 system has been recently proposed to host the first heavily obscured X-ray source at high redshift. We observed PSO167-13 with Chandra/ACIS-S (177 ks), and obtained new spectroscopic observations (7.2 h) with Magellan/FIRE. No significant X-ray emission is detected from the PSO167-13 system, suggesting that the obscured X-ray source previously tentatively detected was either due to a strong background fluctuation or is highly variable. The upper limit (90% confidence level) on the X-ray emission of PSO167-13 ($L_{2-10\,\mathrm{keV}}<8.3\times10^{43}\,\mathrm{erg s^{-1}}$) is the lowest available for a $z>6$ QSO. The ratio between the X-ray and UV luminosity of $\alpha_{ox}<-1.95$ makes PSO167-13 a strong outlier from the $\alpha_{ox}-L_{UV}$ and $L_X-L_{\mathrm{bol}}$ relations. In particular, its X-ray emission is $>6$ times weaker than the expectation based on its UV luminosity. The new Magellan/FIRE spectrum of PSO167-13 is strongly affected by the unfavorable sky conditions, but the tentatively detected C IV and Mg II emission lines appear strongly blueshifted. The most plausible explanations for the X-ray weakness of PSO167-13 are intrinsic weakness or small-scale absorption by Compton-thick material. The possible strong blueshift of its emission lines hints at the presence of nuclear winds, which could be related to its X-ray weakness., Comment: Accepted for publication on A&A

We present deep (265 ks) Chandra X-ray observations of PSO J352.4034$-$15.3373, a quasar at z=5.831 that, with a radio-to-optical flux ratio of R>1000, is one of the radio-loudest quasars in the early universe and is the only quasar with observed extended radio jets of kpc-scale at $z \gtrsim 6$. Modeling the X-ray spectrum of the quasar with a power law, we find a best fit of $\Gamma = 1.99^{+0.29}_{-0.28}$, leading to an X-ray luminosity of $L_{2-10} = 1.26^{+0.45}_{-0.33} \times 10^{45}\ {\rm erg}\ {\rm s}^{-1}$ and an X-ray to UV brightness ratio of $\alpha_{\rm OX} = -1.36 \pm 0.11$. We identify a diffuse structure 50 kpc (${\sim}8^{\prime\prime}$) to the NW of the quasar along the jet axis that corresponds to a $3\sigma$ enhancement in the angular density of emission and can be ruled out as a background fluctuation with a probability of P=0.9985. While with few detected photons the spectral fit of the structure is uncertain, we find that it has a luminosity of $L_{2-10}\sim10^{44}\ {\rm erg}\ {\rm s}^{-1}$. These observations therefore potentially represent the most distant quasar jet yet seen in X-rays. We find no evidence for excess X-ray emission where the previously-reported radio jets are seen (which have an overall linear extent of $0.^{\prime\prime}28$), and a bright X-ray point source located along the jet axis to the SE is revealed by optical and NIR imaging to not be associated with the quasar., Comment: 16 pages, 7 Figures. Accepted for publication the Astrophysical Journal

Radio sources at the highest redshifts can provide unique information on the first massive galaxies and black holes, the densest primordial environments, and the epoch of reionization. The number of astronomical objects identified at z>6 has increased dramatically over the last few years, but previously only three radio-loud (R2500>10) sources had been reported at z>6, with the most distant being a quasar at z=6.18. Here we present the discovery and characterization of P172+18, a radio-loud quasar at z=6.823. This source has an MgII-based black hole mass of ~3x10^8 Msun and is one of the fastest accreting quasars, consistent with super-Eddington accretion. The ionized region around the quasar is among the largest measured at these redshifts, implying an active phase longer than the average lifetime of the z>6 quasar population. From archival data, there is evidence that its 1.4 GHz emission has decreased by a factor of two over the last two decades. The quasar's radio spectrum between 1.4 and 3.0 GHz is steep (alpha=-1.31) and has a radio-loudness parameter R2500~90. A second steep radio source (alpha=-0.83) of comparable brightness to the quasar is only 23.1" away (~120 kpc at z=6.82; projection probability <2%), but shows no optical or near-infrared counterpart. Further follow-up is required to establish whether these two sources are physically associated., Comment: submitted to ApJ on Nov 29, 2020; accepted on Jan 31, 2021. See the companion paper by Momjian et al

Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at $z\ge7.5$, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at $z=7.642$, J0313$-$1806, the most distant quasar yet known. This quasar has a bolometric luminosity of $3.6\times10^{13} L_\odot$. Deep spectroscopic observations reveal a SMBH with a mass of $(1.6\pm0.4) \times10^9M_\odot$ in this quasar. The existence of such a massive SMBH just $\sim$670 million years after the Big Bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in CIV and SiIV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted CIV emission line, indicate that there is a strong active galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum and [CII] emission from the quasar host galaxy, yielding an accurate redshift of $7.6423 \pm 0.0013$ and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of $\rm\sim 200 ~M_\odot ~yr^{-1}$ and a dust mass of $\sim7\times10^7~M_\odot$. Followup observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive galaxies., Comment: 4 figures, 1 table; accepted for publication in ApJL

In 2017, the North Carolina legislature passed House Bill 527 (now State Law 2017-196) in order to foster free, open inquiry in the state's colleges and universities. One of the provisions ordered the University of North Carolina (UNC) system Board of Governors to produce an annual report on two major categories of intellectual freedom: free speech and institutional neutrality. This report is the Martin Center's own analysis of free speech and institutional neutrality on the UNC system's 16 campuses. The report examines universities' choices for summer reading assignments and commencement speakers, whether they make public statements on controversial topics like climate change and diversity, the content of freshman orientation sessions, and universities' mission statements. It also records any barriers to student and faculty speech as well as any known disruptions or incidents that curtailed free speech in the past year. To complete this report, the Martin Center consulted schools' websites and newspaper articles. They also contacted student groups to determine whether they experienced any difficulties exercising their rights to free expression. The results of the Martin Center's analysis of free speech and institutional neutrality are presented in this report, organized by campus.

We present deep Chandra observations of PSO J231.6576$-$20.8335, a quasar at redshift z=6.59 with a nearby (${\sim}8$ proper kpc) companion galaxy. ALMA observed both the quasar and companion to be bright in [C II], and the system has significant extended Ly$\alpha$ emission around the quasar, suggesting that a galaxy merger is ongoing. Unlike previous studies of two similar systems, and despite observing the system with Chandra for 140 ks, we do not detect the companion in X-rays. The quasar itself is detected, but only $13.3^{+4.8}_{-3.7}$ net counts are observed. From a basic spectral analysis, the X-ray spectrum of the quasar is soft (hardness ratio of $\mathcal{HR} = -0.60_{-0.27}^{+0.17}$, power-law index of $\Gamma=2.6^{+1.0}_{-0.9}$), which results in a rest-frame X-ray luminosity comparable to other bright quasars ($L_{2-10} = 1.09^{+2.20}_{-0.70}\times 10^{45}\ \textrm{erg}\ \textrm{s}^{-1}$) despite the faint observed X-ray flux. We highlight two possible interpretations of this result: the quasar has a steep value of $\Gamma$ -- potentially related to observed ongoing Eddington accretion -- thereby pushing much of the emission out of our observed band, or the quasar has a more normal spectrum ($\Gamma{\sim}2$) but is therefore less X-ray luminous ($L_{2-10} \sim 0.6 \times 10^{45}\ \textrm{ erg}\ \textrm{ s}^{-1}$)., Comment: 11 pages, 4 figures. Accepted for publication in ApJ

We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank Telescope (GBT), the Spitzer Space Telescope, and the Dark Energy Camera Legacy Survey (DECaLS). We find that there are five radio sources which are all located within 500 kpc ($\sim$1$^{\prime}$) of the cluster center and have radio luminosities $P_{\mathrm{1.4GHz}}$ > 1.6$\times$10$^{24}$ W Hz$^{-1}$. The typical host galaxies are among the highest stellar mass galaxies in the cluster. The exceptional radio activity among the massive galaxy population appears to be linked to the dynamical state of the cluster. The galaxy distribution suggests an ongoing merger, with a subgroup found to the northwest of the main cluster. Further, two of the five sources are classified as bent-tail sources with one being a potential wide-angle tail (WAT)/hybrid morphology radio source (HyMoRS) indicating a dynamic environment. The cluster also lies in a region of the mass-richness plane occupied by other merging clusters in the Massive and Distant Clusters of WISE Survey (MaDCoWS). The data suggest that during the merger phase radio activity can be dramatically enhanced, which would contribute to the observed trend of increased radio activity in clusters with increasing redshift., Comment: 17 pages and 8 figures. Accepted in ApJ for publication

We present an archival analysis of Chandra X-ray observations for twelve nearby early-type galaxies hosting radio sources with radio power $>10^{23} \, \rm{W}~\rm{Hz}^{-1}$ at 1.4 GHz, similar to the radio power of the radio source in NGC 4261. Previously, in a similar analysis of eight nearby X-ray and optically-bright elliptical galaxies, Werner et al. 2012, found that NGC 4261 exhibited unusually low central gas entropy compared to the full sample. In the central 0.3 kpc of NGC 4261, the ratio of cooling time to freefall time ($t_{\rm{cool}}/t_{\rm ff}$) is less than $10$, indicating that cold clouds may be precipitating out of the hot ambient medium and providing fuel for accretion in the central region. NGC 4261 also hosts the most powerful radio source in the original sample. Because NGC 4261 may represent an important phase during which powerful feedback from a central active galactic nucleus (AGN) is fueled by multiphase condensation in the central kpc, we searched the Chandra archive for analogs to NGC 4261. We present entropy profiles of those galaxies as well as profiles of $t_{\rm{cool}}/t_{\rm ff}$. We find that one of them, IC 4296, exhibits properties similar to NGC 4261, including the presence of only single phase gas outside of $r \sim 2$ kpc and a similar central velocity dispersion. We compare the properties of NGC 4261 and IC 4296 to hydrodynamic simulations of AGN feedback fueled by precipitation. Over the course of those simulations, the single phase galaxy has an entropy gradient that remains similar to the entropy profiles inferred from our observations., Comment: 13 pages, 5 Figures, revised version of paper submitted to ApJ

The properties of galaxy clusters as a function of redshift can be utilized as an important cosmological tool. We present initial results from a program of follow-up observations of the Sunyaev-Zeldovich effect (SZE) in high redshift galaxy clusters detected at infrared wavelengths in the Massive and Distant Clusters of WISE Survey (MaDCoWS). Using typical on-source integration times of 3-4 hours per cluster, MUSTANG2 on the Green Bank Telescope was able to measure strong detections of SZE decrements and statistically significant masses on 14 out of 16 targets. On the remaining two, weaker (3.7 sigma) detections of the SZE signal and strong upper limits on the masses were obtained. In this paper we present masses and pressure profiles of each target and outline the data analysis used to recover these quantities. Of the clusters with strong detections, three show significantly flatter pressure profiles while, from the MUSTANG2 data, five others show signs of disruption at their cores. However, outside of the cores of the clusters, we were unable to detect significant amounts of asymmetry. Finally, there are indications that the relationship between optical richness used by MaDCoWS and SZE-inferred mass may be significantly flatter than indicated in previous studies., Comment: Accepted by ApJ ; 20 pages, 7 figures

We present nearly 500 days of observations of the tidal disruption event ASASSN-18pg, spanning from 54 days before peak light to 441 days after peak light. Our dataset includes X-ray, UV, and optical photometry, optical spectroscopy, radio observations, and the first published spectropolarimetric observations of a TDE. ASASSN-18pg was discovered on 2018 July 11 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d=78.6$ Mpc, and with a peak UV magnitude of $m\simeq14$ it is both one of the nearest and brightest TDEs discovered to-date. The photometric data allow us to track both the rise to peak and the long-term evolution of the TDE. ASASSN-18pg peaked at a luminosity of $L\simeq2.2\times10^{44}$ erg s$^{-1}$, and its late-time evolution is shallower than a flux $\propto t^{-5/3}$ power-law model, similar to what has been seen in other TDEs. ASASSN-18pg exhibited Balmer lines and spectroscopic features consistent with Bowen fluorescence prior to peak which remained detectable for roughly 225 days after peak. Analysis of the two-component H$\alpha$ profile indicates that, if they are the result of reprocessing of emission from the accretion disk, the different spectroscopic lines may be coming from regions between $\sim10$ and $\sim60$ light-days from the black hole. No X-ray emission is detected from the TDE and there is no evidence of a jet or strong outflow detected in the radio. Our spectropolarimetric observations give no strong evidence for significant asphericity in the emission region, with the emission region having an axis ratio of at least $\sim0.65$., Comment: 22 pages, 11 figures, 5 tables. Submitted to ApJ. A machine-readable table containing the host-subtracted photometry presented in this manuscript is included as an ancillary file

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

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

We constrain the evolution of the brightest cluster galaxy plus intracluster light (BCG+ICL) using an ensemble of 42 galaxy groups and clusters that span redshifts of z = 0.05-1.75 and masses of $M_{500,c}=2\times10^{13}-10^{15}$ M$_\odot$ Specifically, we measure the relationship between the BCG+ICL stellar mass $M_\star$ and $M_{500,c}$ at projected radii 10 < r < 100 kpc for three different epochs. At intermediate redshift (z = 0.40), where we have the best data, we find $M_\star\propto M_{500,c}^{0.48\pm0.06}$. Fixing the exponent of this power law for all redshifts, we constrain the normalization of this relation to be $2.08\pm0.21$ times higher at z = 0.40 than at high redshift (z = 1.55). We find no change in the relation from intermediate to low redshift (z = 0.10). In other words, for fixed $M_{500,c}$, $M_\star$ at 10 < r < 100 kpc increases from z = 1.55 to z = 0.40 and not significantly thereafter. Theoretical models predict that the physical mass growth of the cluster from z = 1.5 to z = 0 within $r_{500,c}$ is a factor of 1.4, excluding evolution due to definition of $r_{500,c}$. We find that $M_\star$ within the central 100 kpc increases by a factor of 3.8 over the same period. Thus, the growth of $M_\star$ in this central region is more than a factor of two greater than the physical mass growth of the cluster as a whole. Furthermore, the concentration of the BCG+ICL stellar mass, defined by the ratio of stellar mass within 10 kpc to the total stellar mass within 100 kpc, decreases with increasing $M_{500,c}$ at all redshift. We interpret this result as evidence for inside-out growth of the BCG+ICL over the past ten Gyrs, with stellar mass assembly occuring at larger radii at later times., Comment: Accepted to MNRAS, 9 pages, 6 figures