Open Access
Issue
A&A
Volume 677, September 2023
Article Number A113
Number of page(s) 11
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202346594
Published online 12 September 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model.

Open access funding provided by Max Planck Society.

1. Introduction

Under the current standard cosmological model, galaxies form inside dark-matter halos due to gravitational instability, and grow in mass and size in a bottom-up fashion through mergers (e.g., White & Frenk 1991; Kauffmann et al. 1993). Dominating the massive end of galaxy populations, early-type galaxies (ETGs) are believed to be the end products of this hierarchical merging paradigm. The total mass (i.e., dark matter + baryonic matter) is arguably the single most important property of an ETG, and its distribution across the galaxy contains imprints from various key physical processes that have played a role in the formation and evolution of the galaxy. Dark-matter-only simulations suggest a “universal” density profile with a r−1 inner profile and a r−3 drop-off at large radii, independent of the halo mass (e.g., Navarro et al. 1997). Baryonic physics, on the other hand, is found to have substantially reshaped the mass structure of galaxies. In particular, cooling allows gas to condense and form stars in the central regions of galaxies, which leads to the contraction and dark matter (e.g., Velliscig et al. 2014). Heating due to supernovae or active galactic nucleus feedback, in contrast, can soften the central density concentration (e.g., Duffy et al. 2010). In addition, different types of merging events (e.g., dry versus wet and major versus minor) affect the total mass distribution of ETGs in different fashions (e.g., Ma & Boylan-Kolchin 2004; Naab et al. 2009; Remus et al. 2013). The state-of-the-art cosmological simulations predict that the total-density profile slopes (parameterized by the power-law density model of ρ ∝ rγ) of ETGs will, on average, become shallower from z ∼ 2 to z ∼ 0.8 and remain almost constant afterwards (e.g., Remus et al. 2017; Wang et al. 2019).

Due to the invisible nature of dark matter, observationally measuring the total mass distribution of ETGs turns out to be elusive. So far there have only been two major ways of achieving this, that is through either dynamics or gravitational lensing. In particular, spatially resolved stellar kinematics have been used for measuring the total mass distribution of low- and intermediate-redshift ETGs with z ≤ 0.33 (e.g., Lauer et al. 1995; Faber et al. 1997; Barnabè et al. 2011; Cappellari et al. 2015), but extending it to higher redshifts remains technically challenging due to cosmological dimming and apparent size decrease.

Strong gravitational lensing is the other widely used technique for probing the mass structure of ETGs. As a pure gravitational effect, the results of strong lensing are remarkably accurate and measurements of the total mass are robust against different galaxy models and assumptions on galaxy properties. In particular, the total mass enclosed within the Einstein radius can be determined with a few percent uncertainty (e.g., Koopmans et al. 2006; Auger et al. 2010a; Shu et al. 2016; Sonnenfeld et al. 2013a). More importantly, the constraining power of strong lensing does not degrade with galaxy redshift, which makes it a particularly valuable mass probe at high redshifts.

In addition, strong lensing is a powerful tool for addressing important cosmological and astrophysical questions, such as dark matter, dark energy, and the expansion of the Universe (e.g., Treu & Marshall 2016; Sharma & Linder 2022). For instance, strong gravitational lensing can detect dark matter subhalos (Vegetti et al. 2010; Hezaveh et al. 2016). Also combining total mass measurements from strong lensing with stellar mass estimations from stellar population synthesis or kinematics indicates that the stellar initial mass function (IMF) is not universal but rather varies with galaxy properties (e.g., Auger et al. 2010b; Spiniello et al. 2015). Making use of the time delays in strong lens systems, independent measurements of the Hubble constant have been obtained (e.g., Suyu et al. 2010, 2013; Millon et al. 2020; Wong et al. 2020; Birrer et al. 2020; Van de Vyvere et al. 2022; Shajib et al. 2023), which has been subject to extensive discussions lately (e.g., Abdalla et al. 2022; Vagnozzi et al. 2022; Riess et al. 2022).

In a strong lens system, the distortion of the background source is tightly related to the lensing potential. Modeling the imaging data of a strong lens system can therefore reconstruct the total mass distribution of the lensing galaxy (e.g., Tessore et al. 2016; O’Riordan et al. 2020; Gu et al. 2022). However, due to known degeneracies, most photometric observables remain unchanged under specific transformations of the total mass distribution, which makes the solution nonunique (e.g., Falco et al. 1985; Schneider & Sluse 2013, 2014). To break these degeneracies, stellar kinematics of the lensing galaxies have been introduced to further constrain the total mass distribution (e.g., Grogin & Narayan 1996; Romanowsky & Kochanek 1999; Blandford et al. 2001; Koopmans & Treu 2002, 2003; Warren & Dye 2003; Bolton et al. 2006; van de Ven et al. 2010; Barnabè et al. 2012, 2011; Yıldırım et al. 2020; Shajib et al. 2023). Such a joint strong lensing and dynamical analysis framework has been successfully applied to over 100 strong lens systems, which are mostly conducted on the Strong Lensing Legacy Survey (SL2S), the Sloan Lens ACS Survey (SLACS), the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) – a subsurvey of BOSS –, and BELLS for the GALaxy-Lyα EmitteR sYstems (BELLS GALLERY; e.g., Auger et al. 2010a; Brownstein et al. 2012; Bolton et al. 2012; Sonnenfeld et al. 2013a; Shu et al. 2016; Li et al. 2018) up to z ≈ 1. Most of these analyses to date are based on velocity dispersion measurements in a single aperture or a few bins. Spatially resolved kinematics for the lenses are not yet prevalent, but new facilities such as ground-based adaptive optics assisted integral field units (IFUs) and the James Webb Space Telescope will enlarge the sample of lenses with spatially resolved kinematic data.

Comparing the observed evolution of the total mass distribution of ETGs with cosmological simulations can provide a compelling test on the complex baryonic physics, dark matter-baryon interplay, and galaxy evolution models. However, obtaining accurate evolutionary trends from observations is not straightforward, as it requires identifying an ensemble of observed galaxies that are on the same evolutionary track. To achieve it, the commonly adopted approach so far has been to treat galaxies that satisfy certain observational criteria, which usually do not account for redshift evolution, as from the same population (e.g., Bolton et al. 2012; Sonnenfeld et al. 2013a). Such treatments are generally simplistic and can lead to biased conclusions.

In this paper, we propose a new strategy of tracking the evolution of galaxies by making use of cosmological simulations for selecting galaxies on the same evolutionary track, as an attempt of reaching less biased conclusions. Such a simulation-guided galaxy evolution inference also makes the comparison to simulations more straightforward. This paper is structured the following way: We provide a brief summary of the IllustrisTNG project in Sect. 2, which provides the state-of-the-art cosmological magnetohydrodynamic simulations that we choose to use in this work. Details about our proposed simulation-guided galaxy evolution inference are described in Sect. 3. As an illustration, the evolution of the total mass distribution of a strong lensing ETG at z = 0.884 evaluated following the proposed strategy is presented in Sect. 4. Discussions and conclusions are given in Sects. 5 and 6. Throughout this paper we adopt a flat ΛCDM cosmology with H0 = 67.74 km s−1 Mpc−1, ΩM = 0.3089, split in baryonic and dark matter (with Ωb = 0.0486), and ΩΛ = 0.6911, since these values are used in the IllustrisTNG simulations (Marinacci et al. 2018).

