Cluster Member Galaxies Research Articles

ViCTORIA project: The LOFAR-MeerKAT view of active galactic nuclei in Virgo cluster early-type galaxies

Context. The evolution of active galactic nuclei (AGNs) is closely connected to their host galaxies and surroundings. Via feedback processes, AGNs can counteract the cooling of the intracluster medium (ICM) and suppress star formation in their host galaxies. Radio observations at low frequencies provide a glimpse into the history of AGN activity. The Virgo cluster is a substantial reservoir of nearby galaxies and provides an ideal laboratory for the study of AGNs as well as their feedback mechanisms. Aims. The aim of our work is to characterise the AGN population within the Virgo cluster down to low radio luminosities, constrain the AGN duty cycle, and investigate environmental feedback in cluster member galaxies. Methods. We analysed 144 MHz and 1.3 GHz radio observations of early-type galaxies from the ACS Virgo Cluster Survey (ACSVCS) taken with LOFAR and MeerKAT. Results. We detect 12 of these galaxies at 144 MHz, five of which show clearly extended radio emission. The radio luminosity shows a strong dependence on the stellar mass of the host galaxy, in agreement with previous results. As a notable outlier, the massive elliptical galaxy NGC 4365 (M* = 2.2 × 1011 M⊙) is not detected as compact source in the LOFAR observations. Instead, it is surrounded by diffuse, low-surface brightness emission, which hints towards a past phase of stronger nuclear activity. Furthermore, we find a cavity in NGC 4472 (=M 49) inflated by the wide-angle tail only visible in the LOFAR data, which implies that the cavity was created by a past outburst. The corresponding cavity power is of the same order of magnitude as the jet power in the present duty cycle of the AGN.

Subhalos in Galaxy Clusters: Coherent Accretion and Internal Orbits

Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately 80° between a ∼ 0.1 and the present day. We identify coherent regions in the angular distribution of subhalo accretion, and ∼68% of accreted subhalos enter their host halo through ∼38% of the surface area at the virial radius. The majority of galaxy clusters in the sample have ∼2 such coherent regions. We further measure angular orbits of subhalos with respect to the host major axis and use a clustering algorithm to identify distinct orbit modes with varying oscillation timescales. The orbit modes correlate with subhalo accretion conditions. Subhalos in orbit modes with shorter oscillations tend to have lower peak masses and accretion directions somewhat more aligned with the major axis. One orbit mode, exhibiting the least oscillatory behavior, largely consists of subhalos that accrete near the plane perpendicular to the host halo major axis. Our findings are consistent with expectations from inflow from major filament structures and internal dynamical friction: most subhalos accrete through coherent regions, and more massive subhalos experience fewer orbits after accretion. Our work offers a unique quantification of subhalo dynamics that can be connected to how the intracluster medium strips and quenches cluster galaxies.

THE THREE HUNDRED: Contrasting cluster galaxy density in hydrodynamical and dark matter simulations

Cluster number counts at visible and IR wavelengths will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, the performance of cluster detection algorithms, which at these wavelengths are sensitive to the spatial distributions of the cluster galaxy members and their luminosity functions, need to be accurately characterized. Using THE THREE HUNDRED hydrodynamical and dark-matter-only simulations, we studied a complete sample of massive clusters beyond 7 (5) × 1014 M⊙ at redshift 0 (1) on a (1.48 Gpc)3 volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 × 109 M⊙) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark-matter-only version of these simulations, we demonstrate that baryonic physics preserves significantly low-mass subhalos (galaxies), as has also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark-matter-only simulations of the same objects and taking the same limit in subhalo mass, we find galaxy density profiles that are significantly more cuspy towards the centre of the clusters, where the low-mass substructures tend to concentrate. We conclude that using a dark-matter-only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high-resolution hydro simulations we conclude that a mass resolution of 1.8 × 108h−1 M⊙ will be needed for THE THREE HUNDRED simulations to approach the expected magnitude limits for the Euclid survey. These simulations are currently under way.

ICM-SHOX. I. Methodology Overview and Discovery of a Gas–Dark Matter Velocity Decoupling in the MACS J0018.5+1626 Merger

