## An H -absorption-selected cold rotating disk galaxy at i z ≈ 2 193 .

B. Kaur, 1 N. Kanekar, 1 M. Neeleman, 2 M. Rafelski, 3,4 J. X. Prochaska, 5,6 and R. Dutta 7

1 National Centre for Radio Astrophysics, Tata Institute of Fundamental Research, Pune University, Pune 411007, India

2 National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA

3 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

- 4 Department of Physics &amp; Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA

5 Department of Astronomy &amp; Astrophysics, UCO/Lick Observatory, University of California, 1156 High Street, Santa Cruz, CA 95064,

USA

6 Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), 5-1-5 Kashiwanoha, Kashiwa, 277-8583, Japan 7 Inter-University Centre for Astronomy and Astrophysics, Pune University, Pune 411007, India

## ABSTRACT

We have used the Atacama Large Millimeter/submillimeter Array (ALMA) to map CO(3-2) emission from a galaxy, DLA-B1228g, associated with the high-metallicity damped Lymanα absorber at z ≈ 2 1929 towards the QSO PKS B1228-113. . At an angular resolution of ≈ 0 32 . ′′ × 0 24, DLA-B1228g . ′′ shows extended CO(3-2) emission with a deconvolved size of ≈ 0 78 . ′′ × 0 18, . ′′ i.e. a spatial extent of ≈ 6 4 kpc. . We detect extended stellar emission from DLA-B1228g in a Hubble Space Telescope Wide Field Camera 3 F160W image, and find that H α emission is detected in a Very Large Telescope SINFONI image from only one side of the galaxy. While the clumpy nature of the F160W emission and the offset between the kinematic and physical centers of the CO(3-2) emission are consistent with a merger scenario, this appears unlikely due to the lack of strong H α emission, the symmetric double-peaked CO(3-2) line profile, the high molecular gas depletion timescale, and the similar velocity dispersions in the two halves of the CO(3-2) image. Kinematic modelling reveals that the CO(3-2) emission is consistent with arising from an axisymmetric rotating disk, with an exponential profile, a rotation velocity of v rot = 328 ± 7 km s -1 , and a velocity dispersion of σ v = 62 ± 7 km s -1 . The high value of the ratio v rot /σ v , ≈ 5 3, . implies that DLA-B1228g is a rotation-dominated cold disk galaxy, the second case of a highz H -absorption-selected galaxy identified with a cold rotating disk. i We obtain a dynamical mass of M dyn = (1 5 . ± 0 1) . × 10 11 M ⊙ , similar to the molecular gas mass of ≈ 10 11 M ⊙ inferred from earlier CO(1-0) studies; this implies that the galaxy is baryon-dominated in its inner regions.

Keywords: galaxies: evolution -- galaxies: high-redshift - galaxies: ISM

## 1. INTRODUCTION

Damped Ly α absorbers (DLAs), the highest H i column density absorbers in QSO absorption spectra, provide a unique way of identifying gas-rich galaxies at high redshifts (e.g. Wolfe et al. 2005). In the local Universe, DLA column densities, N HI ≥ 2 × 10 20 cm -2 , arise in or around galaxies (e.g. in galaxy disks, high-velocity clouds, intragroup gas, etc); the presence of a DLA along a QSO sightline has hence long been used to identify highz galaxies via the H i absorption signature (e.g. Wolfe et al. 1986). The importance of such 'H -absorption-selected galaxies' stems i from the fact that their selection is not biased towards high-luminosity objects. In contrast, emission-selected galaxy samples implicitly pick out the brighter galaxies of the population. While studies of emission-selected samples have provided a wealth of information on the properties of high-redshift galaxies (e.g., Madau &amp; Dickinson 2014; Tacconi et al. 2020), the luminosity bias inherent in the selection can affect the interpretation of the results. Addressing this issue critically requires the detection and detailed characterization of the H -absorption-selected galaxies associated i with high-redshift DLAs.

For more than three decades, a range of approaches, including broad-band and narrow-band imaging, multi-angle long-slit spectroscopy, and integral field spectroscopy, have been employed at optical and near-IR wavelengths to search

for the galaxies associated with highz DLAs (e.g., Møller &amp; Warren 1993; Warren et al. 2001; Kulkarni et al. 2006; P´ eroux et al. 2012; Fynbo et al. 2010; Bouch´ et al. 2012; Fumagalli et al. 2014; Wang et al. 2015; Krogager et al. e 2017; Lofthouse et al. 2023; Oyarz´n et al. 2024). u Unfortunately, the success of these studies has been limited by the faintness of the stellar emission of the DLA galaxies. Despite the large number of searches, only ≈ 25 H -absorptioni selected galaxies have been identified via these methods at z ≳ 2 (e.g. Krogager et al. 2017), some at very large impact parameters ( &gt; 100 kpc) to the QSO sightline (e.g. Lofthouse et al. 2023).

In the last few years, mm- and sub-mm searches for the galaxy counterparts of highz DLAs with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimetre Array (NOEMA) have changed the field. ALMA and NOEMA detections of CO and [C ii ] 158 µ m emission in the fields of high-metallicity DLAs have yielded the identification of more than 15 H -absorption-selected galaxies out to i z ≈ 4 3 (e.g. Neeleman et al. 2017, . 2018, 2019; Kanekar et al. 2020; Kaur et al. 2022b; Combes &amp; Gupta 2023). At z ≈ 2, the galaxies identified via CO emission have been found to have relatively low impact parameters ( b ≈ 5 -30 kpc, with a median value of ≈ 11 kpc) to the QSO sightline, or to have a companion at a lower impact parameter; this suggests that the Ly α absorption typically arises in the interstellar medium (ISM) of the DLA galaxy (Neeleman et al. 2018; Fynbo et al. 2018; Kanekar et al. 2020; Kaur et al. 2022b). Conversely, at z ≈ 4, objects identified using [C ii ] 158 µ m emission have higher impact parameters, b ≈ 15 -60 kpc (with a median value of ≈ 27 kpc), suggesting the presence of high-H i column density gas in the circumgalactic medium (CGM) of galaxies at these redshifts (Neeleman et al. 2017, 2019). Surprisingly, a number of the newly-identified H -absorption-selected galaxies, associated with high-metallicity DLAs, have been i found to have high molecular gas masses, ≳ 5 × 10 10 M ⊙ (e.g. Neeleman et al. 2018, 2020; Kanekar et al. 2020; Combes &amp; Gupta 2023). Follow-up ALMA and Karl G. Jansky Very Large Array (VLA) studies have revealed that these highz DLA galaxies also have high, near-thermal, excitation of the mid-J CO rotational levels (Klitsch et al. 2022; Kaur et al. 2022a). This suggests that a significant fraction of high-metallicity DLAs at z ≈ 2 arise in the ISM or CGM of massive, dusty, star-forming galaxies.

Recently, H α , CO, and [C ii ] 158 µ m mapping studies, mostly of emission-selected galaxy samples, have revealed the presence of rotating disk galaxies out to z ≈ 4 5 (e.g. Bouch´ et al. 2013; Genzel et al. 2017; Neeleman et al. 2020; Rizzo . e et al. 2020; Lelli et al. 2021). Interestingly, some of these highz rotating disks are 'cold', i.e. their kinematics are dominated by rotation. In the case of H -absorption-selected galaxies, ALMA mapping studies have so far been carried i out for two systems at z ≈ 4, using the [C ii ] 158 µ m line, finding a major merger for the z ≈ 3 7978 DLA towards . J1201+2117 (Prochaska et al. 2019) and a cold rotating disk galaxy for the z ≈ 4 2546 DLA towards J0817+1351 . (the 'Wolfe disk'; Neeleman et al. 2020). In this Letter , we present the first case of high-resolution CO mapping of a high-redshift H -absorption-selected galaxy, at i z ≈ 2 193 towards the quasar PKS B1228-113. . 1

