MNRAS 000 ,

1-9 (2018)

## A common origin for the circumnuclear disc and the nearby molecular clouds in the Galactic Centre

Alessandro Ballone 1 2 , , Michela Mapelli 1 2 3 4 , , , , Alessandro Alberto Trani 5

- 1 INAF - Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, I-35122 Padova, Italy
- 2 INFN - Padova, Via Marzolo 8, I-35131 Padova, Italy
- 3 Physics and Astronomy Department Galileo Galilei, University of Padova, Vicolo dell'Osservatorio 3, I-35122, Padova, Italy
- 4 Institute f¨r Astro- und Teilchen Physik, Universit¨t Innsbruck, Technikerstrasse 25/8, A-6020 Innsbruck, Austria u a
- 5 Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan

Accepted XXX. Received YYY; in original form ZZZ

## ABSTRACT

The origin of the molecular clouds orbiting SgrA* is one of the most debated questions about our Galactic Centre. Here, we present the hydrodynamic simulation of a molecular cloud infalling towards SgrA*, performed with the adaptive-mesh-refinement code RAMSES. Through such simulation, we propose that the circumnuclear disc and the +20 km/s cloud originated from the same tidal disruption episode, occurred less than 1 Myr ago. We also show that recent star formation is to be expected in the +20 km/s cloud, as also suggested by recent observations.

Key words: black hole physics - Galaxy: centre - ISM: clouds

## 1 INTRODUCTION

The inner few parsecs of the Galactic Centre are an extremely dense environment, hosting a supermassive black hole (SMBH) with mass of around 4 × 10 6 M glyph[circledot] (Boehle et al. 2016; Gillessen et al. 2017), often referred to as SgrA* (i.e., its radio point source counterpart), a large number of late-type stars, few hundreds young and massive stars with very peculiar properties and a widely variegated interstellar medium (for general reviews, see Morris &amp; Serabyn 1996; Genzel et al. 2010; Mapelli &amp; Gualandris 2016).

The late-type stars are in an almost isotropic (though some triaxality is observed in the inner 2 pc; FeldmeierKrause et al. 2017) dense nuclear star cluster, with mass of around 2 × 10 7 M glyph[circledot] in the inner 10 pc (Sch¨del et al. 2007, o 2009, 2014; Feldmeier et al. 2014; Chatzopoulos et al. 2015; Fritz et al. 2016; Feldmeier-Krause et al. 2017; Gallego-Cano et al. 2018).

These young stars are surrounded by a circumnuclear disc (CND) of molecular and neutral gas and dust, sitting at 1.5-4 pc from SgrA* (e.g., Gatley et al. 1986; Guesten et al. 1987; Jackson et al. 1993; Herrnstein &amp; Ho 2002; Shukla et al. 2004; Christopher et al. 2005; Montero-Casta˜o et al. n 2009; Oka et al. 2011; Mart´ ın et al. 2012; Lau et al. 2013; Mills et al. 2013). The CND also contains several tens of M glyph[circledot] of ionized gas in clumpy filaments/arms observable in radio and infrared continuum and emission lines (i.e., the so-called 'minispiral'; see, e.g., Lo &amp; Claussen 1983; Serabyn &amp; Lacy 1985; Paumard et al. 2004; Zhao et al. 2009, 2010; Lau et al. 2013; Tsuboi et al. 2016; Moser et al. 2017).

On the other hand, most young stars reside in the inner pc and a considerable percentage (20-50%) of these lie in one disc (Bartko et al. 2009; Yelda et al. 2014). Most young stars are classified as O-type (many even in a Wolf-Rayet phase) and have an age of few Myr (4-8 Myr: Paumard et al. 2006; 2.5-6 Myr: Lu et al. 2013), hence suggesting that their birth occurred in a recent, single star formation episode. In addition to the disc, a second group of few tens of Btype stars (with age glyph[lessorsimilar] 20 Myr; Habibi et al. 2017), the socalled S-stars, reside in the inner 0.05 pc of the Galaxy and have isotropically oriented, highly eccentric orbits around the SMBH (Ghez et al. 2008; Gillessen et al. 2009, 2017).

Within the inner 20 pc in projected distance from SgrA*, two large molecular clouds, with mass glyph[greaterorsimilar] 10 5 M glyph[circledot] are also observed, often called +50 km/s (or M-0.02-0.07) and +20 km/s (or M-0.13-0.08) (G¨sten &amp; Downes 1980; u Guesten &amp; Henkel 1983; Liszt et al. 1985; Sandqvist 1989; Lis &amp; Carlstrom 1994; Novak et al. 2000). In particular, the +50 km/s lies at a distance of roughly 7 pc from SgrA*, it has a roughly spherical shape and a radius of around 5 pc, though a large portion of the cloud seems to be compressed by the radio non-thermal shell SgrA East (interpreted as a supernova remnant; e.g., Ho et al. 1985; Genzel et al. 1990; Serabyn et al. 1992; McGary et al. 2001; Lee et al. 2003, 2008). The +20 km/s cloud lies around 12 pc away from SgrA* and it is elongated ( ≈ 15 pc × 7.5 pc) in the direction of the SMBH. Recent observations also suggest the presence of few protostars embedded in the +20 km/s cloud (Lu et al. 2015, 2017). According to some authors, these two