Galaxy cluster mergers are rich sources of information to test cluster astrophysics and cosmology. However, cluster mergers produce complex projected signals that are difficult to interpret physically from individual observational probes. Multi-probe constraints on the gas and dark matter (DM) cluster components are necessary to infer merger parameters that are otherwise degenerate. We present Improved Constraints on Mergers with SZ, Hydrodynamical simulations, Optical, and X-ray (ICM-SHOX), a systematic framework to jointly infer multiple merger parameters quantitatively via a pipeline that directly compares a novel combination of multi-probe observables to mock observables derived from hydrodynamical simulations. We report a first application of the ICM-SHOX pipeline to MACS J0018.5+1626, wherein we systematically examine simulated snapshots characterized by a wide range of initial parameters to constrain the MACS J0018.5+1626 merger geometry. We constrain the epoch of MACS J0018.5+1626 to the range 0–60 Myr post-pericenter passage, and the viewing angle is inclined ≈27°–40° from the merger axis. We obtain constraints for the impact parameter (≲250 kpc), mass ratio (≈1.5–3.0), and initial relative velocity when the clusters are separated by 3 Mpc (≈1700–3000 km s−1). The primary and secondary clusters initially (at 3 Mpc) have gas distributions that are moderately and strongly disturbed, respectively. We discover a velocity space decoupling of the DM and gas distributions in MACS J0018.5+1626, traced by cluster-member galaxy velocities and the kinematic Sunyaev–Zel'dovich effect, respectively. Our simulations indicate this decoupling is dependent on the different collisional properties of the two distributions for particular merger epochs, geometries, and viewing angles.

The powerful lens galaxy cluster PLCK G287.0+32.9 (θE ∼ 43″)

Aims. We present a new high-precision strong-lensing model of PLCK G287.0+32.9, a massive lens galaxy cluster at z = 0.383, with the aim of obtaining an accurate estimation of its effective Einstein radius and total mass distribution. We also present a spectroscopic catalog containing accurate redshift measurements for close to 500 objects up to redshift z = 6, including multiply lensed sources and cluster member galaxies. Methods. We exploited high-quality spectroscopic data from the Multi Unit Spectroscopic Explorer (MUSE), covering a central 3 arcmin2 region of the cluster. We supplemented the spectroscopic catalog by including redshift measurements from VIsible MultiObject Spectrograph (VIMOS) and DEep Imaging Multi-Object Spectrograph (DEIMOS). We identified 129 spectroscopic cluster member galaxies with redshift values of 0.360 ≤ z ≤ 0.405, and mF160W ≤ 21. We complemented this galaxy cluster member sample with 24 photometric members identified with a convolutional neural network (CNN) approach. We also identified 114 multiple images from 28 background sources, of which 84 images from 16 sources are new and the remaining ones have already been identified in previous works. From these, we extracted “golden sample” of 47 secure multiple images and used them, together with the selected cluster member, to build and optimize several strong-lensing models with the software lenstool. Results. The best-fitting lens model shows a root mean square (RMS) separation value between the predicted and observed positions of the multiple images of 0.″75. Using its predictive power, we found three new multiple images and we confirm the configuration of three systems of multiple images that were not used for the optimization of the model. For a source at a redshift of zs = 2, we found a cluster with an Einstein radius of θE = 43.4″ ± 0.1″. This value is in agreement with previous estimates and corresponds to a total mass enclosed in the critical curve of ME = 3.33−0.07+0.02 × 1014 M⊙. Conclusions. The combined application of ancillary Hubble Space Telescope (HST) imaging, VIMOS and DEIMOS data, and the new MUSE spectroscopic observations allowed us to build a new lens model of the galaxy cluster PLCK G287.0+32.9, with an improvement in terms of reconstructing the observed positions of the multiple images of a factor of 2.5 with respect to previous models. The derived total mass distribution confirms this cluster to be a very prominent gravitational lens, with an effective Einstein radius of θE ∼ 43″. We were also able to construct an extensive spectroscopic catalog containing 490 objects, of which 153 are bright cluster members with mF160W ≤ 21, and 114 are multiple images.

A Deep Redshift Survey of the Perseus Cluster (A426): Spatial Distribution and Kinematics of Galaxies

