Magnetic field geometry and chemical abundance distribution of the He-strong star CPD 573509
The magnetic field of CPD 3509 was recently detected in the framework of the BOB (B fields in OB stars) collaboration. We acquired low-resolution spectropolarimetric observations of CPD 3509 with FORS 2 and high-resolution UVES observations randomly distributed over a few months to search for periodicity, to study the magnetic field geometry, and to determine the surface distribution of silicon and helium. We also obtained supplementary photometric observations at a timeline similar to the spectroscopic and spectropolarimetric observations. A period of 6.36 d was detected in the measurements of the mean longitudinal magnetic field. A sinusoidal fit to our measurements allowed us to constrain the magnetic field geometry and estimate the dipole strength in the range of 3.9–4.5 kG. Our application of the Doppler imaging technique revealed the presence of He i spots located around the magnetic poles, with a strong concentration at the positive pole and a weaker one around the negative pole. In contrast, high concentration Si iii spots are located close to the magnetic equator. Further, our analysis of the spectral variability of CPD 3509 on short time scales indicates distinct changes in shape and position of line profiles possibly caused by the presence of Cep-like pulsations. A small periodic variability in line with the changes of the magnetic field strength is clearly seen in the photometric data.
keywords:Stars: abundances – Stars: evolution – Stars: magnetic field – Stars: massive – Stars: oscillations – Stars: individual: CPD 3509
Recently, Przybilla et al. (2016) presented a firm detection of a mean longitudinal magnetic field of kG order in the early-B type star CPD 3509, previously studied in a spectroscopic survey of massive stars in NGC 3293 by Evans et al. (2005). In their work, the authors focussed on the investigation of abundances, a model atmosphere, and the evolutionary state of the target in detail. The quantitative spectroscopic analysis of this star with the observed rather low -value of 35 km s yielded an effective temperature and a logarithmic surface gravity of K and , respectively, and a surface helium fraction of by number (see also Maeder et al. 2014). The surface abundances of C, N, O, Ne, S, and Ar were found to be compatible with the cosmic abundance standard (Nieva & Przybilla, 2012), whereas Mg, Al, Si, and Fe were depleted by about a factor of 2. It was suggested that such an abundance pattern can be understood as the consequence of a fractionated stellar wind. Importantly, CPD 3509 with an elapsed main-sequence life time of about 50% has evolved significantly away from the zero-age main sequence and appears to be one of the most evolved He-strong stars known with an independent age constraint due to its cluster membership. Since the evolution of the magnetic field geometry across the main sequence in massive B-type stars is not well sampled in comparison to studies of magnetic fields of late-type Bp and intermediate-mass Ap stars, a detailed study of the magnetic field configuration and the surface chemical inhomogeneities in a significantly evolved magnetic He-strong star is of particular interest. Knowledge of the evolution of the magnetic field geometry, especially of the distribution of the obliquity angle (the orientation of the magnetic axis with respect to the rotation axis) is essential to understand the physical processes taking place in these stars and the origin of their magnetic fields. Further, although it is generally assumed that CPD573509 is a member of the open cluster NGC 3293, no careful study of its membership involving new data from Gaia was carried out yet. We review the membership status involving the available Gaia data in Appendix A.
In order to characterise the properties of CPD 3509 in detail, we obtained time-series spectroscopy and spectropolarimetry with the FOcal Reducer low dispersion Spectrograph (FORS 2; Appenzeller et al. 1998) and the UV-Visual Echelle Spectrograph (UVES; Dekker et al. 2000) to constrain the magnetic field geometry and reconstruct the distribution of silicon and helium on the stellar surface using Doppler Imaging (DI). Further, we carried out photometric observations using the 40 cm Bochum Monitoring Telescope (BMT; Ramolla et al. 2013) of the Cerro Armazones Observatory. In the first part of the paper we report on the period determination using magnetic and photometric data. In the second part we present the results of the application of the Doppler Imaging technique followed by the discussion of the surface element distribution with respect to the magnetic field configuration.
2 Observations and magnetic field measurements
Fifteen FORS 2 spectropolarimetric observations of CPD 573509 were obtained in the framework of the ESO programme 094.D-0355 from 2014 November 17 to 2015 February 18, and further five within the framework of the ESO Large Programme 191.D-0255 during 2014 February and June, and 2015 March. The FORS 2 multi-mode instrument is equipped with polarisation analysing optics comprising super-achromatic half-wave and quarter-wave phase retarder plates, and a Wollaston prism with a beam divergence of 22 in standard resolution mode. We used the GRISM 600B and the narrowest available slit width of 04 to obtain a spectral resolution of . The observed spectral range from 3250 to 6215 Å includes all Balmer lines, apart from H, and numerous helium lines. For the observations, we used a non-standard readout mode with low gain (200kHz,11,low), which provides a broader dynamic range, hence allowing us to reach a higher signal-to-noise ratio (S/N) in the individual spectra. The exposure time for each subexposure accounted for 5 min. Each observation consisted of eight subexposures over approximately one hour including overheads.
Our first description of the assessment of longitudinal magnetic field measurements using FORS 1/2 spectropolarimetric observations was presented in several previous works (e.g. Hubrig et al. 2004a, b, and references therein). To minimize the cross-talk effect and to cancel errors from different transmission properties of the two polarised beams, a sequence of subexposures at the retarder position angles 4545, 4545, 4545, etc. is usually executed during the observations. Moreover, the reversal of the quarter wave plate compensates for fixed errors in the relative wavelength calibrations of the two polarised spectra. According to the FORS User Manual, the spectrum is calculated using:
where and indicate the position angle of the retarder waveplate and and are the ordinary and extraordinary beams, respectively. Rectification of the spectra was performed in the way described by Hubrig, Schöller & Kholtygin (2014a). Null profiles, , are calculated as pairwise differences from all available profiles. From these, 3-outliers are identified and used to clip the profiles. This removes spurious signals, which mostly come from cosmic rays, and also reduces the noise. A full description of the updated data reduction and analysis will be presented in a separate paper (Schöller et al., in preparation; see also Hubrig, Schöller & Kholtygin 2014a). The mean longitudinal magnetic field, , is measured on the rectified and clipped spectra based on the relation following the method suggested by Angel & Landstreet (1970):
where is the Stokes parameter that measures the circular polarization, is the intensity in the unpolarized spectrum, is the effective Landé factor, is the electron charge, is the wavelength, is the electron mass, is the speed of light, is the wavelength derivative of Stokes , and is the mean longitudinal (line-of-sight) magnetic field.
The longitudinal magnetic field was measured in two ways: using the entire spectrum including all available lines, or using exclusively hydrogen lines. Furthermore, we have carried out Monte Carlo bootstrapping tests. These are most often applied with the purpose of deriving robust estimates of standard errors. The measurement uncertainties obtained before and after the Monte Carlo bootstrapping tests were found to be in close agreement, indicating the absence of reduction flaws. The results of our magnetic field measurements, those for the entire spectrum and only the hydrogen lines are presented in Table 1, where we also include in the first rows information about the previous magnetic field measurements presented by Przybilla et al. (2016). The last row shows an additional measurement obtained on 2015 March 18 in the framework of the BOB ESO Large Programme 191.D-0255. The rotation phases presented in the last columns of Table 1 were calculated using the ephemeris determined from our period search described in Sect. 3.
High resolution spectra of CPD 3509 were obtained at the European Southern Observatory with the UVES high resolution spectrograph. The observations with an exposure time of about 40 min were obtained from 2014 December 6 to 2015 March 6 using the standard wavelength setting for dichroic mode (DIC-2 437+760) with the 0.4″ slit in the blue arm giving a spectral resolving power of and the 0.3″ slit in the red arm to achieve a resolution of . The wavelength coverage was 3730–9390 Å with a gap between 5000 and 5700 Å. Initially, we asked for 20 observations in service mode, but, unfortunately, only ten observations were executed. The summary of the UVES spectroscopic observations is presented in Table 2. The table gives the Modified Julian date at the middle of each observation, the observing date, the rotational phase, and the S/N. The S/N is given per resolution element and is measured from the spectral region around 4630 Å. All observations were phased using the ephemeris , referring to the time of the maximum positive magnetic field. Advantageously, the phase distribution of the obtained spectra appears suitable for a Doppler imaging analysis to study the surface distribution of silicon and helium, which exhibit the most distinct spectrum variability in He-strong stars.
Supplementary photometric observations were carried out at the BMT. Taken in 2014 April and May and during 2015 March 31 to April 4, the photometry covers a timeline similar to the spectroscopic and spectropolarimetric observations. Johnson and filters were used and a 33 dithering pattern of observations were obtained each night to improve the sampling and photometric accuracy. Seven photometrically invariable cluster stars (numbers 3, 44, 45, 51, 58, 97, and 110 of Baume et al. 2003) in the vicinity of CPD 3509 were used for comparison in differential photometry. They covered a magnitude range of 94 (m,m) 129 and colors from 006 (B-V) 027, thus well enclosing the target and minimizing filter effects. The rms between the comparison stars were and . The photometric data are summarized in Tables 4 and 5, presented in Appendix B. Our photometric zero values (median of the 2014 data points) for CPD 3509 are and .
3 Period determination
The results of our frequency analysis based on the longitudinal magnetic field measurements presented in Table 1 and performed using a non-linear least squares fit to the multiple harmonics utilizing the Levenberg-Marquardt method (Press et al., 1992) are presented in Fig. 1. To detect the most probable period, we calculated the frequency spectrum and for each trial frequency we performed a statistical F-test of the null hypothesis for the absence of periodicity (Seber, 1977). The resulting F-statistics can be thought of as the total sum including covariances of the ratio of harmonic amplitudes to their standard deviations, i.e. a signal-to-noise ratio. The highest peak in the frequency spectrum for the hydrogen lines not coinciding with the window function corresponds to a period of d and that for the measurements of the entire spectrum corresponds to a period of d. We note that taking into account the precision of the period determination, the difference between the two periods is not significant. On the other hand, as the achieved magnetic field measurement accuracy is higher for the set using all lines, the rotation period of 6.36255 d identified from these measurements is expected to be more reliable and is preferred in the following discussion on the surface element distribution in Sect. 4.
In Fig. 2, we present all measurements, those using the entire spectrum and those using only the hydrogen lines, phased with the corresponding rotation periods and the best sinusoidal fits calculated for these measurements. As already mentioned in the caption of Table 1, the deviating point close to the phases 0.881, respective 0.888, is caused by the unfavourable weather during the observations at that phase.
Since CPD 573509 already finished half of its main-sequence lifetime (Przybilla et al., 2016) and is already passing through the Cep instability strip, we checked the stability of the Stokes spectral lines over the full sequences of sub-exposures obtained on a time scale of tens of minutes. Along with different radial velocity shifts of lines belonging to different elements, we also detect distinct changes in line profiles taking place on time-scales corresponding to the duration of the sub-exposure sequences in the individual observations. In Fig. 3, we present the behaviour of the line profiles in individual spectral lines. The time difference between subexposures accounts for about 20 min. However, with the current data we cannot identify the periodicity of the detected variability, which is probably caused by the presence of Cep-like pulsations. Thus, future observations should focus on the careful search for periodicity and on the identification of the pulsation modes.
The periodicity derived from spectropolarimetry (and the corresponding spottiness, see below) plus the tentative Cep-like pulsations are expected to be detectable in the photometric data. Our differential photometry data in and are displayed in Fig. 4, phased according to the period obtained from the magnetic analysis using only the hydrogen lines – the behaviour is very similar for the case of the phasing based on all lines. A small periodic variability in line with the changes of the magnetic field strength is clearly seen. The dispersion of the higher-cadence data of 2015 (taken within 2/3 of the rotation period) could be interpreted in favour of the presence of Cep pulsations, but this is close to the detection limit and requires dedicated follow-up observations for confirmation. Finally, we want to note a difference in the mean and magnitudes from the present work and that of Baume et al. (2003), which are 007 and 006 fainter than our values. As the comparison stars used for the differential photometry show only an rms of 002 between both studies (data taken 20 years apart), this may imply some long-term variability for CPD 3509 as well.
4 Doppler imaging using UVES observations
As we mentioned above, the rotational phase coverage and S/N of the UVES observations is quite good and enables mapping the chemical element patterns, in particular of the most strongly varying elements He and Si. The phase coverage of the observations is not optimal, though, and four observations have been obtained close in phase. Still, the largest phase gap is 0.22, and occurs between phases 0.79 and 0.01. Other areas on the stellar surface are well covered, and the phase gaps are less than 0.15, posing no problems for the Doppler imaging technique. Tests show that phase gaps as large as about 0.28 do not significantly affect the reconstruction of the features (Rice & Strassmeier, 2000) if the S/N is good. Commonly seen problems caused by non-perfect phase coverage are the blurring of features and not being able to recover small spots (see e.g., Collier-Cameron & Unruh 1994). Chemical spots are large structures and therefore can easily be mapped with the phase coverage of the current observations. Naturally, one should be careful when interpreting the results, especially small features, in the location of the largest phase gap. In Sect. 4.2, the discussion includes the possible impact of the phase gap on the obtained maps.
Since observations with HARPSpol are usually carried out in visitor mode, we do not have high-resolution spectropolarimetric observations of CPD 3509 on our disposal. The magnetic field strength is however only a few kG, and as we show below, it is possible to use the Doppler imaging technique to obtain an approximate elemental distribution without taking into account the impact of the magnetic field on the line profiles. The He mapping is carried out using He i 4713. For Si iii, three lines located close to each other, 4552.622, 4567.840, and 4574.757 Å, were used simultaneously. For the phase calculation we employed the rotation period of d identified in our magnetic field measurements using the entire spectrum. The surface abundance maps were obtained with the INVERS7PD inversion code, which was originally developed by Piskunov, Tuominen & Vilhu (1990) and modified by Hackman, Jetsu & Tuominen (2001). INVERS7PD compares the observations to a grid of synthetic local line-profiles, which were calculated using the SPECTRUM spectral synthesis code (Gray & Corbally, 1994) and ATLAS9 stellar atmospheres by Kurucz (1993). The He i line properties are hard coded into SPECTRUM, and the Si iii line parameters were originally obtained from VALD (Piskunov et al., 1995; Kupka et al., 1999). When using the Si abundance determined by Przybilla et al. (2016), the three Si iii lines were on average well fitted, but the 4552.622 Å line was slightly too weak, while the 4574.7570 Å line appeared slightly too strong. To minimise the influence of the errors in the atomic data, we followed the standard procedure and changed the oscillator strengths of the lines to improve the overall fit (see e.g., Makaganiuk et al. 2011; Korhonen et al. 2013). The changes needed were small and allowed for fitting all three Si iii lines simultaneously. Table 3 gives the parameters used for these lines.
For the model calculations, 20 limb angles were used together with the stellar parameters determined by Przybilla et al. (2016): K, , and micro- and macroturbulence of 2.0 and 10.0 , respectively. From the inversions, we obtained a best fit inclination angle of . Notably, the Doppler imaging technique is very sensitive to the and the inclination angle values in the sense that no converging solution can be found if these values are under- or overestimated.
Przybilla et al. (2016) estimated the value of CPD 3509 to be . In the inversions, the best fit from the Si iii lines was obtained for , in excellent agreement with the results of Przybilla et al. (2016). On the other hand, for the He i line, the best fit was obtained for , which is lower than the other estimates. There are several possible reasons for this discrepancy. First, weaker metal lines are more sensitive to values and broad He i lines dominated by Stark effect are usually avoided in the determination of rotation rates. Si iii and He i lines also have different Landé factors. Further, the star is likely a Cep–like pulsator. Depending on the type of pulsations, this can have more impact on the Si iii lines rather than on the He i lines, causing stronger line broadening in the Si iii lines. We also note that similar discrepancies in the determination of values for different elements were mentioned in other works, which were using imaging methods (e.g. Yakunin et al. 2015). We also remark that SPECTRUM is verified to work at the spectral type range B to mid-M, making CPD 3509 close to the limit of the code’s capabilities. Furthermore, SPECTRUM does not use NLTE in the spectral synthesis, and He i is more sensitive to departures from NLTE than Si iii. For these reasons, the absolute abundance scales in the obtained maps are not necessarily precise. Still, this will not affect the relative abundances of the spots, nor their locations.
4.1 Magnetic field geometry
The simplest model for a magnetic field geometry is based on the assumption that the studied stars are oblique dipole rotators, i.e., their magnetic field can be approximated by a dipole with its magnetic axis inclined with respect to the rotation axis.
From the variation of the phase curve for the field measurements with a mean of G and an amplitude of G, we calculate G and G. Using the definition by Preston (1967)
we find and finally following
and employing obtained from our Doppler imaging inversions, we calculate a magnetic obliquity angle . Assuming a limb-darkening coefficient , typical for stars with K (Przybilla et al., 2016), we estimate a dipole strength of kG using the model by Stibbs (1950), as formulated by Preston (1967):
Using the parameters of the sinusoidal fit to the values resulting from only the hydrogen lines ( G and G), we obtain slightly different values for the magnetic field model: , , and kG. Given the size of the errors of the dipole strength determination, the difference in the derived values of the dipole strengths is not significant.
4.2 Abundance maps and comparison to the magnetic pole location
The He i and Si iii distributions are shown in Fig. 5, and the model fits to the observations are given in Fig. 6. The locations of the chemical spots are very similar in both maps, but the behaviour is opposite. The main concentration of He i occurs at phase 0.0, and a weaker concentration around the phase 0.5. The He i underabundance spots are located at phases 0.4 and 0.7. In contrast, the main concentrations of Si iii occur around the phases 0.4 and 0.7, and the underabundance spots are located around the phases 0.0 and 0.5. In Si iii, also a third smaller overabundance spot is seen around the phase 0.1. Some indication of a corresponding spot could also be seen around phase 0.9, but due to the phase gap in the observations at phases 0.78–0.01, we cannot verify the reality of this potential fourth Si iii spot.
Figure 7 shows the locations of the magnetic poles together with the abundance maps. The He i overabundance spots clearly occur around the magnetic poles, with a stronger concentration at the positive pole and a weaker one around the negative pole. The locations do not coincide completely in phase with the magnetic poles, and the overabundance spots seem to be shifted some 0.1 in phase away from the poles. On the other hand, Si iii shows underabundance around the magnetic poles, and the main overabundance spots fall closer to the magnetic equator. Interestingly, it seems that the Si iii concentrations occur somewhat shifted towards the negative pole, not at the magnetic equator itself. Similarly, the underabundance spots of He i are located closer to the negative magnetic pole.
The presented Doppler maps for Si and He support the dipole-dominated magnetic topology of CPD 3509. On the other hand, the presence of a third Si spot detected around the phase 0.1 and considerable abundance and field strength differences between the two magnetic poles suggest the contribution of a non-negligible higher-order magnetic multipole. Previous studies of upper main-sequence stars showed that inhomogeneous chemical abundance distributions are only observed on the surface of magnetic chemically peculiar Ap and Bp stars with large-scale organised magnetic fields. In these stars, the abundance distribution of certain elements is non-uniform and non-symmetric with respect to the rotation axis. The majority of studies of Ap and Bp stars have revealed a kind of symmetry between the topology of the magnetic field and the element distribution (see e.g. Rice, Whelau & Holmgren 1997; Yakunin et al. 2015; Hubrig et al. 2014b). Thus, the structure of the magnetic field can be studied by producing the surface element maps and measuring the magnetic field using spectral lines of inhomogeneously distributed elements separately. However, due to the low resolution of our FORS 2 spectra, we are not able to study the detailed surface magnetic field distribution. Furthermore, only the availability of high-resolution spectra in all four Stokes parameters would allow us to obtain self-consistent mapping of spots and magnetic fields by means of Zeeman Doppler imaging (ZDI; e.g. Brown et al. 1991).
Our spectropolarimetric monitoring of CPD 3509 using FORS 2 at the VLT shows the presence of an approximately dipolar magnetic field with a polar strength of 4 kG, reversing over the rotation period of 6.36 d. Using the Doppler imaging technique, we were able to constrain the inclination of the rotation axis to the line of sight, , and estimate the obliquity of the magnetic axis, . In the past it was considered neither theoretically nor observationally how the magnetic field geometry in massive B-type stars evolves across the main sequence. Thus, the analysis of the magnetic field configuration of CPD 3509 is of special interest as it is one of the most evolved He-strong stars currently known. Hubrig, North & Schöller (2007) studied the evolution of the magnetic field geometry in late B-type stars with masses with accurate Hipparcos parallaxes and definitely determined longitudinal magnetic fields. They found that the distribution of relative ages peaks at the ZAMS, with two secondary lower peaks at the relative ages around 60% and 80%. The strongest magnetic fields were found in younger stars in terms of the elapsed fraction of their main-sequence life. Further, rotation periods of late-B type stars slightly increase with age, which is consistent with the assumption of conservation of angular momentum during their life on the main sequence, without any hint of a braking mechanism (see also North & Cramer 1984; North 1985). The fact that the strongest magnetic fields are only observed close to the ZAMS can be interpreted as a magnetic field decay in stars at advanced ages. As for the magnetic field geometry, the authors detected a strong hint for an increase of obliquity with elapsed time on the main sequence. Moreover, 21% out of the studied 33 stars have magnetic phase curves fitted by a double wave, indicating that the magnetic topology in late-B type stars is frequently more complex than just a single dipole. Obviously, a comparison of the evolution of magnetic field geometries between the higher-mass He-strong stars and the lower-mass magnetic B-type stars is urgently needed to constrain the mechanism of the magnetic field generation in such stars.
There is also a large dissimilarity between the lower mass magnetic Ap stars and the He-strong stars in respect to the orientation of their magnetic axes and the period lengths. The results of modeling a small sample of Ap stars by Landstreet & Mathys (2000) implied that stars with small obliquity values of the model magnetic axis to the rotation axis, of the order of , have periods longer than 25 days. However, recent studies of the magnetic field geometries in He-strong stars do not confirm this trend: a number of fast rotating He-strong stars with periods below 2 d show low obliquities of their magnetic axes (e.g. Grunhut et al. 2012; Sikora et al. 2015; Hubrig et al. 2017).
The study of the variability of the Si and He lines showed the presence of significant chemical abundance variations across the stellar photosphere. The location of the chemical spots is roughly correlated with the topology of the magnetic field, where the main concentration of He is observed in the vicinity of the positive magnetic pole and Si underabundant spots around the poles. Only three mapping studies were devoted to He-strong stars in the past, all of them using ZDI (Yakunin et al., 2015; Oksala et al., 2015; Kochukhov et al., 2011). Yakunin et al. (2015) studied the He and O abundance distribution on the surface of the He-strong star HD 184927. Similar to our DI result, the He abundance was the highest in the vicinity of the strongest magnetic pole which in that case was the pole of positive polarity. Oksala et al. (2015) studied the distribution of He, C, Si, and Fe on the surface of Ori E. A large overabundant He spot was found to appear at the rotation phase 0.8, but the location of the spot was not correlated with the position of the magnetic poles. The minimum abundances of C, Si, and Fe were found at the same phase where He showed the largest abundance. Finally, Kochukhov et al. (2011) studied the distribution of He on the surface of HD 37776 and concluded that the He concentration is at a maximum in the regions of maximum radial field. As only very few He-strong stars were studied with ZDI, future high-resolution, high signal-to-noise spectropolarimetric observations are needed to obtain a more complete picture on how the distribution of surface chemical spots is related to the magnetic field topology. We note that CPD 3509 has a relatively low projected rotational velocity, and due to the presence of the kG magnetic field and the distinct inhomogeneous element abundance distribution, it can serve as an excellent laboratory to study various atmospheric effects that interact with the magnetic field. Further, the discovered pulsational variability on the time scale of tens of minutes has to be confirmed by future high-resolution spectroscopic time series.
Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programmes 094.D-0355 and 191.D-0255. The observations on Cerro Armazones are supported by the Nordrhein-Westfälische Akademie der Wissenschaften und der Künste in the framework of the academy program of the Federal Republic of Germany and the state Nordrhein-Westfalen. AK acknowledges financial support from RFBR grant 16-02-00604A.
Appendix A The membership of CPD573509 in NGC 3293
The star CPD573509 is located only 66 arcsec from the cluster centre of NGC 3293, whereas the angular radii of the cluster core, of the central part of the cluster, and of the whole cluster were determined as 72 arcsec, 306 arcsec, and 486 arcsec, respectively (Kharchenko et al., 2013). Although this relatively bright star is listed in the Tycho catalogue (ESA, 1997), it has no proper motion measurement in the Tycho-2 catalogue (Høg et al., 2000). Therefore, it was not included in the catalogue of stars in open cluster areas of Kharchenko et al. (2004), where membership probabilities are given for 520 open clusters of the survey of Kharchenko et al. (2005) based on proper motions and optical photometry. In the new Milky Way Star Clusters (MWSC) survey by Kharchenko et al. (2013), based on proper motions from the PPMXL catalogue (Röser, Demleitner & Schilbach, 2010) and near-infrared photometry from the Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006), CPD573509 is however listed with a 95% membership probability from proper motion and 99% from its photometry. Dias et al. (2014) also determined a high (98%) membership probability based only on proper motions from the fourth US Naval Observatory CCD Astrograph Catalog (UCAC4; Zacharias et al. 2013).
One should however note that the proper motion errors of the catalogues used for the above mentioned cluster membership studies are rather large. In case of CPD573509 they are 2.5 mas/yr in the PPMXL and between 2.2 mas/yr and 4.0 mas/yr in the UCAC4 catalogue. A significant improvement of PPMXL proper motions was recently achieved with the Hot Stuff for One Year (HSOY; Altmann et al. 2017) catalogue, which represents a new reduction of the PPMXL including Gaia DR1 data (Gaia Collaboration et al., 2016a, b; Lindegren et al., 2016). The new UCAC5 catalogue (Zacharias, Finch & Frouard, 2017) presents very accurate proper motions combining the UCAC positions with those from Gaia DR1. Because CPD573509 was not in the Tycho-2 catalogue, it was not included in the Tycho-Gaia Astrometric Solution (TGAS) of Gaia DR1. However, its HSOY proper motion, mas/yr, mas/yr, and its UCAC5 proper motion, mas/yr, mas/yr, are in very good agreement with the mean TGAS proper motion of four of the five 1 members from Kharchenko et al. (2004) that can be found in TGAS: mas/yr, mas/yr. Their standard deviation in is still about three times larger than in but was even much larger before we excluded one of the five stars as an outlier.
We do not consider the TGAS parallaxes of the few 1 members from Kharchenko et al. (2004), as they show a large spread, and since CPD573509 is not included in the TGAS. The very large standard deviation of the parallaxes of the five 1 members (0.77 mas) reduces to 0.19 mas, if we again exclude the same one outlier. However, the generally assumed additional systematic error of mas in TGAS parallaxes (Gaia Collaboration et al., 2016b) leads to a % relative distance uncertainty of a cluster at about 2 kpc distance. According to Kharchenko et al. (2005, 2013), the distance to NGC 3293 is about 2400 pc.
Concerning the radial velocity of CPD573509, there is one measurement, km s given without an error estimate (Evans et al., 2005), and slight variability … km s measured by Przybilla et al. (2016) due to the presence of spots. This is in reasonably good agreement with the mean cluster radial velocities determined by Kharchenko et al. (2005, 2013) and Evans et al. (2005), of km s, km s, and km s, respectively.
Thus the membership of CPD573509 in the cluster NGC 3293 is based on the star’s projection to the cluster core, and its available proper motion, photometry, and radial velocity. Gaia DR2, expected for April 2018, will not only provide even more accurate proper motion membership probabilities but also enable a membership study using individual stellar parallaxes of many more stars in the cluster area.
Appendix B Photometric data
- pagerange: Magnetic field geometry and chemical abundance distribution of the He-strong star CPD 573509–B
- pubyear: 2002
- Altmann M., Röser S., Demleitner M., Bastian U., Schilbach E., 2017, A&A, 600, L4
- Angel J. R. P., Landstreet J. D., 1970, ApJ, 160, L147
- Appenzeller I., et al., 1998, The ESO Messenger, 94, 1
- Baume G., Vázquez R. A., Carraro G., Feinstein A., 2003, A&A, 402, 549
- Brown S. F., Donati J.-F., Rees D. E., Semel M., 1991, A&A, 250, 463
- Collier-Cameron A., Unruh Y. C., 1994, MNRAS, 269, 814
- Dekker H., D’Odorico S., Kaufer A., Delabre B., Kotzlowski H., 2000, SPIE, 4008, 534
- Dias W. S., Monteiro H., Caetano T. C., Lépine J. R. D., Assafin M., Oliveira A. F., 2014, A&A, 564, A79
- Gray R. O., Corbally C. J., 1994, AJ, 107, 742
- ESA, 1997, The HIPPARCOS and Tycho catalogues, ESA-SP 1200
- Evans C. J., et al., 2005, A&A, 437, 467
- Gaia Collaboration, et al., 2016a, A&A, 595, A1
- Gaia Collaboration, et al., 2016b, A&A, 595, A2
- Grunhut J. H., et al., 2012, MNRAS, 419, 1610
- Hackman T., Jetsu L., Tuominen I., 2001, A&A, 374, 171
- Høg E., et al., 2000, A&A, 355, L27
- Hubrig S., Kurtz D. W., Bagnulo S., Szeifert T., Schöller M., Mathys G., Dziembowski W. A., 2004a, A&A, 415, 661
- Hubrig S., Szeifert T., Schöller M., Mathys G., Kurtz D. W., 2004b, A&A, 415, 685
- Hubrig S., North P., Schöller M., 2007, Astr. Nachr., 328, 475
- Hubrig S., Schöller M., Kholtygin A. F., 2014a, MNRAS, 440, L6
- Hubrig S., et al. 2014b, MNRAS, 440, L6
- Hubrig S., Kholtygin A. F., Schöller M., Ilyin I., 2017, MNRAS, 467, L81
- Kharchenko N. V., Piskunov A. E., Röser S., Schilbach E., Scholz R.-D., 2004, Astr. Nachr., 325, 740
- Kharchenko N. V., Piskunov A. E., Röser S., Schilbach E., Scholz R.-D., 2005, A&A, 438, 1163
- Kharchenko N. V., Piskunov A. E., Schilbach E., Röser S., Scholz R.-D., 2013, A&A, 558, A53
- Kochukhov O., Lundin A., Romanyuk I., Kudryavtsev D., 2011, ApJ, 726, 24
- Korhonen H., et al., 2013, A&A, 553, A27
- Kupka F., Piskunov N. E., Ryabchikova T. A., Stempels H. C., Weiss W. W., 1999, A&AS, 138, 119
- Kurucz R. L., 1993, CD-ROM No. 13 (Cambridge, Mass: SAO)
- Landstreet J. D., Mathys G., 2000, A&A, 359, 213
- Lindegren L., et al., 2016, A&A, 595, A4
- Makaganiuk V., et al., 2011, A&A, 525, A97
- Maeder A., Przybilla N., Nieva M.-F., Georgy C., Meynet G., Ekström S., Eggenberger P., 2014, A&A, 565, A39
- Nieva M.-F., Przybilla N., 2012, A&A, 539, A143
- North P., Cramer N., 1984, A&AS, 58, 387
- North P., 1985, A&A, 148, 165
- Oksala M. E., et al., 2015, MNRAS, 451, 2015
- Piskunov N. E., Tuominen I., Vilhu O., 1990, A&A, 230, 363
- Piskunov N. E., Kupka F., Ryabchikova T. A., Weiss W. W., Jeffery C. S., 1995, A&AS, 112, 525
- Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 1992, Numerical Recipes, 2nd edn. (Cambridge: Cambridge University Press)
- Preston G. W., 1967, ApJ, 150, 547
- Przybilla N., et al., 2016, A&A, 587, A7
- Ramolla M., et al., 2013, Astr. Nachr., 334, 1115
- Rice J. B., Wehlau W. H., Holmgren D. E., 1997, A&A, 326, 988
- Rice J. B., Strassmeier K. G., 2000, A&AS, 147, 151
- Röser S., Demleitner M., Schilbach E., 2010, AJ, 139, 2440
- Seber G. A. F., 1977, Linear Regression Analysis (New York: Wiley)
- Sikora J., et al., 2015, MNRAS, 451, 1928
- Skrutskie M. F., et al., 2006, AJ, 131, 1163
- Stibbs D. W. N., 1950, MNRAS, 110, 395
- Yakunin I., et al., 2015, MNRAS, 447, 1418
- Zacharias N., Finch C. T., Girard T. M., Henden A., Bartlett J. L., Monet D. G., Zacharias M. I., 2013, AJ, 145, 44
- Zacharias N., Finch C., Frouard J., 2017, AJ, 153, 166