clouds orbit around the Galactic Centre at distances of 50100 pc from the centre of the Galaxy, possibly in a larger belt composed of many molecular clouds (G¨sten &amp; Downes u 1980; G¨sten et al. 1981; Zylka et al. 1990; Kruijssen et al. u 2015; Henshaw et al. 2016; Kruijssen et al. 2019). On the other hand, many authors agree on rather placing them at distances glyph[lessorsimilar] 20 pc, based on their possible interaction with SgrA East (e.g., Genzel et al. 1990; McGary et al. 2001; Herrnstein &amp; Ho 2005; Lee et al. 2008) and on their internal velocity pattern and the presence of several streamers of gas connecting them to the CND (e.g., Okumura et al. 1989; Coil &amp; Ho 1999, 2000; Ferri`re 2012; Oka et al. 2012; Liu e et al. 2012; Hsieh et al. 2017; Takekawa et al. 2017; Tsuboi et al. 2018). In particular, the so-called southern streamer seems to continously extend from the +20 km/s cloud to the CND, feeding it and possibly the inner pc of the Galaxy (in particular, Okumura et al. 1991; Coil &amp; Ho 1999, 2000; Montero-Casta˜o et al. 2009). n

A connection between the inner few parsecs of the galaxy and the close-by molecular clouds is also supported by pure theory. Molecular clouds 'engulfing' the SMBH are often invoked to explain the formation of the disc of young stars (Sanders 1998; Wardle &amp; Yusef-Zadeh 2008; Bonnell &amp; Rice 2008; Hobbs &amp; Nayakshin 2009; Alig et al. 2011; Mapelli et al. 2012; Lucas et al. 2013; Alig et al. 2013) and/or the origin of the CND (Sanders 1998; Wardle &amp; Yusef-Zadeh 2008; Mapelli &amp; Trani 2016; Trani et al. 2018; however, other authors modeled the CND as a long-lived stable structure, 'carved' by feedback from the young stars or by tidal shear, see, e.g., Vollmer &amp; Duschl 2001, 2002; Blank et al. 2016).

Here we draw on this ideas to show that the CND and the +20 km/s cloud might be explained by the same tidal distruption episode.

## 2 METHODS

The simulation we used for this study is part of a set of runs of different molecular clouds in the Galactic Centre, with different masses, radii and initial orbital conditions that were performed to investigate the impact of these on star formation in the inner (tens of) parsecs of our Galaxy (Ballone et al., in preparation). Among this set we found this particular run that seems to qualitatively match some of the properties of the observed CND and +20 km/s cloud. We must then stress that the parameteres of this simulation were not fine-tuned and not all the observed features will be exactly reproduced.

## 2.1 Numerical setup

The simulation has been run with the adpative mesh refinement (AMR) code RAMSES (Teyssier 2002). The domain is cubic, with a size of 100 pc 3 , and it is centred on the position of the SMBH. The resolution of the grid spans 9 levels of refinement, from level 6, corresponding to a resolution of ∆ x 6 = 100 2 / 6 ≈ 1 56 pc, . to level 14, corresponding to a resolution of ∆ x 14 = 100 2 / 14 ≈ 6 1 . × 10 -3 pc. In order to get a good refinement of the dense (and possibly collapsing) material, without getting too high computational cost, we used a so-called 'quasi-Lagrangian' scheme (Teyssier 2015, e.g.,), that forces any cell to be refined if the enclosed mass

Figure 1. Plot of the orbital orientation in the computational domain. The green curve marks the orbit of the initial center of the cloud, from the beginning of the simulation (green square) up to 0.8 Myr (reached at the tip of the arrow). The orbital plane is marked in light blue, with the blue solid line marking the direction opposite to the initial apocentre (we remind the reader that the orbit is precessing). The orange plane is the plane of the sky in our projection (see section 3), the red solid line marks the direction of the Galactic latitude (longitude equal to 0) axis and the red dotted line marks the direction of the line-of-sight. The purple solid line marks the direction of the ascending node and the purple cross marks the position of the SMBH (and the center of the computational domain).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-04716d36-44e9-42a9-b6ea-ad063a629ef1.jpg)

is higher than a certain value m l . For level l = [6 14], we , chose m l = k m l ref , with m ref ≈ 5 82 . × 10 -2 M glyph[circledot] and k l = { 1 1 1 10 , , , 4 , 10 2 , 10 3 1 , , } . Such choice of k l has been based on ensuring that the molecular cloud is refined, at the beginning, at least to level 8; at the same time, it guarantees that only the densest parts of the cloud, after turbulence develops (see Section 2.2), are followed at the highest resolution 1 . We adopted a Lax-Friedrichs Riemann solver, which is the most diffusive, but also the only stable one for our problem, and a MinMod slope limiter for the piecewise linear reconstruction of the Godunov scheme.

## 2.2 Physical setup

The gravity of the SMBH has been simulated by introducing a particle of mass M SMBH = 4 3 . × 10 6 M glyph[circledot] (this is roughly the estimate by Gillessen et al. 2009, 2017) from the beginning of the simulation. Such particle is a 'sink' (hence it can accrete matter, see later) and has been anchored to the centre ( x = y = z = 0) of the simulation domain, to

1 m ref is roughly 1/64 of the jeans mass of the smallest cells, at level 14.

Figure 2. Time evolution of the cloud plunge. The different panels show a density projection along the y axis of the computational domain, at t = 0 15 (panel a), 0.20 (panel b), 0.35 (panel c) and 0.55 (panel d) Myr. The red star marks the position of the SMBH. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-df90d273-657b-47eb-94c2-63c9a61d5a05.jpg)

avoid artificial kicks due to pure numerical errors. The gravitational contribution of the stellar cusp of stars and the nuclear stellar disc (i.e., the so-called 'nuclear bulge'; Launhardt et al. 2002) has also been included as a spherically symmetric rigid potential. The adopted density profile for both components is ρ NB ( r ) = 2 8 . × 10 ( 6 r/r NB ) -γ M glyph[circledot] pc -3 (where r is the distance from the centre of the computational domain), with r NB = 0 22 pc, . γ = 1 2 for . r &lt; r NB and γ = 1 75 for . r &gt; r NB (e.g., Sch¨ odel et al. 2007). For completeness, we also included the Galactic bulge potential, as a Hernquist sphere with density profile ρ GB ( r ) = 2 9 . × 10 10 r GB / [2 πr r ( + r GB ) 3 ] M glyph[circledot] , where r GB = 0 7 kpc . (Mapelli &amp; Trani 2016).

The simulated molecular cloud has mass M cl = 10 5 M glyph[circledot] and radius R cl = 10 pc (corresponding to a density ρ cl ≈ 1 6 . × 10 -21 g cm -3 ) and it has been initially placed at a distance r cl = 30 pc ( x cl = y cl = z cl = 30 / √ 3 pc) from the SMBH, with tangential velocity v t,cl = 10 km s -1 ( v x,cl = v y,cl = 10 / √ 6 km s -1 , v z,cl = -20 / √ 6 km s -1 ). This means that the cloud is initially on one of the diagonals of our computational cubic domain and that the orbit lies on the x=y plane (see Fig. 1). The initial pericenter distance of the cloud center of mass is at a distance of ≈ 0 76 . pc from SgrA*. We gave the cloud a temperature of T cl = 100 K. The cloud is initially given a random Gaussian, divergence-free turbulent velocity field, with power spectrum P k ( ) = || δ 2 v || ∝ k -4 . This is the so-called Burgers power spectrum (Burgers 1948) and it is the one that seems to reproduce, in molecular clouds, the observed trend of the velocity dispersion with the cloud size and the size of its subregions (e.g., Larson 1981). Such turbulent field is needed to produce density peaks in the cloud that later collapse by self-gravity and lead to star formation all over the cloud. At the beginning of the simulation, the cloud is marginally bound (i.e., its turbulent kinetic energy is equal to its potential energy). This also means that the initial turbulent kinetic energy is a factor of 50 larger than the thermal energy of the cloud.

Though the Galactic Centre and its molecular gas are certainly embedded in a complex gaseous environment, for the sake of simplicity we just included in the simulation the molecular gas of the initial cloud, surrounding it with a low-density 'empty atmosphere', with density ρ atm = 1 6 . × 10 -27 g cm -3 and temperature T atm = 10 6 K.

To model the thermodynamics of the simulation we decided to adopt a simple isothermal equation of state (EoS).

This assumption is not capturing the complex heatingcooling interplay affecting the temperature of molecular gas; however, most of the clouds in the Galactic Centre show a temperature of ≈ 50 -150 K (Ao et al. 2013), which is probably the temperature of thermal equilibrium (e.g., Larson 1985; Koyama &amp; Inutsuka 2000) in this region.

Star formation is treated in RAMSES through the sink particle algorithm by Bleuler &amp; Teyssier (2014). This algorithm is quite sophisticated. In extreme synthesis, it checks for cells in the computational domain with density higher than ρ sink , groups these cells into clumps and performs for these a 'collapse', a 'proximity' and a 'virial' 2 check to ensure that the clump will collapse onto a point-like object. For our simulation we adopted ρ sink = 3 × 10 -17 g cm -3 , which is approximately the density for which the Jeans length equals 4 times the minimum adopted cell size (∆ x 14 , see Section 2.1). We performed few tests and found that the value of ρ sink does not significantly impact the result on star formation. The sinks forming during the simulation are all born with an initial 'seed' mass of 3 5 . M glyph[circledot] , which is about the jeans mass at the highest resolution. During the evolution, sink particles can still accrete gas, according to the Bondi-Hoyle accretion algorithm implemented in RAMSES (again, for details, check Bleuler &amp; Teyssier 2014). To avoid spurious accretion due to lack of resolution, we set an upper limit to the sink mass of 150 M glyph[circledot] : if a sink reaches this upper limit during the simulation, no further accretion onto that sink is allowed. This is true for all sinks, except the SMBH, which, as already mentioned, also behaves as a sink in the simulation, hence it is gaining mass with time.

## 3 RESULTS

In Fig. 2, we show the time evolution of the simulation. The cloud gets immediately stretched by the tidal force of the background gravitational potential of the SMBH and the stars, reaching the maximal squeezing at the closest distance from the SMBH (see panel a). The cloud reaches its first apocentre after roughly 0.4 Myr, though the parts of the cloud with lowest angular momentum already start falling back towards the SMBH and whirl around it (panel b and c). At this stage, the apocentre position has significantly changed, compared to the initial one, due to very fast orbital precession induced by the massive stellar cusp (Gualandris et al. 2012; Trani et al. 2018). A circumnuclear ring is already clearly visible after around 0.5 Myr from the beginning of the simulation, while most of the cloud is still elongated towards the SMBH and has not fallen back yet. After this time, the bulk of the cloud has reached its second pericentre passage and its leading part is already receding again, though a fraction of this post-pericentre material is streaming back towards the SMBH (panel d).

In order to try to reproduce the current configuration of the Galactic Centre, we had to define a new reference system corresponding to the actual line of sight (with versor

2 The 'virial' check, as performed in RAMSES' sink particle algorithm, is particularly useful, since it also takes into account the tidal force from the background potential, which is particularly important for our problem.

Figure 3. Density projection along the line-of-sight in the new reference system (see text in section 3.1) for our scenario, at t = 0.550 Myr. The axes represent the Galactic longitude and latitude. The upper panel shows all the gas in the computational domain, the lower panel is just a zoom-in around the inner 16 pc of the domain. The red star labels the position of the SMBH.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-a58c386e-cf82-40a5-b6de-b4c50fda3428.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-10ab1f5e-e713-482d-a6d9-3b20ad029cc5.jpg)

ˆ), Galactic longitude (with versor ˆ ) and latitude (with vers l sor ˆ ). b In this way, we could define our orbit through three angles defining its orientation with respect to the new 'observational' reference system. Specifically, i was defined as the inclination between the orbital plane and the plane of the sky, Ω was defined as the angle, in the sky plane, between the direction of the ascending node and the Galactic latitude (l=0) axis and ω is the angle, in the orbital plane, between the direction of the ascending node and the direction opposite to the initial apocentre (the pericentre direction is changing with time, due to precession; see above). To convert physical lengths to longitude and latitude, we assumed a Galactic Centre distance of 8.32 kpc (Gillessen et al. 2017).

In the following, we will present our main scenario, which is able to explain the formation of the CND and its interaction with the +20 km/s cloud, as well as to reproduce some observational properties of these two molecular structures.

.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-b832389e-4c14-4912-9c3d-ce265da5e8d8.jpg)

0 0

Figure 4. Longitude-velocity diagram for our scenario, at t = 0.550 Myr. The colormap shows the mass in every pixel of the diagram, scaled to the maximum value. At this time our simulation is able to reproduce some of the features of the observed CND and +20 km/s cloud.

## 3.1 Geometry

To simultaneously reproduce the +20 km/s cloud and the CND we had to choose ( x, y, z ) ˆ s ≈ -( 0 050 0 132 . , . , -0 990), . ( x, y, z ) ˆ l ≈ ( -0 661 . , -0 747 . , -0 007) . and ( x, y, z ) ˆ b ≈ (0 749 . , -0 651 . , -0 125), . hence with ˆ b = ˆ s × ˆ . l This corresponds to three orbital orientation angles (see above) i ≈ 97 4 . ◦ , Ω ≈ 267 2 . ◦ and ω ≈ 35 3 . ◦ . Fig. 1 shows the orbital plane and the sky plane, as well as few reference directions, in the adopted projection.

This projection led to the density maps shown in Fig. 3. In this case, the molecular cloud is 'wrapping' around the CND from the south-east as it seems to be suggested by the observations (cf. with the southern streamer/ridge; e.g., Coil &amp; Ho 2000; Ferri`re 2012; Liu et al. 2012; Hsieh et al. 2017). e The zoomed map on the bottom panel of Fig. 3 also shows several smaller streamers in the north-west part of the CND, that quite well resemble the so-called western streamers (W1, W-2, W-3 and W-4; see Fig. 2 in Liu et al. 2012 and Fig. 4 in Hsieh et al. 2017). The inclination of the ring with respect to the sky plane is also similar to the one of the observed CND. On the other hand, the molecular cloud is bigger than the +20 km/s. This problem might potentially be solved by slightly different initial condtions. For example, in this scenario, the +20 km/s cloud can (and does) have a maximum elongation equal to its orbital apocentre distance. Hence, reducing the apocentre distance to the actual length of the +20 km/s cloud (i.e., about 15 pc) would make the simulated cloud less elongated. Furthermore, an initially more compact cloud would be stretched less by the tidal field, also reducing the extension of the simulated cloud.

