Brought to you by:

Gaia-based Isochronal, Kinematic, and Spatial Analysis of the epsilon Cha Association

, , and

Published 2021 January 28 © 2021. The American Astronomical Society. All rights reserved.
, , Citation D. Annie Dickson-Vandervelde et al 2021 AJ 161 87 DOI 10.3847/1538-3881/abd0fd

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

1538-3881/161/2/87

Abstract

The precise parallax, proper motion, and photometric measurements contained in Gaia Data Release 2 (DR2) offer the opportunity to reexamine the membership and ages of nearby young moving groups (NYMGs), i.e., loose groups of stars of age ≲100 Myr in the solar vicinity. Here, we analyze the available DR2 data for members and candidate members of the epsilon Cha Association (epsilonCA) which, at an estimated age of ∼3–5 Myr, has previously been identified as among the youngest NYMGs. The several dozen confirmed members of epsilonCA include MP Mus and T Cha, two of the nearest stars of roughly solar mass that are known to host primordial protoplanetary disks, and the Herbig Ae/Be star HD 104237A. We have used Gaia DR2 data to ascertain the Galactic positions and kinematics and color–magnitude diagram positions of epsilonCA members and candidates so as to reassess their membership status and thereby refine estimates of the distance, age, multiplicity, and disk fraction of the group. Our analysis yields 36 bona fide epsilonCA members, as well as 20 provisional members, including 3 new members identified here as comoving companions to previously known epsilonCA stars. We determine a mean distance to epsilonCA of 101.0 ± 4.6 pc and confirm that, at an age of ${5}_{-2}^{+3}$ Myr, epsilonCA represents the youngest stellar group within ∼100 pc of Earth. We identify several new photometric binary candidates, bringing the overall multiplicity fraction (MF) of epsilonCA to 40%, intermediate between the MFs of young T associations and the field.

Export citation and abstract BibTeX RIS

1. Introduction

Nearby young moving groups (NYMGs) serve as laboratories for the study of stars and their planetary offspring during their first stages of development, i.e., over the first tens of Myr after these young stars and planets have emerged from their birth clouds (Kastner et al. 2019). The age range spanned by coeval NYMGs (a few Myr to several ×100 Myr) provides unique opportunities to study stellar properties over a range of masses and temperatures at specific snapshots of pre-main sequence (pre-MS) and young MS evolution. Thanks to the proximity (<120 pc) of NYMGs, it is possible to observe and characterize an entire stellar population, down to the diminutive brown dwarfs (e.g., Schneider et al. 2019; Phillips et al. 2020) and even massive planets (e.g., Gagné et al. 2018b). The younger NYMGs (ages <25 Myr) feature Sun-like stars orbited by gas-rich circumstellar disks that are likely sites of ongoing planet formation (e.g., Sacco et al. 2014), while dusty debris disks are found among groups at more advanced ages (e.g., Zuckerman et al. 2019).

The precise astrometric and photometric data flowing from the Gaia Space Astrometry mission (Gaia Collaboration et al. 2016, 2018) is particularly useful to the study of NYMGs. Gaia Data Release 1 and the subsequent Gaia Data Release 2 (hereafter DR2) have resulted in a leap in our understanding of the membership of previously known NYMGs and, hence, their fundamental properties, such as ages, mass distributions, internal kinematics, and stellar multiplicity statistics (e.g., Gagné & Faherty 2018; Wright & Mamajek 2018; Lee & Song 2019). Gaia DR2 data have also led to the identification of previously unknown NYMGs (e.g., Gagné et al. 2018a).

Here, we present an analysis of Gaia DR2 data for the epsilon Chamaeleontis Association (hereafter epsilonCA). At an estimated age of 3–5 Myr (Murphy et al. 2013, henceforth M+13), epsilonCA is the youngest NYMG within ∼100 pc of Earth (Kastner et al. 2019 and references therein). Its relative youth and consequent large disk frequency (29%, M+13)—combined with its position well in the foreground of the Cha star formation complex (see Figure 2 of M+13)—make it a key NYMG for purposes of studying planet formation and pre-MS stellar evolution, without the complication of pervasive ambient or intervening molecular cloud material.

The identification and study of stars belonging to epsilonCA began in association with the study of the neighboring η Cha cluster (Mamajek et al. 2000). Chandra X-ray Observatory observations of the Herbig Ae/Be star HD 104237A epsilonCA led to further characterization of the group and discovery of new members (Feigelson et al. 2003). HD 104237A has a close, T Tauri-type companion (HD 104237B), both of which are orbited by a CO-rich circumbinary disk (Hales et al. 2014). An additional three stars are thought to be part of the HD 104237 system, two of which form another close binary (HD 104237D+E); HD 104237D shows evidence of accreting plasma (Testa et al. 2008).

In addition to HD 104237A, the members of epsilonCA include two stars with gas-rich circumstellar disks, MP Mus and T Cha; these stars represent two of only four known examples of solar-mass stars with protoplanetary disks that lie within ∼100 pc of Earth (Kastner et al. 2010; Sacco et al. 2014). Among these four, the disk orbiting MP Mus—a single, 1.2 M star lying a mere 98 pc from Earth—may be most closely analogous to the solar nebula. Meanwhile, the T Cha disk is the only nearby transition disk that is viewed at a high inclination (73°) and furthermore displays evidence for an embedded planet (Hendler et al. 2018).

The most recent comprehensive study of the membership of epsilonCA (M+13) was reliant on pre-Gaia stellar kinematics and photometry. With the benefit of Gaia DR2 data, we have revisited and refined the membership, kinematics, and color–magnitude distribution of epsilonCA, with the primary goal of firmly establishing its position along the NYMG age sequence. We have also reexamined the multiplicity of epsilonCA members and the detailed spatial distribution of disk-bearing stars in the group. In Section 2, we discuss the sample selection and flags and caveats of the available Gaia DR2 data. In Section 3, we describe our methods, including determination of an empirical isochrone and kinematics for the epsilonCA. In Section 4, we present the results from this analysis. In Section 5, we discuss the implications of these results, including the age of epsilonCA relative to other young NYMGs.

2. Sample Selection

Table 1 lists all (65) stars we considered for epsilonCA membership for which Gaia DR2 data is available. The majority of these epsilonCA candidates were drawn from M+13. Those authors gathered all proposed members of epsilonCA from the previous literature (52 in total) and assessed membership via their (pre-Gaia) proper motions and their spectral properties (such as Li absorption-line equivalent widths and infrared excesses). The set of 35 stars initially considered here as bona fide epsilonCA members then consists of M+13's "confirmed members" (M+13 Table 7). Ultimately, we also considered the 6 "provisional members" and 11 "rejected stars" from the M+13 study (Section 4.1) as well as 15 stars from Gagné & Faherty (2018, henceforth GF18) that GF18 designated as potential epsilonCA members (Section 4.1.4).

Table 1. Gaia Data and RVs for Candidate epsilonCA Members

Name a R.A. (J2015.5)Decl. (J2015.5) π σ G σ G − RP σ PMRA σ PMDec σ RV σ RUWE b RV ref c
 (deg)(deg)(mas) (mag) (mag) (mas yr−1) (mas yr−1) (km s−1) 