## 2. ALMA OBSERVATIONS AND DATA ANALYSIS

We used the ALMA Band-3 receivers to observe the redshifted CO(3-2) line from the z ≈ 2 193 DLA galaxy towards . PKS B1228-113. This DLA field was originally searched for CO(3-2) emission by Neeleman et al. (2018), as part of a sample of high-metallicity DLAs (Kanekar et al. 2020), and yielded a detection of strong CO(3-2) emission at low angular resolution (2 . ′′ 6 × 1 7). . ′′ The new ALMA observations were carried out in July 2021 (ALMA proposal ID 2019.1.00373.S; PI M. Neeleman), using the extended C43-7 configuration (longest baseline: 3.6 km) with an onsource integration time of 4.1 hours. PKS B1228-113 was chosen as the pointing centre as it is bright at mm wavelengths, allowing self-calibration of the target field. One of the four ALMA basebands was centred on the redshifted CO(3-2) line frequency of the z = 2 1929 DLA ( . ν obs = 108 291 GHz; Neeleman et al. 2018), with a bandwidth of 1.875 GHz, . 960 channels, and the frequency division mode, yielding a frequency resolution of ≈ 1 95 MHz. . The remaining three basebands were set up in the time division mode, each with a bandwidth of 2.0 GHz and 128 channels, giving a frequency resolution of 15.625 MHz.

Weused the standard ALMA pipeline in the Common Astronomy Software Application package ( casa version 6.5.4.9; CASA Team 2022; Hunter et al. 2023) to carry out the initial data editing and calibration of the visibility dataset. After inspecting the calibrated visibilities and performing manual flagging, we averaged the data to a coarse spectral resolution for the purpose of self-calibration. The self-calibration process involved three rounds of continuum imaging and phase-only self-calibration, followed by one round of amplitude-and-phase self-calibration and imaging, after which no further improvement was seen in either the image or the image/visibility residuals upon further self-calibration.

Table 1. Properties of the spectral cubes, at a channel resolution of 35 km s -1 , with different values of the robust parameter. The columns are (1) the value of the Briggs robust parameter, (2) the synthesized beam of the resulting cube, (3) the RMS noise per channel, in mJy Beam -1 , and (4) the signal-to-noise ratio (S/N) of the peak CO emission in the velocity-integrated CO(3-2) flux density image.

|   Robust | Synthesized beam      |   RMS noise (mJy) |   Peak S/N |
|----------|-----------------------|-------------------|------------|
|      1   | 0 . 32 ′′ × 0 . 24 ′′ |              0.15 |        8.1 |
|      0.5 | 0 . 27 ′′ × 0 . 21 ′′ |              0.16 |        6.7 |
|      0   | 0 . 23 ′′ × 0 . 17 ′′ |              0.19 |        4.7 |
|     -0.5 | 0 . 21 ′′ × 0 . 17 ′′ |              0.23 |        4.7 |
|     -1   | 0 . 20 ′′ × 0 . 16 ′′ |              0.27 |        3.9 |

In the self-calibration procedure, the antenna-based gains were measured using the routine gaincalR (Chowdhury et al. 2020). The continuum imaging in the self-calibration procedure used Briggs weighting (Briggs 1995) with a robust value of -0 5. . For the final continuum image, after self-calibration, we used a robust value of +0 5, obtaining . a synthesized beam of 0 29 . ′′ × 0 22 and an RMS noise of 15 . ′′ µ Jy beam -1 . PKS B1228-113 has a flux density of 24 85 . ± 0 02 mJy in this image. .

We next applied the antenna-based gains obtained from the self-calibration procedure to the CO(3-2) visibility dataset at the original spectral resolution, and subtracted out the continuum emission from the calibrated visibilities. We then used the task tclean to make spectral cubes from the residual visibilities, using a range of velocity resolutions (25, 35, 50 km s -1 ) and values of the robust parameter (-1, -0.5, 0, +0.5, +1). For all cubes, the CO(3-2) line emission was cleaned down to a threshold of 0.5 σ , where σ is the per-channel RMS noise. For cleaning, we used a square mask of size 2 5 . ′′ × 2 5 . ′′ , enclosing the CO emission for all frequency planes. Table 1 summarizes the details of the spectral cubes for the different values of the Briggs robust parameter. We chose a robust value of +1 and a velocity resolution of 35 km s -1 for the final spectral cube used for the kinematic modelling, based on a balance between the per-channel sensitivity and signal-to-noise ratio, and the angular resolution. The final spectral cube has a synthesized beam of 0 32 . ′′ × 0 24, a position angle of = . ′′ -74 4 . ◦ , and an RMS noise of 0 15 mJy beam . -1 at a velocity resolution of 35 km s -1 .

## 3. HST OBSERVATIONS AND DATA ANALYSIS

The field of PKS B1228-113 was observed in 2019 (PID: 15882; PI N. Kanekar) and 2024 (PID: 17483; PI R. Dutta) using the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope , with the F105W and F160W filters 2 respectively. For the z ≈ 2 193 DLA galaxy, this corresponds to rest-frame central wavelengths of . ≈ 3000 ˚ and A ≈ 5000 ˚ , respectively. A A single orbit was obtained with each filter, with a larger dither pattern used to achieve cleaner images by dithering over the size of the IR blobs. The F105W observations use the WIDE-7 dither pattern increased by a factor of 3 over the standard size (ISR 2016-14; Anderson 2016), while the F160W observations used a standard three-point line dither pattern with the size of the dither pattern increased by a factor of 5.

The individual images were combined with AstroDrizzle (Hack et al. 2020) to a scale of 0 . ′′ 06 pixel -1 , and the images astronomically aligned to GAIA DR3 (Gaia Collaboration et al. 2023) with TweakReg, with an uncertainty of ≈ 0 02. . ′′ The effective angular resolution of the HST WFC3 images is ≈ 0 3, based on Gaussian fits to the point spread function . ′′ (PSF) of the quasar.

## 4. RESULTS

Figure 1 shows the velocity-integrated image of the CO(3-2) emission (integrated over the velocity range of -350 km s -1 to +525 km s -1 ) detected in our ALMA spectral cube, at an angular resolution of ≈ 0 32 . ′′ × 0 24 . ′′ (the corresponding spatial resolution at z = 2 193 is 2 6 kpc . . × 2 0 kpc). . The QSO position is indicated by the magenta star, ≈ 3 9 (i.e. . ′′ ≈ 33 kpc) to the south of the DLA galaxy (Neeleman et al. 2018). The CO(3-2) emission has an integrated line flux density of (0 73 . ± 0 12) . Jy km s -1 . This was obtained using the casa task imfit to fit a single-component 2-D Gaussian model to the image of Figure 1; this model was found to yield an acceptable fit, with residuals consistent with noise. The measured integrated CO(3-2) line flux density is consistent with the value of

