Issue |
A&A
Volume 656, December 2021
|
|
---|---|---|
Article Number | A103 | |
Number of page(s) | 16 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/202140896 | |
Published online | 08 December 2021 |
A detailed understanding of the rotation-activity relationship using the 300 Myr old open cluster NGC 3532⋆,⋆⋆,⋆⋆⋆
1
Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
e-mail: dfritzewski@aip.de
2
Space Science Institute, 4750 Walnut St., Boulder, CO 80301, USA
3
Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA
4
Black Hole Initiative at Harvard University, 20 Garden Street, Cambridge, MA 02138, USA
Received:
26
March
2021
Accepted:
8
September
2021
Context. The coeval stars of young open clusters provide insights into the formation of the rotation-activity relationship that elude studies of multi-age field populations.
Aims. We measure the chromospheric activity of cool stars in the 300 Myr old open cluster NGC 3532 in concert with their rotation periods to study the mass-dependent morphology of activity for this transitional coeval population.
Methods. Using multi-object spectra of the Ca II infrared triplet region obtained with the AAOmega spectrograph at the 4 m Anglo-Australian Telescope, we measure the chromospheric emission ratios RIRT′ for 454 FGKM cluster members of NGC 3532.
Results. The morphology of activity against colour appears to be a near-mirror image of the cluster’s rotational behaviour. In particular, we identify a group of ‘desaturated transitional rotators’ that branches off from the main group of unsaturated FGK slow rotators, and from which it is separated by an ‘activity gap’. The few desaturated gap stars are identical to the ones in the rotational gap. Nevertheless, the rotation-activity diagram is completely normal. In fact, the relationship is so tight that it allows us to predict rotation periods for many additional stars. We then precisely determine these periods from our photometric light curves, allowing us to construct an enhanced colour-period diagram that represents 66% of the members in our sample. Our activity measurements show that all fast rotators of near-solar mass (F-G type) have evolved to become slow rotators, demonstrating that the absence of fast rotators in a colour-period diagram is not a detection issue but an astrophysical fact. We also identify a new population of low-activity stars among the early M dwarfs, enabling us to populate the extended slow rotator sequence in the colour-period diagram.
Conclusions. The joint analysis of chromospheric activity and photometric time series data thus enables comprehensive insights into the evolution of the rotation and activity of stars during the transitional phase between the Pleiades and Hyades ages.
Key words: stars: chromospheres / open clusters and associations: individual: NGC 3532 / stars: late-type / stars: activity / techniques: spectroscopic / techniques: photometric
Full Table 1 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/656/A103
© ESO 2021
1. Introduction
The magnetic dynamo of cool stars is believed to be driven jointly by stellar rotation and convection. Although the magnetic field plays a vital role in the creation of various stellar activity phenomena, it is observationally less accessible than these phenomena. Because of this difficulty, magnetic activity is typically accessed, especially in larger samples, through certain indicators, such as ‘excessive’ chromospheric emission. Such emission was first observed by Eberhard & Schwarzschild (1913) in the reversal of the Ca II H and K lines in objective-prism spectra of bright stars and was immediately connected to sunspot-like features on other stars.
For the majority of stars, one observes decreasing activity with decreasing rotation rate. This correlation between the chromospheric emission and stellar rotation rate can be traced back to Kraft (1967) and Skumanich (1972) and similar studies of the time evolution of different activity tracers. Noyes et al. (1984) made an important advance, showing that the chromospheric activity of a sample of field stars is particularly well correlated with a mass-normalised version of the rotation period called the Rossby number, Ro = Prot/τc (with Prot representing the stellar rotation period and τc the convective turnover timescale). The corresponding reduction in the dimensionality of the problem1 allows all cool stars to be considered as part of the same continuum despite differences in their evolutionary timescales. Indeed, this rotation-activity correlation has repeatedly been shown to be fundamental in understanding stellar activity. It is observed not only in chromospheric activity, but also in the coronal X-ray activity of solar-like (Pizzolato et al. 2003) and fully convective stars (Wright & Drake 2016; Newton et al. 2017).
However, the dimensional reduction at the heart of the relationship also conceals some of the underlying activity behaviours of cool stars. Certain stars could potentially be over- or under-active as compared with otherwise similar stars. Likewise, a small change in a single underlying parameter, such as mass or rotation period, could result in a substantial change in activity. Such dependences are best addressed in open clusters, where the coevality and homogeneity of the underlying population tame the parameter dependences, permitting more detailed examination than with field star samples.
A number of tracers are available for probing stellar activity in different layers of outer stellar atmospheres. In the photosphere, one observes starspots (e.g. Strassmeier 2009; Basri & Nguyen 2018) and flares (e.g. Pettersen 1989). In the chromosphere, emission in certain spectral lines (e.g. Hα, Ca II H and K, and the even bluer Mg II H and K) are tracers of heating beyond radiative equilibrium (Hall 2008). Finally, in the corona, the outer layer of the stellar atmosphere, emission in soft X-rays is the most prominent activity tracer (Schmitt 1997; Güdel 2004).
The early work of Eberhard & Schwarzschild (1913) on the Ca II H and K lines was followed up, most famously, by extensive monitoring of the chromospheric activity variability of the same lines by the Mt. Wilson survey (Wilson 1978; Baliunas et al. 1995). However, the well-known Fraunhofer H and K lines are not the only spectral lines of Ca II that are affected by chromospheric activity. In this work, we analyse the Ca II infrared triplet (IRT) lines at 8498 Å, 8542 Å, and 8662 Å.
These lines were first used by Linsky et al. (1979) to probe stellar activity because they are formed by the same upper energy level as the H and K lines (Linsky et al. 1979; Foing et al. 1989) and hence can be used in a similar manner. Their position in the near infrared (photometric I band) enables observations of cool main sequence stars, which have their peak emission in the near infrared, out to larger distances. In fact, the IRT lines have not only been used to study field stars (e.g. Linsky et al. 1979; Žerjal et al. 2013) but also stellar activity in open clusters2.
Yamashita et al. (2020) studied the youngest stars in open clusters and star forming regions and found that most of the T-Tauri stars in their sample show magnetically induced activity similar to that seen in older stars. Marsden et al. (2009) observed the open clusters IC 2391 and IC 2602, which host F- and G-type zero-age main sequence (ZAMS) stars and pre-main sequence stars of lower mass. Soderblom et al. (1993) studied the solar-like stars in the ZAMS Pleiades cluster, and Jackson & Jeffries (2010) observed M dwarfs in NGC 2516. The oldest open cluster to date with stellar activity measured from the IRT lines is the 220 Myr old NGC 6475 (M 7; James & Jeffries 1997). Finally, Soderblom & Mayor (1993) measured stellar activity with the IRT lines for stars in the sparsely populated Ursa Major Moving Group. However, its age is not well constrained (Jones et al. 2015 and references therein). In this work we extend the age sequence with the 300 Myr old open cluster NGC 3532 and analyse the rotation-activity connection in detail using its rich coeval population.
NGC 3532 is a very populous open cluster that is potentially capable of becoming a benchmark object. It is located in Carina, at a distance of 484 pc, and is somewhat embedded in the Galactic disc. For a recent detailed overview of NGC 3532 and related work, we refer the reader to the introduction of Fritzewski et al. (2019, hereafter F19). Our prior efforts on NGC 3532 include a spectroscopic membership study (F19) and, in a companion paper, measurements of the stellar rotation periods of its cool star members (Fritzewski et al. 2021, hereafter F21rot). Here, we present measurements of the chromospheric activity in this cluster, examine their patterns, especially with respect to the rotation periods, and use this knowledge to determine a set of activity-guided photometric rotation periods.
This paper is structured as follows. In Sect. 2 we present our observations, followed by the measurement of the chromospheric emission in Sect. 3. In Sect. 4 we compare the chromospheric activity to the photospheric one. Afterwards, we combine the chromospheric activity with photometric rotation periods of stars in NGC 3532 in Sect. 5 to study the colour-rotation-activity connection in detail. Here, we construct a colour-activity diagram (CAD) and gain insight into the rotational state of stars without photometric rotation periods. In Sect. 6 we construct and analyse the first rotation-activity diagram for NGC 3532. Furthermore, we derive activity-informed photometric rotation periods for members of NGC 3532. The rotational implications of these new periods are analysed in detail in the companion paper (F21rot).
2. Observations, data reduction, and cluster membership
2.1. Spectroscopy
We observed NGC 3532 with the AAOmega multi-object spectrograph on the 3.9 m Anglo-Australian Telescope (ATT) on 10 and 11 March 2017 primarily from the viewpoint of the membership study presented in F19. The spectrograph is fed by the 392-fibre 2dF fibre positioner (Lewis et al. 2002), for which we used three different fibre configurations to observe 1060 stars with the 1700D grating (R = 10 000). Two of the fibre setups targeted the brighter stars (Ic = 12 − 15) with a total exposure time of 30 min each, each split into three sub-exposures per spectrum. The final fibre configuration targeted the fainter stars (Ic = 15 − 17) and was exposed for 120 min.
The data were reduced with the 2dFr pipeline (AAO Software Team 2015, version 6.28), as described in F19; we use the final data products from that work for the present analysis. Before measuring the equivalent widths (EWs), we normalised each spectrum with a fifth-order Chebyshev polynomial. In addition to the spectra, we use the measured radial velocity shifts (in pixel) from F19.
2.2. Photometry
The time series photometry used in this work was obtained between 19 February 2008 and 1 April 2008 from the Cerro Tololo Inter-American Observatory (CTIO) with the Yale 1 m telescope operated by the SMARTS consortium. We observed eight different fields in NGC 3532 with three different exposure settings in both the V (120 s) and Ic (60 s and 600 s) filters. After standard data reduction, we performed point-spread function photometry with DAOPHOT II. The observations, data reduction, photometry, time series analysis, and rotation period determinations are described in the companion paper (F21rot) where these observations are the centrepiece.
For the standardised photometry, we used the optical BV(RI)c cluster survey of Clem et al. (2011) (C11) and combined it with 2MASS (Skrutskie et al. 2006) Ks for the primary colour (V − Ks)0 of our analysis. The photometry was de-reddened with the mean reddening towards NGC 3532 of E(B − V) = 0.034 mag.
2.3. Membership and multiplicity
The spectroscopic observations presented here were originally obtained to establish a radial velocity membership list for NGC 3532. In F19, we presented a set of 660 exclusive members of NGC 3532 based on our own observations as well as survey and literature radial velocities. Here, we analyse the stellar activity of the exclusive members within our spectroscopic dataset. In Fig. 1, we show the spectroscopically observed members of NGC 3532 and highlight different subsets of stars with measured rotation periods (‘rotators’ hereafter).
![]() |
Fig. 1. Colour-magnitude diagram of the spectroscopically observed members of NGC 3532. In orange we show the rotators from Fritzewski et al. (2021) and in light blue from this work. We were unable to determine rotation periods for stars marked in green. Single stars are indicated with circles and candidate binaries from photometry and astrometry with squares. |
Unidentified close binaries could potentially skew our results because binarity can influence a star’s magnetic and rotational history (e.g. Meibom & Mathieu 2005; Meibom et al. 2007). To assist the analysis of stellar activity and rotation, we searched for potential binaries among the cluster members in the companion paper F21rot. There, we classified cluster members as binaries if their distance to the main sequence locus is ΔG > 0.25 mag (photometric criterion) or their re-normalised unit weight error (RUWE) from Gaia early data release 3 (Lindegren et al. 2021) satisfies RUWE > 1.2 (astrometric criterion, Belokurov et al. 2020). In combination, we find 151 potential binaries (22%). These are marked in all relevant figures with squares.
3. Chromospheric emission measurement
Previous studies of stellar chromospheric emission in open clusters in the IRT lines (e.g. Marsden et al. 2009) used spectra of inactive stars and subtracted them from their target spectra to obtain the chromospheric emission. However, this approach is not possible with our dataset because the observations were initially planned as a pure membership survey without the (intentional) inclusion of inactive field stars. The number of true main sequence non-members (in radial velocity, proper motion, and parallax) is too small to build a grid of reference spectra. Stars that pass one or more of the membership criteria could still potentially be members, meaning young and active, and can certainly not be considered inactive. The empirical estimates of the basal flux of inactive stars (Busà et al. 2007; Martin et al. 2017) were also excluded because of their incompleteness with respect to our sample.
Given the lack of suitable reference spectra, we followed a different approach and measured the EW in a 1 Å-wide window in each of the three line cores without prior subtraction of a reference spectrum. We note that a similar approach is often chosen for the Ca II H and K lines (e.g. Baliunas et al. 1995) because the tight 1 Å (triangular) spectral window includes only the line core that contains the signature of the stellar activity – the chromospheric emission.
The resolution of the AAOmega spectra is too low to simply integrate over a 1 Å bandpass. Therefore, we linearly interpolated between the pixels in a 2 Å-wide window and oversampled by a factor of ten. Thereafter, we placed the 1 Å bandpass according to the stellar radial velocity shift (Fig. 2) and integrated over it to obtain the EW in terms of λ and the measured and continuum fluxes
![]() |
Fig. 2. Spectral cutouts for example stars of different chromospheric activity (orange) and their corresponding PHOENIX spectra. Each panel is centred on one IRT line and shifted to the rest frame. The line centre is marked with a dotted line and the limits of the 1 Å wide integration window with dashed lines. The spectra are grouped into low-activity (top) and high-activity stars (bottom). The latter feature prominent fill-ins of the line cores. Earlier spectral types are higher. The topmost spectrum is from a G2 star, and the lowest is from an M 3 star. All six stars are marked with squares in the colour-activity diagram in Fig. 5a. |
We defined the continuum flux Fcont. with an error-weighted least-squares linear fit to the continuum outside of the IRT lines. The fit was evaluated at the position of the integration window, creating a local pseudo-continuum. We followed the approach presented by Cenarro et al. (2001) and used their spectral windows. Cenarro et al. (2001) explain that such a fit improves the continuum determination, in particular in the near infrared with many sky lines and varying signal-to-noise ratio3.
In contrast to the Ca II H and K lines, the Ca II IRT lines have a significant photospheric component. We are unable to remove this flux using inactive stars observed with the same instrument, as explained earlier. Hence, we measured the EW of the IRT lines in the same manner in down-sampled model spectra4 from the PHOENIX spectral library (Husser et al. 2013). These spectra do not contain a chromospheric component. The PHOENIX EWs were linearly interpolated over the effective temperature and act as a substitute for our inactive ‘reference stars’.
To connect this EW scale to our observation, we estimated the effective temperature for each spectrum of the NGC 3532 members from the relations by Mann et al. (2015) through the (V − Ic)0 colour index and subtracted the photospheric flux (in EW) of the corresponding PHOENIX model from the measured EWs. The sum of the EW differences of all three lines is defined as the chromospheric emission in the IRT. It was converted to an excess flux () through the empirical flux scale of Hall (1996)5. In addition to the surface flux in the IRT, we calculated the bolometric flux Fbol, assuming a black body, with the Stefan-Boltzmann law. The flux ratio in the IRT follows as (cf. Marsden et al. 2009):
For each data point, we propagated the resulting uncertainty from the observations. In addition, we also include the uncertainty in the estimated effective temperature (80 K, Mann et al. 2015). When subtracting the photospheric flux, we also calculated the chromospheric emission of an 80 K warmer reference and include the difference in the uncertainty. The temperature uncertainty was also considered in the calculation of the bolometric luminosity.
In our estimate, we ignored the uncertainty associated with the Hall (1996) flux scale for two reasons. Firstly, when considering equally massive stars in the colour-activity diagram, the conversion factor from equivalent width to chromospheric flux is constant. Hence, no matter how large the uncertainty is, the scatter among the sequence is determined by the measurement of the EW and the intrinsic spread. Secondly, we do not find evidence in either the colour-activity diagram or the rotation-activity diagram that the observed spread is an effect of the flux conversion. Hence, we conclude that the uncertainty estimated by Hall (1996) is too large6 and can even be ignored in the case of a coeval group, such as an open cluster7. The IRT excess-EW, the corresponding flux , and the chromospheric emission ratio log
, including their respective uncertainties, are provided in Table 1.
Chromospheric activity measurements, including uncertainties and auxiliary data, for members of NGC 3532.
Among the spectroscopically observed stars in the AAO campaign, we find 454 radial velocity members. For each spectrum, we are able to measure the chromospheric emission, providing one of the largest sets of chromospheric activity measurements in any open cluster to date8. Out of these 454 stars, 387 stars fall within the field of our photometry, meaning that we have light curves for them. We were able to obtain photometric rotation periods for 142 of them, as described exhaustively in the companion paper F21rot.
These periods enable us to probe the rotation-activity relation in this 300 Myr old population and to understand the connection between rotation and activity in a mass-dependent manner. Furthermore, we are able to increase the number of photometric rotation periods through an activity-informed search in the time series data among the 245 (=387−142) remaining stars with light curves but without directly measured prior periods. Prior to these analyses (Sects. 5 and 6), we use the light curves of all stars in the spectroscopic sample to probe the connection between photospheric and chromospheric activity.
4. Chromospheric against photospheric activity
4.1. Photospheric activity
Stellar activity manifests itself in various forms and can be observed in all layers of the stellar atmosphere. In the photospheric layer, starspots are a widely observed magnetic feature of solar-type stars (e.g. Strassmeier 2009), and are often probed using photometric time series observations. The simplest measure of such photospheric activity is the variability amplitude of the light curve, which we can obtain from our own photometric time series. We use the difference between the fifth and the ninety-fifth percentile of the sigma-clipped light curve (Basri et al. 2011). We note that a similar measure is called Rrng in Basri et al. (2011) and measured as a flux difference. In our work we measure it as a magnitude difference and hence call it VP95 to indicate the difference9.
Because the variability amplitude is a proxy for the photospheric stellar activity, it appears to have similar dependences on stellar rotation and magnetic field as chromospheric and coronal activity measures. The left panel in Fig. 3 shows the variability amplitude against the intrinsic colour for NGC 3532. (We note that the (V − Ks)0 colour of a Sun-like G2 star is approximately 1.5). We removed the bluest stars of our sample because they are overexposed. In data optimised for such stars, the activity is observed to be low, in keeping with the thinning of convection zones as the Kraft break is approached. In contrast, the higher variability among the lowest-mass stars is intrinsic, and an expression of their faster rotation.
![]() |
Fig. 3. Left: V band variability amplitude (VP95) for stars in NGC 3532 (orange), compared with those of the zero-age main sequence open cluster NGC 2516 (black) for both stars with periods (squares) and those without (circles). Stars in NGC 3532 generally have lower variability, as seen by the mean amplitudes expressed by the solid lines. Centre: VP90 against Rossby number Ro for the rotators in NGC 3532 and NGC 2516. The distribution can be seen to resemble the activity-rotation relation previously known from chromospheric and coronal activity, and it shows that the photospheric activity amplitude has a similar structure. Right: histogram of the variability amplitudes for the members of the two open clusters, showing that the high-amplitude stars of NGC 3532 are far less numerous in comparison with those of NGC 2516. The dotted grey histogram is restricted to NGC 2516 stars with (V − Ks)0 < 4 and confirms that the high-amplitude peak in NGC 2516 is mostly attributable to the fast-rotating M dwarf members, a population absent from our NGC 3532 data. |
For comparison, we added the variability amplitudes of members in the 150 Myr-old (ZAMS) open cluster NGC 2516 (F20) and show the mean amplitude in 0.5 mag-wide colour bins. We find that the NGC 2516 stars generally have higher photospheric activity levels as compared with stars in NGC 3532 except for the very variable M dwarfs, where no difference is perceptible. This shift towards lower variability and activity with increasing age is expected from the stellar spin-down over time and the corresponding decrease of the magnetic field strength.
Indeed, the lower activity of stars in NGC 3532 is exclusively driven by the slow rotators, as can be clearly seen in the central panel of Fig. 3. Here, we show the variability against the Rossby number for the periodic rotators in both open clusters. This diagram resembles the well-known rotation-activity correlation from Ca II H and K chromospheric- and coronal activity (e.g. Noyes et al. 1984; Pizzolato et al. 2003). The stars with the highest variability amplitudes are found in a flat region among the fastest rotators in a mass-independent way – the saturated regime. In this regime, no difference between stars in NGC 2516 and NGC 3532 is observed, with stars from both clusters occupying the same elevated region across a wide swath in Rossby number. The situation is different in the regime of larger Rossby numbers (Ro ≳ 0.06), where the variability amplitudes are well correlated with the Rossby number, a behaviour seen in other activity indicators such as coronal X-rays. This diagram clearly shows that the lower photometric variability of stars in NGC 3532 is a function of the Rossby number and therefore a sign of the slower rotation compared to NGC 2516.
The lower overall variability can also be seen in the histogram in the rightmost panel of Fig. 3. The distribution for NGC 3532 has the majority of stars between 14 mmag and 40 mmag (median amplitude 32 mmag), with a single-peaked distribution. For the younger NGC 2516, in contrast, we find a double-peaked distribution between 20 mag and 120 mmag (median 67 mmag). The double-peaked structure can be traced to the prevalence of many fast-rotating M dwarfs in NGC 2516. If we restrict the colour range to (V − Ks)0 < 4 (thereby removing the M dwarfs), the peak at the higher amplitudes disappears. However, even in this restricted sample the median amplitude for NGC 2516 is higher than that for NGC 3532.
4.2. Connection between photospheric and chromospheric activity
Because the photospheric activity shows the same structure as is observed in chromospheric activity, we expect the two to be correlated. In Fig. 4, we therefore show the photometric variability amplitude in the V band VP90, V against the chromospheric emission ratio log . A clear trend (correlation coefficient r = 0.665) is visible in this comparison, with higher photometric variability visible at higher chromospheric flux. The large scatter in the trend can potentially be traced to short-term activity observed in the light curves10. Additionally, the inclination angle is expected to influence the variability amplitude. We note that the spectroscopic and photometric data are not even nearly contemporaneous, but were obtained nine years apart. Despite this significant time difference, the chromospheric and photospheric activity are well correlated, not unexpected because of the underlying common origin – the stellar magnetic field and the stellar rotation, both of which evolve on much longer timescales.
![]() |
Fig. 4. Photometric variability amplitude in the V band VP95 against chromospheric emission ratio log |
In related work, Zhang et al. (2020) have compared variability amplitudes of stars in the Kepler field to their Ca H and K chromospheric emission. They found similar results but were able to compare stars of different spectral types because their sample includes stars of various activity levels caused by the inhomogeneous stellar ages in the Kepler field. This kind of comparison is not possible in our coeval dataset because all F and G dwarfs have very similar chromospheric activity. Radick et al. (2018) found a similar correlation between chromospheric activity and photometric activity for young and active stars. Other work on the temporal correlation between photospheric and chromospheric activity found anti-correlations (Morris et al. 2018) and a comparison between two chromospheric activity indicators found a correlation only for very active stars (Rauscher & Marcy 2006).
5. Chromospheric activity and stellar rotation in the colour-activity diagram
Because chromospheric activity is closely linked to the stellar rotation rate, we study both properties jointly in Fig. 5 to emphasise this connection. Our large number of measurements and our provision of rotation periods (instead of v sin i measurements) will be seen to allow more detailed connections between IRT activity and the relevant quantities than was possible in the most notable prior cluster studies to date, those of the Pleiades (Soderblom et al. 1993) and NGC 6475 (James & Jeffries 1997).
![]() |
Fig. 5. Interconnection between chromospheric activity, rotation, and mass. The colour-activity diagram (a) connects the colour-period diagram (b) with the rotation-activity diagram (c). a: chromospheric emission ratio log |
In panel a, we display the measured activity log against the intrinsic (V − Ks)0 colour to investigate its mass dependence11. The second diagram in panel b features the colour-period diagram (CPD) from the companion rotational study. The final panel c shows the rotation-activity diagram for NGC 353212. Stars with common features in panel b have been grouped together and colour-coded so that their locations in all three panels can be recognised. This enables the detailed analysis of the rotation-activity connection even in the colour-activity diagram (CAD), which otherwise does not contain the rotational information.
The CAD in Fig. 5a is analogous to a colour-magnitude diagram (CMD) or a CPD, and it shows the chromospheric emission ratio against the intrinsic colour. Our CAD also contains a large number of NGC 3532 cluster members for which we were unable to derive rotation periods (unfilled symbols). As in a CMD or CPD the (main sequence) open cluster members form narrow bands in the CAD, but notably, what was previously believed to constitute a ‘dispersion’ in activity values here appears to resolve into distinct sequences and gaps that mirror the behaviour in the CPD. The most prominent of these sequences runs diagonally in the CAD from low stellar activity levels in the warmest (F) stars to high levels among the M dwarfs. From this main feature a branch of higher-activity stars splits off at (V − Ks)0 ≈ 2.1. At the red end of the CAD, among the lowest-mass stars, we observe a large spread in the activity levels. A similar spread is observed in the stellar rotation periods. We discuss each group of stars separately below.
5.1. Unsaturated slow rotators
With the help of the colour-coding, the main feature in the CAD can easily be associated with the slow rotator sequence in the CPD (Fig. 5b). There is a very well-defined sequence of FGK dwarfs with log ≲ −4.1 in the CAD that mirrors the clean slow rotator sequence of such stars in the CPD (dark grey). We therefore call these stars the ‘unsaturated slow rotators’ (USRs). It is strengthened considerably by additional stars for which we have been unable to derive rotation periods (open symbols). In agreement with prior work on the rotation-activity connection, the slow rotators clearly correspond to the lowest-activity stars in each mass bin in the open cluster. In keeping with this idea, we identify the M dwarfs with log
≲ −4.2 as the continuation of the sequence of unsaturated stars. These will be found later to lie on the extended slow rotator sequence in the CPD. At this stage, the classification of the M dwarfs marked in purple is ambiguous. We discuss them separately below, and argue below that they should not be considered as USRs because of their increased activity (although short of saturation level) and shorter rotation periods. We refer to them here as ‘desaturated slow M dwarfs’ (DSMs). The significantly lower activity of the warmest (F and G) stars relative to the K dwarfs, is also unsurprising because they have far shallower surface convection zones.
Three stars (C11: 70134, 72069, 121200, marked with black diamonds in Fig. 5a) fall below the well-defined unsaturated sequence. Given the correlation between rotation and activity, those stars might be expected to rotate much slower than other equally massive stars in the cluster, implying that they are likely field stars of an older age. However, all three stars are bona fide cluster members in radial velocity, proper motion, parallax, and photometry. Unfortunately, their light curves and periodograms are too noisy to determine their photometric rotation periods.
We find only one star with an elevated activity level (log = −4.1 at (V − Ks)0 = 1.7) that places it above this well-defined sequence for the F and G dwarfs in the CAD. Our measured rotation period, however, places it in agreement with the other slow rotator stars in the CPD (rather than much faster, for instance). This star is a photometric binary, and it is possible that binary-induced activity is responsible for its elevated position in the CAD. Alternatively, we could have captured it by pure chance during a phase of increased stellar activity.
5.2. Desaturated transitional rotators and activity gap
We have identified an unexpectedly tight sequence of G-K stars with activity above the USR level, but below the saturated level. As can be seen from Fig. 5a, near (V − Ks)0 ≳ 2.0, log ≳ −4.1 a branch of higher-activity stars emerges from the unsaturated sequence. Cross-comparison of this higher-activity branch with the CPD reveals these stars to be associated with the stars in transition from fast to slow rotation (marked in orange). In fact, the association is one-to-one for stars with measured rotation periods. This branch of higher-activity stars reaches the saturation level13 at (V − Ks)0 ≈ 3.0, at which point all of its stars are classified as fast rotators (light grey) in the CPD. We propose the name ‘desaturated transitional rotators’ (DTRs) to capture their transitional rotation-activity behaviour. There are four binaries (squares) in this group, but the majority (seven) of these stars show no evidence of binarity (circles).
In the gap between the two sequences in the CAD discussed above, we find a small number of stars (marked in green) that correspond exactly to the stars in the rotational gap of NGC 3532. With this knowledge, we can clearly define the rotational gap in the CPD separating the sequences. In fact, the correspondence with the CPD, where these stars are in transition from fast to slow rotation, tells us that these stars in the CAD are in transition from the high-activity branch to the low-activity one. In keeping with prior terminology for the rotational gap, and also to reflect the intermediate activity level of these stars, we refer to them as ‘desaturated gap stars’.
We note that the DTRs have a higher binarity fraction compared to USRs of the same colour. Yet, as the rotational analysis in the companion paper shows, these stars do not spin faster due to binary tidal interactions (as for example in Praesepe; Douglas et al. 2019) but are the remainder of the initial fast rotators that emerge from the star formation process. In the companion paper (F21rot), we interpret the higher binarity fraction among the fast rotators with the shorter disc lifetime in binary systems. Even so, only a minority of the DTRs show any sign of binarity, and the higher activity with faster rotation appears not to be directly caused by binarity.
The absence of additional (elevated) stars bluewards of the transitional DTR sequence in the CAD (orange) is also notable, and very important in understanding angular momentum evolution because it clearly shows that even stars without determined photometric rotation periods fall onto the low activity and slow rotator branch in this colour range. Conversely, this correspondence assures us that all fast rotators down to a certain mass limit have already transitioned to the slow rotation phase. The lack of fast-rotating stars in the CPD for solar-mass stars is therefore astrophysical, and not an observational bias caused by small variability amplitudes.
Given the accurate correspondence between the CAD and the CPD, the branching-off point of the DTRs could potentially be used to age-rank young open clusters with high relative sensitivity. As shown in the companion rotational study, the transition from fast- to slow rotation happens on short timescales (tens of megayears), and hence is very age-sensitive. Using a CAD would correspondingly allow potentially accurate age-ranking to be obtained from a single multi-object spectroscopic observation14, rather than a time series, as needed for photometric rotation periods.
5.3. Diversity in activity of low-mass stars
As one approaches the lower-mass stars towards late K and early M dwarfs ((V − Ks)0 ≳ 3.2), the activity of the stars is more diverse, and the structure of the CAD becomes correspondingly more complicated. The diversity is not related to multiplicity as seen from Fig. 5a where the fraction of likely binaries in this regime is seen to be comparable to the USRs. We discuss this mass regime (indicated by the red outlines in Fig. 5a) sequentially from high to low activity.
At the top of the CAD, in the saturated regime (log > −3.8), we find two corresponding groups in the CPD: the fast rotators (light grey) and a group of stars rotating slightly more slowly (yellow). Despite the latter stars rotating more slowly than the former, their chromospheric activity is still saturated. The large number of M dwarfs with saturated activity and as yet undetected periods must therefore all belong to one of these two groups. We infer that we likely did not have the photon sensitivity to directly detect periods for them.
In a band with −3.8 > log > −4.0, we find a group of lower-activity stars that are also significantly slower rotators than the saturated activity ones. These are grouped tightly together in both the CAD and the CPD, where they are marked in purple. We call these the ‘desaturated slow M dwarfs’ (DSMs) to reflect both their reduced activity levels as compared to their saturated counterparts, and their location in the CPD, near the slow rotator sequence, if not quite on it. Because of these joint intermediate behaviours, we suggest that these stars are distinct from the low-activity M dwarfs (discussed immediately below), and that they should not be considered as a continuation of the USR sequence in the CAD. That continuation instead is provided by the sizeable group of low-activity M dwarfs below them in the CAD.
Finally, we find a group of M dwarfs with log < −4.0. They form a sequence of declining activity with declining mass (below the red dotted line). This sequence connects smoothly to the USRs (dark grey), and we identify these stars as the continuation of the USR sequence into the M dwarf regime. At this stage of the discussion, we have no rotation periods for any of them, and thus no counterparts for them in the CPD. We have likely ignored their rotation periods with our demand for a clearly recognisable rotation period. (We reappraise this in the following section and derive their photometric rotation periods, thereby justifying the classification here).
A sparsely populated gap appears between the two sequences traced among the M dwarfs. Our data are not numerous enough to provide sufficient insight into the reality of such a gap between the low-activity M dwarfs and the more active stars. However, we note that a gap is present between the fast rotators and the slow rotators among the M dwarfs in the ZAMS open clusters Pleiades and NGC 2516 (Rebull et al. 2016; Fritzewski et al. 2020), a hint of which is found in the CAD. Ca II H and K and Hα observations of the Pleiades show a similar feature (Fang et al. 2018). It should also be noted that Fang et al. (2018) have questioned the Pleiades membership of the least active stars despite the very slowly rotating stars already known in the Pleiades (Rebull et al. 2016).
A gap between slow and fast rotators also exists in field M dwarfs. The chromospheric emission measurements (LHα/Lbol) presented by Newton et al. (2017) show two distinct groups of high and low activity. They correspond to rapidly and slowly rotating stars that are separated by a gap (Newton et al. 2016). With the gap also present in the field star population, it should be well defined and very sparsely populated in open clusters. However, more observations of mid-M stars are needed to constrain its extent.
6. Rotation-activity relationship and additional rotation periods
Having discussed the mass dependence of the observed stellar activity in the CAD, we move on to investigating the dependence of the IRT activity on rotation. The wide difference in activity between stars with the same rotation periods (for instance the USRs and the DTRs with similar rotation periods of 4 − 6 d) shows that this dependence is not direct, and that a star’s mass (or suitable proxy) has to be considered as well. Such comparisons are therefore typically performed using a mass-normalised version of the rotation period called the Rossby number, to compensate for the variation in activity with stellar mass. The first such diagram was devised and presented by Noyes et al. (1984), and has since been constructed not only for the chromospheric emission of solar-type stars, but also for their coronal emissions (e.g. Pizzolato et al. 2003) and also for the activity of fully convective stars (e.g. Wright et al. 2011).
However, a drawback of the Rossby number is that it depends on the convective turnover timescale, which is not directly accessible from observations, and one has to rely on stellar models or empirically calibrated relations. For this work, we use the relation between Teff and the convective turnover timescale from Barnes & Kim (2010) to derive it for our stars. The effective temperature has already been estimated in Sect. 3 through the relationship given by Mann et al. (2015). The uncertainty in Ro was estimated from the period uncertainty and the corresponding calculation for τc using offsets of ΔTeff = 80 K. As with the uncertainties for log , this approach only includes the measurement and not model-dependent uncertainties. The latter typically differ by only a scaling constant, implying that the underlying power-law dependences are unchanged.
6.1. The rotation-activity relationship
Upon transformation to Ro, we find that all stars with measured rotation periods follow a remarkably tight rotation-activity relationship (Fig. 5c). The separate groups and even sequences of stars in the CPD and CAD collapse onto a single relationship, obscuring their diverse behaviours in those diagrams. Ours echoes similar relationships seen in chromospheric (Noyes et al. 1984) and coronal activity (Pizzolato et al. 2003). Two regimes known from prior work are instantly recognisable. In the saturated regime of low Rossby numbers, the chromospheric emission is independent of the rotation period and has a value of log ∼ −3.7. After a curved transition region (in this log–log representation) in the vicinity of Ro ≈ 0.1, the chromospheric emission declines steeply with Ro to define the unsaturated regime, where activity depends strongly on rotation.
Our prior division of stars in the CPD into distinct groups provides further insight into the different activity regimes. While it is unsurprising that the warmest of the slow rotators are the most evolved from an activity standpoint, and that the fastest M dwarfs are all saturated and thus least evolved, we see from the colour-coded divisions made earlier that there is indeed an order in which the activity evolves. The first stars to evolve after the unsaturated slow rotators (dark grey) are the desaturated slow M dwarfs (purple) and the activity gap stars (green), finally followed by the saturated intermediate rotators (yellow), and then the saturated fast rotators (light grey). Perhaps unsurprisingly, the desaturated transitional rotators (orange) which span a large range of Rossby number, stretch in the rotation-activity diagram all the way from the unsaturated to the saturated stars.
The correlated, unsaturated, part of the rotation activity relation is sometimes described as a power-law (e.g. Pizzolato et al. 2003; Wright et al. 2011), and at other times with a log-linear function (Mamajek & Hillenbrand 2008). We find the log-linear function to follow our data better than a power-law. This suits our goal of effectively locating the stars without measured rotation periods in the activity-rotation space. We fit the correlation in the unsaturated regime with a linear function of the form log = a + b Ro using Orthogonal Distance Regression15 (ODR; Boggs & Rogers 1990) from the SCIPY package. We perform the fits over two ranges, once for the slow rotators alone (defined as having Ro > 0.11) and once for all unsaturated stars (those for which Ro > 0.06). We show both results in the rotation-activity diagram in Fig. 6.
![]() |
Fig. 6. Rotation-activity diagram with two possible log-linear fits to two regions of interest: all unsaturated stars (blue, Ro > 0.06) and only the slow rotators (orange, Ro > 0.11). |
For the slow rotators alone, we find
and including all unsaturated stars the relation is
The root of the differing slopes for both fits might be found among the least active stars. Due to small numbers, we find the scatter among these stars to be smaller compared to other regions of the rotation-activity diagram. Yet the fit to the slow rotators (Ro > 0.11) matches these data points because it has fewer constraints among the more active stars.
6.2. Lower activity among M dwarfs?
All stars, with only one exception, are seen to follow the activity-rotation relation. That star (ID 707560, C11: 210716) falls 0.3 dex (well below) the saturated branch. It has a short rotation period (Prot = 2.76 d) that is associated with a clear peak in the periodograms, and shows no alias periods with significant power. The light curve and the periodograms can be found in Fig. A.1. Given such properties, we would normally expect it to fall into the saturated regime.
However, the star is peculiar in more than its activity; it is also the lowest-mass rotator in our sample. Unfortunately, we were not able to derive rotation periods for other stars of similar mass and activity. Indeed, a large spread in activity can be found among the least-massive stars in the CAD. Given the connection between rotation and activity, we would also expect a correspondingly large spread in periods. However, 707560 is already among the least active stars in the fully convective regime ((V − Ks)0 ≳ 4.6), with a rotation period significantly faster than for other early M dwarfs of similar activity.
In conclusion, we cannot specify whether this star has an unusually low activity or whether such low-mass stars simply do not follow the rotation activity relation of the higher-mass stars. Jackson & Jeffries (2010) have found a decline of the activity, measured in the IRT, with colour for mid-M dwarfs in NGC 2516. In their work mid-M dwarfs populate an unstructured rotation-activity diagram. Such behaviour was not observed by Douglas et al. (2014) who use Hα as their tracer of chromospheric activity. Similarly, again using Hα, Newton et al. (2017) find the rotation-activity relation to be valid for mid- and later-M dwarfs. These stars, however, are field stars and likely not young. Unfortunately, our rotation periods do not extend to lower masses and thus do not allow us to study the activity-rotation relation for mid-M dwarfs or to probe whether the unstructured rotation-activity diagram from Jackson & Jeffries (2010) is intrinsic to the chromospheric emission in the IRT.
6.3. Activity-informed period search
As already pointed out by Noyes et al. (1984), the relation between chromospheric activity indicators and Rossby number in principle allows one to estimate the stellar rotation period (e.g. Wright et al. 2004). We are in the advantageous position of possessing light curves of most of the stars in our activity sample. This allows us not only to estimate the rotation period from the activity, but also subsequently to determine the actual photometric rotation period from the photometry.
The estimated rotation period serves as a prior to our period search, and we calculated it for every star that falls either onto the slow rotator branch of the CAD (log < −4.0) or in the transition and gap region (log
< −3.8 and (V − Ks)0 < 3.6). We excluded the saturated stars (mostly M dwarfs) for which the activity level obviously cannot provide a useful prediction for the rotation period. We calculated the estimated rotation period Pest = Ro τc from Eqs. (3) and (4), respectively, for the two samples. The spread in Rossby number among the stars on the slow rotator sequence is ∼10%. The same relative uncertainty is assigned to Pest as ΔPest.
As described exhaustively in the companion paper, we calculated periodograms for each light curve with five different methods. However, for this particular phase of the work we use only four, discarding the string length method on account of its periodogram becoming noisy for low S/N light curves. Additionally, this method often does not agree with the other four methods, unnecessarily inflating our uncertainties. In the periodograms, we search for the highest peak within 2.5ΔPest of Pest. A manual selection process verifies the existence of a periodic signal in the light curve with this period. A visualisation of this process is given with the example light curves in Fig. A.2. We show the estimated period and the search range in the different periodograms. In keeping with the period determination algorithm in the companion paper, we calculate the final period as the mean of all four possible values, and its uncertainty as the spread among the four initial values. During this process, we restricted the search range for ∼15 light curves to assign the activity-guided period, and not a (chance) higher secondary or alias peak in the periodogram. We confirm all of our previously determined rotation periods, showing that we have included no alias periods among the slow rotators. Additionally, we find 116 new rotation periods among stars previously classified as non-rotators. By definition all of these rotation periods follow the rotation-activity relation. We provide the newly found periods in Table 1 along with the chromospheric activity data.
The success of this method can be attributed to multiple reasons. Firstly, the stellar activity measurement narrows the possible period range and enables us to choose the correct peak in noisy periodograms. It gives us confidence in assigning a period to a noisy light curve and enables us to find the true period in light curves with evolving spots. A second major reason is the breaking of aliases. Since we know that all of our stars are slow rotators (or in the rotational gap), the rotation period cannot be near 1 d for instance, and hence we break the alias with the observational cadence. Finally, we are also able to break the half and double period alias because a factor of two is a very large offset in the rotation-activity diagram. Jointly, these reasons increase the yield of rotation periods for our ground-based observations from 37% to 66% (255/387 stars16 with AAO spectra) in the correlated regime.
However, we are unable to assign a rotation period to every light curve in this set of slow rotators. Among the 184 considered stars, we were not able to assign a rotation period to nine. In some cases the light curve is just too noisy, resulting in periodograms without distinctive peaks. Occasionally, although peaks are found, the different methods provide inconsistent results. In both cases, we do not assign a rotation period.
6.4. Unsaturated M dwarfs on the extended slow rotator sequence
We now turn to an unexpected feature in the CAD – the unsaturated early M dwarfs. Among the stars in that region only one rotator was found in the original analysis of the light curves. We define this area as (V − Ks)0 > 3.7 and log < −4.05. Using the period estimation method described above, we find nine stars with periodic light curve variations. We examine these stars particularly closely because some of these rotation periods are longer than what we would normally have assigned with a 42 d long time series. (Generally, we have required that two complete phases be visible over the observational baseline). However, the stellar activity provides additional information that can justify these long periods. In Table 2, we list all such stars with activity informed rotation periods that can be found in the CAD on the declining branch among the M dwarfs.
Rotation periods for unsaturated M dwarfs and cluster membership criteria.
Despite the good agreement of the new slowly rotating stars in the general activity and the rotational properties of stars in NGC 3532, we reinvestigated their membership status for additional assurance. We applied stronger quality cuts than in our previous membership work (F19) in order to include only high-probability cluster members.
After reconsidering all membership criteria (Table 2), we reclassified three stars as likely non-members because their distances locate them farther from the cluster centre than the others. (We still list them in Table 2 so that they can be reinstated if or when more definitive membership information becomes available). The remaining six stars are cluster members with very high probability. We note, for completeness, that all considered stars are photometric cluster members. However, the non-member 612827 is on the photometric binary sequence and the member 304320 is a potential astrometric binary (RUWE = 1.64).
We show the light curves for the six likely members, together with the generalized Lomb-Scargle (Zechmeister & Kürster 2009), as well as the CLEAN (Roberts et al. 1987; Crane 2001) periodograms in Fig. A.3. The chosen periods do not always coincide with the highest peak in the periodograms for two reasons. Firstly, for the long rotation periods the peaks chosen are near or beyond the typical limit we set for our periods in the companion paper, where we require two cycles in our time series. Secondly, the rotation periods are assigned in concert with the chromospheric activity, which restricts the possible periods. In particular the two peaks that are seen in most periodograms in Fig. A.3 are typically aliases of each other. The chromospheric activity helps to break the half and double period alias effectively.
6.5. New periods in colour-activity-rotation space
With the large number (113) of new rotation periods, we are able to update the central figure of this work (Fig. 5), and show the complete picture, including the slowly rotating M dwarfs, in Fig. 7. We note that the Gaia parallaxes for 13 stars place them far from the cluster centre, despite all other membership criteria being fulfilled. We have conservatively excluded these stars from the rotational analysis despite their being potential cluster members, so as to avoid any contamination of our sample.
![]() |
Fig. 7. Same as Fig. 5 but now including the newly identified (activity-informed) rotation periods; single stars are marked with crosses and possible binaries with plus symbols. Each new period is also colour-coded to match the previously determined rotational groups, enabling cross-identification as in Fig. 5. This diagram includes the unsaturated extended slow rotator sequence (brown), absent in Fig. 5. The light curves and periodograms of the four stars marked with red squares (both in the CAD and CPD) are displayed in Fig. A.2 of the appendix to demonstrate the viability of deriving these activity-informed periods. |
The newly found periods can be observed to strengthen all key regions of the CAD with the exception of the saturated stars. Most notably, we are now able to assign periods to the least active, highest-mass stars of our sample, thus finding many additional periods along the whole slow rotator sequence, to fill it out substantially. In the corresponding rotation-activity diagram (Fig. 7c) all new stars can be seen to follow the same relation as the previously identified rotators. Since the fit to this diagram was the constraint on the new periods, this result is not unexpected, and serves mainly as a check on the results.
The CPD is where we find the most noticeable changes. Not only is the slow rotator sequence much more densely populated than before, but F-type stars are also now present, and the relationship now also extends into the M dwarf regime out to periods of ∼30 d. Slowly rotating M dwarfs in this period regime were also discovered in the Pleiades (Rebull et al. 2016) and NGC 2516 (F20). The stars in NGC 3532 have certainly spun down, and we analyse their rotational evolution in the companion paper (F21rot). The fact that these stars are both located on the rotation-activity relation and extend the slow rotator sequence seamlessly likely indicates that they share a dynamo-regime with the earlier-type stars despite being still on the pre-main sequence. A similar conclusion for main sequence stars was also drawn by Wright & Drake (2016), working with a sample of older, slowly rotating field M dwarfs. An in-depth analysis of the whole CPD is provided in the companion rotational study.
7. Conclusions
We have measured the chromospheric activity of cool stars in the 300 Myr old, rich open cluster NGC 3532 using the Ca II infrared triplet lines and have used these to explore and better define the corresponding colour-activity-rotation relationships. The cluster occupies an important transitional location between the well-studied ZAMS clusters (e.g. Pleiades, NGC 2516) and the 600 Myr old Hyades cluster, and its richness makes it especially valuable. This study follows on the heels of our prior spectroscopic membership study (F19) and a detailed study of the rotation period distribution in the companion paper (F21rot).
We showed how the chromospheric activity is well correlated with the photospheric activity, the latter measured from the starspot-induced variability amplitude of light curves. The correlation holds despite the fact that the two datasets are not contemporaneous, an encouraging sign from the observational perspective, and is ultimately rooted in their shared stellar rotation and its mass and age dependences.
We constructed one of the richest colour-activity diagrams (CAD) available to date for any open cluster. The diagram is remarkably rich in information, with a prominent (main) sequence of lower-activity stars and a side branch of higher-activity stars, the two separated by an ‘activity gap’. The CAD appears to be a near-mirror image of the colour-period diagram (CPD) of the cluster. Colour-coding specific groups of rotators to show the correspondence emphasises that such a correspondence occurs on a star-by-star, rather than statistical, basis. This correspondence suggests that the two quantities (rotation and activity) are potentially predictable from each other to a greater extent than previously thought possible.
In particular, the close correspondence between the activity and the rotation period shows that the absence of higher-mass fast rotators in our photometric work is an astrophysical fact and not a detection issue. We conclude that all stars down to a certain mass seem to have transitioned to the slow rotator sequence, showing the universality and strong mass dependence of this transition. The CAD hosts an unexpected group of low-activity M dwarfs, for which the stellar activity declines with declining mass. In contrast to the majority of young M dwarfs, these stars are clearly in the unsaturated regime.
We also constructed the first rotation-activity diagram for stars in NGC 3532 and, as expected, find the two well-known regimes (saturated and unsaturated) connected by a very wide transition region that can be traced back to the rotational gap in the CPD. The relationship is remarkably tight in relation to prior work, and there is only one outlier.
From this rotation-activity correlation, we predicted rotation periods for stars in the unsaturated regime, searched for, and identified the photometric rotation periods using the photometric light curves obtained for the companion rotational study. We recovered all previously found periods and were able to find 113 additional membership-justified rotation periods with the help of this activity-informed period search. Using the activity as a guideline allows us to break aliases and to assign photometric rotation periods to noisier light curves. The new periods include stars on the extended slow rotator sequence.
In conclusion, the analysis of chromospheric activity in combination with photometric time series analysis enabled us to perform a detailed study of the activity-rotation period correspondence over a wide range of stellar masses among cool stars in NGC 3532, showing that the relationship is tight enough to be predictive over a wide parameter range. In particular, it allows us to substantially increase the number of photometric rotation periods derived. In combination with prior work, this study has contributed towards the establishment of NGC 3532 as a benchmark open cluster and has wider implications for the era of massive multi-object spectroscopic surveys.
The model spectra in Fig. 2 are displayed in full resolution to give the reader a sense of the spectral features; for all calculations, the spectra have been resampled to the resolution of the observed spectra.
The more accurate empirical flux scale by Lorenzo-Oliveira et al. (2016) is optimised for solar-mass stars and is unfortunately not established for M dwarfs.
This is also seen in the good agreement between Hall (1996) and Lorenzo-Oliveira et al. (2016).
The largest set we found in the literature contains 516 members of Praesepe for which Douglas et al. (2014) measured the chromospheric activity in Hα.
Readers might consider comparing this with the CAD in Fig. 5 of James & Jeffries (1997), which shows a similar, but noisier, correlation.
A suitable prior comparison for this panel is Fig. 19 of Soderblom et al. (1993).
Excluding three potential non-members among the M dwarfs (see Sect. 6.4).
Acknowledgments
We thank Marcel Agüeros for the detailed review. SAB acknowledges support from the Deutsche Forschungs Gemeinschaft (DFG) through project number STR645/7-1. SPJ is supported by the German Leibniz-Gemeinschaft, project number P67-2018. Based on data acquired at the Anglo-Australian Telescope under program S/2017A/02. We acknowledge the traditional owners of the land on which the AAT stands, the Gamilaraay people, and pay our respects to elders past and present. Based in part on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory (2008A-0476; S. A. Barnes, SMARTS consortium through Vanderbilt University), which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. 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. This research has made use of NASA’s Astrophysics Data System Bibliographic Services. This research has made use of the SIMBAD database and the VizieR catalogue access tool, operated at CDS, Strasbourg, France. Software: This research made use of ASTROPY, a community-developed core Python package for Astronomy (Astropy Collaboration 2013). This research made use of the following PYTHON packages: PANDAS (McKinney 2010); NUMPY (van der Walt et al. 2011); MATPLOTLIB (Hunter 2007); IPYTHON: (Pérez & Granger 2007).
References
- AAO Software Team 2015, Astrophysics Source Code Library [record ascl:1505.015] [Google Scholar]
- Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58 [Google Scholar]
- Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269 [Google Scholar]
- Barnes, S. A., & Kim, Y.-C. 2010, ApJ, 721, 675 [Google Scholar]
- Basri, G., & Nguyen, H. T. 2018, ApJ, 863, 190 [NASA ADS] [CrossRef] [Google Scholar]
- Basri, G., Walkowicz, L. M., Batalha, N., et al. 2011, AJ, 141, 20 [NASA ADS] [CrossRef] [Google Scholar]
- Belokurov, V., Penoyre, Z., Oh, S., et al. 2020, MNRAS, 496, 1922 [Google Scholar]
- Boggs, P. T., & Rogers, J. E. 1990, Contemp. Math., 112, 186 [Google Scholar]
- Busà, I., Aznar Cuadrado, R., Terranegra, L., Andretta, V., & Gomez, M. T. 2007, A&A, 466, 1089 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cenarro, A. J., Cardiel, N., Gorgas, J., et al. 2001, MNRAS, 326, 959 [NASA ADS] [CrossRef] [Google Scholar]
- Clem, J. L., Landolt, A. U., Hoard, D. W., & Wachter, S. 2011, AJ, 141, 115 [Google Scholar]
- Crane, P. C. 2001, Sol. Phys., 203, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Douglas, S. T., Agüeros, M. A., Covey, K. R., et al. 2014, ApJ, 795, 161 [NASA ADS] [CrossRef] [Google Scholar]
- Douglas, S. T., Curtis, J. L., Agüeros, M. A., et al. 2019, ApJ, 879, 100 [Google Scholar]
- Eberhard, G., & Schwarzschild, K. 1913, ApJ, 38, 292 [Google Scholar]
- Fang, X.-S., Zhao, G., Zhao, J.-K., & Bharat Kumar, Y. 2018, MNRAS, 476, 908 [NASA ADS] [CrossRef] [Google Scholar]
- Foing, B. H., Crivellari, L., Vladilo, G., Rebolo, R., & Beckman, J. E. 1989, A&AS, 80, 189 [Google Scholar]
- Fritzewski, D. J., Barnes, S. A., James, D. J., et al. 2019, A&A, 622, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fritzewski, D. J., Barnes, S. A., James, D. J., & Strassmeier, K. G. 2020, A&A, 641, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fritzewski, D. J., Barnes, S. A., James, D. J., & Strassmeier, K. G. 2021, A&A, 652, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Güdel, M. 2004, A&ARv, 12, 71 [CrossRef] [Google Scholar]
- Hall, J. C. 1996, PASP, 108, 313 [NASA ADS] [CrossRef] [Google Scholar]
- Hall, J. C. 2008, Liv. Rev. Sol. Phys., 5, 2 [Google Scholar]
- Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
- Husser, T. O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jackson, R. J., & Jeffries, R. D. 2010, MNRAS, 407, 465 [NASA ADS] [CrossRef] [Google Scholar]
- James, D. J., & Jeffries, R. D. 1997, MNRAS, 292, 252 [Google Scholar]
- Jones, J., White, R. J., Boyajian, T., et al. 2015, ApJ, 813, 58 [NASA ADS] [CrossRef] [Google Scholar]
- Kraft, R. P. 1967, ApJ, 150, 551 [Google Scholar]
- Lewis, I. J., Cannon, R. D., Taylor, K., et al. 2002, MNRAS, 333, 279 [NASA ADS] [CrossRef] [Google Scholar]
- Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2 [EDP Sciences] [Google Scholar]
- Linsky, J. L., Hunten, D. M., Sowell, R., Glackin, D. L., & Kelch, W. L. 1979, ApJS, 41, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Lorenzo-Oliveira, D., Porto de Mello, G. F., Dutra-Ferreira, L., & Ribas, I. 2016, A&A, 595, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mamajek, E. E., & Hillenbrand, L. A. 2008, ApJ, 687, 1264 [Google Scholar]
- Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64 [Google Scholar]
- Marsden, S. C., Carter, B. D., & Donati, J. F. 2009, MNRAS, 399, 888 [NASA ADS] [CrossRef] [Google Scholar]
- Martin, J., Fuhrmeister, B., Mittag, M., et al. 2017, A&A, 605, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- McKinney, W. 2010, in Proceedings of the 9th Python in Science Conference, eds. S. van der Walt, & J. Millman, 51 [Google Scholar]
- Meibom, S., & Mathieu, R. D. 2005, ApJ, 620, 970 [Google Scholar]
- Meibom, S., Mathieu, R. D., & Stassun, K. G. 2007, ApJ, 665, L155 [NASA ADS] [CrossRef] [Google Scholar]
- Morris, B. M., Curtis, J. L., Douglas, S. T., et al. 2018, AJ, 156, 203 [NASA ADS] [CrossRef] [Google Scholar]
- Newton, E. R., Irwin, J., Charbonneau, D., et al. 2016, ApJ, 821, 93 [Google Scholar]
- Newton, E. R., Irwin, J., Charbonneau, D., et al. 2017, ApJ, 834, 85 [Google Scholar]
- Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763 [Google Scholar]
- Pettersen, B. R. 1989, Sol. Phys., 121, 299 [NASA ADS] [CrossRef] [Google Scholar]
- Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., & Ventura, P. 2003, A&A, 397, 147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pérez, F., & Granger, B. E. 2007, Comput. Sci. Eng., 9, 21 [Google Scholar]
- Radick, R. R., Lockwood, G. W., Henry, G. W., Hall, J. C., & Pevtsov, A. A. 2018, ApJ, 855, 75 [Google Scholar]
- Rauscher, E., & Marcy, G. W. 2006, PASP, 118, 617 [NASA ADS] [CrossRef] [Google Scholar]
- Rebull, L. M., Stauffer, J. R., Bouvier, J., et al. 2016, AJ, 152, 113 [Google Scholar]
- Roberts, D. H., Lehar, J., & Dreher, J. W. 1987, AJ, 93, 968 [Google Scholar]
- Schmitt, J. H. M. M. 1997, A&A, 318, 215 [NASA ADS] [Google Scholar]
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
- Skumanich, A. 1972, ApJ, 171, 565 [Google Scholar]
- Soderblom, D. R., & Mayor, M. 1993, AJ, 105, 226 [NASA ADS] [CrossRef] [Google Scholar]
- Soderblom, D. R., Stauffer, J. R., Hudon, J. D., & Jones, B. F. 1993, ApJS, 85, 315 [NASA ADS] [CrossRef] [Google Scholar]
- Strassmeier, K. G. 2009, A&ARv, 17, 251 [Google Scholar]
- van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22 [Google Scholar]
- Žerjal, M., Zwitter, T., Matijevič, G., et al. 2013, ApJ, 776, 127 [CrossRef] [Google Scholar]
- Wilson, O. C. 1978, ApJ, 226, 379 [Google Scholar]
- Wright, N. J., & Drake, J. J. 2016, Nature, 535, 526 [Google Scholar]
- Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261 [Google Scholar]
- Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48 [Google Scholar]
- Yamashita, M., Itoh, Y., & Takagi, Y. 2020, PASJ, 72, 80 [NASA ADS] [CrossRef] [Google Scholar]
- Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]
- Zhang, J., Bi, S., Li, Y., et al. 2020, ApJS, 247, 9 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Light curves of selected stars
![]() |
Fig. A.1. Light curve in time (top left) and phase space (top right) as well as the Lomb-Scargle periodogram (bottom left) and CLEAN periodogram (bottom right) for the outlier (707560) in the rotation-activity diagram. |
![]() |
Fig. A.2. Four examples of light curves and their associated periodograms for (non-M) stars with newly identified, activity-informed rotation periods. Each column shows the data for one star. The two upper panels show the light curves in both the time and phase domains. We fit a sine with the rotation period and display it along with the data. The left scale is valid for the two leftmost plots, while the right scale is for the two light curves on the right. The lower panels show, from top to bottom, the Lomb-Scargle, CLEAN, phase-dispersion minimization, and Gregory-Loredo periodograms. In each periodogram, we mark both the rotation period (orange line) and the period estimated from the activity and its uncertainty (blue lines). Unlike the other periodograms, in the phase-dispersion minimization the best periods are represented by dips. |
![]() |
Fig. A.3. Left column: Light curves for the early M dwarfs of low stellar activity in NGC 3532 that have newly identified rotation periods. The orange line is a sine fit to the data with the newly identified likely rotation period. The first two light curves (from the top) are from V 120 s exposures, while all others are from Ic 600 s images. Centre column: Corresponding normalised generalized Lomb-Scargle diagram, with the selected rotation period marked by the orange vertical line. The dashed blue line indicates the rotation period estimated from the chromospheric activity measurement, and the dotted lines delineate the search range, which is derived from the period uncertainty. Right column: Same as the centre column but with the CLEAN periodogram. The offset between the peaks and the position of the best period for 102025 in the first row is due to the different peak positions in the phase-dispersion minimization and Gregory-Loredo periodograms (not displayed). |
All Tables
Chromospheric activity measurements, including uncertainties and auxiliary data, for members of NGC 3532.
All Figures
![]() |
Fig. 1. Colour-magnitude diagram of the spectroscopically observed members of NGC 3532. In orange we show the rotators from Fritzewski et al. (2021) and in light blue from this work. We were unable to determine rotation periods for stars marked in green. Single stars are indicated with circles and candidate binaries from photometry and astrometry with squares. |
In the text |
![]() |
Fig. 2. Spectral cutouts for example stars of different chromospheric activity (orange) and their corresponding PHOENIX spectra. Each panel is centred on one IRT line and shifted to the rest frame. The line centre is marked with a dotted line and the limits of the 1 Å wide integration window with dashed lines. The spectra are grouped into low-activity (top) and high-activity stars (bottom). The latter feature prominent fill-ins of the line cores. Earlier spectral types are higher. The topmost spectrum is from a G2 star, and the lowest is from an M 3 star. All six stars are marked with squares in the colour-activity diagram in Fig. 5a. |
In the text |
![]() |
Fig. 3. Left: V band variability amplitude (VP95) for stars in NGC 3532 (orange), compared with those of the zero-age main sequence open cluster NGC 2516 (black) for both stars with periods (squares) and those without (circles). Stars in NGC 3532 generally have lower variability, as seen by the mean amplitudes expressed by the solid lines. Centre: VP90 against Rossby number Ro for the rotators in NGC 3532 and NGC 2516. The distribution can be seen to resemble the activity-rotation relation previously known from chromospheric and coronal activity, and it shows that the photospheric activity amplitude has a similar structure. Right: histogram of the variability amplitudes for the members of the two open clusters, showing that the high-amplitude stars of NGC 3532 are far less numerous in comparison with those of NGC 2516. The dotted grey histogram is restricted to NGC 2516 stars with (V − Ks)0 < 4 and confirms that the high-amplitude peak in NGC 2516 is mostly attributable to the fast-rotating M dwarf members, a population absent from our NGC 3532 data. |
In the text |
![]() |
Fig. 4. Photometric variability amplitude in the V band VP95 against chromospheric emission ratio log |
In the text |
![]() |
Fig. 5. Interconnection between chromospheric activity, rotation, and mass. The colour-activity diagram (a) connects the colour-period diagram (b) with the rotation-activity diagram (c). a: chromospheric emission ratio log |
In the text |
![]() |
Fig. 6. Rotation-activity diagram with two possible log-linear fits to two regions of interest: all unsaturated stars (blue, Ro > 0.06) and only the slow rotators (orange, Ro > 0.11). |
In the text |
![]() |
Fig. 7. Same as Fig. 5 but now including the newly identified (activity-informed) rotation periods; single stars are marked with crosses and possible binaries with plus symbols. Each new period is also colour-coded to match the previously determined rotational groups, enabling cross-identification as in Fig. 5. This diagram includes the unsaturated extended slow rotator sequence (brown), absent in Fig. 5. The light curves and periodograms of the four stars marked with red squares (both in the CAD and CPD) are displayed in Fig. A.2 of the appendix to demonstrate the viability of deriving these activity-informed periods. |
In the text |
![]() |
Fig. A.1. Light curve in time (top left) and phase space (top right) as well as the Lomb-Scargle periodogram (bottom left) and CLEAN periodogram (bottom right) for the outlier (707560) in the rotation-activity diagram. |
In the text |
![]() |
Fig. A.2. Four examples of light curves and their associated periodograms for (non-M) stars with newly identified, activity-informed rotation periods. Each column shows the data for one star. The two upper panels show the light curves in both the time and phase domains. We fit a sine with the rotation period and display it along with the data. The left scale is valid for the two leftmost plots, while the right scale is for the two light curves on the right. The lower panels show, from top to bottom, the Lomb-Scargle, CLEAN, phase-dispersion minimization, and Gregory-Loredo periodograms. In each periodogram, we mark both the rotation period (orange line) and the period estimated from the activity and its uncertainty (blue lines). Unlike the other periodograms, in the phase-dispersion minimization the best periods are represented by dips. |
In the text |
![]() |
Fig. A.3. Left column: Light curves for the early M dwarfs of low stellar activity in NGC 3532 that have newly identified rotation periods. The orange line is a sine fit to the data with the newly identified likely rotation period. The first two light curves (from the top) are from V 120 s exposures, while all others are from Ic 600 s images. Centre column: Corresponding normalised generalized Lomb-Scargle diagram, with the selected rotation period marked by the orange vertical line. The dashed blue line indicates the rotation period estimated from the chromospheric activity measurement, and the dotted lines delineate the search range, which is derived from the period uncertainty. Right column: Same as the centre column but with the CLEAN periodogram. The offset between the peaks and the position of the best period for 102025 in the first row is due to the different peak positions in the phase-dispersion minimization and Gregory-Loredo periodograms (not displayed). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.