# Spatial Segregation of Massive Clusters in Dwarf Galaxies

Bruce Elmegreen, A. Adamo, M. Boquien, F. Bournaud, D. Calzetti, D. Cook, D. Dale, P. -A. Duc, D. Elmegreen, J. Fensch, et al.

► To cite this version:

Bruce Elmegreen, A. Adamo, M. Boquien, F. Bournaud, D. Calzetti, et al. Spatial Segregation of Massive Clusters in Dwarf Galaxies. The Astrophysical Journal Letters, 2020, 888 (2), pp.L27. 10.3847/2041-8213/ab632a. hal-03142652

HAL Id: hal-03142652

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

Submitted on 28 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>

# Spatial Segregation of Massive Clusters in Dwarf Galaxies

Bruce G. Elmegreen<fnref n="1" />, A. Adamo<fnref n="2" />, M. Boquien<fnref n="3" />, F. Bournaud<fnref n="4" />, D. Calzetti<fnref n="5" />, D. O. Cook<fnref n="6" />, D. A. Dale<fnref n="7" />, P.-A. Duc<fnref n="4" />,<fnref n="8" />, D. M. Elmegreen<fnref n="9" />, J. Fensch<fnref n="10" />, K. Grasha<fnref n="11" />,<fnref n="12" />, Hwi Kim<fnref n="13" />, L. Kahre<fnref n="14" />, M. Messa<fnref n="5" />, J. E. Ryon<fnref n="15" />, E. Sabbi<fnref n="15" />, and L. J. Smith<fnref n="16" />

Received 2019 December 2; revised 2019 December 13; accepted 2019 December 16; published 2020 January 16

## Abstract

The relative average minimum projected separations of star clusters in the Legacy ExtraGalactic UV Survey (LEGUS) and in tidal dwarfs around the interacting galaxy NGC 5291 are determined as a function of cluster mass to look for cluster-cluster mass segregation. Class 2 and 3 LEGUS clusters, which have a more irregular internal structure than the compact and symmetric class 1 clusters, are found to be mass-segregated in low-mass galaxies, which means that the more massive clusters are systematically bunched together compared to the lower-mass clusters. This mass segregation is not present in high-mass galaxies or class 1 clusters. We consider possible causes for this segregation, including differences in cluster formation and scattering in the shallow gravitational potentials of low-mass galaxies.

Unified Astronomy Thesaurus concepts: Galaxies (573); Star formation (1569); Star clusters (1567)

## 1. Introduction

Compact star clusters usually form inside more extended associations of young stars (Feitzinger & Galinski 1987; Elmegreen & Efremov 1997; Hunter 1999; Maíz-Apellániz 2001; Lada & Lada 2003; Elmegreen et al. 2006; Elmegreen 2008; Portegies Zwart et al. 2010; Krumholz et al. 2018), as part of a hierarchical structure for star formation that resembles the distribution of dense interstellar clouds (Scalo 1985; Fleck 1996; Elmegreen & Falgarone 1996; Cartwright & Whitworth 2004). The relative positions and ages of these clusters follow power-law correlations (Efremov & Elmegreen 1998; de la Fuente Marcos & de la Fuente Marcos 2009; Grasha et al. 2017b), suggesting this hierarchy is the result of turbulent motions with self-gravity dominating the densest phase.

While this basic structure is well observed, there has been little effort to quantify the spatial correlation as a function of cluster mass. We do not know, for example, if the most massive clusters group together with lower-mass clusters surrounding them. Such mass segregation can be an important constraint on cluster formation models and an indicator of the history of the region, including competitive (e.g., Zinnecker 1982; Bonnell et al. 1997) or cooperative (Elmegreen et al. 2014) accretion of gas into the clusters, or density-dependent cloud masses (Alfaro & Román-Zúñiga 2018). Mutual cluster attraction leading to coalescence in dense regions (Lahén et al. 2019) might also be indicated.

Here we describe a new metric for mass-dependent clustering that has the potential to reveal whether clusters segregate according to mass. We apply this metric to clusters in the Legacy ExtraGalactic UV Survey (LEGUS) database (Calzetti et al. 2015), and, by comparison, to clusters in tidal dwarfs around the interacting galaxy NGC 5291 (Fensch et al. 2019). These tidal dwarfs have *Hubble Space Telescope (HST)* images and the highest number of clusters observed so far for this type of galaxy.

## 2. Method

### 2.1. *The Relative Average Minimum Projected Separation*

To gain insight into mass segregation among high- and low-mass clusters, we consider the average minimum projected separation between clusters as a function of cluster mass. We denote this quantity, corrected for galaxy inclination, by $\bar{D}_{\min}(M)$, where the *bar* denotes the average for all clusters of a particular mass *M*, and the subscript “min” denotes the minimum distance *D* to these other clusters. For a random distribution of cluster positions, this separation is about equal to the inverse square root of the average projected cluster density. In addition, we denote the number of clusters in a logarithmic mass interval by $N(M)d \log M$. For a uniform random distribution of clusters of all masses, $\bar{D}_{\min}(M)$ multiplied by $N(M)^{1/2}$ is independent of *M*. The mass distribution function for clusters and incompleteness at low-mass both enter $\bar{D}_{\min}(M)$ and $N(M)$ in the same way, canceling out.

