# Characterising open clusters in the solar neighbourhood with the *Tycho-Gaia* Astrometric Solution*

T. Cantat-Gaudin<fnref n="1" />, A. Vallenari<fnref n="1" />, R. Sordo<fnref n="1" />, F. Pensabene<fnref n="1" />,<fnref n="2" />, A. Krone-Martins<fnref n="3" />, A. Moitinho<fnref n="3" />, C. Jordi<fnref n="4" />, L. Casamiquela<fnref n="4" />, L. Balaguer-Núnez<fnref n="4" />, C. Soubiran<fnref n="5" />, and N. Brouillet<fnref n="5" />

Received 26 May 2017 / Accepted 29 January 2018

## ABSTRACT

**Context.** The *Tycho-Gaia* Astrometric Solution (TGAS) subset of the first *Gaia* catalogue contains an unprecedented sample of proper motions and parallaxes for two million stars brighter than *G* ~ 12 mag.

**Aims.** We take advantage of the full astrometric solution available for those stars to identify the members of known open clusters and compute mean cluster parameters using either TGAS or the fourth U.S. Naval Observatory CCD Astrograph Catalog (UCAC4) proper motions, and TGAS parallaxes.

**Methods.** We apply an unsupervised membership assignment procedure to select high probability cluster members, we use a Bayesian/Markov Chain Monte Carlo technique to fit stellar isochrones to the observed 2MASS JHKs magnitudes of the member stars and derive cluster parameters (age, metallicity, extinction, distance modulus), and we combine TGAS data with spectroscopic radial velocities to compute full Galactic orbits.

**Results.** We obtain mean astrometric parameters (proper motions and parallaxes) for 128 clusters closer than about 2 kpc, and cluster parameters from isochrone fitting for 26 of them located within a distance of 1 kpc from the Sun. We show the orbital parameters obtained from integrating 36 orbits in a Galactic potential.

**Key words.** galaxies: star clusters: general – methods: numerical

## 1. Introduction

Open clusters (OCs) are basic constituents of our Galaxy. It is believed that most stars, if not all, are born within clusters, which are later disrupted and their population mixes with what is referred to as the “field population” (see e.g. Lada & Lada 2003). The homogeneous age and chemical composition of stars within a cluster makes them valuable tracers of the properties of the Galactic disk (Cantat-Gaudin et al. 2014; Vázquez et al. 2010), as well as testbeds for stellar evolution models. OCs trace the radial metallicity gradient (Cantat-Gaudin et al. 2016; Jacobson et al. 2016); their age, metallicity, and kinematics are important to define how the radial migration process can affect the disk (Minchev 2016).

Although different authors may use slightly different definitions for what constitutes an OC, it is generally considered that stars within a cluster have a common origin (they have similar ages and chemical compositions) and are bound by gravity (they are physically close to each other and share a common orbit through the disk of the Milky Way). Historically, many stellar clusters have been serendipitously discovered as local over-densities of stars (e.g. Messier 1781).

The majority of known OCs are located within 2 kpc of the Sun. The most recent version of the OC catalogue of Dias et al. (2002, hereafter DAML) quotes about 2200 objects, while the Milky Way Star Clusters catalogue of Kharchenko et al. (2013) lists over 3000 objects, many of which are putative or candidate clusters. On one hand, the cluster census is far from being complete, especially at the faint end where small and sparse objects and remnants of disrupted clusters can escape detection (Bica & Bonatto 2011) (even when their stars are visible) because they do not stand as significant over-densities. On the other hand, cluster searches produce false positives and astrometric information is necessary to distinguish a genuine cluster (sharing common proper motions and parallaxes) from a coincidental asterism.

HIPPARCOS (Perryman et al. 1997) was the first space mission dedicated to astrometry. The data it collected yielded parallaxes for about 120 000 stars in our Galaxy (ESA 1997; van Leeuwen 2007) and allowed for the construction of the 2.5 million stars proper motion catalogue *Tycho-2* (Høg et al. 2000). The ongoing *Gaia* space mission (Perryman et al. 2001; Gaia Collaboration 2016b) will deliver proper motions and parallaxes for a billion sources, dwarfing the HIPPARCOS catalogue by four orders of magnitude. *Gaia* is expected to revolutionize the field, discovering new objects (Koposov et al. 2017) and providing accurate astrometrically based characterizations of thousands of clusters. The *Gaia* mission first data release (hereafter GDR1 Gaia Collaboration 2016a) dataset, which contains

<footnote n="1">¹ INAF-Osservatorio Astronomico di Padova, vicolo Osservatorio 5, 35122 Padova, Italy
e-mail: tristan.cantat@oapd.inaf.it</footnote>

<footnote n="2">² Dipartimento di Fisica e Astronomia, Università di Padova, vicolo Osservatorio 3, 35122 Padova, Italy</footnote>

<footnote n="3">³ SIM, Faculdade de Ciências, Universidade de Lisboa, Ed. C8, Campo Grande, 1749-016 Lisboa, Portugal</footnote>

<footnote n="4">⁴ Institut de Ciències del Cosmos, Universitat de Barcelona (IEEC-UB), Martí i Franquès 1, 08028 Barcelona, Spain</footnote>

<footnote n="5">⁵ Laboratoire d'Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, UMR 5804, 33615 Pessac, France</footnote>

<footnote>* Cluster data are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/615/A49</footnote>
positions and G-band magnitudes for a billion sources, also contains a subset called the *Tycho-Gaia* Astrometric Solution (TGAS), consisting of two million sources with a full astrometric solution obtained through the combination of *Gaia* observations and the positions listed in *HIPPARCOS* and *Tycho-2* (Lindegren et al. 2016). The next *Gaia* data releases are expected to deliver astrometric measurements for more than a billion stars. Handling such large, multi-dimensional datasets (whose applications extend well beyond stellar cluster science) requires automated methods in order to identify and select cluster stars and to characterise stellar clusters (see e.g. Robichon et al. 1999, with *HIPPARCOS* data). In this paper we chose to follow the unsupervised photometric membership assignment in stellar clusters approach (UPMASK) introduced by Krone-Martins & Moitinho (2014), applying it to astrometric data for the first time.

The aims of this paper are two-fold: (i) to validate the use of tools that can be applied to multi-dimensional datasets such as the two million sources in TGAS, or the billion sources of the upcoming further *Gaia* releases, and (ii) to update the cluster census in the solar neighbourhood, deriving memberships, and mean parallaxes and proper motions from the stellar clusters that can be clearly identified in the TGAS catalogue.

In Sect. 2 we present the method we applied to select probable cluster members in the TGAS catalogue for 128 OCs, extending the use of the code UPMASK (Krone-Martins & Moitinho 2014) to astrometric data. Section 3 presents the mean astrometric parameters we computed for these clusters, and Sect. 4 presents the ages and cluster parameters we obtained for 34 of them by isochrone fitting using the Bayesian approach of the code (BASE-9 von Hippel et al. 2006). Finally, Sect. 6 discusses the results obtained in the context of the Milky Way disk, and Sect. 7 closes with concluding remarks.

## 2. Membership from TGAS data

In this section, we compare the uncertainties in the proper motions in the TGAS and UCAC4 datasets, and discuss the membership of the clusters present in TGAS.

### 2.1. Proper motions in the TGAS dataset

In this paper the parallaxes used were always the TGAS parallaxes, but the proper motions were either taken from TGAS or from UCAC4 (Zacharias et al. 2012). Because the *Gaia* scanning law in the first 14 months that are included in GDR1 has been favouring more frequent visits to regions closer to the ecliptic poles (*Gaia* Collaboration 2016a), the quality of the TGAS data varies significantly between regions across the sky. In some areas the uncertainties on the proper motions are larger than those of the UCAC4 catalogue. The uncertainties on proper motions of TGAS stars are in the range 0.5–2.6 mas yr⁻¹, with a median value of 1.1 mas yr⁻¹ (Lindegren et al. 2016). UCAC4 proper motion formal uncertainties are in the range 1–10 mas yr⁻¹, with possible systematics of the order of 1–4 mas yr⁻¹ (Zacharias et al. 2013), on average higher than TGAS uncertainties. However, the effect of the *Gaia* scanning law is that in some regions TGAS uncertainties can be on average higher than in UCAC4. This effect is clearly visible in Fig. 1, which shows a comparison between the proper motion uncertainties in UCAC4 and TGAS for the stars present in both catalogues. Figure 2 shows the mean proper motion difference between TGAS and UCAC4 for those stars. Local systematic differences of the order of 1 to 3 mas yr⁻¹ are present all across the sky.

<figure><img src="image_1.png" /><figcaption>**Fig. 1.** Median proper motion error ratio between UCAC4 and TGAS for stars present in both catalogues, in HEALPix level five pixels, in Galactic coordinates (north to the top, increasing longitude to the left, Galactic centre at the centre).</figcaption></figure>

In each field under analysis, in order to see the cluster members as a more compact group in astrometric space, we used the proper motion catalogue, which provided the smallest median uncertainty, thus using TGAS proper motions in the regions where the ratio shown in Fig. 1 is above one (e.g. for NGC 2360, for which the TGAS and UCAC4 proper motions are shown in Fig. 7). In practice, whether UCAC4 or TGAS have better proper motion errors for individual stars does not depend solely on the region, but also on the type of astrometric prior used in TGAS. Even in the regions where UCAC4 provides better median errors than TGAS, the subset of TGAS stars whose proper motion was determined using their *HIPPARCOS* positions as astrometric prior always has smaller errors than UCAC4 (typically under 0.1 mas). Those stars only account for 5% of the whole TGAS dataset and for less than 0.4% of stars fainter than $G = 10$. We therefore based our decision to use one proper motion catalogue or the other based on the cluster position only.

### 2.2. Target selection

We obtained a list of cluster coordinates and parameters from the Milky Way Star Clusters catalogue of Kharchenko et al. (2013, hereafter MWSC). Since in the TGAS data all stars with parallaxes under 0.5 mas have relative errors $\sigma_\omega/\omega$ larger than 50%, we limited our study to OCs with expected parallaxes larger than 0.5 mas (closer than 2000 pc). Since the TGAS stars are all contained in the *Tycho-2* astrometric catalogue, it is affected by the same completeness limit of $V \sim 11.5$. Making use of PARSEC stellar isochrones (Bressan et al. 2012), we rejected all the clusters which, according to their listed age, distance, and extinction, were not expected to contain any star brighter than $V = 12$.

Finally, excluded from this study the 19 nearby OCs studied in *Gaia Collaboration* (2017), namely: the Hyades, the Pleiades, Coma Berenices, Praesepe, Alpha Perseus, Blanco 1, Collinder 140, IC 2391, IC 2602, IC 4665, NGC 2451A, NGC 6475, NGC 6633, NGC 7092, NGC 2516, NGC 2232, NGC 2422, NGC 3532, and NGC 2547.

This left us with a list of 694 OCs to investigate. For these objects, the MWSC quotes tidal radii of up to 15 pc. To ensure we did not cut out potential cluster stars by applying selections
<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Top: mean $\mu_{\alpha}^*$ difference between UCAC4 and TGAS for stars present in both catalogues, in HEALPix level five pixels, in Galactic coordinates (north to the top, increasing longitude to the left, Galactic centre at the centre). Bottom: as for the $\mu_{\delta}$ component of the proper motion.</figcaption></figure>

that were too narrow, we kept all TGAS stars within an angle corresponding to a physical distance of 20 pc at the reference distance of the cluster, and visually examined the sky distribution, parallax and proper motion distribution, and colour-magnitude for the members and field stars for each cluster in order to confirm the reality of the object. When the procedure failed to identify any cluster member, we restricted the field of view to an angle corresponding to 10 pc in order to provide a better contrast between field and cluster stars in the astrometric space. The final search radius used for each cluster is listed in Table A.1.

## 2.3. Membership determination

Our determination of cluster membership relies on proper motions and parallaxes. To eliminate obvious field stars we first performed a broad selection, rejecting all stars with proper motions more distant than 10 mas yr⁻¹ from the literature value. For clusters closer than 500 pc (expected parallax of 2 mas) we also rejected stars with parallaxes under 1 mas. About 90% of the TGAS sources have parallax errors smaller than 0.5 mas. From the parallax distribution error (shown in Fig. 3) we expect that only 1% of TGAS sources have observed parallaxes that are lower than their true value by 1 mas. This means that even for clusters with a true parallax of exactly 2 mas, applying a cut-off at 1 mas introduces no significant bias in the remaining parallax distribution. Although the field of view of some pairs of OCs overlap (e.g. NGC 2451 with NGC 2477), the difference in proper motions allowed us to tell the two clusters apart.

<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Parallax error distribution (in bins of width 0.01 mas) for the stars used in this study.</figcaption></figure>