2. IllustrisTNG simulation

Big volume simulations aim to reproduce the structure of our universe, based on theoretical assumptions. Comparing the outcomes of simulations to observations allows a direct comparison of theoretical predicted structures and observed structures.

In the last years simulations moved to bigger volumes and at the same time better spatial and mass resolutions. The IllustrisTNG simulations are currently one of the biggest simulations at highest volumes available. The IllustrisTNG simulations come in three sizes, the TNG50, TNG100, and TNG300. In each case the name indicates the simulated box side-length in comoving Mpc. All simulations were carried out as dark-matter-only, as well as baryonic-matter simulations (Nelson et al. 2019).

The results of the IllustrisTNG simulations are found to be in good agreement within limits with observations. The simulated structure and galaxy formation, and the overall galaxy morphologies of both elliptical and spiral galaxies agree to a good degree with observations. The galaxy shapes and colors, the stellar populations, the star formation history, the galaxy stellar mass function, the stellar mass to halo mass ratio, and the luminosity of star forming galaxies are found to be close to observations for low and intermediate mass galaxies (e.g., Springel et al. 2018; Pillepich et al. 2018a; Rodriguez-Gomez et al. 2019; Lustig et al. 2023). Further, the simulated metallicities and radio emissions are in agreement with observed metallicities on both galaxy and cluster scale within their uncertainties for low and intermediate mass galaxies (e.g., Marinacci et al. 2018; Cannarozzo et al. 2023). Also on cluster scale the IllustrisTNG simulations shows agreement with observed morthologies (e.g., Pillepich et al. 2018b). Overall the simulations agree well with observations, but particularly for the most massive galaxies, a deviation from observations can be observed. For example, the simulations lack in number of massive galaxies, and their profile shows discrepancies with observations, see the sources above. Further investigating the discrepancies helps to understand how the simulations can be improved to match better with observations.

For galaxy evolution studies of massive ETGs, we want to study these most massive galaxies and use as big simulations as possible to get more massive simulated galaxies, but at the same time we want a good mass resolution. We thus chose the TNG300 baryonic simulation with a baryonic mass resolution of 1.1 × 107M and a dark matter mass resolution of 5.9 × 107M. The softening radius of the TNG300 simulation is 715 pc in comoving coordinates. These mass resolution and softening radius are a good compromise of sufficient resolution and large enough volume for our purposes of studying massive ETGs. The IllustrisTNG simulation data and descriptions are public available on the IllustrisTNG website1.

3. Methodology

3.1. Overview

The ideal way of probing the evolution of an object, is to measure its changes over time. This is, however, infeasible for studying galaxies, because most of the evolutionary processes involved happen on timescales much longer than the human lifetime. As a result, the most common way of quantifying galaxy evolution is by comparing properties of a sample of galaxies at different redshifts. These samples are supposed to be from the same population according to some usually basic and crude criteria. While this “ensemble approximation” is a feasible alternative, how to properly construct the ensemble is nontrivial, especially as the sample size is growing substantially and the call for less bias becomes increasingly more critical.

We hence propose a new strategy of linking observed galaxies at different redshifts by utilizing evolutionary tracks drawn from state-of-the-art cosmological simulations. More specifically, for any given observed galaxy at redshift z0, we first identify a collection of simulated galaxies at the same redshift that can be considered as its analogs. The simulated descendants (at z < z0) of the analog galaxies are then selected according to the individual merger trees. We use the properties of the simulated descendants to match them to observed galaxies at the same redshift. The evolution of the target galaxy can be studied by analyzing the so-selected observed galaxies at z < z0. Furthermore, by comparing the evolution trends of observed galaxies with simulations, all relevant ingredients adopted in the simulations, such as the cosmological model, dark matter properties, and sub-grid physics, can be tested.

3.2. Matching between observed and simulated galaxies

The objects we are studying are massive ETGs, which are found to follow a scaling relation among three basic properties: size, surface brightness, and velocity dispersion, that is the fundamental plane (FP; e.g., Magorrian et al. 1998; Koopmans et al. 2006). Studies have shown that the intrinsic scatter of the FP for massive ETGs is remarkably small (≲0.1 (0.07–0.09) dex) up to z = 0.5 (La Barbera et al. 2000; Bernardi et al. 2001; Jørgensen & Chiboucas 2013; Quenneville et al. 2022), and also in good agreement with the FP spanned by lensing galaxies (Auger et al. 2010a), which suggests that massive ETGs of similar size and surface brightness are likely self-similar in many other aspects. We therefore chose to match between observed and simulated galaxies using size and stellar mass, which is a proxy for the surface brightness and more easily available from the simulations.

To construct the simulated size-stellar mass relation, we selected subhalos, identified by the subfind algorithm (Springel et al. 2001), with the following properties from the TNG300 simulation:

  • total bound stellar mass between 1010.9 and 1012.5M;

  • star formation rate below 5 M yr−1;

  • less than 40% of stellar particles with circularity above 0.7.

The cut on the circularity (the last item) uses a built-in parameter of the simulation to remove clear disk dominated galaxies (described in the data specifications in table C2; Genel et al. 2015). By default, only the three-dimensional sizes calculated by the subfind algorithm are provided for the subhalos. In order to make a better comparison to observations we projected for every subhalo all its stellar particles along the z-axis and compute the projected half-stellar-mass radius. For the stellar mass, we directly used the total stellar mass from the subfind algorithm.

To construct the observed size-stellar mass relation, we used galaxies observed by the Sloan Digital Sky Survey (available catalogs were for galaxies from SDSS-I (Abazajian et al. 2009), and SDSS-III (Ahn & Alexandroff 2014; Alam et al. 2015), but there is no combined catalog for both) for which spectroscopic redshifts are available. The SDSS-I and SDSS-III are complementary in providing a sample of ETGs up to redshift ≈0.9: The SDSS-III BOSS survey focused on higher redshift galaxies with z > 0.45, which results in a lack of highly massive galaxies at lower redshifts (Loureiro et al. 2019), but the SDSS-I galaxy sample does not go to high enough redshifts. To obtain a subsample that can be compared to the simulated galaxies, we needed to apply additional cuts on several properties, including stellar mass. For SDSS galaxies, four different estimates for the stellar mass are available, which are computed by four groups using different methods, that are Granada group (Ahn & Alexandroff 2014), MPA-JHU group3, Portsmouth group (Thomas et al. 2013), and Wisconsin group based on principal component analysis method (Chen et al. 2012) for two different stellar population synthesis models (Bruzual & Charlot 2003; Maraston & Strömbäck 2011). The SL2S and SLACS lens masses use an exponential decay in their star formation history to determine the stellar masses of the galaxies. The only stellar mass catalogs using exponential decay are the MPA-JHU and the Granada group catalogs. To compare the lens galaxies to general ETG galaxies, the star formation history of the catalog should match. Further the flexible stellar population synthesis (FSPS) library used for the Granada catalog is also used for SLACS lenses. The MPA-JHU catalog only contains relatively low redshift galaxies, because it only contains galaxies up to DR8. As a consequence of the mentioned points, we used in this work stellar mass from the Granada group.

