Free Access
Issue
A&A
Volume 654, October 2021
Article Number A134
Number of page(s) 17
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202141122
Published online 22 October 2021

© ESO 2021

1. Introduction

The observational knowledge of fundamental properties of the stars is the basis for most of the astrophysics models explaining the evolution of stars and galaxies. In particular, the mass of the star is the main property that regulates its evolution and eventual fate. The determination of stellar masses with high precision is then essential to constrain these models. Double-lined eclipsing binaries (EBs) provide a unique opportunity to obtain direct measurements of the mass and radius of its stellar components with accuracies of ∼1% (Torres et al. 2010). However, the number of systems where this precision has been achieved is relatively low (Torres et al. 2010; Southworth 2015). This holds true, especially, for young (< 800 Myr) low-mass stars and brown dwarfs, as only a few are known and have accurate mass and radius determinations (David et al. 2019; Lodieu et al. 2020, and references therein). These systems are necessary to constrain pre-main-sequence (PMS) evolutionary models. Most of them have been recently discovered, made possible by the advent of ground-based surveys such as SuperWASP (Pollacco et al. 2006), and space missions such as CoRoT (Baglin et al. 2006), Kepler (Koch et al. 2004; Borucki et al. 2007), or TESS (Ricker et al. 2014), which monitor almost the whole night sky.

The detection of EBs is heavily constrained by the inclination of the orbit of the system (i.e., only nearly edge-on configurations produce eclipses). On the other hand, orbital configurations of double-lined spectroscopic binaries (SB2) are not as strongly limited. Thus, the number of known young SB2s (e.g., Dudorov & Eretnova 2016; Shkolnik et al. 2017; Zúñiga-Fernández et al. 2021, and references therein) is higher than EBs. Most of these young SB2s were discovered, also serendipitously, as part of spectroscopic surveys focusing on relatively few members of young stellar associations, the majority of them solar-type (FGK) stars, while the number of M dwarfs still remains low. Moreover, their orbits, for the most part, have not been characterised to date as only a few spectra for each of them are available. However, SB2s do not provide any information about the radii of the system components and only a minimum limit for their masses, even though they are more abundant than EBs.

The determination of the age of the system is a step needed to test the validity of the predictions of PMS models. This is usually achieved comparing the photometry of the object against isochrones derived from evolutionary models (e.g., Malo et al. 2014b; Herczeg & Hillenbrand 2015; Bell et al. 2015). In open clusters and young moving groups (YMGs), we can assume that all the members were born roughly at the same time (Zuckerman & Song 2004). Under this assumption we can expect that the age derived from the isochrones should be consistent for all of them. This provides a powerful tool for determining the age of a particular open cluster or YMG, and subsequently of any of its newly identified members. For this reason, a proper characterisation of the evolutionary models is also important for determining stellar ages. Nonetheless, the age of the members of open clusters and YMGs can also be estimated using other methods that are mostly independent or do not rely so heavily on evolutionary models, for example traceback simulations to find the convergence point; lithium depletion boundary (LDB) determinations; or age relations with youth indicators such as rotational period, X-ray, or Hα emission (see Soderblom 2010, and references therein)

The relative proximity to the Sun of most YMG members means they are bright enough to achieve a radial velocity (RV) precision of km s−1 with small telescopes (∼1 m) needed to measure dynamical masses. As part of an observing campaign with the Mercator 1.2 m telescope, we observed several early- to mid-type M dwarfs previously classified as members of YMGs in the literature. We confirmed one of these targets, GJ 1284, as an SB2 previously identified by Gizis et al. (2002). Later on, we also characterised the orbit of this system using the Mercator spectra and data from other facilities obtained during this and subsequent campaigns, and data gathered from the literature. The system was initially classified as a member of either the Beta Pictoris Moving Group (βPMG; Barrado y Navascués et al. 1999; Zuckerman et al. 2001) or the Columba Association (COL; Torres et al. 2008) by Malo et al. (2013). Moreover, one other system, GJ 898, had been identified kinematically as a wide co-moving companion to GJ 1284 (Shaya & Olling 2011). This wide companion forms part of a triple system alongside another binary star, GJ 897AB (Kuiper 1943). The three objects combined form a wide quintuple system. In this work we revise the membership of GJ 1284 to a YMG based on its new kinematics and youth indicators, which suggests that it belongs to the Local Association (LA; Eggen 1975; Montes et al. 2001), the same YMG as GJ 898 and GJ 897 (Montes et al. 2001; Makarov et al. 2008).

The LA, also commonly known as the Pleiades moving group, is one of the first classical YMGs identified by Eggen (1975) who initially described it as a stream associated with the Pleiades open cluster. Later on, Asiain et al. (1999) and López-Santiago et al. (2006) revealed that LA is actually made up of several subgroups with ages ranging from 10–300 Myr.

In this work we characterise the orbit of GJ 1284 using spectra collected with four different high-resolution spectrographs, and compare the results against the predictions made by several stellar evolution models. In Sect. 2 we introduce the GJ 1284 system and review the literature about this object and its wide companions. In Sect. 3 we describe the high-resolution spectra collected over several observing campaigns complemented by additional photometric and spectroscopic datasets from the literature. In Sect. 4 we measure the RVs from the spectra and determine the orbital parameters of the system. We also infer the rotational period of GJ 1284 system and its wide companions. Finally, in Sect. 5 we propose an age range for the three systems and compare their photometry and the mass of the GJ 1284 components against PMS evolutionary models and similar studies.

2. GJ 1284

2.1. Description of the system

GJ 1284 is located at a distance d = 15.906 ± 0.019 pc (Gaia early Data Release 3, eDR3; Gaia Collaboration 2016, 2021). It has been extensively analysed as a member of the solar neighbourhood due to its proximity, and has been part of several direct-image searches for substellar companions (e.g., Galicher et al. 2016; Meshkat et al. 2017; Naud et al. 2017; Nielsen et al. 2019). The system was classified as an M3.0V star by Hawley et al. (1996). More recently, Torres et al. (2006) classified it as an M2.0Ve with an uncertainty of one subclass using high-resolution spectra. In this work we estimate the spectral types of the binary components based on the colours derived from the photometric magnitudes.

Gizis et al. (2002) measured short-term variations in the RVs of GJ 1284 of over 46 km s−1 and identified it as an SB2, even though they did not characterise its orbit. Later it was again identified as an SB2 by Torres et al. (2006), who obtained three more RV measurements with a mean value of −5.7 km s−1. Jeffers et al. (2018) also measured the RV of GJ 1284 on seven consecutive nights obtaining a mean value of −4.4 km s−1 with a peak-to-peak variation of 62.63 km s−1. Additionally, Schneider et al. (2019) measured a single value of 8.18 ± 1.07 km s−1. Finally, Sperauskas et al. (2019) obtained two more values of −0.4 ± 1.0 and −14.6 ± 1.0 km s−1. A summary of all the individual measurements is given in Table 1. Nonetheless, the orbital parameters of this binary system have not been derived due to the lack of continuous monitoring.

Table 1.

Compilation of RVs measured for GJ 1284 from the literature.

Torres et al. (2006) first identified GJ 1284 as a young object (< 600 Myr) based on several youth indicators. Malo et al. (2013) proposed it as a member of either COL (99.2%) or βPMG (99.9%), depending on whether or not they included the mean RV value from Torres et al. (2006) in the analysis. Following this work, GJ 1284 has been commonly catalogued as either a member of COL (Malo et al. 2014a; Bell et al. 2015; Naud et al. 2017) or βPMG (Elliott et al. 2016; Shkolnik et al. 2017). More recently, it has been rejected as a member of these two YMGs (Gagné et al. 2018; Schneider et al. 2019) and was classified as a field star (Janson et al. 2017) using the updated astrometry from Gaia. Nevertheless, the analysis of the membership of this object in a YMG is based on the average value of the available RV measurements, not on its systemic velocity.

Malo et al. (2014a) compiled the value of several youth indicators for members of YMGs, including GJ 1284. In particular, four youth indicators were considered: the projected rotational velocity (v sin i), the pseudo-equivalent width (pEW) of the Hα emission line, the X-ray flux, and the pEW of the lithium 6707.8 Å absorption line. They report a value of v sin i = 3.9 ± 1.2 km s−1, in agreement with Torres et al. (2006). However, the authors note that this value might not be reliable because of the binarity of the system. A value of pEW(Hα) = −3.3 Å was obtained by Riaz et al. (2006) from medium-resolution spectra taken at the Cerro Tololo Inter-American Observatory (CTIO) 1.5 m telescope. Again, this value could be affected by the binarity, although it also agrees with the pEW(Hα) = −3.23 ± 0.08  Å obtained by Jeffers et al. (2018) from high-resolution spectra. The log fX = −11.44 erg s−1 cm−2 also comes from Riaz et al. (2006) using ROSAT data (Voges et al. 1999) and the conversion relations described in Schmitt et al. (1995). Finally, Torres et al. (2006) report that the Li 6708 Å line was not detected in their spectra.