Figure 1. The velocity-integrated CO(3-2) image, with contours at ( -3 3 4 2 6) , , . , × σ , where σ is the RMS noise on the image. The magenta star indicates the QSO position. The right panel shows a zoom-in on DLA-B1228g. The x-axis and y-axis are right ascension and declination relative to the kinematic centre of the CO(3-2) emission (see Table 2), in arcseconds. The synthesized beam is shown in the bottom-left corner of each panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00850v1-figure-0864ce71-cf49-40f2-a638-afa28ce55a9c.jpg)

(0 726 . ± 0 031) Jy km s . -1 obtained from the original low-resolution ALMA CO(3-2) spectrum (Neeleman et al. 2018; Kanekar et al. 2020), indicating that no emission is being resolved out in our high-resolution ALMA cube. The 2-D Gaussian fit to the CO(3-2) emission yields a deconvolved size of (0 . ′′ 78 ± 0 15) . ′′ × (0 . ′′ 18 ± 0 07), i.e. . ′′ a physical size of ≈ 6 4 kpc . × 1 5 kpc. . No continuum emission is detected at the location of DLA-B1228g in the 0 29 . ′′ × 0 22 resolution . ′′ continuum image.

Figure 2[A] shows the mean CO(3-2) velocity field which was obtained by fitting a Gaussian function to the spectrum at each pixel of the cube, and then using the mean value of the function at each pixel, via the fitting package qubefit (Neeleman et al. 2021); the image has been blanked at spatial locations where the integrated CO(3-2) flux density has &lt; 3 σ significance. The velocity field appears consistent with that of a rotating disk, as also suggested by the double-horned profile seen in the high-resolution ALMA spectrum. Figure 2[B] shows a position-velocity (PV) slice taken along the major axis 3 of the galaxy disk (indicated by the dashed black line in Figure 2[A]). The velocities along the major axis are seen to increase with increasing distance from the galaxy centre, indicative of emission from a rotating disk galaxy. However, we note that the current ALMA CO(3-2) image has intermediate angular resolution and sensitivity, with only 2 -3 synthesized beams across the CO-emitting region. In such situations, it is possible for a merger system to show a velocity field similar to that of a rotating disk (e.g. Rizzo et al. 2022). In Section 5.5, we will discuss the possibility that DLA-B1228g might be a merger system, rather than a disk galaxy.

Fig 3[A] shows the HST-WFC3 F160W image (in colourscale) overlaid on the ALMA CO(3-2) integrated line flux density image (in black contours). The F160W image shows extended emission (detected at ≈ 8 σ significance) that is approximately coincident with the location of the CO(3-2) emission. The flux-weighted centroid of the HST-WFC3 F160W emission is at RA=12h30m55.510s, Dec.=-11 39 06.48 ◦ ′ ′′ , at an offset of ≈ 0 2 (i.e. . ′′ ≈ 1 5 . kpc) from the kinematic centre of the CO(3-2) emission (see Table 2 and Section 5). The source has an F160W AB magnitude of 24 5 . ± 0 1, . measured over an elliptical aperture with the Kron radius (Kron 1980). We find that the HST F160W emission, corresponding to rest-frame 5000 ˚ , appears clumpy in nature, with no evidence of a disk-like morphology. A

## 5. DISCUSSION

5.1. Stellar properties of DLA-B1228g

3 We note that the major axis was obtained from the kinematic modelling presented in Section 5. Similar results are obtained on taking a PV slice along the major axis of the 2D Gaussian fit to the integrated CO(3-2) intensity.

Relative Declination (arcsec)

0.6

0.3

0.0

-0.3

-0.6

[A]

0.6

-0.6

-0.3

0.0

0.3

Relative Right Acension (arcsec)

Figure 2. [A] The CO(3-2) mean velocity field, after blanking all spatial locations where the integrated CO(3-2) flux density has &lt; 3 σ significance. The velocity field is seen to be similar to that of a rotating disk. The red cross indicates the kinematic centre of the CO(3-2) emission (see Table 2). [B] A PV slice along the major axis of the galaxy, indicated by the dashed line in the left panel. The y-axis contains velocity, in km s -1 , relative to the DLA redshift, z = 2 1929, while the x-axis is distance . from the the peak of the CO(3-2) emission (see Table 2), in arcseconds. The flux density contours start at ± 2 σ and increase in steps of √ 2, where σ is the RMS noise per 35 km s -1 channel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00850v1-figure-80bbdb9b-316a-48c4-99e3-7b1c83852576.jpg)

Figure 3. ALMA CO(3-2) integrated line flux density image (black contours; both panels) overlaid on [A] the HST WFC3 F160W image (colourscale) and [B] the HST WFC3 F105W image (colourscale). The images also show the VLT SINFONI integrated H α intensity image (magenta contours starting at ± 3 σ and increasing in steps of √ 2). The x-axis and y-axis are, respectively, right ascension and declination, relative to the kinematic centre of the CO(3-2) emission (see Table 2). The inset on the bottom left shows the ALMA synthesized beam. In both panels, the yellow cross indicates the possible NUV counterpart of DLA-B1228g identified in the HST-WFC3 F105W image by Kaur et al. (2022a). The new ALMA CO(3-2) image shows that this object is clearly offset from the CO emission, and is thus unlikely to be related to DLA-B1228g. See main text for discussion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00850v1-figure-e189f942-3133-4e4b-a403-a5514d6d86fd.jpg)

Kaur et al. (2022a) presented an HST-WFC3 image of the field of DLA-B1228g in the F105W filter, corresponding to rest-frame near-ultraviolet (NUV) wavelengths at the galaxy redshift of z ≈ 2 193. This yielded the apparent detection . of the stellar emission of DLA-B1228g, at a position consistent with that of the CO emission in the low-resolution ALMA CO(3-2) and VLA CO(1-0) images (Neeleman et al. 2018; Kaur et al. 2022a). Figure 3[B] shows the new

200

100

0

100

200

)

1