Similar to the selection of simulated galaxies, we only selected SDSS galaxies that have

  • Granada stellar mass between 1010.9 and 1012.5M;

  • Granada star formation rate below 5 M yr−1;

  • ellipticity less than 0.3;

  • error on i-band size less than 10%.

For galaxy size, we used the effective radius (in the intermediate-axis convention) of the best-fit de-Vaucouleurs model in the i band that is available in the SDSS imaging catalog.

Figure 1 shows the size and stellar mass distributions for simulated and observed galaxies in eight redshift bins from 0.1 to 0.89. The highest redshift bin was chosen to match the redshift of the target galaxy we are analyzing in the next section. The target observed galaxy has a redshift of z = 0.884, and the closest redshift in IllustrisTNG is at z = 0.89. For simulated galaxies, the redshift values are 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, and 0.89, while for observed galaxies, the redshift values are 0.1 ± 0.02, 0.2 ± 0.02, 0.3 ± 0.02, 0.4 ± 0.02, 0.5 ± 0.02, 0.6 ± 0.02, 0.7 ± 0.02, and 0.75 − 0.90. For the highest redshift bin of simulated galaxies at z = 0.89, we included observed galaxies in the redshift range from 0.75 to 0.90 in order to obtain a large enough sample. In each redshift bin, we fitted a linear relation between log10(Re/kpc) and log10(M*/M) for observed and simulated galaxies separately, and we report the best-fit parameters and the mean scatters in y-direction in log10(Re/kpc) with respect to the best-fit linear relations in Table 1.

thumbnail Fig. 1.

Size-stellar mass relations at eight redshift bins from 0.1 to 0.89. The y-axis is the half light radius in i-band in log10(Re/kpc) for observed galaxies or the projected half stellar mass radius in log10(Re/kpc) for simulated galaxies. The x-axis is the stellar mass in log10(M*/M) for observed galaxies or the total bound stellar mass of a subhalo in log10(M*/M). The observed data points, as well as their best fit, are shown in red, with the mean scatter of the points around the best fit as colored region around it. The median error of the observed data points is shown on the bottom right corner of the plot. The simulated data points and their best fit relation are shown in blue. The cyan colored dots are outliers of the simulated galaxies, identified by 3σ clipping, that were neglected in the fitting process.

thumbnail Fig. 1.

continued.

Table 1.

Best-fit size-mass relations of observed and simulated galaxies at different redshifts.

The slopes for the simulated relations are close to 0.7 in all the redshift bins considered, while there is a trend of flattening with cosmic time in the observed size-stellar mass relation. The typical uncertainties on the slope for the simulated data sizes is 0.007 and on the observed fit 0.02. We note that the shaded bands in Fig. 1 do not show the uncertainty of the best fit, but the mean scatter of the data points around the best fit. We found significant offsets between the observed and simulated size-stellar mass relations, especially in the six lower redshift bins (up to ≈0.25 dex in size). This difference has also been observed in previous studies, which suggests it to be primarily caused by offsets in size, because simulated sizes are measured from all star particles while observed sizes are usually measured from profile fitting that may miss the extended, low surface-brightness component, especially in lower redshifts (e.g., Genel et al. 2018; Rodriguez-Gomez et al. 2019). As a result, in each redshift bin, we chose to use stellar mass and the offset in size (in units of the measured scatter) with respect to the best-fit relation for galaxy matching.

4. Application to strong lens systems

As a proof of concept, we constructed an ensemble of strong lens systems at redshifts from ≈0.9 to ≈0.2 using the proposed strategy and characterized the evolution of their lensing galaxies in terms of the total mass distribution in the central region, particularly in terms of the radial profile of the total mass distribution. In Sect. 4.1 we describe the selection criteria of the lens systems. Section 4.2 gives an overview how we calculated the slope of the radial mass profile of the observed lens systems. In Sect. 4.3 we determine the slope of the radial mass profile of the simulated galaxies, and in Sect. 4.4 we show the results of the application of our comparison method on strong lens systems.

4.1. Lens system selection

The strong lens systems were selected from five samples, that are SLACS (Auger et al. 2010a), BELLS (Brownstein et al. 2012), BELLS GALLERY (Shu et al. 2016), SLACS for the masses (S4TM; Shu et al. 2017), and SL2S (Sonnenfeld et al. 2013b,a, 2015), that contain in total almost 200 strong lens systems with ETG lensing galaxies at redshifts ≈0.05–0.9. Since all the lensing galaxies in the first four samples were observed by the SDSS, we used the effective radius in the i-band from the SDSS catalog as their sizes, and their stellar masses were directly taken from the Granada catalog. For the SL2S lensing galaxies we used the half-light radii from best-fit de-Vaucouleurs models of the surface brightness of the HST images of the lens galaxies from Sonnenfeld et al. (2013b) as their sizes. The stellar masses of the SL2S lensing galaxies were estimated from fitting multiband CFHT photometry using a code implemented by Auger et al. (2009), who used this same code to estimate stellar mass for 49 SLACS lensing galaxies. As shown in Fig. 2, stellar masses estimated by Auger et al. (2009) are systematically lower than the Granada stellar masses by 0.19 dex on average and the scatter is 0.08 dex4. Because this bias does not vary noticeably with stellar mass, we added 0.19 dex to the reported stellar masses and added in quadrature the 0.08 dex scatter to the reported uncertainties for the SL2S lensing galaxies before galaxy matching.

thumbnail Fig. 2.

Comparison of SLACS X stellar mass with Granada stellar mass with Chabrier IMF. The y-axis shows the difference between the SLACS X stellar masses and the Granada masses. The mean difference is found to be −0.19, with a scatter of 0.08. The range of the scatter is indicated by the dotted lines. The errorbar on the left shows the typical uncertainties on the data points. The difference is possibly due to the different assumptions on the stellar initial mass function and the star formation history parameters employed in the models.

We selected the lensing galaxy in SL2SJ021801-080247 (Sonnenfeld et al. 2015) as our target galaxy based on its highest lens redshift among these five samples, which is at z = 0.884. The target galaxy has an estimated stellar mass of log10(M*/M) = 11.46 ± 0.17 (after conversion) and a size of log10(Re/kpc) = 0.91 ± 0.13. To identify analogs of the target galaxy in the simulations, we applied the following criteria:

  1. The simulated galaxies have the same stellar mass as the observed target galaxy within ±1σ uncertainties, that is logarithmic stellar masses in the range of 11.46 ± 0.17;

  2. The simulated galaxies and the observed target galaxy have the same relative offsets in size (after taking into account the scatters) with respect to the best-fit simulated and observed size-stellar mass relation respectively. More specifically,

    1. At the stellar mass of the target galaxy (log10(M*/M) = 11.46), the best-fit observed size-stellar mass relation in the z = 0.89 redshift bin suggests a mean logarithmic size of 0.76 with a scatter of 0.14 dex;

    2. The observed 1σ logarithmic size interval for the target galaxy, that is 0.78–1.04, thus corresponds to 1.13–1.98 times the scatter above the best-fit observed size-stellar mass relation;

    3. The best-fit simulated size-stellar mass relation in the same redshift bin suggests a mean logarithmic size of 0.7306 with a scatter of 0.15 dex at the stellar mass of 1011.46M (see Table 1);

    4. We require simulated galaxies to have logarithmic sizes in the range of 0.7495 (=0.7306 + 1.126 × 0.15) to 1.0281 (=0.7306 + 1.9833 × 0.15).