2.2. Wide companions

Two other stellar systems are identified in the literature as common proper motion companions to GJ 1284: GJ 898 (Shaya & Olling 2011) and GJ 897AB (Kuiper 1943). We describe these two systems below.

GJ 898 is a K6 (Torres et al. 2006) located at a distance of d = 14.549 ± 0.006 pc and classified as a member of the Local Association by Montes et al. (2001) according to its kinematics. The projected separation in the sky to GJ 1284 is 3.59 deg or ∼2 × 106 au (roughly 1 pc) at their mean distance. The physical separation of 1.66 pc together with a difference in their proper motions of over 10σ (Gaia eDR3; Table 2) may suggest this is a chance alignment. Nonetheless, the galactocentric UVW velocities for GJ 898 and GJ 1284 (Sect. 5.1) are in very good agreement (Table 2; Fig. 1). The difference in the astrometry can thus be explained as a projection effect consequence of their wide separation. In Sect. 5.2 we analyse whether the two systems are physically bound.

thumbnail Fig. 1.

Revised galactocentric UVW velocities of GJ 1284 (blue), GJ 898 (orange), and GJ 897AB (green) from this work compared to other YMG members from Montes et al. (2001). The ellipses represent the 2σ values of UVW for each YMG. Similarly, we overplot βPMG and COL as defined by Gagné et al. (2018).

Table 2.

Stellar parameters of the GJ1284 system and its two wide companions.

GJ 897AB was identified as a companion to GJ 898 by Kuiper (1943), and therefore is also a companion to GJ 1284. It is a resolved tight visual binary (ρ = 0.57″; Mason et al. 2002) formed by two components of similar magnitude (VA = 11.14 mag; VB = 11.17 mag; WDS) and mass (0.5 ± 0.1 MHeintz 1986a), which suggests a similar spectral type M2–M3 (Riaz et al. 2006; Torres et al. 2006). The orbit of the system is nearly edge-on (i = 89.30°) with a period of 28.2 years (Heintz 1986b). While Woolley et al. (1970) list a distance for GJ 897 of 14.5 ± 2.7 pc and Gliese & Jahreiss (1991) give 12.9 ± 9.9 pc, the separation of the components combined with the inclination of the orbit affects the quality of the astrometric solutions of more modern catalogues. For instance, HIPPARCOS lists a parallax of 5 ± 34 mas (van Leeuwen 2007), while Gaia does not provide a full astrometric solution in the last eDR3. In Table 2 we list the parallax and proper motions from the dynamical solution provided by Tokovinin (2018) for GJ 897AB, which are in closest agreement with the values for GJ 898. In this work we adopt for GJ 897AB the distance and proper motions of GJ 898 given their higher precision and the proximity of the two objects in the sky, which indicates they are most probably physical companions located at the same distance. Makarov et al. (2008) also noted that the space velocities of GJ 897AB and GJ 898 match those of the LA members from Montes et al. (2001), and supported this assessment with other activity indicators (see Sect. 5.1).

We performed a search for other common proper motion systems in a radius of up to 10 degrees using the Gaia eDR3 data, but it did not produce any additional candidates. All together, GJ 898, GJ 897AB, and GJ 1284 form a group of at least five stellar objects with very similar galactocentric velocities, with GJ 898 being the most massive component. A schematic of the configuration of this quintuple system is shown in Fig. 2. The existence of this common proper motion system has an implication on the determination of the age of GJ 1284. As previously mentioned, both GJ 898 and GJ 897AB have been associated with LA in the literature. In Sect. 5.1 we discuss the membership of the three systems in a YMG based on their kinematics and other youth indicators to infer an age for the quintuplet.

thumbnail Fig. 2.

Schematic view of the spatial configuration of the quintuple system made up of GJ 898, GJ 1284AB, and GJ 897AB. Sizes and distances are not to scale.

3. Observations

We gathered 27 spectra for GJ 1284 with four distinct high-resolution spectrographs (Table 3). Additionally, this object and the two wide companions described in the previous section were observed by TESS, although the individual components of GJ 1284 and GJ 897AB are not resolved in the images.

Table 3.

Instrumental configuration of the spectrographs used in this paper.

3.1. Spectroscopy

3.1.1. HERMES

We first observed and identified GJ 1284 as an SB2 with the High Efficiency and Resolution Mercator Echelle Spectrograph (HERMES; Raskin et al. 2011) as part of observing programme CAT-141-MULTIPLE-2/18B (PI: Cardona Guillén). We also obtained additional spectra with this same instrument through progammes 128-Mercator4/18B (PI: Montes) and ITP18_8 (PI: Tremblay). HERMES is a fibre-fed echelle spectrograph mounted on the Mercator 1.2 m semi-robotic telescope located in the Roque de los Muchachos Observatory (ORM) on La Palma, Spain. A total of 11 spectra were gathered over several runs spanning from 9 November to 22 December 2018 under clear conditions and a seeing of ∼0.8 arcsec. The spectra were taken using the high-resolution mode (HRF) with an exposure time of 1800 s. This mode provides a spectral resolution R ∼ 85 000 and a wavelength coverage from 377 to 900 nm through a 2.5 arcsec fibre equipped with an image slicer. The data were reduced by the automated pipeline and RV toolkit (HermesDRS1). This tool provides one-dimensional wavelength-calibrated spectra using the standard reduction procedure (bias subtraction, flat-field correction, order extraction, and wavelength calibration) from which we derive the RV measurements described in Sect. 4.1. Figure 3 shows the HERMES spectra for each epoch zoomed in around the Na Iλ 8183.25 Å line.

thumbnail Fig. 3.

HERMES spectra of GJ 1284 showing the Na Iλ 8183.25 Å line.

3.1.2. NRES

We obtained five additional spectra on three separate nights (21, 23, and 25 December 2018) using the Network of Robotic Echelle Spectrographs (NRES; Siverd et al. 2018) as part of programme KEY2017AB-002a (PI: Brown). NRES is a set of identical high-resolution spectrographs mounted on the 1 m telescopes of the Las Cumbres Observatory network (LCOGT; Brown et al. 2013). Each spectrograph is fibre-fed by two telescopes simultaneously, and has a resolution R ∼ 53 000 and a wavelength coverage from 380 to 860 nm. The spectra used in this work were obtained using the node located in the South African Astronomical Observatory (SAAO) in Sutherland, South Africa. Currently, three other NRES elements are available in the LCOGT nodes located in Chile, the United States, and Israel. The data were reduced using the automatic NRES pipeline2. The pipeline yields a one-dimensional spectrum calibrated in wavelength similarly to HERMES. In Fig. 4 we plot the NRES spectra for each epoch showing the Na Iλ 8183.25 Å line.

thumbnail Fig. 4.

NRES spectra of GJ 1284 centred around the Na Iλ 8183.25 Å line.

3.1.3. HARPS-N

We also observed the target using the High Accuracy Radial velocity Planet Searcher North (HARPS-N; Cosentino et al. 2012) on four consecutive nights between 9 and 12 January 2019 under programme CAT-141-MULTIPLE-2/18B (PI: Cardona Guillén). HARPS-N is an echelle spectrograph mounted on the 3.58 m Telescopio Nazionale Galileo (TNG) also located in the ORM. The design of the instrument is similar to its predecessor HARPS (Mayor et al. 2003), covering a wavelength range from 383 to 690 nm with a resolution of R ∼ 115 000. We obtained two spectra with exposure times of 900 s on each of the observing nights, except on 12 January 2019 when we obtained only one spectrum. These spectra were reduced automatically at the end of each night using the standard HARPS-N data reduction pipeline (Cosentino et al. 2014; Smareglia et al. 2014). Figure 5 shows the HARPS-N spectra around the Ca II-H λ 3968.47 Å line.

thumbnail Fig. 5.

HARPS-N spectra of GJ 1284 focused on the Ca II-H λ 3968.47 Å line.

3.1.4. CAFE

We also used spectra collected by Jeffers et al. (2018) with the Calar Alto Fiber-fed Echelle spectrograph (CAFE; Aceituno et al. 2013). CAFE is a fibre-fed echelle spectrograph mounted on the Calar Alto 2.2 m telescope in Centro Astronómico Hispano Alemán (CAHA) in Almería, Spain. The spectrograph has a wavelength coverage ranging from 396 to 950 nm spanning 84 orders with a resolution R ∼ 62 000. A total of seven spectra were taken over the same number of consecutive nights from 5 August to 11 August 2014 with exposure times of 1800 s. The data were reduced using a version of the REDUCE package (Piskunov & Valenti 2002). More details about the observations and the data reduction process can be found in Jeffers et al. (2018). Figure 6 shows the CAFE spectra at each epoch around the Na Iλ 8183.25 Å line.

thumbnail Fig. 6.

CAFE spectra of GJ 1284 depicting the Na Iλ 8183.25 Å line.

3.2. Photometry

The Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2014) is a space mission launched on 18 April 2018. During its two-year prime mission, which started on 25 July 2018, TESS monitored more than 200 000 main-sequence stars, covering almost all of the sky using its four cameras with a combined field of view of 24° ×96° in continuous campaigns of 27 days. A pre-selection of targets were monitored with a two-minute cadence.

GJ 1284 and its two co-moving companions GJ 898 and GJ 897AB were observed by TESS as part of Sector 2, from 22 August to 20 September 2018, and Sector 29, from 26 August to 22 September 2020. For GJ 1284 and GJ 897AB there is a two-minute cadence light curve available for each of these sectors, while for GJ 898 only the light curve for Sector 29 can be retrieved. The target input catalogue (TIC) numbers for the three objects are shown in Table 2; their light curves are publicly available and can be downloaded from the Mikulski Archive for Space Telescopes (MAST) Portal3. The two-minute-cadence photometry includes two different sets of data: the single aperture photometry (SAP) light curve, and the pre-search data conditioning (PDC) light curve. The PDC light curve attempts to correct the common instrumental effects present in the SAP light curve. In this work we use the SAP light curve to avoid corrections that may affect the activity-induced photometric variations, especially in relatively slow rotators with low-amplitude variability like ours. The two-minute cadence light curves of the three systems are shown in Fig. 7.

thumbnail Fig. 7.

TESS light curves and periodograms of GJ 1284 (top), GJ 898 (middle) and GJ 897 (bottom). Left: TESS photometric time series. Middle: GLS periodogram of the photometric time series. Right: light curve phase-folded with the selected rotation periods in Table 2. We adjust the ordinate axis range to crop the stellar flares to better visualise the variations in the light curve due to the stellar rotation. The colours in all the panels correspond to data from Sector 2 (blue) or Sector 29 (red) or both combined (black).

4. Data analysis

4.1. Radial velocities determination

The usual approach to calculating the RV of an object with the cross-correlation function (CCF) technique can lead to the determination of less precise values in SBs because of the blending of the spectral lines of the binary components. For this reason we used TODMOR (Zucker et al. 2004) to derive the RV values from all the spectra gathered for this object. TODMOR is a multi-order implementation of the TwO-Dimensional CORrelation technique (TODCOR; Zucker & Mazeh 1994). It uses a combination of two template spectra to simultaneously calculate the CCF for both components of the binary system. This technique is applied to each echelle order, combining the resulting CCFs into a single one to derive the RV.

We used a set of template spectra created using PHOENIX stellar models (Husser et al. 2013) as a reference to calculate the CCF. These templates have varying values of effective temperature (Teff) that can be broadened to match the projected rotational velocity (v sin i) of the star. We fixed the value for the metallicity [Fe/H] = 0 given the proximity of the system to the Sun. We used a range for Teff from 3000 to 4000 K in steps of 100 K and for v sin i from 0 to 10 km s−1 in steps of 1 km s−1. An additional parameter α was included to account for the flux-ratio between the two template spectra to calculate the CCF. We calculated the values that provide the best fit of the CCF using the HERMES spectrum observed closer to quadrature, when the lines of the binary components are most separated. We found that the optimal parameter values for the stellar templates are Teff = 3800 ± 100 K and v sin i = 6 ± 1 km s−1 for the primary component and Teff = 3600 ± 100 K and v sin i = 8 ± 1 km s−1 for the secondary, with a flux ratio α = 0.47 ± 0.10. These templates are used to derive the RV for the two components simultaneously at each epoch (see Zucker & Mazeh 1994, for more details). Taking advantage of the CCF per spectral order provided by TODMOR, we also inspected and excluded from the analysis those orders that are heavily contaminated by telluric lines, and those with low signal-to-noise ratios (S/N), mainly in the blue part of the spectrum where the peaks of the two components are not apparent in the CCF. The final RVs derived by TODMOR are given in Table 4.

Table 4.

Summary of the RV measurements and the associated uncertainties obtained from the spectra along the BJD and barycentric corrections at each epoch.

In addition, we calculated and applied the barycentric correction (BVCOR) to each RV epoch. HARPS-N applies this correction to the spectra themselves so that the extracted RV values are already in the barycentric reference frame. Moreover, it also provides the Barycentric Julian Day (BJD). For HERMES the value for the BVCOR and the BJD of each exposure is calculated as part of the pipeline and provided in the header of each image. Instead, CAFE and NRES do not provide these values as part of the automatic pipeline. For these instruments, we calculated the BJD and BVCOR at each epoch with the routines included in the astropy Python package4. This package uses an implementation of the expressions in Wright & Eastman (2014).

To improve the determination of the RVs we also accounted for relative drifts in the wavelength calibration of the different instruments. To palliate this issue we tracked the position of the telluric lines at each epoch using the IRAF routine fxcor (Tody 1986) to calculate a one-dimensional CCF. We correlated the HERMES spectrum with the highest S/N, which provides the widest wavelength coverage of all the spectrographs, against the rest of the spectra and chose the regions where the presence of the telluric lines is more apparent (e.g., 8119–8179, 8203–8324, 8928–8997 Å). The drift of the telluric lines relative to the reference spectrum is calculated from the highest peak of the resulting CCF. For the HERMES spectra we find that the internal drift is comparable (< 0.5 km s−1) to the uncertainty associated with the individual RV calculated using TODMOR. In the case of the CAFE spectra the internal drift is of the same order as HERMES, but there is a drift of ∼1 km s−1 between the two instruments. For NRES the internal drift appears to be about ∼0.5 km s−1, similar to HERMES and also comparable to the systematic drift between the two instruments. Lastly, the internal drift in the HARPS-N spectra is lower than the precision we can achieve with this method of 0.1 km s−1 (i.e., there is no noticeable internal drift). Meanwhile, there seems to be a systematic drift of ∼ − 0.3 km s−1 with respect to the reference HERMES spectrum, lower than the internal drift of this spectrograph. For the final analysis, we corrected the RV values derived with TODMOR from the drifts measured using the telluric lines at each epoch (Table 4). The internal drifts show that the stability of the four spectrographs is good enough to characterise our binary system, which presents variations of over 60 km s−1.

Table 4 includes the final RVs used in the orbital parameter determination (Sect. 4.2) and their uncertainties together with the BJD, barycentric, and telluric correction at each epoch. In Fig. 8 we plot the RVs as a function of time (top panel) and orbital phase (bottom panel).

thumbnail Fig. 8.

RV values for the four different spectrographs (CAFE, HERMES, HARPS-N, and NRES) and orbital solution for the primary (blue) and secondary (orange) component of GJ 1284 as a function of time (top panel) and orbital phase (bottom panel). The residuals of the Keplerian fit calculated with rvfit are shown at the bottom of each panel using the same scale.

4.2. Orbital parameters determination

We use the IDL code rvfit to determine the orbital parameters of the binary system (Iglesias-Marzoa 2015). This routine uses an adaptive simulated annealing global minimisation method to fit the two components of the system simultaneously with a Keplerian orbit defined by seven parameters: the orbital period (P), the time periastron (Tp), the eccentricity (e), the argument of periastron (ω), the RV amplitude of the two components (K1 and K2), and the systemic radial velocity (γ). From these we can derive lower bounds for the masses and semi-major axes of the binary components as well as the mass ratio q = (M1/M2). We estimate an initial value for the orbital period from the phase-folded RV values. The excellent phase coverage of the orbit and the wide time span relative to the orbital period between the CAFE observations and the rest of the exposures allowed us to obtain a very precise initial value, close to the final value provided later on by rvfit. The priors and final parameters of the best solution (Fig. 8) to the simultaneous fit of the two components are shown in Table 5.

Table 5.

Parameters of the simultaneous Keplerian fit to both components of GJ 1284 from rvfit.

From this analysis, we determine a value for the systemic RV of GJ 1284 γ = 0.84 ± 0.14 km s−1. This value differs from the one used in the literature to determine the membership of this object to a YMG (Sect. 2). Nonetheless, we note that the individual RV measurements derived by various authors all fall within the range of expected RV values for the primary component at different orbital phases.

We also compare the orbital period and eccentricity of this system against other spectroscopic binaries from Pourbaix et al. (2004) in Fig. 9. We note that the orbit of GJ 1284 is highly eccentric for its orbital period, lying close to the upper envelope (solid black line) defined by the binaries of the sample.

thumbnail Fig. 9.

Eccentricity and rotation period of GJ 1284 (in blue) against binary systems from Pourbaix et al. (2004). The black solid line represents the upper envelope for this distribution.

4.3. Rotational period determination

We determined the rotational period of the three systems from the TESS photometric time series (Sect. 3.2). First, we calculated the generalised Lomb-Scargle (GLS) periodogram (Zechmeister & Kürster 2009) of the photometry and identify the most significant peaks (see central panels of Fig. 7). Then, we visually inspected the photometric data folded with the period corresponding to the highest peak. In Fig. 7 we can see that GJ 898 presents a well-defined modulation with a two-peak pattern that is responsible for the peaks at 10.99 and 6.19 days in the periodogram. However, we find visually that the dispersion in the phase-folded light curve is minimised for a period of 11.5 ± 0.1 days, well within the width of the 10.99 days peak, adopting this value as the rotational period of GJ 898. On the contrary, GJ 1284 and GJ 897AB present a more complex pattern that can be attributed to the overlapping of the modulation of each component of these two binaries. Both stars also show intense flare activity, cropped out in Fig. 7 for clarity.

