MNRAS 000 , 1-11 (2020)

## Kinematic study of the association Cyg OB3 with Gaia DR2

Anjali Rao 1 ? , Poshak Gandhi 1 ; Christian Knigge 1 ; John A. Paice 1 2 ; ; Nathan W. C. Leigh 3 4 ; ; and Douglas Boubert 5

- 1 Department of Physics &amp; Astronomy, University of Southampton, Highfield, Southampton SO17 1BJ, UK
- 2 Inter-University Centre for Astronomy and Astrophysics, Pune, Maharashtra 411007, India
- 3 Departamento de Astronom´ ıa, Facultad de Ciencias F´ ısicas y Matem´ticas, Universidad de Concepci´n, Concepci´n, Chile a o o
- 4 Department of Astrophysics, American Museum of Natural History, New York, NY 10024, USA
- 5 Magdalen College, University of Oxford, High Street, Oxford OX1 4AU, UK

Accepted XXX. Received YYY; in original form ZZZ

## ABSTRACT

We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 GLYPH&lt;6&gt; 0.3 kpc. The mean peculiar velocity of the sample after removing Galactic rotation and solar motion is GLYPH&lt;24&gt; 22 km s GLYPH&lt;0&gt; 1 , dominated by the velocity component towards the Galactic center. The relative position and velocity of the black hole X-ray binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most likely its parent association. The peculiar kinematic properties of some of the stars are revealed and are suggestive of past stellar encounters. The sample includes a previously known runaway star HD 227018, and its high peculiar velocity of GLYPH&lt;24&gt; 50 km s GLYPH&lt;0&gt; 1 is confirmed with Gaia . We estimated the velocities of stars relative to the association and the star HD 225577 exhibits peculiar velocity smaller than its velocity relative to the association. The star has lower value of proper motion than the rest of the sample. The results suggest a slowly expanding nature of the association, which is supported by the small relative speeds &lt; 20 km s GLYPH&lt;0&gt; 1 with respect to the association for a majority of the sample stars.

Key words: stars: distances - stars: kinematics and dynamics - parallaxes - proper motions - open clusters and associations: individual: Cyg OB3 - stars: individual: Cyg X-1

## 1 INTRODUCTION

OB associations are regions in the sky with higher densities of O- and B-type stars relative to the field (Ambartsumian 1947). These systems are identified with sites of the most recent star formation regions in the Galaxy and usually host the youngest and most massive stars (see review by Blaauw 1964, and references therein). It has long been known that associations cannot be bound by their own self-gravity and tend to form unbound stellar aggregates spanning tens of parsecs (Ambartsumian 1947, 1955; Mel'nik &amp; Dambis 2017). Many OB associations are known in the Galaxy and catalogued in papers such as Humphreys (1978), Ruprecht et al. (1981), Humphreys &amp; McElroy

(1984), Blaha &amp; Humphreys (1989), Garmany &amp; Stencel (1992), and Mel'Nik &amp; Efremov (1995). OB associations have historically been of interest for several reasons, for instance in the investigation of the motions of loosely bound or unbound systems in the Galactic potential (e.g. Lindblad 1941, 1942), the dynamical interactions of the member stars (e.g Wright &amp; Mamajek 2018), and the formation of stars in giant molecular clouds (Clark et al. 2005). Recently, Gaia DR2 has significantly contributed to the study of OB associations e.g. Cantat-Gaudin et al. (2019), Berlanas et al. (2019), Wright et al. (2019) and Melnik &amp; Dambis (2020).

? E-mail: a.rao@soton.ac.uk

The membership of associations has been assigned and revised by a number of authors, as historical observational advances have been made. Humphreys (1978) provided a list of stars in known associations, but only included stars for which a measurement of the photometric distance was

Figure 1. Image obtained with DSS2 in the R -band shows the Cyg OB3 region. The positions of the O- and B-type stars studied in this work are encircled. Cyg X-1 is encircled in red. There are a total of 41 stars shown in the image and given in Table 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-ea8e5cfc-2011-4dc2-a73a-7abaab456b32.jpg)

possible based on the information from MK spectral types and UBV photometry. Ruprecht et al. (1981) also presented a catalogue of star clusters and associations. Humphreys &amp; McElroy (1984) expanded the catalogue of Humphreys (1978), and it was further updated by Blaha &amp; Humphreys (1989), where the authors catalogued a total of 2263 stars in 91 associations and clusters and another set of 2549 stars in the field. Later, Garmany &amp; Stencel (1992) also provided a list of member stars in associations and computed their distances using the method of cluster fitting of B-type main sequence stars. In the classification scheme discussed in Humphreys (1978), Humphreys &amp; McElroy (1984), Blaha &amp; Humphreys (1989) and Garmany &amp; Stencel (1992), membership assignment is based on a combination of position on the sky, photometric distance and radial velocity.

Mel'Nik &amp; Efremov (1995) adopted a different method, compared to the previous authors, to determine the minimum unitary structures among the stellar associations by implementing an algorithm for cluster analysis. They found the groups and their members and provided a list of new associations in the Galaxy. Clustering depends on the choice of length scale and groups begin to merge for larger length scales. In this scenario, although the presence of different groups and associations is ascertained from the results of their cluster analysis using an independent method, their boundaries are not well defined. This, consequently, affects the membership assignment and identification of parent associations for those stars located at the association boundaries.

Here, we investigate the properties of the well known association Cyg OB3. We primarily focus on the distances, proper motions and three-dimensional velocity components of the members. We study a particular case of a confirmed black hole X-ray binary Cyg X-1 (otherwise known as HD 226868, Bolton 1972; Webster &amp; Murdin 1972) located at the association boundary (Fig. 1). Cyg X-1 is one of the brightest X-ray binaries in the X-ray sky and is composed of a black hole of 14.8 GLYPH&lt;6&gt; 1.0 M GLYPH&lt;12&gt; (Orosz et al. 2011) and a companion of 19.2 GLYPH&lt;6&gt; 1.9 M GLYPH&lt;12&gt; (Orosz et al. 2011) of spectral type O9.7Iab (Walborn 1973). Cyg X-1 has been included in the catalogs of Humphreys (1978), Garmany &amp; Stencel (1992), Ruprecht et al. (1981) and Blaha &amp; Humphreys (1989) as a member of Cyg OB3. The cluster analysis algorithm implemented by Mel'Nik &amp; Efremov (1995) divided the association Cyg OB3 into two subgroups for a length scale of 15 pc viz. Cyg 3A and Cyg 3B, however Cyg X-1 was not included in either of the subgroups and remained as a field star.

Mirabel &amp; Rodrigues (2003) studied the connection between Cyg X-1 and the association Cyg OB3 using the astrometry from the Hipparcos mission. The authors suggested that the parent association of Cyg X-1 is indeed Cyg OB3, and showed that the binary is located at the same distance as the other stars in Cyg OB3 and exhibits a proper motion similar to that of the association. The authors also showed that the peculiar velocity of Cyg X-1 is very small (9 GLYPH&lt;6&gt; 2 km s GLYPH&lt;0&gt; 1 ) with respect to the association, and is consistent with the expected random velocities of stars in expanding associations. These results indicate that Cyg X-1 is most likely a member of Cyg OB3. This was later confirmed with precise measurement of the trigonometric parallax and the proper motion of Cyg X-1 using VLBI radio observations (Reid et al. 2011). In this paper, we revisit the proposed connection between Cyg X-1 and Cyg OB3 using the astrometric data from Gaia DR2, which provides trigonometric parallaxes and proper motions for over 1 billion stars in the Galaxy with the best accuracy and precision available to date (Gaia Collaboration et al. 2018).

We define the sample of stars and their filtering criteria in the next section. The results are discussed in detail in Section 3. Our study of proper motions is discussed in Section 4. Error analysis reveals that the parallax uncertainties are the most significant contributor to the uncertainty when determining the three-dimensional distances between stars (see Section 5). In addition, we find the peculiar velocities and relative velocities of individual stars in the sample, where radial velocity measurements were available from the literature, with respect to the median velocity of the association. A comparison of the two velocities turns out to be an efficient method to filter out stars with peculiar kinematic behavior from the main body of the association. The results are discussed in Section 7 and individual cases are discussed in detail.

## 2 OBSERVATIONS

## 2.1 Sample of Stars

The list of associations published by Humphreys (1978) was restricted to only those stars for which photometric distances were available and included a total of 25 O- and B-type stars in Cyg OB3. The list provided by Garmany &amp; Stencel (1992) includes a larger sample of 37 stars in the association. According to the lists by Humphreys (1978) and Garmany &amp; Stencel (1992), the association Cyg OB3 is located between Cyg OB1 and Cyg OB5 within Galactic coordinates of 71 GLYPH&lt;14&gt; GLYPH&lt;20&gt; l GLYPH&lt;20&gt; 74 , and 0.9 GLYPH&lt;14&gt; GLYPH&lt;14&gt; GLYPH&lt;20&gt; b GLYPH&lt;20&gt; 3.5 GLYPH&lt;14&gt; . Instead of restricting to the list of stars given in the above

