# Molecular Gas Dominated 50 kpc Ram Pressure Stripped Tail of the Coma Galaxy D100

Pavel Jáchym, Ming Sun, J. D. P. Kenney, Luca Cortese, Françoise Combes, Masafumi Yagi, Michitoshi Yoshida, Jan Palous, Elke Roediger

► To cite this version:

Pavel Jáchym, Ming Sun, J. D. P. Kenney, Luca Cortese, Françoise Combes, et al. Molecular Gas Dominated 50 kpc Ram Pressure Stripped Tail of the Coma Galaxy D100. The Astrophysical Journal, 2017, 839, 10.3847/1538-4357/aa6af5. [hal-03730819](https://hal.science/hal-03730819v1)

HAL Id: hal-03730819

https://hal.science/hal-03730819v1

Submitted on 27 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Molecular Gas Dominated 50 kpc Ram Pressure Stripped Tail of the Coma Galaxy D100*

Pavel Jáchym<fnref n="1" />, Ming Sun<fnref n="2" />, Jeffrey D. P. Kenney<fnref n="3" />, Luca Cortese<fnref n="4" />, Françoise Combes<fnref n="5" />, Masafumi Yagi<fnref n="6" />,<fnref n="7" />, Michitoshi Yoshida<fnref n="8" />,
Jan Palouš<fnref n="4" />, and Elke Roediger<fnref n="9" />

Received 2017 February 3; revised 2017 March 20; accepted 2017 March 31; published 2017 April 24

## Abstract

We have discovered large amounts of molecular gas, as traced by CO emission, in the ram pressure stripped gas tail of the Coma cluster galaxy D100 (GMP 2910), out to large distances of about 50 kpc. D100 has a 60 kpc long, strikingly narrow tail, which is bright in X-rays and $H\alpha$. Our observations with the IRAM 30 m telescope reveal in total ~$10^9 M_{\odot} H_2$ (assuming the standard CO-to-$H_2$ conversion) in several regions along the tail, thus indicating that molecular gas may dominate its mass. Along the tail, we measure a smooth gradient in the radial velocity of the CO emission that is offset to lower values from the more diffuse $H\alpha$ gas velocities. Such a dynamic separation of phases may be due to their differential acceleration by ram pressure. D100 is likely being stripped at a high orbital velocity ~$2200$ km s$^{-1}$ by (nearly) peak ram pressure. Combined effects of intra-cluster medium (ICM) viscosity and magnetic fields may be important for the evolution of the stripped interstellar matter. We propose that D100 has reached a continuous mode of stripping of dense gas remaining in its nuclear region. D100 is the second known case of an abundant molecular stripped gas tail, suggesting that conditions in the ICM at the centers of galaxy clusters may be favorable for molecularization. From comparison with other galaxies, we find that there is a good correlation between the CO flux and the $H\alpha$ surface brightness in ram pressure stripped gas tails, over ~2 dex.

**Key words:** galaxies: clusters: individual (Coma) – galaxies: evolution – galaxies: individual (D100) – galaxies: ISM – galaxies: star formation – submillimeter: ISM

## 1. Introduction

The environments of galactic tails formed by ram pressure stripping (RPS) are likely different from typical environments of galaxy disks. The interstellar matter (ISM) stripped from galaxies infalling into clusters by dynamical pressure of the intra-cluster medium (ICM) that fills up the space in between galaxies presumably mixes with the surrounding ICM. Due to various competing thermodynamic processes a spectrum of temperatures and densities develops in the wakes of stripped galaxies. RPS as an efficient hydrodynamical mechanism of ISM removal from galaxies in clusters, as well as the fate of the stripped ISM, has been studied both observationally and theoretically (Gunn & Gott 1972; Cowie & Songaila 1977; Larson et al. 1980; Nulsen 1982; Giovanelli & Haynes 1983; Cayatte et al. 1990; Kenney & Koopmann 1999; Vollmer et al. 2001; Koopmann & Kenney 2004; Roediger & Hensler 2005; Jáchym et al. 2007; Chung et al. 2009; Kapferer et al. 2009; Tonnesen et al. 2011; Tonnesen & Bryan 2012, and others).

In the nearby, $M_{dyn} \sim 10^{14} M_{\odot}$, Virgo cluster, many short off-disk, mostly H I features are known (Kenney et al. 2004, 2014; Chung et al. 2009; Abramson et al. 2011). While there are also examples of longer ram pressure stripped gas tails (length ~ (1–2) × $D_{25}$; Oosterloo & van Gorkom 2005; Chung et al. 2007; Boselli et al. 2016), most of the gas missing in the

Virgo galaxies has not been revealed in the intra-cluster space (e.g., Vollmer & Huchtmeier 2007). In the more distant cluster A1367, with about five times the Virgo mass, several long RPS tails were observed in H I (Scott et al. 2010, 2012).

In more massive clusters ($M_{dyn} \sim 10^{15} M_{\odot}$), such as Coma or Norma, many examples of long (length ~ $3 \times D_{25}$), clearly RPS tails are observed in (1) diffuse $H\alpha$ (Gavazzi et al. 2001; Yoshida et al. 2004, 2008; Cortese et al. 2006, 2007; Sun et al. 2007; Yagi et al. 2007, 2010; Fossati et al. 2012) and (2) X-rays (Finoguenov et al. 2004; Wang et al. 2004; Machacek et al. 2005; Sun & Vikhlinin 2005; Sun et al. 2006, 2010). Some of the tails are observed in multiple wavelengths (mostly X-rays + $H\alpha$), such as ESO 137-001, NGC 4848, and GMP 2910, but also NGC 4388 in Virgo.

Rather surprisingly, regions of young star formation have been revealed in some of the tails, either in $H\alpha$ or UV (Kenney & Koopmann 1999; Cortese et al. 2006; Sun et al. 2007; Yoshida et al. 2008; Hester et al. 2010; Smith et al. 2010; Yagi et al. 2013; Ebeling et al. 2014). This indicates that while ram pressure generally suppresses star formation in the disks, it may induce new star formation in the stripped medium. Dense molecular clouds are the principal sites of star formation in galaxies. Observed ongoing star formation in some of the tails is thus suggestive of the presence of molecular gas.

Indeed, an abundant cold molecular component was discovered for the first time in the tail of the Norma cluster galaxy ESO 137-001 (Jáchym et al. 2014). The detected

<footnote n="1">¹Astronomical Institute, Czech Academy of Sciences, Boční II 1401, 14100, Prague, Czech Republic; jachym@ig.cas.cz</footnote>

<footnote n="2">²Department of Physics, University of Alabama in Huntsville, 301 Sparkman Drive, Huntsville, AL 35899, USA</footnote>

<footnote n="3">³Department of Astronomy, Yale University, 260 Whitney Avenue, New Haven, CT 06511, USA</footnote>

<footnote n="4">⁴International Centre for Radio Astronomy Research, The University of Western Australia, 35 Stirling Hwy, Crawley, WA 6009, Australia</footnote>

<footnote n="5">⁵Observatoire de Paris, LERMA, PSL, CNRS, Sorbonne Univ. UPMC, and College de France, F-75014, Paris, France</footnote>

<footnote n="6">⁶Optical and Infrared Astronomy Division, National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan</footnote>

<footnote n="7">⁷Graduate School of Science and Engineering, Hosei University, 3-7-2, Kajinocho, Koganei, Tokyo, 184-8584, Japan</footnote>

<footnote n="8">⁸Hiroshima Astrophysical Science Center, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan</footnote>

<footnote n="9">⁹Milne Centre for Astrophysics, Department of Physics & Mathematics, University of Hull, Hull, HU6 7RX, UK</footnote>

<footnote>* Based on observations made with the IRAM 30 m Telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain).</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Left: *Chandra* view of the central parts of the Coma cluster. Positions of two large elliptical galaxies NGC 4889 and NGC 4874 are shown with circles, together with the position of D100. The remarkable ~48 kpc long X-ray tail extending to the NE direction from D100 is clearly visible. Image credit: NASA/CXC/ MPE/Sanders et al. (2013). Middle: *HST* view of D100's closest neighborhood: projected to the SE is the S0 galaxy D99 (GMP2897) and to the S a weak galaxy GMP 2913. The Subaru Hα tail is overlaid in red (Yagi et al. 2010). Right: *HST* WFPC2 zoom on D100 showing prominent dust extinction filaments extending from the nucleus, as well as two spiral arms and a central bar. The RGB image was created by combining a *B* (F450W filter, blue), an *I* (F814W filter, red), and a merged *B* + *I* (green) images. A F450W image was already published by Caldwell et al. (1999).</figcaption></figure>

amounts of the cold molecular phase (~$10^9 M_\odot$) are similar to those of the hot ionized phase observed in the tail, suggesting that molecular gas may form a substantial fraction of the multiphase gas. Moreover, in the tail of ESO 137-001, the observed gas components, including the rich molecular phase, together nearly account for the missing gas in the disk. Other tails were searched for molecular content: In the Virgo cluster, only upper limits on CO emission were measured in the tail of IC 3418 (Jáchym et al. 2013), while several regions with molecular gas were recently discovered at large distances in the tail of NGC 4388 (Verdugo et al. 2015), but corresponding only to a small fraction of total mass of the gas tail. Also, in the Virgo cluster, regions of off-disk CO emission were detected in NGC 4522 (Vollmer et al. 2008), as well as in the long Hα trail that connects M86 with NGC 4438 (Dasyra et al. 2012).

Formation and survival of molecular gas in the tails of ram pressure stripped galaxies, in the surroundings of the hot ICM, is an interesting problem worth further observational and numerical efforts. Central regions of galaxy clusters have revealed that cold molecular gas may exist in the long filaments of cooling cores (e.g., Salomé et al. 2011). The central region of the Coma cluster, the most massive and most X-ray luminous cluster at $z < 0.025$, is an ideal laboratory for studying the hydrodynamic effects of the surrounding ICM on galaxies as well as the fate of the stripped gas. The Coma cluster has the richest optical data among nearby massive clusters, already with more than 20 late-type galaxies with one-sided star-forming or ionized gas tails (Smith et al. 2010; Yagi et al. 2010), several of which are also bright in soft X-rays. One of the best galaxies to study is D100, a galaxy near the Coma core with a remarkable ram pressure gas-stripped tail.

### 1.1. D100 (GMP 2910)

D100's tail is the straightest, and has the largest length to width ratio, of any known ram pressure stripped tail. It is bright in multiple wavelengths. Figure 2 shows the Subaru deep image of the galaxy (Yagi et al. 2010). Hα emission is coming from a long (~130′′ ≈ 60 kpc) and extremely narrow (mostly ~4′′5 ≈ 2.1 kpc) area that connects to the core of the galaxy (Yagi et al. 2007). The inner parts of the Hα tail are also shown

in the middle panel of Figure 1, overlaid on the *HST* image of the galaxy (Caldwell et al. 1999). The D100 tail is also bright in soft X-rays (see Figure 1, left panel). *GALEX* observations further revealed a UV component of the tail in the inner ~15 kpc of its length (Smith et al. 2010). Some star formation is thus likely taking place in the stripped material. No HI was detected with VLA in the galaxy or the tail (Bravo-Alfaro et al. 2000, 2001).

D100 is a 0.3 L\* SBab type galaxy (see Table 1) with an estimated stellar mass of $2.1 \times 10^9 M_\odot$ (Yagi et al. 2010).<fnref n="10" /> It exhibits starburst in its core, with a current star-formation rate of ~2.3 $M_\odot$ yr⁻¹ (derived from WISE band 4), and poststarburst characteristics in the rest of the disk. D100 is projected at only ~240 kpc from the Coma cluster center. Figure 1 (left panel) shows its position in a *Chandra* view of the central parts of the Coma cluster. Its radial velocity component relative to the Coma mean is approximately −1570 km s⁻¹ (Yagi et al. 2007).