For GJ 897AB, we find that the shape of the light curve is quite different in the two TESS sectors. In Sector 2 (in blue in Fig. 7) we find an irregular modulation with a two-peak pattern showing two main peaks at 5.32 and 2.43 days in the periodogram. In contrast, the light curve for Sector 29 (in red) shows a clear modulation with a period of 5.02 days. The periodogram of the combined light curves of the two sectors (in black) shows a peak at 5.00 days. This suggests that the periodicity at ∼5 days in the light curve is due either to the individual components or to both of them, assuming they have similar rotational period given their similar spectral types and masses. Thus, we adopted a value of 5.00 ± 0.02 d as the rotational period of this system.

The light curve of GJ 1284 also shows a different behaviour for each sector, similarly to GJ 897. The light curve of Sector 2 (in blue) presents a more complex structure, while in Sector 29 it is more regular. The periodograms of both sectors have a main peak at around ∼7.5 days. We expect a shorter rotation period for the secondary component than for the primary, given its later spectral type (Sect. 5.3.1), which we do not detect in the TESS light curve as it probably appears diluted due to its lower brightness (Sect. 5.3.2). Thus, we determine a rotational period for GJ 1284A of 7.5 ± 0.1 days.

We followed the same procedure described in this section to calculate the rotation period of the LA members from Montes et al. (2001). First, we cross-matched all 120 members with the Gaia eDR3 with a matching radius of 2″ accounting for their proper motions to determine the new positions at the Gaia equinox (J2015.5). We rejected those matches that do not provide a full astrometric solution and radial velocity in Gaia, ending up with a total of 88 stars, for which we derived the UVW velocities and checked that their kinematics are still compatible with a membership to the LA. Of these 88 stars, 86 have two-minute-cadence TESS light curves available, but we can only recover a reliable measure of the rotation period for 64 of them, which are listed in Table A.1 and plotted in yellow in Fig. 10c.

thumbnail Fig. 10.

Comparison of the values of four different youth indicators for the three systems analysed in this work (GJ 1284, GJ 898, GJ 897) against young objects from the literature. Panel a: Pseudo-equivalent width of the Hα line against members of the Pleiades, Praesepe, and Hyades open clusters from Fang et al. (2018). Panel b: X-ray luminosity of the individual components compared with members from four different open clusters from Núñez & Agüeros (2016) presented in a box-plot. Panel c: rotational period derived from the TESS light curves alongside a compilation of members of the Pleiades from Rebull et al. (2016) and of the Hyades from Douglas et al. (2016), and rotational periods for LA members in Montes et al. (2001) derived in this work (see Sect. 5.1). Panel d: upper limits of the EW of the Li 6708 Å feature plotted alongside the upper envelopes of the IC 2602 and Hyades open clusters, as well as the upper and lower envelopes of the Pleiades (Gutiérrez Albarrán et al. 2020, and references therein)

5. Discussion

5.1. Youth assessment

GJ 1284 has been proposed as a member of βPMG (∼25 Myr; Binks & Jeffries 2014; Malo et al. 2014b) and COL (∼40 Myr; Torres et al. 2008; Bell et al. 2015) (Sect. 2). Furthermore, its two wide companions have been kinematically assigned to LA (10–300 Myr; López-Santiago et al. 2009) (Sect. 2.2). The membership of GJ 1284 to one of these groups, and thus its age, is fundamental so it can be compared with stellar evolutionary models.

We calculated the UVW galactocentric velocities from the revised systemic RV γ (Sect. 4.2) and the astrometry from the Gaia eDR3 using the transformations described in Johnson & Soderblom (1987). We also revised the kinematics of GJ 898 using the astrometry and RV from Gaia eDR3. Moreover, we calculated the 3D space motions of GJ 897AB with the astrometry from Tokovinin (2018) and the RV from Torres et al. (2006) (Sect. 2.2). We show the revised UVW velocities in Fig. 1 and list the values in Table 2. The kinematics of all three systems agree and are compatible with other members of LA from Montes et al. (2001). Moreover, our determination of the systemic velocity of GJ 1284 place this object outside βPMG and COL as defined in the BANYAN-Σ code (Gagné et al. 2018).

We employed several youth indicators to further constrain the kinematic age, including the Hα pseudo-equivalent width pEW(Hα), the near-ultraviolet (NUV) excess emission measured with the colour NUV − J, the X-ray luminosity LX, the stellar rotational period, and the pEW of the Li feature at 6707.8 nm.

In the upper left panel of Fig. 10 we show that all three objects present Hα emission compatible with members of the Pleiades (∼110–130 Myr; Stauffer et al. 1998; Dahm 2015), but also with some stars from the Hyades (∼650–800 Myr; Martín et al. 2018; Douglas et al. 2019; Lodieu 2020, and references therein) according to Fang et al. (2018). We set the age range of the systems between the estimated age of these two open clusters.

From the 2MASS (Skrutskie et al. 2006) J and K magnitudes and the GALEX (Martin et al. 2005) NUV magnitude (Table 2), we inferred colours J − K = 0.871, 0.763, 0.854 mag and NUV − J = 11.689, 11.427, 11.395 mag for GJ 1284, GJ 898, and GJ 897, respectively. Findeisen et al. (2011) identified a correlation between age and NUV excess emission for members of YMGs. The NUV − J colours of our objects indicate an excess emission in the NUV compatible with members of YMGs younger than the Hyades (see Figs. 7, 8, and 9 in Findeisen et al. 2011).

We derived a log(LX) = 28.15,  28.74,  28.92 erg s−1 for GJ 898 and the individual components of GJ 1284AB and GJ 897AB, respectively, assuming comparable X-ray luminosities for the two components of the binary systems. These levels of X-ray emission are compatible with objects younger than the Hyades (Núñez & Agüeros 2016), as shown in the upper right panel of Fig. 10.

We compared the rotational periods of GJ 1284A, GJ 898, and GJ 897A against members of the Pleiades and Praesepe open clusters from Rebull et al. (2016) and Douglas et al. (2016, 2019), respectively (Fig. 10). We calculated the G − J colour of GJ 897A assuming that the two components have similar spectral types. As for GJ 1284, the method for calculating the colour of its binary components is described in detail in Sect. 5.3.1.

The lower left panel of Fig. 10 shows that the rotational periods of the three stars are longer that the values expected for members of the Pleiades, but shorter than for members of Praesepe (∼600–750 Myr; Douglas et al. 2019, and references therein). From the rotational period, we conclude that the age of these three systems is most likely bracketed by the age of the two open clusters. We note that the rotational period of GJ 1284 might be affected by the tidal interaction between the two components. Simulations show that the spin and orbit of binary systems with orbital periods of less than 10 days might synchronise on a timescale as short as 200 Myr (Fleming et al. 2019). However, this is not the case for GJ 1284AB as its orbital period is roughly 1.5 times longer than the rotational period of the primary component. In any case, the conclusions of the age determination derived from the rotational period of GJ 898 and GJ 897A still hold.

We found no presence of the Li6708 line in our GJ 1284 spectra similarly to Torres et al. (2006) setting an upper limit of pEWLi6708 < 35 mÅ from our spectra. We also analysed an archival HARPS spectra for GJ 898 (programme 085.C-0019(A); PI: Lo Curto) again finding no sign of the Li6708 line, determining an upper limit for its pEWLi6708 < 17 mÅ. This indicates that their age must be older than IC 2602 (40–50 Myr Dobbie et al. 2010), and that it is compatible with members of the Pleiades (Gutiérrez Albarrán et al. 2020) or older, as shown in the lower right panel of Fig. 10.

We can also put constraints on the age of the GJ 1284 based on the high eccentricity of its orbit (Table 5). Meibom & Mathieu (2005) determined a relation of the circularisation period with age from binary members of several open clusters. Although the Hyades cluster appears as an outlier in the original work, a revision of the binaries used in the analysis showed that the cut-off period for this open cluster is in fact 8.5–11.9 days (see Gillen et al. 2017, for more details). The Hyades upper limit is similar to the orbital period of our system (Table 5), thus we infer that it must be younger than the age of this open cluster, which agrees with the expected age derived from the other youth indicators. Moreover, we note again the high eccentricity of the orbit, which can be seen in Fig. 9 where GJ 1284AB is located near the upper envelope, indicating that it is still far from being circularised and suggesting that it is well below the age of the Hyades.