## 3.2 Longitude-velocity diagram

The longitude-velocity diagram in Fig. 4 is produced including all the gas in our simulation. This plot clearly shows the ring in the inner 2 pc as an oblique line going from positive

1

Figure 5. Density projection, at t = 0.5 Myr, along the line-ofsight for the scenario described in section ?? . This plot shows the position of the sink particles formed in the simulation (red dots). The axes represent the Galactic latitude and longitude, though in physical units (pc).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-6e7a5e3f-0b4f-4bf5-899e-a986cc56e447.jpg)

line-of-sight velocity (around +100 km/s) and positive latitude to negative line-of-sight velocity (around -150 km/s) and negative latitude. This is roughly comparable to what is observed for the CND (cf. with Fig. 1 in Oka et al. 2011, Fig. 8 in Hsieh et al. 2017 and Fig. 1 in Takekawa et al. 2017) and it is coherent with the orbital motion of the cloud (i.e., the CND is corotating with the +20 km/s cloud). The infalling molecular cloud is occupying the right part of the PV diagram, with a median line-of-sight velocity of roughly +20 km/s, though with a moderately larger spread in velocity, compared to the observations. The longitude-velocity diagram also shows the presence of gas with higher velocity, compared to the CND. This is gas that is orbiting the SMBH inside (i.e., faster than) the ring. We know that gas inside the CNDis present but at higher temperatures (i.e., the minispiral, see section 1), most probably ionized by the disc of young massive stars. Due to our simplified assumption of isothermal gas and the lack of such ionization sources, this material is appearing in our simulated longitude-velocity diagram, but it would not if we were including such more sophisticated physical ingredients and produce different longitude-velocity diagrams for different gas phases. We do not detect, in this diagram, any feature similar to the negative-longitude extension (NLE), as visible in the longitude-velocity diagram in Fig. 1 of Takekawa et al. (2017).