HD 8279*142.0874−78.25978.5040.0258.927600.000380.37370.0025−28.2110.04320.3500.04612.760.21.052
CP-68 1388 164.4553−69.23338.5920.02910.02290.00240.67630.0070−34.8420.0513.5520.04415.91.00.9312
VW Cha*167.0058−77.707912.1870.0201.4250.018
TYC 9414-191-1 169.1200−78.42241.5640.02210.5236740.000250.867180.00067−38.3720.0402.4600.0330.91
2MASS J11183572–7935548 169.6478−79.598610.570.1513.86870.00271.37260.0043−41.760.264.980.2419.31.63.35+ 8
RX J1123.2–7924*170.7311−79.41239.3790.02112.86740.00161.00360.0041−31.6700.040−17.4430.0382.72.91.148
HIP 55746*171.3217−84.954511.0770.1387.534350.000440.35240.0021−46.190.1611.920.1620.91.21.463
2MASS J11334926–7618399*173.4550−76.31115.4100.06015.567040.000811.33730.0010−22.410.11−0.5610.0871.07
RX J1137.4–7648*174.3794−76.799711.9880.02113.526530.000791.04940.0022−60.5440.041−8.9350.03214.05.01.168
2MASS J11404967–7459394 175.2063−74.994210.3970.06916.42910.00181.44660.0044−42.740.13−2.280.1010.31.00.948
2MASS J11411722–7315369# 175.3211−73.260210.4160.05215.478040.000841.32290.0024−42.0470.0935.6720.0821.07
TYC 9238-612-1*175.3646−73.78416.1240.02110.59760.00120.56630.0035−35.7470.04110.9380.0361.33
2MASS J11432669–7804454 175.8603−78.07935.540.5614.89310.00721.3210.022−43.980.92−6.520.7515.6110.51+ 8
2MASS J11432968–7418377# 175.8729−74.310510.0630.03914.176330.000661.22130.0022−40.6210.073−3.0260.0671.19
RX J1147.7–7842 176.9496−78.6989.8940.03912.35230.00121.18400.0028−41.6600.075−4.2650.06516.10.91.288
RX J1149.8–7850 177.3818−78.85039.9180.02511.93240.00431.0280.015−41.8760.047−4.2650.04013.41.31.198
RX J1150.4–7704*177.6171−77.07736.5510.03111.52500.00160.77730.0042−42.7440.054−11.0470.0401.26
RX J1150.9–7411 177.6868−74.187110.630.1813.45670.00501.3240.012−39.500.32−8.080.2915.01.26.01+ 8
2MASS J11550485–7919108 178.7692−79.31989.8860.05814.81800.00171.35020.0053−41.180.13−4.3360.0861.08
T Cha 179.3054−79.35889.1220.08312.970.111.610.35−42.000.12−9.2450.0806–300.912,5,10
RX J1158.5–7754B 179.6109−77.91269.6620.03513.20230.00271.17130.0081−39.5730.064−5.6860.06013.02.01.2211
RX J1158.5–7754A 179.6165−77.90839.5180.0359.97980.00110.78880.0033−39.6600.063−12.8440.07014.021.821.662,8,12
HD 104036 179.646−77.82549.5660.0386.692960.000530.16460.0041−41.2860.065−7.7620.07312.60.51.233
CXOU J115908.2–781232 179.7824−78.20899.4250.06215.25630.000631.31390.0022−38.620.14−5.1810.09115.10.21.248
epsilon Cha AB 179.9056−78.22194.78160.0080−0.15050.0086133.711
RX J1159.7–7601 179.9254−76.02410.0250.02310.80950.00230.76940.0058−41.0250.043−6.1900.03813.03.70.962,8
2MASS J12000269–7444068# 180.0105−74.735210.0670.04914.32410.00111.24310.0028−42.0470.079−5.6990.0711.24
HD 104237A 180.0204−78.1939.2260.0586.54270.00180.22130.0074−39.310.11−6.2120.08313.520.391.554
HD 104237D 180.0336−78.19439.8850.06113.02060.00151.19810.0065−38.870.13−3.1950.0981.53
HD 104237E 180.0379−78.19519.7960.03411.8380.0281.160.11−42.9070.074−4.4180.0511.54
2MASS J12005517–7820296 180.2291−78.34159.7130.08215.58260.00121.40050.0040−40.590.14−4.960.1210.71.31.178
HD 104467 180.412−78.988110.140.178.42850.00220.51450.0059−41.110.26−5.410.2612.810.965.172
2MASS J12014343–7835472 180.4302−78.59659.530.1117.1150.0211.1660.069−41.320.19−6.330.1420.00.61.288
USNO-B 120144.7–781926 180.4343−78.32419.8190.06315.27580.00761.3720.020−41.430.11−6.1050.08814.91.11.088
CXOU J120152.8–781840 180.4681−78.31159.7350.06814.98210.00121.33190.0026−40.590.12−6.880.1016.51.11.148
RX J1202.1–7853†180.5145−78.883710.0110.04611.51000.00301.00590.0088−45.1220.081−4.1240.06917.10.21.635
RX J1202.8–7718‡180.7269−77.31069.6060.03513.37490.00261.17620.0080−39.6250.072−6.0440.05914.40.61.217
RX J1204.6–7731†181.1498−77.52639.9220.03412.57140.00181.15130.0053−41.3640.056−6.3970.05210.42.01.206
TYC 9420-676-1*181.2384−79.53463.7000.08010.239470.000630.44430.0016−38.160.16−1.000.182.95
HD 105234 181.7721−78.74129.5660.0357.419470.000410.22470.0028−40.9710.058−9.410.0571.15
2MASS J12074597–7816064*181.9407−78.26859.4540.05814.551440.000731.22930.0021−38.4240.091−6.3190.07715.21.91.338
RX J1207.7–7953 181.9501−79.878510.0160.03013.34030.00221.18600.0057−42.0580.056−7.1130.05215.00.71.128
HIP 59243 182.2811−78.781310.0190.0336.793730.000430.18890.0034−43.3640.062−7.6190.061.10
HD 105923 182.9084−71.17679.4000.0368.836530.000830.53670.0035−38.7210.062−7.4230.05114.341.060.871,2,12
RX J1216.8–7753 184.1906−77.89279.8200.03912.92060.00371.1320.010−39.8250.072−9.0740.07214.02.01.176
RX J1219.7–7403 184.9315−74.06599.8650.02712.14000.00161.01870.0042−40.3230.043−9.2560.03913.861.890.962
RX J1220.4–7407 185.0901−74.12776.710.7011.94060.00201.07060.0050−40.91.2−4.11.412.30.438.28+ 6
2MASS J12203396–7135188# 185.1411−71.588610.7730.05713.219470.000841.24860.0019−42.5730.096−8.5710.0701.28
2MASS J12203619–7353027# 185.1502−73.884210.0070.03513.14090.00101.14240.0024−40.6210.069−9.3260.0501.12
2MASS J12210499–7116493 185.2703−71.280410.0550.02411.20320.00280.88740.0079−40.4160.043−9.6470.03411.442.531.029
2MASS J12220068–7001041# 185.5022−70.017810.0080.04714.02360.00111.23710.0045−41.2290.096−10.4740.0631.23
2MASS J12222238–7137040# 185.5926−71.617810.5160.02712.52990.00811.1050.023−41.8430.044−9.9620.0431.06
2MASS J12224862–7410203# 185.702−74.172410.6050.08916.17510.00241.47740.0068−42.570.15−9.270.131.18
2MASS J12255824–7551116# 186.4921−75.853310.0490.06515.365540.000741.33180.0026−41.230.11−10.460.101.20
2MASS J12324805–7654237# 188.1995−76.906711.690.9112.83460.00121.19600.0033−41.841.67−13.491.3236.83+
2MASS J12332483–6848553# 188.3531−68.815410.0610.08615.71180.00121.36760.0039−41.430.13−12.510.111.16
2MASS J12351540–7043079# 188.8135−70.718910.7050.06615.701480.000681.38440.0026−43.720.11−12.6390.0941.03
RX J1239.4–7502 189.838−75.04439.6460.0279.97510.00170.66790.0048−38.2140.044−12.3680.04213.622.80.979
2MASS J12421315–6943484# 190.5543−69.730210.0620.08315.14440.00111.35510.0039−39.500.14−12.770.101.24
RX J1243.1–7458*190.7219−74.987.990.3713.76520.00441.1860.011−17.250.66−2.630.5313.57.010.71+ 8
2MASS J12425584–7034207# 190.7321−70.572510.0620.07815.974960.000981.40490.0031−40.470.12−13.200.120.99
2MASS J12473611–7031135# 191.9000−70.520610.0120.03113.269200.000651.14980.0016−40.2490.048−14.2410.0451.12
CD-69 1055 194.6061−70.480410.5490.0299.60620.00270.63880.0076−41.0010.047−16.4590.04511.181.671.011,2,8
CM Cha 195.5560−76.63285.1520.02412.55990.00431.0010.012−20.9900.040−9.1440.0411.15
MP Mus 200.531−69.636810.1150.0319.95210.00300.72070.0087−38.2890.044−20.2040.04511.60.20.9812

Notes.

a M+13 bona fide members are indicated by †, M+13 provisional members are indicated by ‡, M+13 rejected members are indicated by *, and GF18 members are indicated by #. b Values marked with + have high astrometric excess noise. c References for radial velocities: (1) Desidera et al. (2015), (2) Gaia Collaboration et al. (2018), (3) Gontcharov (2006), (4) Grady et al. (2004), (5) Guenther et al. (2007), (6) Lopez Martí et al. (2013), (7) Malo et al. (2014), (8) Murphy et al. (2013), (9) Kordopatis et al. (2013), (10) Schisano et al. (2009), (11) Terranegra et al. (1999), (12) Torres et al. (2006).

Download table as:  ASCIITypeset images: 1 2 3

We searched the Gaia archive for DR2 counterparts to these 67 stars from the M+13 and GF18 studies, using a 025 search radius centered at the position listed in M+13. Unique DR2 counterparts were identified for all stars listed in Table 1, comprising a total of 65 of the 67 stars searched, though not all of these counterparts had parallaxes and colors. All Gaia counterparts met or exceeded the minimum number of Gaia visibility periods (nvis > 5) recommended by Lindegren et al. (2018).

To assess the quality of the Gaia data for individual stars, we applied three metrics: the renormalized unit weight error (RUWE; Lindegren 2018), astrometric excess noise, and photometric excess (Gaia Collaboration et al. 2018). The astrometric unit weight error (UWE), in its renormalized form (RUWE), is useful to determine when data are potentially unreliable based on Gaia measurement uncertainties alone. Following Lindegren (2018), we flag data with RUWE > 1.4. Note that this RUWE threshold is implicitly based on a star's ${G}_{{B}_{P}}-{G}_{{R}_{P}}$ color, whereas the color–magnitude diagram analysis carried out here (Section 3.1) uses Gaia $G-{G}_{{R}_{P}}$ colors. We also utilized astrometric_excess_noise to assess the quality of the astrometric data, flagging those stars with astrometric_excess_noise > 0.1 × π as having potentially large parallax uncertainties. We used Gaia's phot_bp_rp_excess_factor, E, to flag stars that did not satisfy $1.0\,+0.015{({G}_{{B}_{P}}-{G}_{{R}_{P}})}^{2}\lt E\lt 1.3+0.06{({G}_{{B}_{P}}-{G}_{{R}_{P}})}^{2}$. As discussed in Gaia Collaboration et al. (2018), stars that fall outside of this range of ${G}_{{B}_{P}}-{G}_{{R}_{P}}$ have colors that may not be trustworthy.

The RUWE values for all 65 stars are listed in Table 1. The results of the astrometric excess noise and photometric excess tests are also noted in Table 1. Two stars, T Cha and HD 104237E, fail the phot_bp_rp_excess_factor test. In Figure 1, we display the Gaia DR2 color–magnitude diagram (CMD) positions of the 30 (of 35) bona fide M+13 epsilonCA members that have viable DR2 parallaxes and photometry (see Section 3.1), with the results of the preceding data-flagging exercise indicated.

Figure 1.

Figure 1. Gaia DR2 color–magnitude diagram (CMD) for the 30 M+13 bona fide members (blue circles) included in empirical single-star locus fitting (SLFR) analysis (Section 3.1). The best-fit empirical isochrone obtained from the SLFR method is represented by the black curve. Three theoretical isochrones (Tognelli et al. 2018) for ages of 3, 5, and 8 Myr (orange, green, and red dashed lines, respectively) are also overlaid on the data. Stars with high RUWE values (low-accuracy astrometry) are marked with black, open circles; stars failing the color excess factor test (bad photometry) are marked with red squares; and stars that are outliers in magnitude offset are marked with crosses. Five stars identified as candidate photometric binaries (via the empirical single-star isochrone fitting) are denoted by gray hexagons. Errors are displayed as horizontal and vertical bars in blue; where no error bar is seen, the errors are smaller than the symbols. The pink line (with arrows) represents the reddening vector inferred for T Cha (the red square and cross at G − GRP  ∼ 1.6), i.e., E(G − GRP ) = 1.1 mag and AG  = 4.5 mag (see Section 4.1.3).

Standard image High-resolution image