<footnote n="1">¹IBM Research Division, T.J. Watson Research Center, Yorktown Hts., NY 10598, USA</footnote>

<footnote n="2">²Department of Astronomy, The Oskar Klein Centre, Stockholm University, Stockholm, Sweden</footnote>

<footnote n="3">³Centro de Astronomía (CITEVA), Universidad de Antofagasta, Avenida Angamos 601, Antofagasta, Chile</footnote>

<footnote n="4">⁴AIM, CEA, CNRS, Université Paris-Saclay, Université Paris Diderot, Sorbonne Paris Cité, F-91191 Gif-sur-Yvette, France</footnote>

<footnote n="5">⁵Department of Astronomy, University of Massachusetts Amherst, Amherst, MA 01003, USA</footnote>

<footnote n="6">⁶California Institute of Technology, Pasadena, CA, USA</footnote>

<footnote n="7">⁷Department of Physics and Astronomy, University of Wyoming, Laramie, WY, USA</footnote>

<footnote n="8">⁸Université de Strasbourg, CNRS, Observatoire astronomique de Strasbourg, UMR 7550, F-67000 Strasbourg, France</footnote>

<footnote n="9">⁹Vassar College, Department of Physics and Astronomy, Poughkeepsie, NY, USA</footnote>

<footnote n="10">¹⁰Univ. Lyon, ENS de Lyon, Univ. Lyon 1, CNRS, Centre de Recherche Astrophysique de Lyon, UMR5574, F-69007 Lyon, France</footnote>

<footnote n="11">¹¹Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia</footnote>

<footnote n="12">¹²ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia</footnote>

<footnote n="13">¹³NSF's National Optical-Infrared Astronomy Research Laboratory/Gemini Observatory, Casilla 603, La Serena, Chile</footnote>

<footnote n="14">¹⁴Department of Astronomy, New Mexico State University, Las Cruces, NM, USA</footnote>

<footnote n="15">¹⁵Space Telescope Science Institute, Baltimore, MD 21218, USA</footnote>

<footnote n="16">¹⁶European Space Agency/Space Telescope Science Institute, Baltimore, MD, USA</footnote>
If high-mass clusters of mass $M_H$ are more bunched together than low-mass clusters of mass $M_L$, then

$$N(M_H)^{1/2} \bar{D}_{\min}(M_H) < N(M_L)^{1/2} \bar{D}_{\min}(M_L). \quad (1)$$

Thus, a plot of $N(M)^{1/2} \bar{D}_{\min}(M)$ versus $M$ indicates the relative segregation of different masses.

For comparisons among different regions in a galaxy or different galaxies, the above quantity should be normalized to the region size, which we take to be the average projected separation (corrected for galaxy inclination) between the clusters in the lowest-mass interval, $\bar{D}(M_{low})$. This interval is chosen because generally it has the most clusters and gives the most accurate region size. Thus, the quantity to consider as a function of mass is the relative average minimum projected separation (RAMPS) corrected for galaxy inclination,

$$\text{RAMPS}(M) = N(M)^{1/2} \bar{D}_{\min}(M) / \bar{D}(M_{\text{low}}). \quad (2)$$

Another measure of relative cluster separation is the two-point correlation function (e.g., Bastian et al. 2009; Grasha et al. 2017a), which determines the relative proportions of all separations. The RAMPS differs because it uses only the nearest distances.

## 2.2. Testing the RAMPS

A fractal hierarchical model shows the trends in RAMPS. This model is made on a 512 × 512 square grid of total size 1 with 8 levels of hierarchy starting with a 4 × 4 grid of cells at the top in level $i = 2$ (see Elmegreen 2018). At each level $i$ there is a probability $p$ of choosing a cell that will be further subdivided into 2 × 2 cells in the next lower level. This probability depends on the fractal dimension $D_f$ and is given by $p = 2^{D_f-2}$. To choose a cell, a random number between 0 and 1 is generated and compared to $p$; if the random number is smaller than $p$, we choose the cell. Note that a fractal dimension $D_f = 2$ causes all cells to be chosen ($p = 1$), filling the square grid completely in two dimensions. For the model we use $D_f = 1.3$ because that matches the observations of interstellar clouds (Elmegreen & Falgarone 1996). Each level sub-divides only the cells chosen at the next higher level. For level $i$ from 2 to 9, the size of the cell is $1/2^i$ and the total number of cells is $2^i \times 2^i$, although only the fraction $p$ are chosen at each level.