## 3.3 Star formation in the +20 km/s cloud

In Fig. 5 we show the position, at t = 0.5 Myr, of all the 32 sink particles formed in the simulation (see section 2.2) at this time. These sink particles should trace (proto)stars, born by the collapse of the gas in local density peaks seeded by the turbulence. All the sinks shown in this plot are younger than 0.3 Myr, i.e., they formed during or after the first pericentre passage of the cloud. This is compatible with the observational indications of ongoing star formation in the +20 km/s cloud, presented by Lu et al. (2017). In par-

Figure 6. Density projection along the line-of-sight, at t = 0.675 Myr, in a different reference system compared to our reference scenario (see Appendix A). The axes represent the Galactic longitude and latitude.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-da5e9888-009f-4adf-b1fb-190f8ef6ce89.jpg)

25

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-6b964374-1e69-4af1-926e-38c5b16021fe.jpg)

1

Figure 7. Density projection along the Galactic latitude in the new reference system (see text in sections 3.1) for our scenario. The axes represent the Galactic longitude (in degrees) and the line-of-sight distance from SgrA* (in pc). The red star labels the position of the SMBH.

ticular, these authors found 13 H 2 O masers coincident with submm continuum detected clumps and further 5 gravitationally bound clumps detected only in the submm continuum, for a total of 18 potential protostars. Such number is roughly of the order (within a factor of 2) of the number of sinks formed in the simulation 3

3 A possible caveat is that our resolution does not allow to count single stars, i.e., our sinks might also represent very small (proto)clusters of stars. Furthermore, we lack stellar feedback, which might also play a role in this sense. Nonetheless, the sink location should capture the regions where collapse has occurred, which is more or less what has been done by Lu et al. (2017), looking at the gravitationally bound clumps.

## 4 DISCUSSION

We have proposed a scenario to explain the simultaneous formation of the CND and the +20 km/s cloud via the tidal disruption of a large molecular cloud. In this scenario, the +20 km/s cloud reaches the Galactic centre from the southeast of the CND (as already suggested by several other authors; e.g., Okumura et al. 1991; Coil &amp; Ho 1999).

In particular, our'picture'is somewhat similar to one of the models presented in Vollmer &amp; Duschl (2002). However, in our simulation the CND is directly formed by the plunge of the cloud, while in Vollmer &amp; Duschl (2002) it is assumed to be formed at earlier times.