3. Analysis

3.1. Empirical Single-star Isochrone

To ascertain the empirical isochrone defined by single-star members of NYMGs like epsilonCA, we have developed the single-star locus fitting routine (SLFR; D. A. Dickson-Vandervelde et al. 2020, in preparation). The SLFR utilizes a recursive method to fit a polynomial to the color–magnitude diagram (CMD) of a set of NYMG members, identifying and rejecting outliers at each iteration. Between iterations, any star with a magnitude greater than 2σ is rejected before refitting the polynomial. The final result is a polynomial that represents the best fit to the single-star locus of the group across well-sampled regions of color–magnitude space, as well as a list of candidate photometric binaries, i.e., stars that lie significantly above, but within ∼0.75 mag of, the single-star locus. The resulting best-fit polynomial then effectively represents the group's empirical single-star isochrone.

Only 32 of the 33 bona fide members of epsilonCA listed in Table 1 have the DR2 parallaxes and photometry necessary for SLFR analysis; we further excluded two stars whose parallaxes appear incompatible with epsilonCA membership (see Section 4.1.2). Figure 1 illustrates the results of the SLFR method as applied to the remaining 30 stars. The SLFR-generated empirical isochrone (polynomial) is only well fit where there is a good sampling of stars in color–magnitude space and is particularly unreliable (and is therefore not plotted) in the bluest region of the CMD (G − GRP  < 0.5). We found a fourth-order polynomial to be sufficient to match the color–magnitude data without introducing artifacts in poorly sampled regions. The final polynomial is given by G = 1.28c4 + 8.44c3 − 28.9c2 + 33.17c + 7.18, where G is the absolute G magnitude and c is the G − GRP color. This fit provided an rms in G of 0.61 mag. Hereafter, we refer to this polynomial as the empirical (5 Myr) isochrone.

The empirical isochrone is evidently a better match to the epsilonCA color–magnitude data than any of the theoretical isochrones (from Tognelli et al. 2018) plotted in Figure 1. While the 5 and 8 Myr theoretical isochrones follow the SLFR-generated empirical isochrone in the blue CMD regions (i.e., G − GRP in the range 0.2–0.8), both of these isochrones fall below the NYMG distribution for redder colors (G − GRP  > 0.8). This general behavior, in which theoretical isochrones fall under the single-star loci of NYMGs in Gaia-based CMDs, has been well documented in the recent literature (e.g., Gagné et al. 2018a). Nonetheless, it is apparent that the 5 Myr isochrone appears to be the best overall match to the SLFR-generated empirical isochrone. We further discuss the implications of the comparison between empirical and theoretical isochrones for the age of the epsilonCA in Section 5.3.

The SLFR polynomial fitting exercise yields five photometric binary candidates among the 30 bona fide M+13 epsilonCA members plotted in Figure 1. Of note, the majority of these binary-star candidates have high RUWE. This is consistent with the results of Belokurov et al. (2020), who found that stars along the binary-star locus within all Gaia data show higher RUWEs than stars along the single-star locus. This correlation is likely a result of marginally resolved, close binaries yielding low-precision astrometric solutions. The SLFR method also flagged two clear CMD outliers among the 32 bona fide members, T Cha and 2MASS J12014343–7835472, both of which are discussed in Section 4.1.3.

3.2. Kinematic Analysis

We calculated the heliocentric space motions and positions of the (35) stars considered as epsilonCA members by M+13 (Section 2) for which proper motions, parallaxes, and radial velocities (RVs) are available (see Table 1). Gaia DR2 provides RVs for eight stars, and the remaining (majority) of the RVs come from M+13 and other RV surveys in the literature (see footnote b of Table 1). Heliocentric velocities (UVW) and positions (XYZ) were calculated using code from BANYAN Σ (Gagné et al. 2018c). In cases where a star had multiple RV measurements, we adopt the mean RV to calculate UVW, except in the case of T Cha (see Section 4.1.3). The resulting UVW were used to reevaluate epsilonCA membership (Section 4.1).

4. Results

4.1. Membership

4.1.1. Kinematic and Photometric Inclusion/Rejection Parameters and Criteria

We employed two quantitative criteria, kinematic offset (Kkin) and magnitude offset (ΔM), to help assess epsilonCA membership. These two metrics, listed in Table 2, correspond to those defined in M+13. Specifically, the kinematic offset is defined as

Equation (1)

where (U0, V0, W0) is the mean space motion for epsilonCA as calculated from the 30 bona fide M+13 members used for the SLFR analysis (Section 3.1). We defined ΔM as

Equation (2)

where c is the color of the star, $M^{\prime} (c)$ is the absolute magnitude of the empirical isochrone at the star's color, and M is the measured absolute G magnitude of the star. Outliers in ΔM are those with ∣ΔM∣ ≥ 2σΔM and in Kkin are ${K}_{\mathrm{kin}}\geqslant 2{\sigma }_{{K}_{\mathrm{kin}}}$, where σΔM  = 0.629 mag and ${\sigma }_{{K}_{\mathrm{kin}}}=3.2$ km s−1. These standard deviations are calculated for the M+13 bona fide membership list with usable data (30 stars), i.e., omitting M+13's (17) provisional and rejected members.

Table 2. Membership Quality Flags

Name Kkin ${\sigma }_{{K}_{\mathrm{kin}}}$ ΔM σΔM StatusNotes
HD 82794.350.211.89110.0066RPreviously rejected (Section 4.1.2)
CP-68 13882.01.0−0.19750.0097M
2MASS J11183572–7935548A4.71.6−1.2000.033M
RX J1123.2–792415.82.90.95090.0065RPreviously rejected (Section 4.1.2)
HIP 557467.11.21.3830.027RPreviously rejected (Section 4.1.2)
RX J1137.4–76484.65.01.87900.0046RPreviously rejected (Section 4.1.2)
2MASS J11404967–74593944.21.00.0250.016M
2MASS J11411722–73153691.0930.012P
2MASS J11432669–780445417.94.1−0.830.24RProbable Cha Cloud (Section 4.1.2)
2MASS J11432968–74183770.86590.0090P
RX J1147.7–78421.810.90−0.65860.0097M
RX J1149.8–78501.01.3−0.00020.0097M
RX J1150.9–74112.71.3−0.8920.042M
2MASS J11550485–7919108−0.0590.015M
T Cha3.01.3−7.720.13MAnomalous CMD position (Section 4.1.3)
RX J1158.5–7754B1.32.00.24540.0105M
RX J1158.5–7754A3.61.8−0.75420.0091M
HD 1040361.870.513.55940.0092M
CXOU J115908.2–7812321.220.250.7740.015M
RX J1159.7–76011.33.70.30290.0073M
2MASS J12000269–74440680.7990.012P
HD 104237A0.690.422.210.016M
HD 104237D−0.1150.015M
HD 104237E−1.0270.036M
2MASS J12005517–78202963.41.3−0.1250.020M
HD 1044671.61.00.0100.038M
2MASS J12014343–78354725.930.654.1680.045MAnomalous CMD position (Section 4.1.3)
USNO-B 120144.7–7819260.831.110.0640.022M
CXOU J120152.8–7818402.51.10.330.0163M
RX J1202.1−78533.350.23−0.2790.013M
RX J1202.8−77180.770.610.3650.010M
RX J1204.6−77313.62.0−0.16950.0091M
HD 1052343.10470.0083P
2MASS J12074597–78160641.51.91.0290.014RPreviously rejected (Section 4.1.2)
RX J1207.7−79531.010.700.33880.0087M
HIP 592433.26970.0075P
HD 1059230.92.00.02310.0091M
RX J1216.8–77531.31.90.3040.012M
RX J1219.7–74039.23.30.24850.0075M
RX J1220.4–74072.11.4−1.090.24RProbable Cha Cloud (Section 4.1.2)
2MASS J12203396–7135188−0.2160.012P
2MASS J12203619–73530270.48540.0086P
2MASS J12210499–71164932.62.50.04820.0080M
2MASS J12220068–70010410.5470.011P
2MASS J12222238–71370400.2480.014P
2MASS J12224862–7410203−0.8220.021P
2MASS J12255824–75511160.7840.015P
2MASS J12324805–76542370.080.18P
2MASS J12332483–68485530.6130.020P
2MASS J12351540–70430790.4730.014P
RX J1239.4–75021.52.80.06720.0077M
2MASS J12421315–69434840.2350.019P
RX J1243.1–745811.17.00.270.11RPreviously rejected (Section 4.1.2)
2MASS J12425584–70342070.2700.018P
2MASS J12473611–70311350.55940.0073P
CD-69 10552.51.90.11760.0087M
MP Mus2.320.21−0.22650.0097M

Note. Key for status column: M = member; P = provisional member; R = rejected.

Download table as:  ASCIITypeset image