We study the global kinematics of the Perseus galaxy cluster (A426) at redshift z = 0.017 using a large sample of galaxies from our new MMT/Hectospec spectroscopic observation for this cluster. The sample includes 1447 galaxies with measured redshifts within 60′ from the cluster center (1148 from this MMT/Hectospec program and 299 from the literature). The resulting spectroscopic completeness is 67% at r-band apparent magnitude r Petro,0 ≤ 18.0 within 60′ from the cluster center. To identify cluster member galaxies in this sample, we develop a new open-source Python package, CausticSNUpy. This code implements the algorithm of the caustic technique and yields 418 member galaxies within 60′ of the cluster. We study the cluster using this sample of member galaxies. The cluster shows no significant signal of global rotation. A statistical test shows that the cluster does not have a noticeable substructure within 30′. We find two central regions where the X-ray-emitting intracluster medium and galaxies show significant velocity differences (>7σ). On a large scale, however, the overall morphology and kinematics between the intracluster medium and galaxies agree well. Our results suggest that the Perseus cluster is a relaxed system and has not experienced a recent merger.

Distribution of Merging and Post-merger Galaxies in Nearby Galaxy Clusters

We study the incidence and spatial distribution of galaxies that are currently undergoing gravitational merging (M) or that have signs of being post-merger (PM) in six galaxy clusters (A754, A2399, A2670, A3558, A3562, and A3716) within the redshift range of 0.05 ≲ z ≲0.08. To this aim, we obtained Dark Energy Camera mosaics in the u′ , g′ , and r′ bands covering up to 3 × R 200 of the clusters, reaching 28 mag arcsec−2 surface brightness limits. We visually inspect u′ g′ r′ color-composite images of volume-limited (M r < −20) cluster member galaxies to identify whether galaxies are of M or PM type. We find 4% M-type and 7% PM-type galaxies in the galaxy clusters studied. By adding spectroscopic data and studying the projected phase-space diagram (PPSD) of the projected clustocentric radius and the line-of-sight velocity, we find that PM-type galaxies are more virialized than M-type galaxies, having a 1%–5% higher fraction within the escape-velocity region, while the fraction of M-type was ∼10% higher than the PM type in the intermediate environment. Similarly, on a substructure analysis, M types were found in groups in the outskirts, while PM-type populated groups were found in ubiquitous regions of the PPSD. Adopting literature-derived dynamical state indicator values, we observed a higher abundance of M types in dynamically relaxed clusters. This finding suggests that galaxies displaying post-merger features within clusters likely merged in low-velocity environments, including clusters outskirts and dynamically relaxed clusters.

The Optical Properties of Galaxy Cluster Abell 2319

The optical photometric and spectroscopic investigations of the massive and merging galaxy cluster Abell 2319 (A2319) are presented here. RTT150 telescope of TÜBİTAK, Antalya, Türkiye used CCD imaging and spectroscopic observations. In this paper, 110 galaxies were determined in A2319 and defined as the magnitudes of the Bessel B and R filters in each cluster member galaxy. Spectral observations were done of the brightest cluster galaxy (BCG) and six additional brilliant galaxies. We estimated the Luminosity Function (LF) of galaxies for each filter. The resulting LF of cluster galaxies for each filter is well-fitted by the Double Schechter function. The best-fit parameter values derived as the characteristic absolute magnitudes are -21.08 ± 0.03, and -20.84 ± 0.01, -21.43 ± 0.02, and -20.54 ± 0.02, and the slopes at the faint end of the LF were -1.34 ± 0.04 and -1.12 ± 0.03, -1.47 ± 0.05 and -1.18 ± 0.03 for B and R filters, respectively.

Augmenting the power of time-delay cosmography in lens galaxy clusters by probing their member galaxies

We present a novel approach to measuring the expansion rate and the geometry of the Universe, which combines time-delay cosmography in lens galaxy clusters with pure samples of ‘cosmic chronometers’ by probing the member galaxies. The former makes use of the measured time delays between the multiple images of time-varying sources strongly lensed by galaxy clusters, while the latter exploits the most massive and passive cluster member galaxies to measure the differential time evolution of the Universe. We applied two different statistical techniques, adopting realistic errors on the measured quantities, to assess the accuracy and the gain in precision on the values of the cosmological parameters. We demonstrate that the proposed combined method allows for a robust and accurate measurement of the value of the Hubble constant. In addition, this provides valuable information on the other cosmological parameters thanks to the complementarity between the two different probes in breaking parameter degeneracies. Finally, we showcased the immediate observational feasibility of the proposed joint method by taking advantage of the existing high-quality spectro-photometric data for several lens galaxy clusters.