We summarise the age range estimated from the lower limit of the younger age and the upper limit from the older age for each activity indicator in Table 6. The upper limit is set by the estimated maximum age of the stellar associations part of LA. In Fig. 10c we see that the rotation periods of the LA members from Table A.1 (in yellow) lie mostly between the sequences defined by the Praesepe and Pleiades stars, which is consistent with the estimated age range of 10–300 Myr from the literature. However, we also note that some of them follow the Praesepe sequence or are located slightly above it. This leads us to believe that the upper limit for the age of this association is not well constrained, although it is indeed young. Given that the main conclusions from this work are drawn from the fact that the objects have arrived to the main sequence (see Sect. 5.3.1), which is firmly based on the lower limit for the age imposed by the rotational period, the upper limit for the age of LA does not affect them. Here we adopt a maximum age for the LA similar to those of Praesepe and Hyades. The estimated age ranges for the different indicators are then compatible with an age for the system of 110–800 Myr.

Table 6.

Summary of the possible age ranges for GJ 1284 according to the different youth indicators considered in this work.

5.2. Gravitational bond of the quintuple system

Throughout Sect. 5.1 and in the following Sect. 5.3.1 we work under the assumption that GJ 1284 and its two wide companions, GJ 898 and GJ 897AB, have a common origin based on their kinematics and thus are roughly the same age. Even though the youth indicators and kinematics support the membership of these three systems to the LA, given the wide age range of this YMG we cannot confirm they have the same age. Moreover, the wide physical separation of 1.66 pc between them might indicate that they are not physically bound.

To explore the different scenarios, we first determine the probability of finding two members of the LA with physical separations and velocities as similar as GJ 1284 and GJ 898 based on the galactocentric positions and velocities of the members from Montes et al. (2001). We estimate that the probability that two of them have UVW velocities as similar as our systems is ∼10–15%, depending on whether we select all the LA members or only the denser region in the UVW velocity space. Furthermore, the probability that two LA members are closer than 1.66 pc is lower than 0.2%, given that members of YMGs are not grouped as in open clusters, but instead are dispersed over a much larger volume. Thus, the combined probability of finding two systems at these separations in the position and velocity space is only 0.05% or less.

We also calculate the disruption times of our system to test the hypothesis that they are physically bound using the expressions from Weinberg et al. (1987) and Binney & Tremaine (1987). For this calculation we adopt a mean mass for the eventual perturber of 1 M, a system density of 0.078–0.084 pc−3 (Henry et al. 2018) and a mean velocity dispersion for the perturbers of 59.8 km s−1 (Zapatero Osorio et al. 2007) in the solar neighbourhood, together with the masses for GJ 898 (∼0.66 M) and GJ 1284AB (∼0.77 M, see Sect. 5.3.2) and their separation (1.66 pc). With these parameters we derive a disruption time tdisr ≃ 100 − 200 Myr, which is comparable to the age estimated from the youth indicators (Sect. 5.1) suggesting that the system might still be physically bound. However, the difference in the W velocities of GJ 1284 and GJ 898, an order of magnitude higher than the escape velocity of the latter of ∼100 m s−1, indicates otherwise. Nonetheless, if we trace back the motion of these two systems using their UVW velocities, we find that they reached their minimum separation (0.6 pc) just 1.3 ± 0.5 Myr ago. Using this distance instead, we infer a disruption time tdisr ≃ 400 − 700 Myr, longer than the time we calculated previously, but still comparable to the age of the system. In summary, although currently the system does not seem to be physically bound, we believe that it is very unlikely to find such relatively close systems in the LA by chance alignment. Additionally, we cannot rule out that the system was physically bound until 1–2 Myr ago when it could have been disrupted by an impact from a perturber.

5.3. Model comparison

5.3.1. Colour-magnitude diagrams

In this section we compare the absolute magnitudes and colours of the three systems with four different PMS evolutionary models: Lyon models (BHAC15; Baraffe et al. 2015); MESA Isochrones & Stellar Tracks (MIST; Dotter 2016; Choi et al. 2016); PAdova and TRieste Stellar Evolution Code (PARSEC; Bressan et al. 2012; Chen et al. 2014); and Stellar Parameters of Tracks with Starspots (SPOTS; Somers et al. 2020). A brief description of the characteristics of each of these models can be found in David et al. (2019). In Fig. 11 we plot the absolute V magnitudes and the V − J colour of each system, and of their individual components together with the isochrones from each theoretical model.

thumbnail Fig. 11.

Comparison of the measured magnitudes for GJ 1284, GJ 898, and GJ 897 and the calculated magnitude for the individual components of the two binary systems against a set of isochrones from four different PMS evolutionary models: BHAC15 (top left), PARSEC (top right), MIST (bottom left), and SPOTS (bottom right). We also include the isomasses (dashed lines) for each of the models for initial mass values ranging from 0.7 to 0.2 M.

To separate the integrated photometry of GJ 1284 into the corresponding magnitudes for the individual components we use the integrated GaiaGBP magnitude and the flux ratio derived by TODMOR from the HARPS-N spectra. The wavelength range of the HARPS-N spectra closely resembles the transmission curve of the GBP filter5, hence justifying this approximation. The absolute GBP is calculated from the parallax of this object (Table 2). Using these values, we derive absolute magnitudes in the GBP filter of 10.73 ± 0.09 and 11.55 ± 0.13 mag for the primary and secondary components, respectively. Additionally, we use the relations for old field stars from Pecaut & Mamajek (2013)6 assuming that the system is close to the main sequence based on its age (Sect. 5.1). We base this approximation on the fact that, for solar-type stars and early- to mid-type M dwarfs, the isochrones of the evolutionary models for ages older than ∼120 Myr, as our objects appear to be, squeeze together (Fig. 11). Moreover, a comparison with empirical masses and luminosities of members of the NGC 1647 (∼150 Myr; Dias et al. 2002), Pleiades, and Praesepe open clusters with similar spectral types (Sect. 5.3.2) show that the mass-luminosity sequences overlap at this range of ages. According to the GBP absolute magnitude, we infer a spectral type of M2–M2.5 and M3–M3.5 respectively for the primary and secondary components of GJ 1284 using these relations. We estimate the difference in magnitudes between the two components for the different filters from these same relations and decompose the other integrated photometry into the individual magnitudes. We take into account the uncertainties in the spectral type estimates to calculate the photometric error bars. For GJ 897AB, the two components of the binary system have similar spectral types (Sect. 2.2), hence we assume identical fluxes.

Figure 11 shows that the BHAC15, MIST, and SPOTS models predict a younger isochronal age for GJ 1284A/B and GJ 897A/B (∼30–50 Myr) than for GJ 898 (∼100 Myr). This effect was already noticed by Herczeg & Hillenbrand (2015) using the BHAC15 models, which predicted younger ages for very low-mass stars than for solar-type stars in other YMGs. In Table 7 we also list the isochronal age calculated minimising the combined distance of the three systems to each isochrone assuming they are coeval. To account for the uncertainties in the magnitude and colour of the individual components we performed a Monte Carlo simulation with 5 000 realisations to estimate the error in the age determination. We see that these ages are younger than those derived from the youth indicators in Sect. 5.1, in particular from the rotational period. The determination of isochronal ages younger than expected for members of other young associations was also noted in other works (e.g., David et al. 2016; Gillen et al. 2017). On the other hand, PARSEC models suggest a common age of ∼40–140 Myr for all of them, consistent with the youth indicators. The version v1.2S of the PARSEC models considered in this work includes an ad hoc correction to reconcile the models with the observational data, particularly for very low-mass stars (Chen et al. 2014). This might explain the differences between PARSEC and the other three sets of models.

Table 7.

Summary of the isochronal ages for each model derived for the three systems analysed in this work, assuming coevality.

The discrepancy in the age derived by the different models indicates that there might be an additional factor that varies for the three systems. For example, the SPOTS models include the effects of magnetic fields in the form of stellar spots, while the other models do not (Somers et al. 2020). Somers & Pinsonneault (2015) explored in detail the effect of spot-related activity in the stellar evolution and found that low-mass stars with starspots have lower luminosities and temperatures compared to those with no spots. As a consequence, ages in spotted stars can be underestimated by a factor of 2–10. However, the spot-filling factor cannot be directly determined from our observations. In Fig. 12 we show isochrones derived from the SPOTS models for different spot-filling factors. We can see for the typical maximum value of 50% in PMS low-mass stars (Somers & Pinsonneault 2015 and references therein) that the isochronal ages of the three systems are in better agreement with the lower limit for the age ∼100 Myr derived from the youth indicators and the values derived by the PARSEC models.

thumbnail Fig. 12.

Isochrones from the SPOTS model for three different ages (10, 50, and 200 Myr) and different spot-filling factors (0, 34, and 51%).

5.3.2. Mass-luminosity relations

In this section we compare the minimum masses of GJ 1284A and B (Sect. 4.2) and other dynamical mass determinations with the predictions from the PMS models in a mass-luminosity diagram (Fig. 13). We use the bolometric corrections from Pecaut & Mamajek (2013) for the spectral types of the individual components and the GBP magnitudes of GJ 1284A and B (Sect. 5.3.1) to derive a bolometric luminosity of log(Lbol/L) = − 1.61 ± 0.05 and −1.77 ± 0.09 for the primary and the secondary components, respectively. These values for the luminosities and the minimum masses are compatible with the theoretical mass-luminosity relations.

