{'fraction_non_alphanumeric': 0.03807039845702012, 'fraction_numerical': 0.030443168368912463, 'mean_word_length': 4.509841806545103, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Multi-phase Circum-Galactic Medium probed with MUSE and ALMA

Galaxy halos appear to be missing a large fraction of their baryons, most probably hiding in the circumgalactic medium (CGM), a diffuse component within the dark matter halo that extends far from the inner regions of the galaxies. A powerful tool to study the CGM gas is offered by absorption lines in the spectra of background quasars. Here, we present optical (MUSE) and mm (ALMA) observations of the field of the quasar Q1130-1449 which includes a log [N(H I)/cm^-2]=21.71+/-0.07 absorber at z=0.313. Ground-based VLT/MUSE 3D spectroscopy shows 11 galaxies at the redshift of the absorber down to a limiting SFR>0.01 M_sun yr^-1 (covering emission lines of [OII], Hbeta, [OIII], [NII] and Halpha), 7 of which are new discoveries. In particular, we report a new emitter with a smaller impact parameter to the quasar line-of-sight (b=10.6 kpc) than the galaxies detected so far. Three of the objects are also detected in CO(1-0) in our ALMA observations indicating long depletion timescales for the molecular gas and kinematics consistent with the ionised gas. We infer from dedicated numerical cosmological RAMSES zoom-in simulations that the physical properties of these objects qualitatively resemble a small group environment, possibly part of a filamentary structure. Based on metallicity and velocity arguments, we conclude that the neutral gas traced in absorption is only partly related to these emitting galaxies while a larger fraction is likely the signature of gas with surface brightness almost four orders of magnitude fainter that current detection limits. Together, these findings challenge a picture where strong-HI quasar absorbers are associated with a single bright galaxy and favour a scenario where the HI gas probed in absorption is related to far more complex galaxy structures.


INTRODUCTION
Baryons are missing from galaxies in what is known as the galaxy halo missing baryon problem (McGaugh 2008). Indeed, galaxy halos appear to be missing approximately 60% of their baryons compared to expectations from the cosmological mass density, suggesting that these are structures nearly devoid of baryons both in mass and spatial extent. It has become apparent that galaxies also exhibit a diffuse baryonic component within the dark matter halo that extends far from the inner regions to the virial radius and beyond. This halo gas or circumgalactic medium (CGM) has been partly detected in the hot gas emitted around galaxies and detected in the X-rays regime (e.g. Anderson & Bregman 2010;Bregman et al. 2018;Nicastro et al. 2018) as well as in Ly-α (e.g. Wisotzki et al. 2016;Leclercq et al. 2017;Wisotzki et al. 2018). On the other hand, a powerful tool to study the cooler gas is offered by absorption lines in the spectra of background quasars as it allows us to study underdense gas regions with a sensitivity independent of redshift (e.g. Morris et al. 1993;Tripp, Lu & Savage 1998). Observations of redshifted UV absorption lines have indeed revealed much lower temperature and density gas in the CGM (Steidel et al. 2010;Rudie et al. 2012;Werk et al. 2013;Turner et al. 2014).
A remaining challenge is to relate the gas traced in absorption to the galaxy properties. A particularly powerful technique to bridge the gas component and stellar content of galaxies has come from integral field spectroscopy (IFS). Over the past few years, this technique has been used to study the CGM of high-redshift galaxies using near-infrared Integral Field Units (IFUs) by successfully detecting the galaxies responsible for strong H i absorbers at redshifts z ∼ 1 and z ∼ 2 in Hα and [N ii] emission (Bouché et al. 2007;Péroux et al. 2011a,b;Rudie, Newman & Murphy 2017). Together with long slit spectroscopy follow-up, these observations have enabled us to map the kinematics, star formation rate, and metallicity of this emission-line gas, and to estimate the dynamical masses of these galaxies Péroux et al. 2012;Bouché et al. 2013;Péroux, Kulkarni & York 2014). As optical IFUs became available, the study of the CGM using absorption techniques has expanded Bielby et al. 2016;Bouché et al. 2016;Fumagalli et al. 2016;Bielby et al. 2018) to also include lower-redshift objects (Péroux et al. 2017;Klitsch et al. 2018;Rahmani et al. 2018a;Péroux et al. 2018).
Here, we present results from new MUSE and ALMA observations of the field of background quasar Q1130−1449 with a strong H i absorber at z abs =0.313. The manuscript is organised as follows: Section 2 presents the observational set-up and data reduction. Section 3 details the analysis and results derived from these data, while Section 4 summarises our findings. Throughout this paper we adopt an H0 = 70 km s −1 Mpc −1 , ΩM = 0.3, and ΩΛ = 0.7 cosmology. At the redshift of the absorber, 1 arcsec corresponds to 4.6 kpc.
From a HST/Faint Object Spectrograph (FOS) UV spectrum, the z abs =0.313 absorber was determined to be a damped Lyα (DLA) system with a large H i column density of neutral gas: log [N (HI)/cm −2 ]=21.71±0.07 (Lane et al. 1998). Lane et al. (1998) also detected 21cm absorption spanning FWHM=42.1±2.7 km s −1 at the redshift of that absorber with the Westerbork Synthesis Radio Telescope (WSRT). This was later confirmed by Kanekar et al. (2009) who derived a spin temperature Tspin=820±145K and covering factor, f=0.9 from Very-long-baseline interferometry (VLBI) observations. Bechtold et al. (2001) additionally reported strong X-ray absorption in excess of the Galactic value in the Chandra quasar spectrum which they interpreted as associated with that absorber. Based on this assumption, they derived the abundance of the oxygen group elements, independently of dust depletion, of 23% solar. Kanekar et al. (2014) computed the zinc metallicity to be [Zn/H]= −0.80 ± 0.16 relative to solar (i.e. the neutral gas metallicity is 16% solar) from a medium-resolution HST/STIS spectrum. Guber, Richter & Wendt (2018) used the VLT/UVES spectrum of the quasar to fit the multicomponent profile in Ti, Mn and Ca which are likely depleted into dust grains. In addition, they estimated the H i column density of each component from the observed Ca ii column density pattern, assuming a constant Ca ii/H i ratio. Interestingly, the thus estimated H i column density for each component, on its own, exceeds the canonical DLA definition (log [N (HI)/cm −2 ] >20.3).
We obtained the VLT/UVES reduced quasar spectra from the ESO advanced data products archives 1 . Fig.A1 of the appendix A presents our independent fit to the metal absorption profiles, which is consistent with findings from Guber, Richter & Wendt (2018), adding lines from Fe ii, Mg i and the saturated Mg ii doublet. The zero velocity in the figure is set to the systemic redshift of galaxy G0, zG0 =0.31305 (see later section). Most of the absorption lies blueward of this systemic redshift.
2.2 Field Imaging: Identifying Absorbing Galaxies