For the remaining star sample, our approach to membership determination is based on the principles of the unsupervised membership assignment method UPMASK (Krone-Martins & Moitinho 2014). This approach does not rely on strong physical assumptions concerning the nature of a cluster (no assumptions on density profile modelling or on the structure in photometric space), except that its stars share common properties and that the spatial distributions of cluster stars is concentrated, while the distribution of field stars is random. In machine learning and data mining it is common to refer to a group of objects as a cluster. In this paper we tried to refer to the output of *k*-means clustering (abstract grouping of datapoints with similar observed properties) as groups, to avoid any possible confusion with stellar clusters (astronomical objects). The core idea of UPMASK is to apply a simple clustering algorithm (for instance, *k*-means clustering) to identify small groups of stars with similar colours and magnitudes in various photometric bands, then check all these small groups individually and determine whether their spatial distribution is more tightly concentrated than a random distribution (this is referred to as the “veto” step).

Although UPMASK was originally built to identify stellar clusters based on photometry (and is used for this, for instance in Costa et al. 2015), its core principle can easily be generalised to other types of quantities, as its only strong assumption is to consider that cluster members must be, in any observable space, more tightly distributed than field stars. This assumption happens to hold even truer for astrometry than photometry, as all stars within a cluster are expected to be located at the same distance from us and moving in the same direction, regardless of their colour and luminosity. In this study, rather than applying the *k*-means clustering to a set of magnitudes in different photometric filters, we applied it to the three-dimensional astrometric space of proper motions and parallaxes ($\mu_{\alpha}^*, \mu_{\delta}, \omega$). As recommended in Krone-Martins & Moitinho (2014), we scale each of these three observables to unit variance. The *k*-means clustering method does not allow the user to impose the number of points in each group, as they can be of variable sizes, but it requires the user to choose the number of desired groups (which is equivalent to setting the mean number of stars in each group). Krone-Martins & Moitinho (2014) report that they obtain the best results using values between ten and 25. In this study, since we expect several clusters to have very few members in the TGAS data, we set this mean number to ten.

To determine whether or not each identified group is spatially more concentrated than a random distribution, we applied
<figure><img src="image_4.png" /><figcaption>Fig. 4. Top left panel: proper motions of stars in field of NGC 752. The colour code corresponds to the groups identified by k-means clustering in the (μα*, μδ, ϖ) space. The sky distribution of the three highlighted groups is shown in the other panels. The quantity Λ is a measurement of spatial clustering, as defined in Eq. (1). Stars are considered potential cluster members when Λ > 1.</figcaption></figure>

the method introduced by Allison et al. (2009) originally used to reveal mass segregation in star clusters (by showing whether high-mass stars are spatially more concentrated than the average of the cluster). The method consists in comparing the total branch length $l_{obs}$ of the minimum spanning tree connecting all stars in that group to the expected value $l$ in a random distribution containing the same number of stars.

To save computation time, for sample sizes from three to 80 we pre-computed and tabulated the expected $l$ and associated standard deviations $\sigma_l$ by generating 2000 random circular distributions. If the $\langle l_{obs} \rangle$ value in a given group is smaller than the tabulated $l$ by at least $\sigma_l$, that is if

$$ \Lambda = \frac{l - l_{obs}}{\sigma_l} > 1, \qquad (1) $$

then all the stars in that group are flagged as possible cluster members. Figure 4 shows an example of k-means clustering applied to the stellar cluster NGC 752 and the sky distribution of three selected groups. In that example, two of the shown groups are considered as containing potential cluster members.

When this veto step has been applied to all stars, the procedure (k-means clustering and spatial veto) is performed again, but instead of using the catalogue value we add to each datapoint a random offset in proper motion and parallax, corresponding to the uncertainties. After a total of 100 iterations, the frequency with which a star was flagged as part of a clustered group is interpreted as its membership probability.

The random nature of the grouping step performed in the heuristic of UPMASK means that small groups of field stars might sometimes be flagged as clustered (in other words, even purely random distribution are expected to satisfy Eq. (1)) and these stars end up with non-zero membership probabilities of a few percent, which can be considered noise level. To obtain cleaner results and a better contrast between field and cluster

<figure><img src="image_5.png" /><figcaption>Fig. 5. Top left: positions of the stars identified as probable cluster members (dots coloured according to membership probability) and field stars (black points). Top right: proper motions of the probable cluster members and of the field stars. Bottom left: JKs photometry for the cluster stars and field stars. Bottom right: parallax ϖ against J-magnitude for cluster and field stars.</figcaption></figure>

stars, we applied the procedure a second time for all OCs, after discarding stars for which the first run yielded probabilities lower than 10%. We consider the final membership probability to be the result of this second run. The individual probabilities of all stars in the investigated fields of 128 clusters (including those with low or zero-probability) are provided as an electronic table at the CDS. In total, we find 4450 potential cluster members (probability > 50%) and 851 secure members (>90%). The result of the member selection procedure is shown in Fig. 5 for the cluster Alessi 3 using TGAS proper motions and parallaxes.

## 2.4. Taking correlations into accounts

Whenever possible (i.e. when the TGAS proper motions were used), the errors we added to the original values in each random drawing took into account the full covariance between the error in three parameters (μα*, μδ, ϖ). For each star, the covariance matrix is

$$ \text{Cov} = \begin{bmatrix} \sigma_{\mu_{\alpha} *}^2 & \sigma_{\mu_{\alpha} * \sigma_{\mu_{\delta}} \rho_{\mu_{\alpha} * \mu_{\delta}}} & \sigma_{\mu_{\alpha} * \sigma_{\omega} \rho_{\mu_{\alpha} * \omega}} \\ \sigma_{\mu_{\alpha} * \sigma_{\mu_{\delta}} \rho_{\mu_{\alpha} * \mu_{\delta}}} & \sigma_{\mu_{\delta}}^2 & \sigma_{\mu_{\delta} \sigma_{\omega} \rho_{\mu_{\delta} \omega}} \\ \sigma_{\mu_{\alpha} * \sigma_{\omega} \rho_{\mu_{\alpha} * \omega}} & \sigma_{\mu_{\delta} \sigma_{\omega} \rho_{\mu_{\delta} \omega}} & \sigma_{\omega}^2 \end{bmatrix}, \quad (2) $$

where $\sigma_{\mu_{\alpha} *}$, $\sigma_{\mu_{\delta}}$, and $\sigma_{\omega}$, and $\rho_{\mu_{\alpha} * \mu_{\delta}}$, $\rho_{\mu_{\alpha} * \omega}$, and $\rho_{\mu_{\delta} \omega}$ (the correlation coefficients) are all provided in the TGAS catalogue. Neglecting the off-diagonal terms is equivalent to considering the correlation coefficients $\rho$ as equal to zero (i.e. that the errors are independent), which is generally not true. This is illustrated in Fig. 6, where $\rho$ for proper motions and parallaxes are shown to be significantly different from zero for a large number of stars used in this study. Indeed, the correlations between the errors of astrometric parameters are expected to vary in different regions of the sky, as explained in Lindegren et al. (2016, see for instance the map in their Fig. 7).
<figure><img src="image_6.png" /><figcaption>**Fig. 6.** Histograms for the value of the correlation coefficient $\rho$ between the three astrometric parameters $\mu_{\alpha}^*$, $\mu_{\delta}$, and $\omega$, for the stars used in this study in the fields where the TGAS proper motions uncertainties are smaller than the UCAC4 uncertainties.</figcaption></figure>

A visual example is shown in Fig. 7, where the bottom left and bottom middle panels show the proper motions of stars in NGC 2360, displaying the uncertainties as if $\mu_{\alpha}^*$ and $\mu_{\delta}$ errors were uncorrelated. The representation in the bottom right of Fig. 7 takes into consideration the correlations between both components of the proper motion, showing the uncertainties as covariant, tilted ellipses, which when correlations are non-zero are always narrower than the non-covariant representation. As a result, some stars that appear as marginally compatible with being cluster members in the bottom left panel of Fig. 7 clearly appear as outliers in the bottom right panel.

When using UCAC4 proper motions, for which the correlations between the errors on the components of the proper motion are not available, we set the correlations coefficients to zero. We also assumed that the error on the TGAS parallax is uncorrelated with the UCAC4 proper motion errors for these stars.

# 3. Mean proper motions and parallaxes

## 3.1. Computation

We computed the mean astrometric parameters of the cluster members after applying a 2-$\sigma$ clipping on the proper motions, parallax, and sky distribution of the member stars. In order to determine a mean proper motion and mean parallax for each cluster taking into account both the membership probability and the uncertainty on the parameters of each star, we performed a thousand random drawings where we picked stars according to their probability of being cluster members. For each redrawing, the mean proper motion and mean parallax was computed as a weighted mean, where the weight for each star corresponds to the inverse variance (inverse uncertainty, squared) and the uncertainty on this mean value is the standard deviation divided by the square root of the number of selected stars. The final value is the average over all redrawings, and the final uncertainty is the quadratic sum of the average the uncertainty and the standard deviation of the mean value over all redrawings.

In this computation of mean parameters we do not take into account the correlations between astrometric parameters for individual stars, as they do not have a significant impact on the mean values and their uncertainties. However, due to the fact that correlations tend to be similar among sources in a given field and therefore do not cancel out, the mean parameters themselves are correlated, and a more complete treatment should be applied if

one was interested in determining the full covariance of the mean parameters. As a result, the volume of parameter space covered by the possible combinations of $\langle\omega\rangle$, $\langle\mu_{\alpha}^*\rangle$, and $\langle\mu_{\delta}\rangle$ is smaller when correlations exist between these quantities.

The final results are listed in Table A.1. Since the number of identified members is rather low for many clusters, we did not update the coordinates of the cluster centre. In all cases, the positions of the member stars were consistent with the literature data.

The astrometric validation of GDR1 revealed local systematic biases of the order of 0.3 mas in the parallax zero-point (Arenou et al. 2017). Although this does not affect our ability to distinguish between background and foreground stars within a small field of view, the uncertainty on the absolute parallax of a given group of stars should include an additional 0.3 mas, which cannot be reduced by averaging over a large number of cluster members.

We compared the mean cluster parallaxes found in this study with the expected values from the DAML and MWSC catalogues. We found two outliers (Collinder 359 and Ruprecht 145) for which the mean parallax is significantly smaller (by more than 3-$\sigma$) than the expected value in DAML, but which are in excellent agreement with the value listed in MWSC. Similarly, one cluster (ASCC 99) has a mean parallax higher by more than 4-$\sigma$ of the MWSC-expected value, but which is in excellent agreement with DAML. Apart from those three outliers, the standard deviation of the difference between mean and expected parallax in units of uncertainty is 0.75 using DAML references and 0.71 using MWSC. The fact that the values derived from TGAS data differ on average by less than one unit of uncertainty with the catalogue value suggests that, in most cases, adding 0.3 mas to the total error budget leads to overestimated uncertainties.

The studies of Stassun & Torres (2016) and Jao et al. (2016) report that the TGAS parallaxes could be underestimated by 0.25 mas, although they note that the effect is larger for stars within 25 pc and seems to vanish when $\omega > 1$ mas. In this present study, the median difference between the mean parallax and the parallax expected from the literature distance is -0.035 mas when using DAML as reference (M.A.D. 0.24) and -0.032 mas (M.A.D. 0.25) using MWSC. These numbers are in agreement with a zero-point of -0.04 ± 0.003 mas reported in Arenou et al. (2017), and with the results of Casertano et al. (2017) and Sesar et al. (2017), who find no significant zero-point offset of the parallaxes in the TGAS data.
NGC 2360 - Proper motions and uncertainties

<figure><img src="image_7.png" /><figcaption>**Fig. 7.** Top: UCAC4 proper motions for NGC 2360. Bottom left: TGAS proper motions for the same set of stars, with nominal uncertainties represented as error bars. Bottom middle: same as bottom left panel, but the uncertainties are represented as non-covariant error ellipses. Bottom right: same data as bottom middle panel, with uncertainties represented as tilted error ellipses, showing the correlations between $\mu_{\alpha}*$ and $\mu_{\delta}$ errors (here of the order of -0.8).</figcaption></figure>

## 3.2. Comparison with other results from TGAS parallaxes

We tested the capability of our approach to reproduce the results obtained in Gaia Collaboration (2017) with TGAS parallaxes, for the six most distant OCs in their list ($\bar{\omega} < 3$ mas). The mean parallaxes we obtain are in perfect agreement for four of those OCs (see Table 1), but differ for the other two (NGC 2516 and NGC 3532), despite our method providing a list of members very similar to that listed in Gaia Collaboration (2017). The origin of this discrepancy is still under investigation, but we remark that very strong correlations are present for the astrometric parameters of the sources in the field of these two objects.

## 3.3. Remarks on individual clusters

Chereul 1: The putative object listed in Chereul et al. (1999) can be clearly identified in proper motion space from TGAS data. Given its estimated distance of ~100 pc, a better calculation of its mean astrometric parameters taking into account the projection effects and variation of proper motions across the sky (such as in e.g. Gaia Collaboration 2017) would provide a more accurate determination. We, however, failed to identify the groups Chereul 2 and Chereul 3, which could be asterisms, and should be further investigated with the next Gaia data release.