Velocity (km s

Table 2. Results for the kinematic modelling of DLA-B1228g

| Parameter estimates for the rotating thin-disk model   | Parameter estimates for the rotating thin-disk model   |
|--------------------------------------------------------|--------------------------------------------------------|
| Right Ascension                                        | 12 h 30 m 55 . 505( ± 1) s                             |
| Declination                                            | - 11 ◦ 39 ′ 06 . 3870( +5 - 2 ) ′′                     |
| Galaxy redshift                                        | 2 . 19326 ± 0 . 00007                                  |
| Inclination ( i , in deg)                              | 66 . 4 ± 2 . 2                                         |
| Position angle (PA, in deg)                            | 305 . 2 ± 1 . 4                                        |
| Peak intensity ( I 0 , in mJy/beam)                    | 3 . 9 +0 . 7 - 0 . 5                                   |
| Scale length (R d , in kpc)                            | 1 . 8 ± 0 . 1                                          |
| Rotation Velocity ( v rot , in km s - 1 )              | 328 ± 7                                                |
| Velocity Dispersion ( σ v , in km s - 1 )              | 62 ± 7                                                 |

Note -The coordinates indicate the kinematic center of the CO(3-2) emission, obtained from the thin disk modelling of DLA-B1228g.

high-resolution ALMA CO(3-2) image (black contours) overlaid on the HST-WFC3 F105W image (colourscale) and a Very Large Telescope (VLT) SINFONI H α image (magenta contours; Neeleman et al. 2018). The object identified in the HST-WFC3 F105W image by Kaur et al. (2022a) is at RA=12 h 30 m 55 44 , Dec.= . s -11 39 05 87 ◦ ′ . ′′ (shown by the yellow cross in Figure 3[B]), and is offset by ≈ 1 02 . ′′ , i.e. ≈ 8 4 kpc, from the peak of the CO(3-2) emission. . This object is also clearly detected in our new F160W image of Fig. 3[A]. The offset between this HST emission feature and the CO(3-2) emission is much larger than the typical spatial offsets of ≲ 4 kpc that have been observed between the stellar and the gas emission in highz dusty star-forming galaxies (e.g. Hodge et al. 2015; Calistro Rivera et al. 2018; Bik et al. 2023). Further, this HST-WFC3 object shows no H α emission at z ≈ 2 193 in the VLT-SINFONI image, . with an upper limit of ≈ 0 9 M . ⊙ yr -1 on its star-formation rate (SFR; P´roux et al. 2011). e This, along with the large spatial offset between the optical and the CO emission, indicates that the HST-WFC3 105W source is unlikely to be directly associated with the CO-emitting DLA galaxy.

Of course, as noted in the previous section, we detect extended emission in the new HST-WFC3 F160W image of Fig. 3[A] which is nearly spatially coincident (offset ≈ 1 4 kpc) with the CO(3-2) emission, and which is not seen . in the HST-WFC3 F105W image. This extended emission, at a rest-frame wavelength of ≈ 5000 ˚ at A z ≈ 2 193, . appears to be the stellar counterpart of the CO(3-2) emission. Dust obscuration, perhaps allied with the Balmer break, is likely to be the main cause of the non-detection of this emission in the F105W image. The F160W emission area is 112 pixels, corresponding to a circular radius of 0.2 arcsec. The 2 σ upper limit on the HST-WFC3 F105W non-detection, at a rest-frame wavelength of ≈ 3000 ˚ A over this area is AB mag of 26.5. This yields an upper limit on the unobscured SFR of ≈ 1 3 M . ⊙ yr -1 , assuming a Chabrier initial mass function (Chabrier 2003) and using the local relation between NUV luminosity and the SFR (Kennicutt &amp; Evans 2012). This is significantly lower than the total SFR of the galaxy measured from the sub-mm continuum, ≈ 100 M ⊙ yr -1 (Neeleman et al. 2018; Klitsch et al. 2022), indicating a high dust obscuration in the DLA galaxy.

Finally, the two panels of Figure 3 show the VLT SINFONI H α image, which displays emission (with &gt; 4 σ statistical significance; Neeleman et al. 2018) at the south-eastern edge of the CO(3-2) emission. Figure 3 of Neeleman et al. (2018) shows that the putative H α emission arises at negative velocities, in good agreement with the CO velocities on the south-eastern side of the disk. It thus appears that the H α emission is detected only from one side of the galaxy. Interestingly, the new HST-WFC3 F160W image shows rest-frame ≈ 5000 ˚ emission at the location of the VLTA SINFONI H α emission. However, H α emission is not detected from the peak of the F160W emission, closest to the centre of DLA-B1228g. This suggests that the brightest F160W emission may arise from an older stellar population, that is not currently undergoing active star formation.

Figure 4. The three rows show the CO(3-2) integrated line flux density (Jy km s -1 ), the CO(3-2) velocity field (km s -1 ), and the CO(3-2) velocity dispersion (km s -1 ), while the three columns show the ALMA image, the best-fit rotating-disk model image, and the residual image obtained after subtracting the model image from the corresponding ALMA image. The x-axis and y-axis are, respectively, right ascension and declination (in arcseconds), relative to the peak of the CO(3-2) emission measured from the best-fit rotating-disk model (see Table 2). In the top left panel, the contours start at ± 3 σ , and increase in steps of √ 2; negative features are shown with dashed contours. We note that the top right panel, showing the residual integrated line flux density image, does not contain any residual emission with ≥ 3 σ or ≤ -3 σ significance, and hence does not have any marked contours. The images of the middle and bottom panels have been blanked at spatial locations where the integrated CO(3-2) flux density has &lt; σ 3 significance. The inset in the left panels shows the ALMA synthesized beam.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00850v1-figure-c99001a5-63f8-4cdf-8140-3a7366898328.jpg)

The CO(3-2) velocity field of DLA-B1228g and the PV slice along the galaxy major axis both indicate that the galaxy is likely to be a rotating disk. We used the Bayesian Markov Chain Monte Carlo (MCMC)-based Python fitting package qubefit (Neeleman et al. 2021) to model the kinematics of the CO(3-2) emission. The fitting routine convolves a user-defined model with the synthesized beam of the observations and then compares the resultant model cube to the actual spectral cube, using a defined likelihood function. It then maximizes the likelihood function, following an MCMC approach in order to estimate the model parameters and their uncertainties. We assumed a

rotating infinitesimally-thin exponential disk model for DLA-B1228g, in which the intensity I of the CO emission at a radius r is given by I ( r ) = I 0 e -r/r d , where I 0 is the central intensity and r d is the disk scale length. The model is described by the following nine parameters: the coordinates of the centre of the CO emission, the galaxy redshift, the galaxy inclination ( ), the position angle (PA) of the major axis, the rotation velocity ( i v rot ), the velocity dispersion ( σ v ), the disk scale length ( r d ), and the peak intensity of the CO(3-2) emission ( I 0 ). For simplicity, we assume that both the rotational velocity and the velocity dispersion are constant across the disk. Finally, we used 400,000 MCMC chains to estimate the probability distribution functions of the nine parameters of the model.

The results of the kinematic modelling with qubefit are shown in Figure 4, where the top row shows the integrated CO(3-2) line flux density, the middle row shows the CO(3-2) velocity field, and the bottom row shows the CO(3-2) velocity dispersion. The three columns contain the data, the model, and the residuals after subtracting the model from the data. It is clear from the top right panel that no residual CO(3-2) emission is detected at a statistically significant ( ≥ 3 σ ) level, implying that the assumed thin rotating disk model provides a good fit to the data. The reduced chi-square of the fit is 1.2, consistent with a good fit.

The best-fit parameters of the thin-disk model are listed in Table 2; note that the errors have been increased appropriately to account for the fact that the reduced chi-square of the fit is slightly higher than unity. We obtain a rotation velocity of v rot = 328 ± 7 km s -1 , and a velocity dispersion of σ v = 62 ± 7 km s -1 .

The ratio of rotation velocity to velocity dispersion provides information on whether a galaxy's kinematics is dominated by rotation ( v rot /σ v &gt; 1) or turbulence ( v rot /σ v &lt; 1). For DLA-B1228g, the kinematic modelling yields the high ratio v rot /σ v = 5 3 . ± 0 6, . indicating that the kinematics is dominated by rotation. This is similar to values of v rot /σ v that have recently been obtained for emission-selected galaxies at similar redshifts, when cold gas tracers such as CO, [C ii ] 158 µ m, or C i emission are used to map the kinematics (e.g. Lelli et al. 2023; Rizzo et al. 2023)

In passing, we note that we used a simple axisymmetric exponential disk model, with a constant rotational velocity and velocity dispersion, for the kinematic modelling as the signal-to-noise ratio of the detected CO(3-2) emission is not very high. Given that this model already provides a good fit, more complicated models (with more parameters) are also likely to yield good fits to the observed CO(3-2) emission.

## 5.3. The dynamical mass of DLA-B1228g