We assign a mass to each of the cells at all levels, considering that the cell center represents the position of a star cluster in the hierarchy of young stellar structures. To give mass segregation, uniformity, or inverse segregation, we let the cluster mass $M$ scale with the level $i$ such that $M = M_0 \mu^i$ where $\mu = 10^{0.5}$, 1, or $10^{-0.5}$ for $i = 2, ..., 6$ in these three cases, respectively. To keep the masses in the range of $10^3 M_\odot-10^5 M_\odot$, we set $M_0 = 10^2 M_\odot$ and $10^6 M_\odot$ for the segregation and inverse segregation cases. For the uniform case, we let $\log M_0$ be a random number uniformly distributed between 3 and 5. Because cells at levels with larger $i$ are on average closer together within their hierarchies than cells at lower $i$ within the same hierarchy, the segregation of high-mass clusters to denser average regions corresponds to a greater proportion of massive clusters at large $i$. This is why $\mu > 1$ corresponds to mass segregation and $\mu < 1$ corresponds to inverse mass segregation, where massive clusters systematically avoid each other compared to low-mass clusters.

Figure 1 shows the RAMPS function for the three values of $\mu$. Mass segregation with $\mu > 1$ has a negative slope and inverse mass segregation has a positive slope. Random mass in

<figure><img src="image_3.png" /><figcaption>Figure 1. (Top) The relative average minimum projected separations (RAMPS) between clusters are plotted vs. the cluster mass for randomly positioned clusters in a hierarchical fractal distribution. The mass units are arbitrary. Increasing slopes indicate that low-mass clusters are preferentially clumped together, while decreasing slopes indicate that high-mass clusters are clumped together. Randomly positioned clusters give a flat slope.</figcaption></figure>

the hierarchy corresponds to a horizontal line in the figure. The way to interpret the negative slope is that the nearest neighbor high-mass cluster to another high-mass cluster is closer than the average cluster spacing would be at that mass if they were randomly distributed over the whole region.

## 3. Data

Catalogs in the HST LEGUS survey were used to obtain the positions, ages, and masses of measured clusters (Calzetti et al. 2015; Adamo et al. 2017), considering Padova stellar evolution models with starburst extinction curves. To keep the sample as free as possible from fading effects with age, we consider only clusters more massive than a distance-dependent limit, $M_{\text{limit}}$, and younger than 125 Myr. For almost all galaxies in LEGUS, the lower limit to the detectable mass at 125 Myr age is $94 D_{\text{Mpc}}^2 M_\odot$ for distance $D_{\text{Mpc}}$ in Mpc. This mass corresponds to an absolute V-band magnitude of $M_v = -6.0$. For a typical distance of 6 Mpc, the typical mass limit is 3400 $M_\odot$. Then, for the entire sample, there are 14 galaxies that have 35 or more such clusters in classes 2 or 3 and which span a factor of $10^{1.5}$ or more in mass. These galaxies and their cluster counts are listed in Table 1 along with the RAMPS slopes. Galaxy distances, star formation rates (SFRs), and stellar masses are from Calzetti et al. (2015), as are inclinations (not listed). Galaxy position angles are from various sources such as the Third Reference Catalog of Bright Galaxies (de Vaucouleurs et al. 1991) and the LITTLE THINGS survey (Hunter et al. 2012) when available, and measured from LEGUS images at the LEGUS website<fnref n="17" /> when not available, all verified by measurements on the Digitized Sky Survey.<fnref n="18" />

<footnote n="17">¹⁷ https://legus.stsci.edu/legus_observations.html</footnote>