New MUSE Observations
We report new MUSE medium-deep (5-hours on source) observations of this field. The observations were carried out in service mode (under programme ESO 96.A-0303, PI: C. Peroux) during the nights of 7-8 Dec 2015 (Texp=4×960s), 5-6 March 2016 (4×1200s) and 19-20 Feb 2017 (8×1200s). The seeing constraint for these observations was < 0.8 arcsec and natural seeing mode was used. The field was rotated by 90 degrees between exposures. These were further divided into two equal sub-exposures, with an additional field rotation of 90 degrees and sub-arcsec dithering offset in 2-step pattern to minimise residuals from the slice pattern. The field of view is 60 arcsec × 60 arcsec, corresponding to a 0.2 arcsec/pixel scale. We used the "nominal mode" resulting in a spectral coverage of ∼4800-9300Å. At the redshift of the target (z abs =0.313), the data cover emission lines from [O ii] λλ3727, 3729 to Hα. The spectral resolution is  ) showing objects with bright continua. The exposure time for the HST image is 4400 s, while for MUSE it is 14400 s. The quasar is at the center of the image. North is up, East to the left. The fields-of-view are ∼1×1'. MUSE observations additionally provide spectroscopic and kinematic information for the majority of the objects in the field. The quasar is shown as a white star and the positions of the galaxies associated with the quasar absorber are marked with light-blue squares. Four of the galaxies previously identified are detected (namely G1, G2, G4 and G6) and their numbering system follows Kacprzak, Murphy & Churchill (2010). The other 7 objects (G0 and >G16, underlined) are new identifications at z abs =0.313. Some objects are bright in the continuum but undetected in pseudo narrow-band images because of weak emission lines (i.e. G4, G20 and G21).
R=1770 at 4800Å and R=3590 at 9300Å resampling the whole spectrum to a spectral sampling of 1.25Å/pixel. The ESO MUSE pipeline (Weilbacher 2015) version v1.6 was used to reduced the data together with additional external routines for sky subtraction and extraction of the 1D spectra. We used master bias, flat field images and arc lamp exposures from observations taken closest in time to the science frames to correct each raw cube. We processed the raw science data with the scibasic and scipost recipes, correcting the wavelength calibration to a heliocentric reference. The wavelength solutions were checked using the known wavelengths of the night-sky OH lines and were found to be accurate within 15 km s −1 . We centered the individual exposures with the exp align recipe using the point sources in the field to ensure accurate relative astrometry. We combined the individual exposures into a single data cube using the exp combine recipe. We measured the seeing of the final combined data from the quasar and other bright point sources in the cube. We calculated that the point spread function (PSF) has a full width at half maximum of 0.76 arcsec. Finally, to estimate the flux uncertainty, we compared the flux calibration with the R magnitude of the quasar and a bright star in the field. We computed the MUSE broadband fluxes in the R-band filter. We determined the flux error to be ±25%.
The removal of OH emission lines from the night sky was accomplished with an additional purpose-developed code. The scipost recipe was performed with sky-removal method turned off. After selecting sky regions in the field, we created PCA components from the spectra which were further applied to the science datacube to remove sky line residuals (Husemann et al. 2016;Péroux et al. 2017). The resulting MUSE white-light image is shown on the right panel of Fig. 1. We retrieved HST/WFPC2 F814W image of the field from the HST archive, and we show it on the left panel of Fig. 1.

