Smooth Accretion Research Articles

JWST/NIRCam Paβ Narrowband Imaging Reveals Ordinary Dust Extinction for Hα Emitters within the Spiderweb Protocluster at z = 2.16

We combine JWST/NIRCam and Subaru/MOIRCS dual Paβ+Hα narrowband imaging to trace the dust attenuation and the star formation activities of a sample of 43 Hα emitters at the core of one of the most massive and best-studied clusters in formation at the cosmic noon: the Spiderweb protocluster at z = 2.16. We find that most Hα emitters display Paβ/Hα ratios compatible with Case B recombination conditions, which translates into nebular extinction values ranging at AV ≈ 0–3 mag, and dust corrected Paβ star formation rates consistent with coeval main sequence field galaxies at fixed stellar mass ( 9.4<logM*/M⊙<11.0 ) during this cosmic epoch. Furthermore, we investigate possible environmental impacts on dust extinction across the protocluster large-scale structure and find no correlation between the dustiness of its members and environmental proxies such as phase-space position, clustercentric radius, or local density. These results support the scenario for which dust production within the main galaxy population of this protocluster is driven by secular star formation activities fueled by smooth gas accretion across its large-scale structure. This downplays the role of gravitational interactions in boosting star formation and dust production within the Spiderweb protocluster, in contrast with observations in higher redshift and less evolved protocluster cores.

Cosmic Sands: The Origin of Dusty, Star-forming Galaxies in the Epoch of Reionization

We present the Cosmic Sands suite of cosmological zoom-in simulations based on the simba galaxy formation model in order to study the buildup of the first massive and dusty galaxies in the early universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the buildup of stellar mass. The galaxies are already appreciably chemically evolved by z = 7, with modeled dust masses comparable to those inferred from observations in the same epoch, except for the most extreme systems. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250 M ⊙ yr−1, but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust RT, we find that, while the infrared luminosities of the most-dust-rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical z = 2 submillimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies.

Gas metallicity distributions in SDSS-IV MaNGA galaxies: what drives gradients and local trends?

ABSTRACTThe gas metallicity distributions across individual galaxies and across galaxy samples can teach us much about how galaxies evolve. Massive galaxies typically possess negative metallicity gradients, and mass and metallicity are tightly correlated on local scales over a wide range of galaxy masses; however, the precise origins of such trends remain elusive. Here, we employ data from SDSS-IV MaNGA to explore how gas metallicity depends on the local stellar mass density and on galactocentric radius within individual galaxies. We also consider how the strengths of these dependencies vary across the galaxy mass-size plane. We find that radius is more predictive of local metallicity than stellar mass density in extended lower-mass galaxies, while we find density and radius to be almost equally predictive in higher-mass and more compact galaxies. Consistent with previous work, we find a mild connection between metallicity gradients and large-scale environment; however, this is insufficient to explain variations in gas metallicity behaviour across the mass-size plane. We argue our results to be consistent with a scenario in which extended galaxies have experienced smooth gas accretion histories, producing negative metallicity gradients over time. We further argue that more compact and more massive systems have experienced increased merging activity that disrupts this process, leading to flatter metallicity gradients and more dominant density-metallicity correlations within individual galaxies.

The buildup of galaxies and their spheroids: The contributions of mergers, disc instabilities, and star formation

