Dancing with the stars: Stirring up extraordinary turbulence in Galactic center clouds

Molecular clouds in the central molecular zone (CMZ) have been observed to feature turbulent line widths that are significantly higher, and scale with cloud size more steeply, than in the rest of the Milky Way. In the same Galactic region, the stellar density is also much higher than in the rest of the Milky Way, and the vertical stellar velocity dispersion is large, meaning that even young stars are likely to cross the entire vertical extent of the CMZ within their lifetimes. Here, we investigate whether interactions of CMZ molecular clouds with crossing stars can account for the extraordinary properties of observed turbulence in this part of the Galaxy. We calculated the rate of energy deposition by stars crossing CMZ clouds due to (a) stellar winds and (b) dynamical friction, and compared it to the rate of turbulence decay. We calculated the predicted scaling of turbulence line width with cloud size in each case. We find that energy deposition by stellar winds of crossing massive stars can account for both the level and the scaling of CMZ cloud turbulence with cloud size. We also find that the mechanism stops being effective at a Galactocentric distance comparable to the CMZ extent. On the other hand, we find that dynamical friction by crossing stars does not constitute a significant driver of turbulence for CMZ clouds.


Introduction
The central molecular zone (CMZ) is an extreme environment in the Milky Way. Turbulence in molecular clouds is extraordinarily high, with line widths of 10-100 km/s, and a scaling with effective cloud size steeper than in the rest of the Milky Way (σ v ∝ L 0.66±0.18 , Kauffmann et al. 2017a). The driver of this excess turbulence is unclear. Cloud-cloud collisions have been suggested as a possible cause, and the shocks associated with such collisions could also explain the shock chemical tracers found in the Galactic center region (Menten et al. 2009;Lis et al. 2001). However, it is not evident that cloud-cloud collisions could also offer a natural explanation for the steep line-widthsize scaling of CMZ turbulence, while the shock chemical tracers could also be the result of hot expanding shells originating in Wolf-Rayet stars (Martín-Pintado et al. 1999).
The stellar density in the Galactic center regions is also very high, including a significant population of massive stars (Bryant & Krabbe 2021;Rosslowe & Crowther 2015;Kauffmann et al. 2017a). Stars in the CMZ show significant vertical velocity dispersion (∼ 50 km/s; e.g., Kim & Morris 2001). Given that an object moving at 1 km/s crosses 1 pc in 1 Myr, these velocities are sufficiently high even for massive stars to cross the entire vertical extent of the CMZ within their lifetimes.
In this paper, we investigate whether or not stars crossing molecular clouds could be responsible for maintaining the high levels of turbulence in the CMZ, and if so, by which mecha-⋆ tassis@physics.uoc.gr ⋆⋆ pavlidou@physics.uoc.gr nism. The interaction between stars and molecular clouds near the Galactic center has been investigated for its effect on the stellar population; Kim & Morris (2001) for example, investigated whether the scattering of stars off of giant molecular clouds in the CMZ can drive the vertical diffusion of newly formed stars in the region. Here, we are looking into the reverse process, that is, whether the presence of a high number of stars affects clouds by driving turbulence within them, via (a) direct kinetic energy deposited by crossing stars that feature fast stellar winds or (b) dynamical friction.

Energetics of turbulence in the CMZ
Turbulence in molecular clouds is supersonic and generates shocks that dissipate the energy of the motions on a timescale of the order of a cloud-crossing time (Stone et al. 1998). To maintain a certain level of turbulence, a driver is therefore needed that will replenish the lost energy at a comparable rate. The decay rate of kinetic energy in turbulent motions of typical velocity v turb in a cloud of average density ρ and size L isĖ turb = 1 2 ρL 3 v 2 turb /t cross , where t cross = L/v turb is the turbulence crossing time. This then givesĖ turb = 1 2 ρL 2 v 3 turb or, by normalizing to values typical for the CMZ, (1) Article number, page 1 of 3 A&A proofs: manuscript no. DwtS