<table><caption>**Table 1.** Mean parallaxes for the most distant clusters in FvL17.</caption><thead><tr><th>OC</th><th>⟨ω⟩<br>Gaia Collaboration (2017)</th><th>⟨ω⟩<br>(this study)</th></tr></thead><tbody><tr><td>NGC 6633</td><td>2.41 ± 0.04</td><td>2.38 ± 0.03</td></tr><tr><td>NGC 3532</td><td>2.42 ± 0.04</td><td>2.26 ± 0.05</td></tr><tr><td>NGC 2547</td><td>2.75 ± 0.08</td><td>2.73 ± 0.10</td></tr><tr><td>IC 4665</td><td>2.83 ± 0.05</td><td>2.71 ± 0.08</td></tr><tr><td>Collinder 140</td><td>2.86 ± 0.11</td><td>2.72 ± 0.09</td></tr><tr><td>NGC 2516</td><td>2.99 ± 0.08</td><td>2.51 ± 0.05</td></tr></tbody></table>

**Notes.** Mean parallax $⟨\omega⟩$ expressed in mas.

The overlapping objects NGC 1746, NGC 1750, and NGC 1758: Although they were historically listed as three distinct objects, the reality of all NGC 1746, NGC 1750, and NGC 1758 has been questioned by different authors. In particular, Galadi-Enriquez et al. (1998b,a) have suggested that NGC 1750 and NGC 1758 might be true physical objects, while NGC 1746 could be an asterism. The more recent study of Landolt & Africano (2010) was not able to solve the issue of the
reality of all three objects. The catalogue of DAML lists them
as three different components, with distinct proper motions
and distances (800 pc, 630 pc, and 760 pc, respectively) while
MWSC lists the whole complex as NGC 1746. In this paper
we only identified the system as one group, whose estimated
distance and proper motions are closest to those of NGC 1750 in
DAML, and therefore list our result under that name.

4. Age determination from photometry: Bayesian
isochrone fitting

We determined cluster parameters for our clusters with the freely
available code Bayesian Analysis for Stellar Evolution with Nine
Parameters (BASE-9 von Hippel et al. 2006). The code derives
posterior distributions for cluster parameters using a Bayesian
approach and performing sampling in parameter space with a
Markov chain Monte Carlo technique (MCMC). Given a set of
stellar isochrones, BASE-9 compares the theoretical photome-
try with the observed one, taking into account the presence of
unresolved binary stars as well as the probability that each indi-
vidual star is a cluster member or a field star. More insight
on the capabilities of BASE-9 can be found in Jeffery et al.
(2016). We set BASE-9 to provide a posterior distribution for
four parameters: age (log *t*), metallicity ([Fe/H]), distance mod-
ulus ((*m* − *M*)<sub>*V*</sub>), and *V*-band extinction (*A*<sub>*V*</sub> = 3.1 × *E*(*B* − *V*)).
We made use of PARSEC isochrones (Bressan et al. 2012),
which are not shipped with BASE-9 but were straightforward to
implement.

Although the *G*-magnitudes contained in the GDR1 cat-
alogue are of exquisite quality, we decided not to make use
of them because the common approximation that all stars are
affected equally by interstellar extinction does not hold for such
precise photometric measurements. In the presence of interstel-
lar extinction, it is common to assume that all stars are affected
in the same way regardless of their spectral type, and that absorp-
tion produces an identical, rigid shift in colour and magnitude for
all stars in a cluster. In reality, stars of different spectral types are
affected differently by extinction (Jordi et al. 2010), with varia-
tions of up to 10–15 mmag in *G*-band, even in cases of moderate
extinction with *A*<sub>*V*</sub> ~ 0.5 (Sordo et al., in prep.). This effect must
be accounted for when working with *Gaia* photometry, which for
the subset of TGAS stars (brighter than *G* = 13) has a median
photometric error of 1.3 mmag. Adapting BASE-9 to the specific
requirements of sub-millimag photometry is beyond the scope of
this study, and we limited the analysis to 2MASS *JHK*<sub>S</sub> photom-
etry using the pre-computed cross-match between the *Gaia* and
2MASS catalogues provided with the GDR1 (Marrese et al. in
prep.).

A Bayesian approach requires the setting of priors on the
parameters we are trying to determine and BASE-9 only allows
for either flat or Gaussian priors. We set a Gaussian prior on
distance modulus, centred on the value found in Sect. 3 using
the mean parallax. We set rather loose Gaussian priors on the
other three parameters, with log *t* centred on the value listed in
the MWSC catalogue and with a dispersion of 0.5. The prior
on the extinction *A*<sub>V</sub> was also centred on the literature value,
with a dispersion of 0.2. The metallicity prior distribution was
centred on [Fe/H] = −0.1 for all stars, with a dispersion of
0.2. When spectroscopic metallicities were available, we cen-
tred the prior on that value, with a dispersion of 0.2 dex as well.
BASE-9 also allows us to take into account prior knowledge of
membership probability, for which we used the results obtained
in Sect. 2.

<figure><img src="image_8.png" /><figcaption>**Fig. 8.** Posterior distribution maps for NGC 2567 combining the outputs of five BASE-9 runs. The dashed contour encircles 68% of the total likelihood and the crossed symbol shows the mean value.</figcaption></figure>

Following the approach of Jeffery et al. (2016), for every
OC we performed five runs of BASE-9 sampling 3000 points
each, using slightly different starting points. The first run started
from the literature value for age, extinction, and metallicity, and
the distance modulus used was the one determined in Sect. 3.
The starting points for the additional four runs were identical,
but we shifted either the initial distance modulus by ±0.3 or the
age by ±0.2 in logarithmic scale. These values correspond to the
average accuracy in the MWSC Catalogue.

In most cases the output of BASE-9 proved sensitive to the
starting values, which we attribute to the low number of mem-
ber stars in most of our clusters. When the inclusion or rejection
of one star can have a strong influence on the choice of the best
fit isochrone, the likelihood space can have a chaotic structure
with deep local minima, and the choice of initial point for the
sampling can have a strong influence on the final results. For
26 OCs, the procedure converged to similar posterior distri-
butions in all five runs and provided a satisfactory fit to the
observed colour-magnitude diagram (CMD). Unsurprisingly,
these clusters tend to be those with a larger number of mem-
bers or featuring red clump stars, which provide good constrains
on the cluster parameters.

For those 26 OCs, we combined the results of all five
runs to compute the final cluster parameters and their uncer-
tainty (effectively using a total of 15 000 samplings). The pos-
terior distributions are generally non-symmetrical and show
correlations between parameters (e.g. age and extinction, see
Fig. 8). For simplicity we only report here the mean and
standard deviation of those posterior distributions. As pointed
out by Jeffery et al. (2016), the uncertainties on the cluster
parameters reflect the internal precision of the procedure (the
certainty with which a certain PARSEC isochrone represents
the data better than a different PARSEC isochrone), rather
than an absolute accuracy. An example of CMD, along with
isochrones sampled from the posterior distribution, is shown
in Fig. 9.

Our age (log *t*) determinations are on average older by
0.26 than the literature values listed in the MWSC. System-
atic differences are expected between studies that make use
<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Colour-magnitude diagram for NGC 2567. The grey lines are five PARSEC isochrones randomly chosen from the posterior distribution returned by BASE-9.</figcaption></figure>

of different sets of models, and the MWSC ages (Kharchenko et al. 2013) were determined from PADOVA isochrones, while this study makes use of PARSEC models. In particular, different choices of solar metallicity reference and mixing-length parameter are known to lead to slight differences in the predicted brightness of red clump stars, which in turn has an incidence on the choice of model that best reproduces the morphology of the CMD (see e.g. Sect. 7.1 in Cantat-Gaudin et al. 2014).

Metallicity is the parameter that is least well constrained by isochrone fitting. The mean metallicity derived from photometry for this sample of OCs is [Fe/H] = -0.23 dex. We find that our determination of the metallicity is systematically lower by Δ[Fe/H] = 0.19 in comparison to MWSC values. Although a significant scatter is present, our sample traces a marginally negative metallicity gradient, with slope -0.04 ± 0.04 dex kpc⁻<fnref n="1" />, as shown in the bottom panel of Fig. 10. The uncertainty on the gradient slope was computed as the standard deviation among 1000 redrawings, performed by picking a metallicity and a mean parallax from a Gaussian distribution representing the values found in this study and their associated errors.

This value is compatible with the results of Netopil et al. (2016) or Jacobson et al. (2016) (-0.10 ± 0.02 and -0.085 ± 0.017 dex kpc⁻<fnref n="1" />, respectively). Our sample is, however, rather small (26 objects) and the precision of metallicites derived from isochrones fitting is not as good as the results that can be obtained from high-resolution spectroscopy.

**5. Computing three-dimensional velocities and full orbits**

In order to compute full, three-dimensional velocities, we combined the proper motions determined in this study with spectroscopic radial velocities listed in Mermilliod et al. (2008) and Mermilliod et al. (2009). After excluding the non-members and the stars flagged as either variables or binaries, we found that

<figure><img src="image_10.png" /><figcaption>**Fig. 10.** *Top:* difference in log t between the ages found in this study and those quoted in the MWSC catalogue. The dashed line indicates the mean value. *Middle:* Difference between the [Fe/H] found in this study and those quoted in the MWSC catalogue (for the OCs with metallicity estimates). The dashed blue line indicates the mean value. *Bottom:* metallicity ([Fe/H]) obtained in this study based on JHKs photometry, against Galactocentric radius. The red cross indicates the solar metallicity and Galactocentric radius (8.34 kpc, Reid et al. 2014).</figcaption></figure>

36 of the OCs in our sample have radial velocities that can be computed from at least two stars.

The velocity of a particle in the Galactic disk can be described as the sum of two components: (i) a circular motion around the Galactic centre, at a velocity depending on the Galactocentric radius and which defines the Regional Standard of Rest (RSR), and (ii) an additional peculiar motion with respect to this RSR. We computed the cylindrical (U_s, V_s, W_s)<fnref n="1" /> components of the peculiar motion as described in Sect. 3.5.1 of Casamiquela et al. (2016), adopting a Galactic rotation curve with Θ₀ = 240 km s⁻<fnref n="1" />, R☉ = 8.34 kpc, and dR/dθ = -0.2 km s⁻<fnref n="1" /> kpc⁻<fnref n="1" />, and taking into account the Sun’s own peculiar motion as (U☉, V☉, W☉) = (10.7, 15.6, 8.9) km s⁻<fnref n="1" /> (Reid et al. 2014). The results of this computation are listed in Table 3 and shown in Fig. 11. We found velocity dispersions of 21.1, 14.8, and 8.8 km s⁻<fnref n="1" /> for U_s, V_s, and W_s (respectively). These values are in good agreement with the values of 20, 15, and 10 km s⁻<fnref n="1" /> found by Holmberg et al. (2009) for stars younger than 2 Gyr.