ABSTRACT We use the GALFORM semi-analytical model of galaxy formation and the Planck-Millennium simulation to investigate the origins of stellar mass in galaxies and their spheroids. We compare the importance of mergers and disc instabilities, as well as the starbursts that they trigger. We find that the fraction of galaxy stellar mass formed ex situ (i.e. through mergers; fex) increases sharply from M* = 1011 M⊙ upwards, reaching 80 per cent at M* = 1011.3 M⊙. The massive end of the fex–M* relation does not evolve with redshift, in disagreement with other models. For low-mass galaxies we find larger ex situ contributions at z = 0 than in other models (7–12 per cent), with a decrease towards higher redshifts. Major mergers contribute roughly half of the ex situ mass, with minor mergers and smooth accretion of satellites both accounting for ≈25 per cent, almost independent of stellar mass and redshift. Mergers dominate in building up high-mass (M*, sph &amp;gt; 1011 M⊙) and low-mass (M*, sph &amp;lt; 108.5 M⊙) spheroids. Disc instabilities and their associated starbursts dominate for intermediate-mass spheroids (108.5 &amp;lt; M*, sph &amp;lt; 1011 M⊙) at z = 0. The mass regime where pseudo-bulges dominate is in agreement with observed pseudo-bulge fractions, but the peak value in the pseudo-bulge fraction predicted by GALFORM is likely too high. Starbursts induced by disc instabilities are the dominant channel for spheroid growth at all redshifts, while merger-induced starbursts are relatively negligible, except at very high redshifts (z &amp;gt; 5).

Study of general relativistic magnetohydrodynamic accretion flow around black holes

ABSTRACT We present a novel approach to study the global structure of steady, axisymmetric, advective, magnetohydrodynamic (MHD) accretion flow around black holes in full general relativity (GR). Considering ideal MHD conditions and relativistic equation of state (REoS), we solve the governing equations to obtain all possible smooth global accretion solutions. We examine the dynamical and thermodynamical properties of accreting matter in terms of the flow parameters, namely energy (${\cal E}$), angular momentum (${\cal L}$), and local magnetic fields. For a vertically integrated GRMHD flow, we observe that toroidal component (bϕ) of the magnetic fields generally dominates over radial component (br) at the disc equatorial plane. This evidently suggests that toroidal magnetic field indeed plays important role in regulating the disc dynamics. We further notice that the disc remains mostly gas pressure (pgas) dominated (β = pgas/pmag &amp;gt; 1, pmag refers magnetic pressure) except at the near horizon region, where magnetic fields become indispensable (β ∼ 1). We observe that Maxwell stress is developed that eventually yields angular momentum transport inside the disc. Towards this, we calculate the viscosity parameter (α) that appears to be radially varying. In addition, we examine the underlying scaling relation between α and β, which clearly distinguishes two domains coexisted along the radial extent of the disc. Finally, we discuss the utility of the present formalism in the realm of GRMHD simulation studies.

Star Formation and Quenching of Central Galaxies from Stacked Hi Measurements

Abstract We quantitatively investigate the dependence of central galaxy Hi mass (M Hi ) on the stellar mass (M *), halo mass (M h ), star formation rate (SFR), and central stellar surface density within 1 kpc (Σ1), taking advantage of the Hi spectra stacking technique using both the Arecibo Fast Legacy ALFA Survey and the Sloan Digital Sky Survey. We find that the shapes of M Hi –M h and M Hi –M * relations are remarkably similar for both star-forming and quenched galaxies, with massive quenched galaxies having constantly lower Hi masses of around 0.6 dex. This similarity strongly suggests that neither halo mass nor stellar mass is the direct cause of quenching, but rather the depletion of the Hi reservoir. While the Hi reservoir for low-mass galaxies of M * &lt; 1010.5 M ☉ strongly increases with M h , more massive galaxies show no significant dependence of M Hi with M h , indicating that the main effect of halo is to determine the smooth cold gas accretion. We find that the star formation and quenching of central galaxies are directly regulated by the available Hi reservoir, with an average relation of SFR ∝ M H I 2.75 / M * 0.40 , implying a quasi-steady state of star formation. We further confirm that galaxies are depleted of their Hi reservoir once they drop off the star formation main sequence and there is a very tight and consistent correlation between M Hi and Σ1 in this phase, with M H I ∝ Σ 1 − 2 . This result is consistent with the compaction-triggered quenching scenario, with galaxies going through three evolutionary phases of cold gas accretion, compaction and post-compaction, and quenching.

Globular cluster numbers in dark matter haloes in a dual formation scenario: an empirical model within emerge