catalogues, we have searched the entire region within the above mentioned Galactic coordinates, and queried the Gaia DR2 archive for sources with 5-parameter astrometric solutions available with position, parallax and proper motion ( astrometric\_params\_solved = 31 ). Since the faintest OB star in the list of Garmany &amp; Stencel (1992) is of magnitude V =10.32, stars fainter than G =10.30 magnitude are excluded in the query. The relation between V band and Gaia's G band is obtained from Gaia data release documentation 1 . This photometric selection restricts the sample to contain stars of early B-type, in fact all the stars in Garmany &amp; Stencel (1992) are of spectral type B5 or earlier. Cyg OB3 adjoins Cyg OB1 and there is no clear boundary between the two associations. Therefore the sample of stars in Cyg OB3 is expected to be contaminated by the stars of Cyg OB1 and vice versa. In order to reduce the contamination of stars from Cyg OB1, the query rejects the sources located at a distance, obtained from parallax inversion, smaller than 1.4 kpc. This distance threshold is chosen because it is the average distance for Cyg OB1 reported in Mel'nik &amp; Dambis (2017). We also removed all the stars from the sample at distances glyph[greaterorsimilar] 3 kpc to remove distant non-members. Since the sample has not been filtered based on proper motion or other kinematic properties, the sample may represent more than one kinematic substructure as has been observed in other OB associations e.g., Cantat-Gaudin et al. (2019).

The resultant list of sources from the query is crossmatched against SIMBAD using a 2 arc second search radius, and the spectral types of matched stars are extracted. A total of 6 stars were removed where more than one potential cross-match was found in order to avoid source confusion while assigning SIMBAD objects to Gaia counterparts. A comparison of magnitudes was not helpful in the identification of sources in these six cases. Four of these systems are binary. Since astrometric solutions of binary systems are not reliable, these were removed from the sample. Stars with spectral type O and B are shortlisted in the sample. All the stars are bright with G -band magnitudes in the range 7.1910.30 with a median value of 9.16.

## 2.2 Filtering Criteria

We followed the three filtering criteria recommended by the Gaia team described in Lindegren et al. (2018), and also in equations (1-3) of Arenou et al. (2018), to extract good astrometric solutions. We also filtered the solutions by re-normalized unit weight error (RUWE) 2 , defined as u norm = u u GLYPH&lt;157&gt; 0 ' G C ; ' , where u = (astromet-ric\_chi2\_al/astrometric\_n\_good\_obs\_al -5) 1 2 GLYPH&lt;157&gt; . RUWE is a quality metric introduced after the release of DR2, based on analysis of reliability of the astrometric solutions by the Gaia team. Solutions with RUWE GLYPH&lt;21&gt; 1.4 were rejected. Finally, another star HD 227573 was removed from the sample, as it has a fractional uncertainty on its parallax &gt; 1. A number of flags viz. astrometric\_gof\_al , astromet-ric\_chi2\_al , astrometric\_excess\_noise and astromet- ric\_excess\_noise\_sig often used to interpret the goodness of an astrometric fit suggested a bad fit for this star. A total of 4 stars were discarded by these filters. There are 22 stars in common between the list provided by Garmany &amp; Stencel (1992) and our sample. Table 1 gives a list of all these stars. The names of the sources obtained from SIMBAD and their Gaia DR2 source IDs are provided in the second and third columns, respectively. The remaining columns list the DR2 values including RA and Dec coordinates, G -band magnitude, parallax, and proper motion. According to the spectral types obtained from SIMBAD , about two-third of the stars in the sample belong to the B GLYPH&lt;0&gt; spectral type.

## 3 DISTANCE MEASUREMENTS

Heliocentric distances can be estimated by inverting the parallaxes provided in Table 1 for individual stars ( r inv = 1/ $ ). Bailer-Jones (2015) and Astraatmadja &amp; Bailer-Jones (2016) emphasized that parallax inversion may not serve as a good estimator to find the distance and its uncertainty interval if the fractional uncertainty in parallax measurement is large ( glyph[greaterorsimilar] 0.2). It is observed that the fractional uncertainty is less than 0.2 for most of our sample (Table 1).

We estimated the distances to individual stars with a Bayesian prior that considers an exponential decrease of space density as in Bailer-Jones (2015) with a length scale of 2 kpc corresponding to the approximate mean distance to Cyg OB3 and compared with the distances obtained from parallax inversion. The two distance estimates are consistent with each other within 2 GLYPH&lt;27&gt; level, which was expected considering the smaller fractional uncertainties in the parallax measurements. We compared the distance estimates from exponential prior with a length scale of 1.35 kpc used in Astraatmadja &amp; Bailer-Jones (2016), which were found to be consistent with those obtained with a length scale of 2 kpc. We also searched for our sources in the catalogue provided by Anders et al. (2019), where the authors estimate accurate distances to a good number of sources using parallaxes from Gaia DR2 and photometric catalogues of Pan-STARRS1, 2MASS, and AllWISE. The distances are found to be consistent with those obtained from Anders et al. (2019) within uncertainties. Distances are provided in Table 2.

A histogram of the distances ( r exp ) is shown in Fig 2; the majority of stars are located at a distance of about 2 kpc, suggesting a higher density of OB-type stars in a region spanning 0.31 kpc equal to the standard deviation of r exp estimates. In order to compare the observed width of the parallax distribution and the expected width entirely due to parallax uncertainties, we generated 500 realizations of parallax distributions for 41 sources by drawing random values from a normal distribution defined by a mean &lt; $ &gt; = 0.5 mas with the respective GLYPH&lt;27&gt; $ . The observed parallax distribution is found to be significantly wider than the ones obtained from simulations, and suggests that the observed scatter of sources along the line of sight is real and is not entirely driven by parallax uncertainties.

1 https://gea.esac.esa.int/archive/documentation/GDR2/Data processing/ chap cu5pho/sec cu5pho calibr/ssec cu5pho PhotTransf.html 2 See the technical note GAIA-C3-TN-LU-LL-124-01 at https://www.cosmos.esa.int/web/gaia/dr2-known-issues

Histograms were once again simulated with parallax uncertainties corresponding to the size of the association (assumed to be 100 pc for a distance of 2 kpc) and the intrinsic variance defined as GLYPH&lt;27&gt; 2 int = GLYPH&lt;27&gt; 2 obs GLYPH&lt;0&gt; GLYPH&lt;27&gt; 2 sim was estimated, where GLYPH&lt;27&gt; obs and GLYPH&lt;27&gt; sim are, respectively, the widths of the observed

Table 1. Gaia data for the 41 stars in our sample. (a) G denotes the G-band magnitudes of the stars obtained from Gaia DR2. (b) pmra and pmdec represent the proper motions of the stars in RA and Dec.