In the *HST* image in Figure 1 (right panel), prominent dust features obscuring the eastern side of the disk are clearly visible (Caldwell et al. 1999). They coincide with the tail direction. The image also reveals a strong two-armed spiral pattern extending out to an ~1.4 kpc radius. While two early-type galaxies are projected close to D100, GMP 2897 (D99) at ~10 kpc and GMP 2852 at ~30 kpc, their radial velocities are substantially larger (by about 4500 and 2000 km s⁻¹, respectively). Another galaxy, a low-surface brightness GMP 2913 occurs at a projected distance of ~9 kpc from D100 (see Figure 1, middle panel). Its radial velocity is only ~130 km s⁻¹ lower than that of D100, which makes their interaction possible (Yagi et al. 2007). However, the optical isophotes of D100 are symmetric and do not indicated any (strong) recent tidal interaction.

In this paper, we report our discovery of abundant molecular gas component in the prominent tail of D100. We study its distribution and kinematics and compare it with the warm ionized component, and discuss the origin of the tail and of the molecular gas in the tail. The paper is organized as follows.

<footnote n="10">¹⁰ The stellar mass estimate comes from the MPA-JHU SDSS catalog in the DR7_v5.2 version.</footnote>
<table><caption>**Table 1**
Parameters of the Galaxy D100 (GMP 2910, PGC 044716, MRK 0060 NED 01)</caption><tr><td>R.A., Decl. (J2000)</td><td>13h00m09s14, +27°51′59″2</td></tr><tr><td>Type</td><td>SBab</td></tr><tr><td>Redshift<sup>a</sup>, V<sub>helio</sub></td><td>0.01784, 5348 km s<sup>-1</sup></td></tr><tr><td>V<sub>Coma</sub></td><td>-1570 km s<sup>-1</sup></td></tr><tr><td>Major diameter (B-band)</td><td>24″</td></tr><tr><td>Major–minor axis ratio</td><td>1.35</td></tr><tr><td>PA, inclination</td><td>178°6, ~43°</td></tr><tr><td>Total B (I) mag</td><td>16.09 ± 0.09 (15.23 ± 0.08)</td></tr><tr><td>Stellar mass<sup>b</sup></td><td>2.1 × 10<sup>9</sup> M<sub>☉</sub></td></tr></table>

**Notes.**

<sup>a</sup> From Yagi et al. (2007) and corrected from observer's to heliocentric velocity frame.

<sup>b</sup> MEDIAN stellar mass in the MPA-JHU SDSS catalog.

Our observations and results are introduced in Sections 2 and 3, kinematics of the tail is studied in Section 4, and the distribution of the CO emission and its correlation with Hα emission in Section 5. Also, the CO-to-H₂ conversion factor is discussed in Section 5. In Section 6, we discuss the origin of the tail and the efficiency of star formation in the stripped gas, and we study the orbit of the galaxy in the Coma. Conclusions are drawn in Section 7. Throughout the paper, we use the Coma distance of 97.5 Mpc, thus 1 arcsec ≈ 0.473 kpc. Cosmological constants in use were $H_0 = 70$, $\Omega_\Lambda = 0.73$, and $\Omega_m = 0.27$.

## 2. Observations

The observations were carried out with the IRAM 30 m antenna operated by the Institut de Radio Astronomie Millimé-trique (IRAM) at Pico Veleta, Spain, in 2014 December. The EMIR receiver in E090 and E230 bands was used to observe simultaneously at the frequencies of the ¹²CO(1–0) ($\nu_{rest}$ = 115.271 GHz) and the ¹²CO(2–1) ($\nu_{rest}$ = 230.538 GHz) lines. Observing conditions were excellent with PWV as low as 1–2 mm and system temperatures typically of about 140–180 K at CO(1-0) and 160–200 K at CO(2-1). The FTS spectrometer with ~200 kHz spectral resolution was connected to both lines. Also the WILMA autocorrelator with a spectral resolution of 2 MHz at both 115 and 230 GHz was used as a back-up. The observations were done in a symmetric Wobbler switching mode with the maximum throw of the secondary reflector of 240" in order to avoid, with OFF positions, the tail if oriented in azimuth.

