MUSE-ALMA Haloes X: The stellar masses of gas-rich absorbing galaxies

The physical processes by which gas is accreted onto galaxies, transformed into stars and then expelled from galaxies are of paramount importance to galaxy evolution studies. Observationally constraining each of these baryonic components in the same systems however, is challenging. Furthermore, simulations indicate that the stellar mass of galaxies is a key factor influencing CGM properties. Indeed, absorption lines detected against background quasars offer the most compelling way to study the cold gas in the circumgalactic medium (CGM). The MUSE-ALMA Haloes survey is composed of quasar fields covered with VLT/MUSE observations, comprising 32 \ion{H}{i} absorbers at 0.2 $<$ $z$ $<$ 1.4 and 79 associated galaxies, with available or upcoming molecular gas measurements from ALMA. We use a dedicated 40-orbit HST UVIS and IR WFC3 broad-band imaging campaign to characterise the stellar content of these galaxies. By fitting their spectral energy distribution, we establish they probe a wide range of stellar masses: 8.1 $<$ log($M_*$/M$_{\odot}$) $<$ 12.4. Given their star-formation rates, most of these objects lie on the main sequence of galaxies. We also confirm a previously reported anti-correlation between the stellar masses and CGM hydrogen column density N(\ion{H}{i}), indicating an evolutionary trend where higher mass galaxies are less likely to host large amounts of \ion{H}{i} gas in their immediate vicinity up to 120 kpc. Together with other studies from the MUSE-ALMA Haloes survey, these data provide stellar masses of absorber hosts, a key component of galaxy formation and evolution, and observational constraints on the relation between galaxies and their surrounding medium.