<footnote n="18">¹⁸ http://archive.stsci.edu/cgi-bin/dss_form</footnote>
<table><caption>**Table 1**
Cluster Sample</caption><thead><tr><th>Galaxy<sup>a</sup></th><th><i>D</i><sup>b</sup></th><th>SFR</th><th><i>M</i><sub>star</sub></th><th><i>N</i><sup>d</sup></th><th>1<sup>c</sup><br><i>B</i></th><th><i>S</i></th><th><i>N</i></th><th>2<br><i>B</i></th><th><i>S</i></th><th><i>N</i></th><th>3<br><i>B</i></th><th><i>S</i></th><th><i>N</i></th><th>2+3<br><i>B</i></th><th><i>S</i></th></tr></thead><tbody><tr><td>628</td><td>9.90</td><td>3.67</td><td>11</td><td>133</td><td>4</td><td>0.02 ± 0.15</td><td>59</td><td>4</td><td>0.05 ± 0.28</td><td>44</td><td>3</td><td>-0.15 ± 0.53</td><td>103</td><td>4</td><td>0.07 ± 0.37</td></tr><tr><td>1249</td><td>6.90</td><td>0.15</td><td>0.55</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>43</td><td>4</td><td>0.28 ± 0.35</td></tr><tr><td>1313</td><td>4.39</td><td>1.15</td><td>2.6</td><td>113</td><td>4</td><td>-0.10 ± 0.13</td><td>250</td><td>4</td><td>0.01 ± 0.13</td><td>286</td><td>4</td><td>0.06 ± 0.16</td><td>536</td><td>4</td><td>-0.01 ± 0.06</td></tr><tr><td>1566</td><td>13.20</td><td>5.67</td><td>27</td><td>212</td><td>4</td><td>0.02 ± 0.26</td><td>119</td><td>3</td><td>0.14 ± 0.63</td><td>93</td><td>3</td><td>0.33 ± 0.47</td><td>212</td><td>3</td><td>0.09 ± 0.47</td></tr><tr><td>1705</td><td>5.10</td><td>0.11</td><td>0.13</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>36</td><td>4</td><td>-0.18 ± 0.92</td></tr><tr><td>3344</td><td>7.00</td><td>0.86</td><td>5.0</td><td>35</td><td>3</td><td>0.22 ± 0.30</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>3351</td><td>10.00</td><td>1.57</td><td>21</td><td>39</td><td>5</td><td>-0.64 ± 0.18</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>3627</td><td>10.10</td><td>4.89</td><td>31</td><td>144</td><td>3</td><td>0.11 ± 0.03</td><td>92</td><td>4</td><td>0.06 ± 0.11</td><td>50</td><td>4</td><td>0.23 ± 0.35</td><td>142</td><td>4</td><td>0.19 ± 0.32</td></tr><tr><td>3738</td><td>4.90</td><td>0.07</td><td>0.24</td><td>...</td><td>...</td><td>...</td><td>60</td><td>4</td><td>-0.43 ± 0.23</td><td>55</td><td>3</td><td>-0.13 ± 0.89</td><td>115</td><td>4</td><td>-0.36 ± 0.23</td></tr><tr><td>4395</td><td>4.30</td><td>0.34</td><td>0.60</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>44</td><td>4</td><td>-0.04 ± 0.35</td></tr><tr><td>4449</td><td>4.31</td><td>0.94</td><td>1.1</td><td>66</td><td>5</td><td>0.04 ± 0.31</td><td>160</td><td>5</td><td>-0.39 ± 0.13</td><td>132</td><td>4</td><td>-0.07 ± 0.20</td><td>292</td><td>5</td><td>-0.38 ± 0.21</td></tr><tr><td>4656</td><td>5.50</td><td>0.50</td><td>0.40</td><td>43</td><td>3</td><td>-0.61 ± 0.17</td><td>61</td><td>3</td><td>-0.31 ± 1.17</td><td>46</td><td>4</td><td>0.27 ± 0.20</td><td>107</td><td>5</td><td>-0.03 ± 0.31</td></tr><tr><td>5194</td><td>7.66</td><td>6.88</td><td>24</td><td>140</td><td>5</td><td>-0.06 ± 0.38</td><td>198</td><td>5</td><td>0.09 ± 0.11</td><td>81</td><td>3</td><td>0.51 ± 0.95</td><td>279</td><td>5</td><td>0.14 ± 0.09</td></tr><tr><td>5194-ML</td><td>7.66</td><td>6.88</td><td>24</td><td>610</td><td>5</td><td>-0.26 ± 0.29</td><td>646</td><td>4</td><td>-0.11 ± 0.16</td><td>89</td><td>3</td><td>0.36 ± 0.62</td><td>735</td><td>4</td><td>-0.08 ± 0.17</td></tr><tr><td>5253</td><td>3.15</td><td>0.10</td><td>0.22</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>93</td><td>5</td><td>-0.18 ± 0.21</td></tr><tr><td>6503</td><td>5.27</td><td>0.32</td><td>1.9</td><td>40</td><td>3</td><td>0.06 ± 0.20</td><td>37</td><td>3</td><td>-0.19 ± 0.05</td><td>59</td><td>3</td><td>0.04 ± 0.09</td><td>96</td><td>3</td><td>0.02 ± 0.13</td></tr><tr><td>7793</td><td>3.44</td><td>0.52</td><td>3.2</td><td>41</td><td>3</td><td>-0.18 ± 0.13</td><td>116</td><td>3</td><td>-0.11 ± 0.44</td><td>85</td><td>4</td><td>0.02 ± 0.30</td><td>201</td><td>4</td><td>0.09 ± 0.24</td></tr><tr><td>5291<sup>e</sup></td><td>63.5</td><td>0.15</td><td>2.00</td><td>106</td><td>6</td><td>-0.65</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SFR:<sup>f</sup></td><td><i>S</i> ± <i>ε</i></td><td>χ<sup>2</sup></td><td>DOF</td><td>0.14 ± 0.35</td><td>37</td><td>9</td><td>0.28 ± 0.11</td><td>7.2</td><td>8</td><td>0.22 ± 0.15</td><td>3.3</td><td>9</td><td>0.14 ± 0.13</td><td>6.5</td><td>12</td></tr><tr><td>Mass:</td><td><i>S</i> ± <i>ε</i></td><td>χ<sup>2</sup></td><td>DOF</td><td>0.10 ± 0.26</td><td>37</td><td>9</td><td>0.25 ± 0.07</td><td>3.3</td><td>8</td><td>0.17 ± 0.14</td><td>3.6</td><td>9</td><td>0.14 ± 0.10</td><td>4.9</td><td>12</td></tr><tr><td>Area:</td><td><i>S</i> ± <i>ε</i></td><td>χ<sup>2</sup></td><td>DOF</td><td>-0.10 ± 0.38</td><td>79</td><td>9</td><td>0.25 ± 0.14</td><td>3.8</td><td>8</td><td>0.24 ± 0.15</td><td>0.9</td><td>9</td><td>0.20 ± 0.10</td><td>4.1</td><td>12</td></tr><tr><td>SFR/A:</td><td><i>S</i> ± <i>ε</i></td><td>χ<sup>2</sup></td><td>DOF</td><td>0.45 ± 0.42</td><td>34</td><td>9</td><td>0.12 ± 0.39</td><td>15</td><td>8</td><td>-0.02 ± 0.42</td><td>3.0</td><td>9</td><td>-0.18 ± 0.22</td><td>14</td><td>12</td></tr><tr><td>SFR:<sup>g</sup></td><td><i>S</i> ± <i>ε</i></td><td>χ<sup>2</sup></td><td>DOF</td><td>0.06 ± 0.39</td><td>14</td><td>9</td><td>0.10 ± 0.08</td><td>2.7</td><td>8</td><td>0.09 ± 0.10</td><td>1.6</td><td>9</td><td>0.08 ± 0.07</td><td>5.4</td><td>12</td></tr><tr><td>Mass:</td><td><i>S</i> ± <i>ε</i></td><td>χ<sup>2</sup></td><td>DOF</td><td>0.05 ± 0.29</td><td>14</td><td>9</td><td>0.10 ± 0.06</td><td>1.9</td><td>8</td><td>0.06 ± 0.09</td><td>1.6</td><td>9</td><td>0.07 ± 0.06</td><td>4.3</td><td>12</td></tr></tbody></table>