Free-form and hybrid lens models for SDSS J1004+4112: substructure and central image time delay constraints

ABSTRACT SDSS J1004+4112 is a well-studied gravitational lens with a recently measured time delay between its first and fourth arriving quasar images. Using this new constraint, we present updated free-form lens reconstructions using the lens inversion method grale, which only uses multiple image and time delay data as inputs. In addition, we obtain hybrid lens reconstructions by including a model of the brightest cluster galaxy as a Sersic lens. For both reconstructions, we use two sets of images as input: one with all identified images, and the other a revised set leaving out images that have been potentially misidentified. We also develop a source position optimization Markov chain Monte Carlo (MCMC) routine, performed on completed grale runs, that allows each model to better match observed image positions and time delays. All the reconstructions produce similar mass distributions, with the hybrid models finding a steeper profile in the centre. Similarly, all the mass distributions are fitted by the Navarro–Frenk–White profile, finding results consistent with previous parametric reconstructions and those derived from Chandra X-ray observations. We identify an ∼5 × 1011 M⊙ substructure apparently unaffiliated with any cluster member galaxy and present in all our models, and study its reality. Using our free-form and hybrid models, we predict a central quasar image time delay of ∼2980 ± 270 and ∼3280 ± 215 d, respectively. A potential future measurement of this time delay will, while being an observational challenge, further constrain the steepness of the central density profile.

The JWST view of the barred galaxy population in the SMACS0723 galaxy cluster

Context. The cosmic evolution of the barred galaxy population provides key information about the secular evolution of galaxies and the settling of rotationally dominated discs. Aims. We study the bar fraction in the SMACSJ0723.37323 (SMACS0723) cluster of galaxies at z = 0.39 using the Early Release Observations obtained with the NIRCam instrument mounted on the JWST telescope. Methods. We visually inspected all cluster member galaxies using the images from the NIRCam F200W filter. We classified the galaxies into ellipticals and discs and determine the presence of a bar. The cluster member selection was based on a combined method using both the available spectroscopy and the color-magnitude relation. Results. As has previously been found in nearby galaxy samples, we find that the bar fraction distribution of SMACS0723 is a strong function of the galaxy stellar luminosity (or stellar mass). The analogy with local clusters, such as Virgo and Coma, reveals a similar distribution among the three clusters for low-mass galaxies (log(M⋆/M⊙)≤9.5). The comparison with a sample of local galaxies in a field environment shows a remarkable lack of bars in this low-mass regime for the SMACS0723 cluster (and, therefore, in Virgo and Coma) with respect to the field. At high masses (log(M⋆/M⊙)≥10.25), galaxies in SMACS0723 show a slightly lower bar fraction than those in Coma. At these high masses, we find a much larger bar fraction in SMACS0723 than previous works on field galaxies at z ∼ 0.4. Nevertheless, the difference is only marginal when we compare with a sample of well-resolved local field galaxies. Thus, we suggest that the improved capabilities of JWST with respect to HST in terms of both spatial resolution and image depth are responsible for the higher bar fraction we obtained. Conclusions. Our results support a scenario where cluster environment affects the formation of bars in a mass-dependent way. At high masses, the mild increase in the bar fraction of local clusters (Coma) with respect to both SMACS0723 and local field galaxies suggests a weak effect coming from the cluster environment possibly triggering bar formation. On the other hand, low-mass galaxies show the same bar fraction in the three clusters (different redshifts) and a significant drop with respect to field galaxies at z = 0, thus suggesting that: (i) the bar fraction of low-mass galaxies in clusters is not evolving during the last 4 Gyr; and (ii) bar formation is severely inhibited in low-mass galaxies residing in clusters.

An excursion into the core of the cluster lens Abell 1689