First of all, as a general remark, we must stress again that the initial conditions of our simulation were not finetuned, hence we do not expect every observational feature to be perfectly matched. For example, the total mass of the +20 km/s cloud and the CND is estimated to be around 4 × 10 5 M glyph[circledot] (Ferri`re 2012; however, this values are controe versial; see, e.g., Tsuboi et al. 2018), while our simulated cloud has an initial mass of only 10 5 M glyph[circledot] . Nonetheless, our scenario depends mostly on the orbit of the cloud, hence we do not expect a significant change of our results under the assumption of an initially higher cloud mass.

We must also acknowledge that, compared to other models, in this scenario the +20 km/s cloud reaches SgrA* with a much closer pericentre. In fact, most of the currently available models connect the +20 km/s and the +50 km/s cloud on the south-east of SgrA*, through the so-called molecular ridge. This hypothesis is indeed supported by the presence of the molecular ridge and by velocity maps showing a smooth velocity gradient between the two clouds (for the most recent velocity maps, see, e.g., Liu et al. 2012; Hsieh et al. 2017). On the other hand, the presence of the southern streamer (e.g., Okumura et al. 1991; Coil &amp; Ho 1999) and of several smaller streamers (such as W-1, W-2, W-3 and W-4 shown in Fig. 2 of Liu et al. 2012 and Fig. 4 of Hsieh et al. 2017) suggest a much 'closer' interaction between the +20 km/s cloud and the CND, which is naturally provided in our scenario.

A connection of the CND and the +20 km/s to the +50 km/s in the framework of a single cloud plunge is possible. In fact, as shown in Fig. 6, when passing pericenter the cloud might appear as two different clouds, connected by a streamer of gas. Fig. 6 might visually remind of the actual projected positions of the CND and the two observed molecular clouds. However, reconciling this configuration with the observations is more problematic (though we try to provide a more through comparison and a discussion of pros and cons of this further possible scenario in Appendix A).

A further constraint from observations could come from the estimated line-of-sight distance from SgrA*. Unfortunately, there is no agreement on the exact position of the +20 km/s cloud along the line-of-sight. For example, Vollmer et al. (2003) put the +20 km/s in front of SgrA*, at a maximumline-of-sight distance of [-50,-25] pc. On the other hand, Ferri`re (2012), based on a collection of studies from differe ent authors, prefers to place the +20 km/s cloud still in front of SgrA*, but at a distance of its center of [-12,-4] pc. In Fig. 7, we show a density map of our scenario, projected along the Galactic latitude. As visible, the +20 km/s cloud lies in

front of SgrA*, at distances that are compatible with the estimates of Ferri`re (2012). e

Finally, we must acknowledge that our scenario does not reproduce all the observed features (e.g., it does not reproduce the NLE found by Takekawa et al. 2017). This is possibly due to some simplified treatment of the involved physics. Nonetheless our model is quite succesful in simultaneously explaining many of them. Furthermore, our conclusions could also be affected by the supernova remnant SgrA East. This supernova remnant has an estimated age of glyph[greaterorsimilar] 10 4 yr, which is much shorter than the orbital time of our cloud. As already shown by, e.g., Coil &amp; Ho (2000) and Herrnstein &amp; Ho (2005), some of the observed spatial and velocity features, close to the CND, might also be explained by the interaction between SgrA East and the observed molecular complex. The inclusion of such supernova explosion in our simulation might be part of a follow-up study.

## 5 SUMMARY

Our simulation of a molecular cloud infalling towards SgrA* matches many features of the observed CND and +20 km/s cloud. In particular, it can

- · explain the formation of the CND;
- · link the CND to the + 20 km/s cloud, as suggested, e.g., by Takekawa et al. (2017);
- · show that star formation in the +20 km/s cloud is possible, as indicated by the recent work by Lu et al. (2015).

Specifically, the + 20 km/s cloud and the CND might both be the result of a single tidal disruption event of a large molecular cloud. The present work is to be thought as a quite effective proof of concept to explain such scenario.

More sophisticated simulations should include additional complex physical ingredients, such as a more realistic equation of state, a complex chemical network to study different molecular species separately and their thermodynamical impact on the gas and the radiative and mechanical feedback from stars. The latter ingredient, in the form of supernova explosion, would be, for example, crucial to try to account for the effect of SgrA East on the molecular gas complex. Nonetheless, the amount of observed features explained by the present model is quite remarkable, considering these limitations and that our initial conditions were not fine-tuned.

## ACKNOWLEDGEMENTS

AB and MM acknowledge financial support from the MERAC Foundation, through grant 'The physics of gas and protoplanetary discs in the Galactic centre'. AAT acknowledges support from the JSPS KAKENHI Grant Number 17F17764. We would like to thank Prof. Romain Teyssier for his help with RAMSES and the whole ForDyS group for useful discussions. Most of the simulation post-processing was carried out with the yt toolkit (Turk et al. 2011).

## REFERENCES

| Alig C., Burkert A., Johansson P. H., Schartmann M., 2011, MN- RAS, 412, 469                                                                                                                                                                                                                                                                                                                |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Alig C., Schartmann M., Burkert A., Dolag K., 2013, ApJ, 771, 119                                                                                                                                                                                                                                                                                                                           |
| Ao Y., et al., 2013, A&A, 550, A135 Bartko H., et al., 2009, ApJ, 697, 1741 Blank M., Morris M. R., Frank A., Carroll-Nellenback J. J., Duschl W. J., 2016, MNRAS, 459, 1721 Bleuler A., Teyssier R., 2014, MNRAS, 445, 4015 Boehle A., et al., 2016, ApJ, 830, 17                                                                                                                          |
| Physics, 82, 3121 Ghez A. M., et al., 2008, ApJ, 689, 1044 Gillessen S., Eisenhauer F., Trippe S., Alexander T., Genzel R., Martins F., Ott T., 2009, ApJ, 692, 1075                                                                                                                                                                                                                        |
| Bonnell I. A., Rice W. K. M., 2008, Science, 321, 1060 Burgers J. M., 1948, Advances in Applied Mechanics, 1, 171 Chatzopoulos S., Fritz T. K., Gerhard O., Gillessen S., Wegg C., Genzel R., Pfuhl O., 2015, MNRAS, 447, 948 Christopher M. H., Scoville N. Z., Stolovy S. R., Yun M. S., ApJ, 622, 346                                                                                    |
| 2005, Coil A. L., Ho P. T. P., 1999, ApJ, 513, 752 Coil A. L., Ho P. T. P., 2000, ApJ, 533, 245 Feldmeier-Krause A., Zhu L., Neumayer N., van de Ven G., de Zeeuw P. T., Sch¨del o R., 2017, MNRAS, 466, 4040 Feldmeier A., et al., 2014, A&A, 570, A2 Ferri`re e K., 2012, A&A, 540, A50 Fritz T. K., et al., 2016, ApJ, 821, 44 Gallego-Cano E., Sch¨del o R., Dong H., Nogueras-Lara F., |
| Gallego-Calvente A. T., Amaro-Seoane P., Baumgardt H., 2018, A&A, 609, A26 Gatley I., Jones T. J., Hyland A. R., Wade R., Geballe T. Krisciunas K., 1986, MNRAS, 222, 299 Genzel R., Stacey G. J., Harris A. I., Townes C. H., Geis N., U. U., Poglitsch A., Stutzki J., 1990, ApJ, 356, 160 Genzel R., Eisenhauer F., Gillessen S., 2010, Reviews of                                       |
| R., Graf Modern                                                                                                                                                                                                                                                                                                                                                                             |
| Gillessen S., et al., 2017, ApJ, 837, 30 Gualandris A., Mapelli M., Perets H. B., 2012, MNRAS, 427, 1793 Guesten R., Henkel C., 1983, A&A, 125, 136                                                                                                                                                                                                                                         |
| Guesten R., Genzel R., Wright M. C. H., Jaffe D. T., Stutzki J., Harris A. I., 1987, ApJ, 318, 124 G¨sten u R., Downes D., 1980, A&A, 87, 6 G¨sten u R., Walmsley C. M., Pauls T., 1981, A&A, 103, 197 Habibi M., et al., 2017, ApJ, 847, 120 Henshaw J. D., et al., 2016, MNRAS, 457, 2675                                                                                                 |
| Herrnstein R. M., Ho P. T. P., 2002, ApJ, 579, L83 Herrnstein R. M., Ho P. T. P., 2005, ApJ, 620, 287 Ho P. T. P., Jackson J. M., Barrett A. H., Armstrong J. T., ApJ, 288, 575 Hobbs A., Nayakshin S., 2009, MNRAS, 394, 191                                                                                                                                                               |
| 1985, Hsieh P.-Y., Koch P. M., Ho P. T. P., Kim W.-T., Tang Y.-W., Wang H.-H., Yen H.-W., Hwang C.-Y., 2017, ApJ, 847, 3 Jackson J. M., Geis N., Genzel R., Harris A. I., Madden Poglitsch A., Stacey G. J., Townes C. H., 1993, ApJ, 173                                                                                                                                                   |
| H., Inutsuka S.-I., 2000, ApJ, 532, 980 J. M. D., Dale J. E., Longmore S. N., 2015, 447, 1059 J. M. D., et al., 2019, arXiv e-prints,                                                                                                                                                                                                                                                       |
| S., 402, MNRAS,                                                                                                                                                                                                                                                                                                                                                                             |
| Kruijssen Kruijssen                                                                                                                                                                                                                                                                                                                                                                         |
| Larson R. B., 1981, MNRAS, 194, 809                                                                                                                                                                                                                                                                                                                                                         |
| Larson R. B., 1985, MNRAS, 214, 379                                                                                                                                                                                                                                                                                                                                                         |
| Koyama                                                                                                                                                                                                                                                                                                                                                                                      |
| Lau R. M., Herter T. L., Morris M. R., Becklin E. E.,                                                                                                                                                                                                                                                                                                                                       |
| Adams                                                                                                                                                                                                                                                                                                                                                                                       |
| J. D., 2013, ApJ, 775, 37                                                                                                                                                                                                                                                                                                                                                                   |

| Lis D. C., Carlstrom J. E., 1994, ApJ, 424, 189 Liszt H. S., Burton W. B., van der Hulst J. M., 1985, A&A, 142,                                                                                                                                                                                                     |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| 237 Liu H. B., Hsieh P.-Y., Ho P. T. P., Su Y.-N., Wright M., Sun A.-L., Minh Y. C., 2012, ApJ, 756, 195 Lo K. Y., Claussen M. J., 1983, Nature, 306, 647 Lu J. R., Do T., Ghez A. M., Morris M. R., Yelda S., K., 2013, ApJ, 764, 155                                                                              |
| Matthews                                                                                                                                                                                                                                                                                                            |
| Lu X., Zhang Q., Kauffmann J., Pillai T., Longmore S. N., Krui- jssen J. M. D., Battersby C., Gu Q., 2015, ApJ, 814, L18 Lu X., et al., 2017, ApJ, 839, 1                                                                                                                                                           |
| Lucas W. E., Bonnell I. A., Davies M. B., Rice W. K. M., 2013, MNRAS, 433, 353 Mapelli M., Gualandris A., 2016, in Haardt F., Gorini V., Moschella U., Treves A., Colpi M., eds, Lecture Notes in Physics, Berlin Springer Verlag Vol. 905, Lecture Notes in                                                        |
| Physics, Berlin Springer Verlag. p. 205 ( arXiv:1505.05473 ), doi:10.1007/978-3-319-19416-5 6 Mapelli M., Trani A. A., 2016, A&A, 585, A161 Mapelli M., Hayfield T., Mayer L., Wadsley J., 2012, ApJ, 749, 168 ı S., Mart´n-Pintado ı J., Montero-Casta˜o n M., Ho P. T. P., Blundell R., 2012, A&A, 539, A29       |
| Mart´n R. S., Coil A. L., Ho P. T. P., 2001, ApJ, 559, 326                                                                                                                                                                                                                                                          |
| McGary Mills E. A. C., G¨sten u R., Requena-Torres M. A., Morris M. R., 2013, ApJ, 779, 47                                                                                                                                                                                                                          |
| Montero-Casta˜o n M., Herrnstein R. M., Ho P. T. P., 2009, ApJ,                                                                                                                                                                                                                                                     |
| Morris M., Serabyn E., 1996, ARA&A, 34, 645 Moser L., et al., 2017, A&A, 603, A68 Novak G., Dotson J. L., Dowell C. D., Hildebrand R. H., Ren- barger T., Schleuning D. A., 2000, ApJ, 529, 241 Oka T., Nagai M., Kamegai K., Tanaka K., 2011, ApJ, 732, 120 Oka T., Onodera Y., Nagai M., Tanaka K., Matsumura S., |
| Kamegai K., 2012, ApJS, 201, 14 Okumura S. K., et al., 1989, ApJ, 347, 240 Okumura S. K., Ishiguro M., Fomalont E. B., Hasegawa T., Ka- suga T., Morita K.-I., Kawabe R., Kobayashi H., 1991, ApJ, 378, 127                                                                                                         |
| Paumard T., Maillard J.-P., Morris M., 2004, A&A, 426, 81                                                                                                                                                                                                                                                           |
| Paumard T., et al., 2006, ApJ, 643, 1011 Sanders R. H., 1998, MNRAS, 294, 35 Sandqvist A., 1989, A&A, 223, 293                                                                                                                                                                                                      |
| Sch¨del o R., et al., 2007, A&A, 469, 125                                                                                                                                                                                                                                                                           |
| Sch¨del o R., Merritt D., Eckart A., 2009, A&A, 502, 91 Sch¨del o R., Feldmeier A., Kunneriath D., Stolovy S., N., Amaro-Seoane P., Nishiyama S., 2014, A&A, 566, A47                                                                                                                                               |
| Serabyn E., Lacy J. H., 1985, ApJ, 293, 445                                                                                                                                                                                                                                                                         |
| Neumayer Serabyn E., Lacy J. H., Achtermann J. M., 1992, ApJ, 395, 166                                                                                                                                                                                                                                              |
| Shukla H., Yun M. S., Scoville N. Z., 2004, ApJ, 616, 231 Takekawa S., Oka T., Tanaka K., 2017, ApJ, 834, 121                                                                                                                                                                                                       |
| Teyssier R., 2002, A&A, 385, 337                                                                                                                                                                                                                                                                                    |
| Teyssier R., 2015, ARA&A, 53, 325                                                                                                                                                                                                                                                                                   |
| Trani A. A., Mapelli M., Ballone A., 2018, ApJ, 864, 17                                                                                                                                                                                                                                                             |
| Tsuboi M., Kitamura Y., Miyoshi M., Uehara K.,                                                                                                                                                                                                                                                                      |
| Tsutsumi Miyazaki A., 2016, PASJ, 68, L7 M., Kitamura Y., Uehara K., Tsutsumi T., Miyawaki                                                                                                                                                                                                                          |
| T., Tsuboi R., Miyoshi M., Miyazaki A., 2018, PASJ, 70, 85 Turk M. J., Smith B. D., Oishi J. S., Skory S., Skillman S.                                                                                                                                                                                              |
| Abel T., Norman M. L., 2011, The Astrophysical Journal                                                                                                                                                                                                                                                              |
| W., Sup- plement Series, 192, 9 B., Duschl W. J., 2001, A&A, 377, 1016                                                                                                                                                                                                                                              |
| B., Duschl W. J., 2002, A&A, 388, 128                                                                                                                                                                                                                                                                               |
| Vollmer Vollmer Vollmer B., Zylka R., Duschl W. J., 2003, A&A, 407, 515                                                                                                                                                                                                                                             |
| Wardle M., Yusef-Zadeh F., 2008, ApJ, 683, L37                                                                                                                                                                                                                                                                      |
| M. R.,                                                                                                                                                                                                                                                                                                              |
| Yelda S., Ghez A. M., Lu J. R., Do T., Meyer L., Morris                                                                                                                                                                                                                                                             |

- Zhao J.-H., Morris M. R., Goss W. M., An T., 2009, ApJ, 699, 186
- Zhao J.-H., Blundell R., Moran J. M., Downes D., Schuster K. F., Marrone D. P., 2010, ApJ, 723, 1097
- Zylka R., Mezger P. G., Wink J. E., 1990, A&amp;A, 234, 133

## APPENDIX A: A FURTHER, TEMPTATIVE CONNECTION WITH THE +50 KM/S CLOUD

As already discussed in section 4 and shown in Fig. 6, when the cloud is passing pericentre it may appear as two distinct clouds, connected close to the CND. In particular, Fig. 6 reminds, in projection, the actual positions and sizes of both the +20 km/s and +50 km/s clouds (cf. with Novak et al. 2000; Ferri`re 2012; Mapelli &amp; Gualane dris 2016). This plot was obtained from the same simulation used for our main scenario, but choosing a different projection. In this case, we had to choose ( x, y, z ) ˆ s ≈ (0 400 . , -0 343 . , -0 850), ( . x, y, z ) ˆ l ≈ ( -0 404 . , -0 898 0 172) . , . and ( x, y, z ) ˆ b ≈ (0 823 . , -0 274 0 498). In this different pro-. , . jection, the three orbital orientation angles (see section 3) are i ≈ 58 3 . ◦ , Ω ≈ 253 1 . ◦ and ω ≈ 36 5 . ◦ .

Such alternative configuration has several interesting features, but also shows inconsistencies with the properties of the observed molecular complex. In this Appendix we will discuss them and we will propose few speculations to possibly reconcile this more complex scenario with the observations.

For example, in this case the +20 km/s cloud has its closest approach to the SMBH in the north of the CND. Recent observations might show indications of this being the case (cf. with the bridge in Takekawa et al. 2017 and the extW streamer in Hsieh et al. 2017). On the other hand, as already discussed in section 4, most of the current models connect the +20 km/s and +50 km/s clouds on the southeast. In principle, this would also be possible in a single plunge scenario. However, compared to the configuration shown in Fig. 6, a different projection should be adopted and the infalling cloud should have a more circular orbit, with a much larger pericentre distance (comparable to the distance of the observed molecular ridge), in order to undergo a southern passage. In this case, it would be harder to explain the southern streamer and the formation of the CND with the plunge of these clouds. Another possibility is that the observed velocity pattern of the two molecular clouds (which seem to be connected by a smooth velocity gradient) is overlapping and confused with the overall motion of the central molecular zone at larger distances (see, e.g., Henshaw et al. 2016).

As also visible in the zoom-in in Fig. A1, in this different configuration, a stream of gas also appears to connect the +50 km/s and the CND. This might, in projection, be consistent with the so-called 'northern ridge' (e.g., McGary et al. 2001; Ferri`re 2012, and references therein), even e though it is not matching its observed kinematical properties. In fact, this gas in the simulation has high ( &gt; +50 km/s) line-of-sight velocities, while the observed northern ridge has a bulk velocity of roughly -10 km/s, even though it seems to connect to the CND with a velocity gradient reaching +60 km/s (McGary et al. 2001). We must stress that the region between the CND and the +50 km/s cloud

Figure A1. Zoom-in of Fig. 6, around the inner 16 pc of the domain. The red star labels the position of the SMBH.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-c6df12be-f5bd-4978-80b6-79e0ac532fd6.jpg)