**Notes.**

a NGC numbers except for 1249, which is a UGC number. NGC 5194 has a first row for clusters identified by eye and a second row ("ML") for cluster identified by machine learning.

b Distance in megaparsecs, star formation rate (SFR) in $M_{\odot}$ yr$^{-1}$, Galaxy stellar mass $M_{\star}$ in $10^9 M_{\odot}$.

c Cluster classes: 1 is compact, class 2 is elongated, and class 3 is multi-core.

d $N$ represents the number of clusters in this class with a mass exceeding a certain distance-dependent limit and age less than 126 Myr (log age = 8.1); $B$ represents the number of mass bins for the clusters, each of width 0.5 dex; $S$ is the slope of the RAMPS, i.e., the derivative of the log of the RAMPS with respect to log of the cluster mass for equal intervals of the log of the cluster mass.

e NGC 5291 is the system of tidal dwarfs, which have a total of 106 usable clusters that have not been classified according to the LEGUS system.

f These four rows give the slopes and errors ($S \pm \epsilon$) of the RAMPS slopes versus log(SFR), log($M_{\star}$), log(Area) (in kpc$^2$), and log(SFR/Area) (in $M_{\odot}$ kpc$^{-2}$ yr$^{-1}$), where area is measured at the radius of 25 mag per square arcsecond in the $B$ band. These rows also show the $\chi^2$ values and the number of degrees of freedom, DOF, which is the number of RAMPS values minus 2.

g These next two rows are for RAMPS with equal numbers of clusters per interval of the cluster mass.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (Top) The RAMPS are plotted vs. the cluster mass *M* (in $M_\odot$) for LEGUS clusters in classes 2 and 3 (blue lines) and for NGC 5291 (red line) clusters. The RAMPS are dimensionless. Error bars are from the student-t distribution at 90% probability; they are symmetric around the linear value of the RAMPS but asymmetric here when plotted logarithmically (points with no error bars have only 2 clusters in that mass interval). Cluster masses are randomly offset from the center of the mass bin for clarity; different galaxies have different symbols. Increasing RAMPS are on the left and decreasing RAMPS are on the right. (middle, bottom) The slopes of the RAMPS are shown vs. the total galaxy star formation rates (in $M_\odot$ yr$^{-1}$) and galaxy stellar masses (in $M_\odot$). The middle panels are for RAMPS calculated using equal log(*M*) intervals, as in the top panels, while the bottom panels are for RAMPS calculated using equal numbers of clusters in each mass interval. Low-mass galaxies tend to have decreasing RAMPS, indicating a greater tendency for high-mass clusters to collect together compared to low-mass clusters. The green line and crosses on the top left and green crosses in the middle and bottom panels are for NGC 5194 clusters determined by machine learning; the blue crosses on the top left are for NGC 5194 visual identifications.</figcaption></figure>