Our kinematic modelling of the CO(3-2) emission of DLA-B1228g allows us to estimate the dynamical mass of the galaxy, following standard procedures. For a spherically symmetric mass distribution, the dynamical mass enclosed within a radius R is given by M dyn ( R ) = v 2 × R/G = 2 33 . × 10 5 × v 2 × R M ⊙ (e.g. Neeleman et al. 2021), where v is the circular velocity, in km s -1 , and R is in kpc.

The circular velocity is related to the rotation velocity and the velocity dispersion by v 2 circ = ( v 2 rot + ησ 2 v ), where the value of η lies in the range ≈ 2 -6 depending on the mass distribution and the galaxy kinematics (e.g. Epinat et al. 2009; Burkert et al. 2010). Since DLA-B1228g is a rotation-dominated system, with v rot /σ v ≈ 5 3, we will neglect the . velocity dispersion in the above expression, and will assume v circ ≈ v rot = 328 ± 7 km s -1 .

Finally, following Neeleman et al. (2021), we estimate the dynamical mass of DLA-B1228g within a radius equal to twice the effective radius ( R e ) of the CO(3-2) emission. For an exponential disk, this radius contains ≈ 85% of the total emission. Further, for such a disk, the effective radius (i.e. the half-light radius) is equal to 1 678 . × r d , where r d is the disk scale length. For DLA-B1228g, r d = 1 8 . ± 0 1 kpc, implying an effective radius of . R e = 3 02 . ± 0 17 kpc. We . note that the effective diameter is similar to the estimate of the deconvolved size of the CO(3-2) emission, 6.4 kpc, as measured from the 2-D Gaussian model. The radius that we will use to estimate the dynamical mass is then R = 2 R e = 6 04 . ± 0 34 kpc. .

We obtain a dynamical mass of (1 5 . ± 0 1) . × 10 11 M ⊙ for DLA-B1228g, within the radius R = 6 04 kpc. We emphasize . that the quoted error only includes the error from the fit. However, the uncertainties in the estimate are dominated by the assumption of a spherical mass distribution. If the actual mass distribution is intermediate between a thin disk and a sphere (e.g. if the baryons are mostly in a disk and the dark matter in a sphere), this estimate would be higher than the true dynamical mass, by ≈ 20% (e.g. Neeleman et al. 2021).

The dynamical mass estimate includes contributions from the dark matter and the baryons (M bary ≈ M mol +M HI + M ) of the galaxy, within a radius of 2 ∗ R e . While the baryonic mass of highz galaxies has been shown to be dominated by the atomic gas (e.g. Chowdhury et al. 2022a), the latter is extended on scales of tens of kpc (e.g. Chowdhury et al. 2022b). The baryonic mass of the inner disk is thus likely to be mainly contributed by the molecular gas and the stars. For DLA-B1228g, Kaur et al. (2022a) used VLA CO(1-0) spectroscopy to obtain a molecular gas mass of

Figure 5. The left panel [A] shows the expected H i disk of DLA-B1228g (large ellipse) and its expected velocity field, based on the velocity field of the ALMA CO(3-2) emission (inner ellipse) and an assumed flat rotation curve (see main text for a description of the model for the H i disk). The magenta star indicates the location of the sightline towards the background QSO. The x and y axes are in kpc, relative to the kinematic centre of the CO(3-2) emission (see Table 2). The top right panel [B] shows the VLT-UVES Ni ii λ 1751 ˚ absorption profile from the A z ≈ 2 1929 DLA, while the bottom right panel [C] shows the . integrated ALMA CO(3-2) flux density spectrum of DLA-B1228g, with the dashed line showing the ± 1 σ error, where σ is the rms noise per 35 km s -1 velocity resolution. The vertical blue line in panel [B] shows the expected absorption velocity if the absorption arises in an extended H i disk with a flat rotation curve (see main text for discussion). In both the right panels, the x-axis is velocity, in km s -1 , relative to the redshift of DLA-B1228g ( z = 2 1933; see Table 2). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00850v1-figure-27674793-548a-4d02-bed5-d080303dcd78.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00850v1-figure-9f67d971-548d-405e-9f44-f9091ed81d15.jpg)

M mol = (1 0 . ± 0 2) . × 10 11 × ( α CO / 4 36) M . ⊙ . We do not currently have an estimate of the stellar mass of DLA-B1228g. The total baryonic mass within 2 R e is thus ≳ 1 0 . × 10 11 M , using the molecular gas mass estimate alone. ⊙ This is comparable to the dynamical mass estimate of ≈ 1 5 . × 10 11 M ⊙ (especially given that the dynamical mass is likely to be over-estimated by ≈ 20%). We thus find that the central regions of DLA-B1228g are baryon-dominated, as has also been found for emission-selected galaxies at similar redshifts (e.g, Genzel et al. 2017, 2020)

## 5.4. The z ≈ 2 1929 . DLA: Disk or CGM?

An important open question in studies of galaxy evolution is whether the high H i column densities characteristic of DLAs arise primarily in the disks of highz galaxies (i.e. in their ISM) or are also commonly found in their CGM. Indeed, one of the main reasons for the choice of the DLA-defining threshold column density of 2 × 10 20 cm -2 is that such high H i column densities are characteristic of galaxy disks in the local Universe and thus provide a signpost for the presence of a gas-rich galaxy along the QSO sightline (Wolfe et al. 1986).

Very interestingly, the H -absorption-selected galaxies associated with DLAs at i z ≳ 4 have impact parameters of ≈ 15 -60 kpc (Neeleman et al. 2017, 2019), with no evidence of a closer companion galaxy found in HST imaging (Kaur et al. 2021). While the number of identified H -absorption-selected galaxies at i z ≳ 4 remains small today, the observational data suggest that DLAs at these redshifts might arise primarily in the CGM.

Conversely, at lower redshifts, z ≈ 2, most of the presently identified H -absorption-selected galaxies have lower i impact parameters, ≲ 15 kpc (e.g. Krogager et al. 2017; Neeleman et al. 2018; Kanekar et al. 2020; Kaur et al. 2022a), consistent with the DLA arising in the ISM of the associated galaxy. We note that some of the optical techniques (e.g. Møller et al. 2004; Fynbo et al. 2010; Wang et al. 2015; Krogager et al. 2017) are biased towards identifying DLA galaxies at low impact parameters.

In this context, the large impact parameter ( ≈ 33 kpc, similar to that of DLA galaxies at z ≈ 4) of DLA-B1228g to the sightline towards PKS B1228-113 is unusual. However, this does not immediately imply that the DLA arises in the CGM of the galaxy, because the H i disk of a galaxy is typically far more extended than the stellar or molecular disks. In the local Universe, the H i diameter of a galaxy of a given H i mass can be obtained from the H i mass-size relation (e.g. Broeils &amp; Rhee 1997; Wang et al. 2016). For example, a galaxy at z ≈ 0 with an H i mass of ≈ 10 10 M ⊙ has a typical H i diameter (at an H i surface density of 1 M ⊙ pc -2 ) of ≈ 60 kpc (Wang et al. 2016). At high redshifts, we do not have measurements of the H i masses of individual galaxies, due to the weakness of the H i 21cm emission line.

However, recent H 21cm stacking studies have shown that the average H i i mass of star-forming galaxies at z ≈ 1 3 is . significantly larger than their molecular gas and stellar masses (Chowdhury et al. 2022a). In fact, atomic gas has been shown to make up ≈ 70% of the baryonic mass of star-forming galaxies at z ≈ 1 3 (Chowdhury et al. 2022a). . This indicates that the atomic gas mass of DLA-B1228g is likely to be larger than the molecular gas mass, i.e. that the H i mass is ≳ 7 3 . × 10 10 × ( α CO / 4 36) M . ⊙ (where we have used a factor of 1.36 to convert the atomic gas mass to the H i mass).