New ALMA Observations
The field of Q1130−1449 was observed with ALMA in Band 3 to cover the CO(1-0) emission lines at the redshift of the z abs = 0.313 absorber (under programme 2016.1.01250.S, PI: C. Peroux) on 4, 8, and 15 Dec 2016. The precipitable water vapour (PWV) varied between 1.5 and 5.4 mm and the total on-source observing time was Texp = 3.6 hrs. A compact antenna configuration (C40-3) resulted in an angular resolution 1.2×1.8". One of the four spectral windows was centred on the redshifted CO(1-0) line frequency of 87.8 GHz and used the high spectral resolution mode, providing 3840 channels of each 0.488 MHz wide. The other three spectral windows were set to low spectral resolution mode and were used for continuum observations of the field. The quasars J1058+0133 and J1139−1350 were used as amplitude and phase calibrators, respectively.
We started the data reduction with the pipelinecalibrated uv-data sets, as delivered by ALMA. Additional data reduction steps were carried out with the Common Astronomy Software Applications (CASA) software package version 4.7.0. Minor manual flags were added to remove some uv-data points with strongly outlying amplitudes. Next, we applied two cycles of self-calibration on the data. The quasar in the centre of our science field is very bright at mm frequencies (∼800 mJy), which makes the field very suitable for self-calibration of both the phases and amplitudes.
First, we used the tclean recipe to Fourier-transform the uv-data into a continuum image. The self-calibration procedure was performed using the tasks gaincal and applycal on individual measurement sets to create corrected uv-data, after which the data were re-imaged to create an improved continuum map. We applied one round of phase self-calibration and one round of amplitude and phase selfcalibration.The next step was to subtract the bright continuum source from the field using the task uvsub. We then built a cube with the tclean task, setting the pixel size to 0.2", so as to oversample the beam sufficiently, and using a 'robust' weighting scheme with a Briggs parameter of 0.5. One of the four spectral windows was centred on the redshifted CO(1-0) line frequency of 87.8 GHz and used the high spectral resolution mode, using a binning of 8, providing 3840 channels of each 3.906 MHz wide. We applied a spectral binning of seven channels to achieve a velocity resolution of 53 km s −1 . We removed remaining continuum signatures around the imperfectly subtracted quasar with the uvcontsub task. Lastly, we produced a cube corrected for the primary beam using the pbcor task. The final rms noise level of the cube is 0.11 mJy/beam per 53 km s −1 channel, and the angular resolution is 1.4" × 1.3", corresponding to 6 kpc at the redshift of the target. We note that the resulting FWHM of the primary beam of ALMA in band 3 is about 60", conveniently matching the MUSE field-of-view.
In addition, we created a small cube around the quasar with the highest spectral resolution (∼7 km s −1 ) to search for absorption lines, but do not detect the CO(1-0) absorption line at the absorber redshift in the spectrum of the background quasar. Given the rms, we compute an optical depth limit of τ =0.0125 km s −1 at 5-σ assuming a 10 km s −1 wide absorption line following Zwaan et al. (2015) calculation. Assuming the excitation temperature equates the TCMB at z=0.313, we derive a stringent limit on the CO column density of N(CO)<2×10 13 cm −2 (Mangum & Shirley 2015). Using the mean CO/H2 column density ratio of 3×10 −6 from Burgh, France & McCandliss (2007), this results in a limit of N(H2)<7×10 18 cm −2 .