INTRODUCTION
One of the key questions in galaxy evolution is how galaxies interact and connect to their immediate environment, called the circumgalactic medium (CGM, Tumlinson et al. 2017).Once stars are formed, galaxies expel ionising photons and heavy elements formed in stars and supernovae into their surrounding environment through galactic winds (Pettini et al. 2008;Shull et al. 2014).Any gas flows into and out of the galaxy traverse through this extended gas halo and ★ E-mail: raugustin@aip.detherefore the CGM is closely connected to the galaxy's evolution (Muratov et al. 2015;Fox & Davé 2017).A detailed investigation of gas inflows and outflows is of paramount importance for understanding these processes.Since gas, stars, and metals are intimately connected, gas flows affect the history of star formation and chemical enrichment in galaxies.Therefore the study of the multi-phase CGM (extending over hundreds of kpc around galaxies - Shull et al. 2014) is crucial for understanding the star formation and galaxy evolution as a whole.
However, determining what drives the physical processes at play in the CGM still remains a complex problem in galaxy evolution, in large part due to the lack of significant observational constraints.Direct observations of this gas halo are challenging as the gas is extended and diffuse and the expected emission from the CGM is faint and mostly below current instrument detection limits (Augustin et al. 2019;Corlies et al. 2020).
In order to effectively explore the diffuse gas in the Universe, absorption-line spectroscopy proves to be a valuable method, as its detection sensitivity remains unaffected by redshift.One particularly compelling approach is to detect absorption lines in the spectra of bright background quasars, which allows for a comprehensive study of the distribution, chemical properties, and kinematics of CGM gas (e.g.Prochaska et al. 2005;Tumlinson et al. 2013).
Typically, strong (log [(H I)/cm −2 ]≥18) absorbers observed in the lines-of-sight to high redshift quasars serve as suitable targets to be followed up with additional imaging and spectroscopy to determine its relation to galaxies.A comprehensive understanding of the baryon cycle can be achieved by combining quantitative data on stellar properties with the already established gaseous measurements.A significant number of studies have successfully identified hosts to so-called DLAs (Damped Lyman- systems) and sub-DLAs (Péroux et al. 2003), which are the strongest H i absorbers (Augustin et al. 2018;Rhodin et al. 2018;Krogager et al. 2017).Due to their high column density these absorbers host a significant amount of neutral gas (Péroux & Howk 2020) and are therefore believed to be closely connected to their host galaxies.This makes them ideal targets to study the relation between galaxies and their CGM.
Recently, 3D Integral Field Spectroscopy (IFS), which produces data cubes where each pixel on the image has a spectrum, has provided a novel technique to examine the gas in absorption against background sources whose lines of sight pass through a galaxy's CGM.Building on early efforts with Near-Infrared (NIR) Integral Field Spectrograph (IFS) SINFONI (Spectrograph for INtegral Field Observations in the Near Infrared) on the Very Large Telescope (VLT) (Bouché et al. 2007;Péroux et al. 2011Péroux et al. , 2013Péroux et al. , 2016;;Augustin et al. 2018), the potential of this technique for studying the CGM with the optical IFU VLT/MUSE (Multi Unit Spectroscopic Explorer) has been demonstrated (Schroetter et al. 2016;Bouché et al. 2016;Zabl et al. 2019;Fumagalli et al. 2016;Muzahid et al. 2020;Dutta et al. 2021;Berg et al. 2022).Using this technique, a novel survey has been designed over the past years, combining IFU data from VLT/MUSE with ALMA and HST data of 19 quasar fields, labeled MUSE-ALMA haloes (Hamanowicz et al. 2020;Péroux et al. 2022;Weng et al. 2023a;Karki et al. 2023).Within these fields, 79 galaxies were identified to be associated within 500 km/s of strong quasar absorbers at  ∼ 1 (Péroux et al. 2022).By probing the neutral, molecular and ionised gas around these galaxies, the physical properties of the CGM have been characterised from the combination of abundance determination and kinematics through Voigt profile fitting of the absorber and emission-line and stellar continuum luminosities of the host galaxies (Augustin et al. 2018;Klitsch et al. 2018;Péroux et al. 2019;Hamanowicz et al. 2020).
The goal of the present work is to compute the stellar mass of galaxies related to gas-rich absorbers and understand what population of galaxies is traced by those absorbers.Particularly, we want to revisit the previously identified stellar mass -H i column density anticorrelation (Augustin et al. 2018) and fit its slope.Additionally, we determine whether the galaxies that are found as absorber hosts have physical properties which differ from the general population of galaxies, whether they are more star forming or quenched than the general galaxy population or whether they are following the same scaling relations as the general galaxy population.Particularly, Péroux et al. (2020) have shown that the metal distribution around galaxies is anisotropic.Their findings indicate that the metallicity trend with azimuthal angle is strongly dependent on the stellar mass of the host, indicating that the stellar mass is a critical component to the structure and evolution of galactic haloes.In this work we determine the stellar masses of galaxies associated with absorbers, providing essential information on these galaxies to better understand their relation to the gas seen in absorption.
The manuscript is organized as follows: Section 2 presents the observations used in this study.Section 3 details the process of fitting Spectral Energy Distribution to the data, while Section 4 relates the stellar properties of these galaxies with their CGM characteristics.We summarize and conclude in Section 5. Here, we adopt an  0 = 67.74km s −1 Mpc −1 , Ω  = 0.3089, and Ω Λ = 0.6911 cosmology.

Survey overview and relevant data
Integral Field Units (IFUs) have opened a new era in establishing the relation between absorption and emission.The optical IFU VLT/MUSE (Bacon et al. 2010) has proven to be a true game-changer in the field.The MUSE-ALMA Haloes survey probes the multi-phase CGM gas of intermediate redshift galaxies.The main goal of the survey is to reveal and understand the physical processes responsible for the transformation of baryons in galaxies (Péroux et al. 2022).The survey is based on a unique selection of known quasar absorbers with measured H i column density log [(H I)/cm −2 ]≥18, from HST UV spectroscopy with resolutions of =20,000-30,000.Weng et al. (2023a) have measured metal lines from optical spectroscopy and analysed the physical and emission-line properties of 79 galaxies associated with 32 H i absorbers at redshift 0.20 ≲  ≲ 1.4.These associated galaxies are selected to have velocities ≤ 500 kms −1 relative to the absorber redshift.They are found at impact parameters of ∼5.7-100 kpc at the lower redshifts and up to 270 kpc at the higher redshifts, where the upper cut-off reflects the field-of-view of MUSE.The star-formation rates (SFRs) of associated galaxies are measured using the H emission line when available.For sources at redshift  ≥ 0.4, where H is not observable with MUSE, Weng et al. (2023a) estimate the SFR using the [O ii] luminosity.For galaxies with dust corrections available, they also provide a dust-corrected SFR.3- SFR limits are calculated for non-detections.In the resulting sample, there are both passive galaxies without detectable emission lines, and star-forming galaxies with SFRs up to 15 M ⊙ yr −1 .