Further, Chowdhury et al. (2022b) show that the stacked H 21cm emission signal from star-forming galaxies at i z ≈ 1 is spatially resolved for synthesized beams &lt; 90 kpc, implying an average H i disk diameter ≳ 50 kpc for galaxies at z ≈ 1 3 with an H . i mass of ≈ 1 4 . × 10 10 M . This is consistent with the expected H ⊙ i mass-size relation at z ≈ 0 (Wang et al. 2016). Assuming that DLA-B1228g follows the local H i mass-size relation and that its H i mass is ≳ 7 3 . × 10 10 M ⊙ yields an H i disk diameter of ≳ 160 kpc. The solid ellipse in Fig. 5[A] shows the expected H i disk geometry for DLA-B1228g, assuming that the H i disk follows the same geometry as the molecular gas observed in CO emission, for an assumed H i mass of ≈ 7 3 . × 10 10 M ⊙ and the local H i mass-size relation. It is clear in this scenario that the sightline to PKS B1228-113 is expected to be covered by the H i disk, indicating that it is possible that the damped Lymanα absorption arises in the H i disk of the galaxy.

Fig. 5[A] also shows the expected H i velocity field of DLA-B1228g, assuming that the H i disk too can be described as a thin exponential disk, with the flat rotation curve of the molecular disk. The expected velocity from the H i disk at the location of the QSO sightline is seen to be ≈ -123 km s -1 . Fig. 5[B] shows the Ni ii λ 1751 ˚ absorption A line detected in a VLT Ultraviolet Echelle Spectrograph (UVES) spectrum from the z = 2 1929 DLA (top panel; . Akerman et al. 2005; Neeleman et al. 2018) and the integrated ALMA CO(3-2) spectrum (bottom panel; this work). Note that zero velocity in this figure is at the systemic redshift of DLA-B1228g, i.e. at z = 2 19326. . It is clear that the bulk of the low-ionization metal line absorption lies redward of the expected velocity of ≈ -123 km s -1 from the H i disk, extending to a velocity of ≈ +100 km s -1 . A simple thin disk model thus fails to explain the observed low-ionization metal line absorption in the z ≈ 2 1929 DLA. Scenarios that explain the observed velocity structure of . the low-ionization metal lines include absorption in outflowing or accreting gas, in gas clumps in the CGM, or in a warp or flare in the outer disk of the galaxy (e.g. Bouch´ et al. 2013). e

## 5.5. Can the CO Kinematics be explained by a Major Merger?

Our kinematic modelling of the CO(3-2) emission has shown that a rotating disk model provides a good fit to the present CO data. However, as emphasized by Rizzo et al. (2022), it is possible for the velocity field of a merger system to appear to be a rotating disk, for observations (such as the present one) with intermediate angular resolution and sensitivity. We hence consider the possibility that DLA-B1228g is actually a merger system, rather than a disk galaxy.

At the outset, we note that the apparently clumpy nature of the HST F160W emission does appear suggestive of a merger. Further, the kinematic centre obtained from the modelling is offset from the physical centre of the CO emission. Conversely, the double-peaked CO(3-2) emission spectrum of Fig. 5[C] is symmetric, which would be surprising for a merger. Further, the merger would have to be of two galaxies of similar CO(3-2) line flux densities and velocity widths, i.e. of approximately similar molecular gas masses, implying a major merger. One of the putative merging galaxies would then give rise to the blueward CO emission of Fig. 2, and the other, the redward CO emission, with a physical separation of ≈ 3 kpc and a velocity separation of ≈ 500 km s -1 . Such a major merger would be very likely to trigger massive star formation, producing strong H α emission (e.g, Alaghband-Zadeh et al. 2012). However, (weak) H α emission is seen from only one of the regions showing F160W emission, with no H α emission seen from the second region with the F160W emission, closer to the centre of the CO emission. This too would be very surprising for a major merger. We note that dust obscuration would affect the F160W rest-frame 5000 ˚ emission more than A the H α emission, and thus cannot explain the non-detection of H α emission from the region with F160W emission. In addition, the velocity dispersions of the blue and red velocity components of the CO emission are seen to be very similar in the bottom-left panel of Fig. 5. In the major-merger scenario, this would be surprising as the velocity dispersions of the two merging galaxies would be expected to be uncorrelated. Further, the DLA sightline lies to the south of DLA-B1228g, i.e. on the side of the blueward CO emission. However, Fig. 5[B] shows that the low-ionization metal line absorption arises roughly midway between the velocities of the two CO peaks. This would be surprising in a merger scenario, as one would expect the absorption to arise from gas that is associated with the galaxy producing the blueward CO emission. Finally, combining the SFR estimate of ≈ 100 M ⊙ yr -1 (Neeleman et al. 2018; Klitsch et al. 2022) with the molecular gas mass of ≈ 10 11 M ⊙ (Kaur et al. 2022b) yields a molecular gas depletion timescale

of ≈ 1 Gyr for DLA-B1228g, far larger than typical molecular gas depletion timescales of high-redshift starburst galaxies ( ≲ 0 1 Gyr; . Tacconi et al. 2020). Assuming a significantly lower value of the CO-to-H 2 conversion factor, α CO ≈ 1 M ⊙ (K km s -1 pc ) 2 -1 (e.g. Bolatto et al. 2013), would reduce this discrepancy, but would still yield a molecular gas depletion timescale of ≈ 0 25 Gyr, higher than that in typical major mergers at . z ≈ 2. Overall, it appears very unlikely that DLA-B1228g is a major merger system, although a definitive test of this scenario would require a high-resolution, high-sensitivity CO imaging study (e.g. Rizzo et al. 2022).

In passing, we note that the DLA metallicity of -0 22 (Akerman et al. 2005; Neeleman et al. 2018) appears quite . high, given that the DLA sightline has a high impact parameter, ≈ 33 kpc. However, studies of metallicity gradients in large samples of massive emission-selected galaxies at z ≈ 2 have mostly yielded flat or only slightly negative metallicity gradients, out to radii of ≈ 10 kpc, and even mildly positive gradients in some galaxies (e.g. Wuyts et al. 2016; F¨rster o Schreiber et al. 2018). In the case of absorption-selected galaxies, Christensen et al. (2014) compared emission and absorption metallicities for two DLAs at z ≈ 2 3 . -2 6 and their associated galaxies (at impact parameters of . ≈ 6 3 kpc . and ≈ 16 kpc), obtaining flat metallicity gradients in both cases. Overall, the present data on metallicity gradients in both emission-selected massive galaxies and DLA galaxies at z ≈ 2 do not suggest any inconsistency with the observed metallicity of -0 22 at an impact parameter of . ≈ 33 kpc for the z ≈ 2 193 DLA towards QSO B1228-113. .

## 6. SUMMARY