thumbnail Fig. 13.

Minimum masses and luminosities of GJ 1284 A (blue) and B (orange) against the mass-luminosity relations derived from the BHAC15 (top left), PARSEC (top right), MIST (bottom left), and SPOTS (bottom right) models. The blue and orange horizontal bars represent the expected masses at different inclination angles for the primary and the secondary components, respectively. The dashed vertical lines indicate the masses derived from the mass-luminosity relations for the minimum and maximum ages derived from the youth indicators (shaded in grey). For comparison, we overplot the dynamical masses of the EB members of NGC 1647, the Pleiades (red) and Praesepe (green) open clusters from Lodieu et al. (2020).

Given the bolometric luminosity of the individual components of GJ 1284, we can estimate their masses using the mass-luminosity relations from the models assuming an age range for the system of 110–800 Myr (Sect. 5.1) obtaining a mean mass of 0.432 − 0.462 M and 0.365 − 0.390 M for the primary and secondary, respectively. The intervals of masses for each model (Table 8) take into account the estimated range of ages. We compare these results with dynamical masses measured for EBs with similar ages and spectral types (Fig. 13) from a compilation by Lodieu et al. (2020), namely with members of the Pleiades (in red) from David et al. (2016) and Hebb et al. (2006) and Praesepe (in green) from Kraus et al. (2017) and Gillen et al. (2017) with spectral types M2–M4 for at least one of their components. We see that the masses derived for GJ 1284A and B from their luminosity and age are consistent with the dynamical masses for these EBs, and in particular they seem to be delimited between the masses of 2MASSJ0446+19A (M2.5V; 0.47 ± 0.05 M; Hebb et al. 2006) and PTFEB132.707+19.810 (M3.5V; 0.3953 ± 0.0020 M; Kraus et al. 2017), whose spectral types are also compatible with the GJ 1284 components (Sect. 5.3.1). Assuming that the masses of each component derived from the photometry are correct, they imply an inclination angle for the orbit of the binary system of 43.6 − 42.8° or 43.5 − 42.3° according to the minimum mass of primary and the secondary, respectively.

Table 8.

Summary of maximum and minimum masses (in M) of the two binary components derived from the isochrones of the different models.

We can also see, with the exception of PARSEC, that the masses derived from the mass-luminosity relations for each model are higher than the values inferred from the CMD for the same model. Moreover, the masses derived from the CMD are lower than the dynamical masses from the literature. The colour and magnitude corrections needed to reconcile the masses derived from the mass-luminosity relations and the CMDs are listed in Table 9. This suggests that the mass-luminosity relations seem to be more reliable to derive masses from the photometry than CMDs. One of the possible explanations for this could be that, although the evolutionary models may correctly predict the values of the luminosities, they appear to underestimate the effective temperatures. This was already noted in other works (David et al. 2016; Gillen et al. 2017), where models predicted effective temperatures 200–300 K hotter for early to mid M dwarfs. Alternatively, model atmospheres and spectral synthesis for these objects may predict magnitudes and colours for these very cool objects that are not fully correct. Furthermore, this also indicates that, unless the age of the system is younger than the value derived from the youth indicators in this work, only the PARSEC models seem to predict masses consistent with those derived from mass-luminosity relations and CMDs and also compatible with empirical masses of EBs in open clusters older than 100 Myr.

Table 9.

Colour and magnitudes differences needed to reconcile the mass estimates derived from the CMD and the mass-luminosity relations for each model.

6. Conclusions

In this work we analyse a young SB2 system, GJ 1284, alongside its two wide companion, GJ 898 and GJ 897AB, and compare the results of this analysis with several PMS evolutionary models. The main results and conclusions from this work are as follows:

  • We characterised the orbit of GJ 1284, a young SB2 system, using RVs derived using high-resolution spectra from HERMES, HARPS-N, CAFE, and NRES. We find that this system is composed of two stars of similar mass (q = 0.845 ± 0.012) with estimated spectral types of M2–M2.5 for the primary component and M3–M3.5 for the secondary in an eccentric orbit (e = 0.505 ± 0.005) with a period of P = 11.838033 ± 0.000016 d. We also derived a minimum mass of 0.141 ± 0.003, 0.119 ± 0.003 M for the primary and secondary components, respectively.

  • We re-analysed the kinematics of GJ 1284AB using the revised systemic velocity of 0.84 ± 0.14 km s−1, which suggests a membership in the LA, consistent GJ 898 and GJ 897AB, in the same YMG. Based on several youth indicators, we adopted an age range for the three systems of 110–800 Myr, assuming they have a common origin.

  • We find that the probability of finding a system such as the one made up of GJ 1284 and its two relatively close co-moving companions in the LA by chance alignment is very low, suggesting they might have a common origin. Moreover, although currently they do not seem to be physically bound, we cannot rule out that they were in the past before being disrupted just 1–2 Myr ago.

  • We find that the isochronal ages derived from the comparison of the photometry with four PMS evolutionary models (BHAC15, PARSEC, MIST, and SPOTS) derived from the photometry of these three systems are younger than those estimated from the youth indicators, except for the PARSEC models. Moreover, the isochrones of the models, again excluding PARSEC, suggest a different age for the early M dwarf binaries than for the solar-type star.

  • Finally, we see that the masses for GJ 1284 A and B from the mass-luminosity relations for the different evolutionary models, except PARSEC, are higher than the values derived from the photometry alone, but are consistent with the dynamical masses of EBs with similar spectral type and age from the literature. This suggests that the determination of masses using mass-luminosity relations from evolutionary models may be more reliable than a comparison with the photometry in CMDs. Alternatively, the system might be younger than we estimated from the youth indicators.


6

See updated version including Sloan and Gaia filters in https://www.pas.rochester.edu/~emamajek/EEM_dwarf_UBVIJHK_colors_Teff.txt

Acknowledgments