|   Nr. | Source          |   Source Id ( Gaia DR2) |   RA GLYPH<11> (deg) |   Dec. GLYPH<14> (deg) |   G a Gmag | Parallax $ (mas)   |   f GLYPH<27> $ GLYPH<157> $ | pmra b GLYPH<22> GLYPH<11> cos GLYPH<14> (mas y GLYPH<0> 1 )   | pmdec b GLYPH<22> GLYPH<14> (mas y GLYPH<0> 1 )   |
|-------|-----------------|-------------------------|----------------------|------------------------|------------|--------------------|------------------------------|----------------------------------------------------------------|---------------------------------------------------|
|     1 | HD 227415       |     2058706399118040960 |              300.97  |                35.5661 |       9.39 | 0.57 GLYPH<6> 0.03 |                         0.05 | -3.85 GLYPH<6> 0.04                                            | -6.53 GLYPH<6> 0.05                               |
|     2 | HD 227943       |     2058738972149344896 |              302.267 |                34.7291 |       9.82 | 0.65 GLYPH<6> 0.04 |                         0.06 | -0.73 GLYPH<6> 0.06                                            | -2.66 GLYPH<6> 0.07                               |
|     3 | HD 228022       |     2058830510807237248 |              302.46  |                35.4299 |      10.1  | 0.59 GLYPH<6> 0.04 |                         0.07 | -3.03 GLYPH<6> 0.05                                            | -7.12 GLYPH<6> 0.06                               |
|     4 | HD 228041       |     2058842807280883584 |              302.493 |                35.496  |       8.9  | 0.48 GLYPH<6> 0.04 |                         0.08 | -1.43 GLYPH<6> 0.06                                            | -2.81 GLYPH<6> 0.05                               |
|     5 | HD 191495       |     2058845629092327040 |              302.223 |                35.5129 |       8.35 | 0.59 GLYPH<6> 0.05 |                         0.09 | -3.12 GLYPH<6> 0.07                                            | -7.08 GLYPH<6> 0.09                               |
|     6 | HD 191567       |     2058844903225261568 |              302.308 |                35.4847 |       8.64 | 0.63 GLYPH<6> 0.04 |                         0.07 | -4.21 GLYPH<6> 0.06                                            | -6.65 GLYPH<6> 0.06                               |
|     7 | HD 228104       |     2058950628143309312 |              302.66  |                35.8744 |       8.96 | 0.51 GLYPH<6> 0.03 |                         0.06 | -1.19 GLYPH<6> 0.05                                            | -3.15 GLYPH<6> 0.05                               |
|     8 | HD 191917       |     2058952728407702272 |              302.738 |                35.9531 |       7.71 | 0.44 GLYPH<6> 0.04 |                         0.08 | -3.41 GLYPH<6> 0.06                                            | -6.51 GLYPH<6> 0.06                               |
|     9 | HD 228282       |     2058988217688303488 |              303.087 |                36.3469 |      10.28 | 0.60 GLYPH<6> 0.06 |                         0.1  | -3.30 GLYPH<6> 0.11                                            | -5.81 GLYPH<6> 0.08                               |
|    10 | HD 190967       |     2059011376158159488 |              301.541 |                35.386  |       7.92 | 0.49 GLYPH<6> 0.03 |                         0.07 | -3.44 GLYPH<6> 0.05                                            | -6.43 GLYPH<6> 0.05                               |
|    11 | HD 227722       |     2059003168499159168 |              301.724 |                35.3061 |       9.53 | 0.59 GLYPH<6> 0.03 |                         0.05 | -3.52 GLYPH<6> 0.05                                            | -6.84 GLYPH<6> 0.05                               |
|    12 | BD+35 3976      |     2059033370710817152 |              301.986 |                35.6615 |       9.87 | 0.58 GLYPH<6> 0.03 |                         0.05 | -2.97 GLYPH<6> 0.04                                            | -6.19 GLYPH<6> 0.05                               |
|    13 | HD 190919       |     2059072197216667776 |              301.484 |                35.6721 |       7.19 | 0.43 GLYPH<6> 0.03 |                         0.07 | -3.46 GLYPH<6> 0.05                                            | -7.00 GLYPH<6> 0.05                               |
|    14 | HD 227634       |     2059075873709364864 |              301.506 |                35.7655 |       7.85 | 0.54 GLYPH<6> 0.03 |                         0.06 | -3.15 GLYPH<6> 0.05                                            | -6.56 GLYPH<6> 0.05                               |
|    15 | LS II +35 32    |     2059076148587293696 |              301.485 |                35.7899 |       9.66 | 0.59 GLYPH<6> 0.04 |                         0.06 | -3.09 GLYPH<6> 0.05                                            | -6.45 GLYPH<6> 0.06                               |
|    16 | BD+35 3955      |     2059076251666505728 |              301.495 |                35.7972 |       7.28 | 0.53 GLYPH<6> 0.03 |                         0.06 | -3.14 GLYPH<6> 0.05                                            | -6.29 GLYPH<6> 0.05                               |
|    17 | HD 227586       |     2059073159271061632 |              301.401 |                35.6232 |       8.76 | 0.38 GLYPH<6> 0.04 |                         0.11 | -3.14 GLYPH<6> 0.07                                            | -6.37 GLYPH<6> 0.07                               |
|    18 | HD 227621       |     2059075255233453824 |              301.469 |                35.7068 |      10.24 | 0.54 GLYPH<6> 0.03 |                         0.05 | -3.02 GLYPH<6> 0.04                                            | -6.33 GLYPH<6> 0.04                               |
|    19 | BD+35 3956      |     2059075839349023104 |              301.5   |                35.7623 |       8.82 | 0.54 GLYPH<6> 0.03 |                         0.06 | -2.93 GLYPH<6> 0.05                                            | -6.37 GLYPH<6> 0.06                               |
|    20 | HD 190864       |     2059070135632404992 |              301.416 |                35.6077 |       7.68 | 0.49 GLYPH<6> 0.03 |                         0.06 | -3.07 GLYPH<6> 0.05                                            | -6.49 GLYPH<6> 0.05                               |
|    21 | HD 227696       |     2059095424401407232 |              301.646 |                35.7406 |       8.21 | 0.66 GLYPH<6> 0.04 |                         0.06 | -3.21 GLYPH<6> 0.06                                            | -6.42 GLYPH<6> 0.06                               |
|    22 | HD 227680       |     2059129303105142144 |              301.607 |                36.3291 |       9.59 | 0.39 GLYPH<6> 0.03 |                         0.08 | -2.75 GLYPH<6> 0.05                                            | -5.71 GLYPH<6> 0.05                               |
|    23 | HD 227611       |     2059112879149465600 |              301.438 |                35.9008 |       8.55 | 0.48 GLYPH<6> 0.03 |                         0.06 | -3.15 GLYPH<6> 0.05                                            | -6.54 GLYPH<6> 0.06                               |
|    24 | HD 191612       |     2059130368252069888 |              302.369 |                35.7337 |       7.7  | 0.47 GLYPH<6> 0.04 |                         0.08 | -3.59 GLYPH<6> 0.05                                            | -5.76 GLYPH<6> 0.06                               |
|    25 | HD 227883       |     2059154007751836672 |              302.116 |                36.094  |      10.02 | 0.46 GLYPH<6> 0.03 |                         0.06 | -2.98 GLYPH<6> 0.05                                            | -5.92 GLYPH<6> 0.04                               |
|    26 | HD 227960       |     2059150640497117952 |              302.299 |                36.0485 |       9.38 | 0.47 GLYPH<6> 0.03 |                         0.07 | -2.78 GLYPH<6> 0.05                                            | -5.99 GLYPH<6> 0.05                               |
|    27 | HD 191611       |     2059236196250413696 |              302.359 |                36.4888 |       8.47 | 0.37 GLYPH<6> 0.04 |                         0.1  | -2.57 GLYPH<6> 0.05                                            | -4.98 GLYPH<6> 0.06                               |
|    28 | HD 227757       |     2059219875373718016 |              301.803 |                36.3593 |       9.16 | 0.54 GLYPH<6> 0.03 |                         0.05 | 0.37 GLYPH<6> 0.04                                             | -3.63 GLYPH<6> 0.04                               |
|    29 | HD 191139       |     2059223002094535168 |              301.74  |                36.3966 |       7.9  | 0.46 GLYPH<6> 0.03 |                         0.07 | -3.67 GLYPH<6> 0.05                                            | -6.50 GLYPH<6> 0.06                               |
|    30 | BD+36 3882      |     2059225819608656768 |              301.855 |                36.5545 |       9.8  | 0.54 GLYPH<6> 0.02 |                         0.05 | -3.02 GLYPH<6> 0.04                                            | -5.77 GLYPH<6> 0.04                               |
|    31 | HD 226868       |     2059383668236814720 |              299.59  |                35.2016 |       8.52 | 0.42 GLYPH<6> 0.03 |                         0.08 | -3.88 GLYPH<6> 0.05                                            | -6.17 GLYPH<6> 0.05                               |
|    32 | HD 227245       |     2059455613211582080 |              300.59  |                35.6749 |       9.47 | 0.63 GLYPH<6> 0.03 |                         0.04 | -3.96 GLYPH<6> 0.04                                            | -6.42 GLYPH<6> 0.05                               |
|    33 | HD 227018       |     2059434898612310656 |              299.955 |                35.3093 |       8.86 | 0.47 GLYPH<6> 0.03 |                         0.07 | -5.59 GLYPH<6> 0.05                                            | -7.99 GLYPH<6> 0.05                               |
|    34 | LS II +35 21    |     2059458778631320064 |              300.869 |                35.6928 |      10.3  | 0.64 GLYPH<6> 0.02 |                         0.04 | -3.54 GLYPH<6> 0.04                                            | -6.34 GLYPH<6> 0.04                               |
|    35 | HD 227132       |     2059524543174097664 |              300.312 |                35.9566 |      10.02 | 0.40 GLYPH<6> 0.02 |                         0.06 | -2.47 GLYPH<6> 0.04                                            | -5.29 GLYPH<6> 0.04                               |
|    36 | HD 226951       |     2059706203121650432 |              299.798 |                36.1153 |       9.08 | 0.50 GLYPH<6> 0.03 |                         0.06 | -4.22 GLYPH<6> 0.05                                            | -6.11 GLYPH<6> 0.05                               |
|    37 | HD 227070       |     2059744720391006720 |              300.118 |                36.5328 |      10    | 0.40 GLYPH<6> 0.03 |                         0.08 | -0.83 GLYPH<6> 0.05                                            | -6.79 GLYPH<6> 0.06                               |
|    38 | TYC 2682-2524-1 |     2059878001807789952 |              301.482 |                36.2742 |      10.1  | 0.43 GLYPH<6> 0.03 |                         0.07 | -3.07 GLYPH<6> 0.05                                            | -5.95 GLYPH<6> 0.05                               |
|    39 | BD+35 3929      |     2059854121788578432 |              300.865 |                36.035  |       9.39 | 0.54 GLYPH<6> 0.03 |                         0.06 | -2.88 GLYPH<6> 0.05                                            | -5.87 GLYPH<6> 0.05                               |
|    40 | HD 227460       |     2059863536357618816 |              301.068 |                36.2654 |       9.46 | 0.56 GLYPH<6> 0.03 |                         0.06 | -3.03 GLYPH<6> 0.05                                            | -6.13 GLYPH<6> 0.05                               |
|    41 | HD 227607       |     2059885801468671744 |              301.44  |                36.52   |       9.88 | 0.52 GLYPH<6> 0.03 |                         0.05 | -4.36 GLYPH<6> 0.04                                            | -6.26 GLYPH<6> 0.05                               |