Figure 2, like Figure 8 in M+13, illustrates the kinematic and magnitude offsets plotted against each other. The figure includes all stars from M+13 (i.e., M+13's bona fide, provisional, and rejected members) with Gaia parallaxes and G − GRP colors. Stars with Kkin or ΔM values that place them within both 2σ boundaries, within errors, are henceforth considered to be high-probability members of epsilonCA. Stars outside these bounds were reconsidered for membership on a case by case basis. For instance, the two magnitude outliers beyond 1σK , HD 104036 and HD 104237A, are both A stars in the bluest region of the CMD. This region is not well fit by the empirical isochrone so the large ΔM is not enough to disqualify their membership. Additional cases, including stars whose DR2 parallaxes render their epsilonCA membership doubtful, are discussed in Sections 4.1.2 and 4.1.3.

Figure 2.

Figure 2. Comparison of the offsets in CMD (G vs. GRp ) space and kinematic (UVW) space used to assess epsilonCA membership (Section 4.1). Blue circles, purple triangles, and red squares represent the M+13 epsilonCA study's bona fide members, provisional members, and rejected stars, respectively. The horizontal dashed lines correspond to 2σΔM  = 1.25 mag; the vertical dashed lines correspond to ${\sigma }_{{K}_{\mathrm{kin}}}$ (3.2 km s−1) and $2{\sigma }_{{K}_{\mathrm{kin}}}$. Points appearing at 0.0 on the x-axis are stars lacking radial velocity measurements, for which kinematic offsets cannot be determined. Four outliers discussed in the text are labeled: RX J1220.4–7407 and 2MASS J11432669–7804454 in Section 4.1.2, T Cha and 2MASS J12014343–7835472 in Section 4.1.3, and RX J1202.8–7718 and CXOU J1159082–781232 in Section 4.1.4.

Standard image High-resolution image

Figure 3 presents Gaia CMDs for the full sample of M+13 stars initially considered epsilonCA members, i.e., for the stars listed in Table 1 (left panel) and for the final membership list, which is presented and discussed in Section 4.1.5 (right panel). The empirical isochrone in both panels is that found for the original 30 star data set (i.e., is the same as the SLFR curve in Figure 1). In the following, we describe the reasoning behind retaining or rejecting the >2σ outliers plotted in Figure 2.

Figure 3.

Figure 3. Gaia DR2 CMDs for the initial sample of stars considered for epsilonCA membership (Table 1; left panel) and for our final member lists, both bona fide and provisional members (Tables 3, 4); center and right panels). The empirical isochrone obtained from SLFR analysis of the original 30 bona fide members from M+13 (Figure 1) is overlaid in all three panels (blue curves). In the left panel, the blue circles, purple triangles, and red squares indicate M+13 membership candidacy (as in Figure 2); new candidate epsilonCA members identified via our wide-separation companion search (Section 4.2.2) are represented by orange diamonds; and candidates originating from GF18 are represented with green pentagons. In the center panel, our final confirmed epsilonCA members (Table 3 are indicated with brown circles and our provisional epsilonCA members (Table 4) are indicated with pink circles. This panel and the right panel also include theoretical isochrones (dashed curves) from Tognelli et al. (2018) for ages of 3.0 Myr (orange), 5.0 Myr (green), and 8.0 Myr (red). In the right panel, both provisional and confirmed epsilonCA members (Tables 3 and 4) are represented with brown circles; black dashes indicates stars with disks and black open circles indicate Gaia unresolved binaries and photometric binary candidates. Three stars of particular interest, discussed in the text, are labeled: T Cha (Section 4.1.3), 2MASS J12014343–7835472 (Section 4.1.3), and 2MASS J11550336–7919147 (Section 4.2.2).

Standard image High-resolution image

4.1.2. Candidates Rejected on the Basis of Kinematics, Photometry, and/or Parallaxes

We corroborate the rejected status of all of the eleven epsilonCA members rejected by M+13. Five of the M+13 rejects clearly maintain their nonmember status on the basis of the large Gaia-based kinematic and magnitude offsets determined here (Table 2). One rejected M+13 member with offsets <2σ, 2MASS J12074597–7816064, is unlikely to be a member given its position in the Gaia CMD (Figure 1). This star falls below our empirical isochrone, despite being a suspected spectroscopic binary (M+13). Its kinematic offset is also at the edge of acceptability (Kkin = 3.03 km s−1). We hence exclude 2MASS J12074597–7816064 from our final epsilonCA membership list. Another five M+13 rejected stars now convincingly lose membership status based on Gaia parallaxes (2MASS J11334926–7618399, RX J1243.1–7458, RX J1150.4–7704, TYC 9238-612-1, and TYC 9420-676-1). VW Cha, which was previously rejected on the basis of a pre-Gaia distance estimate, was not reconsidered in this work because no Gaia parallax was reported in DR2.

Two stars previously considered bona fide epsilonCA members by M+13, 2MASS J11432669–7804454 and RX J1220.4–7407, have CMD positions consistent with young ages and are marginally consistent with epsilonCA membership (in terms of ΔM) given our SLFR analysis, although both were flagged as photometric binaries due to their elevation above the single-star locus. Kinematically, RX J1220.4–7407 lies within the distribution of epsilonCA members (Kkin = 2.11 km s−1), while 2MASS J11432669–7804454 has a large kinematic offset (17.9 km s−1). However, these stars have Gaia parallaxes of 5.54 and 6.71 mas, respectively, a factor of 1.5–2 smaller than the mean for epsilonCA (9.81 mas). Notably, both stars have high RUWE and astrometric_excess_noise values, suggesting that their parallax solutions may not be reliable. If their parallaxes are in fact accurate, their relatively large distances would indicate they are younger, background T Tauri stars associated with the Cha star-forming region (see, e.g., Kastner et al. 2012). Their positions above the single-star locus would be consistent with their belonging to this (generally younger) Cha cloud population. Given these uncertainties, we do not include these two stars in our final roster of bona fide epsilonCA members, although their status is worth revisiting in future Gaia data releases.

Of the six provisional members identified by M+13, two can be rejected outright based on their Gaia parallaxes: TYC 9414-191-1 (π = 1.56 mas) and CM Cha (π = 5.15 mas). TYC 9414-191-1 is likely a background star with proper motions similar to those of epsilonCA members, and CM Cha (like 2MASS J11432669–7804454 and RX J1220.4–7407) appears to belong to the more distant Cha cloud T Tauri star population. The other four M+13 provisional members are discussed in Section 4.1.4.

4.1.3. Outliers Retained as Likely Members

T Cha: T Cha was previously established as a kinematic member of epsilonCA (Torres et al. 2008). The kinematic analysis of T Cha is more complicated than other epsilonCA members, however, due to its variable RV; measurements range from ∼6 to ∼30 km s−1 (Schisano et al. 2009). We adopt the RV measured by Guenther et al. (2007), 14.0 ± 1.3, which is near the median of the values reported by Schisano et al. (2009); with this RV, the UVW of T Cha (presented in Section 4.1.5) yield a Kkin of 3.0 ± 1.3 km s−1, well within the 2σ region. Adopting a larger value for its RV (e.g., the Gaia DR2 RV, 25.52 ± 4.24) would make T Cha a kinematic outlier. However, its XYZ (52.7, −90.1, −31.4 pc) positions T Cha near the median of epsilonCA. We conclude that the Gaia DR2 astrometric data for T Cha support its membership in epsilonCA, although the variable nature of its RV still remains to be characterized (see discussion in Schisano et al. 2009).

The Gaia DR2 photometric data for T Cha (Figure 1) show it to be far redder and fainter than expected for an epsilonCA star of its spectral type and mass (G8 and 1.3 M, respectively; Schisano et al. 2009). In DR2, T Cha was flagged as having a suspect phot_bp_rp_excess_factor factor (i.e., E = 2.2), indicative of poor-quality photometry. This is perhaps due to the star's notable variable behavior: T Cha is known to exhibit variable extinction in the optical regime, as a consequence of its highly inclined, dusty disk (i = 73°;  Hendler et al. 2018). Schisano et al. (2009) found that T Cha shows a visual extinction (corrected for ISM extinction) of ∼0.5 mag on average, but with large excursions, sometimes reaching a maximum of ∼3 mag. Variable extinction in young stars with disks, such as that displayed by T Cha, has been hypothesized to result from a nested inner/outer disk structure wherein a warped inner disk causes quasi-periodic occultations of the photosphere (Alencar et al. 2010). Modeling of the disk around T Cha shows that instead, its variability could be caused by an asymmetric, puffed-up inner disk rim (Olofsson et al. 2013).

Based on its position in the Gaia CMD (Figure 3, right) and assuming a spectral type of G8 (Schisano et al. 2009)—which suggests an absolute G magnitude of 3.3, given the SLFR-derived empirical isochrone (Figure 1)—and applying a standard ISM reddening law (Cardelli et al. 1989), we infer an AV of 6.0 mag at the time of Gaia observations. This AV is a factor of ∼10 larger than the typical AV reported by Schisano et al. (2009)—a surprising result, given that (as of DR2) Gaia data included nine viable visibility periods. Alternatively, if we were to accept the G − GRP of T Cha at face value, then the star displays a color excess E(G − GRP ) = 1.11, which would imply AV  ∼ 3.5 assuming T Cha's obscuration follows the same (Cardelli et al. 1989) standard reddening law. This is still a factor of ∼5 larger than typical for T Cha, and is furthermore discrepant with AV as determined from its spectral-type-based absolute magnitude. Thus, while T Cha's anomalous position in the epsilonCA CMD is likely due in part to obscuration by its disk, and suggests the disk dust exhibits non-ISM-like reddening, we caution that its position (low and red) in the CMD may also be a consequence of suspect Gaia photometry.

2MASS J120143437835472: This early-M star (M2.25; Luhman 2004), also known as epsilonCA 11, is notably underluminous for a star in this spectral type regime (ΔM ∼ 4 mag; Figure 1). Previous studies of epsilonCA have proposed that 2MASS J12014343–7835472 is orbited by a nearly edge-on disk (Luhman 2004), such that the cold, outer parts of the disk flare to block the stellar photospheric emission, which is then only detected in the form of scattering off the inner disk (Fang et al. 2013). The star also appears to be actively accreting disk material, given its large Hα equivalent width (Luhman 2004).

Although 2MASS J12014343–7835472 has a large Kkin value (5.7 km s−1), this is mainly due to its low measured RV (11.44 ± 2.53 km s−1); its Gaia DR2 proper motions are generally consistent with those of other epsilonCA members. Furthermore, its lithium absorption-line (λ6708) equivalent width is consistent with the young age of the epsilonCA (M+13). Its membership in epsilonCA is hence supported by the available data, provided its photospheric emission is indeed strongly attenuated by a nearly edge-on disk. Unlike T Cha, the $G-{G}_{{R}_{P}}$ color of 2MASS J12014343–7835472 appears consistent with its spectral type (Figure 1), suggesting that the occulting disk dust has a significant component of large grains even as the Gaia photometry is dominated by starlight scattered off the disk surface.

4.1.4. Other epsilonCA Candidates

Provisional M+13 candidates: Two stars designated as provisional epsilonCA members in the M+13 epsilonCA study, RX J1202.8–7718 and CXOU J115908.2–781232, are upgraded here to bona fide members. RX J1202.8–7718 was identified as a kinematic member in M+13 but its status was uncertain, with the possibility remaining that it could belong to the Lower Centaurus Crux (LCC) subgroup of the Scorpius-Centaurus OB association. In our analysis, RX J1202.8–7718 has both low Kkin (0.31 km s−1) and ΔM (0.39 mag), suggesting that it is unlikely to be an interloper from the LCC. Likewise, on the basis of its pre-Gaia proper motion and distance, M+13 suspected that CXOU J115908.2–781232 may be associated with the Cha I cloud. However, Gaia DR2 places the star at ∼105 pc, consistent with epsilonCA and much nearer to Earth than Cha I. Given its low Kkin (0.31 km s−1) and ΔM (0.82 mag), we include CXOU J115908.2–781232 among our bona fide epsilonCA member list.

Two other M+13 provisional members, HD 105234 and HIP 59243, lack the requisite RVs to determine their UVW and, hence, kinematic offsets, and are thus also classed as provisional members in our final membership list. Therefore, of the six provisional members presented by M+13, two are rejected, two are upgraded to bona fide, and two maintain their provisional status.

Gagné & Faherty ( 2018 ) candidates: On the basis of a statistical analysis of Gaia DR2 data (utilizing BANYAN; Gagné et al. 2018c), GF18 identified a total of 15 stars that are candidate members of epsilonCA. These stars were designated as high-probability members (five stars), possible members (five), and low-probability members (five), based on BANYAN probabilities of membership in epsilonCA as well as in other young moving groups. Their XZ distribution (Figure 4, left panel) makes apparent that the GF18 high-probability members are found within the spatial locus of stars we designate as bona fide epsilonCA members—albeit outside the core epsilonCA group (centered at X ∼ 49 pc, Z ∼ −28 pc), because the GF18 study excluded stars with distances >100 pc, whereas the lower-probability members are found at the upper XZ periphery of the group.

Figure 4.

Figure 4. Individual components of heliocentric positions (XZ) and the Gaia CMD for both the final membership list of epsilonCA (Table 3; brown circles) and GF18 epsilonCA candidates. Yellow squares, pink triangles, and green pentagons indicate GF18's high-probability, possible, and low-probability members, respectively. In the left panel, the heliocentric median and mean positions for epsilonCA members (Table 7) are indicated as green and orange plus signs, respectively. The green shaded region represents the area within the inferred tidal radius (∼4 pc; Section 5.2) as centered on the median position.

Standard image High-resolution image

In contrast, there is no such clear distinction between the CMD positions of the three categories of GF18 candidates. All 15 stars lie along the single-star CMD locus of epsilonCA determined via our SLFR analysis, although most (13 of 15) lie below the empirical isochrone, by up to a magnitude (Figure 4, right panel). While this might cast doubt on their epsilonCA membership, we caution that most of the GF18 stars lie in a region of CMD space that is poorly sampled by kinematically verified epsilonCA members. Hence, the final membership status of all 15 stars is not easily assessed and must await measurements of, e.g., RVs and Li absorption-line strengths.

To assign spectral types to the provisional epsilonCA members gleaned from the GF18 study, as well as other provisional members (Section 4.1.5), we have determined an empirical relationship between spectral subtype and Gaia color for bona fide K- and M-type epsilonCA members. For this purpose, we selected those K and M stars from our final membership list (Section 4.1.5) for which spectral types were determined from optical spectroscopy and for which E(B − V) ≤ 0.05, based on data presented in M+13. The resulting empirical relationship is illustrated in Figure 5. The best-fit second-order polynomial shown in the Figure is given by

Equation (3)

where Sp is the spectral index of the star, defined such that (...−1, 0, +1, ...) = (...K7, M0, M1, ...), and c is the G − GRP color. Figure 5 demonstrates that Equation (3) provides an accuracy of roughly a subtype for K and M stars in the epsilonCA. The spectral types determined from this empirical color–subtype relationship for the provisional members from GF18 are, in most cases, 0.5–2.0 subtype later than those reported by GF18. Our determinations are likely to be more reliable, however, given that they are based on the Gaia colors of bona fide epsilonCA members of known spectral type, and the color–subtype relationship is known to be age-dependent (e.g., Pecaut & Mamajek 2013).

Figure 5.

Figure 5. Spectral index vs. Gaia G − GRP color for all M+13 members with E(B − V) ≤ 0.05. Spectral index is here defined such that (...−1, 0, +1, ...) = (...K7, M0, M1, ...). Most stars have inferred E(B − V) = 0.0; two stars with E(B − V) = 0.05 are marked with red circles. The black curve represents the best-fit second-order polynomial, and the dark and light gray curves represent deviations of 0.5 and 1.0 subtype, respectively.

Standard image High-resolution image

4.1.5. Final epsilonCA Membership; Provisional Members

Our final epsilonCA membership list, including spectral types, distances, Gaia photometry, and 2MASS photometry (Skrutskie et al. 2006), is presented in Table 3. Most epsilonCA candidates without RVs are here designated as provisional members—the exceptions being three systems lacking RVs that were classified as bona fide members by M+13, for which we find Gaia DR2-based XYZ and CMD positions that are consistent with membership. The provisional members are presented in Table 4. The latter list includes three new members revealed via a search for wide-separation companions (Table 5; see Section 4.2.2).

Table 3. Final epsilon Cha Association Membership List a

NameSpTDistance (pc) b GBP G GRP J H K Multiplicity c Disk? d
CP-68 1388K1 ${115.66}_{-0.38}^{+0.39}$ 10.56810.0239.3478.488.017.79
2MASS J11183572−7935548AM4.5 ${94.14}_{-1.30}^{+1.34}$ 15.51313.86912.5010.509.899.62P
2MASS J11183572−7935548B* ${93.04}_{-3.35}^{+3.61}$ 16.189
2MASS J11404967−7459394M5.5 ${95.68}_{-0.64}^{+0.64}$ 18.83316.42914.98312.6812.1511.77
RX J1147.7−7842M3.5 ${100.52}_{-0.39}^{+0.39}$ 13.81412.35211.1689.528.868.59P
RX J1149.8−7850M0 ${100.28}_{-0.25}^{+0.25}$ 13.041911.93210.9049.458.728.49Y
RX 1150.9−7411M3.7 ${93.61}_{-1.56}^{+1.62}$ 14.85213.45712.13310.389.789.48VP
2MASS J11550485−7919108M3 ${100.60}_{-0.59}^{+0.60}$ 16.82614.81813.46811.2210.4710.08CY
T ChaG8 ${108.98}_{-0.97}^{+0.99}$ 13.68212.97411.3688.967.866.95Y
RX J1158.5−7754BM3 ${102.92}_{-0.37}^{+0.37}$ 14.62913.20212.03110.349.729.44
RX J1158.5−7754AK4 ${104.47}_{-0.38}^{+0.39}$ 10.6889.9809.1918.227.567.40VP
HD 104036A7 ${103.95}_{-0.41}^{+0.42}$ 6.8066.6936.5286.296.226.11V
CXOU J115908.2−781232*M4.75 ${105.50}_{-0.68}^{+0.69}$ 17.11815.25613.94212.0111.4511.17Y
epsilon Cha AaAbBB9111H4.8624.7824.9315.025.044.98V
RX J1159.7−7601K4 ${99.21}_{-0.23}^{+0.23}$ 11.49110.80910.0409.148.478.30
HD 104237CM/L
HD 104237BK/M
HD 104237AA7.75 ${107.76}_{-0.67}^{+0.68}$ 6.7106.5436.3215.815.254.59Y
HD 104237DM3.5 ${100.61}_{-0.62}^{+0.62}$ 14.44013.02111.8239.628.748.12?
HD 104237EK5.5 ${101.52}_{-0.34}^{+0.35}$ 12.65311.83810.675≥9.10≥8.257.49PY
2MASS J12005517−7820296M5.75 ${102.39}_{-0.85}^{+0.87}$ 17.86515.58314.18211.9611.4011.01SY
HD 104467G3 ${98.11}_{-1.56}^{+1.62}$ 8.7968.4287.9147.266.976.85C
2MASS J12014343−7835472M2.25 ${104.334}_{-1.13}^{+1.16}$ 18.20517.11515.94914.3613.3812.81Y
USNO-B 120144.7−781926M5 ${101.29}_{-0.64}^{+0.65}$ 17.24815.27613.90411.6811.1210.78Y
CXOU J120152.8−781840M4.75 ${102.16}_{-0.70}^{+0.71}$ 16.93914.98213.65011.6311.0410.77
RX J1202.1−7853M0 ${99.36}_{-0.45}^{+0.46}$ 12.54711.51010.5049.2158.468.31V
RX J1202.8−7718*M3.5 ${103.52}_{-0.37}^{+0.37}$ 14.79513.37512.19910.519.839.59
RX J1204.6−7731M3 ${100.24}_{-0.33}^{+0.34}$ 13.9412.57111.4209.779.138.88
RX J1207.7−7953M3.5 ${99.30}_{-0.30}^{+0.30}$ 14.79313.34012.15410.439.769.57
HD 105923G8 ${105.78}_{-0.40}^{+0.41}$ 9.2458.8378.3007.677.317.18VC
RX J1216.8−7753M4 ${101.27}_{-0.40}^{+0.40}$ 14.23612.92111.78910.099.479.24
RX J1219.7−7403M0 ${100.82}_{-0.27}^{+0.27}$ 13.20112.14011.1219.759.058.86
2MASS J12210499−7116493K7 ${98.92}_{-0.23}^{+0.24}$ 12.04411.20310.3169.098.428.24
RX J1239.4−7502K3 ${103.10}_{-0.28}^{+0.29}$ 10.5319.9759.3078.437.957.78
CD-69 1055K0 ${94.32}_{-0.26}^{+0.26}$ 10.1329.6068.9678.187.707.55
MP MusK1 ${98.34}_{-0.30}^{+0.30}$ 10.5479.9529.2318.287.647.29Y

Notes.

a * = previously not considered a bona fide member (M+13). b Distances have been calculated using the inverse parallax method with zero-point corrections (Lindegren et al. 2018). c S = suspected spectroscopic binary (M+13); V = resolved binary (Briceño & Tokovinin 2017); P = possible photometric binary (this work); C = potential wide-separation companion (this work). d IR excess and/or other evidence indicative of presence of circumstellar disk (M+13 and references therein).

Download table as:  ASCIITypeset image

Table 4. Provisional epsilon Cha Association Membership List

NameSpT a Distance (pc) b GBP G GRP J H K Multiplicity c Disk? d
2MASS J11411722−7315369M4.5 ${95.51}_{-0.47}^{+0.47}$ 17.38215.47814.15512.1611.5611.31
2MASS J11432968−7418377M3.5 ${98.84}_{-0.38}^{+0.38}$ 15.74214.17612.95511.2510.6010.37
2MASS J11550336−7919147M6 ${101.26}_{-5.12}^{+5.70}$ 20.47619.92518.18915.8515.03 ≥ 12.64Y
2MASS J12000269−7444068M4 ${98.80}_{+0.48}^{-0.48}$ 15.95914.32413.08111.3710.6710.42
2MASS J12011981−7859057M5 ${101.55}_{-0.77}^{+0.78}$ 17.25115.17913.82311.7511.2010.89
HD 105234A9‡ ${103.96}_{-0.37}^{+0.38}$ 7.5707.4197.1956.876.766.68VY
HIP 59243A6‡ ${99.27}_{-0.32}^{+0.32}$ 6.9346.7946.6056.356.236.17V
2MASS J12115619−7108143M3 ${105.28}_{-0.75}^{+0.76}$ 14.53013.02111.82810.099.519.24
2MASS J12203396−7135188M4 ${92.36}_{-0.48}^{+0.49}$ 14.88013.21911.97110.149.579.27
2MASS J12203619−7353027M2.5 ${99.40}_{-0.34}^{+0.35}$ 14.48513.14111.99810.419.719.48
2MASS J12220068−7001041M4 ${99.39}_{-0.47}^{+0.47}$ 15.40214.02412.787 ≥ 10.40 ≥ 10.00 ≥ 9.90
2MASS J12222238−7137040M2 ${94.61}_{-0.24}^{+0.25}$ 13.76912.53011.4259.909.278.99
2MASS J12224862−7410203M5.5 ${93.82}_{-0.78}^{+0.79}$ 18.68816.17514.69812.2611.6911.29P
2MASS J12255824−7551116M4.5 ${98.98}_{-0.63}^{+0.64}$ 17.32415.36614.03412.0211.4911.19
2MASS J12324805−7654237M3 ${85.17}_{-6.11}^{+7.14}$ 14.31012.83511.6399.889.268.96
2MASS J12332483−6848553M5 ${98.86}_{-0.83}^{+0.85}$ 17.82915.71214.34412.2811.6811.37
2MASS J12351540−7043079M5 ${92.94}_{-0.56}^{+0.57}$ 17.85315.70114.31712.1811.6611.33
2MASS J12421315−6943484M5 ${98.85}_{-0.80}^{+0.82}$ 17.16715.14413.78911.7611.1510.87
2MASS J12425584−7034207M5 ${98.86}_{-0.75}^{+0.77}$ 18.26115.97514.57012.3211.7911.46
2MASS J12473611−7031135M3 ${99.35}_{-0.30}^{+0.31}$ 14.62813.26912.11910.509.869.71

Notes.

a Spectral types determined. b Distances have been calculated using the inverse parallax method with zero-point corrections (Lindegren et al. 2018). c S = suspected spectroscopic binary (M+13); V = resolved binary (Briceño & Tokovinin 2017); P = possible photometric binary (this work). d IR excess and/or other evidence indicative of presence of circumstellar disk (M+13 and references therein).

Download table as:  ASCIITypeset image

Table 5. New Wide-separation Companions to epsilonCA Members a

NameSpTR.A.Decl. π G G − GRP PMRAPMDec
  (deg)(deg)(mas)(mag)(mag)(mas yr−1)(mas yr−1)
2MASS J11550485−7919108M3178.7692−79.31989.886 ± 0.05814.8181.35−41.18 ± 0.13−4.34 ± 0.09
2MASS J11550336−7919147M6178.7628−79.32089.82 ± 0.5319.921.74−39.74 ± 1.22−4.66 ± 0.68
HD 104467G3180.4120−78.988110.14 ± 0.178.4280.515−41.11 ± 0.26−5.41 ± 0.26
2MASS J12011981−7859057M5180.3316−78.98499.794 ± 0.07515.1791.3562−41.99 ± 0.12−5.35 ± 0.10
HD 105923G8182.9084−71.17679.400 ± 0.0368.8370.537−38.720 ± 0.060−7.42 ± 0.05
2MASS J12115619−7108143M3182.9836−71.13749.44 ± 0.6813.0211.1928−38.64 ± 0.11−8.11 ± 0.09

Note.

a Candidate wide-separation comoving systems are listed as pairs, with the previously identified epsilonCA member listed first and its candidate wide-separation companion listed second.

Download table as:  ASCIITypeset image

For each star in Tables 3 and 4, we also indicate stellar multiplicity and the presence/absence of evidence of a circumstellar disk. Multiplicity is further subdivided into photometric binaries (P), visual binaries from Briceño & Tokovinin (2017; V), suspected spectroscopic binaries from M+13 (S), and candidate wide-separation companion systems (C). Binarity and multiplicity within the final membership list are further discussed in Section 4.2.

Table 6 lists the heliocentric space positions (XYZ) and velocities (UVW) of the final epsilonCA membership. The medians and (uncertainty-weighted) means of these spatial and kinematic coordinates, as calculated from the positions and velocities of bona fide members in Table 3, are tabulated in Table 7.

Table 6. Bona Fide epsilonCA Members: Heliocentric Positions and Velocities

Name X σ Y σ Z σ U σ V σ W σ
    (pc)     (km s−1)  
CP-68 138844.890.15−105.200.35−17.2090.057−10.860.39−20.690.91−8.720.15
2MASS J11183572−7935548A43.0490.605−78.81.1−28.320.40−7.800.77−23.61.3−10.370.50
2MASS J11183572−7935548B42.51.6−77.92.9−28.01.0
RX J1147.7−784247.860.19−83.820.33−28.090.11−9.540.44−20.990.75−11.230.25
RX J1149.8−785047.890.12−83.460.21−28.2300.071−10.870.62−18.81.1−10.390.37
RX 1150.9−741144.080.75−80.31.4−19.160.33−7.880.64−19.31.0−10.570.29
2MASS J11550485−791910848.510.29−83.230.49−28.980.17
T Cha52.7120.474−90.060.81−31.430.28−11.980.65−19.41.1−13.030.38
RX J1158.5−7754B49.750.18−85.890.31−27.2200.098−10.450.97−18.51.7−10.000.53
RX J1158.5−7754A50.500.19−87.180.32−27.620.10−10.060.88−18.31.5−13.710.48
HD 10403650.250.20−86.800.35−27.340.11−11.490.25−18.260.41−11.070.14
CXOU J115908.2−78123251.080.33−87.830.57−28.400.19−9.450.16−20.270.18−10.3980.076
RX J1159.7−760147.870.11−83.770.19−23.1160.053−10.41.8−18.73.1−9.650.86
HD 104237A52.250.33−89.680.56−28.960.18−10.850.23−19.150.33−10.600.12
HD 104237D48.790.30−83.730.52−27.040.17
HD 104237E49.230.17−84.490.29−27.2900.093
2MASS J12005517−782029649.720.42−85.090.71−27.760.23−11.920.65−16.91.1−9.010.36
HD 10446747.740.77−81.11.3−27.630.45−10.390.55−18.320.81−9.650.30
2MASS J12014343−783547250.750.56−86.520.95−28.710.32−8.020.36−24.720.51−12.370.19
USNO-B 120144.7−78192649.260.32−84.160.54−27.420.18−10.020.55−20.340.92−10.610.30
CXOU J120152.8−78184049.690.34−84.880.59−27.620.19−9.020.55−21.440.92−11.360.30
RX J1202.1−785348.370.22−82.220.38−27.800.13−10.180.13−23.050.17−10.6610.067
RX J1202.8−771850.370.18−86.550.31−26.2650.094−9.850.30−19.890.50−10.110.16
RX J1204.6−773148.920.17−83.610.28−25.7610.087−11.980.98−16.671.67−9.110.52
RX J1207.7−795348.770.15−81.380.25−29.3250.089−9.880.35−20.290.58−10.920.21
HD 10592351.830.20−90.860.35−15.7360.060−9.640.52−20.670.91−8.760.16
RX J1216.8−775350.430.20−83.740.33−26.460.10−9.71.0−19.51.6−10.390.52
RX J1219.7−740350.500.14−84.980.23−19.8080.053−9.710.95−20.021.59−9.400.37
2MASS J12210499−711649349.610.12−84.310.20−14.7130.035−10.51.3−18.162.16−8.380.38
RX J1239.4−750253.600.15−85.330.24−21.7720.060−9.31.5−19.82.3−9.660.59
CD-69 105551.630.14−77.940.22−12.5010.034−9.81.0−18.71.6−8.280.25
MP Mus56.850.17−79.360.24−11.8570.036−9.130.13−19.640.17−8.5660.043

Download table as:  ASCIITypeset image

Table 7. Mean and Median Heliocentric Positions and Velocities of epsilonCA Members

 Mean a Median σb
Distance (pc)100.99100.814.62
X (pc)49.66049.6502.854
Y (pc)−84.328−84.2325.185
Z (pc)−18.905−27.3785.143
U (km s−1)−9.847−9.9541.089
V (km s−1)−20.667−19.5571.788
W (km s−1)−9.682−10.3821.305

Notes.

a Weighted means and medians as calculated from values listed in Table 6. b Standard deviation of the corresponding weighted mean.

Download table as:  ASCIITypeset image

Figure 6 illustrates the individual and mean positions and velocities. Based on our final membership, the mean distance to the epsilonCA is 100.99 pc, with a standard deviation σ = 4.62 pc. The structure and spectral type distribution of the association are discussed in Section 5.2.

Figure 6.

Figure 6. Individual components of heliocentric positions (XYZ) and velocities (UVW) plotted against one another for the final membership list (Table 3, cyan circles) and the provisional members (Tables 4, pink circles), with open black circles marking stars that are unresolved binaries and black Xs marking stars with disks. Table 3 stars lacking measured RVs (see Table 1) are omitted from the UVW plots. The mean and median values of UVW and XYZ are indicated by orange and green crosses, respectively, in each plot; the previously obtained mean values (from M+13) are represented by red crosses. The region within the calculated tidal radius (4 pc, Section 5.2) is indicated as a green circle in the XYZ plots, centered on the median. When no error bar can be seen, the error lies within the marker.

Standard image High-resolution image

4.2. Multiplicity

4.2.1. Photometric Binaries

On the basis of our SLFR analysis, we identify five photometric binary candidates, i.e., stars that lie above the single-star locus at positions consistent with their being double or perhaps triple systems: RX J1150.9−7411, RX J1158.5−775A, HD 104237E, RX J1147.7−7842, and 2MASS J11183572−7935548. These stars have ΔM values of −1.211, −0.670, −0.903, −0.767, and −1.038, respectively. Two of these systems, RX J1150.9−7411 and RX J1158.5−775A, have been resolved visually (separations 0875 and 0073, respectively; Köhler 2001; Briceño & Tokovinin 2017), but evidently, even the former is unresolved in Gaia DR2.

HD 104237E is already known to have a companion, HD 104237D (at separation 424;  Briceño & Tokovinin 2017), and would be a hierarchical triple system if HD 104237E is confirmed as a binary (see Section 4.2.3). The DR2-based binary candidacy of RX J1147.7−7842 is also novel. Both systems warrant spectroscopic follow-up to search for RV variability, since close binaries can have suspect Gaia DR2 astrometry (e.g., Kastner 2018). We note that Briceño & Tokovinin (2017) failed to detect companions at separations as small as ∼01 (∼10 au) for these two stars.

Briceño & Tokovinin (2017) found that the 2MASS J11183572−7935548 system (henceforth J1118AB) is a 092 separation binary consisting of an M4.5 primary and a lower-luminosity companion that may be related to the transition-disk nature of the object. The Gaia DR2 data confirm the angular separation (09) and establish that this corresponds to a projected physical separation of 85 au. The parallaxes for J1118AB agree, within the errors, confirming that they constitute a physical pair. Unfortunately, the only Gaia photometry for J1118B is in the G band, so we could not ascertain the CMD location of this companion. J1118A was flagged as a spectroscopic binary in M+13 and, via our SLFR method, we find it is also a candidate photometric binary. In light of the presence of the faint visual companion 2MASS 1118B (Briceño & Tokovinin 2017), it appears that 2MASS 1118AB is a possible hierarchical triple system.

4.2.2. Wide-separation Binaries

We searched for potential wide-separation companions to all epsilonCA members retained after applying the criteria described in Section 4.1. Specifically, we searched the Gaia DR2 catalog for equidistant, comoving stars by querying the catalog for all Gaia sources within a 500'' radius of the position of each of these epsilonCA members and then reordering the resulting source list by parallax. This search radius corresponds to ∼50 kau, or ∼0.25 pc, at the mean distance to epsilonCA. While a search of epsilonCA member fields within this radius typically returns ∼2000 stars, once ordered by parallax, potential (equidistant) companions to the star originally queried rise to the top of the list and are hence conspicuous. All stars so identified have parallaxes and proper motions within a few percent of the star searched and so were accepted as wide-separation companion candidates.

Applying this method, we have identified three new candidate members of epsilonCA. Gaia DR2 data for these candidates are listed in Table 5. The faintest and reddest of these three, 2MASS J11550336−7919147, is described below. The other two comoving companion candidates, 2MASS J12011981−7859057 and 2MASS J12115619−7108143, are mid-M stars that fall along the empirical single-star isochrone, with ΔM ≤ 2σΔM . These two candidate comoving companions have projected physical separations from their primaries of 5.7 kau and 17.7 kau, respectively. We estimated their spectral types (M5 and M3, respectively) from the empirical relationship between Gaia G − GRP color and spectral type for epsilonCA members described in Section 4.1.4. Neither star displays evidence of a dusty disk in the form of an IR excess, i.e., both have 2MASS/WISE colors consistent with those of "diskless" young-M stars of similar spectral type (Pecaut & Mamajek 2013).

This comoving companion search also recovered multiple stars already considered epsilonCA members as components of possible wide binary systems. The RX J1158.5−7754 system was matched with HD 104036. Three stars, 2MASS J12005517−7820296, CXOU J120152.8−781840, and USNO-B 120144.7−781926, were identified as a three-component comoving system. We also identified CXOU J115908.2−781232, a provisional M+13 member, as another possible component of the HD 104237 multiple system (see Section 4.2.3).

2MASS J115503367919147: The third, newly identified wide-comoving companion, 2MASS J11550336−7919147 (henceforth 2MASS J1155−7919B), at an absolute G magnitude of ∼15 and G − GRP  ∼ 1.75, is both the faintest and the reddest object thus far identified in epsilonCA (Dickson-Vandervelde et al. 2020). As described in Dickson-Vandervelde et al. (2020), we find 2MASS J1155−7919B is comoving with 2MASS 11550485−7919108 (hereafter 2MASS J1155−7919A). The star 2MASS J1155−7919A was itself initially thought to be a wide-separation comoving companion to T Cha Kastner et al. (2012), before Gaia DR2 data established that T Cha and J1155−7919A are neither equidistant nor precisely comoving (Kastner 2018). The 2MASS J1155−7919AB pair has a projected separation of 575, corresponding to a projected physical separation of 566 au. In Dickson-Vandervelde et al. (2020), we suggested that the position of 2MASS J1155−7919B at the extreme faint, red end of the single-star locus of the epsilonCA CMD reflected its likely status as a substellar object, with a bolometric luminosity ($\mathrm{log}{L}_{\mathrm{bol}}/{L}_{\odot }=-3.2$) that would imply its mass is a mere ∼10 MJup. This would make 2MASS J1155−7919B the lowest-mass epsilonCA member presently known—even less massive than WISEA J120037.79−784508.3, a recently identified brown dwarf candidate and possible epsilonCA member (Schutte et al. 2020) 5 However, a potential alternative model to explain the large absolute G magnitude and red color of 2MASS J1155−7919B is that the object is in fact a mid-M star—possibly a near-twin to host 2MASS J1155−7919A—that is obscured by a large column density of gray dust in a nearly edge-on disk. Under this interpretation, 2MASS J1155−7919B would be analogous to 2MASS J12014343−7835472 (=epsilonCA 11; 4.1.3), but even more highly obscured by its disk. Such an alternative explanation is motivated by the fact that 2MASS J1155−7919B is consistently 5 mag dimmer than its host in G, Bp , Rp , J, and H (Table 5), and that its spectral type, as obtained from its color (via the relationship in Figure 5), is M6. In a forthcoming paper (D. A. Dickson-Vandervelde et al. 2020, in preparation), we further explore these two possible scenarios for the nature of 2MASS J1155−7919B.

4.2.3. HD 104237

HD 104237 is a proposed quintuplet system within epsilonCA, consisting of a triple system dominated by the bright Herbig Ae/Be star HD 104237A and including the close binary HD 104237DE (Feigelson et al. 2003; Grady et al. 2004). Our wide-separation companion search (Section 4.2.2) also suggests that the previously identified epsilonCA member CXOU J115908.2−781232 is a possible additional companion to HD 104237A. We also flagged HD 104237E as a possible photometric binary, as noted in Section 4.2.1. While stellar components A, D, and E were all resolved within Gaia DR2, B and C were not, most likely being incorporated into the Gaia point-spread function of HD 104237A.

The majority of stars in the HD 104237 system lie near the spatial median of epsilonCA, the densest region of epsilonCA, and the system is well within the group's tidal radius (Section 5.2). However, the Gaia DR2 parallaxes of HD 104237D and E place these two stars 6–7 pc closer to Earth than HD 104237A. Indeed, during our wide-separation companion search, HD 104237D+E flagged each other, while HD 104237A only matched (in terms of parallax and proper motion tolerances) with the star CXOU J115908.2−781232. If HD 104237A+E is in fact a bound system, this would imply that the Gaia DR2 parallaxes for components D and E are spurious. On the other hand, as in the case of T Cha and J1155−7919A, the star originally designated T Cha B (Kastner et al. 2012; Kastner 2018), it is possible that the apparent D and E components of the HD 104237 system are in fact merely epsilonCA members that are seen in projection near the bright primary star. This would not be surprising, given the system's position near the median XYZ of epsilonCA. If the potential new components of HD 104237, CXOU J115908.2−781232 and HD 104237Eb, are, in fact, bound to HD 104237A, this could be a seven-star system. Alternatively, the bound components may in fact consist of a triple system comprising HD 104237Ea+Eb+D plus a quadruple system comprising HD 104237A+B+C and CXOU J115908.2−781232. Gaia DR3 should help resolve some of this uncertainty concerning the composition of the HD 104237 system by confirming and/or refining the DR2 parallaxes to its individual components.

5. Discussion

5.1. Multiplicity Fraction and Spectral Type Distribution

In light of the preceding, we can now revisit the epsilonCA multiplicity fraction (MF), or the number of systems consisting of more than one star. Previous to this work, it was found that higher mass stars of epsilonCA all have companions, while the lower-mass stars have a low companion frequency (Briceño & Tokovinin 2017). The binary fraction, including suspected binaries, was reported in M+13 as ${36}_{-8}^{+10}$%. The MF of our final membership list (Table 3) is 40%, with 12 of 30 systems in multiple systems, which is consistent within the errors with the MF determined by M+13.

We further separate the MF of epsilonCA into mass groups, bearing in mind small number statistics. For the population of intermediate-mass (A and B, n = 5), solar-type (F, G, and K, n = 11), and low-mass (M, n = 19) star systems in epsilonCA we find MF = 100%, MF = 36%, and MF = 42%, respectively. This can be compared with the MFs of intermediate-mass, solar-type mass, and low-mass stars on the main sequence, for which the MFs are on average MF ≥ 50%, MF = 44% ± 2%, and MF = 26% ± 3%, respectively (Duchêne & Kraus 2013). In contrast, the MF of T associations lies in the range ∼66%–75% (Duchêne & Kraus 2013). The MF of epsilonCA is hence less than that of T associations and greater than that of MS stars. This is consistent with the ∼5 Myr age of epsilonCA (Section 5.3), which is intermediate between T associations and the MS field. The epsilonCA hence appears to be a key group in which to study the dissolution of young binary systems as T associations evolve toward the field population.

Figure 7 presents the spectral type distribution obtained from our final membership list (Table 3). For purposes of this histogram, K stars of all subtypes are grouped together. The figure demonstrates that the epsilonCA spectral type distribution, which peaks at M4, resembles that of the χ1 For cluster (see Figure 6 in Zuckerman et al. 2019) and other nearby associations like TWA, BPMG, and Columba, many of which have far better statistics (see, e.g., Figure 9 in Lee & Song 2019).

Figure 7.

Figure 7. Spectral type histogram compiled from the final epsilonCA membership lists (Table 3 and 4). Unfilled portions of the histogram bars indicate provisional members. Non-M spectral types are grouped across subtypes, while the M subtypes are plotted individually. The color scheme follows that of Figure 9 of Lee & Song (2019).

Standard image High-resolution image

5.2. Structure of the Association

With accurate heliocentric positions, we can analyze the structure of the epsilonCA NYMG in the Galactic context. Following Zuckerman et al. (2019), the tidal radius of the group can be roughly estimated as $r=R{({M}_{c}/3{M}_{g})}^{1/3}$ (King 1962), where R is the distance between Earth and the Galactic Center (∼8200 pc), Mc is the mass of the stellar group, and Mg is the Galactic mass interior to the Sun (∼1011 M). In order to obtain a rough estimate for Mc , we adopt masses of 3.0, 2.0, 1.0, 0.7, and 0.3 M for B, A, G, K, and M stars, respectively. We then obtain an estimate for the group mass, based on the spectral types of individual bona fide members of the Association (Table 3), of Mc  ∼ 28 M. 6 Adopting this value of Mc , the tidal shredding radius of the group is found to be 4 pc, with a large uncertainty (given the large uncertainties in the estimates of the epsilonCA and Galactic masses). Given the estimated group mass and this tidal radius, and assuming (for simplicity) a spherical stellar distribution, the mass density is ∼0.13 M pc−3 and the stellar density is ∼0.18 pc−3. This mass density is within a factor of 2 of the estimated local Galactic disk density ( ∼ 0.1M pc−3; Mamajek 2016), reflecting the fact that epsilonCA is a diffuse association (as opposed to cluster). Inclusion of provisional members would increase the Association mass estimate to 36M but negligibly change the aforementioned tidal radius and density calculations.

The spatial and kinematic structure of the final membership of the epsilonCA is illustrated in Figure 6. The vast majority of group members (∼80%) fall within the tidal shredding radius of the median group XYZ position, shown as the green shaded circle and red cross, respectively. The binaries in the group fall both inside and outside of the central region of the association, with no obvious correlation with position. Figure 8 illustrates the disk fraction versus distance from the median XYZ position of epsilonCA. The figure indicates that 9 of the 11 epsilonCA systems known to host disks lie within 5 pc of the median position, a volume roughly coincident with that defined by the group's tidal shredding radius; the disk fraction falls from ∼50% within 5 pc of the median position to ∼15% beyond 5 pc of the median. This centrally concentrated distribution suggests that stars in the core region of the epsilonCA are more likely to retain dusty disks, hinting at the possibility that the nascent planetary systems orbiting stars in this region are subject to more frequent and/or more catastrophic dust-generating collisions. The distribution in Figure 8 stands in contrast to that of the (older) χ1 For cluster (age ∼40 Myr), for which the majority of stars with IR excesses (hence dusty disks) lie in a narrow, spherical annulus just outside the tidal shredding radius of the cluster (Zuckerman et al. 2019). However, as noted by those authors, χ1 For also represents a particularly striking and unusual case of a large disk frequency in a ∼40 Myr old cluster.

Figure 8.

Figure 8. Histogram of the distance from the median XYZ position for the final list of epsilonCA members (Table 3). The filled orange bars indicate stars with disks, and the disk fractions are represented as horizontal black lines. The tidal shredding radius (4 pc; Section 5.2) is denoted by a vertical gray dashed line.

Standard image High-resolution image

5.3. Age of the Association

The SLFR analysis applied to epsilonCA (Section 3.1) was additionally performed for two other NYMGs, the TW Hya Association (TWA) and β Pic Moving Group (BPMG). These NYMGs are both slightly older than epsilonCA, at 8 Myr (Donaldson et al. 2016) and 24 Myr (Bell et al. 2015), respectively, and hence provide good references for the relative age of epsilonCA as determined from their respective SLFR-based empirical isochrones. For our SLFR analysis, we used lists of stars considered bona fide members of each group (Lee & Song 2019), totaling 30 stars for TWA and 113 stars for BPMG.

The empirical isochrones for the three NYMGs are compared in the right panel of Figure 9. With the exception of BPMG, the blue ends of the empirical isochrones are not well fit, due to the small sample sizes in these regions, and are hence not plotted in Figure 9 (and are excluded from this discussion). All three of the empirical isochrones fit their respective NYMG data well in the regions 0.2 ≤ G − Rp  ≤ 1.2. In this CMD region, the empirical isochrones show the expected hierarchical pattern—i.e., with epsilonCA highest, BMPG lowest, and TWA intermediate between the two—reflecting the relative ages of these three NYMGs. At redder colors, all three groups drop off in population and appear to display a larger degree of scatter, and hence the empirical isochrones (fourth-order polynomials) do not provide as good a fit in these regions.

Figure 9.

Figure 9. Gaia DR2 CMDs for epsilonCA (Table 3; blue circles), the TW Hya Association (green squares), and the β Pic Moving Group (red diamonds), where data for the latter two groups are based on for the membership lists in Lee & Song (2019). In the left panel, the data for the three groups are overlaid with the corresponding empirical isochrones obtained from SLFR analysis (Sections 3.1 and 5.3). In the right panel, the data are overlaid with theoretical isochrones from Tognelli et al. (2018) for ages of 5, 10, and 20 Myr.

Standard image High-resolution image

Three representative isochrones from Tognelli et al. (2018) are shown in the right panel of Figure 9. 7 While the theoretical isochrones well match the Gaia NYMG data blueward of G − Rp  ∼ 0.8 (apart from the TWA, which lacks stars in this region), these curves fall below the Gaia data at regions redder than  ≃ 0.9 in G − Rp . This divergence between theoretical isochrones and data is commonly observed in Gaia CMDs, as previously noted (Section 3.1), and can likely be ascribed to the high levels of magnetic activity of late K and M stars (e.g., Gagné et al. 2018a).

Based on the close correspondence of the theoretical 5 Myr isochrone and the empirical (SLFR-generated) single-star isochrone in the region 0.5 ≤ G − Rp  ≤ 1.0, we estimate an age of ${5}_{-2}^{+3}$ Myr for epsilonCA, where the range of uncertainty is based on the vertical offset of the 3 and 8 Myr theoretical isochrones in this same domain of G − Rp . This is consistent with the age range determined via theoretical isochrone analysis in the (pre-Gaia) M+13 study, i.e., a median of 3–5 M for lower-mass stars (albeit with somewhat older inferred ages for solar-mass stars), as well as with the age obtained from a reanalysis of the M+13 data by (${3.7}_{-1.4}^{+4.6}$ Myr; Schutte et al. 2020). Future age estimates for the epsilonCA should be informed by spectroscopic observations designed to confirm the membership status of provisional members (Table 4), a DR3-based search for additional members, and application of isochrones obtained from "magnetic" pre-MS evolution models (e.g., Simon et al. 2019).

6. Summary

We have used Gaia DR2 astrometric and photometric data to refine the membership of epsilonCA and have thereby established the distance (mean D = 100.99 pc, σD  = 4.62 pc), age (${5}_{-2}^{+3}$ Myr), and spatial and kinematic distributions (Figure 6; Table 7) of the association. We confirm that epsilonCA is significantly younger than both the TW Hya Association and β Pic Moving Group (Figure 9) and, hence, that epsilonCA represents the youngest NYMG within ∼100 pc of Earth. Our analysis includes the determination of an empirical relationship between Gaia G − GRP color and spectral type that should provide an accurate means to determine the spectral subtypes of ∼5 Myr old K and M stars from dereddened Gaia photometry.

Our analysis yields a final epsilonCA membership list consisting of 36 bona fide members and 20 provisional members (Table 3; Table 4). The provisional members require follow-up study (in particular, RV determinations) in order to confirm their membership status kinematically. These include three new members of the epsilonCA that we have identified via a search of Gaia DR2 for wide-separation, comoving companions to previously identified members. One of these newly identified members is either a substellar object or an M star viewed through (and hence highly obscured by) a nearly edge-on circumstellar disk. We identified five photometric binaries among the group members, three of which are new binary-star candidates.

Like other nearby stellar groups, the presently known membership of epsilonCA is dominated by M stars, and the spectral type distribution peaks in the mid-M range. The MF of epsilonCA (40%) is intermediate between those of the field star population and T associations. We find an overall circumstellar disk fraction of 30% for epsilonCA, with the vast majority of disk-bearing stars lying within ∼5 pc (Figures 6 and 8).

In providing a well-defined region of XYZ and UVW space encompassing this NYMG, this work sets the stage for a more complete search of Gaia DR2 (and eventual DR3) data for new members of epsilonCA. This work thereby provides the framework for future investigations of the initial mass function as well as multiplicity and disk fractions of epsilonCA.

We thank the anonymous referee for numerous comments that improved this manuscript. We acknowledge helpful comments from Jonathan Gagné, Beth Klein, and Eric Mamajek. This research is supported by NASA Exoplanets Research Program grant 80NSSC19K0292 to Rochester Institute of Technology. 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 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.

Footnotes

  • 5  

    WISEA J120037.79−784508.3, which has an absolute G magnitude of 11.3, is not included in Table 4.

  • 6  

    Note that this estimate ignores the potential contribution to Mc of molecular gas possibly associated with epsilonCA—or at least with T Cha—as evidenced by detection of interstellar CO at coincident RV (Sacco et al. 2014).

  • 7  

    We found MIST isochrones (Paxton et al. 2011, 2013, 2015; Choi et al. 2016; Dotter 2016) to be indistinguishable from the Tognelli et al. (2018) curves, so only the latter are used as representative of the expected temporal behavior.

Please wait… references are loading.
10.3847/1538-3881/abd0fd