Table 1 also lists the 11 galaxies that have 35 or more class 1 clusters spanning a factor of $10^{1.5}$ or more in mass within the same age and mass limits. Similarly, the table lists the cluster counts and RAMPS slopes considering only class 2 types alone and only class

3 types alone. Class 1 clusters are compact, class 2 are somewhat elongated, and class 3 are multi-core (Adamo et al. 2017). The galaxies span a factor of ~100 in SFR and stellar mass, so they represent a fair sample of spiral and dwarf galaxy types.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (Left) Galaxy stellar mass in $M_{\odot}$ is shown as a function of galaxy distance in Mpc. (Right) The slope of the RAMPS for class 2 and 3 clusters is shown as a function of distance. The blue points are LEGUS galaxies and the red point is for tidal dwarfs around NGC 5291. For the LEGUS galaxies, there is a correlation between galaxy mass and distance, so a resulting correlation appears between RAMPS slope and distance, but the RAMPS does not appear to be biased for distance by itself.</figcaption></figure>

The clusters were first divided into logarithmic mass intervals of 0.5 dex; the number of bins in Table 1 represents the number of these intervals spanned by the cluster masses. As for the random trial discussed above, the projected separations between each cluster and all the other clusters in the same mass interval (corrected for galaxy inclination) were determined and the minimum of these separations was noted. This minimum represents the distance between each particular cluster and its nearest cluster of the same mass. The average of these minimum distances was then determined for each mass interval. After multiplication by the square root of the number of clusters in the mass interval and division by the average separation in the lowest-mass interval, we obtain the RAMPS.

Figure 2 shows the RAMPS for the LEGUS clusters in classes 2 and 3, divided into those with increasing or nearly constant average slopes on the left and those with decreasing average slopes on the right. These slopes are determined from the whole mass range, which, e.g., for 4 bins, corresponds to a factor of ~100.

The middle panels in Figure 2 show the RAMPS slope versus the SFR (left) and the galaxy stellar mass $M_{\star}$ (right). There is an increasing trend in both panels. The slopes of these trends were determined by least-squares linear fits with uncertainties given by the student-t distribution at 90% probability. Corresponding $\chi^2$ values are from the sum of the squared differences between the RAMPS slopes and the linear fits, normalized to the uncertainties in the slopes. The result is a slope of $0.14 \pm 0.13$, with $\chi^2 = 6.5$ for 12 degrees of freedom (DOF) in the plot of RAMPS slope versus log(SFR), and $0.14 \pm 0.10$ with $\chi^2 = 4.9$ for 12 DOF in the plot of RAMPS slope versus log $M_{\star}$. These $\chi^2$ values are smaller than the number of DOF, which is number of RAMPS values minus 2, indicating reasonably good fits. A summary is in Table 1.

The bottom panels in Figure 2 show the RAMPS slopes calculated with an equal number of clusters in each of 6 cluster mass intervals (rather than equal log($M$) intervals). The upward trends are present for this binning too, although slightly

smaller. Still, their slopes (Table 1) are within the error bars of the slopes in the middle panels.

The correlations between the slope of the RAMPS and the galaxy SFR or stellar mass imply that low-mass galaxies have their most massive class 2–3 clusters closer together than average, whereas high-mass galaxies have all cluster masses randomly distributed.

Inclination errors could affect the results, but a recalculation of the class 2+3 case with zero inclinations gave about the same slopes: $0.16 \pm 0.14$ versus logSFR and $0.16 \pm 0.11$ versus log($M_{\star}$). The $\chi^2$ values were much higher without inclination corrections, however, 31 and 18, respectively.

The clusters were identified by eye for all galaxies, but for NGC 5194, clusters were also identified with machine learning (ML) techniques, using the visual identifications as a training set (Messa et al. 2018; Grasha et al. 2019). There are more ML clusters than visual clusters for this galaxy, but the slope of the RAMPS is about the same in both cases. The green points in the bottom of Figure 2 and the green crosses and line in the top left are for clusters identified by ML in NGC 5194, compared to the blue points at the same SFR and stellar mass on the bottom and the blue crosses on the top left.