ABSTRACT We present an empirical model for the number of globular clusters (GCs) in galaxies based on recent data showing a tight relationship between dark matter halo virial masses and GC numbers. While a simple base model forming GCs in low-mass haloes reproduces this relation, we show that a second formation pathway for GCs is needed to account for observed younger GC populations. We confirm previous works that reported the observed linear correlation as being a consequence of hierarchical merging and its insensitivity to the exact GC formation processes at higher virial masses, even for a dual formation scenario. We find that the scatter of the linear relation is strongly correlated with the relative amount of smooth accretion: the more dark matter is smoothly accreted, the fewer GCs a halo has compared to other haloes of the same mass. This scatter is smaller than that introduced by halo mass measurements, indicating that the number of GCs in a galaxy is a good tracer for its dark matter mass. Smooth accretion is also the reason for a lower average dark matter mass per GC in low-mass haloes. Finally, we successfully reproduce the observed general trend of GCs being old and the tendency of more massive haloes hosting older GC systems. Including the second GC formation mechanism through gas-rich mergers leads to a more realistic variety of GC age distributions and also introduces an age inversion in the halo virial mass range log Mvir/M⊙ = 11–13.

Green valley galaxies in the cosmic web: internal versus environmental quenching

We analyze the SDSS data to classify the galaxies based on their colour using a fuzzy set-theoretic method and quantify their environments using the local dimension. We find that the fraction of the green galaxies does not depend on the environment and 10%–20% of the galaxies at each environment are in the green valley depending on the stellar mass range chosen. Approximately 10% of the green galaxies at each environment host an AGN. Combining data from the Galaxy Zoo, we find that ∼ 95% of the green galaxies are spirals and ∼ 5% are ellipticals at each environment. Only ∼ 8% of green galaxies exhibit signs of interactions and mergers, ∼ 1% have dominant bulge, and ∼ 6% host a bar. We show that the stellar mass distributions for the red and green galaxies are quite similar at each environment. Our analysis suggests that the majority of the green galaxies must curtail their star formation using physical mechanism(s) other than interactions, mergers, and those driven by bulge, bar and AGN activity. We speculate that these are the massive galaxies that have grown only via smooth accretion and suppressed the star formation primarily through mass driven quenching. Using a Kolmogorov-Smirnov test, we do not find any statistically significant difference between the properties of green galaxies in different environments. We conclude that the environmental factors play a minor role and the internal processes play the dominant role in quenching star formation in the green valley galaxies.

Evolution of splashback boundaries and gaseous outskirts: insights from mergers of self-similar galaxy clusters

ABSTRACT A self-similar spherical collapse model predicts a dark matter (DM) splashback and accretion shock in the outskirts of galaxy clusters while missing a key ingredient of structure formation – processes associated with mergers. To fill this gap, we perform simulations of merging self-similar clusters and investigate their DM and gas evolution in an idealized cosmological context. Our simulations show that the cluster rapidly contracts during the major merger and the splashback radius rsp decreases, approaching the virial radius rvir. While in the self-similar model rsp depends on a smooth mass accretion rate parameter Γs, our simulations show that in the presence of mergers, rsp responds to the changes in the total mass accretion rate Γvir, which accounts for both mergers and smooth accretion. The scatter of the Γvir − rsp/rvir relation indicates a generally low Γs ∼ 1 in clusters in cosmological simulations. In contrast to the DM, the hot gaseous atmospheres significantly expand by the merger-accelerated (MA-) shocks formed when the runaway merger shocks overtake the outer accretion shock. After a major merger, the MA-shock radius is larger than rsp by a factor of up to ∼1.7 for Γs ≲ 1 and is ∼rsp for Γs ≳ 3. This implies that (1) mergers could easily generate the MA-shock-splashback offset measured in cosmological simulations, and (2) the smooth mass accretion rate is small in regions away from filaments where MA-shocks reside. We further discuss the shapes of the DM haloes, various shocks, and contact discontinuities formed at different epochs of the merger, and the ram-pressure stripping in cluster outskirts.

Star–Gas Misalignment in Galaxies. II. Origins Found from the Horizon-AGN Simulation