and simulated histograms. It is found that the parallax uncertainties in the calculation of GLYPH&lt;27&gt; sim would have to be increased by a factor of 3 to get zero intrinsic variance. This is unlikely to be produced by systematic effects and therefore, is suggestive of either a real extension of the stars or contamination by non-members. However, if we adopt the maximum global systematic uncertainty of 0.1 mas stated by the Gaia team (Luri et al. 2018), the result is consistent with zero intrinsic variance.

The mean and median distance for the sample are found to be, respectively, 1.99 kpc and 1.92 kpc, with a standard deviation of 0.31 kpc. The distance of 1.92 GLYPH&lt;6&gt; 0.31 kpc for the association Cyg OB3 is in agreement with the findings of Dambis et al. (2001), Mel'Nik et al. (2001), Mirabel &amp; Rodrigues (2003), Mel'Nik &amp; Efremov (1995), and Mel'nik

&amp; Dambis (2017). The distance to Cyg X-1 is found to be 2.38 GLYPH&lt;6&gt; 0.19 kpc (Gandhi et al. 2019), consistent with the distance to the association. The nearest stars are at a distance of about GLYPH&lt;24&gt; 1.5 kpc, as constrained by the minimum distance in the archive query.

## 4 PROPER MOTIONS

Gaia DR2 provides proper motion measurements for a large number of stars in the region covered by Cyg OB3. Proper motions are defined as GLYPH&lt;22&gt; = q GLYPH&lt;22&gt; 2 GLYPH&lt;11&gt; GLYPH&lt;3&gt; + GLYPH&lt;22&gt; 2 GLYPH&lt;14&gt; , where GLYPH&lt;22&gt; GLYPH&lt;11&gt; GLYPH&lt;3&gt; ( = GLYPH&lt;22&gt; GLYPH&lt;11&gt; cos GLYPH&lt;14&gt; ) and GLYPH&lt;22&gt; GLYPH&lt;14&gt; are proper motions along RA ( GLYPH&lt;11&gt; ) and Dec ( GLYPH&lt;14&gt; ), respectively. The proper motions are shown as a function of parallax in Fig 3. Since nearby stars usually exhibit

Table 2. Spectral types of the stars mentioned in the third column are obtained from SIMBAD . The distances, r exp , in the fourth column are the heliocentric distances to the stars calculated using exponential prior. v pec mentioned in the fifth column are the peculiar velocities of stars after removing the components of Solar motion and Galactic rotation. v rel are the relative velocities of individual stars with respect to the association. Radial velocities v r given in the last column are obtained from Kharchenko et al. (2007).

|   Nr. | Source          | Spectral type   | r exp (kpc)        | v pec (km s GLYPH<0> 1 )   | v rel (km s GLYPH<0> 1 )   | v r (km s GLYPH<0> 1 )   |
|-------|-----------------|-----------------|--------------------|----------------------------|----------------------------|--------------------------|
|     1 | HD 227415       | B3              | 1.76 GLYPH<6> 0.08 | 28.3 GLYPH<6> 7.8          | 21.4 GLYPH<6> 10.0         | -25.0 GLYPH<6> 10.0      |
|     2 | HD 227943       | B8              | 1.54 GLYPH<6> 0.09 | -                          | -                          | -                        |
|     3 | HD 228022       | B3:III:         | 1.70 GLYPH<6> 0.12 | -                          | -                          | -                        |
|     4 | HD 228041       | B0.5V:e         | 2.09 GLYPH<6> 0.17 | -                          | -                          | -                        |
|     5 | HD 191495       | B0IV-V(n)       | 1.70 GLYPH<6> 0.15 | 23.2 GLYPH<6> 4.9          | 17.0 GLYPH<6> 5.5          | 5.0 GLYPH<6> 3.7         |
|     6 | HD 191567       | B1V             | 1.59 GLYPH<6> 0.11 | 30.5 GLYPH<6> 5.0          | 24.6 GLYPH<6> 5.1          | -29.0 GLYPH<6> 1.8       |
|     7 | HD 228104       | B1:IV:pe        | 1.98 GLYPH<6> 0.12 | -                          | -                          | -                        |
|     8 | HD 191917       | B1III           | 2.26 GLYPH<6> 0.19 | 22.7 GLYPH<6> 5.3          | 13.6 GLYPH<6> 5.3          | -18.0 GLYPH<6> 3.7       |
|     9 | HD 228282       | B8              | 1.69 GLYPH<6> 0.17 | -                          | -                          | -                        |
|    10 | HD 190967       | O9.5V+B1Ib      | 2.03 GLYPH<6> 0.14 | 20.8 GLYPH<6> 4.8          | 8.8 GLYPH<6> 4.2           | -16.1 GLYPH<6> 1.8       |
|    11 | HD 227722       | B1III           | 1.70 GLYPH<6> 0.08 | 21.8 GLYPH<6> 4.9          | 13.0 GLYPH<6> 7.0          | -1.0 GLYPH<6> 7.4        |
|    12 | BD+35 3976      | OB-             | 1.74 GLYPH<6> 0.09 | -                          | -                          | -                        |
|    13 | HD 190919       | B0.7Ib          | 2.32 GLYPH<6> 0.15 | 26.2 GLYPH<6> 5.2          | 11.3 GLYPH<6> 4.4          | -15.0 GLYPH<6> 3.7       |
|    14 | HD 227634       | B0.2II          | 1.85 GLYPH<6> 0.11 | 17.5 GLYPH<6> 4.5          | 7.7 GLYPH<6> 3.6           | -10.0 GLYPH<6> 3.7       |
|    15 | LS II +35 32    | B1V             | 1.71 GLYPH<6> 0.11 | -                          | -                          | -                        |
|    16 | BD+35 3955      | B0.7Iab         | 1.88 GLYPH<6> 0.11 | 15.3 GLYPH<6> 4.3          | 9.1 GLYPH<6> 4.3           | -6.0 GLYPH<6> 1.8        |
|    17 | HD 227586       | B0.5IVp         | 2.63 GLYPH<6> 0.32 | 19.6 GLYPH<6> 5.3          | 12.0 GLYPH<6> 5.8          | -6.0 GLYPH<6> 7.4        |
|    18 | HD 227621       | B1.5IV          | 1.87 GLYPH<6> 0.09 | -                          | -                          | -                        |
|    19 | BD+35 3956      | B0.5Vne         | 1.85 GLYPH<6> 0.11 | 19.7 GLYPH<6> 6.5          | 13.0 GLYPH<6> 6.3          | -19.0 GLYPH<6> 7.4       |
|    20 | HD 190864       | O6.5III(f)      | 2.05 GLYPH<6> 0.13 | 17.6 GLYPH<6> 4.6          | 11.5 GLYPH<6> 4.9          | -2.8 GLYPH<6> 3.6        |
|    21 | HD 227696       | B1III           | 1.52 GLYPH<6> 0.09 | 17.0 GLYPH<6> 4.7          | 8.8 GLYPH<6> 4.3           | -15.0 GLYPH<6> 3.7       |
|    22 | HD 227680       | B3II            | 2.56 GLYPH<6> 0.21 | -                          | -                          | -                        |
|    23 | HD 227611       | B1:III/Ve       | 2.08 GLYPH<6> 0.13 | 25.9 GLYPH<6> 8.2          | 15.1 GLYPH<6> 8.3          | -24.0 GLYPH<6> 10.0      |
|    24 | HD 191612       | O8fpe           | 2.14 GLYPH<6> 0.17 | 26.3 GLYPH<6> 5.4          | 17.3 GLYPH<6> 5.2          | -27.6 GLYPH<6> 3.3       |
|    25 | HD 227883       | B               | 2.19 GLYPH<6> 0.14 | -                          | -                          | -                        |
|    26 | HD 227960       | O8.5            | 2.14 GLYPH<6> 0.15 | 12.3 GLYPH<6> 4.4          | 11.2 GLYPH<6> 4.3          | -8.0 GLYPH<6> 3.7        |
|    27 | HD 191611       | B0.5III         | 2.71 GLYPH<6> 0.27 | 13.4 GLYPH<6> 5.5          | 26.6 GLYPH<6> 5.8          | 4.0 GLYPH<6> 3.7         |
|    28 | HD 227757       | O9.5V           | 1.87 GLYPH<6> 0.09 | 28.3 GLYPH<6> 4.9          | 43.8 GLYPH<6> 4.7          | -16.0 GLYPH<6> 7.0       |
|    29 | HD 191139       | B0.5III         | 2.19 GLYPH<6> 0.16 | 22.2 GLYPH<6> 4.9          | 8.2 GLYPH<6> 3.9           | -13.0 GLYPH<6> 3.7       |
|    30 | BD+36 3882      | B1III:          | 1.85 GLYPH<6> 0.08 | -                          | -                          | -                        |
|    31 | HD 226868       | O9.7Iabpvar     | 2.38 GLYPH<6> 0.19 | 22.2 GLYPH<6> 4.3          | 11.4 GLYPH<6> 3.0          | -5.1 GLYPH<6> 0.5        |
|    32 | HD 227245       | O7              | 1.60 GLYPH<6> 0.07 | 21.1 GLYPH<6> 4.2          | 9.7 GLYPH<6> 3.2           | -13.0 GLYPH<6> 3.7       |
|    33 | HD 227018       | O6.5III         | 2.12 GLYPH<6> 0.15 | 52.2 GLYPH<6> 5.2          | 42.4 GLYPH<6> 5.6          | 20.0 GLYPH<6> 4.4        |
|    34 | LS II +35 21    | B               | 1.55 GLYPH<6> 0.06 | -                          | -                          | -                        |
|    35 | HD 227132       | B2III           | 2.48 GLYPH<6> 0.15 | -                          | -                          | -                        |
|    36 | HD 226951       | B0.5III         | 1.99 GLYPH<6> 0.13 | 23.4 GLYPH<6> 4.0          | 11.7 GLYPH<6> 2.4          | -9.0 GLYPH<6> 2.9        |
|    37 | HD 227070       | B2              | 2.54 GLYPH<6> 0.20 | -                          | -                          | -                        |
|    38 | TYC 2682-2524-1 | B5?             | 2.36 GLYPH<6> 0.16 | -                          | -                          | -                        |
|    39 | BD+35 3929      | B1III           | 1.85 GLYPH<6> 0.11 | -                          | -                          | -                        |
|    40 | HD 227460       | B0.5:V          | 1.78 GLYPH<6> 0.10 | 14.2 GLYPH<6> 4.4          | 8.3 GLYPH<6> 3.9           | -10.0 GLYPH<6> 3.7       |
|    41 | HD 227607       | B1:Ib:          | 1.92 GLYPH<6> 0.09 | 27.6 GLYPH<6> 5.2          | 19.8 GLYPH<6> 6.7          | 5.0 GLYPH<6> 7.4         |