ABSTRACT Abell 1689 is a well-studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the most notable feature of free-form methods like grale. We used two different sets of multiple image systems from the available strong lensing data – one containing only the secure systems (107 images) and the other containing all available systems, only excluding some very non-secure systems (151 images). For the very well-constrained central ∼100 kpc region of the cluster, we made a detailed comparison of the grale reconstructed lensing mass and stellar mass retrieved by the Spectral Energy Distribution (SED) fitting software fast++. We found a light-unaccompanied mass peak in this region, whose existence, while tentative, is favoured by the distribution of nearby images that are local maxima in the Fermat potential. However, further tests, using different methodologies are needed to confirm the reality of this feature. If it is shown to be real, this light-unaccompanied mass peak is consistent with dark matter self-interaction cross-section σ ≲ 1 cm2 g−1, while being in tension with larger cross-sections.

VICTORIA project: The LOFAR HBA Virgo Cluster Survey

Context. The Virgo cluster is the nearest (d = 16.5 Mpc) massive (M ≥ 1014 M⊙) galaxy cluster and is therefore a prime target for studying astrophysical processes in dense large-scale environments. In the radio band, we can probe the nonthermal components of the interstellar medium, the intracluster medium (ICM), and of active galactic nuclei (AGN). This allows an investigation of the impact of the environment on the evolution of galaxies and the contribution of AGN to ICM heating. With the ViCTORIA (VIrgo Cluster multi-Telescope Observations in Radio of Interacting galaxies and AGN) project, we are carrying out multiple wide-field surveys of the Virgo cluster at different frequencies. Aims. We aim to investigate the impact of the environment on the evolution of galaxies and the contribution of AGN to ICM heating – from the inner cluster regions out to beyond the virial radius. Methods. We performed a survey of the cluster at 120–168 MHz using the LOw-Frequency ARray (LOFAR). We imaged a 132 deg2 region of the cluster, reaching an order-of-magnitude greater sensitivity than existing wide-field radio surveys of this field at three times higher spatial resolution compared to other low-frequency observations. We developed a tailored data processing strategy to subtract the bright central radio galaxy M 87 from the data. This allowed us to correct for systematic effects due to ionospheric variation as a function of time and direction. Results. In the final mosaic, which has a resolution of 9″ × 5″, we reach a median noise level of 140 μJy beam−1 inside the virial radius and 280 μJy beam−1 for the full area. We detect 112 Virgo member galaxies and 114 background galaxies. In at least 18 cases, the radio morphology of the cluster member galaxies shows clear signs of ram-pressure stripping. This includes three previously unreported candidates. In addition, we reveal previously undiscovered tails of 150 kpc in length from a previous epoch of AGN activity for NGC 4472 (M 49). While no cluster-scale diffuse radio sources are discovered, we find the presence of an extended radio signature of the W′ group. This feature is coincident with an X-ray filament detected with SRG/eROSITA in the outskirts of the cluster. We speculate that this emission is synchrotron radiation, which could be related to shocks or turbulence from accretion processes. Conclusions. The data published in this paper serve as a valuable resource for future studies. In the follow-up work of the ViCTORIA project, we will use these data for an analysis of environmental effects on the radio properties of star-forming galaxies in Virgo.

Are Brightest Cluster Galaxies Special?

A longstanding question when studying Brightest Cluster Galaxies (BCGs) and their evolution is whether their luminosities are statistical extremes of the luminosity distribution of other cluster member galaxies, or whether their luminosities follow a different,“special” distribution. To explore this question, we use the two Tremaine–Richstone (TR) statistics that sometimes lead to different conclusions about BCG specialness. We run Monte Carlo simulations of galaxy clusters with varying galaxy luminosity functions and BCG properties, and explore the sensitivity of the TR statistics to the adjusted BCG properties. We identify cases where the two statistics lead to different conclusions. The first of the TR statistics, which relies on the spread of BCG luminosities, is more reliable for determining specialness than the second, which utilizes the spread in the luminosity gap between the BCG and second ranked galaxy.

A Large-scale Kinematic Study of Molecular Gas in High-z Cluster Galaxies: Evidence for High Levels of Kinematic Asymmetry