<footnote n="1">¹ The velocity component *U*<sub>S</sub> is positive towards the Galactic centre, *V*<sub>S</sub> along the Galactic rotation, and *W*<sub>S</sub> towards the Galactic north pole.</footnote>
<table><caption>**Table 2.** Cluster parameters derived from *JHK* photometry for 26 OCs.</caption><thead><tr><th>OC</th><th>log t</th><th>[Fe/H]</th><th>A<sub>V</sub></th><th>dist. mod.</th></tr></thead><tbody><tr><td>Alessi 2</td><td>8.96 ± 0.08</td><td>-0.27 ± 0.17</td><td>0.45 ± 0.13</td><td>8.52 ± 0.20</td></tr><tr><td>IC 4756</td><td>9.09 ± 0.01</td><td>-0.09 ± 0.03</td><td>0.31 ± 0.04</td><td>8.24 ± 0.05</td></tr><tr><td>NGC 0752</td><td>9.17 ± 0.02</td><td>0.03 ± 0.05</td><td>0.11 ± 0.05</td><td>8.19 ± 0.07</td></tr><tr><td>NGC 1528</td><td>8.94 ± 0.04</td><td>-0.09 ± 0.11</td><td>0.30 ± 0.10</td><td>9.49 ± 0.14</td></tr><tr><td>NGC 1662</td><td>9.03 ± 0.03</td><td>-0.41 ± 0.07</td><td>0.66 ± 0.06</td><td>7.60 ± 0.08</td></tr><tr><td>NGC 1750</td><td>8.58 ± 0.09</td><td>-0.43 ± 0.05</td><td>0.94 ± 0.07</td><td>9.09 ± 0.09</td></tr><tr><td>NGC 1912</td><td>8.85 ± 0.10</td><td>-0.36 ± 0.10</td><td>0.69 ± 0.10</td><td>9.77 ± 0.13</td></tr><tr><td>NGC 2099</td><td>8.95 ± 0.05</td><td>-0.19 ± 0.08</td><td>0.50 ± 0.08</td><td>9.69 ± 0.14</td></tr><tr><td>NGC 2281</td><td>8.85 ± 0.04</td><td>-0.14 ± 0.11</td><td>0.29 ± 0.06</td><td>8.40 ± 0.12</td></tr><tr><td>NGC 2482</td><td>8.88 ± 0.06</td><td>-0.20 ± 0.15</td><td>0.15 ± 0.08</td><td>9.76 ± 0.16</td></tr><tr><td>NGC 2527</td><td>9.03 ± 0.04</td><td>-0.26 ± 0.11</td><td>0.24 ± 0.05</td><td>8.73 ± 0.10</td></tr><tr><td>NGC 2539</td><td>8.96 ± 0.15</td><td>-0.18 ± 0.16</td><td>0.29 ± 0.15</td><td>9.97 ± 0.22</td></tr><tr><td>NGC 2548</td><td>8.90 ± 0.04</td><td>-0.38 ± 0.09</td><td>0.28 ± 0.05</td><td>9.07 ± 0.08</td></tr><tr><td>NGC 2567</td><td>8.87 ± 0.07</td><td>-0.23 ± 0.13</td><td>0.41 ± 0.09</td><td>10.68 ± 0.18</td></tr><tr><td>NGC 4852</td><td>8.89 ± 0.08</td><td>-0.13 ± 0.15</td><td>0.26 ± 0.12</td><td>9.37 ± 0.19</td></tr><tr><td>NGC 5822</td><td>9.15 ± 0.02</td><td>-0.24 ± 0.08</td><td>0.41 ± 0.10</td><td>9.32 ± 0.09</td></tr><tr><td>NGC 6152</td><td>8.75 ± 0.09</td><td>-0.03 ± 0.15</td><td>0.69 ± 0.14</td><td>10.34 ± 0.27</td></tr><tr><td>NGC 6281</td><td>8.80 ± 0.09</td><td>-0.38 ± 0.07</td><td>0.63 ± 0.07</td><td>8.19 ± 0.11</td></tr><tr><td>NGC 6793</td><td>9.07 ± 0.11</td><td>-0.27 ± 0.19</td><td>0.47 ± 0.10</td><td>8.56 ± 0.20</td></tr><tr><td>NGC 6811</td><td>9.16 ± 0.03</td><td>-0.38 ± 0.06</td><td>0.36 ± 0.07</td><td>9.75 ± 0.10</td></tr><tr><td>NGC 6991</td><td>9.15 ± 0.02</td><td>-0.00 ± 0.06</td><td>0.25 ± 0.06</td><td>8.82 ± 0.11</td></tr><tr><td>NGC 7209</td><td>9.01 ± 0.06</td><td>-0.41 ± 0.07</td><td>0.34 ± 0.10</td><td>9.45 ± 0.15</td></tr><tr><td>Platais 3</td><td>8.90 ± 0.15</td><td>-0.22 ± 0.14</td><td>0.17 ± 0.08</td><td>6.23 ± 0.11</td></tr><tr><td>Ruprecht 98</td><td>8.96 ± 0.05</td><td>-0.29 ± 0.11</td><td>0.53 ± 0.09</td><td>8.01 ± 0.12</td></tr><tr><td>Stock 1</td><td>8.77 ± 0.04</td><td>-0.12 ± 0.11</td><td>0.27 ± 0.07</td><td>7.76 ± 0.10</td></tr><tr><td>Turner 5</td><td>8.80 ± 0.08</td><td>-0.19 ± 0.15</td><td>0.11 ± 0.08</td><td>7.58 ± 0.17</td></tr></tbody></table>

<figure><img src="image_11.png" /><figcaption>**Fig. 11.** Distribution of the three components of the peculiar motion for the OCs listed in Table 3. NGC 5617 was discarded because of the large uncertainty affecting its distance determination.</figcaption></figure>

From the three-dimensional position and three-dimensional velocity of each cluster, we computed orbits using the software galpy and the static, axisymmetric potential MWPotential2014 (Bovy 2015). Figure 12 shows the current distance (|z|) and maximum altitude above the Galactic plane (z<sub>max</sub>) for those integrated orbits. Both these quantities correlate with age, with OCs younger than 300 Myr being all contained within 180 pc of the plane, while half the older clusters have orbits that extend beyond this limit. The oldest OC in our sample is NGC 2682, one of the oldest known open clusters, and its orbit strays more than 400 pc from the Galactic plane. We also computed the eccentricity of each orbit ($e = r_a/r_p$ where $r_p$ and $r_a$ are the perigalacticon and apogalacticon of the

orbit), and found no apparent correlation of eccentricity with age (bottom panel of Fig. 12). The orbital parameters of the integrated orbits are listed in Table 3.

## 6. Discussion

The two-million-star TGAS sample allows us to identify members and derive mean parallaxes in clusters located further than 1000 pc from us, thus enabling studies of the Galactic disk on a larger scale than the works based on the results of the HIPPARCOS catalogue (e.g. Robichon et al. 1999; Bouy & Alves 2015), which is severely incomplete beyond 100 pc and for which the relative parallax error reaches 50% at 350 pc.
<table><caption>**Table 3.** Selected parameters for 36 integrated orbits.</caption><thead><tr><th>OC</th><th>Log age</th><th>U<sub>s</sub></th><th>V<sub>s</sub></th><th>W<sub>s</sub></th><th>z<br>[pc]</th><th>z<sub>max</sub><br>[pc]</th><th>e</th></tr></thead><tbody><tr><td>IC 4725</td><td>7.97</td><td>13.7 ± 0.5</td><td>-4.4 ± 1.4</td><td>10.4 ± 1.1</td><td>-25 ± 9</td><td>92 ± 9</td><td>0.086 ± 0.021</td></tr><tr><td>IC 4756</td><td>9.09</td><td>-15.3 ± 0.3</td><td>-3.4 ± 0.3</td><td>2.3 ± 0.3</td><td>72 ± 5</td><td>72 ± 5</td><td>0.087 ± 0.004</td></tr><tr><td>NGC 0752</td><td>9.17</td><td>-13.1 ± 0.3</td><td>-6.0 ± 0.5</td><td>-9.5 ± 0.3</td><td>-141 ± 14</td><td>229 ± 30</td><td>0.069 ± 0.009</td></tr><tr><td>NGC 1342</td><td>8.6</td><td>11.9 ± 0.4</td><td>8.5 ± 0.6</td><td>8.1 ± 0.5</td><td>-170 ± 19</td><td>196 ± 17</td><td>0.035 ± 0.003</td></tr><tr><td>NGC 1647</td><td>8.3</td><td>18.8 ± 0.3</td><td>14.0 ± 0.3</td><td>6.5 ± 0.3</td><td>-125 ± 11</td><td>148 ± 8</td><td>0.06 ± 0.0</td></tr><tr><td>NGC 1662</td><td>9.03</td><td>26.1 ± 0.4</td><td>15.5 ± 0.2</td><td>10.2 ± 0.3</td><td>-121 ± 10</td><td>177 ± 7</td><td>0.084 ± 0.001</td></tr><tr><td>NGC 2099</td><td>8.95</td><td>-1.1 ± 0.3</td><td>-30.6 ± 4.0</td><td>1.2 ± 1.8</td><td>81 ± 13</td><td>93 ± 24</td><td>0.138 ± 0.042</td></tr><tr><td>NGC 2168</td><td>8.26</td><td>23.1 ± 0.3</td><td>0.1 ± 0.6</td><td>2.8 ± 0.5</td><td>59 ± 7</td><td>64 ± 5</td><td>0.081 ± 0.01</td></tr><tr><td>NGC 2281</td><td>8.85</td><td>-15.1 ± 0.2</td><td>4.8 ± 0.5</td><td>-3.1 ± 0.6</td><td>176 ± 18</td><td>188 ± 25</td><td>0.059 ± 0.006</td></tr><tr><td>NGC 2360</td><td>8.8</td><td>-4.8 ± 1.1</td><td>17.9 ± 1.1</td><td>24.6 ± 1.2</td><td>-9 ± 8</td><td>415 ± 120</td><td>0.055 ± 0.027</td></tr><tr><td>NGC 2423</td><td>9.03</td><td>26.4 ± 0.5</td><td>-7.9 ± 0.5</td><td>-0.8 ± 0.6</td><td>88 ± 10</td><td>104 ± 20</td><td>0.12 ± 0.013</td></tr><tr><td>NGC 2447</td><td>8.68</td><td>-1.4 ± 0.8</td><td>10.8 ± 0.6</td><td>5.6 ± 0.6</td><td>27 ± 0</td><td>68 ± 9</td><td>0.023 ± 0.01</td></tr><tr><td>NGC 2477</td><td>8.91</td><td>35.7 ± 4.4</td><td>10.7 ± 2.4</td><td>9.6 ± 5.2</td><td>-115 ± 49</td><td>198 ± 59</td><td>0.098 ± 0.032</td></tr><tr><td>NGC 2527</td><td>9.03</td><td>-15.1 ± 0.3</td><td>-8.6 ± 0.2</td><td>6.8 ± 0.2</td><td>42 ± 1</td><td>97 ± 1</td><td>0.107 ± 0.004</td></tr><tr><td>NGC 2539</td><td>8.96</td><td>8.3 ± 1.0</td><td>-9.3 ± 0.9</td><td>2.5 ± 1.7</td><td>196 ± 54</td><td>211 ± 62</td><td>0.106 ± 0.009</td></tr><tr><td>NGC 2546</td><td>8.13</td><td>4.2 ± 1.3</td><td>10.9 ± 0.8</td><td>-1.9 ± 0.8</td><td>-13 ± 9</td><td>27 ± 12</td><td>0.022 ± 0.005</td></tr><tr><td>NGC 2548</td><td>8.9</td><td>14.4 ± 0.6</td><td>15.8 ± 0.6</td><td>14.2 ± 0.7</td><td>211 ± 27</td><td>297 ± 27</td><td>0.05 ± 0.005</td></tr><tr><td>NGC 2567</td><td>8.87</td><td>10.4 ± 1.7</td><td>-5.5 ± 1.0</td><td>-2.3 ± 1.2</td><td>123 ± 48</td><td>156 ± 68</td><td>0.097 ± 0.016</td></tr><tr><td>NGC 2682</td><td>9.54</td><td>-19.0 ± 1.7</td><td>-11.0 ± 1.3</td><td>-3.0 ± 2.9</td><td>392 ± 80</td><td>403 ± 93</td><td>0.111 ± 0.006</td></tr><tr><td>NGC 3680</td><td>9.2</td><td>-1.8 ± 2.5</td><td>1.8 ± 0.7</td><td>-0.3 ± 1.0</td><td>264 ± 34</td><td>266 ± 33</td><td>0.036 ± 0.003</td></tr><tr><td>NGC 5138</td><td>7.55</td><td>0.4 ± 6.5</td><td>-7.8 ± 3.5</td><td>5.2 ± 1.6</td><td>98 ± 21</td><td>108 ± 27</td><td>0.06 ± 0.021</td></tr><tr><td>NGC 5316</td><td>8.23</td><td>-39.5 ± 5.0</td><td>-26.0 ± 3.1</td><td>21.2 ± 0.9</td><td>27 ± 0</td><td>84 ± 3</td><td>0.05 ± 0.028</td></tr><tr><td>NGC 5617</td><td>8.25</td><td>-62.5 ± 20.4</td><td>-35.2 ± 11.1</td><td>22.8 ± 2.9</td><td>21 ± 1</td><td>64 ± 30</td><td>0.073 ± 0.029</td></tr><tr><td>NGC 5662</td><td>8.28</td><td>-22.0 ± 1.0</td><td>-3.1 ± 0.9</td><td>0.9 ± 0.4</td><td>69 ± 5</td><td>128 ± 19</td><td>0.04 ± 0.008</td></tr><tr><td>NGC 5822</td><td>9.15</td><td>-31.0 ± 0.8</td><td>-6.2 ± 0.8</td><td>13.5 ± 0.4</td><td>78 ± 10</td><td>81 ± 9</td><td>0.066 ± 0.01</td></tr><tr><td>NGC 6067</td><td>7.97</td><td>-0.7 ± 4.5</td><td>-7.6 ± 4.9</td><td>14.4 ± 2.5</td><td>-88 ± 77</td><td>117 ± 61</td><td>0.1 ± 0.071</td></tr><tr><td>NGC 6124</td><td>8.29</td><td>-3.5 ± 0.4</td><td>20.3 ± 0.7</td><td>7.0 ± 0.7</td><td>88 ± 7</td><td>107 ± 6</td><td>0.037 ± 0.002</td></tr><tr><td>NGC 6134</td><td>9.02</td><td>-2.8 ± 1.1</td><td>17.8 ± 2.1</td><td>-8.3 ± 2.6</td><td>21 ± 0</td><td>99 ± 50</td><td>0.025 ± 0.007</td></tr><tr><td>NGC 6281</td><td>8.8</td><td>6.9 ± 0.4</td><td>6.2 ± 0.9</td><td>9.7 ± 0.9</td><td>47 ± 2</td><td>89 ± 8</td><td>0.038 ± 0.005</td></tr><tr><td>NGC 6494</td><td>8.52</td><td>-2.4 ± 0.2</td><td>16.7 ± 0.9</td><td>4.5 ± 1.0</td><td>71 ± 7</td><td>76 ± 3</td><td>0.019 ± 0.005</td></tr><tr><td>NGC 6705</td><td>8.51</td><td>35.5 ± 4.7</td><td>-4.2 ± 7.0</td><td>4.4 ± 4.0</td><td>-45 ± 50</td><td>102 ± 70</td><td>0.157 ± 0.039</td></tr><tr><td>NGC 6811</td><td>8.8</td><td>28.3 ± 3.5</td><td>16.4 ± 0.5</td><td>22.1 ± 1.4</td><td>206 ± 38</td><td>281 ± 47</td><td>0.077 ± 0.004</td></tr><tr><td>NGC 6866</td><td>8.64</td><td>7.7 ± 4.6</td><td>27.1 ± 0.6</td><td>8.8 ± 3.8</td><td>175 ± 46</td><td>212 ± 57</td><td>0.078 ± 0.001</td></tr><tr><td>NGC 6940</td><td>8.98</td><td>24.8 ± 1.6</td><td>9.5 ± 0.5</td><td>-10.5 ± 0.9</td><td>-99 ± 13</td><td>205 ± 36</td><td>0.075 ± 0.005</td></tr><tr><td>NGC 7209</td><td>9.01</td><td>-22.4 ± 1.2</td><td>-5.4 ± 0.4</td><td>5.1 ± 0.9</td><td>-81 ± 23</td><td>114 ± 11</td><td>0.09 ± 0.013</td></tr><tr><td>Trumpler 3</td><td>8.01</td><td>18.2 ± 1.6</td><td>31.2 ± 1.8</td><td>-8.0 ± 1.1</td><td>78 ± 9</td><td>84 ± 12</td><td>0.054 ± 0.005</td></tr></tbody></table>