higher proper motions, we expect to find a correlation between parallax and proper motion, as observed in the figure. However, a significantly higher proper motion for HD 227018 [ 33 ] and lower values for four of the stars viz. HD 227943 [ 2 ], HD 228041 [ 4 ], HD 228104 [ 7 ] and HD 227757 [ 28 ] marked in the figure are evident outliers. The number mentioned in parenthesis indicates the index of the star according to the serial number in Tables 1 and 2. Cyg X-1 does not exhibit any peculiarity in its proper motion, when compared with the other stars. We investigated the presence of kinematic substructures in the sample based on their proper motion using Shapiro-Wilk test. The test suggests that the distributions of GLYPH&lt;22&gt; GLYPH&lt;11&gt; GLYPH&lt;3&gt; and GLYPH&lt;22&gt; GLYPH&lt;14&gt; do not follow single normal distributions, and hence indicate the presence of kinetic subgroups. The distributions are consistent with single normal distributions when stars 2 , 4 , 7 , 28 , 33 , and 37 are removed from the test.

Fig 4 shows the proper motions of the stars in the l GLYPH&lt;0&gt; b plane. The starting positions of the arrows correspond to their locations extrapolated back in time by 0.1 Myr, with the arrowheads shown at their present positions. The length of each arrow is proportional to the proper motion of the respective star. The serial numbers of some of the stars from Tables 1 and 2 are shown in the figure. Not all stars are marked for the sake of clarity. Although the proper motions of most of the stars point in a similar direction and have similar magnitudes, some anomalies are observed. For example, the smaller proper motion of HD 227943 [ 2 ] is noticeable in the figure.

We also calculate the proper motion of the association

Figure 2. Histogram of heliocentric distance estimates, r exp , obtained from exponential prior.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-ace4bab5-adf3-408f-a25e-612a339b0129.jpg)

Figure 3. Parallax as a function of proper motion for the stars in our sample. Cyg X-1 is shown in red. The serial numbers for some of the stars from Table (1) are shown in the figure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-55fe95f5-e344-4da9-81c8-e8b45a3d2ab6.jpg)

by determining the median l GLYPH&lt;0&gt; b coordinates of the stars and extrapolating into the past using the median proper motion of the stars. The proper motion of the association is presented with a yellow arrow in the figure. Note that the yellow arrow is very similar to that of Cyg X-1, shown in red, both in terms of direction and magnitude. This result, obtained from Gaia DR2 data for 41 stars, is consistent with the findings of Mirabel &amp; Rodrigues (2003), who studied the association using Hipparcos data for 22 stars.

Galactic rotation and Solar motion significantly contribute to the observed proper motions of stars in the sky. Although the study of the proper motions of any sample provides important clues about the peculiarity in the kinematics of stars, it is important to exclude these two strong contributors to the proper motion in order to find the intrinsic motion of the stars. Therefore, next we calculate the

Figure 4. Stars are shown in the Galactic l GLYPH&lt;0&gt; b plane. The grey arrows represent their proper motions over the past 0.1 Myr. Serial numbers for some of the stars from Table (1) are marked in the figure. The red and green arrows represent GLYPH&lt;22&gt; for Cyg X-1 [ 31 ] and HD 227018 [ 33 ], respectively. The proper motion of the center of the association is shown via the yellow arrow.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-0a04d8ae-d67c-4b99-84fa-25eb7c82d201.jpg)

peculiar velocities of the stars in our sample, accounting for the full three-dimensional motions of the stars, including the line of sight, radial velocity.

## 5 PECULIAR VELOCITIES

The estimation of the peculiar velocity relative to Galactic rotation and Solar motion for any star requires a full astrometric solution, i.e. parallax, proper motion and radial velocity ( v r ), as input. Gaia has a spectrometer on board that can measure radial velocities, however none of the stars studied in this work have measured DR2 v r values. Therefore, these measurements were extracted from the catalogue published by Kharchenko et al. (2007). The catalogue lists 25 of the stars given in Table 1 for which peculiar velocities ( v pec ) and subsequent inferences are discussed herein. The v r for Cyg X-1 is taken from Gies et al. (2008). The formalism of Johnson &amp; Soderblom (1987) is used to obtain the heliocentric space velocities ( U V W ; ; ). Peculiar velocities were calculated by removing Solar motion and Galactic rotation according to the formalism and constants given in Reid et al. (2009, 2014) and is defined as the space velocity after removing Galactic rotation and Solar motion, v pec =

q U 2 s + V 2 s + W 2 s , where , U s , V s and W s are the velocity components towards the Galactic centre, along the direction of Galactic rotation, and towards the North Galactic pole.

The mean and standard deviation of 10,000 realizations of v pec have been quoted in Table 2 by drawing random values of $ GLYPH&lt;22&gt; , and v r , accounting for the covariance between $ and GLYPH&lt;22&gt; quoted in the Gaia DR2 astrometric solution. The calculation of v pec also accounts for the uncertainties of the circular rotation speed of the Galaxy at the location of the Sun ( GLYPH&lt;2&gt; 0 ), Solar motion towards the Galactic centre ( U GLYPH&lt;12&gt; ),

Figure 5. Acomparison between the velocities of stars relative to the association and their peculiar velocities is shown. The dotted line corresponds to equality between the two quantities. The serial numbers for some of the stars from Table (1) are shown in the figure. Cyg X-1 is shown in red.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-91a33729-1395-4656-a3b9-74a656339cb0.jpg)

Solar motion in the direction of Galactic rotation ( V GLYPH&lt;12&gt; ), Solar motion towards the north Galactic pole ( W GLYPH&lt;12&gt; ) and the distance of the Sun from the Galactic centre ( R GLYPH&lt;12&gt; ) from Reid et al. (2014).

A majority of the stars exhibit v pec of GLYPH&lt;24&gt; 22 km s GLYPH&lt;0&gt; 1 , with a median value of 22 km s GLYPH&lt;0&gt; 1 and a range of 10-30 km s GLYPH&lt;0&gt; 1 , except HD 227018 [ 33 ] that has a peculiar velocity of GLYPH&lt;24&gt; 52 km s GLYPH&lt;0&gt; 1 . The velocity dispersion, measured as the v pec standard deviation is equal to 7.7 km s GLYPH&lt;0&gt; 1 . This velocity dispersion is small considering the typical v pec uncertainties of individual stars given in Table 2 of about 6 km s GLYPH&lt;0&gt; 1 . This is consistent with Mel'nik &amp; Dambis (2017) and Mel'Nik &amp; Efremov (1995).

The peculiar velocity of Cyg X-1 is found to be 22.2 GLYPH&lt;6&gt; 4.3 km s GLYPH&lt;0&gt; 1 , consistent with the mean of the association. This result is consistent with the findings of Mirabel &amp; Rodrigues (2003) that the two systems have low peculiar velocities (see also Gandhi et al. 2019). The match also indicates that Cyg X-1 is a member of the common group formed by the stars in the region. This is consistent with the inferences drawn by Mirabel (2017a) and Mirabel &amp; Rodrigues (2003) regarding the membership of Cyg X-1 in Cyg OB3.

We also calculate the relative three-dimensional velocities v rel of individual stars with respect to the association defined as v rel = p ' U s GLYPH&lt;0&gt; U m ' 2 + ' V s GLYPH&lt;0&gt; V m ' 2 + ' W s GLYPH&lt;0&gt; W m ' 2 , where U m , V m , W m are the median velocity components of the association. A majority of the stars exhibit v rel of about 15 km s GLYPH&lt;0&gt; 1 . The mean and standard deviation of 10,000 v rel realizations for each star are shown in Table 2. Astrometric covariances are taken into account. The mean v rel is 15.9 GLYPH&lt;6&gt; 9.4 km s GLYPH&lt;0&gt; 1 , where the quoted uncertainty is the standard deviation of the relative velocities of 25 stars.