HST Broad-band Imaging Data
MUSE-ALMA Haloes also includes broad-band imaging of all but one field in the sample.This includes HST WFC3 imaging data in selected near-UV and optical filters, combined with archival Ultra-Violet, optical and Infra-Red WFPC2 and WFC3 imaging data to investigate the immediate surroundings of the quasar absorbers (Péroux et al. 2022).The filters for the observations are carefully designed to surround the 4000 Å break at each absorber's redshift.For the majority (16/19) of the fields images in three broad-band filters are available.Using additional data available in the HST archive, we cover a total of four filters for some of the targets.The MUSE IFU observations of the quasar fields are used to design the HST observations.We identified [OII] emission from galaxies at the redshift of the known absorbers as potential host galaxies and determined their location with respect to the quasar.The locations of these galaxies are used to optimize the HST pointings such that the diffraction spikes from the HST optics would not spatially overlap with the galaxies, while simultaneously keeping the whole MUSE FOV within the HST field of view.We then performed a careful quasar PSF subtraction, following Augustin et al. (2018), and photometry measurements (as described in detail in Péroux et al. 2022;Karki et al. 2023).The Astropy package Photutils is used on the processed HST broadband images in each filter for each field to search for all objects and to perform photometry of the detected objects.The details are described in Karki et al. (2023).In cases where an object detected in one filter was not detected in another filter, a 3- magnitude limit was calculated for the filter with the non-detection by measuring the 3- noise level within the same aperture as for the detection.The resulting magnitudes are publicly released with Péroux et al. (2022) and are used for the SED fitting described below.

SPECTRAL ENERGY DISTRIBUTION FITTING
Building on the successes of earlier multi-band approaches (e.g.Christensen et al. 2014;Augustin et al. 2018), we perform Spectral Energy Distribution (SED) fits of the galaxies associated with the absorbers.The wavelength coverage and SNR of the continuum emission of the galaxies in question in the VLT/MUSE observations alone does not provide sufficient data for state-of-the-art SED algorithms to converge for the redshifts we are probing.Instead, the new HST data cover wavelengths around the 4000 Å break and are designed to detect the stellar continuum of absorber counterparts in order to probe the stellar population of these galaxies.The coverage of this break significantly improves the accuracy of the stellar mass estimates, so we ensure that we obtain observations on both sides of the break.Additionally, the strength of the break directly affects the ability to measure the stellar masses effectively.We note that, highly star-forming galaxies have a less defined break than quenched galaxies.Therefore the masses of quenched galaxies are determined more robustly by the fitting procedure.
We use the Le Phare (Arnouts et al. 1999;Ilbert et al. 2006) SED fitting code to determine the physical properties of the observed galaxy sample.We assume the Bruzual & Charlot (2003) SED library, a Calzetti et al. (1994) dust extinction law and a Chabrier (2003) initial mass function.We focus here on the galaxies whose redshifts have been measured with VLT/MUSE spectroscopy and are known to be within 500 kms −1 from the strong H i absorbers.The resulting sample comprises 79 galaxies for 31 absorbers (see Table 1).To perform the SED fitting, we use the reliably determined magnitudes based on the high spatial resolution HST imaging.Some of the faint objects are detected in VLT/MUSE through their line-emission signatures, but not in continuum in the broad-band HST images, particularly in the bluest filters.In cases where galaxies had no detection in any HST filter, we make use of the photometry catalogs from the Dark Energy Legacy Survey (DECaLS, Dey et al. 2019) to perform the SED fitting.These systems are labeled as "LEGACY" in the last column of table 1.For galaxies with detections in one of the three bands, the SED fitting provides less robust stellar mass estimates, dubbed "probable".For illustration, we show in the Appendix examples of SED fits for a robust and a probable system (Fig. A1).These estimates are included in our sample for completeness, while we highlight the more reliable stellar mass measurements.We therefore focus here on the 34 galaxies with the most robust fits.
Figure 1 displays the distribution of stellar masses resulting from the SED fits.The grey histogram indicates the possible  * determinations where fewer than two photometry measurements are available while the red histogram shows the more robust determinations where 2-4 photometry measurements are available.The derived stellar masses cover a broad range from log [ * /M ⊙ ]= 7.8 -12.4,which illustrates the diversity of these systems.We see that the galaxies with stellar masses above log[ * /M ⊙ ]= 10 have their SED fit typically robustly determined while the lower mass galaxies often have non-detections in the bluer filters, limiting their stellar mass determination and potentially biasing our detection numbers and robust mass determinations towards higher masses.
We additionally run the Le Phare SED fitting code on all galaxies in the HST fields regardless of whether or not they have spectroscopic redshifts, in order to assess the potential of determining photometric redshifts with SED fitting.We derive uncertainties of the order  ≥0.5 in redshift, so we deem the results for such objects not robust enough to identify further absorber hosts beyond the ones detected through spectroscopy in MUSE cubes.