Notes. The velocity $U_s$ is positive towards the Galactic centre, $V_s$ along the Galactic rotation, and $W_s$ towards the Galactic north pole. The uncertainties were derived from 10 000 random realisations, varying the cluster distance, proper motion, and radial velocity according to their uncertainties.

The distribution of OCs identified in this paper is shown in Fig. 13, along with the schematic location of the spiral arms of the Milky Way in the model of Reid et al. (2014). We remark that almost all OCs in this study are located in the inter-arm region. In order to trace the spiral structure of the thin disk, our sample would require a significant number of objects younger than 20 Myr (log t < 7.3), as clusters drift away from their birthplace over this timescale (see e.g. Dias & Lépine 2005). These young clusters and associations are often sparsely populated and/or embedded in their progenitor molecular cloud (thus requiring the use of infrared photometry). Due to the magnitude limits of the TGAS solution, even nearby clusters such as Chamaeleon I (2 Myr old, 160 pc from us, Whittet et al. 1997) are not present in the catalogue. Hints of objects such as IC 2395 (6 Myr, 800 pc Clariá et al. 2003) can be seen, but not with a sufficient number of stars for our method to establish secure membership.

Another common tracer of the spiral structure of the disk are the so-called OB associations, which are even sparser groups of bright, short-lived blue stars, whose identification requires the use of multi-band photometry (and ideally spectroscopy). Such a search is well beyond the scope of this paper, aimed at identifying stellar clusters from astrometric data.

Out of the 128 OCs for which we derived mean astrometric parameters, 67 have mean proper motions obtained from TGAS data (with a median proper motion error of 0.23 mas yr⁻¹), while UCAC4 data was used for 61 of them (these OCs have a median proper motion error of 0.32 mas yr⁻¹).

The difference in scale height observed between young and old objects is traditionally attributed to disk heating: stars and stellar clusters are formed on flat and circular orbits, which might be disturbed during the course of their life, in particular by interactions with giant molecular clouds and non-axisymmetric
<figure><img src="image_12.png" /><figcaption>**Fig. 12.** *Top:* current distance $|z|$ from the Galactic plane as a function of age for the 36 OCs for which we computed full orbits. *Middle:* maximum altitude above the Galactic plane for the integrated orbits of those OCs. *Bottom:* eccentricity of the integrated orbits against age of the cluster.</figcaption></figure>

features of the disk (Carlberg 1987; Aumer et al. 2016; Grand et al. 2016). The clear correlation between age and altitude above the plane shown in Fig. 12, where clusters younger than 300 Myr have a mean value of $z_{max}$ of 100 pc while the older clusters have a mean $z_{max}$ of 190 pc, illustrates the fact that their vertical motion can be affected during the first hundred million years of their life.

The absence of an apparent correlation between age and eccentricity seems to indicate that radial heating has not affected the objects in our sample as much as vertical heating, and the timescale for radial heating is longer than that of vertical heating. The sample of old outer disk clusters studied by Cantat-Gaudin et al. (2016) contains several OCs older than 2 Gyr and eccentricities larger than 0.2. The study of Vande Putte et al. (2010) also shows that OCs older than 1 Gyr tend to have larger eccentricities, but contains a significant number of young objects that appear to follow eccentric orbits as well. Both these studies might be affected by an observational bias, since old stellar

<figure><img src="image_13.png" /><figcaption>**Fig. 13.** Position of the OCs studied in this paper in Galactic rectangular XYZ coordinates. Black dots: OCs from this study with relative parallax errors under 50% (considering an additional 0.3 mas calibration error for all OCs). Crosses: OCs from this study with relative parallax error over 50%. Open symbols: OCs from Gaia Collaboration (2017). The spiral arms are traced according to the model of Reid et al. (2014).</figcaption></figure>

clusters are more easily detected at high Galactic latitudes, towards the less crowded anticentre direction, and therefore more likely to be those with perturbed orbits. The uncertainties on proper motions for these distant objects also translate into large uncertainties on their physical motion. The unprecedented quality of the Gaia-DR2 data will allow for a more accurate characterisation of the kinematics for a much larger sample, giving us a deeper insight into the age-kinematics relations over a large portion of the Galactic disk.

## 7. Conclusion and future prospects

In this paper, we make use of a suitable combination of Gaia DR1 TGAS parallaxes and proper motions and UCAC4 proper motions to derive a list of high probability cluster members for 128 OCs. For those objects, we compute mean proper motions and parallaxes. For 26 clusters, we obtain parameters such as age, extinction, [Fe/H], and distance modulus from comparison with isochrones using a Bayesian/MCMC method.

The issue of dealing with a large astronomical dataset is not unique to Gaia, as other current (Chambers et al. 2016, e.g. the Panoramic Survey Telescope and Rapid Response System) or future (Ivezic et al. 2008, e.g. the Large Synoptic Survey Telescope) observational campaigns will deliver multidimensional catalogues of unprecedented size. The work conducted in this study shows the strength of automatic approaches
to cluster membership selection such as UPMASK when dealing with large datasets. While the two-million-star TGAS subset of GDR1 is larger than the HIPPARCOS catalogue by an order of magnitude, and therefore already too rich to be investigated by hand, the second Gaia data release (GDR2) will contain a full astrometric solution for one billion star, dwarfing TGAS by a factor 500. The upcoming ground-based LSST is expected to reach four magnitudes deeper than Gaia, and to deliver better astrometry than Gaia for sources fainter than $G \sim 20$ as well as multi-band photometry for an end-of-mission catalogue totalling many billions of objects.

Systematic cluster studies conducted with automated methods should also attempt to provide general cluster parameters such as ages, using stellar evolution models. Ideally, these determinations should be performed for as many clusters as possible in order to build large homogeneous samples and avoid the additional dispersion in results introduced by compiling age determinations originating from various studies making use of various sets of models. Automated tools provide an objective estimate of the cluster parameters, a convenient alternative to fitting colour-magnitude diagrams by eye (which yields non-reproducible results and is completely impractical when dealing with samples of hundreds of clusters), and allows us to consider independent measurements such as distance estimates obtained from trigonometric parallaxes.

This preparatory study shows the strength of our tools, applied to astrometric data, to identify and study stellar clusters. The upcoming GDR2 dataset will be much deeper than TGAS (with magnitudes reaching down to $G \sim 21$, against $G \sim 13$ in TGAS), with a better astrometric precision for both proper motions and parallaxes, allowing us to identify and characterise clusters at much larger distances. The data of GDR2 will also contain Gaia $G_{BP}$ and $G_{RP}$ magnitudes, allowing us to perform photometric studies using Gaia data alone. Efforts should be made towards developing tools to take into account the effect of interstellar extinction on stars of different spectral types (see Sordo et al., in prep.), in order to take full advantage of the sub-millimag precision of Gaia magnitudes.

In addition to providing better characterisations of known OCs, the GDR2 data is expected to contain many as yet unknown clusters and associations. Completing the cluster census in the Galactic disk and estimating ages is crucial to understanding cluster formation and disruption, and for tracing the structures and the processes that shape the disk of the Milky Way.

**Acknowledgements.** Part of the work presented in this paper was support by the European Science Foundation under ESF grant number 5004. This work has made use of data from the European Space Agency (ESA) mission Gaia (www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. T. C.-G. and A. V. acknowledge the support of ASI through grant 2014-025-R.1.2015. A. M. and A. K.-M. acknowledge support from the Portugese Fundação para a Ciência e a Tecnologia (FCT) through grants: SFRH/BPD/74697/2010 and the Strategic Programme UID/FIS/00099/2013. This work was supported by the MINECO (Spanish Ministry of Economy) - FEDER through grant ESP2016-80079-C2-1-R, ESP2014-55996-C2-1-R, and MDM-2014-0369 of ICCUB (Unidad de Excelencia 'María de Maeztu'), as well as the Juan de la Cierva program 2015. This work has made extensive use of Topcat (Taylor 2005), and of NASA's Astrophysics Data System. The figures in this paper were produced with Matplotlib

(Hunter 2007) and Healpy, a Python implementation of HEALPix (Górski et al. 2005).

## References

Allison, R. J., Goodwin, S. P., Parker, R. J., et al. 2009, MNRAS, 395, 1449

Arenou, F., Luri, X., Babusiaux, C., et al. 2017, A&A, 599, A50

Aumer, M., Binney, J., & Schönrich, R. 2016, MNRAS, 462, 1697

Bica, E., & Bonatto, C. 2011, A&A, 530, A32

Bouy, H., & Alves, J. 2015, A&A, 584, A26

Bovy, J. 2015, ApJS, 216, 29

Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127

Cantat-Gaudin, T., Vallenari, A., Zaggia, S., et al. 2014, A&A, 569, A17

Cantat-Gaudin, T., Donati, P., Vallenari, A., et al. 2016, A&A, 588, A120

Carlberg, R. G. 1987, ApJ, 322, 59

Casamiquela, L., Carrera, R., Jordi, C., et al. 2016, MNRAS, 458, 3150

Casertano, S., Riess, A. G., Bucciarelli, B., et al. 2017, A&A, 599, A67

Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560]

Chereul, E., Crézé, M., & Bienaymé, O. 1999, A&AS, 135, 5

Clariá, J. J., Lapasset, E., Piatti, A. E., et al. 2003, A&A, 409, 541

Costa, E., Moitinho, A., Radiszc, M., et al. 2015, A&A, 580, A4

Dias, W. S., Alessi, B. S., Moitinho, A., et al. 2002, A&A, 389, 871

Dias, W. S., & Lépine, J. R. D. 2005, ApJ, 629, 825

ESA (ed.) 1997, The HIPPARCOS and TYCHO catalogues. Astrometric and photometric star catalogues derived from the ESA HIPPARCOS Space Astrometry Mission, ESA Special Publication, 1200

Gaia Collaboration (Brown, A. G. A., et al.) 2016a, A&A, 595, A2

Gaia Collaboration (Prusti, T., et al.) 2016b, A&A, 595, A1

Gaia Collaboration (van Leeuwen, F., et al.) 2017, A&A, 601, A19

Galadi-Enriquez, D., Jordi, C., Trullols, E., et al. 1998a, A&AS, 131, 239

Galadi-Enriquez, D., Jordi, C., Trullols, E., et al. 1998b, A&A, 333, 471

Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759

Grand, R. J. J., Springel, V., Gómez, F. A., et al. 2016, MNRAS, 459, 199

Hög, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27

Holmberg, J., Nordström, B., & Andersen, J. 2009, A&A, 501, 941

Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90

Ivezic, Z., Tyson, J. A., Abel, B., et al. 2008, ArXiv e-prints [arXiv:0805.2366]

Jacobson, H. R., Friel, E. D., Jílková, L., et al. 2016, A&A, 591, A37

Jao, W.-C., Henry, T. J., Riedel, A. R., et al. 2016, ApJ, 832, L18

Jeffery, E. J., von Hippel, T. van Dyk, D. A., et al. 2016, ApJ, 828, 79

Jordi, C., Gebran, M., Carrasco, J. M., et al. 2010, A&A, 523, A48

Kharchenko, N. V., Piskunov, A. E., Schilbach, E., et al. 2013, A&A, 558, A53

Koposov, S. E., Belokurov, V., & Torrealba, G. 2017, MNRAS, 470, 2702

Krone-Martins, A., & Moitinho, A. 2014, A&A, 561, A57

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

Landolt, A. U., & Africano, III, J. L. 2010, PASP, 122, 1008

Lindgren, L., Lammers, U., Bastian, U., et al. 2016, A&A, 595, A4

Mermilliod, J. C., Mayor, M., & Udry, S. 2008, A&A, 485, 303

Mermilliod, J.-C., Mayor, M., & Udry, S. 2009, A&A, 498, 949