.

Figure A2. Longitude-velocity diagram for the configuration in Fig. 6. The colormap shows the mass in every pixel of the diagram, scaled to the maximum value. At this time our simulation is able to reproduce some of the features of the observed CND, +20 km/s cloud and +50 km/s complex.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-ee5e8e78-f7fd-4f11-a9a5-b958f3d18034.jpg)

is believed to be strongly affected by the supernova remnant SgrA East (see section 1), so it is also difficult (and maybe deceptive) to strictly compare it to our simulation, which does not include such supernova explosion.

Fig. A2 shows a longitude-velocity diagram for this different configuration. The black line shows the connection between the pre-pericentre and the post-pericentre parts of the cloud. Due to the fast orbital precession, the material that has passed pericentre keeps on moving with positive line-of-sight velocity. The post-pericentre gas has a bulk velocity that is a bit higher than +50 km/s and a large spread. A lot of gas is actually at velocities of around 150 km/s, which is much more than observed. Interestingly enough, this longitude-velocity diagram also shows motion of gas that is not fully coherent with the motion of the cloud and the CND (see the magenta dashed lines). These might be compared with the longitude-velocity diagrams of the extW and eastern streamers in Hsieh et al. (2017) (Fig. 5, 6 and 8 of their paper; see also the longitude-velocity diagrams in Fig. 2 of Oka et al. 2011). Hence, in this more complex configuration, the extW streamer could be part of the +20 km/s cloud reaching the CND from 'above', while (as in the model of Hsieh et al. 2017) the eastern streamer could be material of the +50 km/s cloud, falling back towards the CND.