Energy deposition by stellar winds in the CMZ
The Galactic center is a region of high stellar density, with a high abundance of stars that produce significant winds, of the order of 10 −3 M ⊙ /yr within the central parsec (Coker 2001). Clouds in the CMZ will interact with these stars, both young and older. Stars will have significant relative velocities with respect to molecular clouds, and will be crossing these clouds, stirring turbulence in them in two ways: dynamical friction, and cloudwind interaction. In this section, we calculate the latter. Between the central parsec and a Galactocentric distance R ∼ 100 − 300 pc typical of the CMZ , the self-consistent model of Sormani et al. (2021) for the inner Milky Way bulge predicts a decline of the stellar population by a factor of ∼ 10 3 − 10 4 . Assuming the integrated wind mass-loss rate will decline by the same factor, and adopting the conservative estimate of a ×10 4 decline, we expect a typical mass deposition due to stellar winds of 10 −7 M ⊙ yr −1 pc −3 . For stars that happen to be crossing clouds in the CMZ, these winds will deposit a kinetic energy of ∼ 5 × 10 42 erg yr −1 pc −3 (using a conservative average energy deposition rate over the lifetimes of large stars of ∼ 10 49 erg per M ⊙ of lost mass 1 from the numerically integrated models of Fierlinger et al. 2016, who used tabulated mass-loss rates from Ekström et al. 2012 and wind velocities from Voss et al. 2009). Normalizing again, as in Eq. 1, to a cloud of size ∼ 1pc, this implies an energy deposition rate by winds oḟ which, very suggestively, lies at exactly the right level to replenish turbulence decaying at the rate of Eq. (1). The least secure assumption in deriving Eq.
(2) is that the number of massive stars producing strong winds scales in the same way as the total stellar mass between the central parsec and the typical CMZ radii. For this reason, we also present an alternative calculation ofĖ kin,winds with different, although still large, uncertainties.
From Fierlinger et al. (2016) (their Fig. 1), the (averaged over stellar lifetime) kinetic energy injection rate of a single star as a function of stellar mass can be approximately described -for stars more massive than 20 M ⊙ -byĖ kin,1 (m) ≈ (3.5 × 10 37 erg/yr) m 4 , where m is the stellar mass in M ⊙ . From the same figure, stellar lifetimes in the same mass range (including the post-main sequence wind stage) can be approximated by If stars are produced in the CMZ at a rate ofṀ * (in M ⊙ /yr, assumed constant over the past few tens of millions of years) with a power-law initial mass function 2 (IMF) dn/dm = Cm −s , then the number of (massive) stars with masses between m and m+dm will be T (m)(dn/dm)dm and their kinetic energy injection rate over the entire CMZ will beĖ kin,1 T (m)(dn/dm)dm. Integrating over all stars with m > 20, we obtain a totalĖ kin,winds from such stars of m max m minĖ kin,1 T (m)(dn/dm)dm. Adopting standard parameters for a Salpeter IMF (s = 2.35, m max = 120, and m min = 1 From a comparison of top and bottom panels of Fig. 1 of Fierlinger et al. (2016) it follows that kinetic energy injection rate and mass-loss rate scale with stellar mass in approximately the same way for massive stars (> 20 M ⊙ ), implying a wind velocity roughly constant with mass for these stars; this is why in Eq. (2) it is sufficient to scale the kinetic energy injection rate with the mass wind loss rate alone. 2 We normalize the IMF so that 0.15 so that m ∼ 0.5), and aṀ * ∼ 0.1 (e.g. Henshaw et al. 2022 and references therein) we obtainĖ kin,winds ∼ 10 40 erg/s for the entire volume of the CMZ to which the star formation estimate refers. Converting this to a per-volume estimate is not trivial. First, the star formation is not uniformly distributed over the CMZ volume, as is clear from Henshaw et al. (2022), who compare different estimates ofṀ * in the CMZ using different tracers and studying different volumes: as can be seen in their Table 1, there is no strong correlation between the volume of the CMZ studied and the derived value ofṀ * , which is to be expected if the star formation, massive stars, and molecular clouds are all spatially correlated and only cover a small fraction of the CMZ volume. If we assume that most of this star formation takes place in the same region of the central CMZ covered by the clouds studied in Kauffmann et al. (2017a) (radius of ∼ 100 pc, height of ∼ ±5 pc from the Galactic plane), then we obtaiṅ E kin,winds ∼ 4 × 10 34 erg s −1 pc −3 (corresponding to a population of about 1500 O and B stars of all evolutionary stages within the CMZ). However, this estimate is very uncertain. It is most sensitive to our assumptions on the IMF and the volume within most of the star formation is confined. For example, if we assume a top-heavy IMF similar to that preferred for the central parsec (e.g. s = 0.85, m min = 0.7 as in Model 8 of Maness et al. 2007), we obtainĖ kin,winds ∼ 6 × 10 35 erg s −1 pc −3 . Conversely, for every factor-two increase in the Galactrocentric radius of the CMZ zone, where most of the star formation takes place, the estimate ofĖ kin,winds is decreased by a factor of four. However, it is encouraging that, even with these large uncertainties, this first-principles estimate is of the same order of magnitude as the simple calculation of Eq. (2).