We have used ALMA in the extended C43-7 configuration to carry out the first high-resolution CO mapping of a high-redshift, H -absorption-selected galaxy, DLA-B1228g at i z ≈ 2 193. . We find that the CO(3-2) emission from DLA-B1228g is consistent with arising from a rotating disk, with a spatial extent of ≈ 6 4 kpc. . It appears unlikely that the galaxy is a major merger system, although this scenario cannot be formally ruled out by the present data. Kinematic modelling of the CO(3-2) emission as a rotating thin disk, with an exponential surface density profile, yields a rotation velocity of v rot = 328 ± 7 km s -1 and a velocity dispersion of σ v = 62 ± 7 km s -1 . The high value of the ratio v rot /σ v , 5 3 . ± 0 6, implies that DLA-B1228g is a cold, rotation-dominated disk galaxy. We obtain a dynamical mass of . ≲ 1 5 . × 10 11 M ⊙ from the kinematical modelling, within a radius of ≈ 6 kpc. This is similar to the molecular gas mass inferred from the VLA CO(1-0) emission, indicating that the inner regions of the galaxy are baryon-dominated. The non-detection of DLA-B1228g in an HST-WFC3 F105W image yields an upper limit of 1 3 M . ⊙ yr -1 on its unobscured SFR, far lower than the total SFR measured from the dust continuum emission, and indicating that the object is a highly dusty galaxy. We present a new HST-WFC3 F160W image that identifies extended stellar emission (at a rest-frame wavelength of ≈ 5000 ˚ at A z ≈ 2 193) from DLA-B1228g, close to the CO(3-2) emission, and overlapping . with the H α emission detected in an earlier VLT-SINFONI image. The DLA galaxy is likely to have a high H i mass, ≳ 7 3 . × 10 10 M , with the damped Lyman-⊙ α absorption arising from an extended H i disk. However, the velocity offset between the metal-line absorption detected in the VLT-UVES absorption spectrum and the expected disk velocity at the QSO location suggests that a simple thin disk model is not sufficient to explain the observed metal-line absorption.