We also considered clusters in the tidal dwarf galaxies connected with the interacting galaxy NGC 5291. These were obtained from the study by Fensch et al. (2019) with distances between the clusters determined by assuming zero inclination and position angle, as these are highly irregular galaxies. Among their sample of 272 clusters with masses less than $3 \times 10^5 M_{\odot}$ and not category 0 (which are excluded because they have more than two HST passbands with only upper limits on the flux), we include all 106 clusters with masses larger than $10^4 M_{\odot}$ and ages less than 100 Myr. These limits avoid the loss of clusters from fading. There are several dwarfs in the collision debris but we can treat all of them as one large distribution to derive the slope of the RAMPS because the nearest cluster to any given cluster of the same mass is likely to be inside the same tidal dwarf. The other factors in RAMPS, $N(M)$ and $D_{\text{low}}$, are constant and do not affect the slope. The
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (Top) The RAMPS are plotted vs. the cluster mass ($M_{\odot}$) for LEGUS clusters in class 1. (Middle, bottom) The slopes of the RAMPS are shown vs. the total galaxy star formation rates (in $M_{\odot}$ yr$^{-1}$) and galaxy stellar masses (in $M_{\odot}$). Unlike the class 2 and 3 clusters in Figure 2, the class 1 clusters appear to be randomly distributed for all galaxy masses. The green line and crosses on the top right and the green crosses in the middle and bottom panels are for NGC 5194 clusters discovered by machine learning. As in Figure 2, the middle panels are for RAMPS calculated with equal log($M$) intervals of cluster mass, and the bottom panels are for RAMPS calculated with equal numbers of clusters in each mass interval.</figcaption></figure>

red segmented line and red points in Figure 2 show the results for NGC 5291. The RAMPS for NGC 5291 has a negative slope like the other dwarf galaxies.

Figure 3 shows $M_{\star\text{star}}$ versus distance, with the expected trend for lower-mass galaxies, which are more common, to be closer. Because of this, the closer galaxies also have more negative

slopes in the RAMPS, as shown in the right panel. These distance trends are not the cause of the varying slopes for RAMPS, however. The cluster separations are well resolved for all of the galaxy distances, and the cluster mass ranges are about the same scaled for distance. Moreover, the distance to NGC 5291 (red points in Figure 3) is larger than even the
massive galaxies in LEGUS, yet the slope of its RAMPS is negative, like the closer dwarfs.

The RAMPS for class 1 clusters in LEGUS is shown in Figure 4. Again the positive and negative slopes are separated in the two top panels, but now the slopes are all around 0, as also shown in the middle and bottom panels. Evidently the class 1 clusters have different grouping properties than the class 2 and 3 clusters. Separate plots for class 2 clusters and for class 3 clusters alone (not shown) repeat the correlation in Figure 2, which was for the combined classes. The numbers of clusters and RAMPS slopes for these separate classes are given in Table 1 for completeness.

We also determined the RAMPS for only the young class 1 clusters with ages less than 30 Myr (not shown); these had no obvious trends with galaxy stellar mass either, although there were only six galaxies with enough young clusters to plot. The slope of the RAMPS slope versus log(SFR) linear fit for young class 1 clusters is $0.11 \pm 0.28$ ($\chi^2 = 0.8$, DOF = 4), but the range in SFR is only a factor of 7. Versus log($M_{\text{star}}$), the slope is $0.05 \pm 0.18$ ($\chi^2 = 0.9$, DOF = 4) with a factor of 28 range in $M_{\text{star}}$.

We also checked whether class 1 and class 2+3 clusters show a significant correlation between RAMPS and the SFR per unit area. The areas were determined from the distances and from $D_{25}$ in de Vaucouleurs et al. (1991) or the NASA/IPAC Extragalactic Database. No trends were found. There is a correlation with area alone, however, similar to that in Figures 2 and 4, such that more massive clusters are more clumped together in smaller galaxies, which are also the lower-mass galaxies in the previous figures.

Table 1 summarizes for all cluster classes the slopes of the linear fits, $S \pm \epsilon$, the $\chi^2$ values, and the number of DOF for the RAMPS slopes versus log(SFR), log($M_{\text{star}}$), log(Area), and log(SFR/Area).

## 4. Discussion

The RAMPS method suggests that class 2 and 3 star clusters in LEGUS, which are somewhat elongated or irregular in shape, are mass-segregated in low-mass galaxies, whereas class 1 clusters, which are compact, are not. Clusters in the tidal dwarfs around NGC 5291 are also mass-segregated.

Dwarf galaxies differ from high-mass galaxies in ways that could account for the trends. For example, dwarfs have weaker gravitational potentials than massive galaxies, so the mutual attraction between massive clouds and the clusters they form is larger in proportion to background tidal forces for a dwarf galaxy. The tidal dwarfs around NGC 5291 could also be devoid of dark matter, which makes their background gravitational potential even weaker. This implies that massive clouds and clusters can move closer to each other, or accrete more interstellar gas, with less of an influence from galactic shear in lower-mass galaxies or tidal dwarfs. This explanation does not obviously account for the lack of mass segregation by class 1 clusters, however.

Alternatively, high-mass clusters could scatter away their low-mass neighbors more effectively when Coriolis forces are low, leaving the high-mass clusters more concentrated in each star-forming region than average. This explanation could include the observed difference between cluster classes because class 1 clusters might scatter better than class 2–3 clusters, which could break apart during the process. Class 1 clusters in LEGUS are also older on average than class 2 and 3 clusters