ANALYSIS AND RESULTS
3.1 The Environment of z abs ∼0.3 H i-rich Absorbers 3.1.1 Galaxies at z abs =0.313 in the Field of Q1130−1449 The field is part of an observing campaign targeting quasar absorbers whose redshifts allow detection of Hα with MUSE, with measured N (H I) from HST UV spectroscopy and a known associated galaxy. On top of this, the field of Q1130−1449 has been the target of a number of observing programs in the past which we summarise here. We use the numbering system introduced by Kacprzak, Murphy & Churchill (2010) when relevant and add to it for the new detections. Fig. 1 illustrates this naming scheme. Bergeron & Boisse (1991) first reported emission (and absorption) lines at the redshift of the absorber in two objects (G2 and G4). Deharveng, Buat & Bergeron (1995) had successfully detected Hβ (but no Ly-α) in G2 from HST/FOS spectroscopy. Lane et al. (1998) later obtained a spectroscopic redshift of another galaxy (G1) which is also consistent with the absorption redshift. Kacprzak, Murphy & Churchill (2010) spectroscopically identified two additional galaxies (G6 and G14), the latter of which is outside our MUSE field. Here, we used the HST images as a prior to extract spectra from the MUSE cube. To this end, we used Sextractor to identify and deblend faint objects in the broad band WFPC2 image. Using the mpdaf 2 routine (Piqueras et al. 2017), we extract the pixels at these coordinates above a 0.75-σ threshold. We then reviewed the extracted spectra to identify those matching z abs =0.313. In some cases, no emission lines were detected, but strong Ca ii absorption lines and a Balmer break matching z abs =0.313 were identified (e.g. G20). In addition, we made pseudo narrow-band images of the MUSE cube at the expected position of the emission lines. An example of this is shown in Fig. 2 Fig. 4 (top panel) for Hα. This approach reveals additional objects which remain undetected in the HST image because of their faint continuum level (e.g. G0, G16 and G17).
In particular, we discovered a faint object, G0, closer to the quasar at impact parameter b = 2.3 arcsec (10.6 kpc). We note that it is not matched by any object detected in Gemini AO-images by Chun et al. (2006) because of the different wavelength coverage (optical versus H-band) and different sensitivity (emitters versus continuum objects). While the object lies under the quasar PSF, several emission lines are clearly identified in the MUSE spectrum. In the following analysis, we use the systemic redshift of G0 as reference for the analysis of the absorption profile plotted in velocity (see appendix A) simply because it has the smallest impact parameter to the quasar of all objects spectroscopically detected in the field. In essence, this system is like the "Galaxies with background QSOs" searches which have found quasars shining through low redshift, foreground galaxies at small impact parameters (<10 kpc) within the SDSS fibers Straka et al. 2013;Straka et al. 2015;Joshi et al. 2017). These, together with recent findings of low impact parameter objects in high spatial resolution HST imaging (i.e. Augustin et al. 2018), provide observational evidence to the postulate of York et al. (1986) that "some quasar absorption-line systems may arise when a quasar sight line intersects a Magellanic-type irregular galaxy (i.e., a gas-rich dwarf)". Fig. 1 (right panel) summarises the numbering system of the new identifications. The corresponding MUSE spectra are displayed in Fig. 5. In total, we report 11 objects at the redshift of the absorber in our MUSE data, 7 of which are new discoveries.
Similarly, we identified the galaxies at the absorber redshift in the ALMA observations by running the Duchamp 3 source finder (version v1.6.2) which performs a threedimensional search via wavelet reconstruction with a detection threshold of 5-σ on the cube before primary beam cor-  rection. We find three secure emitters associated with galaxies G2, G4 and G6, respectively. Fig. 3 shows the CO(1-0) contours overlaid on the background HST image. In addition, we performed a visual inspection of the cube at the redshifted frequency of the CO(1-0) line. We detect a fourth emission line, 26.4 arcsec south of the quasar. The MUSE spectrum shows a featureless continuum object consistent with either CO(2-1) at z=1.614 or CO(3-2) at z=2.921 (a higher CO transition would have matching Ly-α line in the MUSE spectrum). Interestingly, one of the continuum spectral window also reveals an emission line at 101.4 GHz which is identified as CO(2-1) at z=1.274 thanks to the detection of the [O ii] doublet in the corresponding MUSE spectrum. These findings illustrate the power of combining ALMA with MUSE observations to securely assess the redshift of single line mm detections. Fig. 2 (bottom panel) shows the 0 th moment map of the field. The spectra are extracted within the CASA software and shown in Fig. B2. We also built a continuum image of aggregate bandwidth to check for possible detection of the galaxies in continuum. None of the objects are detected in the map with rms=6 µJy, corresponding to a 5-σ continuum flux limit of 30 µJy.

Cosmological zoom-in Simulations
In order to reproduce the typical environment of z∼0.3 quasar absorbers, we turn to cosmological hydrosimulations. Our calculations are based on dedicated RAM-SES Adaptive Mesh Refinement (AMR) simulations with a total run time of 1.3 million CPU hrs (Frank et al. 2012). We zoomed on a region around the most massive halo with a box of size of ∼14 Mpc/h. The higher-level of refinement allows us to reach a spatial resolution of ∼380pc/h (comoving) at z=0. In short, the simulations include non-thermal supernova feedback with "on-the-fly" self-shielding option (Teyssier et al. 2013). With the goal of reproducing the expected surface brightness in the CGM regions of galaxies, we have post-processed the simulations taking into account gravitational cooling due to collisional ionisation of accreting gas, photo-ionisation by external UV sources and scattering from star forming regions. Details of the simulation are presented in an upcoming paper (Augustin et al. submitted).
Current observations find that the star formation rate (SFR) of strong H i-absorbers at z∼0.3 is just below 1 M /yr (Rahmani et al. 2016  . To best match these observable constraints, we picked a typical z=0.3 halo in the zoom-in simulations which best mimics the properties of H i absorbers: a SFR=1 M /yr), Mstar=10 11 M and MDM=10 12 M (Rvir=200 kpc/h comoving). Fig. 4 presents a MUSE-size FoV Hα emission surface brightness map of that simulated halo in units of erg s −1 cm −2 arcsec −2 . These simulations show the presence of multiple galaxies at the same redshift, typical of a small group. The group members have a broad range of masses so that some of these objects will remain undetected in current state-of-the-art observations. While the simulated objects appear more extended, the predicted surface brightnesses from the simulations are comparable to our Hα MUSE observations. The robustness of the surface brightnesses prediction is a remarkable achievement which indicates that collisional and photo-ionisation in combination can be the physical processes at the origin of the emission seen in the observations. In addition, this demonstrates that the brightest components of the CGM of low-redshift galaxies can be probed with current observations. Another important feature of the simulated cube is the presence of faint low surface-brightness gas with large sky cross-section in between galaxies.