The relative velocity of Cyg X-1 is estimated to be 11.4 GLYPH&lt;6&gt; 3.0 km s GLYPH&lt;0&gt; 1 , consistent with an estimate of 9 GLYPH&lt;6&gt; 2 km s GLYPH&lt;0&gt; 1 by Mirabel &amp; Rodrigues (2003). We compare the peculiar velocities of stars with their relative velocities and the results are shown in Fig 5. A positive correlation is observed between the two quantities. Most of the stars show v rel &lt; v pec , however a few exceptions are apparent viz. HD 191611 [ 27 ], and HD 227757 [ 28 ].

Figure 6. The figure shows components of space velocities U s , V s and W s . Cyg X-1 is shown in red. Median values of the components and their standard deviation are shown in blue. Serial numbers marked for some of the stars are from Table 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-2bb40e1b-950c-40ac-b915-4ab5c5a03d35.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-1c517d80-1753-43e1-b2ee-a6e493bd9002.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-fd5a8bd7-6659-40e7-b420-40b7c1576ebc.jpg)

## 5.1 Peculiar velocities ( U s , V s and W s )

We study the components of v pec viz. U s , V s and W s for our 25 stars, shown in Fig 6. The median values of the components that can be considered as the velocity components of the association are found to be ( U m , V m , W m ) = (16.6 GLYPH&lt;6&gt; 10.8, GLYPH&lt;0&gt; 4.6 GLYPH&lt;6&gt; 10.8, 2.6 GLYPH&lt;6&gt; 5.5) km s GLYPH&lt;0&gt; 1 , where the quoted uncertainties are the standard deviation for the sample. v pec is dominated

by the U s component along the direction to the Galactic centre, which has a median value of 16.6 GLYPH&lt;6&gt; 10.8 km s GLYPH&lt;0&gt; 1 for the sample. The first panel in Fig. 6 shows that the stars move as a coherent group in this direction, excluding three stars marked in the figure with their serial numbers. The standard deviation in U s reduces to 3.2 km s GLYPH&lt;0&gt; 1 , when the three stars are not considered. A majority of the stars are observed moving towards the North Galactic pole with a median velocity of 2.6 GLYPH&lt;6&gt; 5.5 km s GLYPH&lt;0&gt; 1 for the sample. The two stars HD 191611 [ 27 ] and HD 227757 [ 28 ], exhibiting higher relative velocities in Fig (5), have a velocity component U s that deviates from the rest of the stars. The velocity components for the runaway star HD 227018 [ 33 ] are the highest in all three directions. The velocity components for Cyg X-1 are found to be ( U s , V s , W s ) = (19.1 GLYPH&lt;6&gt; 4.9, GLYPH&lt;0&gt; 1.4 GLYPH&lt;6&gt; 4.9, 9.8 GLYPH&lt;6&gt; 1.2) km s GLYPH&lt;0&gt; 1 and it appears to be moving with the coherent structure formed by other stars in the association.

## 6 EVOLUTION OF THE ASSOCIATION

A publicly available tool galpy (Bovy 2015) allows the calculation of stellar orbits under a variety of Galactic potentials, and extrapolation of orbits in time. We make use of this facility to study the distance of individual stars with respect to the center of the association ( r rel ) at present, and also to study the projected RA and Dec positions after 10 Myr. The tool requires a six-parameter input as initial conditions to integrate the orbits viz. GLYPH&lt;11&gt; , GLYPH&lt;14&gt; , r , GLYPH&lt;22&gt; GLYPH&lt;11&gt; cos GLYPH&lt;14&gt; , GLYPH&lt;22&gt; GLYPH&lt;14&gt; , and v r . The calculation accounts for the Galactic potential, but does not consider the gravitational potential of the association or surrounding interstellar medium as these are expected to be unimportant (see below). The orbits of individual stars are integrated 1000 times by drawing random values of $ GLYPH&lt;22&gt; , and v r . The covariance between $ and GLYPH&lt;22&gt; has been taken into account. Most of the stars are found to be within GLYPH&lt;24&gt; 500 pc from the centre of the association at present, though the uncertainties on r rel are large. The results of our error analysis reveal that the parallax is the most significant contributor to the uncertainties in r rel estimates.

We study the current and projected RA and Dec positions over the next 10 Myr on the sky plane on parsec scales using galpy assuming all stars lie at the median association distance. The current positions for 41 stars and the projected positions for 25 of the stars are shown in Fig 7 in black and grey, respectively. r GLYPH&lt;11&gt; and r GLYPH&lt;14&gt; are the distances from the centre of the association along RA and Dec. The light grey arrows connect the corresponding stars at T=0 and T=10 Myr. The figure suggests an expansion of the association in the plane of the sky with time when only the Galactic potential is considered, although the association is found to be moving towards the Galactic center (Fig 6). We estimated the gravitational force on a given star due to the total mass of the association and compared it to the force exerted on the star by the Galactic potential. It was found that the binding force due to the association on a given member is highly insignificant when compared with the force due to the Galactic potential, rendering the association as an unbound entity that is slowly expanding with time.

Figure 7. Extrapolated expansion of the association in the sky plane when only the Galactic potential is considered. r GLYPH&lt;11&gt; and r GLYPH&lt;14&gt; are the distances from the centre of the association along the RA and Dec directions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00810v2-figure-7ff3de76-abbd-4598-ab3f-4c2d27aaf974.jpg)

## 7 DISCUSSION

The results obtained from Gaia DR2 for the parallax, proper motion and peculiar velocities of Cyg X-1 and Cyg OB3 are consistent with those obtained using the Hipparcos mission. This further supports the identification of Cyg OB3 as the parent association of Cyg X-1. We find that, although the parallaxes obtained with Gaia are the most accurate and precise measurements available so far, the measurement uncertainties are large enough to introduce scatter in the membership determination of individual stars, as well as a measurement of the physical scale of the association along the line of sight.

We studied the association Cyg OB3 using the full astrometric solution provided by Gaia DR2 and radial velocities obtained from the literature. The distribution of distances (Fig 2) shows that the majority of the stars lie at a distance of 1.92 GLYPH&lt;6&gt; 0.31 kpc.

The star HD 227018 [ 33 ] shows the highest proper motion and HD 227943 [ 2 ], HD 228041 [ 4 ], HD 228104 [ 7 ] and HD 227757 [ 28 ] are the stars with the smallest proper motions in the sample (Fig 3). Peculiar velocities could be calculated only for HD 227018 [ 33 ] and HD 227757 [ 28 ]. The study of v pec values further corroborates the cohesive nature of Cyg OB3 as the majority of stars have peculiar velocities of about 22 km s GLYPH&lt;0&gt; 1 , with a small velocity dispersion. A comparison of v pec and v rel of individual stars in the association has provided important evidence about the kinematics of stars in the association, and helped to identify those sources with peculiar kinematic behavior. Some of the stars identified in this way are discussed in detail below.

## 7.1 Cyg X-1 / HD 226868 [ 33 ]

The black hole binary Cyg X-1 is located at the boundary of Cyg OB3. It is important to identify its parent association in order to understand the formation of black holes and natal kicks (see Blaauw 1961; Lyne &amp; Lorimer 1994; van Paradijs &amp; White 1995; White &amp; van Paradijs 1996; Hansen &amp; Phinney 1997; Jonker &amp; Nelemans 2004; Repetto et al. 2012)

imparted during supernovae explosions. Since a majority of known black hole binaries host low mass companions that can be long lived, they could have travelled large distances from their parent associations or site of birth. Therefore it is difficult to identify the parent association of a system with an evolved star and high peculiar velocity. Cyg X-1 is an important black hole in this regard as it is a nearby and young binary with a massive companion and it has been found to exhibit low peculiar velocities (Mirabel &amp; Rodrigues 2003; Mirabel 2017a,b).

It has been suggested by Mirabel (2017a) that Cyg X-1 is at the same distance as the stellar association Cyg OB3 using data from Hipparcos and VLBI radio observations. There are three distance estimates for Cyg X-1 from geometric parallaxes available in the literature. Hipparcos' new reduction places the source at a distance of 0.60 GLYPH&lt;6&gt; 0.32 kpc, while the distance measured with VLBI radio observations is 1.86 + 0 12 : GLYPH&lt;0&gt; 0 11 : kpc (Reid et al. 2011). The latest distance estimate using geometric parallax measurements from Gaia DR2 is 2.37 GLYPH&lt;6&gt; 0.18 kpc (Gandhi et al. 2019). Massey et al. (1995) measure the distance to the black hole binary to be 2.14 GLYPH&lt;6&gt; 0.07 kpc using photometric and geometric methods. Though more work is needed to understand systematic differences, all of the above mentioned distance estimates, excluding the distance from Hipparcos' new reduction, are in agreement at the 2 GLYPH&lt;27&gt; level of confidence, and hence are consistent with the inferred distance to the centre of the association.