Abstract There have been many studies aiming to reveal the origins of the star–gas misalignment found in galaxies, but there still is a lack of understanding of the contribution from each formation channel candidate. We aim to answer the question by investigating the misaligned galaxies in the Horizon-AGN simulation. There are 27,903 galaxies of stellar mass M * &gt; 1010 M ⊙ in our sample, of which 5984 are in a group in the halo mass of M 200 &gt; 1012 M ⊙. We have identified four main formation channels of misalignment and quantified their levels of contribution: mergers (35%), interaction with nearby galaxies (23%), interaction with dense environments or their central galaxies (21%), and secular evolution, including smooth accretion from neighboring filaments (21%). We found in the simulation that the gas, rather than stars, is typically more vulnerable to dynamical disturbances; hence, misalignment formation is mainly due to the change in the rotational axis of the gas rather than stars, regardless of the origin. We have also inspected the lifetime (duration) of the misalignment. The decay timescale of the misalignment shows a strong anticorrelation with the kinematic morphology (V/σ) and the cold gas fraction of the galaxy. The misalignment has a longer lifetime in denser regions, which is linked with the environmental impact on the host galaxy. There is a substantial difference in the length of the misalignment lifetime depending on the origin, and it can be explained by the magnitude of the initial position angle offset and the physical properties of the galaxies.

Evolution of subhalo orbits in a smoothly growing host halo potential

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

FIR-luminous [C ii] Emitters in the ALMA-SCUBA-2 COSMOS Survey (AS2COSMOS): The Nature of Submillimeter Galaxies in a 10 Comoving Megaparsec-scale Structure at z ∼ 4.6