The Typical Environment of low-redshift Absorbers
Over the years, there has been mounting evidence that bright galaxies can be found at the redshift of quasar absorbers (see Krogager et al. 2017, for a recent review). Tremendous progress has come from IFU observations in particular which have provided a robust means to remove the bright quasar contamination and reach low impact parameters Péroux et al. 2012;Schroetter et al. 2016). At low redshifts, however, it is becoming clear that the picture is more complex than a simple one-to-one correspondence between a single galaxy and a quasar absorber. Thanks to the MUSE wide FoV, it is now possible to survey a large portion of sky in one observational set-up, with physical sizes (300 kpc across) typical of the CGM of galaxies (Steidel et al. 2010). Clearly, the case presented here does not contain the typical isolated galaxy-absorber match often reported in past searches for absorber counterparts. In part, this is due to deeper data as we are reaching limits of SFR>0.01 M yr −1 , i.e. an order of magnitude deeper than similar z∼1 studies. In addition, we report the detection of passive galaxies whose redshift estimate comes solely from Ca H&K absorption features and Balmer breaks detected in their continuum. While galaxies at these low redshifts will be more clustered, finding simlar objects at higher redshifts is challenging as the continuum flux gets fainter.
While we note the possible alignments of the objects in what could be a filament with a north-west to south-east  In addition, the predicted surface brightness from the simulations is comparable to our Hα MUSE observations indicating that collisional and photo-ionisation in combination can be the physical processes at the origin of the emission seen in the observations. These results demonstrate that the brightest components of the CGM of low-redshift galaxies can be probed with current observations of the emission, while the diffuse gas in between objects is at present best traced in absorption. These findings add to the paradigm shift where our former view of strong-N (H I) quasar absorbers being associated with a single bright galaxy changes towards a picture where the H i gas probed in absorption is related to far more The panels for G1 and G2 show key emission and absorption lines (at redshifted wavelengths), which lines can also be noted in galaxies plotted below. In some cases, no emission lines were detected, but strong Ca ii absorption lines and a Balmer break matching z abs =0.313 were identified (e.g. G20). The spectrum of G21 is smoothed with a 3-pixel boxcar. In other cases, the galaxies have strong emission lines but faint continuum level (e.g. G0, G16 and G17) so that such objects would not appear in broad band imaging. In particular, we discovered a faint galaxy closer to the quasar G0 at impact parameter b = 2.3 arcsec (10.6 kpc). While the object lies under the quasar PSF, several emission lines are clearly identified in the MUSE spectrum. In total, we report 11 objects at the redshift of the absorber in our MUSE data, 7 of which are new discoveries. The systemic redshift of galaxy G0 is used as zero velocity reference. The galaxies are ordered from bottom to top as a function of impact parameter to the quasar line-of-sight (see right y-axis for values in kpc). The light blue horizontal lines indicate the Vmax extent for galaxies exhibiting ordered rotation. The overlap in velocity space with the closest galaxies is somewhat limited, and much of the absorption bluewards of G0 could be related to the low surface brightness tidal gas predicted in our simulations to be 4 orders of magnitude below the current detection threshold.
orientation (Møller & Fynbo 2001;Fumagalli et al. 2016;Péroux et al. 2017), it is also possible that these objects are part of a galaxy group. Assuming G2 as the central member galaxy, we note that G6, G4, G19 and G21 (lying to the north) are at velocity positions blueshifted with respect to G2 while G0, G16 and G17 (lying to the south) are redshifted with respect to G2. This alignment shows an overall direction of rotation in the halo of the group. The velocity dispersion of the system is ∼170 km s −1 which is in the range of small galaxy groups. Based on such dispersion and assuming a virialised spherical system, we obtain a virial radius of 420 kpc so that the MUSE FoV probes 1/3 rd of the virial radius of the group. The virial mass for such system will be 2. . The cosmological zoomin simulations presented in the section 3.1.2 are specifically evolved to redshift zero with the goal to predict the emissivity of CGM gas around low-redshift galaxies thus enabling a more direct comparison with the observations presented here. They indicate that strong-H i absorbers could be associated with faint low surface-brightness intra-group gas. In the simulations, this gas is typically remnant tidal debris from previous interactions between the main galaxy and other smaller satellite galaxies. Indeed, the selection of absorbing systems on the basis of their sky coverage will favour interacting systems that have spread their gas around. The MUSE white light image (Fig. 1) shows clear indications of extended gas around G2 and G4 in particular.
This lay-out bears resemblance to cold gas streams similar to the Magellanic Stream (Richter et al. 2014). Lowredshift analogues of such groups include the M81 (NGC 3031)/M82 (NGC 3034)/NGC 3077 complex where optical observations tracing starlight indicate well-separated objects. Twenty-one cm maps however show H i gas distributed in between these objects dominated by filamentary structures, clearly demonstrating the violent disruption of this system by tidal interaction (Yun, Ho & Lo 1993. Similarly, de Blok et al. (2014) report the detection of a cloud of neutral gas with admittedly lower column density outside the main H i disk of NGC2403 which they argue could be either accreting from the intergalactic medium, or the result of a minor interaction with a neighbouring dwarf galaxy. Interestingly, they note that the velocities of the H i in NGC 2366 in the same group completely overlap with those of the H i in NGC 2403, which would make their absorption signature indistinguishable in the spectrum of a background object. Twenty five years ago, Morris & van den Bergh (1994) had already calculated that the space density of Ly-α absorbers and local small groups matched, thus proposing that quasar absorbers can be produced by tidal debris in groups of galaxies.
These results pose a new challenge to the interpretation of the neutral gas probed in absorption. Indeed, most studies of the absorbing galaxies have concentrated on the brightest objects which are likely to be the main group member. Others have reasonably argued that the absorption is likely to arise from the galaxy with the smallest impact parameter to the quasar line-of-sight. Indeed, we note that lower-mass galaxies such as G0 revealed here by the new MUSE observations of Q1130−1449 are expected to be one of the main contributors to galactic winds, given that the gas can more easily escape the galaxies' smaller potential wells. Fig. 6 shows the galaxies' positions with respect to the absorption profile in velocity space. The galaxies are ordered from bottom to top as a function of impact parameter to the quasar lineof-sight (see right y-axis for values in kpc). The light blue horizontal lines indicate the Vmax extent for galaxies indicating rotation. The overlap in velocity space with the closest galaxies is somewhat limited, and much of the absorption bluewards of G0 could be related to the low surface brightness tidal gas predicted in our simulations to be four orders of magnitude below the current detection threshold. These results strenghten previous reports of low-redshift strong-H i absorbers associated with groups (Péroux et al. 2017;Bielby et al. 2017;Rahmani et al. 2018a;Klitsch et al. 2018;Rahmani et al. 2018b;Klitsch et al. 2019;Borthakur et al. 2018). Altogether, these findings add to the paradigm shift where our former view of strong-N (H I) quasar absorbers being associated to a single bright galaxy changes towards a picture where the H i gas probed in absorption is related to far more complex galaxy structures. Since both winds and interactions are predicted to be enhanced in the distant Universe, these factors will likely become more important at higher redshifts but harder to discern.

Physical Properties of the Galaxy Group
The MUSE observations provide a wealth of information on the stellar properties of the galaxies in the group. Accurate redshifts are determined from multiple lines detected in the 1D MUSE spectra calibrated in air. We measured the emission fluxes from a Gaussian fit to the nebular lines. The 3σ upper limits for non-detection are computed for an unresolved source spread over 6 spatial pixels and spectral FWHM = 2.4 pixels = 3Å. We followed the method described in Péroux, Kulkarni & York (2014) based on Hα and Hβ to estimate moderate E(B-V) and further correct the observed fluxes for dust attenuation. The resulting fluxes are tabulated in Table 1. The SFR with 1σ uncertainties were estimated from these Hα fluxes following the relation of Kennicutt (1998). The emission line galaxies have values of SFR ranging from 0.1 to a few solar masses per year (see Table 2).
Following the methodology of Péroux et al. (2017), the HII region metallicities were measured from the strong-line indices (N2, O3N2 and R23 most relevant branch) following prescription by Kobulnicky, Kennicutt & Pizagno (1999) and the mean value of each of these was computed. Our estimates agreed within the errors with earlier measurements of G2 and G6 metallicities by Kacprzak, Murphy & Churchill (2010). The metallicity of the galaxies in the group ranges from [X/H]=−0.57±0.17 to 0.09±0.16. The total metallicity of the neutral gas probed in absorption, [Zn/H] = −0.80±0.16, is therefore consistent within the error bars with the lowest of the metallicities found in the galaxies in the MUSE image.
While metallicity gradients in galaxies could partly reconcile these measurements (Carton et al. 2018), we note that the low dust-corrected abundances measured in absorption resemble the properties of the extragalactic diffuse gas produced in our numerical simulations. Furthermore, the nebular emission line ratios place all but one galaxy into the star-forming region of the BPT diagram (Baldwin, Phillips & Terlevich 1981 Table 2. We used the ALMA observations to constrain the molecular gas content of the galaxy group. We measured the flux density in the lines by integrating the CO(1-0) emission lines in the extracted spectra (Fig. B2). We then computed the corresponding luminosities using the following formula: where dL is the luminosity distance, Fint the flux density and f obs =87.8GHz is the observed frequency (Solomon, Downes & Radford 1992). Because we observe the CO(1-0) transition directly, our data naturally overcome un-certainties related to unknown CO line ratios. Estimates of the corresponding H2 molecular masses however rely on the rather uncertain CO-H2 conversion factor (Bolatto, Wolfire & Leroy 2013). Previous works have used Galactic CO-H2 conversion factor of αCO=4.6 M (K km/s pc 2 ) −1 . Recent results by Klitsch et al. (2019) find that this is not always appropriate for absorbing-galaxies so that here we opt to express the molecular masses as a function of αCO. The resulting molecular masses range from M mol = (1.00±0.05) to (3.16±0.07) × (αCO/1 K km/s pc 2 ) × 10 9 M and are listed Table 3.
The ALMA observations indicate substantial cold gas reservoirs in three of the galaxies in the group. These molecular masses are large in comparison with the typical stellar mass of absorbing galaxies . In fact, a number of recent works have reported CO associated with quasar absorbers from ALMA observations with a high rate of incidence and derived equally high molecular masses (Neeleman et al. 2016;Møller et al. 2018;Neeleman et al. 2018;Kanekar et al. 2018;Klitsch et al. 2018;D'Odorico et al. 2018;Klitsch et al. 2019). However, it is unclear why H i-rich absorbing galaxies should have large reservoir of H2 gas. At first sight, this might appear at odds with searches for H2 rotational and vibrational transitions at UV restframe wavelengths which have only reported both low incident rates (10%) and low molecular gas content (f mol 0.01) with respect to the neutral gas (Ledoux, Petitjean & Srianand 2003;Noterdaeme et al. 2008). While the most abundant molecule after H2, CO is known to be less abundant with a CO/H2 column density ratio ranging from 10 −7 to 10 −5 and a mean value of 3×10 −6 (Burgh, France & Mc-Candliss 2007). However, it must be borne in mind that the neutral atomic phase of the gas traces significantly lower gas densities than the molecular gas traced by CO (Snow & McCall 2006). Indeed the observations of nearby galaxies indicate that H i and H2 are not co-spatial (Schruba et al. 2011), with the molecular gas concentrated in the central part of galaxies so that the detection of H2 in absorption at an impact parameter distance at the edges of galaxies is less likely (Zwaan & Prochaska 2006;Obreschkow & Rawlings 2009).
In the field of Q1130−1449, we detected three objects in CO at the absorber redshift, two of which have large impact parameters (b>50 kpc). The objects close to the quasar (i.e. G0 and G1) do not show CO(1-0) emission, although G1 is the most actively star-forming member of the group. We calculated the molecular depletion times for the three objects detected with ALMA following: τ depl =MH2/SFR=1/SFE, where SFE is the star formation efficiency. To compute the depletion time scale, we assumed a Galactic CO-H2 conversion factor of αCO=4.6 M (K km/s pc 2 ) −1 . We derive τ depl ranging from 1.7±0.5 to 12±5 Gyr (see Table 3). The gas depletion times are long compared to the timescales of the processes driving the evolution of the interstellar medium state ). Semenov, Kravtsov & Gnedin (2017) argue that star-forming gas converts only a small fraction (∼1%) of its mass into stars while most of it (∼90%) is dispersed by dynamical and feedback processes. Selection based on strong H i-absorbers thus possibly preferentially finds galaxies with large molecular gas reservoirs (at given SFR) which inefficiently convert their gas into stars . Table 1. Dust-corrected nebular emission line fluxes. The fluxes are expressed in units of erg s −1 cm −2 . The quoted errors are 1σ uncertainties. Galaxies G1 to G6 (shown in italics) were previously identified, while all the remaining are new identifications. Note: a Galaxies showing NaID absorption features in their spectra. b The flux limits are not corrected for dust. Table 2. Ionised gas properties of the galaxies detected at the redshift of the z abs =0.313 absorber in the MUSE observations. We use the numbering system introduced by Kacprzak, Murphy & Churchill (2010) when relevant and extend it for the seven new detections. δ is the angular distance from the quasar in arcsec and b is the impact parameter in kpc. SFR estimates are dust-corrected. G20 and G21 have no, or weak, emission lines. Galaxies G1 to G6 (shown in italics) were previously identified.

Morphological and Kinematical Properties of Group Members
MUSE and ALMA observations of the field are also complementary in that they both provide information on the kinematics of the stars and molecular gas respectively. We performed a 3-dimensional morpho-kinematic analysis of the objects detected in MUSE and ALMA, using the method described in Péroux et al. (2017). In brief, we used the GalPak 3D algorithm  which compares directly the data cube with a parametric model mapped in x, y, λ coordinates. The algorithm uses a Markov Chain Monte Carlo (MCMC) approach with a proposed distribution of the parameters in order to efficiently probe the parameter space. The code has the advantages that it fits the galaxy in 3-d space and provides a robust description of the morpho-kinematics of the data. We assessed the quality of the fit from a χ 2 minimisation as well as the inspection of the residual maps. We mostly used the brightest [O iii] 5007 A line to perform the analysis. Our results are tabulated in Table 4 and illustrated in Fig. B1 and B2 in appendix B for the MUSE and ALMA detections respectively. The CO velocity fields are well resolved in our ALMA observations. As expected, the stellar sizes as measured by r 1/2 are larger than the molecular gas sizes. Our analysis further provides a determination of the inclination and position angle (PA) of the galaxies. In cases where GalPak 3D did not converge (namely for G0, G1, G6, G20 and G21), we use Sextractor to derive these parameters from the broad band HST/WFPC2 image.
The analysis of galaxy G0 is hampered by the proximity to the bright quasar (2.3 arcsec). G1 shows no indication of a velocity shear notwithstanding the bright emission lines in the object and a favourable inclination parameter (i=51 deg) as measured in the HST/WFPC2 broad band image. Indeed, the nebular emission line ratios in the object show evidence for ionisation by an AGN with Hα emission produced in the central engine (Baldwin, Phillips & Terlevich 1981). G2 is a large galaxy (half-light radius r 1/2 =14±2 kpc) dominated Table 3. CO emitters at z abs =0.313 in the ALMA observations. Flux densities, velocity widths, CO luminosities, molecular masses expressed in units of α CO and depletion timescales are provided for each detected galaxy.  Table 4. Nebular emission lines and cold molecular gas (bold) morpho-kinematic properties. The results of the 3-dimensional MCMC forward modelling of the morphological and kinematical parameters listed are all inclination corrected. Inclination and PA (NoE) are derived from a Sextractor fit to the broad band HST/WFPC2 image when the algorithm did not converge. When the object is detected in both MUSE and ALMA, i.e. G2, G4 and G6, the kinematic models (bold) are sometimes found to be in good agreement indicating a relation between the stellar content and the molecular gas in some of these galaxies . The maximum velocity of the disks are used to derive the dynamical and halo masses. Galaxies G1 to G6 (shown in italics) were previously identified.  (Fig. 1) also shows extended gas around the galaxy. For galaxy G4, we used the Hα emission line rather than the fainter [O iii] 5007Å. The ALMA observations of G4 indicate two, spectrally and spatially well-separated components (Fig. 3) clearly seen in the spectrum (Fig. B2) with wide FWHM=530±50 km s −1 . This is reflected in the kinematics analysis where we derive a large Vmax=290±4 km s −1 for this system. We tentatively interpret these characteristics as a signature of merging and note extended diffuse gas is also present in the MUSE white light image (Fig. B1). Galaxy G6 is almost face-on (i=18 deg) hindering the kinematical analysis of this system. An inspection of the velocity field of G16 reveals that it has yet to converge to its Vmax value, which should be taken as an upper limit. The faint galaxy G17 on the other hand has a well-converged velocity field as illustrated in Fig. B1. Our detection of galaxy G19 indicates that our observations have not fully reached Vmax either, although rotation can be securely reported. Galaxy G20 is a continuum-detected object while G21 has faint emission lines thus precluding further kinematic analysis. Overall, the model converges for most galaxies. While not a unique solution, these findings indicate that the emission lines are well described by disk rotation. However, such an analysis might be rather insensitive to interactions of large objects with lower mass galaxies (say, as in the case of G2 and G17).
We derived the dynamical mass of each galaxy from the enclosed mass: M dyn = V 2 max r 1/2 / G (Péroux et al. When an object is detected from its nebular emission lines and in CO(1-0), i.e. G2, G4 and G6, both the ionised gas traced by nebular lines and the CO molecular line widths (FWHM=200-530 km s −1 ) are significantly broader than the neutral gas traced by absorption at some distance from the galaxies (∼150 km s −1 ) and especially the 21cm line absorption (FWHM=42.1±2.7 km s −1 ). The kinematics of the stars and molecular gas in terms of orientation (PA) deviate by 7-σ (G2) or rather similar (G4). The maximum velocities differ significantly: 264±3 km s −1 (Hα) vs 134±2 km s −1 (CO) for G2 and 231±4 km s −1 (Hα) vs 290±4 km s −1 (CO) for G4. Levy et al. (2018) made a detailed comparison of the molecular and ionised gas (traced by Hα) kinematics in a sample of local galaxies. They find that ∼75% of their sample galaxies have smaller ionised gas rotation velocities than the molecular gas in the outer part of the rotation curve. They report no case where the molecular gas rotation velocity is measurably lower than that of the ionised gas unlike what we observe in the case of galaxy G2.

CONCLUSION
In this paper, we have presented deep MUSE and ALMA observations of the field of the QSO Q1130−1449 which shows a log [N (HI)/cm −2 ]=21.71±0.07 damped Lyα system at z abs = 0.313. Our main findings can be summarised as follow: • Our MUSE observations cover 11 galaxies at the redshift of the absorber, 7 of which are new discoveries. In particular, we report a new object with the smallest impact parameter to the quasar line-of-sight (b=10.6 kpc). Three of the objects are also detected in CO(1-0) in our ALMA observations.
• Using dedicated numerical cosmological simulations, we infer that the physical properties of these galaxies qualitatively resemble a small group environment, possibly part of a filamentary structure.
• Based on metallicity and velocity arguments, we conclude that the neutral gas traced by strong H i absorbers is only partly related to these emitting galaxies while a larger fraction is likely the signature of low surface brightness emitting gas four orders of magnitude fainter than the current detection limits outside the Local Group.
• We report large molecular gas reservoirs with long depletion timescales in the three galaxies detected with ALMA. These results together with other reports of large molecular masses in strong H i-absorption systems indicate that selection based on absorption preferentially picks galaxies which inefficiently convert their gas into stars.
• Detailed kinematics analysis of both the ionised and molecular component of these galaxies shows signatures of past interactions and possible merging between various members of the group. While the stellar component is spatially more extended, the resolved molecular lines are broader in velocity space.
• Together with other similar reports, our findings challenge a picture where strong-N (H I) quasar absorbers are associated with a single bright galaxy and favour a scenario where the H i gas probed in absorption is related to far more complex galaxy structures.  North is up, east to the left. The signal in galaxy G0 panel is dominated by the bright nearby quasar (slightly off-center to the West). While the object lies under the quasar PSF, several emission lines are clearly identified in the MUSE spectrum after quasar spectral-PFS subtraction (see Fig. 5). Some object are bright in continuum but have weak or undetected emission lines (e.g. G20) so that they are apparent in the HST data but weak in MUSE narrow-band. Conversely, some objects show strong emission lines with faint continuum (e.g. G0 and G16).