The relative velocity of Cyg X-1 is very small with respect to the association, as has been discussed by Mirabel &amp; Rodrigues (2003). The results for r , GLYPH&lt;22&gt; , v pec and v rel , obtained with the astrometric data from Gaia DR2, confirm the findings of Mirabel &amp; Rodrigues (2003) and support the conjecture that Cyg X-1 is a member of Cyg OB3. The implication of this is that it formed in situ without a strong natal kick. Cyg X-1 is one of few dynamically measured black holes in high mass X-ray binaries, and its low kick is consistent with weak trends of mass dependent natal kicks observed in black hole binaries (e.g. Gandhi et al. 2019).

## 7.2 HD 227943 [ ], HD 228041 [ ] and HD 228104 2 4 [ 7 ]

These three stars show significantly smaller values of GLYPH&lt;22&gt; as compared to the rest of the sample. All three are located towards the low end of the range of Galactic latitudes within Cyg OB3. HD 227943 [ ] is located at the boundary of the 3 region (Fig 4) selected in this work. HD 227943 [ ] 3 and HD 228041 [ 5 ] are not in the list of stars of Cyg OB3 published by Garmany &amp; Stencel (1992), however HD 228104 [ 8 ] is included. The peculiar velocities for all three stars could not be studied as their radial velocities are not available. Assuming v r is equal to the median value of the association, the v pec and v rel of the three stars are GLYPH&lt;24&gt; 22 km s GLYPH&lt;0&gt; 1 and GLYPH&lt;24&gt; 38 km s GLYPH&lt;0&gt; 1 , respectively.

## 7.3 HD 227757 [ 28 ]

HD 227757 [ 28 ] shows a smaller magnitude of proper motion relative to the mean of the association, as shown in Fig 3. In addition, the direction of proper motion differs in Galactic latitude from the rest of the stars in the sample as shown by the blue arrow in Fig 4. The estimated v pec of the object is small ( GLYPH&lt;24&gt; 28 km s GLYPH&lt;0&gt; 1 ), and therefore the object cannot be classified as a runaway star ( v pec &gt; 30 km s GLYPH&lt;0&gt; 1 ; Blaauw 1961). Xu et al. (2018) measure its velocity with respect to the local standard of rest to be 22.9 GLYPH&lt;6&gt; 20.0 km s GLYPH&lt;0&gt; 1 . Mahy et al. (2013) found no significant variation in the radial velocity measurements for this star from spectroscopic studies and suggested the star as 'presumably' single. A different direction of proper motion from the majority of the stars in the sample may suggest a possible interaction or encounter in the past. This is further supported by the higher v rel = 43.8 GLYPH&lt;6&gt; 4.7 km s GLYPH&lt;0&gt; 1 (see Fig 5).

## 7.4 HD 227018 [ 33 ]

HD 227018 has manifested its distinct kinematic behavior in all diagnostics studied in this work. The highest proper motion of this star amongst our sample is evident from Fig 4. It shows v pec of GLYPH&lt;24&gt; 50 km s GLYPH&lt;0&gt; 1 (Fig 5) and it is a known runaway star moving with a velocity of &gt; 30 km s GLYPH&lt;0&gt; 1 (Blaauw 1961; van Buren et al. 1995). Its relative velocity v rel is also higher than for the other stars ( GLYPH&lt;24&gt; 40 km s GLYPH&lt;0&gt; 1 ). Tetzlaff et al. (2011) measure a peculiar spatial velocity of 37.4 + 4 8 : GLYPH&lt;0&gt; 7 2 : km s GLYPH&lt;0&gt; 1 and classify the source as a runaway star using Hipparcos data. The star was not included in Cyg OB3 (Cyg 3A and Cyg 3B) by Mel'Nik &amp; Efremov (1995), but was considered as a member of NGC 6871, the core of the Cyg OB3 association, by Massey et al. (1995). As discussed in Gvaramadze et al. (2009) and Gvaramadze &amp; Gualandris (2011), the high velocity of runaway stars can be a consequence of either the disruption of a massive binary following a supernova explosion (Blaauw 1961; Stone 1991; Leonard &amp; Dewey 1993; Iben &amp; Tutukov 1996) or dynamical three or four-body encounters in dense stellar systems (Poveda et al. 1967; Aarseth 1974; Gies &amp; Bolton 1986; Leonard &amp; Duncan 1990; Leigh &amp; Sills 2011; Leigh &amp; Geller 2012; Leigh et al. 2016; Leigh &amp; Wegsman 2018; Ryu et al. 2017a,b). Our analysis reveals that the star HD 227611 [ 23 ] came to a distance of closest approach of GLYPH&lt;24&gt; 20 pc about 1 Myr ago, however this result in not conclusive of a possible interaction between the two stars resulting in the high velocity of the runaway star.

## 8 SUMMARY