We investigate the resolved kinematics of the molecular gas, as traced by the Atacama Large Millimeter/submillimeter Array in CO (2−1), of 25 cluster member galaxies across three different clusters at a redshift of z ∼ 1.6. This is the first large-scale analysis of the molecular gas kinematics of cluster galaxies at this redshift. By separately estimating the rotation curve of the approaching and receding sides of each galaxy via kinematic modeling, we quantify the difference in total circular velocity to characterize the overall kinematic asymmetry of each galaxy. 3/14 of the galaxies in our sample that we are able to model have similar degrees of asymmetry as that observed in galaxies in the field at similar redshift based on observations of mainly ionized gas. However, this leaves 11/14 galaxies in our sample with significantly higher asymmetry, and some of these galaxies have degrees of asymmetry of up to ∼50 times higher than field galaxies observed at similar redshift. Some of these extreme cases also have one-sided tail-like morphology seen in the molecular gas, supporting a scenario of tidal and/or ram pressure interaction. Such stark differences in the kinematic asymmetry in clusters versus the field suggest the evolutionary influence of dense environments, established as being a major driver of galaxy evolution at low redshift, is also active in the high-redshift universe.

Testing general relativity: New measurements of gravitational redshift in galaxy clusters

Context. The peculiar velocity distribution of cluster member galaxies provides a powerful tool to directly investigate the gravitational potentials within galaxy clusters and to test the gravity theory on megaparsec scales. Aims. We exploit spectroscopic galaxy and galaxy cluster samples extracted from the latest releases of the Sloan Digital Sky Survey (SDSS) to derive new constraints on the gravity theory. Methods. We considered a spectroscopic sample of 3058 galaxy clusters, with a maximum redshift of 0.5 and masses between 1014 − 1015 M⊙. We analysed the velocity distribution of the cluster member galaxies to make new measurements of the gravitational redshift effect inside galaxy clusters. We accurately estimated the cluster centres, computing them as the average of angular positions and redshifts of the closest galaxies to the brightest cluster galaxies. We find that this centre definition provides a better estimation of the centre of the cluster gravitational potential wells, relative to simply assuming the brightest cluster galaxies as the cluster centres, as done in past literature works. We compared our measurements with the theoretical predictions of three different gravity theories: general relativity (GR), the f(R) model, and the Dvali–Gabadadze–Porrati (DGP) model. A new statistical procedure was used to fit the measured gravitational redshift signal, and thus to discriminate among the considered gravity theories. Finally, we investigated the systematic uncertainties that possibly affect the analysis. Results. We clearly detect the gravitational redshift effect in the exploited cluster member catalogue. We recover an integrated gravitational redshift signal of −11.4 ± 3.3 km s−1, which is in agreement, within the errors, with past literature works. Conclusions. Overall, our results are consistent with both GR and DGP predictions, while they are in marginal disagreement with the predictions of the considered f(R) strong field model.

High-precision strong lensing models of galaxy clusters in the JWST era

AbstractWe present high-precision strong lensing models for the galaxy clusters MACS J0416.1–0403 at z=0.396 and Abell 2744 at z=0.307. The models are constrained by two of the largest data-sets of secure multiple images ever used in lensing. These are identified from the photometric images observed by the Hubble space telescope and JWST in combination with spectroscopic data obtained by the Multi-Unit Spectroscopic Explorer at the Very Large Telescope. The same spectro-photometric data are used to create pure and complete samples of cluster member galaxies. Our models allow an extremely precise estimation of the cluster total mass distribution and produce accurate magnification maps that are fundamental to study the physical properties (mass, size, luminosity, etc.) of the lensed high-redshift galaxies.

First JWST observations of a gravitational lens