Messier, C. 1781, Catalogue des Nébuleuses & des amas d'Étoiles (Catalog of Nebulae and Star Clusters), Technical report

Minchev, I. 2016, Astron. Nachr., 337, 703

Netopil, M., Paunzen, E., Heiter, U., et al. 2016, A&A, 585, A150

Perryman, M. A. C., Lindgren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49

Perryman, M. A. C., de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130

Robichon, N., Arenou, F., Mermilliod, J.-C., et al. 1999, A&A, 345, 471

Sesar, B., Fouesneau, M., Price-Whelan, A. M., et al. 2017, ApJ, 838, 107

Stassun, K. G., Torres, G. 2016, ApJ, 831, L6

Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, ASP Conf. Ser., 347, 29

van Leeuwen, F. 2007, A&A, 474, 653

Vande Putte, D., Garnier, T. P., Ferreras, I., et al. 2010, MNRAS, 407, 2109

Vázquez, R. A., Moitinho, A., Carraro, G., et al. 2010, A&A, 511, A38

von Hippel, T., Jefferys, W. H., Scott, J., et al. 2006, ApJ, 645, 1436

Whittet, D. C. B., Prusti, T., Franco, G. A. P., et al. 1997, A&A, 327, 1194

Zacharias, N., Finch, C. T., Girard, T. M., et al. 2012, VizieR Online Data Catalog, I/322A

Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44
Appendix A

<table><caption>**Table A.1.** Mean astrometric parameters computed for 128 OCs.</caption><thead><tr><th>OC</th><th>α<br>[deg]</th><th>δ<br>[deg]</th><th>⟨ω⟩<br>[mas]</th><th>ω<sub>DAML</sub><br>[mas]</th><th>ω<sub>MWSC</sub><br>[mas]</th><th>⟨μ<sub>α</sub>⟩<br>[mas yr<sup>-1</sup>]</th><th>⟨μ<sub>δ</sub>⟩<br>[mas yr<sup>-1</sup>]</th><th>PM†</th><th>Search radius<br>[arcmin]</th><th>nb</th></tr></thead><tbody><tr><td>ASCC 10</td><td>51.765</td><td>35.030</td><td>1.92 ± 0.20</td><td>–</td><td>1.43</td><td>-3.38 ± 0.16</td><td>-1.13 ± 0.20</td><td>U</td><td>97</td><td>21</td></tr><tr><td>ASCC 112</td><td>304.095</td><td>52.098</td><td>1.52 ± 0.27</td><td>–</td><td>1.76</td><td>-0.96 ± 0.31</td><td>-0.45 ± 0.33</td><td>U</td><td>61</td><td>4</td></tr><tr><td>ASCC 113</td><td>318.007</td><td>38.610</td><td>1.74 ± 0.07</td><td>2.22</td><td>2.0</td><td>0.78 ± 0.13</td><td>-3.86 ± 0.09</td><td>T</td><td>74</td><td>35</td></tr><tr><td>ASCC 123</td><td>340.597</td><td>54.250</td><td>4.58 ± 0.16</td><td>4.0</td><td>4.0</td><td>11.42 ± 0.28</td><td>-1.54 ± 0.29</td><td>U</td><td>136</td><td>9</td></tr><tr><td>ASCC 124</td><td>341.895</td><td>46.310</td><td>1.27 ± 0.13</td><td>–</td><td>1.43</td><td>0.36 ± 0.10</td><td>-1.85 ± 0.10</td><td>T</td><td>54</td><td>11</td></tr><tr><td>ASCC 16</td><td>81.037</td><td>1.640</td><td>2.89 ± 0.12</td><td>–</td><td>2.52</td><td>0.92 ± 0.28</td><td>0.56 ± 0.21</td><td>U</td><td>70</td><td>27</td></tr><tr><td>ASCC 18</td><td>81.562</td><td>0.900</td><td>2.89 ± 0.10</td><td>2.0</td><td>3.19</td><td>0.54 ± 0.24</td><td>0.50 ± 0.14</td><td>U</td><td>133</td><td>45</td></tr><tr><td>ASCC 19</td><td>81.922</td><td>-1.965</td><td>2.69 ± 0.12</td><td>2.86</td><td>3.18</td><td>0.97 ± 0.18</td><td>-0.97 ± 0.15</td><td>T</td><td>96</td><td>25</td></tr><tr><td>ASCC 21</td><td>82.260</td><td>3.560</td><td>3.11 ± 0.15</td><td>2.0</td><td>2.64</td><td>0.60 ± 0.23</td><td>0.79 ± 0.28</td><td>U</td><td>66</td><td>13</td></tr><tr><td>ASCC 23</td><td>95.089</td><td>46.675</td><td>1.59 ± 0.14</td><td>1.67</td><td>1.67</td><td>0.23 ± 0.24</td><td>-0.48 ± 0.11</td><td>U</td><td>112</td><td>18</td></tr><tr><td>ASCC 32</td><td>105.525</td><td>-26.560</td><td>1.23 ± 0.08</td><td>–</td><td>1.3</td><td>-3.36 ± 0.10</td><td>3.68 ± 0.08</td><td>T</td><td>89</td><td>33</td></tr><tr><td>ASCC 41</td><td>116.745</td><td>0.020</td><td>3.30 ± 0.17</td><td>–</td><td>2.52</td><td>1.13 ± 0.18</td><td>-4.45 ± 0.63</td><td>U</td><td>89</td><td>12</td></tr><tr><td>ASCC 51</td><td>139.522</td><td>-69.690</td><td>1.99 ± 0.29</td><td>2.0</td><td>2.02</td><td>-8.12 ± 0.35</td><td>8.30 ± 0.29</td><td>T</td><td>66</td><td>6</td></tr><tr><td>ASCC 99</td><td>282.277</td><td>-18.740</td><td>3.37 ± 0.15</td><td>3.57</td><td>1.47</td><td>5.72 ± 0.48</td><td>-0.49 ± 0.33</td><td>U</td><td>119</td><td>15</td></tr><tr><td>Alessi 10</td><td>301.182</td><td>-10.485</td><td>1.92 ± 0.19</td><td>1.95</td><td>1.95</td><td>1.90 ± 0.28</td><td>-7.10 ± 0.37</td><td>T</td><td>63</td><td>11</td></tr><tr><td>Alessi 12</td><td>310.845</td><td>23.785</td><td>1.90 ± 0.13</td><td>1.86</td><td>1.84</td><td>3.66 ± 0.20</td><td>-5.21 ± 0.15</td><td>U</td><td>126</td><td>21</td></tr><tr><td>Alessi 13</td><td>52.020</td><td>-35.870</td><td>9.50 ± 0.44</td><td>10.0</td><td>8.93</td><td>36.43 ± 0.56</td><td>-4.41 ± 0.33</td><td>T</td><td>624</td><td>15</td></tr><tr><td>Alessi 21</td><td>107.670</td><td>-9.300</td><td>1.72 ± 0.09</td><td>2.0</td><td>1.68</td><td>-5.41 ± 0.19</td><td>2.53 ± 0.20</td><td>T</td><td>66</td><td>31</td></tr><tr><td>Alessi 2</td><td>71.565</td><td>55.255</td><td>1.69 ± 0.09</td><td>2.0</td><td>1.92</td><td>-0.86 ± 0.27</td><td>-0.80 ± 0.18</td><td>U</td><td>62</td><td>21</td></tr><tr><td>Alessi 3</td><td>109.050</td><td>-46.617</td><td>3.71 ± 0.08</td><td>3.47</td><td>3.22</td><td>-9.70 ± 0.11</td><td>12.21 ± 0.12</td><td>T</td><td>237</td><td>29</td></tr><tr><td>Alessi 5</td><td>160.830</td><td>-61.170</td><td>2.51 ± 0.04</td><td>2.51</td><td>2.6</td><td>-15.50 ± 0.08</td><td>2.62 ± 0.11</td><td>T</td><td>84</td><td>16</td></tr><tr><td>Alessi 6</td><td>220.132</td><td>-66.120</td><td>1.23 ± 0.03</td><td>–</td><td>1.63</td><td>-10.36 ± 0.09</td><td>-5.67 ± 0.08</td><td>T</td><td>76</td><td>21</td></tr><tr><td>Alessi 9</td><td>266.295</td><td>-47.185</td><td>4.89 ± 0.13</td><td>4.74</td><td>4.59</td><td>11.10 ± 0.33</td><td>-8.21 ± 0.41</td><td>U</td><td>351</td><td>25</td></tr><tr><td>BH 99</td><td>159.465</td><td>-59.125</td><td>2.39 ± 0.15</td><td>1.97</td><td>1.93</td><td>-14.48 ± 0.12</td><td>1.30 ± 0.13</td><td>T</td><td>67</td><td>12</td></tr><tr><td>Chereul 1</td><td>217.267</td><td>55.392</td><td>10.15 ± 0.31</td><td>–</td><td>–</td><td>-17.89 ± 0.32</td><td>-5.01 ± 0.38</td><td>T</td><td>653</td><td>8</td></tr><tr><td>Collinder 135</td><td>109.485</td><td>-37.035</td><td>3.46 ± 0.12</td><td>3.16</td><td>2.77</td><td>-9.28 ± 0.27</td><td>5.57 ± 0.22</td><td>U</td><td>217</td><td>36</td></tr><tr><td>Collinder 350</td><td>267.052</td><td>1.360</td><td>2.71 ± 0.13</td><td>3.57</td><td>3.31</td><td>-3.68 ± 0.35</td><td>0.29 ± 0.29</td><td>U</td><td>242</td><td>19</td></tr><tr><td>Collinder 359</td><td>270.120</td><td>2.875</td><td>1.93 ± 0.10</td><td>4.02</td><td>1.56</td><td>1.98 ± 0.23</td><td>-8.19 ± 0.26</td><td>U</td><td>271</td><td>54</td></tr><tr><td>Collinder 394</td><td>283.080</td><td>-20.200</td><td>1.38 ± 0.10</td><td>1.45</td><td>1.5</td><td>-0.76 ± 0.28</td><td>-4.81 ± 0.32</td><td>U</td><td>45</td><td>32</td></tr><tr><td>Collinder 463</td><td>27.300</td><td>71.780</td><td>1.36 ± 0.07</td><td>1.42</td><td>1.25</td><td>-1.83 ± 0.07</td><td>-0.29 ± 0.15</td><td>T</td><td>96</td><td>41</td></tr><tr><td>IC 4725</td><td>277.927</td><td>-19.110</td><td>1.52 ± 0.12</td><td>1.61</td><td>1.78</td><td>-3.43 ± 0.35</td><td>-5.87 ± 0.33</td><td>U</td><td>50</td><td>44</td></tr><tr><td>IC 4756</td><td>279.735</td><td>5.450</td><td>1.94 ± 0.03</td><td>2.07</td><td>2.07</td><td>0.66 ± 0.13</td><td>-2.58 ± 0.14</td><td>U</td><td>141</td><td>99</td></tr><tr><td>Lynga 2</td><td>216.112</td><td>-61.330</td><td>0.97 ± 0.12</td><td>1.11</td><td>1.07</td><td>-6.76 ± 0.19</td><td>-4.65 ± 0.23</td><td>T</td><td>32</td><td>9</td></tr><tr><td>Melotte 101</td><td>160.567</td><td>-65.095</td><td>0.41 ± 0.18</td><td>0.5</td><td>0.5</td><td>-6.91 ± 0.24</td><td>2.47 ± 0.18</td><td>T</td><td>31</td><td>10</td></tr><tr><td>NGC 0752</td><td>29.272</td><td>37.790</td><td>2.34 ± 0.05</td><td>–</td><td>2.22</td><td>8.19 ± 0.09</td><td>-12.12 ± 0.07</td><td>U</td><td>148</td><td>73</td></tr><tr><td>NGC 1027</td><td>40.665</td><td>61.620</td><td>0.74 ± 0.09</td><td>0.97</td><td>1.05</td><td>-2.91 ± 0.22</td><td>2.11 ± 0.16</td><td>U</td><td>42</td><td>21</td></tr><tr><td>NGC 1039</td><td>40.549</td><td>42.790</td><td>1.94 ± 0.10</td><td>2.0</td><td>1.96</td><td>-0.51 ± 0.17</td><td>-6.39 ± 0.13</td><td>U</td><td>134</td><td>41</td></tr><tr><td>NGC 1342</td><td>53.004</td><td>37.357</td><td>1.47 ± 0.06</td><td>1.5</td><td>1.5</td><td>0.08 ± 0.19</td><td>-1.54 ± 0.08</td><td>T</td><td>49</td><td>22</td></tr><tr><td>NGC 1528</td><td>63.847</td><td>51.190</td><td>0.99 ± 0.09</td><td>0.92</td><td>1.05</td><td>2.08 ± 0.19</td><td>-1.71 ± 0.16</td><td>U</td><td>40</td><td>21</td></tr><tr><td>NGC 1545</td><td>65.295</td><td>50.252</td><td>1.38 ± 0.14</td><td>1.41</td><td>1.16</td><td>-2.26 ± 0.17</td><td>0.45 ± 0.18</td><td>U</td><td>45</td><td>10</td></tr><tr><td>NGC 1647</td><td>71.497</td><td>19.170</td><td>1.88 ± 0.04</td><td>1.85</td><td>1.75</td><td>-0.97 ± 0.09</td><td>-1.66 ± 0.12</td><td>U</td><td>122</td><td>57</td></tr><tr><td>NGC 1662</td><td>72.112</td><td>10.925</td><td>2.56 ± 0.07</td><td>2.29</td><td>2.29</td><td>-1.05 ± 0.14</td><td>-2.02 ± 0.11</td><td>T</td><td>154</td><td>36</td></tr><tr><td>NGC 1750</td><td>76.140</td><td>23.830</td><td>1.48 ± 0.08</td><td>1.59</td><td>–</td><td>-3.39 ± 0.14</td><td>-3.71 ± 0.13</td><td>U</td><td>107</td><td>57</td></tr><tr><td>NGC 1778</td><td>76.995</td><td>37.028</td><td>0.52 ± 0.11</td><td>0.68</td><td>0.67</td><td>-0.37 ± 0.41</td><td>-4.42 ± 0.38</td><td>U</td><td>18</td><td>5</td></tr><tr><td>NGC 1901</td><td>79.560</td><td>-68.440</td><td>2.51 ± 0.11</td><td>2.17</td><td>2.46</td><td>0.78 ± 0.31</td><td>13.11 ± 0.24</td><td>T</td><td>81</td><td>6</td></tr><tr><td>NGC 1912</td><td>82.215</td><td>35.800</td><td>0.88 ± 0.07</td><td>0.71</td><td>0.87</td><td>-0.01 ± 0.23</td><td>-4.51 ± 0.25</td><td>U</td><td>29</td><td>39</td></tr><tr><td>NGC 1960</td><td>84.082</td><td>34.158</td><td>1.05 ± 0.07</td><td>0.75</td><td>0.83</td><td>-0.50 ± 0.17</td><td>-5.29 ± 0.14</td><td>U</td><td>21</td><td>26</td></tr><tr><td>NGC 1977</td><td>83.887</td><td>-4.830</td><td>2.31 ± 0.15</td><td>2.0</td><td>2.79</td><td>-0.25 ± 0.23</td><td>-0.17 ± 0.32</td><td>T</td><td>64</td><td>14</td></tr><tr><td>NGC 2099</td><td>88.087</td><td>32.570</td><td>0.76 ± 0.11</td><td>0.72</td><td>0.71</td><td>2.25 ± 0.24</td><td>-7.12 ± 0.20</td><td>U</td><td>48</td><td>30</td></tr><tr><td>NGC 2168</td><td>92.302</td><td>24.360</td><td>1.12 ± 0.05</td><td>1.1</td><td>1.2</td><td>0.62 ± 0.10</td><td>-4.06 ± 0.08</td><td>U</td><td>83</td><td>103</td></tr><tr><td>NGC 2215</td><td>95.220</td><td>-7.295</td><td>1.31 ± 0.26</td><td>1.27</td><td>0.88</td><td>0.19 ± 0.29</td><td>-3.12 ± 0.34</td><td>T</td><td>52</td><td>16</td></tr><tr><td>NGC 2244</td><td>97.980</td><td>4.940</td><td>0.86 ± 0.13</td><td>0.6</td><td>0.65</td><td>-0.68 ± 0.51</td><td>0.56 ± 0.52</td><td>U</td><td>20</td><td>8</td></tr></tbody></table>