Astrometric data from Gaia DR2 for 41 stars located in the region covered by Cyg OB3 are studied in this work. The stars in the sample are found to be concentrated at a heliocentric distance of about 2 kpc. Most of the stars are observed to be moving with peculiar velocity, v pec of about 22 km s GLYPH&lt;0&gt; 1 . One kinematic subgroup has been identified consisting of majority of the sample stars. A small dispersion in v pec conforms to the results suggesting that the association forms a coherent structure in velocity space (e.g. Brown et al. 1999; Tian et al. 1996; Mathieu 1986). OB associations are unbound systems formed from single molecular clouds, wherein the stars begin to drift apart after the escape of dust and gas (Humphreys 1978; Mel'nik &amp; Dambis 2017). We have calculated the relative velocities of individual stars with respect to the association, which can be considered as a measure of this drift. The small relative velocities of &lt; 20

km s GLYPH&lt;0&gt; 1 for a majority of the stars suggest that Cyg OB3 is a slowly expanding association. Since the relative velocities of the stars are very close to or only marginally different from their respective peculiar velocities, we conclude that expansion is a significant contributor to the peculiar velocities of individual stars.

The kinematic characteristics of Cyg X-1 are found to be consistent with those of the majority of stars in the sample. The membership of Cyg X-1 in the association is not ruled out, despite the location of the binary at the boundary of the association. This result is in agreement with Mirabel (2017a), Mirabel (2017b), and Mirabel &amp; Rodrigues (2003) that Cyg X-1 is most likely a member of the association and was formed in situ.

Our analysis reveals interesting kinematic behavior for some of the stars; e.g. HD 227757 [ 28 ], which exhibits a proper motion in a direction different to the majority of the stars in the association. HD 191611 [ 27 ] is found to have v rel &gt; v pec . The peculiar kinematic properties of these stars might be remnants of past encounters of the stars in a dense environment. The status of HD 227018 as a runaway star is further supported by data from Gaia DR2. The rest of the stars in the sample display v pec glyph[lessorsimilar] 30 km s GLYPH&lt;0&gt; 1 and are not classified as runaway stars.

Finally, projected positions of the sample stars on the Galactic xy-plane reveals an extended structure along the line of sight. This is likely an artifact of the uncertainties in the parallax measurements. At the present level of uncertainties, we restricted ourselves to the study of the kinematic properties of the stars and refrained from addressing spatial structure and from identifying members and non-members in the association. Future Gaia data releases with improved uncertainties will make such studies more feasible.

## ACKNOWLEDGMENT

The authors thank the anonymous reviewer for useful comments that helped to improve the paper. This work has made use of data from the European Space Agency (ESA) mission Gaia ( https://www.cosmos.esa.int/gaia ), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://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. AR acknowledges a Commonwealth Rutherford Fellowship. JP is in part supported by funding from a University of Southampton Central VC Scholarship. PG thanks STFC for support (ST/R000506/1). NWCL gratefully acknowledges a Fondecyt Iniciac´ ıon grant (#11180005). DB thanks Magdalen College for his fellowship and the Rudolf Peierls Centre for Theoretical Physics for providing office space and travel funds. This work is supported by a UGC-UKIERI Phase 3 Thematic Partnership. This research has made use of observatory archival image from DSS2.

## REFERENCES

Ambartsumian V. A., 1947, in Stellar Evolution and Astrophysics

| (Armenian Academy of Science [German translation, Ab- handl. Sowjetischen. Astron. 1, 33, (1951)]). Ambartsumian V. A., 1955, The Observatory, 75, 72 Anders F., et al., 2019, arXiv e-prints, p. arXiv:1904.11302 Arenou F., et al., 2018, A&A, 616, A17                                                                                                                                           |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Gaia Collaboration et al., 2018, A&A, 616, A1 Gandhi P., Rao A., Johnson M. A. C., Paice J. A., Maccarone T. J., 2019, MNRAS, 485, 2642 Garmany C. D., Stencel R. E., 1992, A&AS, 94, 211                                                                                                                                                                                                           |
| RAS, 396, 570 Hansen B. M. S., Phinney E. S., 1997, MNRAS, 291, 569 Humphreys R. M., 1978, ApJS, 38, 309 Humphreys R. M., McElroy D. B., 1984, ApJ, 284, 565 Iben Jr. I., Tutukov A. V., 1996, ApJ, 456, 738                                                                                                                                                                                        |
| Johnson D. R. H., Soderblom D. R., 1987, AJ, 93, 864                                                                                                                                                                                                                                                                                                                                                |
| Jonker P. G., Nelemans G., 2004, MNRAS, 354, 355                                                                                                                                                                                                                                                                                                                                                    |
| Berlanas S. R., Wright N. J., Herrero A., Drew J. E., Lennon D. J., 2019, MNRAS, 484, 1838 Blaauw A., 1961, Bull. Astron. Inst. Netherlands, 15, 265 Blaauw A., 1964, ARA&A, 2, 213 Blaha C., Humphreys R. M., 1989, AJ, 98, 1598 Bolton C. T., 1972, Nature, 235, 271                                                                                                                              |
| Bovy J., 2015, ApJS, 216, 29                                                                                                                                                                                                                                                                                                                                                                        |
| Brown A. G. A., Blaauw A., Hoogerwerf R., de Bruijne J. H. de Zeeuw P. T., 1999, in Lada C. J., Kylafis N. D.,                                                                                                                                                                                                                                                                                      |
| NATO Advanced Science Institutes (ASI) Series C Vol. NATO Advanced Science Institutes (ASI) Series C. p. ( arXiv:astro-ph/9902234 ) Cantat-Gaudin T., Mapelli M., Balaguer-N´˜ez un L., Jordi Sacco G., Vallenari A., 2019, A&A, 621, A115 Clark P. C., Bonnell I. A., Zinnecker H., Bate M. R., 2005, RAS, 359, 809 Dambis A. K., Mel'Nik A. M., Rastorguev A. S., 2001, Astronomy Letters, 27, 58 |
| J., eds, 540, 411 C., MN-                                                                                                                                                                                                                                                                                                                                                                           |
| Gies D. R., Bolton C. T., 1986, ApJS, 61, 419 Gies D. R., et al., 2008, ApJ, 678, 1237 Gvaramadze V. V., Gualandris A., 2011, MNRAS, 410, 304 Gvaramadze V. V., Gualandris A., Portegies Zwart S., 2009,                                                                                                                                                                                            |
| MN- N. V., Scholz R.-D., Piskunov A. E., Roeser S., Schilbach E., 2007, VizieR Online Data Catalog, 3254                                                                                                                                                                                                                                                                                            |
| Kharchenko                                                                                                                                                                                                                                                                                                                                                                                          |
| Leigh N., Geller A. M., 2012, MNRAS, 425, 2369                                                                                                                                                                                                                                                                                                                                                      |
| Leigh N., Sills A., 2011, MNRAS, 410, 2370                                                                                                                                                                                                                                                                                                                                                          |
| Leigh N. W. C., Wegsman S., 2018, MNRAS, 476, 336                                                                                                                                                                                                                                                                                                                                                   |
| Leigh N. W. C., Stone N. C., Geller A. M., Shara M. M.,                                                                                                                                                                                                                                                                                                                                             |
| H., Solano-Oropeza D., Thomas Y., 2016, MNRAS, 463,                                                                                                                                                                                                                                                                                                                                                 |
| Muddu Leonard P. J. T., Dewey R. J., 1993, in Sasselov D. D., ed.,                                                                                                                                                                                                                                                                                                                                  |
| tronomical Society of the Pacific Conference Series Vol.                                                                                                                                                                                                                                                                                                                                            |
| Luminous High-Latitude Stars. p. 239                                                                                                                                                                                                                                                                                                                                                                |
| Leonard P. J. T., Duncan M. J., 1990, AJ, 99, 608                                                                                                                                                                                                                                                                                                                                                   |
| B., 1941, Stockholms Observatoriums Annaler, 13,                                                                                                                                                                                                                                                                                                                                                    |
| Lindblad Lindblad B., 1942, Stockholms Observatoriums Annaler, 14,                                                                                                                                                                                                                                                                                                                                  |
| Lindegren L., et al., 2018, A&A, 616, A2 Luri X., et al., 2018, A&A, 616, A9 Lyne A. G., Lorimer D. R., 1994, Nature, 369, 127 Mahy L., Rauw G., De Becker M., Eenens P., Flores C. A.,                                                                                                                                                                                                             |
| A&A, 550, A27                                                                                                                                                                                                                                                                                                                                                                                       |
| Massey P., Johnson K. E., Degioia-Eastwood K., 1995, ApJ,                                                                                                                                                                                                                                                                                                                                           |
| 151                                                                                                                                                                                                                                                                                                                                                                                                 |
| Mathieu R. D., 1986, Highlights of Astronomy, 7, 481                                                                                                                                                                                                                                                                                                                                                |
| A. M., Efremov Y. N., 1995, Astronomy Letters, 21,                                                                                                                                                                                                                                                                                                                                                  |
| Mel'Nik Mel'Nik A. M., Dambis A. K., Rastorguev A. S., 2001,                                                                                                                                                                                                                                                                                                                                        |
| 3311                                                                                                                                                                                                                                                                                                                                                                                                |
| As-                                                                                                                                                                                                                                                                                                                                                                                                 |
| 45,                                                                                                                                                                                                                                                                                                                                                                                                 |
| 8.1                                                                                                                                                                                                                                                                                                                                                                                                 |
| 2013, 454, 10 Astronomy                                                                                                                                                                                                                                                                                                                                                                             |
| 3.1                                                                                                                                                                                                                                                                                                                                                                                                 |
| Letters, 27, 521                                                                                                                                                                                                                                                                                                                                                                                    |
| Mel'nik A. M., Dambis A. K., 2017, MNRAS, 472, 3887                                                                                                                                                                                                                                                                                                                                                 |

MNRAS 000 , 1-11 (2020)

| Mirabel I. F., 2017b, in Gomboc A., ed., IAU Symposium Vol. 324, New Frontiers in Black Hole Astrophysics. pp 303-306 ( arXiv:1611.09266 ), doi:10.1017/S1743921316012904 Mirabel I. F., Rodrigues I., 2003, Science, 300, 1119                            |
|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Orosz J. A., McClintock J. E., Aufdenberg J. P., Remillard R. A., Reid M. J., Narayan R., Gou L., 2011, ApJ, 742, 84                                                                                                                                       |
| Poveda A., Ruiz J., Allen C., 1967, Boletin de los Observatorios Tonantzintla y Tacubaya, 4, 86                                                                                                                                                            |
| Reid M. J., et al., 2014, ApJ, 783, 130                                                                                                                                                                                                                    |
| S., Davies M. B., Sigurdsson S., 2012, MNRAS, 425,                                                                                                                                                                                                         |
| Repetto                                                                                                                                                                                                                                                    |
| 2799                                                                                                                                                                                                                                                       |
| Ruprecht J., Bal´zs a B., White R. E., 1981, in Catalogue of Star Clusters and Associations. Supplement 1., by Ruprecht, J.; Bal´zs, a B.; White, R. E.. Budapest (Hungary): Akad´miai e Ki- ad´, o Publishing House of the Hungarian Academy of Sciences, |
| 732 p.. Ryu T., Leigh N. W. C., Perna R., 2017a, MNRAS, 467, 4447 Ryu T., Leigh N. W. C., Perna R., 2017b, MNRAS, 470, 2 Stone R. C., 1991, AJ, 102, 333                                                                                                   |
| Tetzlaff N., Neuh¨user a R., Hohle M. M., 2011, MNRAS, 410, 190                                                                                                                                                                                            |
| K. P., van Leeuwen F., Zhao J. L., Su C. G., 1996, 118, 503                                                                                                                                                                                                |
| A&AS,                                                                                                                                                                                                                                                      |
| Walborn N. R., 1973, ApJ, 179, L123                                                                                                                                                                                                                        |
| Webster B. L., Murdin P., 1972, Nature, 235, 37 White N. E., van Paradijs J., 1996, ApJ, 473, L25 Wright N. J., Mamajek E. E., 2018, MNRAS, 476, 381                                                                                                       |
| Wright N. J., et al., 2019, MNRAS, 486, 2477                                                                                                                                                                                                               |
| Xu Y., et al., 2018, A&A, 616, L15 van Buren D., Noriega-Crespo A., Dgani R., 1995, AJ, 110, van Paradijs J., White N., 1995, ApJ, 447, L33                                                                                                                |
| 2914                                                                                                                                                                                                                                                       |
| Tian                                                                                                                                                                                                                                                       |
Anjali Rao Poshak Gandhi Christian Knigge John A. Paice Nathan W. C. Leigh Douglas Boubert
Kinematic study of the association Cyg OB3 with Gaia DR2
We study the stellar kinematic properties and spatial distribution of the
association Cyg OB3 using precise astrometric data from Gaia DR2. All known O-
and B-type stars in Cyg OB3 region with positions, parallaxes and proper
motions available are included, comprising a total of 41 stars. The majority of
stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3
kpc. The mean peculiar velocity of the sample after removing Galactic rotation
and solar motion is ~22 km/s, dominated by the velocity component towards the
Galactic center. The relative position and velocity of the black hole X-ray
binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most
likely its parent association. The peculiar kinematic properties of some of the
stars are revealed and are suggestive of past stellar encounters. The sample
includes a previously known runaway star HD 227018, and its high peculiar
velocity of ~50 km/s is confirmed with Gaia. We estimated the velocities of
stars relative to the association and the star HD 225577 exhibits peculiar
velocity smaller than its velocity relative to the association. The star has
lower value of proper motion than the rest of the sample. The results suggest a
slowly expanding nature of the association, which is supported by the small
relative speeds <20 km/s with respect to the association for a majority of the
sample stars.