We thank the anonymous referee for his/her comments that helped to improve this work. CCG is supported by the Ministerio de Economía y Competitividad del Gobierno de España (MINECO) under project SEV-2015-0548-16-3. We acknowledge financial support from the Agencia Estatal de Investigación of the Ministerio de Ciencia, Innovación y Universidades through projects PID2019-109522GB-C5[3, 4] and AEI/10.13039/501100011033. Based on observations made with the Mercator Telescope, operated on the island of La Palma by the Flemish Community, at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Based on observations obtained with the HERMES spectrograph, which is supported by the Research Foundation - Flanders (FWO), Belgium, the Research Council of KU Leuven, Belgium, the Fonds National de la Recherche Scientifique (F.R.S.-FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de Genève, Switzerland and the Thüringer Landessternwarte Tautenburg, Germany. 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 work makes use of observations from the NRES spectrograph on the LCOGT 1m telescope at the SAAO node of the Las Cumbres Observatory global telescope network. Based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundación Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. Based on observations collected at the Centro Astronómico Hispano-Alemán (CAHA) at Calar Alto, operated jointly by Junta de Andalucía and Consejo Superior de Investigaciones Científicas (IAA-CSIC). This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA’s Science Mission Directorate. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

References

  1. Aceituno, J., Sánchez, S. F., Grupp, F., et al. 2013, A&A, 552, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Asiain, R., Figueras, F., Torra, J., & Chen, B. 1999, A&A, 341, 427 [NASA ADS] [Google Scholar]
  3. Baglin, A., Auvergne, M., Boisnard, L., et al. 2006, in 36th COSPAR Sci. Assem., 3749 [Google Scholar]
  4. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Barrado y Navascués, D., Stauffer, J. R., Song, I., & Caillault, J. P. 1999, ApJ, 520, L123 [CrossRef] [Google Scholar]
  6. Bell, C. P. M., Mamajek, E. E., & Naylor, T. 2015, MNRAS, 454, 593 [Google Scholar]
  7. Binks, A. S., & Jeffries, R. D. 2014, MNRAS, 438, L11 [NASA ADS] [CrossRef] [Google Scholar]
  8. Binney, J., & Tremaine, S. 1987, Galactic dynamics (Princeton, N.J.: Princeton University Press), c1987 [Google Scholar]
  9. Borucki, W., Koch, D., Basri, G., et al. 2007, Proc. Int. Astron. Union, 3, 17 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [Google Scholar]
  11. Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031 [Google Scholar]
  12. Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525 [Google Scholar]
  13. Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 [Google Scholar]
  14. Cosentino, R., Lovis, C., Pepe, F., et al. 2012, in Ground-based and Airborne Instrumentation for Astronomy IV, eds. S. K. Ramsay, I. S. McLean, & H. Takami, 84461V [Google Scholar]
  15. Cosentino, R., Lovis, C., Pepe, F., et al. 2014, in Ground-based and Airborne Instrumentation for Astronomy V, eds. S. K. Ramsay, I. S. McLean, & H. Takami, 91478C [Google Scholar]
  16. Dahm, S. E. 2015, ApJ, 813, 108 [Google Scholar]
  17. David, T. J., Conroy, K. E., Hillenbrand, L. A., et al. 2016, AJ, 151, 112 [NASA ADS] [CrossRef] [Google Scholar]
  18. David, T. J., Hillenbrand, L. A., Gillen, E., et al. 2019, ApJ, 872, 161 [Google Scholar]
  19. Dias, W. S., Alessi, B. S., Moitinho, A., & Lépine, J. R. D. 2002, A&A, 389, 871 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Dobbie, P. D., Lodieu, N., & Sharp, R. G. 2010, MNRAS, 409, 1002 [Google Scholar]
  21. Dotter, A. 2016, ApJS, 222, 8 [Google Scholar]
  22. Douglas, S. T., Agüeros, M. A., Covey, K. R., et al. 2016, ApJ, 822, 47 [NASA ADS] [CrossRef] [Google Scholar]
  23. Douglas, S. T., Curtis, J. L., Agüeros, M. A., et al. 2019, ApJ, 879, 100 [Google Scholar]
  24. Dudorov, A. E., & Eretnova, O. V. 2016, Astron. Astrophys. Trans., 29, 437 [NASA ADS] [Google Scholar]
  25. Eggen, O. J. 1975, PASP, 87, 37 [Google Scholar]
  26. Elliott, P., Bayo, A., Melo, C. H. F., et al. 2016, A&A, 590, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Fang, X.-S., Zhao, G., Zhao, J.-K., & Bharat Kumar, Y. 2018, MNRAS, 476, 908 [NASA ADS] [CrossRef] [Google Scholar]
  28. Findeisen, K., Hillenbrand, L., & Soderblom, D. 2011, AJ, 142, 23 [Google Scholar]
  29. Fleming, D. P., Barnes, R., Davenport, J. R. A., & Luger, R. 2019, ApJ, 881, 88 [Google Scholar]
  30. Gagné, J., Mamajek, E. E., Malo, L., et al. 2018, ApJ, 856, 23 [Google Scholar]
  31. Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Galicher, R., Marois, C., Macintosh, B., et al. 2016, A&A, 594, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Gillen, E., Hillenbrand, L. A., David, T. J., et al. 2017, ApJ, 849, 11 [Google Scholar]
  35. Gizis, J. E., Reid, I. N., & Hawley, S. L. 2002, AJ, 123, 3356 [Google Scholar]
  36. Gliese, W., & Jahreiss, H. 1991, NASA STI/Recon Tech. Rep. A, 224, 161 [NASA ADS] [Google Scholar]
  37. Gutiérrez Albarrán, M. L., Montes, D., Gómez Garrido, M., et al. 2020, A&A, 643, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Hawley, S. L., Gizis, J. E., & Reid, I. N. 1996, AJ, 112, 2799 [Google Scholar]
  39. Hebb, L., Wyse, R. F. G., Gilmore, G., & Holtzman, J. 2006, AJ, 131, 555 [NASA ADS] [CrossRef] [Google Scholar]
  40. Heintz, W. D. 1986a, A&AS, 65, 411 [NASA ADS] [Google Scholar]
  41. Heintz, W. D. 1986b, AJ, 92, 446 [NASA ADS] [CrossRef] [Google Scholar]
  42. Henden, A. A., Templeton, M., Terrell, D., et al. 2016, VizieR Online Data Catalog:II/336 [Google Scholar]
  43. Henry, T. J., Jao, W.-C., Winters, J. G., et al. 2018, AJ, 155, 265 [NASA ADS] [CrossRef] [Google Scholar]
  44. Herczeg, G. J., & Hillenbrand, L. A. 2015, ApJ, 808, 23 [Google Scholar]
  45. Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Janson, M., Durkan, S., Hippler, S., et al. 2017, A&A, 599, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Jeffers, S. V., Schöfer, P., Lamert, A., et al. 2018, A&A, 614, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [Google Scholar]
  49. Jones, D. O., & West, A. A. 2016, ApJ, 817, 1 [NASA ADS] [CrossRef] [Google Scholar]
  50. Koch, D. G., Borucki, W., Dunham, E., et al. 2004, Overview and status of the Kepler Mission, 1491 [Google Scholar]
  51. Kraus, A. L., Douglas, S. T., Mann, A. W., et al. 2017, ApJ, 845, 72 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kuiper, V. G. P. 1943, ApJ, 97, 275 [NASA ADS] [CrossRef] [Google Scholar]
  53. Lodieu, N. 2020, Mem. della Soc. Astron. Ital., 91, 84 [Google Scholar]
  54. Lodieu, N., Paunzen, E., & Zejda, M. 2020, in Rev. Front. Mod. Astrophys. (Cham: Springer International Publishing), 213 [CrossRef] [Google Scholar]
  55. López-Santiago, J., Montes, D., Crespo-Chacon, I., & Fernandez-Figueroa, M. J. 2006, ApJ, 643, 1160 [CrossRef] [Google Scholar]
  56. López-Santiago, J., Micela, G., & Montes, D. 2009, A&A, 499, 129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. López-Santiago, J., Montes, D., Gálvez-Ortiz, M. C., et al. 2010, A&A, 514, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Makarov, V. V., Zacharias, N., & Hennessy, G. S. 2008, ApJ, 687, 566 [NASA ADS] [CrossRef] [Google Scholar]
  59. Malo, L., Doyon, R., Lafrenière, D., et al. 2013, ApJ, 762, 88 [Google Scholar]
  60. Malo, L., Artigau, É., Doyon, R., et al. 2014a, ApJ, 788, 81 [NASA ADS] [CrossRef] [Google Scholar]
  61. Malo, L., Doyon, R., Feiden, G. A., et al. 2014b, ApJ, 792, 37 [NASA ADS] [CrossRef] [Google Scholar]
  62. Martin, D. C., Fanson, J., Schiminovich, D., et al. 2005, ApJ, 619, L1 [Google Scholar]
  63. Martín, E. L., Lodieu, N., Pavlenko, Y., & Béjar, V. J. S. 2018, ApJ, 856, 40 [CrossRef] [Google Scholar]
  64. Mason, B. D., Hartkopf, W. I., Urban, S. E., et al. 2002, AJ, 124, 2254 [NASA ADS] [CrossRef] [Google Scholar]
  65. Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
  66. Meibom, S., & Mathieu, R. D. 2005, ApJ, 620, 970 [Google Scholar]
  67. Meshkat, T., Mawet, D., Bryan, M. L., et al. 2017, AJ, 154, 245 [Google Scholar]
  68. Montes, D., López-Santiago, J., Gálvez, M., et al. 2001, MNRAS, 328, 45 [NASA ADS] [CrossRef] [Google Scholar]
  69. Myers, J. R., Sande, C. B., Miller, A. C., Warren, W. H. J., & Tracewell, D. A. 2015, VizieR Online Data Catalog:V/145 [Google Scholar]
  70. Naud, M.-E., Artigau, É., Doyon, R., et al. 2017, AJ, 154, 129 [Google Scholar]
  71. Nielsen, E. L., De Rosa, R. J., Macintosh, B., et al. 2019, AJ, 158, 13 [Google Scholar]
  72. Núñez, A., & Agüeros, M. A. 2016, ApJ, 830, 44 [Google Scholar]
  73. Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 [Google Scholar]
  74. Piskunov, N. E., & Valenti, J. A. 2002, A&A, 385, 1095 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Pollacco, D. L., Skillen, I., Cameron, A. C., et al. 2006, PASP, 118, 1407 [Google Scholar]
  76. Pourbaix, D., Tokovinin, A. A., Batten, A. H., et al. 2004, A&A, 424, 727 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Raskin, G., Van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69 [CrossRef] [EDP Sciences] [Google Scholar]
  78. Rebull, L. M., Stauffer, J. R., Bouvier, J., et al. 2016, AJ, 152, 113 [Google Scholar]
  79. Riaz, B., Gizis, J. E., & Harvin, J. 2006, AJ, 132, 866 [Google Scholar]
  80. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, Proc. SPIE, 9143, 914320 [Google Scholar]
  81. Riedel, A. R., Blunt, S. C., Lambrides, E. L., et al. 2017, AJ, 153, 95 [NASA ADS] [CrossRef] [Google Scholar]
  82. Schmitt, J. H. M. M., Fleming, T. A., & Giampapa, M. S. 1995, ApJ, 450, 392 [Google Scholar]
  83. Schneider, A. C., Shkolnik, E. L., Allers, K. N., et al. 2019, AJ, 157, 234 [Google Scholar]
  84. Shaya, E. J., & Olling, R. P. 2011, ApJS, 192, 2 [CrossRef] [Google Scholar]
  85. Shkolnik, E. L., Allers, K. N., Kraus, A. L., Liu, M. C., & Flagg, L. 2017, AJ, 154, 69 [Google Scholar]
  86. Siverd, R. J., Brown, T. M., Barnes, S., et al. 2018, in Ground-based and Airborne Instrumentation for Astronomy VII, eds. H. Takami, C. J. Evans, & L. Simard, Proc. SPIE, 10702, 107026C [NASA ADS] [Google Scholar]
  87. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [Google Scholar]
  88. Smareglia, R., Bignamini, A., Knapic, C., Molinaro, M., & GAPS Collaboration 2014, ASP Conf. Ser., 485, 435 [NASA ADS] [Google Scholar]
  89. Soderblom, D. R. 2010, ARA&A, 48, 581 [Google Scholar]
  90. Somers, G., & Pinsonneault, M. H. 2015, ApJ, 807, 174 [Google Scholar]
  91. Somers, G., Cao, L., & Pinsonneault, M. H. 2020, ApJ, 891, 29 [Google Scholar]
  92. Southworth, J. 2015, DEBCat: A Catalog of Detached Eclipsing. Binary Stars, 496, 164 [NASA ADS] [Google Scholar]
  93. Sperauskas, J., Deveikis, V., & Tokovinin, A. 2019, A&A, 626, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Stauffer, J. R., Schultz, G., & Kirkpatrick, J. D. 1998, ApJ, 499, L199 [NASA ADS] [CrossRef] [Google Scholar]
  95. Tody, D. 1986, The Iraf Data Reduction And Analysis System, 733 [Google Scholar]
  96. Tokovinin, A. 2018, ApJS, 235, 6 [NASA ADS] [CrossRef] [Google Scholar]
  97. Torres, C. A. O., Quast, G. R., da Silva, L., et al. 2006, A&A, 460, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Torres, C. A. O., Quast, G. R., Melo, C. H. F., & Sterzik, M. F. 2008, in Handb. Star Form. Reg. Vol. II, ed. B. Reipurth, 5, 757 [NASA ADS] [Google Scholar]
  99. Torres, G., Andersen, J., & Giménez, A. 2010, A&ARv, 18, 67 [Google Scholar]
  100. van Leeuwen, F. 2007, A&A, 474, 653 [CrossRef] [EDP Sciences] [Google Scholar]
  101. Voges, W., Aschenbach, B., Boller, T., et al. 1999, ApJ, 349, 389 [Google Scholar]
  102. Weinberg, M. D., Shapiro, S. L., & Wasserman, I. 1987, ApJ, 312, 367 [NASA ADS] [CrossRef] [Google Scholar]
  103. Woolley, R., Epps, E. A., Penston, M. J., & Pocock, S. B. 1970, R. Obs. Ann., 5, 227 [Google Scholar]
  104. Wright, J. T., & Eastman, J. D. 2014, PASP, 126, 838 [Google Scholar]
  105. Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44 [Google Scholar]
  106. Zapatero Osorio, M. R., Martin, E. L., Bejar, V. J. S., et al. 2007, ApJ, 666, 1205 [NASA ADS] [CrossRef] [Google Scholar]
  107. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]
  108. Zucker, S., & Mazeh, T. 1994, ApJ, 420, 806 [NASA ADS] [CrossRef] [Google Scholar]
  109. Zucker, S., Mazeh, T., Santos, N. C., Udry, S., & Mayor, M. 2004, A&A, 426, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  110. Zuckerman, B., & Song, I. 2004, ARA&A, 42, 685 [Google Scholar]
  111. Zuckerman, B., Song, I., Bessell, M. S., & Webb, R. A. 2001, ApJ, 562, L87 [NASA ADS] [CrossRef] [Google Scholar]
  112. Zúñiga-Fernández, S., Bayo, A., Elliott, P., et al. 2021, A&A, 645, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Additional table