**Notes.** Flag †U indicates that the UCAC4 proper motions were used, T indicates TGAS proper motions. $\omega_{DAML}$ and $\omega_{MWSC}$ are the expected parallaxes assuming no error on the catalogue value.
<table><caption>Table A.1. continued</caption><thead><tr><th>OC</th><th>α<br>[deg]</th><th>δ<br>[deg]</th><th>⟨ω⟩<br>[mas]</th><th>ω<sub>DAML</sub><br>[mas]</th><th>ω<sub>MWSC</sub><br>[mas]</th><th>⟨μ<sub>α</sub>⟩<br>[mas yr<sup>-1</sup>]</th><th>⟨μ<sub>δ</sub>⟩<br>[mas yr<sup>-1</sup>]</th><th>PM<sup>†</sup></th><th>Search radius<br>[arcmin]</th><th>nb</th></tr></thead><tbody><tr><td>NGC 2264</td><td>100.245</td><td>9.880</td><td>1.38 ± 0.11</td><td>1.5</td><td>1.54</td><td>-0.51 ± 0.30</td><td>-3.68 ± 0.28</td><td>U</td><td>48</td><td>6</td></tr><tr><td>NGC 2281</td><td>102.075</td><td>41.080</td><td>2.03 ± 0.08</td><td>1.79</td><td>2.0</td><td>-3.96 ± 0.14</td><td>-8.05 ± 0.14</td><td>U</td><td>121</td><td>35</td></tr><tr><td>NGC 2287</td><td>101.512</td><td>-20.745</td><td>1.29 ± 0.03</td><td>1.41</td><td>1.3</td><td>-4.40 ± 0.06</td><td>-1.39 ± 0.04</td><td>T</td><td>97</td><td>77</td></tr><tr><td>NGC 2323</td><td>105.679</td><td>-8.370</td><td>1.01 ± 0.04</td><td>1.05</td><td>1.11</td><td>-0.69 ± 0.11</td><td>-0.35 ± 0.10</td><td>T</td><td>72</td><td>110</td></tr><tr><td>NGC 2353</td><td>108.627</td><td>-10.250</td><td>0.67 ± 0.16</td><td>0.85</td><td>0.85</td><td>-1.03 ± 0.25</td><td>0.81 ± 0.24</td><td>T</td><td>29</td><td>18</td></tr><tr><td>NGC 2360</td><td>109.425</td><td>-15.640</td><td>0.77 ± 0.04</td><td>0.89</td><td>0.64</td><td>-0.14 ± 0.16</td><td>5.99 ± 0.11</td><td>T</td><td>35</td><td>40</td></tr><tr><td>NGC 2423</td><td>114.277</td><td>-13.850</td><td>1.07 ± 0.06</td><td>1.31</td><td>1.21</td><td>-0.82 ± 0.11</td><td>-3.54 ± 0.09</td><td>T</td><td>78</td><td>78</td></tr><tr><td>NGC 2437</td><td>115.455</td><td>-14.805</td><td>0.66 ± 0.03</td><td>0.66</td><td>0.73</td><td>-3.78 ± 0.08</td><td>0.33 ± 0.06</td><td>T</td><td>49</td><td>120</td></tr><tr><td>NGC 2447</td><td>116.145</td><td>-23.860</td><td>1.13 ± 0.06</td><td>0.96</td><td>0.96</td><td>-3.51 ± 0.15</td><td>5.06 ± 0.06</td><td>T</td><td>65</td><td>62</td></tr><tr><td>NGC 2451B</td><td>116.115</td><td>-37.670</td><td>2.64 ± 0.10</td><td>3.31</td><td>2.04</td><td>-9.70 ± 0.26</td><td>4.33 ± 0.13</td><td>T</td><td>227</td><td>34</td></tr><tr><td>NGC 2477</td><td>118.035</td><td>-38.530</td><td>0.64 ± 0.12</td><td>0.75</td><td>0.69</td><td>-0.70 ± 0.57</td><td>1.94 ± 0.52</td><td>U</td><td>25</td><td>18</td></tr><tr><td>NGC 2482</td><td>118.792</td><td>-24.275</td><td>0.75 ± 0.11</td><td>0.84</td><td>0.74</td><td>-4.47 ± 0.24</td><td>2.27 ± 0.12</td><td>T</td><td>23</td><td>35</td></tr><tr><td>NGC 2527</td><td>121.275</td><td>-28.170</td><td>1.71 ± 0.04</td><td>1.66</td><td>1.56</td><td>-5.38 ± 0.08</td><td>7.37 ± 0.08</td><td>T</td><td>113</td><td>83</td></tr><tr><td>NGC 2539</td><td>122.670</td><td>-12.840</td><td>1.17 ± 0.18</td><td>0.73</td><td>0.81</td><td>-2.48 ± 0.25</td><td>-2.33 ± 0.29</td><td>U</td><td>47</td><td>25</td></tr><tr><td>NGC 2546</td><td>122.962</td><td>-37.610</td><td>0.98 ± 0.07</td><td>1.09</td><td>1.07</td><td>-4.30 ± 0.11</td><td>4.41 ± 0.10</td><td>T</td><td>35</td><td>46</td></tr><tr><td>NGC 2548</td><td>123.435</td><td>-5.770</td><td>1.44 ± 0.06</td><td>1.3</td><td>1.27</td><td>-0.39 ± 0.20</td><td>2.68 ± 0.23</td><td>U</td><td>87</td><td>63</td></tr><tr><td>NGC 2567</td><td>124.642</td><td>-30.650</td><td>0.56 ± 0.06</td><td>0.6</td><td>0.6</td><td>-2.97 ± 0.11</td><td>2.34 ± 0.08</td><td>T</td><td>19</td><td>21</td></tr><tr><td>NGC 2571</td><td>124.735</td><td>-29.735</td><td>1.02 ± 0.14</td><td>0.75</td><td>0.75</td><td>-4.50 ± 0.21</td><td>4.36 ± 0.14</td><td>T</td><td>23</td><td>14</td></tr><tr><td>NGC 2669</td><td>131.617</td><td>-52.940</td><td>1.06 ± 0.16</td><td>0.96</td><td>0.95</td><td>-4.09 ± 0.35</td><td>4.55 ± 0.40</td><td>T</td><td>30</td><td>11</td></tr><tr><td>NGC 2670</td><td>131.392</td><td>-48.820</td><td>0.59 ± 0.09</td><td>0.84</td><td>0.77</td><td>-5.68 ± 0.20</td><td>3.70 ± 0.23</td><td>T</td><td>24</td><td>14</td></tr><tr><td>NGC 2682</td><td>132.847</td><td>11.815</td><td>1.42 ± 0.20</td><td>1.24</td><td>1.12</td><td>-9.38 ± 0.35</td><td>-4.35 ± 0.33</td><td>U</td><td>31</td><td>17</td></tr><tr><td>NGC 3228</td><td>155.355</td><td>-51.730</td><td>2.27 ± 0.11</td><td>1.84</td><td>1.75</td><td>-15.63 ± 0.48</td><td>-0.34 ± 0.38</td><td>T</td><td>125</td><td>19</td></tr><tr><td>NGC 3330</td><td>159.690</td><td>-54.115</td><td>0.65 ± 0.07</td><td>1.12</td><td>0.8</td><td>-7.20 ± 0.13</td><td>0.97 ± 0.14</td><td>T</td><td>36</td><td>21</td></tr><tr><td>NGC 3680</td><td>171.405</td><td>-43.235</td><td>1.18 ± 0.10</td><td>1.07</td><td>1.06</td><td>-7.15 ± 0.22</td><td>0.90 ± 0.18</td><td>T</td><td>34</td><td>17</td></tr><tr><td>NGC 4103</td><td>181.665</td><td>-61.250</td><td>0.38 ± 0.27</td><td>0.61</td><td>0.6</td><td>-6.56 ± 0.31</td><td>-0.60 ± 0.24</td><td>T</td><td>17</td><td>6</td></tr><tr><td>NGC 4609</td><td>190.575</td><td>-62.990</td><td>0.60 ± 0.17</td><td>0.76</td><td>0.76</td><td>-4.97 ± 0.27</td><td>-1.20 ± 0.19</td><td>T</td><td>25</td><td>8</td></tr><tr><td>NGC 4852</td><td>195.045</td><td>-59.590</td><td>0.63 ± 0.14</td><td>0.91</td><td>0.9</td><td>-8.03 ± 0.14</td><td>-2.10 ± 0.19</td><td>T</td><td>29</td><td>21</td></tr><tr><td>NGC 5138</td><td>201.795</td><td>-59.030</td><td>0.62 ± 0.11</td><td>0.5</td><td>0.55</td><td>-3.67 ± 0.10</td><td>-1.41 ± 0.20</td><td>T</td><td>16</td><td>7</td></tr><tr><td>NGC 5316</td><td>208.485</td><td>-61.870</td><td>0.67 ± 0.05</td><td>0.82</td><td>0.83</td><td>-6.34 ± 0.07</td><td>-1.53 ± 0.07</td><td>T</td><td>25</td><td>29</td></tr><tr><td>NGC 5460</td><td>211.890</td><td>-48.330</td><td>1.49 ± 0.18</td><td>1.43</td><td>1.53</td><td>-3.74 ± 0.50</td><td>-0.51 ± 0.30</td><td>U</td><td>48</td><td>19</td></tr><tr><td>NGC 5617</td><td>217.447</td><td>-60.715</td><td>0.38 ± 0.08</td><td>0.5</td><td>0.57</td><td>-5.38 ± 0.17</td><td>-3.20 ± 0.08</td><td>T</td><td>21</td><td>24</td></tr><tr><td>NGC 5662</td><td>218.895</td><td>-56.615</td><td>1.37 ± 0.06</td><td>1.5</td><td>1.6</td><td>-6.54 ± 0.14</td><td>-7.18 ± 0.07</td><td>T</td><td>49</td><td>31</td></tr><tr><td>NGC 5822</td><td>226.117</td><td>-54.390</td><td>1.19 ± 0.04</td><td>1.07</td><td>1.26</td><td>-7.44 ± 0.09</td><td>-5.30 ± 0.05</td><td>T</td><td>74</td><td>96</td></tr><tr><td>NGC 6025</td><td>240.817</td><td>-60.400</td><td>1.21 ± 0.05</td><td>1.32</td><td>1.31</td><td>-3.05 ± 0.12</td><td>-3.00 ± 0.10</td><td>T</td><td>44</td><td>35</td></tr><tr><td>NGC 6067</td><td>243.300</td><td>-54.220</td><td>0.37 ± 0.06</td><td>0.71</td><td>0.56</td><td>-1.98 ± 0.15</td><td>-2.78 ± 0.08</td><td>T</td><td>22</td><td>46</td></tr><tr><td>NGC 6087</td><td>244.710</td><td>-57.940</td><td>0.93 ± 0.05</td><td>1.12</td><td>1.12</td><td>-1.67 ± 0.08</td><td>-2.47 ± 0.08</td><td>T</td><td>36</td><td>15</td></tr><tr><td>NGC 6124</td><td>246.315</td><td>-40.650</td><td>1.66 ± 0.05</td><td>1.95</td><td>1.79</td><td>-0.53 ± 0.20</td><td>-0.53 ± 0.21</td><td>U</td><td>130</td><td>48</td></tr><tr><td>NGC 6134</td><td>246.952</td><td>-49.155</td><td>1.03 ± 0.11</td><td>0.79</td><td>1.12</td><td>0.77 ± 0.38</td><td>-4.18 ± 0.41</td><td>T</td><td>36</td><td>20</td></tr><tr><td>NGC 6152</td><td>248.175</td><td>-52.640</td><td>0.63 ± 0.05</td><td>0.97</td><td>0.97</td><td>-3.14 ± 0.20</td><td>-4.88 ± 0.20</td><td>T</td><td>30</td><td>30</td></tr><tr><td>NGC 6281</td><td>256.170</td><td>-37.980</td><td>1.70 ± 0.08</td><td>2.09</td><td>1.95</td><td>-2.04 ± 0.26</td><td>-2.81 ± 0.32</td><td>U</td><td>70</td><td>57</td></tr><tr><td>NGC 6405</td><td>265.080</td><td>-32.215</td><td>2.00 ± 0.14</td><td>2.05</td><td>2.81</td><td>-1.44 ± 0.31</td><td>-4.76 ± 0.31</td><td>U</td><td>68</td><td>41</td></tr><tr><td>NGC 6416</td><td>266.092</td><td>-32.365</td><td>1.06 ± 0.12</td><td>1.35</td><td>1.3</td><td>-1.09 ± 0.46</td><td>-0.47 ± 0.36</td><td>U</td><td>43</td><td>28</td></tr><tr><td>NGC 6494</td><td>269.227</td><td>-19.000</td><td>1.23 ± 0.08</td><td>1.59</td><td>1.54</td><td>1.17 ± 0.24</td><td>0.09 ± 0.22</td><td>U</td><td>49</td><td>53</td></tr><tr><td>NGC 6604</td><td>274.507</td><td>-12.238</td><td>0.54 ± 0.15</td><td>0.59</td><td>0.53</td><td>-0.72 ± 0.36</td><td>-0.97 ± 0.46</td><td>U</td><td>37</td><td>23</td></tr><tr><td>NGC 6694</td><td>281.317</td><td>-9.380</td><td>0.51 ± 0.05</td><td>0.6</td><td>0.57</td><td>0.28 ± 0.28</td><td>-1.01 ± 0.24</td><td>T</td><td>17</td><td>12</td></tr><tr><td>NGC 6705</td><td>282.747</td><td>-6.280</td><td>0.57 ± 0.13</td><td>0.53</td><td>0.57</td><td>-2.19 ± 0.46</td><td>-5.07 ± 0.48</td><td>T</td><td>15</td><td>7</td></tr><tr><td>NGC 6716</td><td>283.642</td><td>-19.885</td><td>1.38 ± 0.12</td><td>1.27</td><td>1.49</td><td>-0.56 ± 0.32</td><td>-4.14 ± 0.39</td><td>U</td><td>85</td><td>56</td></tr><tr><td>NGC 6793</td><td>290.824</td><td>22.140</td><td>1.66 ± 0.16</td><td>0.91</td><td>1.38</td><td>3.77 ± 0.30</td><td>3.51 ± 0.26</td><td>T</td><td>60</td><td>27</td></tr><tr><td>NGC 6811</td><td>294.342</td><td>46.395</td><td>0.98 ± 0.11</td><td>0.82</td><td>0.81</td><td>-4.60 ± 0.18</td><td>-7.69 ± 0.17</td><td>U</td><td>54</td><td>27</td></tr><tr><td>NGC 6866</td><td>300.987</td><td>44.160</td><td>0.96 ± 0.25</td><td>0.84</td><td>0.75</td><td>-1.90 ± 0.63</td><td>-4.97 ± 0.27</td><td>U</td><td>9</td><td>4</td></tr><tr><td>NGC 6913</td><td>306.000</td><td>38.510</td><td>0.52 ± 0.06</td><td>0.65</td><td>0.85</td><td>-3.97 ± 0.13</td><td>-5.51 ± 0.13</td><td>U</td><td>58</td><td>51</td></tr><tr><td>NGC 6940</td><td>308.610</td><td>28.280</td><td>0.93 ± 0.05</td><td>1.3</td><td>1.18</td><td>-2.10 ± 0.11</td><td>-9.45 ± 0.12</td><td>T</td><td>87</td><td>74</td></tr><tr><td>NGC 6991</td><td>313.614</td><td>47.466</td><td>1.85 ± 0.09</td><td>1.43</td><td>1.77</td><td>5.22 ± 0.16</td><td>8.79 ± 0.08</td><td>U</td><td>97</td><td>42</td></tr><tr><td>NGC 7209</td><td>331.225</td><td>46.552</td><td>1.14 ± 0.13</td><td>0.86</td><td>0.8</td><td>2.05 ± 0.19</td><td>0.38 ± 0.14</td><td>T</td><td>57</td><td>28</td></tr><tr><td>NGC 7243</td><td>333.795</td><td>49.875</td><td>1.30 ± 0.11</td><td>1.24</td><td>1.18</td><td>0.10 ± 0.10</td><td>-2.33 ± 0.09</td><td>U</td><td>40</td><td>29</td></tr><tr><td>Platais 10</td><td>205.290</td><td>-59.225</td><td>3.84 ± 0.18</td><td>4.07</td><td>4.07</td><td>-30.22 ± 0.30</td><td>-10.59 ± 0.44</td><td>T</td><td>138</td><td>10</td></tr><tr><td>Platais 3</td><td>69.450</td><td>71.470</td><td>5.60 ± 0.05</td><td>5.0</td><td>5.88</td><td>3.51 ± 0.26</td><td>-21.10 ± 0.20</td><td>T</td><td>343</td><td>15</td></tr><tr><td>Platais 8</td><td>136.875</td><td>-59.160</td><td>7.59 ± 0.11</td><td>7.58</td><td>7.09</td><td>-15.82 ± 0.29</td><td>14.93 ± 0.32</td><td>T</td><td>457</td><td>19</td></tr></tbody></table>
<table><caption>**Table A.1.** continued</caption><thead><tr><th>OC</th><th>α<br>[deg]</th><th>δ<br>[deg]</th><th>⟨ω⟩<br>[mas]</th><th>ω<sub>DAML</sub><br>[mas]</th><th>ω<sub>MWSC</sub><br>[mas]</th><th>⟨μ<sub>α</sub>⟩<br>[mas yr<sup>-1</sup>]</th><th>⟨μ<sub>δ</sub>⟩<br>[mas yr<sup>-1</sup>]</th><th>PM<sup>†</sup></th><th>Search radius<br>[arcmin]</th><th>nb</th></tr></thead><tbody><tr><td>Platais 9</td><td>137.955</td><td>-43.530</td><td>5.59 ± 0.26</td><td>5.75</td><td>5.0</td><td>-25.47 ± 0.38</td><td>13.54 ± 0.30</td><td>T</td><td>342</td><td>21</td></tr><tr><td>Roslund 3</td><td>299.662</td><td>20.514</td><td>0.77 ± 0.15</td><td>0.68</td><td>0.63</td><td>-1.36 ± 0.32</td><td>-4.18 ± 0.26</td><td>U</td><td>19</td><td>14</td></tr><tr><td>Roslund 6</td><td>307.140</td><td>39.205</td><td>2.68 ± 0.10</td><td>2.22</td><td>1.84</td><td>5.18 ± 0.21</td><td>1.75 ± 0.20</td><td>U</td><td>152</td><td>41</td></tr><tr><td>Ruprecht 145</td><td>282.660</td><td>-18.227</td><td>1.47 ± 0.08</td><td>3.13</td><td>1.3</td><td>7.51 ± 0.33</td><td>-2.60 ± 0.23</td><td>U</td><td>85</td><td>30</td></tr><tr><td>Ruprecht 147</td><td>289.092</td><td>-16.250</td><td>3.26 ± 0.09</td><td>3.39</td><td>3.7</td><td>-1.48 ± 0.26</td><td>-26.85 ± 0.20</td><td>U</td><td>194</td><td>63</td></tr><tr><td>Ruprecht 1</td><td>99.135</td><td>-14.193</td><td>0.87 ± 0.15</td><td>0.68</td><td>0.83</td><td>-1.26 ± 0.30</td><td>0.16 ± 0.18</td><td>T</td><td>28</td><td>16</td></tr><tr><td>Ruprecht 98</td><td>179.715</td><td>-64.588</td><td>2.08 ± 0.10</td><td>2.02</td><td>1.64</td><td>-4.30 ± 0.16</td><td>-8.66 ± 0.13</td><td>T</td><td>138</td><td>20</td></tr><tr><td>Stock 10</td><td>84.810</td><td>37.805</td><td>2.90 ± 0.20</td><td>2.63</td><td>2.0</td><td>-4.16 ± 0.27</td><td>-1.08 ± 0.29</td><td>U</td><td>87</td><td>13</td></tr><tr><td>Stock 12</td><td>353.895</td><td>52.685</td><td>2.17 ± 0.10</td><td>2.08</td><td>2.33</td><td>8.22 ± 0.23</td><td>-2.32 ± 0.13</td><td>U</td><td>84</td><td>22</td></tr><tr><td>Stock 1</td><td>293.955</td><td>25.175</td><td>2.31 ± 0.05</td><td>3.14</td><td>2.86</td><td>5.92 ± 0.09</td><td>0.49 ± 0.12</td><td>T</td><td>213</td><td>50</td></tr><tr><td>Stock 2</td><td>33.660</td><td>59.440</td><td>2.78 ± 0.03</td><td>3.3</td><td>2.5</td><td>16.26 ± 0.06</td><td>-13.79 ± 0.05</td><td>U</td><td>226</td><td>149</td></tr><tr><td>Stock 7</td><td>37.432</td><td>60.675</td><td>1.44 ± 0.26</td><td>1.43</td><td>1.18</td><td>-4.57 ± 0.19</td><td>1.11 ± 0.16</td><td>U</td><td>48</td><td>12</td></tr><tr><td>Trumpler 10</td><td>131.872</td><td>-42.400</td><td>2.28 ± 0.07</td><td>2.36</td><td>2.4</td><td>-12.38 ± 0.13</td><td>6.84 ± 0.12</td><td>T</td><td>161</td><td>59</td></tr><tr><td>Trumpler 2</td><td>39.367</td><td>55.950</td><td>1.49 ± 0.10</td><td>1.38</td><td>1.49</td><td>0.14 ± 0.16</td><td>-6.28 ± 0.13</td><td>U</td><td>50</td><td>14</td></tr><tr><td>Trumpler 33</td><td>276.172</td><td>-19.720</td><td>0.59 ± 0.13</td><td>0.57</td><td>0.78</td><td>-0.81 ± 0.76</td><td>-1.01 ± 1.41</td><td>U</td><td>10</td><td>7</td></tr><tr><td>Trumpler 3</td><td>48.052</td><td>63.150</td><td>1.36 ± 0.10</td><td>2.17</td><td>1.61</td><td>-3.90 ± 0.18</td><td>-0.25 ± 0.10</td><td>U</td><td>46</td><td>20</td></tr><tr><td>Turner 5</td><td>143.317</td><td>-36.410</td><td>2.49 ± 0.08</td><td>2.5</td><td>2.5</td><td>0.19 ± 0.29</td><td>-2.67 ± 0.24</td><td>T</td><td>83</td><td>11</td></tr><tr><td>vdBergh 92</td><td>106.020</td><td>-11.530</td><td>0.86 ± 0.21</td><td>0.67</td><td>0.64</td><td>-5.28 ± 0.41</td><td>2.29 ± 0.30</td><td>T</td><td>20</td><td>10</td></tr></tbody></table>