We considered the 158 simulated galaxies that satisfy these criteria as the simulated analogs for the target galaxy, and tracked every individual of them down to z = 0 to construct the simulated evolutionary trends of the target galaxy. The region of the 158 simulated analogs at z = 0.89 is indicated by the blue shaded box in Fig. 3.

thumbnail Fig. 3.

Matching simulated and observed galaxies. Left: size-stellar mass distribution of the simulated galaxies at z = 0.89 (blue) and the target galaxy at redshift z = 0.884 (red error bar). The 158 simulated galaxies in the blue shaded box are considered as the counterparts of the target galaxy. Right: blue points represent the size-stellar mass distribution of the 158 simulated analogs at z = 0.7. Their size and stellar mass have evolved, as indicated by the blue shaded box that brackets their distributions at z = 0.89. The red error bar indicates an observed lensing galaxy that was identified as a descendant of the target galaxy at z = 0.7. All the 158 simulated analogs and the identified observed descendant(s) are used to estimate the simulated and observed property (the total mass slope in this work) of the target galaxy at z = 0.7, respectively.

The same 158 simulated analogs were used to identify observed descendants of the target galaxy at progressively lower redshift bins. We first followed forward the merger trees of the simulated analogs to a desired redshift bin and updated their sizes and stellar masses. We note that the stellar mass of some simulated galaxies decreased at lower redshifts, likely because they experienced substantial interactions with other objects and lost mass due to tidal forces. We took the simulated galaxies that decrease in mass also into account in our comparison, because the same effect of tidal stripping can happen to observed ETGs. Using the z = 0.7 bin as an example, the observed descendants were selected from the lensing galaxies in the five strong lens samples as follows:

  1. For every simulated analog, we calculate its relative size offset with respect to the best-fit simulated size-stellar mass relation at z = 0.7;

  2. For every lensing galaxies within Δz = ±0.03 from the desired redshift bin, we calculate its relative size offset interval with respect to the best-fit observed size-stellar mass relation at z = 0.7 using its observed ±1σ logarithmic size interval;

  3. A lensing galaxy is considered as an observed descendant if its (observed) relative size offset interval and ±1σ logarithmic stellar mass interval overlap with at least one simulated analog.