(Grasha et al. 2017a), giving them more time to scatter. However, the young class 1 clusters, less than 30 Myr, did not show a correlation with galaxy stellar mass. Class 1 clusters could be a random selection of clusters that are dense.

Other models for cluster mass segregation in small galaxies could be developed around a possibly larger Jeans mass for fragments near the center of a star-forming gas complex, or a higher gas density in higher-mass clouds. Numerical simulations need to address these possibilities.

The authors are grateful to the referee for comments on the manuscript. This work is based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with programs #13364 (LEGUS) and #14727 (NGC 5291), including grant HST-GO-14727.004-A to B.G.E. A.A. acknowledges the support of the Swedish Research Council, Vetenskapsrådet, and the Swedish National Space Agency (SNSA). This research has made use of the NASA/IPAC Extragalactic Database (NED), which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology.

## ORCID iDs

Bruce G. Elmegreen https://orcid.org/0000-0002-1723-6330

A. Adamo https://orcid.org/0000-0002-8192-8091

F. Bournaud https://orcid.org/0000-0002-5743-0250

D. Calzetti https://orcid.org/0000-0002-5189-8004

D. O. Cook https://orcid.org/0000-0002-6877-7655

D. A. Dale https://orcid.org/0000-0002-5782-9093

D. M. Elmegreen https://orcid.org/0000-0002-1392-3520

J. Fensch https://orcid.org/0000-0002-4383-0202

K. Grasha https://orcid.org/0000-0002-3247-5321

Hwi Kim https://orcid.org/0000-0003-4770-688X

L. Kahre https://orcid.org/0000-0002-2037-3096

E. Sabbi https://orcid.org/0000-0003-2954-7643

L. J. Smith https://orcid.org/0000-0002-0806-168X

## References

Adamo, A., Ryon, J. E., Messa, M., et al. 2017, ApJ, 841, 131

Alfaro, E. J., & Román-Zúñiga, C. G. 2018, MNRAS, 478, L110

Bastian, N., Gieles, M., Ercolano, B., & Guttermuth, R. 2009, MNRAS, 392, 868

Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 1997, MNRAS, 285, 201

Calzetti, D., Johnson, K. E., Adamo, A., et al. 2015, ApJ, 811, 75

Cartwright, A., & Whitworth, A. P. 2004, MNRAS, 348, 589

de la Fuente Marcos, R., & de la Fuente Marcos, C. 2009, ApJ, 700, 436

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H., et al. 1991, Third Reference Catalogue of Bright Galaxies (New York: Springer)

Efremov, Y. N., & Elmegreen, B. G. 1998, MNRAS, 299, 588

Elmegreen, B. G. 2008, ApJ, 672, 1006

Elmegreen, B. G. 2018, ApJ, 853, 88

Elmegreen, B. G., & Efremov, Y. N. 1997, ApJ, 480, 235

Elmegreen, B. G., Elmegreen, D. M., Chandar, R., Whitmore, B., & Regan, M. 2006, ApJ, 644, 879

Elmegreen, B. G., & Falgarone, E. 1996, ApJ, 471, 816

Elmegreen, B. G., Hurst, R., & Koenig, X. 2014, ApJL, 782, L1

Feitzinger, J. V., & Galinski, T. 1987, A&A, 179, 249

Fensch, J., Duc, P.-A., Boquien, M., et al. 2019, A&A, 628, A60

Fleck, R. C., Jr. 1996, ApJ, 458, 739

Grasha, K., Calzetti, D., Adamo, A., et al. 2017a, ApJ, 840, 113
Grasha, K., Calzetti, D., Adamo, A., et al. 2019, MNRAS, 483, 4707

Grasha, K., Elmegreen, B. G., Calzetti, D., et al. 2017b, ApJ, 842, 25

Hunter, D. A. 1999, in Wolf-Rayet Phenomena in Massive Stars and Starburst Galaxies, ed. K. A. van der Hucht, G. Koenigsberger, & P. R. J. Eenens (San Francisco, CA: ASP), 418

Hunter, D. A., Ficut-Vicas, D., Ashley, T., et al. 2012, AJ, 144, 134

Krumholz, M. R., McKee, C. F., & Bland-Hawthorn, J. 2018, ARA&A, 57, 227

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57

Lahén, N., Naab, T., Johansson, P. H., et al. 2019, ApJL, 879, L18

Maíz-Apellániz, J. 2001, ApJ, 563, 151

Messa, M., Adamo, A., Östilin, G., et al. 2018, MNRAS, 473, 996

Portegies Zwart, S. F., McMillan, S. L. W., & Gieles, M. 2010, ARA&A, 48, 431

Scalo, J. S. 1985, in Protostars and Planets II, ed. D. C. Black & M. S. Matthews (Tucson, AZ: Univ. Arizona Press), 201

Zinnecker, H. 1982, in Symposium on the Orion Nebula to Honor Henry Draper, ed. A. E. Glassgold et al. (New York: New York Academy of Sciences), 226