The half power beamwidth (HPBW) of the IRAM 30 m main beam is described to a good accuracy by HPBW (") = 2460/$\nu$(GHz). Thus at the CO(1-0) and CO(2-1) sky frequencies, HPBW ~21.7" ≈ 10.3 kpc and 10.9" ≈ 5.2 kpc, respectively. The corresponding main beam projected area $\Omega_B$ ≃ 533 arcsec² = 120 kpc² and ≃135 arcsec² = 31 kpc², respectively, including a factor of 1/ln2 of a Gaussian beamshape correction.

Seven integration points were selected to cover the main body of D100, as well as most of its Hα bright tail, out to a projected distance of ~1.6 ≈ 45 kpc from the galaxy (see the scheme in Figure 2). During the observation, we first focused on Hα bright regions in the tail (pointings T1–T4) and, only after revealing strong CO emission, we moved on to intermediate parts of the tail with less (or less clumpy) Hα emission ("complementary" pointings TC2 and TC3) in the remaining

time. The list of observed positions is given in Table 2, together with information on actual on-source observing times.

The FTS backend often suffered from “platforming” between individual units. There is a correction CLASS script FtsPlatformingCorrection5.class<fnref n="11" /> that can subtract baselines individually from the affected FTS units. However, our expected lines were placed in the middle of the central FTS unit, thus platforming was not a big issue.

The data were reduced in the standard manner using CLASS from the GILDAS<fnref n="12" /> software package developed at IRAM. Bad scans were flagged and emission line-free channels in the total width of about 1000 km s⁻¹ were used to subtract (mostly) first-order baselines. The corrected antenna temperatures, $T_A^*$, provided by the IRAM 30 m calibration pipeline, were converted to main-beam brightness temperature by $T_{mb} = T_A^*F_{eff}/\eta_{mb}$, using a main beam efficiency of about $\eta_{mb} = 0.78$ at 115 GHz and 0.59 at 230 GHz, and the forward efficiencies $F_{eff}$ of 0.94 and 0.92, respectively. The rms noise levels typically of 1–1.5 mK per 10.6 km s⁻¹ channels were obtained. Gaussian fits were used to measure peak $T_{mb}$, width, and position of the detected CO lines. The flux density-to-antenna temperature conversion factor is $S_\nu/T_A^* = 3514/(\eta_a D^2)$, where $\eta_a$ is the telescope aperture efficiency and D is the diameter of the telescope in meters. The IRAM 30 m aperture efficiency is ~0.6 at CO(1-0) frequency and ~0.41 at CO(2-1) frequency. The $S_\nu/T_{mb}$ conversion is thus ~5 Jy beam⁻¹ K⁻¹ for both bands.

## 3. Results

### 3.1. Main Body of D100

The CO(1-0) and CO(2-1) spectra measured in the main body D100 pointing are shown in Figure 2 (top-right panel). Both lines are clearly detected, the CO(2-1) line is about two times stronger than CO(1-0). We calculate the CO luminosity from the standard relation of Solomon & Vanden Bout (2005)

$$L'_{CO} = 3.25 \times 10^7 S_{CO} \Delta\nu \nu_{obs}^{-2} D_L^2 (1+z)^{-3}, \quad (1)$$

where $L'_{CO}$ is the CO line luminosity in K km s⁻¹ pc², $S_{CO}$ Δ$\nu$ is the CO velocity integrated line flux in Jy km s⁻¹, $\nu_{obs}$ is the observed CO line frequency in GHz, and $D_L$ is the distance in Mpc. The CO(1-0) luminosity is ~1.1 × 10⁸ K km s⁻¹ pc². Following

$$M_{H_2} [M_{\odot}] = 4.5 L'_{CO} [\text{K km s}^{-1} \text{ pc}^2], \quad (2)$$

where we assume a CO/H₂ conversion factor of $2 \times 10^{20}$ cm⁻²(K km s⁻¹)⁻¹ that is standard under Milky Way disk conditions (e.g., Bolatto et al. 2013). However, the actual value of the conversion factor is uncertain, especially in the tail of D100. This is discussed later in Section 5.3. The luminosity corresponds to a molecular gas mass of ~4.8 × 10⁸ M☉ (including a factor of 1.36 to account for helium).

The CO lines are fitted with Gaussians in Figure 2. Their parameters are given in Table 3. The CO(1-0) central velocity is ~5248 km s⁻¹, assuming the radio definition of the velocity ($\nu_{rad}/c = 1 - \nu_{sky}/\nu_0 = z/(z+1)$). This is consistent with the optical radial velocity $v_{opt} = c(\nu_0/\nu_{sky} - 1) = c z = 5348$ km s⁻¹ measured from Hα, [N II], and [S II] emission in

<footnote n="11">¹¹ www.iram.es/IRAMES/mainWiki/AstronomerOfDutyChecklist</footnote>

<footnote n="12">¹² http://www.iram.fr/IRAMFR/GILDAS</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** CO(1-0) and CO(2-1) spectra measured in D100 (top right panel, smoothed to 10.3 km s⁻¹ channels) and the tail regions T1–T4 (smoothed to 20 km s⁻¹). The IRAM 30 m $^{12}$CO(1-0) and $^{15}$CO(2-1) beams (HPBW of 21′′7 ≈ 10.3 kpc and 10′′9 ≈ 5.2 kpc, respectively) are shown for the observed pointings overlaid on the optical Hα image of D100 (Subaru Telescope, Yagi et al. 2010). Parameters of the CO line Gaussian fits (red curves) are given in Table 3. Dashed vertical lines in the tail spectra indicate the central velocity of the CO(1-0) line in the main body. The strongest CO emission in the tail is clearly in the region T3, at about a 1′′ ≈ 32 kpc distance from the galaxy. The spectra from the two complementary regions TC2 and TC3 are shown in Figure 5. The CO velocity scale is LSR; for the sky position of D100, $V_{hel} = V_{LSR} - 8.5$ km s⁻¹.</figcaption></figure>

<table><caption>**Table 2**
List of Observed Positions</caption><thead><tr><th></th><th>R.A.<br>(J2000)</th><th>Decl.<br>(J2000)</th><th>d<sub>D100</sub><br>(kpc)</th><th>T<sub>ON</sub><br>(min)</th></tr></thead><tbody><tr><td>D100</td><td>13:00:09.14</td><td>+27:51:59.2</td><td>...</td><td>99</td></tr><tr><td>T1</td><td>13:00:09.88</td><td>+27:52:04.1</td><td>5.2</td><td>196</td></tr><tr><td>T2</td><td>13:00:11.84</td><td>+27:52:19.4</td><td>19.4</td><td>112</td></tr><tr><td>TC2</td><td>13:00:12.77</td><td>+27:52:26.2</td><td>26.1</td><td>83</td></tr><tr><td>T3</td><td>13:00:13.64</td><td>+27:52:32.4</td><td>32.3</td><td>127</td></tr><tr><td>TC3</td><td>13:00:14.53</td><td>+27:52:39.0</td><td>38.7</td><td>106</td></tr><tr><td>T4</td><td>13:00:15.45</td><td>+27:52:45.3</td><td>45.2</td><td>177</td></tr></tbody></table>

the galaxy center (Yagi et al. 2007). The CO(1-0) line width is ~120 km s⁻¹, contrary to ~80 km s⁻¹ for the CO(2-1) line. The difference may reflect detecting emission from larger galactocentric radii with increasing rotation velocity within the larger CO(1-0) beam. From the Tully–Fisher relation (e.g., Gnedin et al. 2007), a disk galaxy with D100’s stellar mass has a typical circular velocity of ~80 km s⁻¹. Given the inclination of the galaxy of ~43° (see Section 6.3), this corresponds to the linewidth of $2 \times v_{circ} \sin 43° \approx 110$ km s⁻¹, thus consistent

with the measured CO(1-0) linewidth. The larger CO(1-0) beam is also likely contaminated by the emission from the inner tail (see below).

The CO (2-1)/(1-0) integrated luminosity ratio measured in the D100 pointing is ~1.5 (and the ratio of the peak temperatures is ~2.2), not corrected for the different beam sizes. Assuming a typical CO (2-1)/(1-0) ratio of 0.8–1, the measured value is consistent with emission coming from an extended source that is larger than the CO(2-1) beam.

## 3.2. Molecular Gas in the Tail of D100

Our new CO observations of the D100 tail reveal for the first time the presence of abundant molecular gas coexisting with the previously observed Hα and X-ray components. Figure 2 depicts CO(1-0) and CO(2-1) spectra measured in the T1–T4 pointings. The complementary spectra from the two less-sensitive regions TC2 and TC3 are shown in Figure 5. The parameters of the Gaussian fits to the detected lines are summarized in Table 3. Further in this section, we divide the tail into three parts (inner, intermediate, and outer). As the morphology of the tail does not change substantially along
<table><caption>**Table 3**
Properties of the Detected CO Lines in D100</caption><thead><tr><th>Source</th><th>Line</th><th>rms (mK)</th><th>Velocity (km s<sup>-1</sup>)</th><th>Rel. Vel. (km s<sup>-1</sup>)</th><th>FWHM (km s<sup>-1</sup>)</th><th>T<sub>peak</sub> (mK)</th><th>I<sub>CO,fit</sub> (K km s<sup>-1</sup>)</th><th>L<sub>CO</sub> (10<sup>7</sup> K km s<sup>-1</sup> pc<sup>2</sup>)</th><th>M<sub>mol</sub><sup>a</sup> (10<sup>8</sup> M<sub>☉</sub>)</th></tr></thead><tbody><tr><td rowspan="2">D100</td><td>CO(1-0)</td><td>1.7</td><td>5247.9 ± 5.6</td><td>0</td><td>118.0 ± 14.2</td><td>7.4</td><td>0.93 ± 0.09</td><td>10.62 ± 1.03</td><td>4.78 ± 0.46</td></tr><tr><td>CO(2-1)</td><td>1.7</td><td>5242.7 ± 2.5</td><td>0</td><td>83.9 ± 6.5</td><td>16.</td><td>1.38 ± 0.09</td><td>3.94 ± 0.26</td><td>2.17 ± 0.14</td></tr><tr><td rowspan="4">T1</td><td>CO(1-0)</td><td>1.0</td><td>5280.7 ± 5.6</td><td>33</td><td>122.9 ± 13.7</td><td>4.5</td><td>0.59 ± 0.05</td><td>6.74 ± 0.57</td><td>3.03 ± 0.26</td></tr><tr><td>CO(2-1)</td><td>1.1</td><td>5272.4 ± 10.0</td><td>30</td><td>145.5 ± 25.5</td><td>2.7</td><td>0.43 ± 0.06</td><td>1.23 ± 0.14</td><td>0.68 ± 0.08</td></tr><tr><td>CO(2-1) #1</td><td></td><td>5193.2 ± 9.1</td><td>-55</td><td>41.6 ± 19.7</td><td>1.9</td><td>0.08 ± 0.04</td><td>0.23 ± 0.11</td><td>0.13 ± 0.06</td></tr><tr><td>CO(2-1) #2</td><td></td><td>5282.0 ± 6.8</td><td>39</td><td>82.5 ± 25.9</td><td>3.3</td><td>0.29 ± 0.06</td><td>0.83 ± 0.14</td><td>0.46 ± 0.08</td></tr><tr><td rowspan="6">T2</td><td>CO(1-0)</td><td>1.3</td><td>5328.7 ± 12.4</td><td>81</td><td>165.3 ± 27.4</td><td>3.1</td><td>0.55 ± 0.08</td><td>6.28 ± 0.91</td><td>2.82 ± 0.41</td></tr><tr><td>CO(1-0) #1</td><td></td><td>5287.4 ± 9.5</td><td>40</td><td>66.0 ± 19.5</td><td>3.4</td><td>0.24 ± 0.07</td><td>2.74 ± 0.80</td><td>1.23 ± 0.36</td></tr><tr><td>CO(1-0) #2</td><td></td><td>5379.1 ± 11.6</td><td>131</td><td>77.2 ± 25.7</td><td>3.0</td><td>0.25 ± 0.08</td><td>2.86 ± 0.91</td><td>1.29 ± 0.41</td></tr><tr><td>CO(2-1)</td><td>1.3</td><td>5350.0 ± 13.8</td><td>107</td><td>151.1 ± 30.4</td><td>2.9</td><td>0.46 ± 0.08</td><td>1.31 ± 0.23</td><td>0.72 ± 0.13</td></tr><tr><td>CO(2-1) #1</td><td></td><td>5289.4 ± 6.3</td><td>47</td><td>33.3 ± 15.6</td><td>3.1</td><td>0.11 ± 0.04</td><td>0.31 ± 0.11</td><td>0.17 ± 0.06</td></tr><tr><td>CO(2-1) #2</td><td></td><td>5374.5 ± 7.1</td><td>132</td><td>72.2 ± 17.6</td><td>3.7</td><td>0.29 ± 0.06</td><td>0.83 ± 0.14</td><td>0.46 ± 0.08</td></tr><tr><td rowspan="2">T3</td><td>CO(1-0)</td><td>1.3</td><td>5331.3 ± 3.0</td><td>83</td><td>52.8 ± 7.4</td><td>6.4</td><td>0.36 ± 0.04</td><td>4.11 ± 0.46</td><td>1.85 ± 0.21</td></tr><tr><td>CO(2-1)</td><td>1.3</td><td>5328.4 ± 1.7</td><td>86</td><td>37.7 ± 4.3</td><td>10.7</td><td>0.43 ± 0.04</td><td>1.23 ± 0.11</td><td>0.68 ± 0.06</td></tr><tr><td rowspan="2">T4</td><td>CO(1-0)</td><td>1.0</td><td>5379.5 ± 8.6</td><td>132</td><td>31.3 ± 18.5</td><td>1.4</td><td>0.05 ± 0.03</td><td>0.57 ± 0.34</td><td>0.26 ± 0.15</td></tr><tr><td>CO(2-1)</td><td>1.0</td><td>5377.5 ± 3.6</td><td>135</td><td>44.6 ± 8.1</td><td>3.7</td><td>0.18 ± 0.03</td><td>0.51 ± 0.09</td><td>0.28 ± 0.05</td></tr><tr><td rowspan="2">TC2</td><td>CO(1-0)</td><td>1.5</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>CO(2-1)</td><td>1.3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td rowspan="2">TC3</td><td>CO(1-0)</td><td>1.2</td><td>5350.7 ± 10.6</td><td>103</td><td>91.1 ± 26.4</td><td>2.8</td><td>0.27 ± 0.06</td><td>3.08 ± 0.69</td><td>1.39 ± 0.31</td></tr><tr><td>CO(2-1)</td><td>1.5</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

**Notes.** The table gives the $1\sigma$ rms in $10.2 \text{ km s}^{-1}$ channels, parameters of (multiple) Gaussian fits (the line heliocentric central velocity, the FWHM, the peak temperature, and the integrated intensity), the measured integrated intensity, and the molecular gas mass. The temperatures are given in $T_{\text{mb}}$ scale. Baselines were subtracted in the velocity range of 4200–6200 km s$^{-1}$. Radio definition of velocity was used to convert the sky frequency of the source.

<sup>a</sup> There is an additional systematic uncertainty of a factor of several due to the CO-to-H<sub>2</sub> conversion relation.

its length, we base the division only on the projected distance from the galaxy.

As noted above, the CO/H<sub>2</sub> conversion relation is uncertain in the tail of D100. It may introduce some systematic errors into the values of H<sub>2</sub> masses derived in the following subsections from the observed CO(1-0) luminosities. While it is hard to quantify the uncertainty in the conversion factor, it may be a factor of a few, but it is likely less than an order of magnitude. For more details, see the discussion in Section 5.3.

### 3.2.1. Inner Tail

In the innermost tail pointing T1, CO(1-0) emission is strong and comparable to that coming from the D100 disk itself. Following Equations (1) and (2), the corresponding H<sub>2</sub> mass is $\sim 3 \times 10^8 M_{\odot}$. CO(2-1) emission is substantially weaker than in the disk—the ratio of integrated intensities is low, $I_{\text{CO(2-1)}/I_{\text{CO(1-0)}}} \sim 0.7$. This is very likely due to contribution from the disk covered by the outer parts of the larger CO(1-0) beam. The values of the CO (2-1)/(1-0) line-ratios are for the observed regions depicted in Figure 3.

As compared to the main body pointing, the spectral lines are in T1 shifted by $\sim 40 \text{ km s}^{-1}$ to higher velocities. The CO(2-1) line profile moreover suggests that the emission is coming from two substructures at different radial velocities, with a separation of $\Delta v \approx 90 \text{ km s}^{-1}$. We fit the substructures with a double Gaussian (see Table 3). The lower-velocity peak is at $\sim 50 \text{ km s}^{-1}$ below the D100 central velocity and the higher-velocity peak at $\sim 40 \text{ km s}^{-1}$ above it (however, the latter component is much stronger than the former one). This could be an imprint of the galactic rotation.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** CO (2-1)/(1-0) integrated intensity and peak temperature line ratios measured in D100 and the tail regions T1 to T4. The values are not corrected for different beam sizes. With four times the beam area, a compact source will experience four times the beam dilution in the CO(1-0) beam.</figcaption></figure>

### 3.2.2. Intermediate Tail

At the projected distance of $\sim 20 \text{ kpc}$ from the galaxy, in the second closest region T2, both CO(1-0) and CO(2-1) emission
are again bright. The CO intensity corresponds to a H₂ mass of ~2.8 × 10⁸ M☉. The linewidths of ~150 km s⁻¹ are larger than in the main body and in the T1 region, but the line profiles clearly suggest that there are (at least) two distinct peaks with Δv ≈ 90 km s⁻¹. Presumably, the two peaks could correspond to the substructures (bifurcation) weakly visible in the Hα image. The FWHMs of the T1 and T2 sub-peaks are ~30–80 km s⁻¹, which corresponds to the velocity dispersion range of ~13–34 km s⁻¹. Much better spatial and spectral resolution are needed to further resolve smaller entities, such as individual molecular clouds whose velocity dispersions are smaller, typically ~10 km s⁻¹.

Assuming an NFW profile for the Coma mass distribution (Navarro et al. 1996; Kubo et al. 2007), and the DM halo mass of D100 of ~2 × 10¹¹ M☉ (following the relation of Behroozi et al. 2010), we estimate the tidal truncation radius $R_t \approx R(M_{D100}/2 \times M_{\text{Coma}})^{1/3}$ at the projected galactocentric distance $R = 240$ kpc of D100 in Coma to be $\lesssim 20$ kpc. Thus, the material occurring in the pointing T2 (and farther out) is beyond the tidal truncation radius and will likely contribute to the ICM. However, this does not take into account the gas velocities relative to the escape speed from the galaxy’s potential. Thus, even gas inside the tidal truncation will likely join the ICM since it will exceed the escape speed due to ram pressure acceleration.

The strongest emission (in terms of the S/N ratio) comes from the ~30 kpc distant region T3. The corresponding molecular gas mass is ~1.9 × 10⁸ M☉. The lines are substantially narrower than in the inner tail—the linewidth drops to ~40–50 km s⁻¹, with no substructures suggested. This could be due to leaving behind the circular velocity component and/or due to a compact size of the CO-emitting region, possibly associated with the bright compact Hα (probably H II) region clearly visible in the center of the T3 beam.

The CO (2-1)/(1-0) integrated intensity line ratios (not corrected for the beam sizes) are in the T2 and T3 regions again rather low, ~0.8 and ~1.2, respectively (see Figure 3). This is probably due to more CO emission coming from the outer parts of the larger CO(1-0) beam, outside the CO(2-1) beam. However, taking the distribution of Hα emission as a proxy for the CO distribution, most of the emission may be in the T3 region enclosed by the smaller CO(2-1) beam. This would then suggest a low (corrected) line ratio. In the perpendicular direction, we assume the distribution of the CO emission is confined to the width of the Hα (and X-ray) tail.

### 3.2.3. Outer Tail

In the most distant observed pointing T4, at ~45 kpc from the galaxy, CO(2-1) emission is reliably detected, while there is only a hint of CO(1-0) at the same velocity. The line intensity corresponds to a H₂ mass of ~2.6 × 10⁷ M☉. The CO(2-1) line is again narrow (FWHM ~ 40 km s⁻¹) and shifted toward higher velocities, by about 150 km s⁻¹ relative to the main body. Contrary to the previous pointings, the integrated intensity line ratio in T4 is ~3.6, thus consistent with a point-source emission diluted by the four times larger CO(1-0) beam (and a typical $R_{21}$ line ratio of ~0.8–1). Alternatively, it could suggest that molecular gas is in T4 more extended, and heated to higher temperatures, indicating optically thin gas (see, e.g., Crosthwaite & Turner 2007). Figure 3 clearly illustrates that along the tail the integrated intensity line ratio (not corrected for the different beam sizes) increases.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** CO(1-0) and CO(2-1) peak temperatures and integrated intensities in observed positions along the tail as a function of projected distance from D100.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** CO(1-0) and CO(2-1) spectra measured in the complementary regions TC2 and TC3 in the D100 tail (see the scheme in Figure 2). The rms sensitivity in CO(1-0) and CO(2-1) is 1.5 mK and 1.3 mK, respectively, for the TC2 region, and 1.2 mK and 1.5 mK, respectively, for the TC3 region. Dashed vertical lines indicate the central velocity of the CO(1-0) line in the main body spectrum.</figcaption></figure>

In Figure 4, the peak CO main beam temperatures are shown as a function of the downstream distance in the tail, together with the CO integrated intensities.

### 3.2.4. Complementary, Less-sensitive Pointings

In addition to the Hα bright regions T1–T4 that revealed strong CO emission, less-sensitive observations in regions TC2 and TC3 covering the parts of the tail with less (or less clumpy) Hα emission were performed. These regions would need more observing time to either confirm detection or decrease rms. It is still interesting to inspect their spectra in Figure 5. They suggest that CO(1-0) emission is there stronger than CO(2-1), contrary to the neighboring regions T3 and T4. In TC2, at a projected distance of ~26 kpc from the galaxy (see the scheme in Figure 2), there is a hint of CO(1-0) emission, with an S/N ratio of ~3–4, while CO(2-1) emission is absent. In TC3, at ~39 kpc, CO(1-0) emission is (marginally) detected with an S/N of ~4.6, but CO(2-1) is not detected. The corresponding H₂ mass in TC3 is ~1.4 × 10⁸ M☉.

The presumably low CO (2-1)/(1-0) line ratios in the TC2 and TC3 regions may correspond to emission coming mainly
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Velocity offset between denser, more compact CO and Hα sources and more diffuse Hα warm ionized gas along the D100 tail. Values derived from fitting Hα and [N II] spectral lines (revised from Yagi et al. 2007, converted to radio notation), and CO lines (single fits from Table 3, converted from LSR to heliocentric frame). Dashed lines show linear regression fits to both measurement sets (at distances $\ge$5 kpc). Filled circles indicate compact (possibly H II) regions in the Hα tail, as measured by Hα surface brightness ($>2.1 \times 10^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$). To estimate Hα velocity errors, 1000 Monte Carlo simulations adding different noise levels to the spectra were used. Consequently, lower S/N measurements have larger errors.</figcaption></figure>

from the outer parts of the larger CO(1-0) beams. There indeed seems to be less Hα emission, especially in the location of TC3, which could also point to less molecular emission, assuming there is a correlation between Hα and CO emission. Moreover, in both TC2 and TC3, the Hα emission is probably smoother, with no obvious compact regions, as compared to the neighboring parts of the tail (see Figure 2). The low line ratios in TC2 and TC3 could thus also suggest that conditions in the stripped gas are locally different: the gas may be more extended and have a lower temperature (since CO(1-0) usually traces a more extended component and has a lower excitation temperature than CO(2-1)).

## 4. Kinematic Separation of Dense Gas Component in the Tail

Our CO observations have revealed a radial velocity gradient of ~130 km s⁻¹ along ~50 kpc of the D100’s tail. This is visible in Figure 6, which depicts the central velocities of the Gaussian fits of the detected CO(1-0) and CO(2-1) lines as a function of their (projected) distance from the galaxy. For comparison, the plot also shows radial velocities measured from optical Hα slit spectra in 22 regions along the D100 tail (Yagi et al. 2007, with a couple of previously incorrect values revised). While the ionized gas velocity field is not smooth but shows small-scale variations, a clear trend is visible: the total Hα velocity difference between the galaxy center and outer tail is ~190 km s⁻¹. It can be described by the sum of a sharp rise in the first 4 kpc, and a more gradual rise between r = 4 kpc and r = 60 kpc of ~140 km s⁻¹. We note that the angular resolution of the two measurement sets is different—the CO(2-1) beamsize is ~5 kpc, while the optical slits are typically 0.5–1 kpc long.

Figure 6 reveals that radial velocities of the CO emission lines are systematically lower than the mean optical velocities by ~30 km s⁻¹, while having a similar slope. The CO velocities are mostly consistent with the lower edge of the velocity span of the Hα regions. We identified slit regions with large Hα surface brightnesses (possible H II regions) and marked them in the plot with filled circles. All of them are below (or in one case at) the fit.

This suggests that denser, more compact gas clumps, including CO and compact Hα sources, are less accelerated by ram pressure due to their large momentum, and lag behind more diffuse Hα-emitting gas. Such gas lumps decouple from the surrounding diffuse gas and may fall back to the galaxy potential. The same phenomenon, though observed spatially rather than kinematically, has been observed in high-resolution optical images of other cluster spirals, which show dense decoupled clouds spatially offset from more diffuse dust, which is further downstream (Crowl et al. 2005; Abramson & Kenney 2014; Kenney et al. 2015; Abramson et al. 2016). The velocity field of the stripped ISM may have a large range at any distance. This was clearly shown in numerical simulations by Tonnesen & Bryan (2010). Velocity gradients have been previously reported in ram pressure stripped gas tails of IC 4040, IC 3418, and NGC 4388 (Oosterloo & van Gorkom 2005; Yoshida et al. 2012; Kenney et al. 2014). This is the first time that effects of differential acceleration of individual phases of the ISM by ram pressure are indicated from observation.

The radial velocity gradients of ~130 km s⁻¹ in CO and nearly 200 km s⁻¹ in Hα emission form only about one-tenth of the (projected) ICM wind speed. Simulations by Tonnesen & Bryan (2012) show that even after a few hundred megayears of stripping, most of the gas does not reach the ICM wind velocity but has velocities one-third to one-half of the way to it. The large momentum of dense molecular gas prevents ram pressure to accelerate it quickly. Similarly, in the RPS Virgo cluster galaxy IC 3418, H II regions in the tail extend kinematically only to ~15% of the ICM wind velocity (Kenney et al. 2014).

Interestingly enough, the increase in velocity per unit tail length is constant along the tail (see the linear fits to both the Hα and CO measurements in Figure 6). In numerical simulations, Roediger & Brüggen (2008b) measured the velocity of the stripped gas along the tail to increase linearly with increasing distance from the galaxy (despite higher-velocity width of the tail). They found that the slope changes for different ram pressures, i.e., it is larger for stronger pressures and shallower for weaker pressures (see their Figures 12–14). Also Yoshida et al. (2012) measured almost monotonically increasing velocity of Hα regions in the tail of IC4040 with the distance from the galaxy. They measured the overall acceleration rate of ~10 km s⁻¹ kpc⁻¹ in a range of 5–80 kpc. For comparison, in the tail of D100, the average acceleration is (in projection) about 2.5 km s⁻¹ kpc⁻¹.

In order to illustrate the effects of differential ram pressure acceleration, in Figure 7, we use simple semi-analytic calculations: a set of test particles that represent ISM parcels with a range of column densities are initially distributed in a thin (2D) disk in a model galaxy potential that experiences an increasing ram pressure with a Lorentzian-shaped time profile. Both the model potential and the ram pressure pulse are set to approximately correspond to the case of D100 in Coma. Equation of motion of the particles under the combined effects of the galaxy’s potential and (face-on) ram pressure is
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Effects of differential acceleration by ram pressure as shown from semi-analytic calculations. The plot shows distances and velocities relative to the disk that can reach ISM parcels with different column densities that are stripped from within a radius of 5 kpc from a model galaxy disk. Situation after 1 Gyr, when ram pressure reaches its maximum, is depicted. Some parcels are allowed to multiply their column densities once during the calculation (by a factor of 2.5) in order to mimic condensation in situ in the tail (green triangles). Two groups of particles develop in the tail—those retaining their original density (magenta dots), and the “condensed” ones (triangles). The initial column density of the parcels is in the calculation limited to $100 M_{\odot} \text{ pc}^{-2}$. While some denser parcels could also be stripped, they would reach only small distances from the galaxy.</figcaption></figure>

calculated. Figure 7 shows the distribution of velocities and distances (relative to the galaxy) that the stripped ISM parcels reach after 1 Gyr, at the time when the ram pressure reaches its peak. The results are qualitatively consistent with the observations in Figure 6—a velocity range forms in the distribution of the stripped gas parcels as denser clumps are less accelerated. In the calculations, we moreover mimic simply in situ condensation of denser clumps by multiplying suddenly column density of randomly selected parcels. This can reproduce the offset in velocity (triangles versus dots in Figure 6)—the “condensed” parcels decouple from the main cloud of stripped parcels (they shift down in the plot). With the present calculations we only aim to illustrate the effects of differential ram pressure acceleration of the stripped gas components with different column densities in the tail, not to accurately model the tail of D100.

## 5. H₂-dominated Multiphase Gas Tail

D100 is after ESO 137-001 (Jáchym et al. 2014) only the second known example of a ram pressure stripped galaxy whose tail shows up simultaneously in X-ray, Hα, and CO emission.<fnref n="13" /> It is therefore of great interest to study how different phases co-exist and evolve in this special environment.

The hot ionized gas is expected to fill up the volume of the tail ($f \sim 1$). Its total mass is estimated to $\sim 10^8 M_{\odot}$, which corresponds to the observed (bolometric) luminosity of $\sim 1 \times 10^{40} \text{ erg s}^{-1}$ of the soft X-ray emitting gas at a

temperature of $\sim 1.0$ keV (assuming a single-component model with metallicity $0.3 Z_{\odot}$; Sanders et al. 2014). The electron density of the material is $\sim 8 \times 10^{-3} \text{ cm}^{-3}$ assuming the emission comes from a cylinder of 2 kpc radius and 36 kpc in length.

The warm ionized, Hα emitting gas is, on the other hand, expected to have a very low volume filling factor ($f \sim 0.05$; see Jáchym et al. 2013, for details). Its total mass is $\sim 2 \times 10^8 f^{-1/2} M_{\odot} \approx 4.5 \times 10^7 M_{\odot}$, assuming the emission fills a 60 kpc × 1 kpc cylinder (Yagi et al. 2007).

No H I was detected in D100 with a $3\sigma$ upper limit of $\sim 3 \times 10^8 M_{\odot}$ for the 30" beam (Bravo-Alfaro et al. 2000, 2001), which nearly corresponds to the area of the tail. Re-reduction of the original VLA data decreased the $3\sigma$ limit to $\sim 0.5 \times 10^8 M_{\odot}$ (H. Bravo-Alfaro 2017, private communication). New sensitive VLA observations are forthcoming (PI: M. Sun).

Surprisingly, the tail of D100 may be dominated by a cold molecular component, despite the hot ICM surroundings. The total H₂ mass that we detected is $\sim 1 \times 10^9 M_{\odot}$, thus by a factor of about 5–10 higher than the mass of ionized hot and warm phases together. Given the uncertainty of the X-factor (see the discussion in Section 5.3), it is possible that the molecular gas fraction is somewhat smaller, however, very likely still exceeding the fraction of the other gas components. For comparison, in the tail of the Norma galaxy ESO 137-001, the amounts of cold molecular and hot ionized gas are similar ($\sim 1 \times 10^9 M_{\odot}$; Sun et al. 2006; Jáchym et al. 2014).

In numerical simulations of ram pressure stripped gas tails (Tonnesen et al. 2011), the ratio of warm ($10^4 K < T < 10^5 K$) + hot ($7 \times 10^5 K < T < 4 \times 10^7 K$) to cold ($300 K < T < 10^4 K$) gas components varies in the range of $\sim 2–9$ for different simulation runs (after $\sim 80–110$ Myr of stripping). In the tail of D100, an analogous ratio of warm +hot ionized to cold (molecular+H I) components is thus likely much lower, $\sim 0.15$ (for the standard X-factor value). Though, more suited simulations would be needed to directly compare with the D100 observations, because the CO-emitting gas has a temperature lower than the range covered in Tonnesen et al. (2011), and the X-ray emitting gas in the D100 tail has a temperature of $\sim 1.2 \times 10^7 K$ (Sanders et al. 2014).

### 5.1. Mass Gradients along the Tail

The large-scale distribution of the stripped gas along the tail is determined by the mass-loss per orbital length, which is given by a combination of ram pressure that determines how much gas the galaxy loses in a given time, and the orbital velocity that determines over which volume the lost gas is spread (Roediger & Brüggen 2008b). In numerical simulations, the mass per orbital length decreases behind the galaxy, but at later times saturates (becomes constant) typically at about 70 kpc from the galaxy (e.g., Roediger & Brüggen 2008b). Local physical processes, mixing with the surrounding ICM, as well as separation of phases drive the evolution of the stripped gas and shape the composition of the gas in the tail.

Table 4 gives X-ray and Hα fluxes and corresponding masses encompassed by our IRAM 30 m apertures. In Figure 8, hot ionized, warm ionized, and molecular gas masses within CO(2-1) apertures T1 to T4 are plotted as a function of radial (projected) distance from the galaxy. The plot clearly suggests that the balance of the different gas phases does not change

<footnote n="13">¹³ See the Discussion section for information about the multiphase tail of the Virgo cluster galaxy NGC 4388.</footnote>
<table><caption>**Table 4**
Hα and X-Ray Fluxes and Corresponding Masses in IRAM 30 m Beams</caption><thead><tr><th>Source</th><th>Beam</th><th>F<sub>Hα</sub></th><th>M<sub>Hα</sub></th><th>F<sub>X-Ray</sub></th><th>M<sub>X-Ray</sub></th></tr></thead><tbody><tr><td rowspan="2">D100</td><td>CO10</td><td>130.5 ± 3.7</td><td>...</td><td>47.4 ± 15.0</td><td>...</td></tr><tr><td>CO21</td><td>129.1 ± 1.1</td><td>...</td><td>46.6 ± 9.1</td><td>...</td></tr><tr><td rowspan="2">T1</td><td>CO10</td><td>101.9 ± 3.7</td><td>14.86</td><td>49.1 ± 14.1</td><td>2.3</td></tr><tr><td>CO21</td><td>7.6 ± 0.3</td><td>4.05</td><td>15.8 ± 8.3</td><td>0.9</td></tr><tr><td rowspan="2">T2</td><td>CO10</td><td>11.9 ± 0.3</td><td>5.09</td><td>42.4 ± 15.8</td><td>2.1</td></tr><tr><td>CO21</td><td>6.2 ± 0.3</td><td>3.66</td><td>29.9 ± 10.0</td><td>1.3</td></tr><tr><td rowspan="2">TC2</td><td>CO10</td><td>9.4 ± 0.3</td><td>4.52</td><td>69.0 ± 11.6</td><td>2.7</td></tr><tr><td>CO21</td><td>4.6 ± 0.2</td><td>3.17</td><td>33.3 ± 8.3</td><td>1.3</td></tr><tr><td rowspan="2">T3</td><td>CO10</td><td>7.9 ± 0.3</td><td>4.14</td><td>62.4 ± 21.6</td><td>2.6</td></tr><tr><td>CO21</td><td>4.8 ± 0.2</td><td>3.22</td><td>12.5 ± 8.3</td><td>0.8</td></tr><tr><td rowspan="2">TC3</td><td>CO10</td><td>5.6 ± 0.2</td><td>3.47</td><td>...</td><td>...</td></tr><tr><td>CO21</td><td>2.0 ± 0.1</td><td>2.07</td><td>...</td><td>...</td></tr><tr><td rowspan="2">T4</td><td>CO10</td><td>5.7 ± 0.2</td><td>3.51</td><td>3.3 ± 1.7</td><td>0.6</td></tr><tr><td>CO21</td><td>4.5 ± 0.2</td><td>3.13</td><td>6.7 ± 3.3</td><td>0.6</td></tr></tbody></table>

**Note.** Hα and X-rays fluxes are in $10^{-16}$ erg s$^{-1}$ cm$^{-2}$ units, masses in $10^7 f^{1/2} M_\odot$ units, where f is the respective filling factor. CO(2-1) beam area is ~93.3 arcsec$^2$; CO(1-0) beam area is ~369.9 arcsec$^2$ (1/ln 2 factor is not included).

strongly (within a factor of two) along most of the tail, except its outer parts, where the molecular gas mass drops by a factor of ~10.

To characterize the evolution of the stripped gas, we can measure mass ratios of individual gas phases. The ratio of molecular to warm + hot ionized components decreases from about 17 (in T1) and 13 (in T2, T3) to ~2 (in T4), assuming a constant X-factor along the tail. The ratio of hot-to-warm ionized (X-ray-to-Hα) masses is quite similar in all tail regions and close to ~1–2. For comparison, in the tail of ESO 137-001, the hot-to-warm gas ratio is also nearly constant, ~2, and the ratio of cold-to-(warm+hot) masses also decreases along the tail but in a much less extreme—from ~3 in the innermost region to ~1 at 40 kpc distance downstream (Jáchym et al. 2014). The difference could correspond to different evolutionary states of stripping of D100 as compared to ESO 137-001.

A possible origin of the observed elevated molecular-to-ionized gas ratio could be heating from shocks induced by ram pressure interaction. This can create an abundant warm, CO-emitting molecular component in the tail. The contribution from shocks is also suggested from optical observations that revealed high values of the [N II]/Hα ratio of ~0.5–0.9 (Yagi et al. 2007). Also in the tail of ESO 137-001 were measured elevated values of [N II]/Hα ~ 0.4 and [O I]/Hα ≳ 0.2 (VLT/MUSE observations; Fossati et al. 2016), as well as in other Coma ram pressure stripped galaxies, IC 4040 and RB 199 (Yoshida et al. 2012).

## 5.2. Correlation of Molecular and Warm Ionized Gas Phases

Studying the correlation of individual gas phases in ram pressure stripped tails can reveal details about the evolution of the stripped ISM and its mixing with the surrounding ICM. A detailed comparison of the distributions of different gas phases will be the subject of a future study. With the limited level of spatial resolution of our current CO observations, we focus on the relationship between molecular and warm ionized gas in the tail of D100 and in tails of other known ram pressure stripped galaxies.

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Masses of individual gas components in the observed regions (IRAM 30 m CO(2-1) beams) as a function of projected distance from the galaxy. A volume filling factor of ~0.05 (see e.g., Jáchym et al. 2013) was applied to the ionized gas (Hα) masses. We note that the lines connecting individual points are intended to lead the eye and do not correspond to the actual amount of molecular gas in between the observed regions.</figcaption></figure>

In Figure 9, the integrated CO fluxes measured in the tails of D100 (present work), ESO 137-001 (Jáchym et al. 2014) and ESO 137-002 (P. Jáchym et al. 2017, in preparation) are plotted as a function of the Hα surface brightness in the respective tail areas encompassed by IRAM 30 m or APEX apertures. The tails of the two Norma galaxies are by a factor of ~3–10 brighter than that of D100, both in Hα and in CO emission. The three data points from ESO 137-001 correspond to APEX integrations at ~8, 16, and 40 kpc projected distances downstream in the tail, while those from ESO 137-002 are from regions in the outer disk and the inner tail (upper limit). There is a difference in spatial resolution as the HPBW of the IRAM 30 m beam at Coma distance is ~5.2 kpc, while the HPBW of the APEX beam at Norma distance is ~8.4 kpc.

The plot in Figure 9 reveals a rather tight correlation between cold molecular and warm ionized phase over about 2 dex, both in CO fluxes and Hα surface brightnesses. Molecular gas is the primary fuel of star formation. However, the amount of star formation in CO-abundant tails is different: there are over 30 H II regions in the tail of ESO 137-001 (Sun et al. 2007), possibly only a few H II regions in the tail of D100 (this work), but none were revealed in the tail of ESO 137-002 (Zhang et al. 2013). The observed correlation thus indicates that different excitation processes are active in the stripped gas. It is possible that stronger ram pressure (either in a richer cluster or in a galaxy that gets closer to the cluster center, or in a smaller galaxy) strips more dense gas that can form more molecular gas, as well as it produces brighter Hα emission, not only by star formation (denser, more compact component) but also due to stronger shock ionization (more diffuse component). Consequently, stronger RPS produces data points at the top-right end of the correlation.

In Figure 9, we further include points corresponding to the CO upper limits measured in three Hα bright regions in the
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Correlation between CO and Hα emission in ram pressure stripped gas tails: CO(2-1) integrated flux vs. Hα surface brightness measured in IRAM 30 m apertures (HPBW ~ 5.2 kpc) over the D100 tail and APEX apertures (HPBW ~ 8.4 kpc) over the ESO 137-001 and ESO 137-002 tails (Jáchym et al. 2014, P. Jáchym et al. 2017, in preparation). We assumed the widths of the tails of ~2 kpc (D100), ~4.8 kpc (ESO 137-001; or 5.5 kpc at the 001-C region location), and ~3 kpc (ESO 137-002). CO upper limits are marked with open symbols. CO(1-0) flux is shown for the TC3 region. Data points from the Virgo cluster galaxies IC3418 (diamonds; Jáchym et al. 2013) and NGC 4388 (squares; Verdugo et al. 2015) are also shown. The arrows on the left edge indicate CO detected regions in NGC 4388 with only weak associated Hα emission.</figcaption></figure>

outer tail of the ram pressure stripped Virgo galaxy IC3418 (Fumagalli et al. 2011; Jáchym et al. 2013), and one CO detection and two upper limits in the tail of another ram pressure stripped Virgo galaxy NGC 4388 (Verdugo et al. 2015). While these data points are basically well consistent with the revealed correlation, the spatial scale is very different—at the Virgo distance the IRAM 30 m CO(2-1) beam HPBW is ~0.9 kpc (or ~1.7 kpc at CO(1-0)), thus sampling more local conditions in the tails. Also, IC3418 is a dwarf galaxy with sub-solar metallicity, which might be the primary reason for weak CO emission, that moreover is in a late stripping state when most of molecular gas in the tail would have been consumed by star formation (Jáchym et al. 2013; Kenney et al. 2014). In NGC 4388, on the other hand, CO emission was detected in several more regions not associated with strong Hα emission ($L_{H\alpha} < 8 \times 10^{35}$ erg s⁻¹; arrows at the left edge of the plot in Figure 9).

## 5.3. Uncertainty of X-factor

The value of the CO-to-H₂ conversion factor is in the D100 tail uncertain. Local conditions in the ISM including metallicity, radiation field, temperature, density, pressure, or dust properties determine the value of the conversion between measured CO integrated intensities and H₂ column densities.

The environment of the gas-stripped tail is likely distinct from typical galactic disks in a number of parameters. For example, the volume density of the stripped gas may be lower due to ram pressure preferentially stripping lower density gas from the galaxy that moreover distributes in an extended three-

dimensional tail. Also, the relative lack of young stars and thus of strong UV emission suggests different chemistry in the tail, possibly driven by shocks and turbulence induced by the ram pressure interaction (e.g., Godard et al. 2009, and references therein).

In the inner tail of ESO 137-001, elevated ratios of warm-to-cold molecular gas of $\ge 0.1$ were observed (Sivanandam et al. 2010; Jáchym et al. 2014), indicating that some part of CO emission is coming from a warmer, more diffuse molecular component. Also, in several Virgo spirals, higher ratios of warm H₂/PAH were found (Wong et al. 2014). Recent works that analyzed the CO-to-H₂ relation in the diffuse component of the ISM in the Milky Way (Liszt et al. 2010; Liszt & Pety 2012; Roman-Duval et al. 2016) and in nearby galaxies (Sandstrom et al. 2013) surprisingly found that the value of the X-factor is in the diffuse gas similar to that in dense gas. The reason may be that on small spatial scales (arcmins in the Galaxy) there are likely large differences in the value of the X-factor in the diffuse gas between bright CO (low $X_{CO}$) and faint CO (high $X_{CO}$) regions, however, on large scales they average out (see also Bolatto et al. 2013, for comprehensive review).

While some warm, subthermally excited diffuse component producing enhanced CO emission could exist in the D100 tail, there likely are also dense, compact regions that re-formed in situ, copying the distribution of compact Hα regions (see Figure 2). For example, in the farthest region T4, the measured CO (2-1)/(1-0) line ratio is consistent with a point-source distribution.

We can also estimate the (stellar) metallicity of D100 from the stellar mass-metallicity relation of Tremonti et al. (2004). For $M_* = 2.1 \times 10^9 M_\odot$, the predicted median metallicity of D100 is rather high, $12 + \log(O/H) \sim 8.7$, i.e., about solar (we assume a solar value of 8.69; Asplund et al. 2009). While there is no data available on the gas metallicity, the above estimate suggests that the $X_{CO}$ factor is not much affected by metallicity in D100. We expect the metallicity of the stripped gas is consistent with the value in the main body (Fossati et al. 2016). Since D100 is likely in a late stripping phase (see Section 6), gas originating from outer disk regions with potentially lower metallicity due to a radial metallicity gradient currently does not contribute to the tail.

In summary, while there probably is an uncertainty in the CO/H₂ relation of a factor of a few, there is no strong evidence that the $X_{CO}$ factor is in the D100 tail systematically lower than the standard value. The detected bright CO emission then suggests that molecular gas mass dominates the tail by a factor of $\lesssim 5-10$.

## 6. Discussion

Our detection of the abundant molecular gas in the tail of D100 has revealed a substantial component of the stripped gas that until now was hidden. However, the detection has also raised new questions about the evolution of the stripped gas, regarding the origin of molecular gas in the special environment of ram pressure stripped gas tails, as well as how common the phenomena of molecular gas-rich ram pressure stripped tails are.
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Comparison of inner parts of the tails of D100 (top; Subaru; Yagi et al. 2007) and ESO 137-001 (bottom; MUSE; Fumagalli et al. 2014) as seen in Hα emission. The images are scaled so that the optical diameters of both galaxies are roughly the same (24" in D100 vs. 75" in ESO 137-001; marked with dashed ellipses). The striking difference in morphology between the two tails is clearly visible, mainly due to the presence of a broad component in the tail of ESO 137-001. The displayed lengths of the tails correspond to ~20 kpc (D100) and ~40 kpc (ESO 137-001). The dashed areas in the ESO 137-001 image, as well as the continuation of the tail to the left-hand side were not covered by the MUSE observations.</figcaption></figure>

## 6.1. Ram Pressure Stripping Origin of the Tail

The 60 kpc multi-wavelength gas tail of D100 has all the characteristics of an RPS origin because it is gaseous, one-sided, straight, and well centered on the disk of the galaxy. Moreover, the galaxy is occurring close to the Coma cluster center, where it is likely currently experiencing a peak ram pressure.

What makes it clearly stand out of other known cases is its remarkable narrowness along its length. In ESO 137-001, for comparison, a wider tail of Hα regions orphaned from X-ray emission, together with the secondary X-ray tail, accompany the main tail. This is clearly visible in Figure 10, where we show the inner parts of the Hα tails of the two galaxies. The images are scaled so that the optical diameters of the galaxies match each other. The comparison reveals that the tail of D100 has simpler morphology, with much less substructure perpendicular to its length, and it is also straighter. On the other hand, the two tails bear some similarities—they are brightest in the innermost parts and there are local enhancements along the tail with Hα peaks separated by relatively diffuse emission regions.

RPS proceeds in a galaxy from the outside in, which was confirmed by many observations of truncated H I, Hα, and/or dust disks (Koopmann & Kenney 2004). Consequently, the gas stripped at later stages is expected to form a narrower tail. Numerical simulations have confirmed that small cross-sections lead to less flaring in the tail (see, e.g., Roediger et al. 2006, their Figure 8).

The narrow profile of the D100's tail thus suggests that the galaxy is in an advanced stripping stage when its stripping radius is small, and all ISM from outer disk radii was lost to the intra-cluster space far behind the galaxy. Dense (molecular) gas is not expected to be directly stripped from the galaxy

(see Section 6.3.3). A galaxy can also have a smaller cross-section if it is moving (near) edge-on through the ICM. As we will show later in Section 6.3, this may be the case for D100, however, this effect would probably be more important in an earlier stripping phase. Radiative cooling processes were also identified to make for much narrower tails with significantly less flaring by reducing the pressure in the tail (Tonnesen & Bryan 2010). We also note that the extreme linearity of the tail may point to high viscosity of the surrounding ICM that can suppress hydrodynamic instabilities that otherwise cause vortices and turbulence (Roediger & Brüggen 2008a). Magnetic fields if aligned with the tail can further inhibit instabilities and transport processes with the surroundings making the tail smoother than it would be in the absence of magnetic fields. However, magnetic fields may lead to suppression of viscosity, thus the combined effects of the ICM viscosity and magnetic fields in the Coma center are likely important for the structure of the D100 tail.

### 6.1.1. Continuous Stripping of Nuclear Gas?

The fact that the narrow tail of D100 is not (much) widening suggests that the galaxy is heavily stripped and that the stripping radius has not changed much over the lifetime of the observed tail. At the same time, the CO luminosity of the tail is about twice that of the galaxy, so a lot of gas has been removed from the galaxy without changing the gas stripping radius much. This may be compatible with the stripping of dense gas from the (circum)nuclear region of the galaxy. D100 contains a bar (see the HST image in Figure 1, right panel) that could have driven gas from the outer disk into the inner kiloparsec via gravitational torques. This behavior has been observed in barred galaxies as indicated from their larger central molecular gas concentrations compared to unbarred galaxies (Sakamoto et al. 1999).

The stripping radius of D100 is only ~600 pc, which is about the same as the diameter of the bar. When stripping reaches the nuclear region ($r < r_{bar}$), the stripping radius does not change quickly, both due to large gas concentrations, and the gravitational potential that becomes significantly steeper there. In many galaxies, especially those with central starbursts, the gas surface density increases sharply in the nuclear region, with a scale length much smaller than that of the outer disk (Jogee et al. 2005). D100 shows an ongoing nuclear starburst (Caldwell et al. 1999; Yagi et al. 2007).

Assuming that the stripped gas (~$10^9 M_\odot$) now forming the narrow tail originated from an annular region $R = 500–1000$ pc (which corresponds to an area ~$10^6 \text{ pc}^2$), we can get a rough estimate on the surface density of the original gas in the disk of ~500 $M_\odot \text{ pc}^{-2}$. The gas surface density profiles in the circumnuclear regions are typically ~100–300 $M_\odot \text{ pc}^{-2}$ or less at $R \sim 1$ kpc, and >1000 $M_\odot \text{ pc}^{-2}$ inside $R \sim 500$ pc (Jogee et al. 2005). It is probable that an analogous galaxy to D100 with the same stellar mass but no CN gas concentration would be completely stripped at the current D100's location in the cluster, possessing only an old detached gas-stripped tail.

### 6.1.2. In Situ Origin of Molecular Gas

Formation and survival of large quantities of (cold) molecular gas in the tail is rather surprising given the presence of the surrounding hot ICM. It suggests that processes able to heat and disperse the stripped cool gas, such as heat
conduction, ionizing soft X-ray radiation, cosmic rays or turbulence, are not efficient enough to prevent the gas from cooling and condensing. For example, assuming the saturated flux equations, Tonnesen et al. (2011) estimated the efficiency of heat conduction to be low in the tails, of the order of 10%–20%, otherwise cold clouds would quickly evaporate on a timescale of ~10 Myr. The viscosity of the ICM is also important for the evolution of the stripped gas. With increasing viscosity, the stripped galactic gas mixes less readily with the ambient ICM. Viscously stripped galaxies thus are expected to have unmixed, cool wakes that are also X-ray bright (Roediger et al. 2015).

A key precondition for H₂ formation is the presence of dust in the tail, it means its stripping and survival. The HST image of D100 (see Figure 1, right panel) clearly shows strong dust extinction filaments extending from the disk central regions in the direction of the tail. This indicates that dust has been stripped to the tail; however, its survival to large distances from the galaxy is yet to be explored. The D100 tail was also covered by Herschel observations, but only upper limits on FIR emission in the tail can be obtained (will be done elsewhere; S. Sivanandam et al. 2017, in preparation). HST images of the Virgo galaxies NGC 4522 and NGC 4402 that are being stripped clearly show dust being removed from the disk along with the gas (Abramson et al. 2016).

Given the large amounts of detected molecular gas along the tail, including large distances from the galaxy, together with the small measured velocity gradients (measured in Section 4), in situ molecularization of the stripped material appears to be a probable scenario of the origin of the detected gas. As we will estimate in Section 6.3.3, ram pressure operating on D100 could have directly stripped gas components as dense as ≈50 M☉ pc⁻² down to about a 1 kpc disk radius. However, due to strong internal density gradients in real GMCs reaching values of hundreds of M☉ pc⁻², gradual ablation of dense clouds due to ram pressure is expected rather than pushing them as a whole. The molecular component of the ISM is also expected to be dissociated during stripping by shock introduced by ram pressure interaction. Moreover, the lifetime of dense molecular clouds is typically 10⁶–10⁷ Myr, which likely is by at least a factor of 10 shorter than the presumed tail age.

Numerical hydrodynamical simulations suggested that overdensities may form from the less dense stripped gas that has cooled and condensed in the tail (e.g., Tonnesen & Bryan 2010, 2012). The timescale for condensation and H₂ formation is determined by the gas density following an inverse relation (Guillard et al. 2009). The densest gas lumps in the tail were likely formed from the densest gas lumps that were stripped. Rather dense clumps stripped from the circumnuclear gas concentration of D100, that further were compressed by the surrounding ICM, started to cool down efficiently (the cooling time in the tail is expected to be rather short). Lack of ambient UV photo-dissociating radiation field in the tail further can favor H₂ formation on dust grains.

## 6.2. Star Formation in Extreme Environments

It is of great interest to study star formation in ram pressure stripped gas tails since these are distinct environments from typical star-forming regions in disks of spiral galaxies.

### 6.2.1. Weak Star Formation in the Tail

The detailed Subaru Telescope Hα image of the tail of D100 shows a great deal of substructure in the ionized gas distribution (despite its narrowness) with several compact regions visible at different distances from the galaxy (Yagi et al. 2010, see the image in Figure 2). The three most compact Hα sources also have the lowest optical [N II]/Hα ratios in the D100 tail (in the range of 0.44–0.54; Yagi et al. 2007, their Table 2). However, these values are too large for typical H II regions (generally <0.3; e.g., Fossati et al. 2016, their Figure 7). Thus, they do not look like typical H II regions formed around star-forming clouds.

We can obtain a rough estimate on the star-formation rate (SFR) in the tail if we assumed the most compact sources were H II regions. From Hα fluxes in the three most compact regions identified from the Hα Subaru image in Yagi et al. (2007).<fnref n="14" /> The corresponding fluxes are ~ (3.4, 13.0, and 8.8) × 10⁻¹⁷ erg s⁻¹ cm⁻². Using the SFR–L<sub>Hα</sub> formula of Kennicutt & Evans (2012)

$$ \log \text{SFR} (M_{\odot} \text{ yr}^{-1}) = \log L_{\text{H}\alpha} (\text{erg s}^{-1}) - 41.27, \quad (3) $$

the corresponding summed SFR is ~3.9 × 10⁻³ M☉ yr⁻¹ assuming a typical 1 mag extinction correction for Hα. For comparison, in the observed regions in the tail of ESO 137-001, the summed SFR is by a factor of 10 stronger, ~4 × 10⁻² M☉ yr⁻¹ (Jáchym et al. 2014).

A weak UV emission is also seen in the inner ~15 kpc of the tail with GALEX (Smith et al. 2010). The presence of both Hα and UV emission could suggest star formation has been ongoing in the tail within ~10 Myr (which is the timescales of the Hα emission as a tracer of star formation). However, much UV emission seems to come from background galaxies that are around the tail and the small residual is hard to be isolated for GALEX's large PSF. The current GALEX data thus cannot be used to constrain the SFR robustly. Smith et al. (2010) have also shown a MegaCam archival u-band deep image with higher resolution than in the GALEX image. The u-band continuum shows the tail clearly up to ~30 kpc. While the u-band versus Hα overlay does show a lot of correspondence between these two tracers, most of the u-band emission may not be from stars.

We thus conclude that some star formation may be present in D100's tail, though it is weak. Given the wealth of discovered molecular gas, the efficiency of star formation must be very low (characterized by molecular gas depletion times τ<sub>dep,H₂</sub> = M<sub>H₂</sub>/SFR, that are larger than the Hubble time). Previous observations of ESO 137-001 and NGC 4388 indeed suggested that most of molecular gas in the stripped gas tails does not form stars and ultimately joins the ICM (Jáchym et al. 2014; Verdugo et al. 2015).

### 6.2.2. Starburst in the Disk

Our proposed scenario of RPS of nuclear gas concentration is consistent with the starburst observed in the center of D100: the nucleus of the galaxy (within ~2′′) exhibits starburst characteristics, such as strong emission and strong underlying Balmer absorption. Larger radii in the disk (~3′′) show a poststarburst feature (only Balmer absorption) with a quenching time between

<footnote n="14">¹⁴ The sources #8, #17, and #19, located at (13:00:10.53, +27:52:10.88) at the edge of the CO(1-0) T1 region, (13:00:13.64, +27:52:33.44) in the T3 region, and (13:00:15.4, +27:52:45.6)—in the T4 region, respectively.</footnote>
<figure><img src="image_13.png" /></figure>

**Figure 11.** Scheme of the D100 tail geometry: the disk inclination angle is $i \sim 43^{\circ}$; the tail extends away from us; the E side of the disk is the far side. Consequently, the plane-of-sky orbital velocity component is equal to or larger than the l-o-s orbital velocity component. Thus, the wind angle (w) is close to edge-on (for $v_{sky} \approx v_{los}$) or $\lesssim 43^{\circ}$ (for $v_{sky} \gg v_{los}$), assuming that the tail follows the past orbital path of the galaxy.

0 and 250 Myr, depending on the burst strength (Caldwell et al. 1999; Yagi et al. 2007).

We can check the depletion timescale of the presumed nuclear region in D100. In the main body pointing, we detected $\sim 5 \times 10^8 M_{\odot}$ of molecular gas. The estimated star-formation rate in the nucleus is $\sim 2.3 M_{\odot} \text{ yr}^{-1}$ (derived from WISE band 4), thus $\tau_{\text{dep}} \approx 0.22 \text{ Gyr}$, which is short and indeed consistent with a burst.

## 6.3. Galaxy's Orbit in Coma

To better understand the timescale of stripping of D100, we analyze the available parameters of its orbit. The galaxy is projected close to ($8.5 \approx 240 \text{ kpc}$) the cluster center and its tail points nearly perpendicularly to the direction to the cluster center. This suggests that D100 is likely near the pericenter of its orbit in Coma.<fnref n="15" /> However, the details of the geometry of the ram pressure interaction are not immediately clear.

### 6.3.1. Constraining the Wind Angle

We can get an interesting constraint on the ram pressure interaction geometry and the 3D velocity of the interaction from the HST image of the galaxy (Figure 1, right panel). Apparent in the image are elongated filaments of dust extinction, which can be seen against the Eastern disk side. The filaments are elongated in the tail direction, and are likely part of the base of the tail. We also know that the stripped gas tail must extend in the direction away from us, since the radial velocity of D100 is approximately $-1580 \text{ km s}^{-1}$ relative to the Coma mean ($\sim 6925 \text{ km s}^{-1}$; Yagi et al. 2007). Since the kinematics indicate that the tail must extend away from us, and the extinction features indicate that the base of the tail is in front of the Eastern side of the stellar disk, the Eastern side of the disk must be the far side. The geometry is shown in Figure 11 in the plane-of-the-sky-line-of-sight view.

To constrain the wind direction relative to the disk plane (wind angle), we first derive the disk inclination angle by inspecting isophotes from the HST image. While slight lopsidedness of the disk is apparent, consistent values of the major-to-minor axes ratio are obtained at several different radii. The corresponding average value of the disk inclination is $\sim (43 \pm 3)^{\circ}$. Since the tail is extending away from us, and at the same time it is in front of the E disk side, the tail angle relative to the disk plane must be in the range of $\sim 0^{\circ}-43^{\circ}$. This is illustrated in Figure 11. If we assume that the tail follows the past orbital path, the wind angle is equal to the tail angle. Then the above analysis prescribes that the plane-of-sky (tangential) orbital velocity component is equal to or larger than the radial velocity component of D100. This means that the total orbital velocity of D100 is $\gtrsim 1580 \sqrt{2} \approx 2200 \text{ km s}^{-1}$.

The system orientation implies that the larger the total orbital velocity (and thus the tangential component), the larger the wind angle, but it is always smaller than $\sim 43^{\circ}$. However, for reasonable values of the velocity, the interaction is close to edge-on: for the total orbital velocity of $3000 \text{ km s}^{-1}$ (or $4000 \text{ km s}^{-1}$), the wind angle would be only $\sim 11^{\circ}$ (or $\sim 21^{\circ}$).

Numerical simulations have suggested that (close-to) edge-on ram pressure stripping is less efficient than face-on or only slightly inclined stripping (e.g., Roediger et al. 2006; Jáchym et al. 2009). D100, however, has likely been experiencing very strong (near pericenter) ram pressure in which case stripping saturates and the stripping efficiency difference between edge-on and face-on stripping disappears (see Jáchym et al. 2009).

### 6.3.2. Estimated Age of the D100 Tail

Knowing the disk inclination angle and constraining the tail angle relative to the disk allows us to deproject the tail length ($\sim 60 \text{ kpc in } H\alpha$), as well as the velocity gradient along the tail ($\sim 150 \text{ km s}^{-1}$ in $H\alpha$). For the 3D orbital velocity of $3000 \text{ km s}^{-1}$ (or $4000 \text{ km s}^{-1}$), the deprojected tail length would be $\sim 71 \text{ kpc}$ (or $\sim 65 \text{ kpc}$), and the deprojected velocity gradient along the tail would be $\sim 283 \text{ km s}^{-1}$ (or $\sim 400 \text{ km s}^{-1}$).

We can do a zeroth-order estimate on the age of the tail: to cross the length of 71 kpc (or 65 kpc) at the velocity of $283 \text{ km s}^{-1}$ (or $400 \text{ km s}^{-1}$), the duration of $\sim 245 \text{ Myr}$ (or $\sim 159 \text{ Myr}$) would be needed. These are probably lower limits because we have neglected the fact that gas parcels are from the galaxy ram pressure accelerated only gradually. On the other hand, ram pressure changes steeply with time as the galaxy orbits through the cluster.

From the simple estimate of the tail age, we can also calculate a zeroth-order estimate on the average mass-loss rate due to stripping. The total gas mass of $\sim 1.2 \times 10^9 M_{\odot}$ detected in the tail would be stripped from the galaxy over 245 Myr (or 159 Myr) at a high rate of $5-8 M_{\odot} \text{ yr}^{-1}$.

It is interesting to realize that within the time the tail of the observed length has formed, the galaxy has traveled a much larger distance through the cluster, possibly exceeding $\sim 0.7-0.8 \text{ Mpc}$ (neglecting the fact that the orbital velocity peaks at the pericenter). Thus the environment in which the tail had started to form is likely completely disconnected from the galaxy's current location.

### 6.3.3. Current Ram Pressure Estimate

Modeling the ICM distribution in the Coma cluster with a $\beta$-profile with parameters given by Mohr et al. (1999) and Fossati et al. (2012), we can get an upper limit estimate on the

<footnote n="15">¹⁵ The location and orientation of D100 in Coma is similar to that of the completely stripped dwarf IC3418 in Virgo (Jáchym et al. 2013; Kenney et al. 2014).</footnote>
ICM density at the projected location of D100¹⁶ of ~3.3 × 10⁻²⁷ g cm⁻³. Assuming the orbital velocity of 3000 km s⁻¹ (or 4000 km s⁻¹), the current ram pressure is ~3.0 × 10⁻¹⁰ dyne cm⁻² (or ~5.3 × 10⁻¹⁰ dyne cm⁻²). The value is lower by a factor of several if the (deprojected) orbit has larger 1-o-s distance from the mid-point of the cluster.

For comparison, the current ram pressure acting on D100 may be about two times (or 3.5 times) higher than the estimated current ram pressure on ESO 137-001 in the Norma cluster (Jáchym et al. 2014). Moreover, the effects of ram pressure likely have been much stronger on D100 due to its approximately five times lower stellar mass. Solving semi-analytically the equation of motion of gas parcels in a D100-like model galaxy (see Jáchym et al. 2013, 2014, assuming a ~2 × 10¹¹ $M_\odot$ dark matter halo and a radial orbit), it is probable that the gas parcels with column densities up to ~50 $M_\odot$ pc⁻² could be completely (directly) stripped from the galaxy, down to ~1 kpc radius.

## 6.4. Galaxy Transformation Watched Live

Several studies of S0 galaxies have suggested that their bulges contain younger stellar populations that are also richer in metallicity than the surrounding disks (e.g., Cortesi et al. 2013; Johnston et al. 2014). This indicates that the last episode of star formation occurs in the central regions of galaxies that are being transformed by cluster environments. The starburst observed in the circumnuclear region of D100, that otherwise is gas-truncated, is thus consistent with the image of the galaxy occurring at the late stage of its transformation toward passive S0-type (Johnson et al. 2016). We are lucky enough to witness the ongoing transformation, moreover accompanied by the spectacular multiphase gas tail.

## 7. Conclusions

With the IRAM 30 m telescope we have discovered large amounts of molecular gas in the ram pressure stripped gas tail of the Coma cluster galaxy D100 that up to now was known to be bright in $H\alpha$ and X-rays. After ESO 137-001 (Jáchym et al. 2014), this is only the second known example of a molecular gas-rich RPS tail. While it is currently not clear how common such tails are, it is possible that RPS galaxies may be an important source of cold gas (and sometimes stars) to intra-cluster space. The main results of our analyses are as follows.

1. Bright CO emission was detected in several regions along the tail of D100 out to a nearly 50 kpc distance from the galaxy, with the corresponding total $H_2$ mass ~$10^9 M_\odot$. While the value of the X-factor in the tail is uncertain and may be locally lower by a factor of a few than the standard Galactic value, our observations indicate that molecular gas is the dominant mass component of the D100 tail. In situ formation of molecular gas in the stripped gas is a preferable scenario of its origin.

2. The extremely narrow morphology of the tail that does not change with length indicates that D100 is currently at a late evolutionary stage when strong ram pressure has been continuously stripping the nuclear dense gas. Along the tail, at a spatial scale of ~5 kpc, the distribution of the stripped gas, including the molecular, warm and hot

ionized components, is nearly flat, indicating that the balance between different gas components does not change much. The elevated ratio of molecular-to-ionized gas mass measured in the D100 tail may be due to excitation and heating from shocks induced by ram pressure interaction.

3. The kinematics of the stripped gas in the D100 tail reveals an offset of ~30 km s⁻¹ between denser, more compact gas, including CO and compact $H\alpha$ sources, and more diffuse $H\alpha$-emitting gas. Such a dynamic separation may be due to differential acceleration of the ISM gas phases by ram pressure. The stripped gas is along the tail accelerated to only a fraction (~10%) of the presumed ram pressure wind speed that is likely above 3000 km s⁻¹. The age of the (visible) tail is ~200 Myr. Over that period, the galaxy has traveled ~800 kpc through the cluster and thus the formation of the tail has started in a completely disconnected environment.

4. Some star formation is likely present in the D100 tail, but it is weak. Given the detected large amounts of molecular gas, the efficiency of star formation is very low. This may be due to a presumably mostly diffuse morphology of the molecular gas. While there are several compact regions visible in the $H\alpha$ image of the tail, their optical line ratios may be too large for typical H II regions.

5. The new CO data for D100 in comparison with a limited sample of other galaxies with known CO components in ram pressure stripped gas tails indicate a rather tight correlation between the CO integrated intensity and the $H\alpha$ surface brightness that holds over ~2 dex.

Forthcoming interferometric CO observations of D100 tail will reveal the distribution of molecular gas in the tail as well as more details about the kinematics of the cold component and allow for a better comparison with the kinematics of other gas phases. The present results emphasize the importance of multi-wavelength observations for better understanding of the fate and evolution of ram pressure stripped gas.

This work has been supported by the project LG 14013 of the Ministry of Education, Youth and Sports of the Czech republic, the grant project 15-06012S of the Czech Science Foundation, and the institutional research project RVO:67985815. M.S. thanks the support provided by the National Aeronautics and Space Administration through Chandra Award Number GO6-17111X issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. We thank Hector Bravo-Alfaro for re-reducing existing VLA data. Based in part on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We further acknowledge the usage of the HyperLeda database (http://leda.univ-lyon1.fr).

Facility: IRAM 30 m.

## References

<footnote n="16">¹⁶ For fully ionized gas with primordial abundances (mass of helium is 0.25 of the total mass of the gas), electron and the proton number densities relates to the mass density of the ICM as $n_e = 1.167 n_p$ and $\rho = 1.143 n_e$ amu (e.g., Pavlovski et al. 2008).</footnote>

<footnote>Abramson, A., Kenney, J., Crowl, H., & Tal, T. 2016, AJ, 152, 32</footnote>

<footnote>Abramson, A., & Kenney, J. D. P. 2014, AJ, 147, 63</footnote>

<footnote>Abramson, A., Kenney, J. D. P., Crowl, H. H., et al. 2011, AJ, 141, 164</footnote>
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

Behroozi, P. S., Conroy, C., & Wechsler, R. H. 2010, ApJ, 717, 379

Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207

Boselli, A., Cuillandre, J. C., Fossati, M., et al. 2016, A&A, 587, A68

Bravo-Alfaro, H., Cayatte, V., van Gorkom, J. H., & Balkowski, C. 2000, AJ, 119, 580

Bravo-Alfaro, H., Cayatte, V., van Gorkom, J. H., & Balkowski, C. 2001, A&A, 379, 347

Caldwell, N., Rose, J. A., & Dendy, K. 1999, AJ, 117, 140

Cayatte, V., van Gorkom, J. H., Balkowski, C., & Kotanyi, C. 1990, AJ, 100, 604

Chung, A., van Gorkom, J. H., Kenney, J. D. P., Crowl, H., & Vollmer, B. 2009, AJ, 138, 1741

Chung, A., van Gorkom, J. H., Kenney, J. D. P., & Vollmer, B. 2007, ApJL, 659, L115

Cortese, L., Gavazzi, G., Boselli, A., et al. 2006, A&A, 453, 847

Cortese, L., Marcillac, D., Richard, J., et al. 2007, MNRAS, 376, 157

Cortesi, A., Merrifield, M. R., Coccato, L., et al. 2013, MNRAS, 432, 1010

Cowie, L. L., & Songaila, A. 1977, Natur, 266, 501

Crosthwaite, L. P., & Turner, J. L. 2007, AJ, 134, 1827

Crowl, H. H., Kenney, J. D. P., van Gorkom, J. H., & Vollmer, B. 2005, AJ, 130, 65

Dasyra, K. M., Combes, F., Salomé, P., & Braine, J. 2012, A&A, 540, A112

Ebeling, H., Stephenson, L. N., & Edge, A. C. 2014, ApJL, 781, L40

Finoguenov, A., Briel, U. G., Henry, J. P., et al. 2004, A&A, 419, 47

Fossati, M., Fumagalli, M., Boselli, A., et al. 2016, MNRAS, 455, 2028

Fossati, M., Gavazzi, G., Boselli, A., & Fumagalli, M. 2012, A&A, 544, A128

Fumagalli, M., Fossati, M., Hau, G. K. T., et al. 2014, MNRAS, 445, 4335

Fumagalli, M., Gavazzi, G., Scaramella, R., & Franzetti, P. 2011, A&A, 528, A46

Gavazzi, G., Boselli, A., Mayer, L., et al. 2001, ApJL, 563, L23

Giovanelli, R., & Haynes, M. P. 1983, AJ, 88, 881

Gnedin, O. Y., Weinberg, D. H., Pizagno, J., Prada, F., & Rix, H.-W. 2007, ApJ, 671, 1115

Godard, P., Falgarone, E., & Pineau Des Forêts, G. 2009, A&A, 495, 847

Guillard, P., Boulanger, F., Pineau Des Forêts, G., & Appleton, P. N. 2009, A&A, 502, 515

Gunn, J. E., & Gott, J. R., III 1972, ApJ, 176, 1

Hester, J. A., Seibert, M.,'Neill, J. D., et al. 2010, ApJL, 716, L14

Jáchym, P., Combes, F., Cortese, L., Sun, M., & Kenney, J. D. P. 2014, ApJ, 792, 11

Jáchym, P., Kenney, J. D. P., Růžička, A., et al. 2013, A&A, 556, A99

Jáchym, P., Köppen, J., Palouš, J., & Combes, F. 2009, A&A, 500, 693

Jáchym, P., Palouš, J., Köppen, J., & Combes, F. 2007, A&A, 472, 5

Jogee, S., Scoville, N., & Kenney, J. D. P. 2005, ApJ, 630, 837

Johnson, H. L., Harrison, C. M., Swinbank, A. M., et al. 2016, MNRAS, 460, 1059

Johnston, E. J., Aragón-Salamanca, A., & Merrifield, M. R. 2014, MNRAS, 441, 333

Kapferer, W., Sluka, C., Schindler, S., Ferrari, C., & Ziegler, B. 2009, A&A, 499, 87

Kenney, J. D. P., Abramson, A., & Bravo-Alfaro, H. 2015, AJ, 150, 59

Kenney, J. D. P., Geha, M., Jáchym, P., et al. 2014, ApJ, 780, 119

Kenney, J. D. P., & Koopmann, R. A. 1999, AJ, 117, 181

Kenney, J. D. P., van Gorkom, J. H., & Vollmer, B. 2004, AJ, 127, 3361

Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531

Koopmann, R. A., & Kenney, J. D. P. 2004, ApJ, 613, 866

Kubo, J. M., Stebbins, A., Annis, J., et al. 2007, ApJ, 671, 1466

Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692

Liszt, H. S., & Pety, J. 2012, A&A, 541, A58

Liszt, H. S., Pety, J., & Lucas, R. 2010, A&A, 518, A45

Machacek, M., Dosaj, A., Forman, W., et al. 2005, ApJ, 621, 663

Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517, 627

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563

Nulsen, P. E. J. 1982, MNRAS, 198, 1007

Oosterloo, T., & van Gorkom, J. 2005, A&A, 437, L19

Pavlovski, G., Kaiser, C. R., Pope, E. C. D., & Fangohr, H. 2008, MNRAS, 384, 1377

Roediger, E., & Brüggen, M. 2008a, MNRAS, 388, L89

Roediger, E., & Brüggen, M. 2008b, MNRAS, 388, 465

Roediger, E., Brüggen, M., & Hoeft, M. 2006, MNRAS, 371, 609

Roediger, E., & Hensler, G. 2005, A&A, 433, 875

Roediger, E., Kraft, R. P., Nulsen, P. E. J., et al. 2015, ApJ, 806, 104

Roman-Duval, J., Heyer, M., Brunt, C. M., et al. 2016, ApJ, 818, 144

Sakamoto, K., Scoville, N. Z., Yun, M. S., et al. 1999, ApJ, 514, 68

Salomé, P., Combes, F., Revaz, Y., et al. 2011, A&A, 531, A85

Sanders, J. S., Fabian, A. C., Churazov, E., et al. 2013, Sci, 341, 1365

Sanders, J. S., Fabian, A. C., Sun, M., et al. 2014, MNRAS, 439, 1182

Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5

Scott, T. C., Bravo-Alfaro, H., Brinks, E., et al. 2010, MNRAS, 403, 1175

Scott, T. C., Cortese, L., Brinks, E., et al. 2012, MNRAS, 419, L19

Sivanandam, S., Rieke, M. J., & Rieke, G. H. 2010, ApJ, 717, 147

Smith, R. J., Lucey, J. R., Hammer, D., et al. 2010, MNRAS, 408, 1417

Solomon, P. M., & Vanden Bout, P. A. 2005, ARA&A, 43, 677

Sun, M., Donahue, M., Roediger, E., et al. 2010, ApJ, 708, 946

Sun, M., Donahue, M., & Voit, G. M. 2007, ApJ, 671, 190

Sun, M., Jones, C., Forman, W., et al. 2006, ApJL, 637, L81

Sun, M., & Vikhlinin, A. 2005, ApJ, 621, 718

Tonnesen, S., & Bryan, G. L. 2010, ApJ, 709, 1203

Tonnesen, S., & Bryan, G. L. 2012, MNRAS, 422, 1609

Tonnesen, S., Bryan, G. L., & Chen, R. 2011, ApJ, 731, 98

Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898

Verdugo, C., Combes, F., Dasyra, K., Salomé, P., & Braine, J. 2015, A&A, 582, A6

Vollmer, B., Braine, J., Pappalardo, C., & Hily-Blant, P. 2008, A&A, 491, 455

Vollmer, B., Cayatte, V., Balkowski, C., & Duschl, W. J. 2001, ApJ, 561, 708

Vollmer, B., & Huchtmeier, W. 2007, A&A, 462, 93

Wang, Q. D., Owen, F., & Ledlow, M. 2004, ApJ, 611, 821

Wong, O. I., Kenney, J. D. P., Murphy, E. J., & Helou, G. 2014, ApJ, 783, 109

Yagi, M., Gu, L., Fujita, Y., et al. 2013, ApJ, 778, 91

Yagi, M., Komiyama, Y., Yoshida, M., et al. 2007, ApJ, 660, 1209

Yagi, M., Yoshida, M., Komiyama, Y., et al. 2010, AJ, 140, 1814

Yoshida, M., Ohyama, Y., Iye, M., et al. 2004, AJ, 127, 90

Yoshida, M., Yagi, M., Komiyama, Y., et al. 2008, ApJ, 688, 918

Yoshida, M., Yagi, M., Komiyama, Y., et al. 2012, ApJ, 749, 43

Zhang, B., Sun, M., Ji, L., et al. 2013, ApJ, 777, 122