Although all lensing galaxies satisfy the above requirements can be considered as the observed descendants of the target galaxy, for the purpose of this test case, we further required the selected observed descendants to have similar half-light radius to Einstein radius ratios as SL2SJ021801-080247 (i.e., 1.02) with a tolerance of ±20%. This assures that we compare the mass slopes within similar radii of the lensing galaxies. A much larger Einstein radius than half-light radius would measure the profile slope of more dark component than an Einstein radius at the same size of the half-light radius. The tolerance of 20% was found to be a suitable compromise between a large variety in the half-light radius to Einstein radius ratio and a sufficient number of observed analogs. With many more lens systems discovered in the near future, a tighter bound will be possible. For the considered seven lower redshift bins (see Table 1, 0–4 lensing galaxies satisfy the above requirements, which are used to construct the observed evolutionary trends of the target galaxy (Table 2).

Table 2.

Lens systems selected as possible lower redshift counter parts of system SL2SJ021801-080247.

4.2. Mass slope calculation of observed galaxies

We assumed that the total mass distributions of the lensing galaxies follow a power law profile, that is ρ(r) = ρ0rγ, which is found to be a good approximation for massive ETGs that were studied in this work (e.g., Kormann et al. 1994; Suyu et al. 2006; Koopmans et al. 2006; Auger et al. 2009; Sonnenfeld et al. 2013a; Shu et al. 2017). We then inferred γ (and ρ0) by comparing the lensing mass constraint from imaging data and the dynamical mass constraint from stellar velocity dispersion measurement. Following Kormann et al. (1994), the projected total mass within the Einstein radius (i.e., the lensing mass) is given by

M lensing = π θ E 2 D L 2 crit $$ \begin{aligned} M_{\rm lensing} = \pi \theta _{\rm E}^2 D_{\rm L}^2 \textstyle \sum _{\rm crit} \end{aligned} $$(1)

where θE is the Einstein radius inferred from lens modeling, DL is the angular diameter distance from the observer to the lens and ∑crit is the critical density defined as

crit = c 2 4 π G D S D L D LS $$ \begin{aligned} \textstyle \sum _{\rm crit} = \frac{c^2}{4\pi G} \frac{D_{\rm S}}{D_{\rm L} D_{\rm LS}} \end{aligned} $$(2)

with c the speed of light, G the gravitational constant, DS the angular diameter distance from the observer to the source, and DLS the angular diameter distance between lens and source.

The observed velocity dispersion is the projected and luminosity-weighted average of the three-dimensional velocity dispersion profile, which is determined by the total mass distribution. Assuming spherical symmetry and following the spherical Jeans equation, one can get the expression for the radial velocity dispersion profile of luminous matter σ2(r)

σ 2 ( r ) = G r d r v ( r ) M ( r ) r 2 β 2 r 2 β v ( r ) $$ \begin{aligned} \sigma ^2(r) = \frac{G \int _r^\infty \mathrm{d}r^{\prime } v(r^{\prime }) M(r^{\prime }) r^{\prime 2\beta -2}}{r^{2\beta } v(r)} \end{aligned} $$(3)

with G the gravitational constant, r the radius, v(r) the luminosity density, M(r) the mass enclosed at radius r, and β the velocity anisotropy parameter (see Binney 1980). The equation is valid under the assumption that the ratio of stellar number density to stellar luminosity density is spatially constant. The assumption usually holds within the half-light radius of ETGs. The luminosity density v(r) can be inferred from the observed surface brightness distribution via the Abel integral, and the anisotropy parameter is assumed to be zero.

The luminosity weighting of the velocity dispersion over the aperture can be expressed as

σ , | | 2 = A I ( R ) σ 2 ( R ) PSF R d R d θ A I ( R ) PSF R d R d θ $$ \begin{aligned} \langle \sigma ^2 _{*, ||} \rangle = \frac{ \int _\mathcal{A} I(R) \sigma _\parallel ^2(R) * \mathrm{PSF} \, R \mathrm{d}R \mathrm{d}\theta }{ \int _\mathcal{A} I(R) * \mathrm{PSF} \, R \mathrm{d}R \mathrm{d}\theta } \end{aligned} $$(4)

where I(R) is the surface brightness distribution (in our case a de-Vaucouleur profile), σ(R) is the projected velocity dispersion profile, PSF is the point spread function, and denotes convolution (see Suyu et al. 2010). The velocity dispersions of SL2S lenses were measured from spectra extracted within rectangular slits, so the above integration was done in Cartesian coordinates over the slit width and length. For the remaining lenses, whose velocity dispersions were measured from fiber-fed spectra, a circular aperture was used. The slit width and length for the SL2S lenses, the fiber radius Rfib for the other lenses, and the atmospheric seeing (full width half maximum) FWHM are reported in Table 2. We note that the predicted σ*, || is now a function of γ alone.

The best-fit γ value can be inferred by minimizing the χ2 defined as

χ 2 = ( σ , | | ( γ ) σ ap ) 2 σ unc 2 $$ \begin{aligned} \chi ^2 = \frac{(\sigma _{*,||} (\gamma ) - \sigma _{\rm ap} )^2}{\sigma _{\rm unc}^2} \end{aligned} $$(5)

with σ*, ||(γ) as the predicted velocity dispersion for different mass slopes γ, and σap as the measured velocity dispersion with uncertainty σunc. Since there is only one degree of freedom in the χ2 function, the 1-σ uncertainty on γ can be estimated by looking at γ values that correspond to Δχ2 = 1 relative to the minimum χ2 (e.g., see Chapter 7 of the book “Statistical data analysis” by Cowan 1997). The derived best-fit γ values and their uncertainties are reported in Table 2. As a consistency check, we show a comparison between our estimated γ values and the values estimated by Sonnenfeld et al. (2013a) for the five SL2S lenses involved with reported mass slope in Fig. 4, which suggests an almost perfect agreement. The system SL2SJ021801-080247 has no mass slope value provided in literature.

thumbnail Fig. 4.

Comparison of the γ values of this work with γ values from Sonnenfeld et al. (2013a). Our mass slopes agree well with the literature values, so we conclude that our code for determining the mass slope values works well.

4.3. Mass slope calculation of simulated galaxies

To determine the mass slopes of simulated galaxies, we fitted a power-law profile to their projected mass density distributions. The inner bound of the fit was set to the softening radius of dark matter particles of the Illustris TNG300-1 simulation, which is 1/h kpc = 1.48 kpc with h = 0.6774 in co-moving coordinates (see Springel et al. 2018). The outer bound of the fit was chosen to be the half stellar-mass radius multiplied by the Einstein-to-effective radius ratio of SL2SJ021801-080247 (i.e., 1.00/1.02). An example for the projected mass density and the power law fit for simulated galaxies is shown in Fig. 5. The plot shows the projected dark matter distribution (in blue), the projected stellar mass distribution (in yellow), the projected total mass distribution (in green), and the fit to the total mass distribution between the outer and inner bounds (in black, dotted). The plot is in logarithmic scales for both axes, so the power law fit corresponds to a linear fit. The relation between the mass slope and the total mass distribution is given by

Σ ¯ ( < R ) R 1 γ $$ \begin{aligned} \bar{\Sigma } ({ < } R) \propto R^{1-\gamma } \end{aligned} $$(6)

thumbnail Fig. 5.

Example for the projected mass density of a simulated galaxy at redshift 0.88 and the power-law fit between the outer and inner radial limit (the red lines). The slope of the fitted line is −1.146 and thus the power law mass profile slope γ = 2.146.

with γ the mass slope of the profile and R the radius.

4.4. Test case results

Figure 6 shows the measured total mass-density slopes of the 11 observed lensing galaxies and 158 simulated galaxies that are considered as the descendants of the lensing galaxy in SL2SJ021801-080247. By fitting a line to the measured γ values and uncertainties of the 11 observed galaxies, we found that the observed γ-z relation shows a slightly increasing trend from redshift 0.9 to redshift 0.2 with a slope of δγ/δz = −0.22 ± 0.18. For the simulated γ-z relation, we fitted a line to the mean γ values in individual redshift bins using the 1σ derivation of the mean as uncertainties. We found that the simulated γ-z relations shows a close-to-constant trend from redshift 0.9 to redshift 0 with a slope of δγ/δz = −0.01 ± 0.01. Nevertheless, given the small number of galaxies in each redshift bin, especially on the observation side, these results are by no means conclusive.

thumbnail Fig. 6.

Mass slope comparison of simulated (blue) and observed galaxies (red). The red line is a linear fit to the observed mass slopes over redshift. The blue shaded regions correspond to the 1-σ, 2-σ, and 3-σ bands inferred from simulated galaxies.

Interestingly, the above evolutionary trends appear to be qualitatively similar with previous findings based on larger samples of galaxies that were usually mass-selected. For example, Bolton et al. (2012) found δγ/δz = −0.60 ± 0.15 using ≈80 ETG strong lenses at z ≈ 0.1–0.7. Sonnenfeld et al. (2013a) found, by taking radius and mass dependence into consideration, the mean evolution for an individual galaxy dγ/dz to be on the level of −0.1, and Li et al. (2018) found with similar considerations ∂γ/∂z to be on the level of −0.3 using different samples of ETG strong lenses up to z ≈ 1.0. Remus et al. (2017) studied galaxies in the Magneticum Pathfinder simulations and found close-to-constant redshift evolution on the total mass-density slope since z ≈ 1 for ETGs. Wang et al. (2019) examined the IllustrisTNG simulations and found similar redshift evolution on the total mass-density slope since z ≈ 1 for massive ETGs (> 1011M).

We further compared the stellar mass evolution of the simulated and observed galaxies in Fig. 7. According to the 158 simulated analogs, the target galaxy is expected to become slightly more massive over time. From z = 0.88 to z = 0.1, the average stellar mass of the simulated analogs increases by 0.1 dex (from a logarithmic stellar mass of 11.45 to 11.55). The stellar mass evolution of the observed descendants agrees with the evolution proposed by the simulation. We chose to not compare the size evolution of observed and simulated galaxies, because the simulated and observed sizes are not directly comparable, as shown in Fig. 1.

thumbnail Fig. 7.

Total bound stellar mass of the simulated galaxies in blue, and the observed stellar mass in red plotted over the redshift. The plot suggests that there is no bias between the two samples.

5. Discussions

The evolution study presented above is clearly limited by the small sample size in both observation and simulation. From the ≈200 strong lens systems used in this work, we only found 0–4 galaxies in each redshift bin (and on average 1.6 galaxies per redshift bin) that are consistent with the descendants of the target galaxy, albeit the relatively loose selection criteria adopted. The situation is not expected to improve much, even if we extend to all currently confirmed strong lens systems with ETGs acting as lenses, which amount to about 4005. We suggest that an order of magnitude more strong lens systems will be needed in order to obtain a statistically significant result. Thanks to various large imaging surveys and advanced machine learning techniques, a few thousand promising strong lens candidates have recently been discovered (e.g., Chan et al. 2015; Huang et al. 2020, 2021; Lemon et al. 2020; Cañameras et al. 2021; Shu et al. 2022; Tran et al. 2022; Rojas et al. 2022). Upcoming surveys such as Euclid, Vera C. Rubin Observatory, and Chinese Space Station Telescope will further deliver ∼105 new galaxy-scale strong lens systems with high-resolution imaging data in the coming decade (e.g., Serjeant 2014; Collett 2015). It is equally crucial to carry out follow-up spectroscopic observations that are necessary for confirming the lensing nature and measuring stellar kinematics.

Ideally we want to not only compare lenses that are consistent with the descendants, but a set of lenses that have the same distribution in the observables. We have built that in by the requirement of having a similar half-light to Einstein radius ratio. As the number of observed strong lens systems continues growing in the near future, more stringent cuts on these two quantities can be used and additional requirements on properties, such as ellipticity and environment, can also be included to allow an even cleaner selection.

On the simulation side, we found 158 galaxies that can be considered as the analogs of the target galaxy in the TNG300 simulation, which is the largest we have access to. Although this appears much better than the situation on the observation side, it is mainly due to the large uncertainties on the stellar mass and size for the target galaxy. If the uncertainties shrink by a factor of 2 (to approximately 0.08 dex in stellar mass and 0.06 dex in size), we could only find 31 analogs from the simulation. This is because known strong lensing galaxies are predominantly massive ETGs (M* ≳ 1011M), which are simply rare. The Illustris team has recently completed the MillenniumTNG simulation6 (e.g., Barrera et al. 2022; Contreras et al. 2023; Hadzhiyska et al. 2023; Pakmor et al. 2023), which, compared to TNG300, has a similar mass resolution but a much larger box with cubic side length of approximately 500 Mpc/h. Given the factor of almost 15 increase in volume, much better statistics can be obtained when the MillenniumTNG simulation can be used.

Regarding the density slope inference, the uncertainty on γ is currently dominated by the velocity dispersion uncertainty, which is typically 10–20%. More precise inference can be obtained when higher-quality spectroscopic data are available. In addition, more accurate inference is expected when the single velocity dispersion measurement is superseded by two-dimensional stellar kinematics data (Barnabè et al. 2011).

Obviously, any observed trend constructed by our proposed strategy is simulation specific, since the observed and simulated galaxy selections depend on the actual merging history, as well as the implemented physics. In other words, the observed trend can be different when a different simulation is used for galaxy matching. As a result, comparing the observed trend constructed this way with the simulated trend can be used to test all relevant ingredients involved in the simulation, such as cosmological models, dark matter properties, and sub-grid physics.

6. Conclusion

In this work, we presented a new strategy of probing galaxy evolution by linking observed galaxies at different redshifts according to merger trees drawn from cosmological simulations. The so-linked observed galaxies were considered to be on the same evolutionary path and are therefore used to track how galaxies evolve. Such a simulation-guided galaxy evolution inference is supposed to be less biased than the commonly used approach, that relies on a sample of galaxies that are considered as a single population according to some usually basic and crude observational criteria.

As an illustration, we examined the evolution of the total mass-density slope of massive ETGs acting as strong lenses following this new strategy. We chose the strong lensing galaxy at z = 0.884 in SL2SJ021801−080247 as a starting point and identify 158 simulated galaxies at the same redshift from the TNG300 simulation as analogs of the target galaxy according to the stellar mass and size. We then followed the merger trees of the 158 simulated galaxies and successfully found, with the evolution of stellar mass and size taken into account, 11 observed strong lensing galaxies at redshifts of 0.7–0.2 that are analogs of the simulated galaxies and therefore considered descendants of the target galaxy.

The total mass-density slopes of the target galaxy and its 11 descendants were estimated from strong lensing and stellar kinematics data, and the slopes of the 158 simulated galaxies at eight redshift slices from 0.1 to 0.884 were computed directly from the total mass distributions. The observed results suggest the total mass-density slope of the target galaxy will likely increase with time, while the simulated results suggest a mild decreasing trend. Since the current observed results are heavily dominated by the Poisson fluctuation, we choose not to further investigate the apparent different results, which are nonetheless consistent with each other within 1.5σ.

We would like to emphasize that the proposed simulation-guided galaxy evolution inference is already feasible with currently available strong lens systems and cosmological simulations. In the coming decade, the number of strong lens systems will increase by two orders of magnitude. Cosmological simulations larger than the used TNG300 have already been developed and will keep upgrading. The proposed new strategy will take the advantage of this upcoming data explosion and play an even more important role in studying galaxy evolution.


4

Possibly due to the different assumptions on the stellar initial mass function and the star formation history parameters employed by Auger et al. (2009) and Granada.

5

According to the Master Lens Database:

https://test.masterlens.org/index.php

Acknowledgments

A.F. and S.H.S. thank the Max Planck Society for support through the Max Planck Research Group and the Max Planck Fellowship for S.H.S. A.F. was supported in part from NSERC discovery grant number RGPIN-2020-05102. Y.S. acknowledges the support from the Max Planck Society and the Alexander von Humboldt Foundation in the framework of the Max Planck-Humboldt Research Award endowed by the Federal Ministry of Education and Research and the China Manned Spaced Project (No. CMS-CSST-2021-A07 and CMS-CSST-2021-A12). X.H. acknowledges the University of San Francisco Faculty Development Fund. We thank Alessandro Sonnenfeld for helpful discussions, and the referee for helpful comments. A.F. thanks Laurence Perreault-Levasseur and Yashar Hezaveh for their support.

References

  1. Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543 [Google Scholar]
  2. Abdalla, E., Abellán, G. F., Aboubrahim, A., et al. 2022, J. High Energy Astrophys., 34, 49 [NASA ADS] [CrossRef] [Google Scholar]
  3. Ahn, C. P., & Alexandroff, R. C. A. P., et al. 2014, ApJS, 211, 17 [NASA ADS] [CrossRef] [Google Scholar]
  4. Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [Google Scholar]
  5. Auger, M. W., Treu, T., Bolton, A. S., et al. 2009, ApJ, 705, 1099 [Google Scholar]
  6. Auger, M. W., Treu, T., Bolton, A. S., et al. 2010a, ApJ, 724, 511 [NASA ADS] [CrossRef] [Google Scholar]
  7. Auger, M. W., Treu, T., Gavazzi, R., et al. 2010b, ApJ, 721, L163 [Google Scholar]
  8. Barnabè, M., Czoske, O., Koopmans, L. V. E., Treu, T., & Bolton, A. S. 2011, MNRAS, 415, 2215 [Google Scholar]
  9. Barnabè, M., Dutton, A. A., Marshall, P. J., et al. 2012, MNRAS, 423, 1073 [Google Scholar]
  10. Barrera, M., Springel, V., White, S., et al. 2022, MNRAS, submitted [arXiv:2210.10419] [Google Scholar]
  11. Bernardi, M., Sheth, R. K., Annis, J., et al. 2001, ArXiv e-prints [arXiv:astro-ph/0110344] [Google Scholar]
  12. Binney, J. 1980, MNRAS, 190, 873 [NASA ADS] [CrossRef] [Google Scholar]
  13. Birrer, S., Shajib, A. J., Galan, A., et al. 2020, A&A, 643, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Blandford, R., Surpi, G., & Kundić, T. 2001, ASP Conf. Ser., 237, 65 [NASA ADS] [Google Scholar]
  15. Bolton, A. S., Burles, S., Koopmans, L. V. E., Treu, T., & Moustakas, L. A. 2006, ApJ, 638, 703 [NASA ADS] [CrossRef] [Google Scholar]
  16. Bolton, A. S., Brownstein, J. R., Kochanek, C. S., et al. 2012, ApJ, 757, 82 [NASA ADS] [CrossRef] [Google Scholar]
  17. Brownstein, J. R., Bolton, A. S., Schlegel, D. J., et al. 2012, ApJ, 744, 41 [NASA ADS] [CrossRef] [Google Scholar]
  18. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cañameras, R., Schuldt, S., Shu, Y., et al. 2021, A&A, 653, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Cannarozzo, C., Leauthaud, A., Oyarzún, G. A., et al. 2023, MNRAS, 520, 5651 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cappellari, M., Romanowsky, A. J., Brodie, J. P., et al. 2015, ApJ, 804, L21 [NASA ADS] [CrossRef] [Google Scholar]
  22. Chan, J. H. H., Suyu, S. H., Chiueh, T., et al. 2015, ApJ, 807, 138 [NASA ADS] [CrossRef] [Google Scholar]
  23. Chen, Y.-M., Kauffmann, G., Tremonti, C. A., et al. 2012, MNRAS, 421, 314 [NASA ADS] [Google Scholar]
  24. Chen, Y., Li, R., Shu, Y., & Cao, X. 2019, MNRAS, 488, 3745 [NASA ADS] [CrossRef] [Google Scholar]
  25. Collett, T. E. 2015, ApJ, 811, 20 [NASA ADS] [CrossRef] [Google Scholar]
  26. Contreras, S., Angulo, R. E., Springel, V., et al. 2023, MNRAS, 524, 2489 [NASA ADS] [CrossRef] [Google Scholar]
  27. Cowan, G. 1997, Statistical Data Analysis (Oxford: Clarendon Press) [Google Scholar]
  28. Duffy, A. R., Schaye, J., Kay, S. T., et al. 2010, MNRAS, 405, 2161 [NASA ADS] [Google Scholar]
  29. Faber, S. M., Tremaine, S., Ajhar, E. A., et al. 1997, AJ, 114, 1771 [NASA ADS] [CrossRef] [Google Scholar]
  30. Falco, E. E., Gorenstein, M. V., & Shapiro, I. I. 1985, ApJ, 289, L1 [Google Scholar]
  31. Genel, S., Fall, S. M., Hernquist, L., et al. 2015, ApJ, 804, L40 [NASA ADS] [CrossRef] [Google Scholar]
  32. Genel, S., Nelson, D., Pillepich, A., et al. 2018, MNRAS, 474, 3976 [Google Scholar]
  33. Grogin, N. A., & Narayan, R. 1996, ApJ, 464, 92 [NASA ADS] [CrossRef] [Google Scholar]
  34. Gu, A., Huang, X., Sheu, W., et al. 2022, ApJ, 935, 49 [NASA ADS] [CrossRef] [Google Scholar]
  35. Hadzhiyska, B., Eisenstein, D., Hernquist, L., et al. 2023, MNRAS, 524, 2507 [NASA ADS] [CrossRef] [Google Scholar]
  36. Hezaveh, Y. D., Dalal, N., Marrone, D. P., et al. 2016, ApJ, 823, 37 [Google Scholar]
  37. Huang, X., Storfer, C., Ravi, V., et al. 2020, ApJ, 894, 78 [NASA ADS] [CrossRef] [Google Scholar]
  38. Huang, X., Storfer, C., Gu, A., et al. 2021, ApJ, 909, 27 [NASA ADS] [CrossRef] [Google Scholar]
  39. Jørgensen, I., & Chiboucas, K. 2013, AJ, 145, 77 [CrossRef] [Google Scholar]
  40. Kauffmann, G., White, S. D. M., & Guiderdoni, B. 1993, MNRAS, 264, 201 [Google Scholar]
  41. Koopmans, L. V. E., & Treu, T. 2002, ApJ, 568, L5 [NASA ADS] [CrossRef] [Google Scholar]
  42. Koopmans, L. V. E., & Treu, T. 2003, ApJ, 583, 606 [NASA ADS] [CrossRef] [Google Scholar]
  43. Koopmans, L. V. E., Treu, T., Bolton, A. S., Burles, S., & Moustakas, L. A. 2006, ApJ, 649, 599 [Google Scholar]
  44. Kormann, R., Schneider, P., & Bartelmann, M. 1994, A&A, 284, 285 [NASA ADS] [Google Scholar]
  45. La Barbera, F., Busarello, G., & Capaccioli, M. 2000, A&A, 362, 851 [NASA ADS] [Google Scholar]
  46. Lauer, T. R., Ajhar, E. A., Byun, Y. I., et al. 1995, AJ, 110, 2622 [NASA ADS] [CrossRef] [Google Scholar]
  47. Lemon, C., Auger, M. W., McMahon, R., et al. 2020, MNRAS, 494, 3491 [NASA ADS] [CrossRef] [Google Scholar]
  48. Li, R., Shu, Y., & Wang, J. 2018, MNRAS, 480, 431 [NASA ADS] [CrossRef] [Google Scholar]
  49. Loureiro, A., Moraes, B., Abdalla, F., et al. 2019, MNRAS, 485 [Google Scholar]
  50. Lustig, P., Strazzullo, V., Remus, R.-S., et al. 2023, MNRAS, 518, 5953 [Google Scholar]
  51. Ma, C.-P., & Boylan-Kolchin, M. 2004, Phys. Rev. Lett., 93, 021301 [NASA ADS] [CrossRef] [Google Scholar]
  52. Magorrian, J., Tremaine, S., Richstone, D., et al. 1998, AJ, 115, 2285 [Google Scholar]
  53. Maraston, C., & Strömbäck, G. 2011, MNRAS, 418, 2785 [Google Scholar]
  54. Marinacci, F., Vogelsberger, M., Pakmor, R., et al. 2018, MNRAS, 480, 5113 [NASA ADS] [Google Scholar]
  55. Millon, M., Galan, A., Courbin, F., et al. 2020, A&A, 639, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Naab, T., Johansson, P. H., & Ostriker, J. P. 2009, ApJ, 699, L178 [Google Scholar]
  57. Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [Google Scholar]
  58. Nelson, D., Springel, V., Pillepich, A., et al. 2019, Comput. Astrophys. Cosmol., 6, 2 [Google Scholar]
  59. O’Riordan, C. M., Warren, S. J., & Mortlock, D. J. 2020, MNRAS, 496, 3424 [CrossRef] [Google Scholar]
  60. Pakmor, R., Springel, V., Coles, J. P., et al. 2023, MNRAS, 524, 2539 [NASA ADS] [CrossRef] [Google Scholar]
  61. Pillepich, A., Nelson, D., Hernquist, L., et al. 2018a, MNRAS, 475, 648 [Google Scholar]
  62. Pillepich, A., Springel, V., Nelson, D., et al. 2018b, MNRAS, 473, 4077 [Google Scholar]
  63. Quenneville, M. E., Blakeslee, J. P., Ma, C.-P., et al. 2022, MNRAS, submitted [arXiv:2210.08043] [Google Scholar]
  64. Remus, R.-S., Burkert, A., Dolag, K., et al. 2013, ApJ, 766, 71 [NASA ADS] [CrossRef] [Google Scholar]
  65. Remus, R.-S., Dolag, K., Naab, T., et al. 2017, MNRAS, 464, 3742 [Google Scholar]
  66. Riess, A. G., Yuan, W., Macri, L. M., et al. 2022, ApJ, 934, L7 [NASA ADS] [CrossRef] [Google Scholar]
  67. Rodriguez-Gomez, V., Snyder, G. F., Lotz, J. M., et al. 2019, MNRAS, 483, 4140 [NASA ADS] [CrossRef] [Google Scholar]
  68. Rojas, K., Savary, E., Clément, B., et al. 2022, A&A, 668, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Romanowsky, A. J., & Kochanek, C. S. 1999, ApJ, 516, 18 [NASA ADS] [CrossRef] [Google Scholar]
  70. Schneider, P., & Sluse, D. 2013, A&A, 559, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Schneider, P., & Sluse, D. 2014, A&A, 564, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Serjeant, S. 2014, ApJ, 793, L10 [NASA ADS] [CrossRef] [Google Scholar]
  73. Shajib, A. J., Mozumdar, P., Chen, G. C. F., et al. 2023, A&A, 673, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Sharma, D., & Linder, E. V. 2022, JCAP, 2022, 033 [CrossRef] [Google Scholar]
  75. Shu, Y., Bolton, A. S., Mao, S., et al. 2016, ApJ, 833, 264 [Google Scholar]
  76. Shu, Y., Brownstein, J. R., Bolton, A. S., et al. 2017, ApJ, 851, 48 [Google Scholar]
  77. Shu, Y., Cañameras, R., Schuldt, S., et al. 2022, A&A, 662, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Sonnenfeld, A., Gavazzi, R., Suyu, S. H., Treu, T., & Marshall, P. J. 2013a, ApJ, 777, 97 [Google Scholar]
  79. Sonnenfeld, A., Treu, T., Gavazzi, R., et al. 2013b, ApJ, 777, 98 [Google Scholar]
  80. Sonnenfeld, A., Treu, T., Marshall, P. J., et al. 2015, ApJ, 800, 94 [Google Scholar]
  81. Spiniello, C., Barnabè, M., Koopmans, L. V. E., & Trager, S. C. 2015, MNRAS, 452, L21 [Google Scholar]
  82. Springel, V., Yoshida, N., & White, S. D. M. 2001, New A, 6, 79 [NASA ADS] [CrossRef] [Google Scholar]
  83. Springel, V., Pakmor, R., Pillepich, A., et al. 2018, MNRAS, 475, 676 [Google Scholar]
  84. Suyu, S. H., Marshall, P. J., Hobson, M. P., & Blandford, R. D. 2006, MNRAS, 371, 983 [Google Scholar]
  85. Suyu, S. H., Marshall, P. J., Auger, M. W., et al. 2010, ApJ, 711, 201 [Google Scholar]
  86. Suyu, S. H., Auger, M. W., Hilbert, S., et al. 2013, ApJ, 766, 70 [Google Scholar]
  87. Tessore, N., Bellagamba, F., & Metcalf, R. B. 2016, MNRAS, 463, 3115 [NASA ADS] [CrossRef] [Google Scholar]
  88. Thomas, D., Steele, O., Maraston, C., et al. 2013, MNRAS, 431, 1383 [NASA ADS] [CrossRef] [Google Scholar]
  89. Tran, K.-V. H., Harshan, A., Glazebrook, K., et al. 2022, AJ, 164, 148 [NASA ADS] [CrossRef] [Google Scholar]
  90. Treu, T., & Marshall, P. J. 2016, A&ARv, 24, 11 [Google Scholar]
  91. Vagnozzi, S., Pacucci, F., & Loeb, A. 2022, J. High Energy Astrophys., 36, 27 [NASA ADS] [CrossRef] [Google Scholar]
  92. van de Ven, G., Falcón-Barroso, J., McDermid, R. M., et al. 2010, ApJ, 719, 1481 [Google Scholar]
  93. Van de Vyvere, L., Gomer, M. R., Sluse, D., et al. 2022, A&A, 659, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Vegetti, S., Koopmans, L. V. E., Bolton, A., Treu, T., & Gavazzi, R. 2010, MNRAS, 408, 1969 [Google Scholar]
  95. Velliscig, M., van Daalen, M. P., Schaye, J., et al. 2014, MNRAS, 442, 2641 [Google Scholar]
  96. Wang, Y., Vogelsberger, M., Xu, D., et al. 2019, MNRAS, 490, 5722 [NASA ADS] [CrossRef] [Google Scholar]
  97. Warren, S. J., & Dye, S. 2003, ApJ, 590, 673 [Google Scholar]
  98. White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 [Google Scholar]
  99. Wong, K. C., Suyu, S. H., Chen, G. C. F., et al. 2020, MNRAS, 498, 1420 [Google Scholar]
  100. Yıldırım, A., Suyu, S. H., & Halkola, A. 2020, MNRAS, 493, 4783 [Google Scholar]

All Tables

Table 1.

Best-fit size-mass relations of observed and simulated galaxies at different redshifts.

Table 2.

Lens systems selected as possible lower redshift counter parts of system SL2SJ021801-080247.

All Figures

thumbnail Fig. 1.

Size-stellar mass relations at eight redshift bins from 0.1 to 0.89. The y-axis is the half light radius in i-band in log10(Re/kpc) for observed galaxies or the projected half stellar mass radius in log10(Re/kpc) for simulated galaxies. The x-axis is the stellar mass in log10(M*/M) for observed galaxies or the total bound stellar mass of a subhalo in log10(M*/M). The observed data points, as well as their best fit, are shown in red, with the mean scatter of the points around the best fit as colored region around it. The median error of the observed data points is shown on the bottom right corner of the plot. The simulated data points and their best fit relation are shown in blue. The cyan colored dots are outliers of the simulated galaxies, identified by 3σ clipping, that were neglected in the fitting process.

In the text
thumbnail Fig. 1.

continued.

In the text
thumbnail Fig. 2.

Comparison of SLACS X stellar mass with Granada stellar mass with Chabrier IMF. The y-axis shows the difference between the SLACS X stellar masses and the Granada masses. The mean difference is found to be −0.19, with a scatter of 0.08. The range of the scatter is indicated by the dotted lines. The errorbar on the left shows the typical uncertainties on the data points. The difference is possibly due to the different assumptions on the stellar initial mass function and the star formation history parameters employed in the models.

In the text
thumbnail Fig. 3.

Matching simulated and observed galaxies. Left: size-stellar mass distribution of the simulated galaxies at z = 0.89 (blue) and the target galaxy at redshift z = 0.884 (red error bar). The 158 simulated galaxies in the blue shaded box are considered as the counterparts of the target galaxy. Right: blue points represent the size-stellar mass distribution of the 158 simulated analogs at z = 0.7. Their size and stellar mass have evolved, as indicated by the blue shaded box that brackets their distributions at z = 0.89. The red error bar indicates an observed lensing galaxy that was identified as a descendant of the target galaxy at z = 0.7. All the 158 simulated analogs and the identified observed descendant(s) are used to estimate the simulated and observed property (the total mass slope in this work) of the target galaxy at z = 0.7, respectively.

In the text
thumbnail Fig. 4.

Comparison of the γ values of this work with γ values from Sonnenfeld et al. (2013a). Our mass slopes agree well with the literature values, so we conclude that our code for determining the mass slope values works well.

In the text
thumbnail Fig. 5.

Example for the projected mass density of a simulated galaxy at redshift 0.88 and the power-law fit between the outer and inner radial limit (the red lines). The slope of the fitted line is −1.146 and thus the power law mass profile slope γ = 2.146.

In the text
thumbnail Fig. 6.

Mass slope comparison of simulated (blue) and observed galaxies (red). The red line is a linear fit to the observed mass slopes over redshift. The blue shaded regions correspond to the 1-σ, 2-σ, and 3-σ bands inferred from simulated galaxies.

In the text
thumbnail Fig. 7.

Total bound stellar mass of the simulated galaxies in blue, and the observed stellar mass in red plotted over the redshift. The plot suggests that there is no bias between the two samples.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.