Table A.1.

Rotation period derived from the TESS light curves of the LA members compiled from Montes et al. (2001).

All Tables

Table 1.

Compilation of RVs measured for GJ 1284 from the literature.

Table 2.

Stellar parameters of the GJ1284 system and its two wide companions.

Table 3.

Instrumental configuration of the spectrographs used in this paper.

Table 4.

Summary of the RV measurements and the associated uncertainties obtained from the spectra along the BJD and barycentric corrections at each epoch.

Table 5.

Parameters of the simultaneous Keplerian fit to both components of GJ 1284 from rvfit.

Table 6.

Summary of the possible age ranges for GJ 1284 according to the different youth indicators considered in this work.

Table 7.

Summary of the isochronal ages for each model derived for the three systems analysed in this work, assuming coevality.

Table 8.

Summary of maximum and minimum masses (in M) of the two binary components derived from the isochrones of the different models.

Table 9.

Colour and magnitudes differences needed to reconcile the mass estimates derived from the CMD and the mass-luminosity relations for each model.

Table A.1.

Rotation period derived from the TESS light curves of the LA members compiled from Montes et al. (2001).

All Figures

thumbnail Fig. 1.

Revised galactocentric UVW velocities of GJ 1284 (blue), GJ 898 (orange), and GJ 897AB (green) from this work compared to other YMG members from Montes et al. (2001). The ellipses represent the 2σ values of UVW for each YMG. Similarly, we overplot βPMG and COL as defined by Gagné et al. (2018).

In the text
thumbnail Fig. 2.

Schematic view of the spatial configuration of the quintuple system made up of GJ 898, GJ 1284AB, and GJ 897AB. Sizes and distances are not to scale.

In the text
thumbnail Fig. 3.

HERMES spectra of GJ 1284 showing the Na Iλ 8183.25 Å line.

In the text
thumbnail Fig. 4.

NRES spectra of GJ 1284 centred around the Na Iλ 8183.25 Å line.

In the text
thumbnail Fig. 5.

HARPS-N spectra of GJ 1284 focused on the Ca II-H λ 3968.47 Å line.

In the text
thumbnail Fig. 6.

CAFE spectra of GJ 1284 depicting the Na Iλ 8183.25 Å line.

In the text
thumbnail Fig. 7.

TESS light curves and periodograms of GJ 1284 (top), GJ 898 (middle) and GJ 897 (bottom). Left: TESS photometric time series. Middle: GLS periodogram of the photometric time series. Right: light curve phase-folded with the selected rotation periods in Table 2. We adjust the ordinate axis range to crop the stellar flares to better visualise the variations in the light curve due to the stellar rotation. The colours in all the panels correspond to data from Sector 2 (blue) or Sector 29 (red) or both combined (black).

In the text
thumbnail Fig. 8.

RV values for the four different spectrographs (CAFE, HERMES, HARPS-N, and NRES) and orbital solution for the primary (blue) and secondary (orange) component of GJ 1284 as a function of time (top panel) and orbital phase (bottom panel). The residuals of the Keplerian fit calculated with rvfit are shown at the bottom of each panel using the same scale.

In the text
thumbnail Fig. 9.

Eccentricity and rotation period of GJ 1284 (in blue) against binary systems from Pourbaix et al. (2004). The black solid line represents the upper envelope for this distribution.

In the text
thumbnail Fig. 10.

Comparison of the values of four different youth indicators for the three systems analysed in this work (GJ 1284, GJ 898, GJ 897) against young objects from the literature. Panel a: Pseudo-equivalent width of the Hα line against members of the Pleiades, Praesepe, and Hyades open clusters from Fang et al. (2018). Panel b: X-ray luminosity of the individual components compared with members from four different open clusters from Núñez & Agüeros (2016) presented in a box-plot. Panel c: rotational period derived from the TESS light curves alongside a compilation of members of the Pleiades from Rebull et al. (2016) and of the Hyades from Douglas et al. (2016), and rotational periods for LA members in Montes et al. (2001) derived in this work (see Sect. 5.1). Panel d: upper limits of the EW of the Li 6708 Å feature plotted alongside the upper envelopes of the IC 2602 and Hyades open clusters, as well as the upper and lower envelopes of the Pleiades (Gutiérrez Albarrán et al. 2020, and references therein)

In the text
thumbnail Fig. 11.

Comparison of the measured magnitudes for GJ 1284, GJ 898, and GJ 897 and the calculated magnitude for the individual components of the two binary systems against a set of isochrones from four different PMS evolutionary models: BHAC15 (top left), PARSEC (top right), MIST (bottom left), and SPOTS (bottom right). We also include the isomasses (dashed lines) for each of the models for initial mass values ranging from 0.7 to 0.2 M.

In the text
thumbnail Fig. 12.

Isochrones from the SPOTS model for three different ages (10, 50, and 200 Myr) and different spot-filling factors (0, 34, and 51%).

In the text
thumbnail Fig. 13.

Minimum masses and luminosities of GJ 1284 A (blue) and B (orange) against the mass-luminosity relations derived from the BHAC15 (top left), PARSEC (top right), MIST (bottom left), and SPOTS (bottom right) models. The blue and orange horizontal bars represent the expected masses at different inclination angles for the primary and the secondary components, respectively. The dashed vertical lines indicate the masses derived from the mass-luminosity relations for the minimum and maximum ages derived from the youth indicators (shaded in grey). For comparison, we overplot the dynamical masses of the EB members of NGC 1647, the Pleiades (red) and Praesepe (green) open clusters from Lodieu et al. (2020).

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.