Dependence of SFR, redshift and impact parameter on M *
First we investigate how the determined stellar masses depend on other properties of the galaxy that are connected to its evolution.
In the upper left panel of Figure 2 we show the dependency of the determined stellar mass on the redshift of the galaxy.Given that galaxies grow and evolve over time we expect to find generally more massive galaxies at lower redshifts and less massive galaxies at higher redshifts.However, at the same time, due to surface brightness dimming, we expect to find preferentially bright and therefore massive galaxies at higher redshift and a range over all stellar masses at lower redshifts.What we find in our sample is a mix of those effects.Overall we find a large scatter of stellar masses at all redshifts, however we do notice a lack of lower mass systems at redshifts higher than  ≥ 0.7.
In the upper right panel of Figure 2, we investigate the dependency of the associated host stellar mass on the impact parameter to the absorber.We find no correlation between these two quantities, indicating that galaxies of any mass can be found at any distance to a known absorber.The diamonds indicate the sample of galaxies with lowest impact parameter for each absorption system, which are thought to be the main absorber hosts sample.Here we see that the In all these figures we denote robust stellar mass measurements in red and probable measurements in grey.The closest impact parameter galaxy of any system is marked with a diamond instead of a circle.Upper left: The redshift distribution of the determined stellar masses.We probe a range of stellar masses within 0.2≤z≤1.2,although there is a slight trend of finding the lower mass systems preferentially at lower redshifts.Upper right: The impact parameter range probed by the absorber hosts.We find no correlation between stellar masses and impact parameters.We consider the galaxy closest to the quasar sightline to be the main absorber host (diamond).The impact parameter range for these is up to ∼ 120 kpc.Bottom left: We show how the absorber associated galaxies in our sample fall into the stellar mass -SFR plane in comparison to the main sequence at a matching redshift range (Belfiore et al. 2018).We establish that most of these galaxies are normal star-forming galaxies and follow the main sequence.Some SFR measurements are found below the main sequence and could indicate quenching, particularly at the higher mass end.Bottom right: Same as the Figure on the left, but showing the specific SFR (SFR over stellar mass) instead of absolute SFR, hightlighing the deviation of the higher mass galaxies from the star forming main sequence.
maximum impact parameter we probe with this sample of main absorber hosts is ∼120 kpc.Following Behroozi et al. (2010) and Read et al. (2017), we approximate the halo mass of the galaxies in our sample as M halo ≈ 100 × M ★ .The dashed line in Figure 2 shows the resulting virial radius, assuming a standard Navarro-Frenk-White (NFW) profile (Navarro et al. 1996).We find that particularly at the lower mass end, a significant fraction of galaxies associated with the absorber lie outside the virial radius.As in any absorber-hostconnection study, there is a possibility that the closest galaxy contributing most the absorbing gas lies directly close to the projected position of the background quasar and/or is too faint to be detected.Weng et al. (2024) have recently quantifies this effect using hydrodynamical cosmological TNG50 simulations.Therefore, we cannot fully exclude the possibility of missing faint galaxies close to the quasar sightline.
In the bottom row (left and right panels) of Figure 2, we show the measurements for stellar mass from SED fitting and SFR (left) and specific SFR (right) from emission lines in MUSE (H and [OII], Weng et al. 2023a) compared to the expected main-sequence relation at =0.5 (Belfiore et al. 2018).We find that on average our sample of galaxies follows the main sequence, consistent with results of previous studies (e.g., Kulkarni et al. 2022).However, at higher masses the galaxies fall slightly below the main sequence.We recall that a fraction (17 systems) of the SFR measurements is dust corrected, specifically those with H and H measurements, corresponding to systems with  < 0.4 (see Weng et al. 2023a).Therefore estimates of SFR in the low-resdhift galaxies represent a lower limit on the total SFR.In particular, if the highest mass galaxies with log([ * /[M ⊙ ])> 11 were dust corrected, they might then well follow the main sequence.We furthermore note that the combined sample of H and [OII] determined SFRs may introduce a larger scatter in the SFR distributions than a more homogeneous sample would.We note that the main sequence should also flatten at the higher mass end of the mass distribution (Tomczak et al. 2016).Therefore we conclude that the sample of galaxies is representative of normal star-forming galaxies.Langan et al. (2023) and Guha & Srianand (2023) find a similar result for Mg ii absorption selected systems, indicating that they follow the main sequence.Together with our findings, this indicates that absorber host galaxies follow the trends of the typical galaxies population.2018) and Christensen et al. (2014).We find and quantify for the first time the anticorrelation between H i absorption strength and stellar mass of the main host galaxy (dotted line).The fit was performed on the stars and red diamond data points in the plot.Right: Summary of all measured quantities in this work.We show the impact parameter scaled with virial radius of each system versus the neutral gas H i column density.The color of each circle indicates the stellar mass of that object and the size of the circle is proportional to the star formation rate of this galaxy.We find that the galaxies associated with high column density H i absorbers span a range of SFR and stellar mass.As explored in the panel on the left, the highest column density systems seem to be hosted in groups of galaxies with smaller stellar masses, while the higher stellar mass objects are found at lower column densities.Simultaneously, at the higher column density end the host galaxies show relatively low SFRs compared to the galaxies found at the lower column density end of our sample.