Nonetheless, the major problem of this simulated diagram is the too high velocity observed on the side of the +50 km/s cloud. This might be due to drag of part of the CND (see the dotted branch of the red line in Fig. A2). We must remind the reader that this is a line-of-sight velocity, thus it also depends on the actual trajectory followed by the cloud. The two observed clouds have relatively small line-of-sight velocities, despite their distance from the Galactic Centre. So, for such an eccentric orbit, lower line-of-sight velocities for the +50 km/s cloud could be obtained if the cloud trajectory lay closer to the plane of the sky. A faster precessing orbit could bring the second apocentre position away from the line-of-sight. This might be obtained by a higher stellar mass enclosed by the orbit, compared to our assumptions. The current precession angle is roughly 120 ◦ and, on a 1st order approximation, the precession angle scales linearly with the enclosed mass (see, e.g., Gualandris et al. 2012). Hence, the enclosed stellar mass should be higher within a factor ≈ 1 5, which, nonetheless, seems to be within the uncertain-. ties of the currently available observations (Chatzopoulos et al. 2015; Fritz et al. 2016; Feldmeier-Krause et al. 2017; Gallego-Cano et al. 2018).

In terms of line-of-sight distance from SgrA*, Vollmer et al. (2003) put the +50 km/s cloud in front of SgrA*, at a maximum line-of-sight distance of -5 pc, while Ferri`re e (2012) prefers to place the +50 km/s cloud behind SgrA*, at a distance of 3 pc. As visible in Fig. A3, the +20 km/s cloud lies roughly at the same distance from SgrA* as for our main scenario (cfr. Fig. 7), while the +50 km/s lies behind it, but it extends along the line-of-sight up to 15-20 pc from SgrA*, which is in disagreement with both the aforementioned estimates. Again, a faster precessing orbit might bring the post-pericentre trajectory of the infalling cloud closer to the plane of the sky, eventually giving the +50 km/s cloud a smaller line-of-sight extent and distance from SgrA*.

This paper has been typeset from a T E X/LT E X file prepared by A the author.

## 10 A. Ballone et al.

Figure A3. Density projection along the Galactic latitude in the new reference system for the configuration in Fig. 6. The axes represent the Galactic longitude (in degrees) and the line-of-sight distance from SgrA* (in pc). The red star labels the position of the SMBH.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00549v1-figure-6904b05e-8344-4a02-8648-627ca33a10dd.jpg)
Alessandro Ballone Michela Mapelli Alessandro Alberto Trani
A common origin for the circumnuclear disc and the nearby molecular clouds in the Galactic Centre
The origin of the molecular clouds orbiting SgrA* is one of the most debated
questions about our Galactic Centre. Here, we present the hydrodynamic
simulation of a molecular cloud infalling towards SgrA*, performed with the
adaptive-mesh-refinement code RAMSES. Through such simulation, we propose that
the circumnuclear disc and the +20 km/s cloud originated from the same tidal
disruption episode, occurred less than 1 Myr ago. We also show that recent star
formation is to be expected in the +20 km/s cloud, as also suggested by recent
observations.