Abstract We report the discovery of a 10 comoving megaparsec (cMpc)-scale structure traced by massive submillimeter galaxies (SMGs) at z ∼ 4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources (S 850μm ≥ 6.2 mJy) across 1.6 degrees2 in the COSMOS field. We identify four [C ii] emitting SMGs and two probable [C ii] emitting SMG candidates at z = 4.60–4.64 with velocity-integrated signal-to-noise ratio of S/N &gt; 8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330–550 km s−1 and the inferred host dark-matter halo masses are ∼2–8 × 1012 M ⊙. From their estimated halo masses and [C ii] luminosity function, we suggest that these galaxies have a high (50%–100%) duty cycle and high (∼0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute ≃2% to the total star formation rate density at z = 4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ∼1015 M ⊙ cluster at z = 0. This survey reveals a synchronized dusty starburst in massive halos at z &gt; 4, which could be driven by mergers or fed by smooth gas accretion.

The distribution and properties of DLAs at z ≤ 2 in the EAGLE simulations

ABSTRACT Determining the spatial distribution and intrinsic physical properties of neutral hydrogen on cosmological scales is one of the key goals of next-generation radio surveys. We use the EAGLE galaxy formation simulations to assess the properties of damped Lyman α absorbers (DLAs) that are associated with galaxies and their underlying dark matter haloes between 0 ≤ z ≤ 2. We find that the covering fraction of DLAs increases at higher redshift; a significant fraction of neutral atomic hydrogen (H i) resides in the outskirts of galaxies with stellar mass ≥1010 M⊙; and the covering fraction of DLAs in the circumgalactic medium (CGM) is enhanced relative to that of the interstellar medium (ISM) with increasing halo mass. Moreover, we find that the mean density of the H i in galaxies increases with increasing stellar mass, while the DLAs in high- and low-halo mass systems have higher column densities than those in galaxies with intermediate halo masses (∼1012 M⊙ at z = 0). These high-impact CGM DLAs in high-stellar mass systems tend to be metal poor, likely tracing smooth accretion. Overall, our results point to the CGM playing an important role in DLA studies at high redshift (z ≥ 1). However, their properties are impacted both by numerical resolution and the detailed feedback prescriptions employed in cosmological simulations, particularly that of active galactic nuclei.

On the accretion history of galaxy clusters: temporal and spatial distribution

ABSTRACT We analyse the results of an Eulerian adaptive mesh refinement cosmological simulation in order to quantify the mass growth of galaxy clusters, exploring the differences between dark matter and baryons. We have determined the mass assembly histories (MAHs) of each of the mass components and computed several proxies for the instantaneous mass accretion rate (MAR). The mass growth of both components is clearly dominated by the contribution of major mergers, but high MARs can also occur during smooth accretion periods. We explored the correlations between MARs, merger events, and clusters’ environments, finding the mean densities in 1 ≤ r/R200m ≤ 1.5 to correlate strongly with Γ200m in massive clusters that undergo major mergers through their MAH. From the study of the dark matter velocity profiles, we find a strong anticorrelation between the MAR proxies Γ200m and α200m. Last, we present a novel approach to study the angularly resolved distribution of gas accretion flows in simulations, which allows to extract and interpret the main contributions to the accretion picture and to assess systematic differences between the thermodynamical properties of each of these contributions using multipolar analysis. We have preliminarily applied the method to the best numerically resolved cluster in our simulation. Amongst the most remarkable results, we find that the gas infalling through the cosmic filaments has systematically lower entropy compared to the isotropic component, but we do not find a clear distinction in temperature.

New empirical constraints on the cosmological evolution of gas and stars in galaxies

ABSTRACT We combine the latest observationally motivated constraints on stellar properties in dark matter haloes, along with data-driven predictions for the atomic (H I) and molecular (H2) gas evolution in galaxies, to derive empirical relationships between the build-up of galactic components and their evolution over cosmic time. At high redshift (z ≳ 4), the frameworks imply that galaxies acquire their cold gas (both atomic and molecular) mostly by accretion, with the fraction of cold gas reaching about 20 per cent of the cosmic baryon fraction. We infer a strong dependence of the star formation rate on the H2 mass, suggesting a near-universal depletion time-scale of 0.1–1 Gyr in Milky Way-sized haloes (of masses 1012 M⊙ at z = 0). There is also evidence for a near-universality of the Kennicutt–Schmidt relation across redshifts, with very little dependence on stellar mass, if a constant conversion factor (αCO) of CO luminosity to molecular gas mass is assumed. Combining the atomic and molecular gas observations with the stellar build-up illustrates that galactic mass assembly in Milky Way-sized haloes proceeds from smooth accretion at high redshifts towards becoming merger-dominated at late times (z ≲ 0.6). Our results can be used to constrain numerical simulations of the dominant growth and accretion processes of galaxies over cosmic history.

Where did the globular clusters of the Milky Way form? Insights from the E-MOSAICS simulations

ABSTRACT Globular clusters (GCs) are typically old, with most having formed at z ≳ 2. This makes understanding their birth environments difficult, as they are typically too distant to observe with sufficient angular resolution to resolve GC birth sites. Using 25 cosmological zoom-in simulations of Milky Way-like galaxies from the E-MOSAICS project, with physically motivated models for star formation, feedback, and the formation, evolution, and disruption of GCs, we identify the birth environments of present-day GCs. We find roughly half of GCs in these galaxies formed in situ (52.0 ± 1.0 per cent) between z ≈ 2–4, in turbulent, high-pressure discs fed by gas that was accreted without ever being strongly heated through a virial shock or feedback. A minority of GCs form during mergers (12.6 ± 0.6 per cent in major mergers, and 7.2 ± 0.5 per cent in minor mergers), but we find that mergers are important for preserving the GCs seen today by ejecting them from their natal, high density interstellar medium (ISM), where proto-GCs are rapidly destroyed due to tidal shocks from ISM substructure. This chaotic history of hierarchical galaxy assembly acts to mix the spatial and kinematic distribution of GCs formed through different channels, making it difficult to use observable GC properties to distinguish GCs formed in mergers from ones formed by smooth accretion, and similarly GCs formed in situ from those formed ex situ. These results suggest a simple picture of GC formation, in which GCs are a natural outcome of normal star formation in the typical, gas-rich galaxies that are the progenitors of present-day galaxies.

Dynamical heating of the X-ray emitting intracluster medium: the roles of merger shocks and turbulence dissipation

ABSTRACT The diffuse plasma inside clusters of galaxies has X-ray emitting temperatures of a few keV. The physical mechanisms that heat this intracluster medium (ICM) to such temperatures include the accretion shock at the periphery of a galaxy cluster, the shocks driven by merger events, as well as a somewhat overlooked mechanism – the dissipation of intracluster turbulent motions. We study the relative role of these heating mechanisms using galaxy clusters in Lagrangian tracer particle re-simulations of the Omega500 cosmological simulation. We adopt a novel analysis method of decomposing the temperature increase at each time-step into the contribution from dissipative heating and that from adiabatic heating. In the high-resolution spatial–temporal map of these heating rates, merger tracks are clearly visible, demonstrating the dominant role of merger events in heating the ICM. The dissipative heating contributed by each merger event is extended in time and also occurs in the rarefaction regions, suggesting the importance of heating by the dissipation of merger-induced turbulence. Quantitative analysis shows that turbulence heating, rather than direct heating at merger shocks, dominates the temperature increase of the ICM especially at inner radii r &amp;lt; r500c. In addition, we find that many merger shocks can propagate with almost constant velocity to very large radii r ≫ r500c, some even reach and join with the accretion shock and becoming the outer boundary of the ICM. Altogether, these results suggest that the ICM is heated more in an ‘inside–out’ fashion rather than ‘outside–in’ as depicted in the classical smooth accretion picture.

Direct Measurement of the H i-halo Mass Relation through Stacking

Abstract We present accurate measurements of the total H i mass in dark matter halos of different masses at z ∼ 0, by stacking the H i spectra of entire groups from the Arecibo Fast Legacy ALFA Survey. The halos are selected from the optical galaxy group catalog constructed from the Sloan Digital Sky Survey DR7 Main Galaxy sample, with reliable measurements of halo mass and halo membership. We find that the H i-halo mass relation is not a simple monotonic function, as assumed in several theoretical models. In addition to the dependence of halo mass, the total H i gas mass shows a strong dependence on the halo richness, with larger H i masses in groups with more members at fixed halo masses. Moreover, halos with at least three member galaxies in the group catalog have a sharp decrease of the H i mass, potentially caused by the virial halo shock-heating and the active galactic nucleus (AGN) feedback. The dominant contribution of the H i gas comes from the central galaxies for halos of , while the satellite galaxies dominate over more massive halos. Our measurements are consistent with a three-phase formation scenario of the H i-rich galaxies. The smooth cold gas accretion is driving the H i mass growth in halos of , with late-forming halos having more H i accreted. The virial halo shock-heating and AGN feedback will take effect to reduce the H i supply in halos of . The H i mass in halos more massive than generally grows by mergers, with the dependence on halo richness becoming much weaker.

Growing a ‘cosmic beast’: observations and simulations of MACS J0717.5+3745

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

LoCuSS: The infall of X-ray groups on to massive clusters

Galaxy clusters are expected to form hierarchically in a LCDM universe, growing primarily through mergers with lower mass clusters and the continual accretion of group-mass halos. Galaxy clusters assemble late, doubling their masses since z~0.5, and so the outer regions of clusters should be replete with infalling group-mass systems. We present an XMM-Newton survey to search for X-ray groups in the infall regions of 23 massive galaxy clusters at z~0.2, identifying 39 X-ray groups that have been spectroscopically confirmed to lie at the cluster redshift. These groups have mass estimates in the range 2x10^13-7x10^14Msun, and group-to-cluster mass ratios as low as 0.02. The comoving number density of X-ray groups in the infall regions is ~25x higher than that seen for isolated X-ray groups from the XXL survey. The average mass per cluster contained within these X-ray groups is 2.2x10^14Msun, or 19% of the mass within the primary cluster itself. We estimate that ~10^15Msun clusters increase their masses by 16% between z=0.223 and the present day due to the accretion of groups with M200>10^13.2Msun. This represents about half of the expected mass growth rate of clusters at these late epochs. The other half is likely to come from smooth accretion of matter not bound in halos. The mass function of the infalling X-ray groups appears significantly top-heavy with respect to that of field X-ray systems, consistent with expectations from numerical simulations, and the basic consequences of collapsed massive dark matter halos being biased tracers of the underlying large-scale density distribution.