We thank an anonymous referee for detailed comments and suggestions that have improved this paper. BK and NK acknowledge support from the Department of Atomic Energy, under project 12-R&amp;D-TFR-5.02-0700; NK also acknowledges support from the Science and Engineering Research Board of the Department of Science and Technology via a J. C. Bose Fellowship (JCB/2023/000030). ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The data reported in this paper are available though the ALMA archive (https://almascience.nrao.edu/alma-data/archive) with project code: ADS/JAO.ALMA #2019.1.00373.S. This research is based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. Support for Program numbers 15882 and 17483 was provided through a grant from the STScI under NASA contract NAS5-26555. This material is based upon work supported by the National Science Foundation under Grant No. 2107989.

Facilities: ALMA, HST

Software: Qubefit (Neeleman et al. 2021), Astropy (Astropy Collaboration et al. 2018), CASA (CASA Team 2022; Hunter et al. 2023).

## REFERENCES

| Akerman, C. J., Ellison, S. L., Pettini, M., & Steidel, C. C. 2005, A&A, 440, 499, doi: 10.1051/0004-6361:20052947                                        |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------|
| Alaghband-Zadeh, S., Chapman, S. C., Swinbank, A. M., et al. 2012, MNRAS, f424, 2232, doi: 10.1111/j.1365-2966.2012.21386.x                               |
| Anderson, J. 2016, Supplemental Dither Patterns for WFC3/IR, Space Telescope WFC Instrument Science Report                                                |
| Astropy Collaboration, Price-Whelan, A. M., Sip˝cz, o B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f                                     |
| Bik, A., ´ Alvarez-M´rquez, a J., Colina, L., et al. 2023, arXiv e-prints, arXiv:2312.03074, doi: 10.48550/arXiv.2312.03074                               |
| Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207, doi: 10.1146/annurev-astro-082812-140944                                                |
| Bouch´, e N., Murphy, M. T., Kacprzak, G. G., et al. 2013, Science, 341, 50, doi: 10.1126/science.1234209                                                 |
| Bouch´, e N., Murphy, M. T., P´roux, e C., et al. 2012, MNRAS, 419, 2, doi: 10.1111/j.1365-2966.2011.19500.x                                              |
| Briggs, D. S. 1995, in American Astronomical Society Meeting Abstracts, Vol. 187, American Astronomical Society Meeting Abstracts, 112.02                 |
| Broeils, A. H., & Rhee, M. H. 1997, A&A, 324, 877 Burkert, A., Genzel, R., Bouch´, e N., et al. 2010, ApJ,                                                |
| 725, 2324, doi: 10.1088/0004-637X/725/2/2324                                                                                                              |
| Calistro Rivera, G., Hodge, J. A., Smail, I., et al. 2018, ApJ, 863, 56, doi: 10.3847/1538-4357/aacffa                                                    |
| doi: 10.1088/1538-3873/ac9642 Chabrier, G. 2003, ApJL, 586, L133, doi: Chowdhury, A., Kanekar, N., & Chengalur, J. N. 2022a,                              |
| 10.1086/374879 ApJL, 935, L5, doi: 10.3847/2041-8213/ac8150                                                                                               |
| S., Dwarkhanath, K. S. 2020, Nature, 586, 369                                                                                                             |
| Christensen, L., Møller, P., Fynbo, J. P. U., & Zafar, T. 2014, MNRAS, 445, 225, doi: 10.1093/mnras/stu1726 Combes, F., & Gupta, N. 2023, arXiv e-prints, |
| Epinat, B., Contini, T., Le F`vre, e O., et al. 2009,                                                                                                     |
| F¨rster o Schreiber, N. M., Renzini, A., Mancini, C., et al. 2018, ApJS, 238, 21, doi: 10.3847/1538-4365/aadd49                                           |
| 504, 789, doi: 10.1051/0004-6361/200911995                                                                                                                |
| arXiv:2310.17204, doi:                                                                                                                                    |
| Team. 2022, PASP, 134, 114501,                                                                                                                            |
| CASA                                                                                                                                                      |
| 2022b, ApJ, 937, 103, doi:                                                                                                                                |
| -.                                                                                                                                                        |
| A., Kanekar, N., Chengalur, J. N., Sethi,                                                                                                                 |
| 10.3847/1538-4357/ac7d52 &                                                                                                                                |
| Chowdhury,                                                                                                                                                |
| 10.48550/arXiv.2310.17204                                                                                                                                 |
| A&A,                                                                                                                                                      |

| Fumagalli, M., O'Meara, J. M., Prochaska, J. X., Kanekar, N., & Wolfe, A. M. 2014, MNRAS, 444, 1282, doi: 10.1093/mnras/stu1512                       |
|-------------------------------------------------------------------------------------------------------------------------------------------------------|
| Fynbo, J. P. U., Laursen, P., Ledoux, C., et al. 2010, MNRAS, 408, 2128, doi: 10.1111/j.1365-2966.2010.17294.x                                        |
| Fynbo, J. P. U., Heintz, K. E., Neeleman, M., et al. 2018, MNRAS, 479, 2126, doi: 10.1093/mnras/sty1520                                               |
| Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1, doi: 10.1051/0004-6361/202243940                                       |
| Genzel, R., F¨rster o Schreiber, N. M., ¨ Ubler, H., et al. 2017, Nature, 543, 397, doi: 10.1038/nature21685                                          |
| Genzel, R., Price, S. H., ¨ Ubler, H., et al. 2020, ApJ, 902, 98, doi: 10.3847/1538-4357/abb0ea                                                       |
| Hack, W. J., Cara, M., Sosey, M., et al. 2020, doi: 10.5281/zenodo.3743274                                                                            |
| Hodge, J. A., Riechers, D., Decarli, R., et al. 2015, ApJL, 798, L18, doi: 10.1088/2041-8205/798/1/L18                                                |
| Hunter, T. R., Indebetouw, R., Brogan, C. L., et al. 2023, PASP, 135, 074501, doi: 10.1088/1538-3873/ace216                                           |
| Kanekar, N., Prochaska, J., Neeleman, M., et al. 2020, ApJL, 901, L5. https://arxiv.org/abs/2009.04001                                                |
| Kaur, B., Kanekar, N., Rafelski, M., et al. 2022a, ApJL, 933, L42, doi: 10.3847/2041-8213/ac7bdd -. 2021, ApJ, 921, 68, doi: 10.3847/1538-4357/ac12d2 |
| Kaur, B., Kanekar, N., Revalski, M., et al. 2022b, ApJ, 934, 87, doi: 10.3847/1538-4357/ac7b2c                                                        |
| Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531, doi: 10.1146/annurev-astro-081811-125610                                                       |
| Klitsch, A., Christensen, L., Valentino, F., et al. 2022, MNRAS, 514, 2346, doi: 10.1093/mnras/stac1190                                               |
| Krogager, J. K., Møller, P., Fynbo, J. P. U., & Noterdaeme, P. 2017, MNRAS, 469, 2959, doi: 10.1093/mnras/stx1011                                     |
| Kron, R. G. 1980, ApJS, 43, 305, doi: 10.1086/190669                                                                                                  |
| Lelli, F., Di Teodoro, E. M., Fraternali, F., et al. 2021, Science, 371, 713, doi: 10.1126/science.abc1893                                            |
| Lelli, F., Zhang, Z.-Y., Bisbas, T. G., et al. 2023, A&A, 672, A106, doi: 10.1051/0004-6361/202245105                                                 |
| Lofthouse, E. K., Fumagalli, M., Fossati, M., et al. 2023, MNRAS, 518, 305, doi: 10.1093/mnras/stac3089                                               |
| Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415, doi: 10.1146/annurev-astro-081811-125615                                                             |

| Møller, P., Fynbo, J. P. U., & Fall, S. M. 2004, A&A, 422, L33, doi: 10.1051/0004-6361:20040194                                                                                                        |
|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Neeleman, M., Kanekar, N., Prochaska, J. X., et al. ApJL, 856, L12, doi: 10.3847/2041-8213/aab5b1                                                                                                      |
| 2018, -. 2017, Science, 355, 1285, doi: 10.1126/science.aal1737                                                                                                                                        |
| Neeleman, M., Kanekar, N., Prochaska, J. X., Rafelski, M. A., & Carilli, C. L. 2019, ApJL, 870, L19, doi: 10.3847/2041-8213/aaf871                                                                     |
| Neeleman, M., Prochaska, J. X., Kanekar, N., & Rafelski, M. 2020, Nature, 581, 269, doi: 10.1038/s41586-020-2276-y                                                                                     |
| Neeleman, M., Novak, M., Venemans, B. P., et al. 2021, ApJ, 911, 141, doi: 10.3847/1538-4357/abe70f                                                                                                    |
| Oyarz´n, u G. A., Rafelski, M., Kanekar, N., et al. 2024, ApJ, 962, 72, doi: 10.3847/1538-4357/ad1182                                                                                                  |
| P´roux, e C., Bouch´, e N., Kulkarni, V. P., York, D. G., & Vladilo, G. 2011, MNRAS, 410, 2237, doi: 10.1111/j.1365-2966.2010.17598.x -. 2012, MNRAS, 419, 3060, doi: 10.1111/j.1365-2966.2011.19947.x |

| Prochaska, J. X., Neeleman, M., Kanekar, N., & Rafelski, M. 2019, ApJL, 886, L35, doi: 10.3847/2041-8213/ab55eb     |
|---------------------------------------------------------------------------------------------------------------------|
| Rizzo, F., Kohandel, M., Pallottini, A., et al. 2022, A&A, 667, A5, doi: 10.1051/0004-6361/202243582                |
| Rizzo, F., Vegetti, S., Powell, D., et al. 2020, Nature, 584, 201, doi: 10.1038/s41586-020-2572-6                   |
| Rizzo, F., Roman-Oliveira, F., Fraternali, F., et al. 2023, A&A, 679, A129, doi: 10.1051/0004-6361/202346444        |
| Tacconi, L. J., Genzel, R., & Sternberg, A. 2020, ARA&A, 58, 157, doi: 10.1146/annurev-astro-082812-141034          |
| Wang, J., Koribalski, B. S., Serra, P., et al. 2016, MNRAS, 460, 2143, doi: 10.1093/mnras/stw1099                   |
| Wang, W.-H., Kanekar, N., & Prochaska, J. X. 2015, MNRAS, 448, 2832, doi: 10.1093/mnras/stv171                      |
| Warren, S. J., Møller, P., Fall, S. M., & Jakobsen, P. 2001, MNRAS, 326, 759, doi: 10.1046/j.1365-8711.2001.04629.x |
| Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 861, doi: 10.1146/annurev.astro.42.053102.133950     |
| Wolfe, A. M., Turnshek, D. A., Smith, H. E., & Cohen, R. D. 1986, ApJS, 61, 249, doi: 10.1086/191114                |
| Wuyts, E., Wisnioski, E., Fossati, M., et al. 2016, ApJ, 827, 74, doi: 10.3847/0004-637X/827/1/74                   |
B. Kaur N. Kanekar M. Neeleman M. Rafelski J. X. Prochaska R. Dutta
An HI-absorption-selected cold rotating disk galaxy at $z\approx2.193$
We have used the Atacama Large Millimeter/submillimeter Array (ALMA) to map
CO(3-2) emission from a galaxy, DLA-B1228g, associated with the
high-metallicity damped Lyman-$\alpha$ absorber at $z \approx 2.1929$ towards
the QSO PKS B1228-113. At an angular resolution of $\approx0.32''\times0.24''$,
DLA-B1228g shows extended CO(3-2) emission with a deconvolved size of
$\approx0.78''\times0.18''$, i.e. a spatial extent of $\approx6.4$ kpc. We
detect extended stellar emission from DLA-B1228g in a Hubble Space Telescope
Wide Field Camera 3 F160W image, and find that H$\alpha$ emission is detected
in a Very Large Telescope SINFONI image from only one side of the galaxy. While
the clumpy nature of the F160W emission and the offset between the kinematic
and physical centers of the CO(3-2) emission are consistent with a merger
scenario, this appears unlikely due to the lack of strong H$\alpha$ emission,
the symmetric double-peaked CO(3-2) line profile, the high molecular gas
depletion timescale, and the similar velocity dispersions in the two halves of
the CO(3-2) image. Kinematic modelling reveals that the CO(3-2) emission is
consistent with arising from an axisymmetric rotating disk, with an exponential
profile, a rotation velocity of $v_{rot}=328\pm7$ km s$^{-1}$, and a velocity
dispersion of $\sigma_{v}=62\pm7$ km s$^{-1}$. The high value of the ratio
$v_{rot}/\sigma_v$, $\approx5.3$, implies that DLA-B1228g is a
rotation-dominated cold disk galaxy, the second case of a high-$z$
HI-absorption-selected galaxy identified with a cold rotating disk. We obtain a
dynamical mass of $M_{dyn}= (1.5\pm0.1)\times10^{11}~M_\odot$, similar to the
molecular gas mass of $\approx10^{11} M_\odot$ inferred from earlier CO(1-0)
studies; this implies that the galaxy is baryon-dominated in its inner regions.