Correlation between M * and CGM atomic gas content
One of the early findings of relating the H i absorption properties with the host galaxy properties was an apparent anti-correlation between stellar mass of the host galaxy and the absorption H i column density (Augustin et al. 2018).Such an anti-correlation seems counterintuititive and may be caused by additional factors such as measuring the impact parameter between the galaxy and the absorber, star formation rate of the host, or the environment of the galaxy.While previous studies were characterising somewhat limited sample sizes, we make use here of the currently complete MUSE-ALMA Haloes survey sample, covering quasar fields with 32 known absorbers and 79 associated galaxies with those absorbers, out of which we have robust masses for 34 galaxies.We are therefore testing this relation over a larger redshift and column density range and in more environments than previously possible.
In the left panel of Figure 3 we show for our sample the relation between the stellar mass of the host and the detected H i column density in its CGM.While previous studies that found this anti-correlation (Augustin et al. 2018;Rhodin et al. 2018) considered objects where each absorber was associated with only one -typically the closest -galaxy, here we are considering more complex systems, including groups of galaxies associated with a strong H i absorber.The reason for the single detections in the literature data, particularly the ones using SINFONI (Augustin et al. 2018) is the somewhat limited field-of-view around the quasar.To simplify and homogenize this investigation and make it comparable to previous studies we consider here the galaxy at the closest impact parameter, marked with a diamond, to be the main absorber host galaxy.The reason behind choosing the closest galaxy in impact parameter is the assumption that the physically closest galaxy to the gas we are probing in absorption will have the largest contribution to the absorption signal of all the galaxies within a given group (?).Now, we investigate this previously identified anti-correlation between stellar mass and circumgalactic H i column density for the larger sample.We also add the data points from Christensen et al. (2014) for comparison and note that all of these previous DLA host galaxy studies, have a larger average and greater spread in redshift, whereas our sample is at a slightly lower redshift on average.
Assuming an anti-correlation power law between stellar mass and H i column density, we fit the relation.We note that the upper limits were not considered in the fit.The data points used for the fit are marked with black dots in the left panel of Figure 3.For the fit we consider robust  * measurements of the closest impact parameter galaxies (red data points with diamond) as well as the literature samples from Christensen et al. (2014), Augustin et al. (2018) and Rhodin et al. (2018).These literature samples consist of single counterparts detections of DLAs and subDLAs and showed a tentative trend for this anti-correlation, which we now have enough data points to actually measure.Our choice of considering the closest galaxy in line of sight to the quasar homogenizes the chosen sample for the fit of the anticorrelation.We perform a Pearson correlation test and find a value of −0.57thus confirming the anticorrelation.Interestingly, we note that the trend is driven by the literature subsample.When focussing on the new measurements from this work, we find a Pearson coefficient of −0.36, indicating little anticorrelation.We fit a simple power-law to the relation in logarithmic space: The lack of high-mass high-N HI systems could be due to an observational bias or selection effect where high-N HI , high stellar mass systems are artificially excluded.Indeed, a high-N HI , high stellar mass system could arise in the center of a massive, potentially dusty galaxy, causing a reddening of the background quasar, which we would miss from typical quasar optical selection (Vladilo & Péroux 2005).
However, assuming this anti-correlation, or rather, upper envelope in detected N HI -stellar mass values has a physical origin, it holds information on the composition and evolution of galactic haloes.The lack of high-mass high-N HI systems can intuitively be explained by the ambient temperature of higher mass galaxies being typically higher than in lower mass galaxies, due to stronger feedback heating up the halo (see e.g.Suresh et al. 2017).This hotter gas may prevent the survival of cooler gas clouds hosting HI, at least within the inner ∼ 120 kpc as probed by our sample of closest impact parameter hosts.Likewise, the higher mass galaxies have undergone a rapid star forming process and already used up their cool H i gas reservoir, whereas lower mass galaxies take longer to use their gas supplies and still maintain more of it in their CGM.As we probe lower column density systems, the spread in stellar mass of associated galaxies increases, with the mean shifting towards higher masses.This means that we can find lower column density systems around a variety of galaxies but preferentially around higher mass systems as compared to the high N HI absorbers.In a lower column density regime, Berg et al. (2022) find a scatter of stellar masses log([ * /[M ⊙ ])∼ 8.5 -12 for absorber systems with  H i < 10 17 cm −2 without any clear trend.Similarly Chen et al. ( 2019) find column densities of ∼ 10 18 cm −2 around galaxies with log([ * /[M ⊙ ])∼ 7.5 -11.5.A reason for this phenomenon could be that more massive galaxies reside in larger haloes, providing a larger cross section for low column density absorption systems than lower mass galaxies.Our results are also in line with Khare et al. (2007) and Kulkarni et al. (2010) who argue that sub-DLAs hosts are more massive than DLA hosts.Finally, Dutta et al. (2023) recently found that Ly rest-frame equivalent widths show a peak at log([ * /[M ⊙ ])∼ 9 and decline at both high and low masses.They attributed this trend to high virial temperature and efficient feedback in high mass halos.Since most of the galaxies in the present study have log([ * /[M ⊙ ])> 9, and while we considered column densities rather than equivalent widths, our observed anticorrelation is aligned with their results.
Taken together, these findings suggest an evolutionary trend of the CGM composition with stellar mass where lower mass systems host haloes that are abundant with cool and dense H i gas and higher mass systems are more depleted of cool gas.However, even the more massive galaxies are fairly normally star forming, given the depletion of fuel in their immediate surroundings.More detailed studies of the gas content around galaxies and the galaxy properties such as stellar mass and star-formation rate but also metallicity, can measure the accretion timescales.These timescales describe the gas flows from the halo onto the disk and their transformation into stars versus the timescales for heating up the halo through feedback from a growing galaxy and depleting the gas fuel.

CONCLUSIONS
In this work, we have performed stellar mass measurements of 79 0.2<z<1.4galaxies known to have velocities within 500 kms −1 from a strong N(H i) absorber.To this end, we have used Le Phare SED fitting on multi-broad-band HST imaging of the fields to robustly determine the stellar masses of the subset of 34 galaxies that were detected in at least two imaging bands.
Our main results are: • We find that galaxies associated with high H i column density absorbers span a large range in stellar masses of log([ * /[M ⊙ ])= 8.1 -12.4 at 0.2<z<1.4.There is no siginificant trend with redshift or impact parameter.All these galaxies are normal star-forming galaxies following the expected main-sequence relation.
• A previously tentative detection of an anti-correlation between the stellar mass and CGM H i column density has been confirmed by the larger sample of absorber-host pairs, indicating an evolutionary trend of cool gas depletion in the CGM with stellar mass.
Ultimately, connecting the CGM gas properties of galaxies with their stellar content provides fresh clues on the baryon cycle, a key component of galaxy formation and evolution.In particular, knowledge of the stellar masses of galaxies allows to perform studies of the azimuthal metallicity distribution (Weng et al. 2023b).Further studies of the detailed chemical composition of absorber host galaxies may reveal the timescales involved in gas accretion and recycling.Cosmological simulations are now challenged to resolve both the cool gas in the CGM as well as the properties of the host galaxies in order to provide physical ground to our findings and make predictions to the galactic baryon cycle that can be tested by such observations.

Figure 1 .
Figure 1.Distribution of stellar masses of absorber host galaxies: We find a large spread of stellar masses for absorber-associated galaxies, with log[ * /M ⊙ ]= 7.8 -12.4.

Figure 2 .
Figure 2. Physical properties of the galaxies associated with gas-rich absorbers.In all these figures we denote robust stellar mass measurements in red and probable measurements in grey.The closest impact parameter galaxy of any system is marked with a diamond instead of a circle.Upper left: The redshift distribution of the determined stellar masses.We probe a range of stellar masses within 0.2≤z≤1.2,although there is a slight trend of finding the lower mass systems preferentially at lower redshifts.Upper right: The impact parameter range probed by the absorber hosts.We find no correlation between stellar masses and impact parameters.We consider the galaxy closest to the quasar sightline to be the main absorber host (diamond).The impact parameter range for these is up to ∼ 120 kpc.Bottom left: We show how the absorber associated galaxies in our sample fall into the stellar mass -SFR plane in comparison to the main sequence at a matching redshift range(Belfiore et al. 2018).We establish that most of these galaxies are normal star-forming galaxies and follow the main sequence.Some SFR measurements are found below the main sequence and could indicate quenching, particularly at the higher mass end.Bottom right: Same as the Figure on the left, but showing the specific SFR (SFR over stellar mass) instead of absolute SFR, hightlighing the deviation of the higher mass galaxies from the star forming main sequence.

Figure 3 .
Figure 3. Left: The relation between H i column density of the absorber versus the stellar masses of galaxies associated with those absorbers.The literature data are taken from Augustin et al. (2018); Rhodin et al. (2018) andChristensen et al. (2014).We find and quantify for the first time the anticorrelation between H i absorption strength and stellar mass of the main host galaxy (dotted line).The fit was performed on the stars and red diamond data points in the plot.Right: Summary of all measured quantities in this work.We show the impact parameter scaled with virial radius of each system versus the neutral gas H i column density.The color of each circle indicates the stellar mass of that object and the size of the circle is proportional to the star formation rate of this galaxy.We find that the galaxies associated with high column density H i absorbers span a range of SFR and stellar mass.As explored in the panel on the left, the highest column density systems seem to be hosted in groups of galaxies with smaller stellar masses, while the higher stellar mass objects are found at lower column densities.Simultaneously, at the higher column density end the host galaxies show relatively low SFRs compared to the galaxies found at the lower column density end of our sample.

Figure A1 .
Figure A1.The red points display the AB magnitudes measurements from the HST observations.The blue line represents the best fit galaxy spectrum derived from Le Phare algorithm.Left: An example of a fit based on detections of the galaxy in three HST bands and resulting in a fit providing a robust estimate of the stellar mass of the object.Right: An example of a fit based on a single magnitude measurement and providing a probable estimate of the stellar mass of that galaxy.