Scaling of turbulent line widths with cloud size
If the turbulence driving of Eq. (2) balances the decay of Eq. (1), then a specific scaling of turbulent velocities, v turb , with cloud or clump size, L, is predicted. Indeed, if the cloud density scales with radius as ρ ∝ r −k , then the average number density within a cloud or clump of typical size L will also scale as n H 2 ∝ L −k . ThenĖ kin,winds ∼Ė turb implies L 2−k v 3 turb ∝ L 3 or v turb ∝ L (1+k)/3 .
(3) Kauffmann et al. (2017b) find, for the majority of structures in the CMZ, k ∼ 1.3, so Eq. (3) gives v turb ∝ L 0.77 , in excellent agreement with the observed scaling in the CMZ reported by Shetty et al. (2012) (most-probable scaling through Bayesian inference of v turb ∝ L 0.78 with a 95% confidence interval on the slope from 0.41 to 1.13) and by Kauffmann et al. (2017a) (v turb ∝ L 0.66±0.18 ).

Galactocentric distance cutoff of enhanced turbulence due to stellar winds
In our discussion, we assume that the number and distribution of massive, fast-wind stars scales with the total stellar density between the central parsec and the outer edge of the CMZ. Interestingly, young stars are observed throughout the CMZ (Bryant & Krabbe 2021). Kauffmann et al. (2017a) do observe O-stars and/or masers in several of the CMZ clouds that they have investigated, but not in all of them. It is therefore interesting to ask whether or not it is necessary to have at least one highmass star currently present in every cloud observed to be abnormally turbulent, if the excess turbulence is indeed produced by winds from crossing stars. The answer is no: the present level K. Tassis and V. Pavlidou 2022: Cloud turbulence from stars at Galactic Center of turbulence in a cloud could have been driven by a crossing massive star in the past. Stars near the Galactic center have vertical velocity dispersions that are significantly higher (∼ ×5) than typical turbulence velocities in clouds, and therefore the star crossing time will be shorter than the turbulence crossing time (on which turbulence decays) by the same factor. Clouds with extraordinary turbulence driven by this mechanism could spend most of their lifetimes devoid of high-mass, high-wind stars, which is consistent with the lack of evidence for the presence of massive stars in several CMZ clouds. We can further quantify this argument, to predict the Galactocentric distance beyond which this mechanism will no longer be efficient in stirring up excess turbulence in molecular clouds. This will occur once the number of stars with significant winds making a cloud crossing over a turbulence decay time L/v turb drops below 1. The number of crossings through a cloud of size L by wind-producing stars of number density n ⋆ and crossing velocity v ⋆ in this time is n ⋆ v ⋆ L 2 (L/v turb ). For the mechanism to remain effective, we need this to be at least 1, that is, we demand that n ⋆ L 3 ≥ v turb /v ⋆ . If the suppression factor of stellar density compared to the central parsec as a function of Galactocentric distance R is f (R), and the central parsec has ∼ 10 2 stars producing significant fast winds, then this mechanism will stop being effective at a Galactocentric distance R such that f (R) ≥ (L/0.2pc) −3 (v turb /50 km s −1 ). For clouds of L ∼ 3pc in size with v turb ∼ 10 km s −1 (Kauffmann et al. 2017a), we require f (R) 10 −4 , which corresponds to a Galactocentric distance of ∼ 250pc (Sormani et al. 2021), which is comparable to the extent of the CMZ.