We present our lens mass model of SMACS J0723.3−7327, the first strong gravitational lens observed by the James Webb Space Telescope (JWST). We use data from the Hubble Space Telescope and the Multi Unit Spectroscopic Explorer (MUSE) to build our ‘pre-JWST’ lens model and then refine it with newly available JWST near-infrared imaging in our JWST model. To reproduce the positions of all multiple lensed images with good accuracy, the adopted mass parameterisation consists of one cluster-scale component, accounting mainly for the dark matter distribution, the galaxy cluster members, and an external shear component. The pre-JWST model has, as constraints, 19 multiple images from six background sources, of which four have secure spectroscopic redshift measurements from this work. The JWST model has more than twice the number of constraints: 30 additional multiple images from another 11 lensed sources. Both models can reproduce the multiple image positions very well, with a δrms of 0.​​″39 and 0.​​″51 for the pre-JWST and JWST models, respectively. The total mass estimates within a radius of 128 kpc (roughly the Einstein radius) are 7.9−0.2+0.3 × 1013 M⊙ and 8.7−0.2+0.2 × 1013 M⊙ for the pre-JWST and JWST models, respectively. We predict with our mass models the redshifts of the newly detected JWST sources, which is crucial information, especially for systems without spectroscopic measurements, for further studies and follow-up observations. Interestingly, one family detected with JWST is found to be at a very high redshift, z &gt; 7.5 (68% confidence level), and with one image that has a lensing magnification of |μ| = 9.5−0.8+0.9, making it an interesting case for future studies. The lens models, including magnification maps and redshifts estimated from the model, are made publicly available, along with the full spectroscopic redshift catalogue from MUSE.

Galaxies in the central regions of simulated galaxy clusters

Context. Recent observations found that observed cluster member galaxies are more compact than their counterparts in ΛCDM hydrodynamic simulations, as indicated by the difference in their strong gravitational lensing properties, and they reported that measured and simulated galaxy–galaxy strong lensing events on small scales are discrepant by one order of magnitude. Among the possible explanations for this discrepancy, some studies suggest that simulations with better resolution and implementing different schemes for galaxy formation could produce simulations that are in better agreement with the observations. Aims. In this work, we aim to assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. Methods. We compared several zoom-in re-simulations of a sub-sample of cluster-sized haloes obtained by varying mass resolution and softening the length and AGN energy feedback scheme. We studied the impact of these different setups on the sub-halo (SH) abundances, their radial distribution, their density and mass profiles, and the relation between the maximum circular velocity, which is a proxy for SH compactness Results. Regardless of the adopted numerical resolution and feedback model, SHs with masses of MSH ≲ 1011 h−1 M⊙, the most relevant mass range for galaxy–galaxy strong lensing, have maximum circular velocities ∼30% smaller than those measured from strong lensing observations. We also find that simulations with less effective AGN energy feedback produce massive SHs (MSH ≳ 1011 h−1 M⊙) with higher maximum circular velocity and that their Vmax − MSH relation approaches the observed one. However, the stellar-mass number count of these objects exceeds the one found in observations, and we find that the compactness of these simulated SHs is the result of an extremely over-efficient star formation in their cores, also leading to larger than observed SH stellar mass. Conclusions. Regardless of the resolution and galaxy formation model adopted, simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged previous works. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the ΛCDM paradigm.

Dark matter in galaxy clusters: Parametric strong-lensing approach

We present a parametric strong-lensing analysis of three massive galaxy clusters for which Hubble Space Telescope imaging is available, as well as spectroscopy of multiply imaged systems and galaxy cluster members. Our aim is to probe the inner shape of dark matter haloes, in particular the existence of a core. We adopted the following working hypothesis: any group- or cluster-scale dark matter clump introduced in the modelling should be associated with a luminous counterpart. We also adopted some additional well-motivated priors in the analysis, even when this degraded the quality of the fit, quantified using the root mean square between the observed and model-generated images. In particular, in order to alleviate the degeneracy between the smooth underlying component and the galaxy-scale perturbers, we used the results from previous spectroscopic campaigns, which allowed us to fix the mass of the galaxy-scale component. In the unimodal galaxy cluster AS 1063, a core mass model is favoured over a non-core mass model, and this is also the case in the multimodal cluster MACS J0416. In the unimodal cluster MACS J1206, we fail to reproduce the strong-lensing constraints using a parametric approach within the adopted working hypothesis. We then successfully added a mild perturbation in the form of a superposition of B-spline potentials, which allowed us to obtain a decent fit (root mean square = 0.5″), and finally find that a core mass model is favoured. Overall, our analysis suggest evidence for core cluster-scale dark matter haloes in these three clusters. These findings may be useful for the interpretation within alternative dark matter scenario, such as self-interacting dark matter. We propose a working hypothesis for parametric strong-lensing modelling in which the quest for the best-fit model is balanced by the quest for presenting a physically motivated mass model, in particular by imposing priors.