Energy deposition through dynamical friction
That the level and scaling of the kinetic energy injection by stellar winds matches the observations of turbulence in CMZ clouds is a non-trivial result. To emphasise this point, we present here a calculation of an alternative source of kinetic energy injection.
As stars cross through molecular clouds, they experience gravitational drag, that is, a loss of momentum and energy due to gravitational interaction with the surrounding matter. The phenomenon, known as dynamical friction, was first studied by Chandrasekhar (1943). The rate of change of the (relative to the cloud) velocity vector of a star, v ⋆ , due to dynamical friction, is (see Binney et al. 1987, chapter 8 where M ⋆ is the mass of the star; ρ is the gas density; G is the gravitational constant; X = v ⋆ /( √ 2σ), where σ is the velocity dispersion of the gas; and Λ = b max v 2 0 /GM ⋆ , where b max is the maximum distance that needs to be considered (of the order of the size of the cloud, and so we take b max ≈ L/2) and v 0 ≈ v ⋆ is the initial relative velocity vector of the star. The rate of energy deposit in the cloud by a single star isĖ ⋆,1 = M ⋆ v ⋆ dv ⋆ dt and the total energy deposition rate by stars due to dynamical friction isĖ ⋆ =Ė ⋆,1 n ⋆ L 3 where n ⋆ is the number density of stars. For L between 1 and 100 pc, and a typical v ⋆ of 50 km/s (e.g. Kim & Morris 2001), ln Λ will vary between 12 and 17; we adopt ln Λ ∼ 15. As v ⋆ is much higher than σ, X will be large, and so the term in the square brackets will be of order unity. For Galactocentric distances 100 pc R 300 pc, Sormani et al. (2021) modeled a local stellar density of 10 2 M ⊙ pc −3 ρ ⋆ 10 M ⊙ pc −3 . Normalizing to a 1 pc cloud, as in Eqs. (1) and (2), and adopting an average stellar mass of 0.5M ⊙ , we finḋ E ⋆ = 8 ×10 29 erg/s ρ ⋆ 10 2 M ⊙ pc −3 n H 2 10 4 cm −3 L 1 pc over five orders of magnitude too low to be relevant for maintaining the observed levels of turbulence in CMZ clouds. Even if the level of energy deposition rate by dynamical friction were comparable to the turbulence decay rate, the predicted scaling in this case would be incompatible with observations: from comparing Eqs. (1) and (5) we obtain L 2 n H 2 v 3 turb ∝ n H 2 L 3 , implying v turb ∝ L 1/3 , much shallower than the observed size-linewidth relation in the CMZ (Shetty et al. 2012;Kauffmann et al. 2017a). We conclude that dynamical friction of clouds on stars does not contribute appreciably to the observed level of turbulence in the CMZ.

Conclusions
We show that fast stellar winds from massive stars crossing CMZ clouds can nicely reproduce the observed properties of turbulence in CMZ clouds, including its magnitude and scaling with cloud size. We also show that this mechanism ceases to be effective at a Galactocentric radius comparable to the extent of the CMZ. In contrast, we find that dynamical friction is not important as a turbulence driver for CMZ clouds in magnitude, nor does it reproduce the observed scaling between line width and size.