Issue |
A&A
Volume 558, October 2013
|
|
---|---|---|
Article Number | A77 | |
Number of page(s) | 21 | |
Section | Galactic structure, stellar clusters and populations | |
DOI | https://doi.org/10.1051/0004-6361/201220704 | |
Published online | 07 October 2013 |
A kinematic study and membership analysis of the Lupus star-forming region ⋆,⋆⋆
1
Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de
São Paulo, Rua do Matão 1226, Cidade Universitária,
05508-900
São Paulo
SP
Brazil
e-mail:
galli@astro.iag.usp.br
2
Institut d’Astrophysique, 98bis Bd Arago, 75014
Paris,
France
3
Observatoire Aquitain des Sciences de l’Univers, CNRS-UMR 5804,
BP 89,
Floirac,
France
Received:
8
November
2012
Accepted:
5
August
2013
Aims. A precise determination of the distance to individual stars is required to reliably determine the fundamental parameters (mass and age) of young stellar objects. This paper is dedicated to investigating the kinematic properties of the Lupus moving group of young stars with the primary objective of deriving individual parallaxes for each group member.
Methods. We identify those stars in the Lupus star-forming region that define the comoving association of young stars by utilizing our new and improved convergent point search method that allows us to derive the precise position of the convergent point of the comoving association from the stars’ proper motions. We used published proper motion catalogs and searched the literature for radial velocities, which are needed to compute individual parallaxes. We supplemented the radial velocity data with new measurements from spectroscopic observations performed with the FEROS spectrograph mounted on the MPG/ESO 2.2 m telescope at La Silla.
Results. We identify a comoving group with 109 pre-main sequence stars and candidates that define the kinematic properties of the Lupus low-mass star-forming region. We derive individual parallaxes for stars with known radial velocity and tentative parallaxes for the remaining group members by assuming that all stars share the same space motion. The convergent point method, combined with the k-NN algorithm, makes it possible to distinguish the Lupus and Upper Centaurus Lupus stars from the adjacent Scorpius-Centaurus association. We find significant depth effects in this region and show that the classical T Tauri stars, located in the close vicinity of the Lupus molecular clouds, form a background population, while the weak-emission line T Tauri stars are dispersed not only in angular extent but also in depth.
Conclusions. The newly derived individual parallaxes will be used in a forthcoming paper to refine the masses and ages of Lupus T Tauri stars, with the aim of better constraining the lifetimes of their circumstellar, protoplanetary disks.
Key words: stars: distances / stars: kinematics and dynamics / stars: pre-main sequence / proper motions / parallaxes / open clusters and associations: individual: Lupus
Based partly on observations collected at the European Southern Observatory, Chile (ESO Programme 087.C-0315).
Tables 4, 7, and Appendices A, B are available in electronic form at http://www.aanda.org
© ESO, 2013
1. Introduction
The inferred fundamental parameters of T Tauri stars (TTSs), the young solar-type pre-main sequence stars first discovered by Joy (1945), are sensitive to their assumed distance. TTSs are usually associated with molecular clouds, the distances of which can be estimated, say, from the photometry of a few bright stars enshrouded in reflection nebulosity (see, for example, Racine 1968) or from 2MASS extinction maps combined with known stellar parallaxes (Lombardi et al. 2008). The distances to the nearby star-forming regions (SFRs) are thus known to relatively good accuracy, providing a first estimate of the distances to the members of the stellar population associated with them. Although these average parallaxes provide valuable information, we must know the distances to the individual members of the young associations more precisely to better constrain their ages and masses by comparing observed stellar properties to evolutionary models.
TTSs are late-type objects and usually faint in the visible, suffering typically from 1 to 2 mag of extinction in that range, while the nearby SFRs to which they belong are typically in the distance range 0.1–1 kpc. Parallax determination from the ground is usually impossible for these objects, and they even presented a challenge for the Hipparcos mission (ESA 1997), which observed only a few of them. Bertout et al. (1999) computed new astrometric solutions using the Hipparcos data for groups of TTSs in various SFRs, thus providing post-Hipparcos average distances to these groups.
Some progress in determining the distance of individual TTSs has been made in recent years, using two different methods. Loinard and collaborators used the Very Long Baseline Array (VLBA) to determine parallaxes of a few selected objects, mainly in the Taurus-Auriga SFR (Loinard et al. 2007; Torres et al. 2007, 2009; Dzib et al. 2010, 2011; Torres et al. 2012). The VLBA data allow for a very precise parallax determination, but the method is observationally intensive, limiting its application to a few remarkable stars. Another approach rests on the fact that the members of young associations share the same spatial motion and uses their proper motions to determine individual parallaxes for members of some nearby associations. This was done, for example, for the TW Hydrae association by Mamajek (2005) and for the Taurus-Auriga SFR by Bertout & Genova (2006). The derived parallaxes using this method are not as precise as those obtained by VLBA observations, but it yields usable parallaxes for all members of the moving group that have measured radial velocities.
Thus, while we have progressed, we are still far from knowing the distances to the large sample of pre-main sequence stars (PMSs), which is required to study such timely topics as the effect of differing protostellar environments on disk lifetimes and the timescales of planet formation. The situation will change considerably with the launch of the Gaia mission, because its instruments will measure the parallaxes and proper motions of millions of faint stars. Although the satellite’s launch date is 2013, its catalog will be published several years later, so one would appreciate making some progress in determining the distances of TTSs in the meantime. Also, the astrometric methods that are developed and tested for this purpose will certainly be useful in the Gaia era.
In this paper, we study the largest southern SFR, located in Lupus, using the recently developed new version of the convergent point (CP) method presented by Galli et al. (2012). Section 2 is devoted to a brief presentation of the Lupus SFR and its previous distance determinations, while Sect. 3 discusses the sample of Lupus candidate stars found in recent catalogs and proper motion information needed for the CP analysis. Section 4 discusses our search for radial velocity data, including observations carried out for a number of stars in our sample with FEROS mounted on the 2.2 m MPG/ESO telescope at La Silla, Chile. The CP and membership analysis are discussed in Sect. 5, while Sect. 6 presents the parallax computations for the moving group members. We discuss the results of this investigation in Sect. 7. Finally, our conclusions are given in Sect. 8.
2. The Lupus association of young stars
The Lupus dark cloud complex is a low-mass SFR that contains four main star-forming clouds (Lupus 1 to 4) and constitutes one of the richest nearby associations of TTSs. Many of the Lupus TTSs have been identified from ROSAT X-ray observations and from the spectroscopic surveys conducted by Krautter et al. (1997) and Wichmann et al. (1997a,b). These surveys showed that the Lupus molecular clouds are surrounded by an extended halo of late-type, X-ray active stars, which bear resemblance to the so-called weak emission-line TTSs (WTTSs), the X-ray active TTSs without circumstellar accretion disks (Walter 1986). If the ROSAT-detected objects were all WTTSs belonging to the Lupus association, then they would greatly exceed the number of so-called classical TTSs (CTTSs), the TTS subgroup showing evidence of circumstellar accretion disks (Bertout 1987; Bertout et al. 1988). Lupus CTTSs were identified by Schwartz (1977) on the basis of their association with the molecular clouds and their Hα emission. This early census of Lupus young stars has recently been expanded by infrared observations with the Spitzer Space Observatory (Merín et al. 2008), and a proper motion study of these objects was performed by López Martí et al. (2011).
The focus of most kinematic studies over the past century was not the Lupus clouds, but the adjacent Scorpius-Centaurus (Sco-Cen) association that contains the nearest OB association. It is divided into the three subgroups Upper Scorpius (US), Upper Centaurus-Lupus (UCL) and Lower Centaurus-Crux (LCC). The Lupus molecular clouds, which are probably signposts of a more recent episode of star formation in that region (see Comerón 2008), occupy a gap between US and UCL. The distance of the Sco-Cen subgroups was investigated by de Zeeuw et al. (1999), who assessed association membership again using the Hipparcos data and derived a mean distance of 140 pc for the UCL subgroup. The Lupus SFR was first estimated to be at the same distance as UCL (Hughes et al. 1993). The close proximity of the UCL subgroup makes it difficult to distinguish between members of the Lupus association and of the UCL association. Mamajek et al. (2002) conducted a spectroscopic survey of UCL candidate members and selected 56 stars that had the same proper motion as UCL association members. A more complete list containing 81 UCL candidate stars is given in the recent review on the Sco-Cen association by Preibisch & Mamajek (2008). As it turned out, many of these stars were previously reported as WTTS members of the Lupus association by Krautter et al. (1997) and Wichmann et al. (1997a,b). Since the membership status of these objects remains unclear, a more detailed study is clearly necessary.
In the past decade, several works have cast doubt on the distance to the Lupus SFR. From the angular extent of the molecular clouds (about 15deg), and by assuming that the depth of the molecular region is comparable, the distances of individual association members are expected to range from about 110 pc to 190 pc. Two bright TTSs observed by Hipparcos, RY Lup, and V856 Sco are apparently associated with clouds of the Lupus SFR, and their parallaxes are compatible with the above range, but the cloud distances reported in the literature range from 100 pc (Knude & Hog 1998) to 360 pc for Lupus 2 (Knude & Nielsen 2001).
Bertout et al. (1999) used five stars connected with
the Lupus complex that had been observed by Hipparcos to compute the average
parallax of the Lupus SFR. They found a distance equal to
206
pc, which is a high value compared to the previous estimates, but noticed that the group
parallax determination was dominated by the brightest stars HIP 79080 and 79081, two
components of a Herbig Ae/Be (HAeBe) system located in Lupus 3. Considering only the three
fainter stars HIP 77157, 78094, and 78317, located in Lupus 1, 2, and 4, respectively, they
obtained a distance of 147
pc, in agreement with previous determinations of the Lupus association distance, but with a
large uncertainty, and a value of 228
pc for HIP 79080 and 79081. Bertout et al. (1999)
concluded that either HIP 79080 and 79081 are not members of Lupus 3, which appears
unlikely, or that this cloud is farther away than the other subgroups. Lombardi et al. (2008) conclude from their detailed study of 2MASS
extinction maps that Lupus has a depth of 51
pc and that this “might be the result of different Lupus subclouds being at
different distances”. The kinematic study presented hereafter aims to shed light
on the structure of the Lupus clouds, as well as on the membership of the Lupus association
of PMS stars.
Rejected stars from our sample of Lupus candidate members.
3. Sample of Lupus candidate members with known proper motions
3.1. Proper motions properties
The Ducourant et al. (2005, hereafter D05) proper motion catalog for PMS stars contains 197 stars in the general area of the Lupus SFR, which approximately ranges from 325° ≤ l ≤ 342° in Galactic longitude and 0° ≤ b ≤ 25° in Galactic latitude. The stars listed in this catalog were mostly identified by Herbig & Bell (1988), Krautter et al. (1997), and Wichmann et al. (1997a,b). We then included in the sample 37 stars of the comprehensive review performed by Comerón (2008) that were not considered in D05, and another 24 stars from the c2d Spitzer Legacy Program detected in the recent paper of López Martí et al. (2011, only group A sources of their study). Our first list of presumed members of the Lupus association therefore consists of 258 stars.
To access the more recent measurements, we searched for proper motion data in the PPMXL (Roeser et al. 2010), SPM4 (Girard et al. 2011), and UCAC4 (Zacharias et al. 2012) catalogs for all the 258 stars in our sample. The main source of proper motions is the SPM4 catalog, which represents the best present-day compromise between proper motion precision and target coverage for the stars in our sample. For the brightest stars (V ≤ 12) in the sample we also searched the Tycho2 (Høg et al. 2000) catalog and used these values when the proper motions given in SPM4 were of lower precision. Doing so, we found proper motion information for 241 stars in our initial sample, and for the remaining ones we kept the proper motions given in D05. We used the proper motions of the PPMXL and UCAC4 catalogs for two stars, since these are the only values available in the literature.
Whenever dealing with the kinematics of SFRs, the high fraction of binaries and multiple systems plays an important role. The overall binary fraction in Lupus is expected to reach about 30% to 40% (Ghez et al. 1997; Merín et al. 2008). An advantage of the D05 catalog is that these systems are clearly identified with a mention AB indicating that the given proper motion is representative of the binary system. In these cases we decided to use the proper motion given in D05 for both resolved and unresolved binaries rather than taking the value provided by SPM4, where the existence of close companions is not mentioned.
The CP search method selects cluster members and defines the CP of a moving group based on proper motion data. Therefore, before starting our analysis, it is necessary to reject all stars with proper motion that carry poor information because of measurement errors. From the sample we reject 33 stars whose proper motion is dominated by errors (i.e., σμα,δ ≥ μα,δ) in both components. After a 3σ elimination in both proper motion components we reject another eight stars that are in obvious disagreement with the common streaming motion of the Lupus moving group. These stars are given in Table 1. The remaining 217 PMS stars will be used in this paper to investigate the kinematic properties of this SFR and discuss their membership status based on the CP analysis. Their average proper motion is (μαcosδ,μδ) = (−16, −21) mas/yr with an average precision of about 2 mas/yr in each component.
Figure 1 displays all proper motion values and their associated uncertainties. We note that a small group of these stars displays small proper motions. Because we considered only stars related to the Lupus SFR whose PMS status had already been established in previous works, it appears unlikely that field stars pollute our sample. We therefore assume for the moment that these stars form a background population belonging to the Lupus complex, and we come back to this point in Sect. 7. Figure 2 shows the proper motion vectors for our sample of Lupus stars. Although the position of the CP is not clearly apparent in this plot, one notices that most proper motion vectors point in a common direction. The precise coordinates of the CP will be derived in Sect. 5.
![]() |
Fig. 1 Proper motion and associated errors for the 217 stars considered in this work. |
![]() |
Fig. 2 Proper motion vectors of the 217 Lupus candidate stars. |
Median positions, proper motions, radial velocities, and the number of stars for the various subgroups in Lupus.
3.2. Refining the Lupus pre-main sequence star sample
While the majority of the stars in our sample lie very close to the Lupus clouds, a fraction of the stars are scattered over a wider region (see Fig. 3). Most of these stars were identified and classified as WTTSs by Krautter et al. (1997) and Wichmann et al. (1997a,b) from ROSAT X-ray pointed observations and the ROSAT All Sky Survey (RASS). However, we mentioned that the Lupus SFR is located near the UCL subgroup of the Sco-Cen association, so one suspects that some of these stars might not be related to the Lupus SFR but might instead be part of the older Sco-Cen association. Indeed, a list of 81 low-mass candidate members of the UCL subgroup is given in the review by Preibisch & Mamajek (2008), and our sample of Lupus association candidate members includes 25 stars of their list. Preibisch & Mamajek (2008) also claim that other stars previously identified as Lupus PMS stars by Krautter et al. (1997) and Wichmann et al. (1997b) appear to be UCL members because of their proper motions and positions in the HR-diagram.
It appears (Preibisch & Mamajek 2008) that the UCL stars are spread over a large extent on the sky, while the spatial distribution of Lupus stars exhibits two components: (i) an on-cloud population concentrated in the immediate vicinity of the molecular clouds and (ii) a more dispersed off-cloud population surrounding the clouds. The corresponding region spans the range of Galactic coordinates 334° ≤ l ≤ 342° and 5° ≤ b ≤ 25° and contains 160 stars, while the off-cloud population in our sample is located in the region with l < 334°, which contains 57 stars. There is some arbitrariness involved in defining these regions, but the idea of separating our sample into two groups turned out to be necessary to our analysis, as is seen below, and has also been suggested by Preibisch & Mamajek (2008), who wrote “... it is probably wise for astrophysical studies to separate the on- or near-cloud Lupus members from the off-cloud UCL/US members.”. We show that the off-cloud population is indeed a mix of Lupus and UCL stars in Sect. 5.
In Table 2 we present the median positions, proper motions, and radial velocities for the Lupus subgroups. We find good agreement between the proper motions of the various star-forming clouds (Lupus 1–4). When comparing the proper motions of on-cloud and off-cloud stars, we note a difference of about 7 mas/yr in right ascension, although the median values for both populations are perfectly compatible in declination. One possibility to explain this difference is the existence of field stars or UCL members in the off-cloud sample as mentioned before. However, the reported difference of 7 mas/yr is still consistent e.g. with the observed extreme proper motion values for the various subgroups of the Taurus complex where a single UVW velocity value can be adopted assuming a velocity dispersion of about 1 km s-1 among the different subgroups (see Luhman et al. 2009). In the case of radial velocities, we consider both radial velocities from the literature and additional measurements derived in this paper (to be discussed in Sect. 4). Binaries and stars with insignificant measurements (i.e., σVr ≥ Vr) are excluded from the analysis presented in Table 2. We note that stars in the off-cloud region exhibit radial velocities that are slightly higher than the values observed for the on-cloud population, which can be explained by geometrical effects, because both populations are at different angular separations from the CP. In the following, we assume that the Lupus stars spread over the various subgroups are comoving and use the CP search method to identify a moving group of PMS stars in this SFR.
![]() |
Fig. 3 Location of the 217 Lupus PMS stars overlaid on the 12CO intensity map from Tachihara et al. (2001). Stars from US and UCL subgroups of the Sco-Cen association as given in de Zeeuw et al. (1999) and Preibisch & Mamajek (2008) are also included in this figure. Black filled circles denote the Lupus candidate stars and open symbols indicate US and UCL members. Stars marked with green filled squares were classified as both Lupus and UCL members in the literature. The blue dashed box encloses the Lupus star-forming clouds. |
4. Radial velocities
We mentioned that stellar radial velocities (RVs) are needed to determine individual parallaxes. We summarize here our search for RVs in the literature and some additional observations for stars with unknown RVs.
4.1. Radial velocities from the literature
We searched the CDS databases to access RV information for the stars in our sample. The search made use of the data mining tools available on the CDS site. We also looked for published RV data that are not available via the web-based CDS service. Our search for RVs, made as exhaustive as possible, is based on Herbig & Bell (1988), Gregorio-Hetem et al. (1992), Barbier-Brossat et al. (1994), Duflot et al. (1995), Dubath et al. (1996), Grenier et al. (1999), Wichmann et al. (1999), Barbier-Brossat & Figon (2000), Madsen et al. (2002), Melo (2003), Nordström et al. (2004), Bobylev (2006), Gontcharov (2006), James et al. (2006), Malaroda et al. (2006), Torres et al. (2006), Guenther et al. (2007), Kharchenko et al. (2007), and White et al. (2007).
We found RVs for only 108 stars of the full sample (on-cloud and off-cloud populations), which reflects the scarcity of this measurement in the literature. This is because previous spectroscopic investigations in Lupus have often focused on the few bright stars of this region. Figure 4 displays the RV distribution in our sample.
![]() |
Fig. 4 Histogram of RVs found in the literature for Lupus stars. |
4.2. Additional radial velocity observations
The scarcity of measured RVs is the main limitation in deriving individual distances in this work. To increase the number of Lupus stars with known RV information, we thus performed spectroscopic observations with the high-resolution (R = 48 000) FEROS (Kaufer et al. 1999) échelle spectrograph mounted at the ESO/MPG 2.2 m telescope operated at La Silla (Chile). In addition to its high performance, FEROS provides a full wavelength coverage in the optical region (between 3500 Å and 9200 Å) over 39 spectral orders. The observations were taken in object calibration mode, which allows acquiring simultaneous spectra of the object and of the ThAr cell. Exposure times ranged from 5 min to 60 min, so that a S/N of about 20–30 was achieved. We observed 52 stars spread over the Lupus and Ophiuchus1 SFRs during the nights of April 29 to May 05, 2011. Binaries and stars belonging to multiple systems as given in the literature were not included in our list of targets. The observed spectra were reduced with the standard FEROS data reduction pipeline, which performs bias subtraction, flat-fielding, scattered-light removal, échelle-order extraction, barycentric velocity correction, and wavelength calibration of the spectra. The extracted wavelength-calibrated spectra for each spectral order were not merged into a single spectrum, but used as 39 separate spectra (one per order) to determine the RV of the target. This procedure enabled us to eliminate the largest noise contributions that come from the orders in the red (orders ≤7) and blue (orders ≥30) regions of the spectra (see Setiawan et al. 2003).
We derived RVs by cross-correlating the reduced spectra of program stars with template spectra. To maximize the quality of our RV measurements we used both a standard-star spectrum and a numerical mask as template spectrum.
The cross-correlation with a stellar template utilized the standard star HD 82106 (V = 7.2 mag, K3V), which was observed every night and whose RV is known to a high accuracy (Vr = + 29.75 ± 0.05 km s-1, Udry et al. 1999). We computed the RVs of program stars by cross-correlating their spectrograms order-by-order with one reference spectrogram of HD 82106 (taken on May 4, 2011) using the IRAF task fxcor. We obtained the RV for each order separately and then averaged these values. The errors of our RV measurements were determined from the variance of all orders considered. Orders with discrepant values (due to a lower S/N) were obviously not considered. To investigate the accuracy of our results we calculated the RV of HD 82106 as derived from our observations. We cross-correlated the reference spectrogram with all spectrograms of HD 82106 spread over the whole observing campaign. The absolute RV of the reference spectrum was determined by cross-correlating it with the solar spectrum. The mean RV derived from our observations is Vr = + 29.71 ± 0.05 km s-1, in good accordance with Udry et al. (1999).
The other alternative to computing RVs consisted of a cross-correlation with a box-shaped binary template. The procedure follows the method outlined in Baranne et al. (1996), which fits a Gaussian to the cross-correlation function of each order. The center of the Gaussian then gives the RV of the target (see Weise et al. 2010, for more details). To better compare the RVs derived for HD 82106 using both techniques, we used a binary template of a G2V star. Our choice of the template was motivated by the use of the solar spectrum to derive the absolute RV of HD 82106 in our first approach (see above). For program stars, we used a K0V star template that is more representative of our targets’ spectral types. The mean RV derived with this alternative strategy is Vr = + 29.74 ± 0.04 km s-1, which is perfectly consistent with the values mentioned above.
Using both a standard-star spectrum and a numerical template is valuable for gaining confidence in the derived results. To spot possible errors, the RVs of program stars were computed using only the common (not rejected) orders in both procedures. The final RV of our targets is the average of two independent values (one for each method). The uncertainties are calculated by propagating the individual errors and they are of a few hundred m/s. When one of the methods failed to return a RV value (due to low S/N in some orders) we only considered the result derived by the other method. A more realistic idea of our external precision is given in Table 3 by comparing the RVs derived in this work with published results for a control group of 3 stars in our sample with previously known RVs. We conclude that our results are fully compatible with the values found in the literature, with a rms uncertainty of about 500 m/s. We present the RVs of 52 stars belonging to the Lupus and Ophiuchus SFRs in Table 4, together with the Li I and Hα equivalent widths (EWs) derived from our observations.
Comparison of RVs derived in this paper with those published in the literature.
Equivalent widths of the Li I and Hα lines were measured from our spectra using IRAF splot routine. The presence of Li absorption in late-type stars is one of the primary criteria for stellar youth (see Basri et al. 1991). The Li I λ6708 Å resonance doublet is blended by Fe lines and the lithium isotope 6Li. It is not possible to separate individual lines at our spectral resolution (~48 000), so we see no evidence of these features. The lithium EWs were measured using both a Gaussian fit and direct integration. The difference between these values is smaller than 25 mÅ, which we consider to be the upper limit of our measurement errors. The contribution from the neighboring blending lines is expected to be smaller than the uncertainty of our results, which we estimated by varying the location of the continuum adjacent to the line.
The Hα emission line is one the most prominent spectroscopic features in the visible spectra of CTTSs (Joy 1945; Herbig 1962). The line profile is often complex and takes different shapes, which were studied and classified by Reipurth et al. (1996). We measured the Hα EWs by direct integration and by fitting a Voigt profile (when the target exhibited a Hα profile with a single peak). Our measurement errors are mainly caused by the uncertainty on the continuum level in the vicinity of the line. We use the standard limit of 10 Å (see, e.g., Appenzeller & Mundt 1989) to distinguish between CTTSs (EW(Hα) ≥ 10 Å) and WTTSs (EW(Hα) < 10 Å) and find that our sample contains 10 CTTSs and 42 WTTSs. A more detailed inspection of the spectra in our sample revealed two PMS spectroscopic binaries (SB2), for which no evidence of binarity could be found in the literature prior to our observations. We present the results of our observations in Table 4.
5. Analysis and results
In the following we derive the CP of the comoving group of Lupus stars considered in this paper. We first use the k-nearest neighbor (k-NN, Fix & Hodges 1951; Venables & Ripley 2002) algorithm to distinguish between Lupus and UCL stars based on their position, then we apply the CP search method to the sample of Lupus candidate members and perform a membership analysis. The technique that we use to find the CP position of the Lupus moving group is the new CP search method that was recently developed by our team. We refer the reader to the original paper (Galli et al. 2012) for more details on the implementation of this method.
5.1. Preliminary analysis
The velocity dispersion in Lupus is a first input parameter that we must determine for performing the CP analysis. The velocity dispersion of young moving groups is expected to be low, only a few km s-1 (Mathieu 1986). Typical values of velocity dispersion in nearby SFRs are 1–2 km s-1 (Jones & Herbig 1979; Dubath et al. 1996; Makarov 2007; Luhman et al. 2009). We adopt σv = 1 km s-1 as the one-dimensional velocity dispersion of the Lupus moving group, and we come back to this point in Sect. 6. Another input parameter in the CP method is the mean distance to the moving group. As discussed in Sect. 2, the distance to the Lupus SFR has undergone substantial revision in the literature over the past decade. Here we use the value of d = 150 pc that seems to be adequate for most of the clouds in the complex (Comerón 2008). The CP search method is rather insensitive to small variations in both parameters.
When we apply the CP search method to the 217 stars in the sample, we end up with a CP
position given by with
125 moving group members. In a recent paper, Makarov
(2007) investigated the kinematics of the Lupus sample of 93 PMS stars identified
by Krautter et al. (1997) using the UCAC2 catalog
(Zacharias et al. 2004) and derived a velocity
dispersion of 1.3 km s-1. The CP position presented in his work is
(αcp,δcp) = (92.8°, −28.1°) ± (3.1°,5.0°).
We note that there is good agreement between both solutions, and it is probable that the
small differences are caused by the different samples of Lupus stars, the different proper
motion sources, and the different CP search methods. However, whether the above result is
the most appropriate for the sample of TTSs associated with the Lupus molecular clouds is
questionable since the derived CP solution might be affected by UCL stars misidentified as
Lupus WTTSs. We investigate this question below.
5.2. The k-NN algorithm applied to the sample of Lupus stars
The k-NN method is a nonparametric machine learning algorithm used in pattern recognition for classifying objects. Although the k-NN method is not yet widely used in astronomy, it has already proven to be effective, such as for the photometric search of brown dwarfs (Marengo & Sanchez 2009) and estimations of photometric redshifts for quasars (Ball et al. 2007). It is used in this work to segregate potential Lupus and UCL members in our test sample of 217 stars based on their position with respect to both the star-forming clouds and confirmed members of each association.
In the classic version of the k-NN method, a test element is classified by the majority
vote of its k nearest neighbors from the training set.
To begin with, we define the training set
T = { (x1,y1),(x2,y2),...,(xN,yN) },
where xi is the vector of attributes (stellar
positions),
xi = (αi,δi),
and
yi ∈ { c1,c2 }
denote the class membership (Lupus or UCL). The training set used in this work includes
both Lupus and UCL stars, and it is constructed as follows. Since most members of the
Lupus SFR are already included in our test sample (and many of them lack membership
confirmation) we use the position of the 105 molecular clouds detected in 12CO
by Tachihara et al. (2001) for this purpose, and
the 81 stars given by Preibisch & Mamajek
(2008) as UCL members2. Stars in our test
sample are classified as potential Lupus or UCL members based on the majority vote of
their k nearest neighbors where the value of k is
optimized to our specific problem (see below). In this context, the probability
p (y|x) for a given star with
attributes (x,y) is given by
(1)where
kC denotes the number of nearest neighbors
in the training set labeled with class c (Lupus or UCL).
One important point to be considered in the k-NN algorithm is the best choice of k. While higher values of k would make the boundaries between Lupus and UCL less clear, a low value, on the other hand, can lead to noisy classification. An odd number for k is preferable and avoids tied votes since we only have to distinguish between two classes. To gain confidence in the derived results and determine the optimum value for k, we construct a total of 1000 synthetic test samples to be classified by the k-NN method using the procedure described above. To do so, we use the well-known UCL members from the Hipparcos catalog given by de Zeeuw et al. (1999), and simulate synthetic Lupus training stars by generating stellar positions from a Gaussian distribution using the position of the 12CO-peak-integrated intensity and the radius of the molecular clouds as given in Table 1 of Tachihara et al. (2001). At each run our routine randomly chooses Lupus and UCL test stars to construct samples with 200 stars and classify them based on our training set. The fraction of Lupus and UCL stars in the synthetic test samples is not fixed and varies at each iteration to avoid any bias when evaluating the accuracy of the procedure.
The first results of our simulations, however, revealed that many Lupus synthetic stars
located in the core of the star-forming clouds were misclassified. This is because some
data points in our training set refer to single stars, while others in principle represent
the center position of molecular clouds. Neither the size of the molecular clouds nor the
distance to nearer neighbors are taken into account, showing that the majority vote
criteria is not suited to our specific case. One way to solve this problem is to not give
the same weight (vote) to all the neighbors. In this context, we implement a
weighted k-NN algorithm that uses a distance weighted function. We
distinguish whether the positions
(αi,δi)
in the training set refer to molecular clouds (Lupus) or stars (UCL) by assigning a
different weight. The weight function ωi is
given by where
r is the radius of the molecular cloud and d is
calculated from the position of the 12CO-peak-integrated intensity of the Lupus
molecular clouds (see Tachihara et al. 2001) or the
position of the UCL star in the training set. The angular distance d
between a test star with given position
(αj,δj)
and a training star
(αi,δi)
is given by
(2)In this case the
posterior probability for a given star in the test sample takes the weight function into
account and it is defined as
(3)where the sum in the
numerator runs only over the nearest neighbors in the training set with class
c (Lupus or UCL). We use these probabilities to classify the stars in
our test sample as potential Lupus or UCL candidate members.
Figure 5 compares the performance of our refurbished weighted k-NN method with its classic version for various values of k. We conclude that our new strategy indeed yields better results and exhibits an accuracy of ~80%. The accuracy is defined by the fraction of test stars in our simulated test samples that are correctly classified by the k-NN algorithm. We conclude that any value of k between 7 and 15 can be used in this work to segregate potential Lupus and UCL stars. We adopt k = 15 under the assumption that more neighbors will help better define both groups. Doing so, we run the weighted k-NN method in our original test sample of 217 stars and identify 183 stars as potential Lupus members and 34 stars as potential UCL members. The distribution of Lupus and UCL stars is illustrated in Fig. 6 and shows that the majority of Lupus stars is located in the vicinity of the molecular clouds (as expected).
![]() |
Fig. 5 Comparison of performance between the classic and weighted k-NN methods as a function of k derived from our simulations. Each point represents an average value of 1000 iterations. |
5.3. Convergent point analysis for the Lupus moving group
The procedure described in the previous section made it possible to separate some
potential UCL members included in our initial sample of Lupus PMS population. However, in
order to search for association members, proper motions must be considered at this stage.
Our final analysis consists in running the CP search method on the sample of 183 stars
classified by the k-NN method as Lupus stars (see Sect. 5.2). We identify a moving group
with 109 members and CP located at with
chi-squared statistics
(i.e.,
χ2/ν = 118.8/107).
We note that the rejection of some likely UCL members in our sample dramatically changed
the CP position of the moving group as compared to our first solution in Sect. 5.1, which
now appears to be a mixed CP solution of Lupus and UCL stars, and is therefore not valid
for either moving group. Tables 5–7 present the 109 Lupus moving group members selected
by the CP search method, together with their parallaxes (to be discussed in Sect. 6). We
note that six stars that were classified as UCL members by Preibisch & Mamajek (2008) and previously regarded as Lupus stars (Krautter et al. 1997; Wichmann et al. 1997a,b) have been
accepted as Lupus members in our CP analysis. These stars are marked with the symbol “*”.
On the other hand, we found that 19 stars previously regarded as WTTSs following their
discovery in X-ray have been rejected in this analysis, and they are instead probable UCL
members, as discussed by Preibisch & Mamajek
(2008).
![]() |
Fig. 6 Location of the 217 PMS stars overlaid on the 12CO intensity map of Tachihara et al. (2001). Different symbols and colors indicate the membership classification (Lupus or UCL) that results from our k-NN analysis. |
Proper motion and RV for the 19 stars that define the Lupus core moving group.
5.4. Result check via Monte Carlo simulations
To assess the validity of the CP location presented above, we performed Monte Carlo
simulations of the 109 Lupus moving group members. We constructed 1000 samples of moving
groups by resampling the stellar proper motions from a Gaussian distribution where the
mean and variance are equal to the individual stellar proper motions and its uncertainty
(in both components). We ran the CP search method for each set of simulated stars and
computed the CP of each moving group. The results of this study are presented in Fig.
7. The CP solution derived in this paper is fully
consistent with the centroid distribution of the Monte Carlo realizations, located at
and
we therefore conclude that it is representative of the Lupus moving group.
Individual parallax and velocity components for Lupus stars with known RVs.
![]() |
Fig. 7 Realizations of CP coordinates for 1000 Monte Carlo simulations (red dots) overlaid on the X2 contours (solid lines) for the CP solution derived in Sect. 5.3. The blue square denotes the CP coordinates for the Lupus moving group derived in that section, and the green triangle denotes the centroid distribution of simulated CPs (see Sect. 5.4). The dashed lines indicate the 1σ, 2σ, and 3σ contour levels of our CP solution. |
![]() |
Fig. 8 Histograms of the Galactic velocity components for members of the Lupus moving group with known RVs. The red dashed line denotes the average values given in Sect. 6.1. |
![]() |
Fig. 9 Spatial velocities of the 19 group members of the Lupus moving group with known RVs projected on the XY, YZ and ZX planes. |
6. Kinematic parallaxes
6.1. Parallaxes and space velocities of group members with known radial velocity
Once the moving group is defined, it is possible to derive individual kinematic
parallaxes πind for group members if their RVs
Vr are known. Individual parallaxes are given by
(4)where
A = 4.74047 km yr/s is the ratio of one astronomical unit in km to the
number of seconds in one Julian year, λ is the angular distance from the
CP position to a given star in the moving group, and μ∥ the
stellar proper motion component that points towards the CP (see Galli et al. 2012, for more details). The parallax uncertainty is
derived by error propagation of this equation and takes the error budget of proper
motions, RVs, and the CP into account (see Appendix A for more details).
We found RVs for only 60 stars in the sample of 109 moving group members. We reject
(resolved) binaries since it will not be possible to derive their parallaxes from a single
RV measurement. Stars that exhibit poor RVs because of their errors are also excluded from
this analysis. The observed RV for Lupus stars is expected to be low and a small variation
accounts for a more significant shift in the parallax and space velocity (see Appendix B
for more details). To spot possible errors on parallaxes and velocities, we define lower
and upper limits for the space velocity of Lupus stars. To do so, we use five stars3 with known trigonometric parallax in the Hipparcos
catalog that have been selected as moving group members in our CP analysis. The space
velocities derived using Hipparcos parallaxes range from
Vlower ≃ 17 km s-1 and
Vupper ≃ 33 km s-1, which define the limits for
the space velocity of Lupus stars. This leaves us with a sample of 19 stars (see Table
5) that we define here as the Lupus core
moving group. The mean parallax is mas, which agrees well with the
estimated distance of 150 pc assumed in the CP analysis (see Sect. 5.1).
We computed the Galactic velocities for each star using the procedure described in Johnson & Soderblom (1987). Figure 8 displays the distribution of the UVW Galactic velocity components while Fig. 9 shows the velocity vectors of the 19 stars with individual parallaxes in a XYZ grid defined as follows. This reference system has its origin at the Sun where X points to the Galactic center, Y points in the direction of Galactic rotation, and Z points to the Galactic north pole.
The average space velocity for the Lupus moving group derived in this paper is
The
space motion for the UCL subgroup was recently revised by Chen et al. (2011) to take into account the parallaxes from Hipparcos
new reduction (van Leeuwen 2007) and modern RV
compilations. They report a space velocity of
(U,V,W) = (−5.1, −19.7, −4.6) ± (0.6,0.4,0.3)
km s-1. The relative space motion between Lupus stars and the UCL subgroup is
(ΔU,ΔV,ΔW) = (0.0, −1.0, −1.6) ± (0.8,1.2,0.6)
km s-1. These results suggest that Lupus is moving at 1.9 ± 1.6
km s-1 with respect to UCL and that their velocities are statistically
indistinguishable at the 1–2 km s-1 level which roughly corresponds to the
velocity dispersion in each group. A similar conclusion was reached regarding Ophiuchus
stars and the US subgroup of the Sco-Cen association (see Mamajek 2008).
One particular point of the Lupus association is that the observed RVs of young stars are expected to be low and exhibit both positive and negative values (see e.g. Table 4 and Fig. 4). Among the moving group members with individual parallaxes presented in Table 5, a total of six stars exhibit negative values for their RVs. In such cases we use the absolute value of the RVs to compute the individual parallaxes using Eq. (4). In this context, Makarov (2007) reports a significant mismatch between the observed spectroscopic radial velocities and the value inferred from his CP solution, implying a moderate degree of expansion. That the Lupus association of young stars is undergoing expansion with a velocity of ≃1 km s-1 (see Makarov 2007) and the RVs are also near zero possibly explains the existence of the two RV populations. This situation contrasts with other SFRs, e.g. Taurus, where the observed RV is higher, Vr ≃ + 16 km s-1 (see Luhman et al. 2009), and the existence of group members with RVs changing signs (i.e., Vr < 0) cannot be tolerated if assuming a one-dimensional velocity dispersion of only a few km s-1.
When we run the CP search method on the sample of stars with known
Vr < 0, we identify a moving group
with 17 stars and CP located at
(αcp,δcp) = (99.9°, −36.1°) ± (7.3°,9.6°).
The large error bars in this CP solution arise from the small number of stars used for
deriving the CP coordinates (see Galli et al.
2012). As one should notice, the effect of linear expansion changes the CP position
(see also discussion in Makarov 2007), but the CP
mentioned above is still compatible with the one given in Sect. 5.3. In the following we
seek to learn whether the use of either CP solutions lead to significant differences in
the derived parallaxes. We computed the individual parallaxes for those six stars again
with Vr < 0 using the above mentioned
CP and the procedure outlined in Appendix A. The average difference between the recomputed
parallaxes and those given in Table 6 is
Δπ = 0.3 ± 0.1 mas with rms of 1.1 mas, and they are statistically
compatible within their error bars. The re-evaluated UVW velocity components for these
stars has a negligible effect on the average space motion of the Lupus moving group,
yielding instead There
is therefore no reason to reject those stars with
Vr < 0 from our analysis. Owing to
the negligible impact to our results and for clarity of presentation we have decided not
to present the re-evaluated parallaxes and space velocities.
6.2. Approximate parallaxes for other moving group members
The hypothesis that all members of a moving group share the same space motion allows us
to compute tentative parallaxes for group members with unknown RVs. First we derive the
average spatial velocity Vspace from the Galactic velocity of
the stars with known RVs using the 19 stars that define the Lupus core moving group (see
Sect. 6.1). Then we compute an approximate parallax πapp as
(5)The uncertainty
on this tentative parallax is again derived by error propagation and considers in this
case the error budget of proper motions, space velocity, and the CP errors. We present in
Fig. 10 a comparison between individual parallaxes
(see Table 5) and approximate parallaxes for the Lupus core moving group. Both procedures
return similar results within the admittedly large error bars, which tends to justify the
assumption of a common space motion.
![]() |
Fig. 10 Comparison of parallaxes computed with RVs (πind) and the spatial velocity (πapp) for the Lupus core moving group of 19 stars. The red solid line indicates perfect correlation. The mean difference between parallaxes computed with both strategies is 0.2 mas, and the rms is 1.4 mas. |
The velocity dispersion of the cluster prevents all stars from having exactly the same space velocity. However, the procedure described above for deriving approximate parallaxes considers a single value of the space velocity for all stars in the group. We performed Monte Carlo simulations by resampling the spatial velocity from a Gaussian distribution where the mean and variance correspond to the values given in Sect. 6.1 for the average space velocity of the Lupus core moving group. We constructed a total of 1000 realizations. In each run we assigned a different value of space velocity and derived the approximate parallax for each star by using Eq. (5). The average value of the computed parallaxes gives our final result, and the standard error of the mean is propagated into the parallax uncertainty. This strategy not only allows us to reproduce the effect of velocity dispersion, but also makes our parallax results less dependent on a single value used for the spatial velocity of the group. We present in Table 7 the approximate parallaxes derived for the remaining 90 Lupus stars with unknown RVs.
We do stress, however, that the individual approximate parallaxes derived in this way should only be seen as tentative values that will be useful for preliminary statistical analyses of the Lupus association. They should be superseded by more precise values when RV measurements become available for these association members.
6.3. Comparison with Hipparcos parallaxes
As a final check of our results we compare the parallaxes derived in this work with Hipparcos trigonometric parallaxes. We consider both versions of the catalog: the original catalog (ESA 1997, hereafter HIP97), and the new reduction of Hipparcos data (van Leeuwen 2007, hereafter HIP07). We found only six Hipparcos stars among the 109 moving group members and used the group spatial velocity to compute approximate parallaxes. The results of this comparison are presented in Fig. 11. The rms with respect to HIP97 is 2.7 mas and 2.5 mas for HIP07. The mean difference between the parallaxes derived in this work and the ones in HIP97 and HIP07 are –0.9 mas and +0.8 mas, respectively. We conclude that our results are in good agreement with the trigonometric parallaxes given in Hipparcos.
![]() |
Fig. 11 Parallaxes derived in this paper compared with the trigonometric parallaxes given in HIP97 (upper panel) and HIP07 (lower panel). The red solid line represents the expected distribution for equal results. |
6.4. Comparison with Makarov (2007) results
Makarov (2007) estimated kinematic distances for Lupus stars in a similar manner to the procedure described in Sect. 6.2. To do so, he assumed a mean spatial velocity of 22 km s-1 and used the derived distances to infer the depth of the Lupus association. In the following we discuss these findings in light of our own results derived in this study.
The sample of 109 moving group members defined in our CP analysis is cross-correlated with the 93 stars used by Makarov (2007) resulting in 44 stars. When comparing the parallaxes derived in both studies we find a mean difference of Δπ = −0.9 ± 0.2 mas (in the sense “this work” minus “Makarov 2007”) with rms of 1.5 mas (see Fig. 12). In Sect. 6.1 we derived a mean space velocity of Vspace = 22.5 ± 1.1 km s-1 that confirms the value of 22 km s-1 used by Makarov (2007) in his analysis. Since the difference between both values amounts to only 2%, it seems unlikely that the reported systematic effect between the two parallax sets can only be explained by the adopted value of the group spatial velocity. As discussed in Sect. 5.3 the CP position derived in both studies differs significantly by Δλ = 25.4° ± 8.3°.
It is also important to note that the proper motions used in each work play an important role not only for determining the CP position, but also for deriving the parallaxes of group members. When we use proper motions from the UCAC2 catalog as done by Makarov (2007) with the CP derived in this paper, the mean difference between the parallaxes presented in both studies drops to Δπ = −0.5 ± 0.2 mas. This result is illustrated in Fig. 12. However, given the more recent and precise astrometric catalogs available now, such as the ones used in this paper (see Sect. 3.1), there is no good reason to use the UCAC2 catalog where the proper motions for some stars are poorly defined with only two observational data points. That Table 2 of Makarov (2007) does not include the uncertainty of the computed distances makes it difficult to compare the quality of both parallax results. The average error of the proper motions used in this paper and by Makarov (2007) is, respectively, 1.9 mas/yr and 3.3 mas/yr in each component. Given the proper motion and CP errors (see Sects. 5.1 and 5.3) presented in both studies we conclude that the parallaxes derived in this paper are more precise and accurate, because our results take the various sources of errors included in the parallax computation into account.
To investigate the depth of the Lupus moving group we first determine the so-assumed center of the association, (XC,YC,ZC) = (142, −61,38) pc, as defined by the 19 stars in Table 5, and then compute the distance of all 109 group members to the center of the association. Figure 13 displays the histogram of the computed distances. We conclude that the Lupus association shows a large depth of at least 100 pc, exceeding the value of 80 pc reported by Makarov (2007). However, we stress that a more detailed study about the spatial distribution of Lupus stars and the size of the association is necessary when more RVs become available in the future.
![]() |
Fig. 12 Comparison between the parallaxes derived in this work and by Makarov (2007). Parallaxes are calculated using the proper motion catalogs included in this work (see Sect. 3.1) and UCAC2 following Makarov (2007). The red solid line indicates perfect correlation. |
![]() |
Fig. 13 Histogram of the individual distances of the 109 Lupus moving group members to the center of the association. |
![]() |
Fig. 14 Location of the moving group members overlaid on the 12CO intensity map from Tachihara et al. (2001). Different symbols and colors mark the various YSO subclasses. |
7. Discussion
7.1. Positions and parallaxes for pre-main sequence subclasses
The Lupus moving group identified in this paper contains 39 CTTSs, 68 WTTSs, and 2 HAeBes. Among the WTTSs we have 49 on-cloud and 19 off-cloud stars. Figure 14 displays the location of the various PMS subclasses. We note that the Lupus 4 population consists mainly of CTTSs, while in the remaining clouds we find both CTTSs and WTTSs, which may represent a selection effect of our input list of Lupus candidate stars (see Sect. 3). As already predicted, the off-cloud population contains only WTTSs while the CTTSs are located in the immediate vicinity of the molecular clouds.
An interesting point arises when we compare the distances of these various TTS populations. Indeed, we find that the off-cloud WTTSs tend to be closer to us and the CTTSs more distant. We present in Table 8 the average parallaxes (distances) for each subgroup in the Lupus SFR.
Properties of the TTSs in Lupus.
7.2. Notes on the parallaxes of Lupus subgroups
The parallaxes derived in this paper allow us to investigate the properties of the various subgroups in this cloud complex. Figure 15 displays the location of Lupus members and corresponding clouds. The histogram of parallaxes is presented in Fig. 16. In the following we discuss the parallaxes derived for these subgroups based on the results presented in Table 9. Using the information provided in Table 2 we calculate the mean UVW velocities for each subgroup and confirm that they exhibit coherent motions with a one-dimensional velocity dispersion of ~1 km s-1 among the various subgroups (see Table 9).
Properties of the various subgroups in Lupus.
![]() |
Fig. 15 Location of the moving group members overlaid on the 12CO intensity map from Tachihara et al. (2001). Different symbols and colors mark the various clouds of the Lupus complex. |
Full sample:
the average error on parallaxes is
mas, yielding an average precision of about 25% in our distance results derived in this
paper. We note the existence of background and foreground populations (as anticipated in
Sect. 3), confirming that the Lupus complex occupies a large volume in space.
Off-cloud population:
the WTTSs that form this population are scattered not only in angular extent but also in depth (see Fig. 16). Different scenarios have been suggested to explain the dispersed population of WTTSs: the stars may have been formed in the vicinity of the clouds and have reached their present location because of the velocity dispersion in the group (Wichmann et al. 1997a); they might have been ejected with high velocities by dynamical interactions in multiple systems (Sterzik & Durisen 1995); star formation could have occurred in small cloudlets that have now dispersed (Feigelson 1996). One important parameter in this discussion is the stellar age that can now be accurately determined for those stars with individual parallaxes. Age determination goes beyond the scope of the present work and will be discussed in a forthcoming paper.
Lupus 1, 2, and 4:
all together they contain 37% of moving group members. While Lupus 1 and 2 are somewhat closer, Lupus 4 is more distant (see Table 9). We derived individual parallaxes for four stars in Lupus 1 and one star in Lupus 2. Our results for Lupus 4 are less accurate, because we only computed tentative parallaxes, and our sample is biased towards CTTSs as mentioned before. However, Hip 78092 is an HAeBe star projected in the direction of Lupus 4, and Hipparcos results for this star, πHIP97 = 5.04 ± 1.18 mas and πHIP07 = 4.29 ± 0.98 mas, confirm the approximate distance to the cloud derived in this paper.
Lupus 3:
we find significant depth effects in Lupus 3, while the stars in the remaining clouds
are less dispersed along the line of sight (see, e.g., Fig. 16). One possibility for explaining this result is the existence of
various components along the line of sight that lie at different distances. We computed
the average parallax of the CTTSs (26 stars) and WTTSs (23 stars) in Lupus 3 and find
The
depth of Lupus 3 derived from the closest and remotest parts of this cloud is in good
agreement with Lombardi et al. (2008, see
Sect. 2). This preliminary result must be confirmed by further investigations
since we derived accurate parallaxes for only a few stars.
![]() |
Fig. 16 Histogram of parallaxes for the various subgroups of the Lupus complex. |
7.3. Expansion of the Lupus association
In Sect. 6.1 we mentioned that the Lupus association is undergoing expansion (see also Makarov 2007). From the proper motions alone, it is impossible to distinguish between the state of linear expansion and parallel space motions, so RVs are needed to distinguish between them. In the following we use the Blaauw’s expansion model (Blaauw 1964), combined with the radial velocities and individual parallaxes given in Tables 5 and 6, to investigate whether the expansion rate of the association can be detected.
In the Blaauw model for linear expansion the stellar RVs are given by (6)where
κ is the expansion term, K a systematic error in the
radial velocities, and d the individual distance of the star. A linear
expansion exists when κ > 0 in this equation. To
derive κ, we plot the difference between the observed spectroscopic RVs
and the predicted RVs (
) inferred from the CP coordinates
as a function of the individual distances. Then, we solve for the slope κ
in Eq. (6) (see Fig. 17). The best fit to the data yields κ = 0.021 ± 0.004
km s-1 pc-1. We reject GQ Lup, RX J1608.5-3847, and Sz121 from
this analysis, because of the large errors in distance (see Table 6), and RX J1613.0-4004 because of its large deviation, presumably due
to a poor RV. Although the derived expansion term κ is small when
compared to, say, the TW Hydrae association (see Mamajek
2005), it is positive within 3σ of the computed uncertainty.
This result confirms that the linear expansion in the Lupus association of young stars is
real but small.
![]() |
Fig. 17 Blaauw’s linear expansion model applied to the Lupus association of young stars. The red dashed line indicates the weighted least-squares fit to the data. |
8. Conclusions
We have identified a moving group of 109 stars in the Lupus SFR by applying our new CP search method combined with the k-NN algorithm, which made it possible to perform a membership analysis and to distinguish between Lupus and UCL stars. We applied the CP search method to different subsets of association members and confirmed our solution with Monte Carlo simulations. Because of the close proximity to the UCL subgroup and the difficulties encountered in separating these groups, we claim that we have detected a minimum moving group in our CP analysis that may not contain all stars kinematically associated to Lupus. We derived accurate parallaxes for members with known RVs that we define as the Lupus core moving group and used the group spatial velocity to tentatively calculate approximate parallaxes for the remaining stars.
Determination of individual parallaxes is restricted to moving group members with known RVs. The RVs of many stars in the group either have never been measured or are of low quality so cannot be used to derive accurate parallaxes. We presented new RV measurements for 52 PMS stars of our starting sample of Lupus candidate stars based on spectroscopic observations performed with FEROS. We encourage observers to employ the available precise spectrographs to perform spectroscopic surveys of this complex SFR. Additional observations will also allow one to detect new binaries and confirm the youth of many PMS candidates.
We show that the CTTSs are located close to or within the molecular clouds, while the WTTSs are dispersed not only in angular extent but also in depth. We find evidence of depth effects in Lupus 3 that must be confirmed by further investigations. Based on the individual distances and radial velocities derived in this paper, we also confirm that the Lupus association is undergoing expansion. These results represent a first step toward better understanding of the structure of the Lupus cloud complex. The individual distances derived in this paper will be used in a forthcoming paper to investigate the physical properties of Lupus stars.
Online material
Results from our observations.
Approximate parallax derived in this paper for Lupus stars with unknown RVs.
Appendix A: Error propagation of parallaxes
The parallax uncertainty for those stars with known RVs is obtained by error
propagation of Eq. (4) and it is given by
(A.1)While
σμ∥ and
σVr refer to errors on
observational quantities (i.e., proper motions and RVs), the use of
σλ in Eq. (A.1) is not straightforward. Galli et al. (2012) performed extensive simulations
that convincingly demonstrate that the precision of the CP position is influenced by
several parameters. The increasing velocity dispersion of the cluster and the existence
of possible interlopers in the moving group also play an important role when evaluating
the precision and accuracy of the CP position. These sources of errors can be roughly
divided into two parts: (i) observational errors that arise mainly from proper motion
errors (since errors on stellar position can be neglected to a first-order
approximation); and (ii) geometric effects (e.g., angular distance from the CP to the
moving group, cluster concentration, distance, and number of group members). To
correctly introduce the CP errors in the parallax uncertainty, one needs to separate the
contribution of these effects, otherwise the error budget due to proper motions will be
considered twice in Eq. (A.1). To do so,
we set the proper motion errors to zero (i.e.,
σμα,δ = 0).
Then, we use Eq. (19) of de Bruijne (1999) to
estimate the velocity dispersion of the Lupus moving group that results from the
peculiar motion of the stars given by the μ⊥ statistics.
This procedure yields
km s-1 and confirms the adopted value of 1 km s-1 as the
one-dimensional velocity dispersion of the moving group used in our CP analysis (see
Sect. 5.1). By adopting σv = 0.7
km s-1 and
σμα,δ = 0,
the uncertainties on the CP position decrease to
(σαcp,σδcp) = (1.6°,1.2°) yielding σλ ≃ 2.0°. This result represents a more realistic estimate of the CP error budget due to geometric effects, not including proper motions errors, to be used in the error propagation of parallaxes.
Appendix B: Notes on radial velocities
One important point is that the observed RVs for Lupus stars are expected to be lower
than to other SFRs, such as Taurus-Auriga and Chamaeleon (see Bertout & Genova 2006; James
et al. 2006). The average value and standard deviation of the RVs presented in
Table 4 for Lupus stars is
km s-1, and the average error is
σVr = 0.5 km s-1.
Here we investigate the errors on parallaxes and space velocities caused by
uncertainties in RVs. A small variation ΔVr in RVs accounts
for the variation Δπ in parallaxes that can be approximated by
(B.1)This shifts the space
velocity of the star by ΔVspace that is given as
(B.2)Using the values of
position, proper motion, and radial velocities given in Table 2 for the on-cloud population, and the CP solution derived in Sect.
5.3, we note that a small shift ΔVr = 0.5 km s-1
in the measured RVs accounts for Δπ ≃ 1.4 mas, and consequently,
ΔVspace ≃ 4 km s-1. We thus emphasize that
high-precision RVs are needed to derive reliable kinematic parallaxes of Lupus stars.
Although a larger sample of UCL members from the Hipparcos catalog exists (see e.g. de Zeeuw et al. 1999), we prefer at this stage to use only those provided by Preibisch & Mamajek (2008) in order to keep approximately the same fraction of Lupus and UCL data points in the training set to avoid a classification bias. The remaining UCL stars will be used later in this section to investigate the accuracy of our procedure and determine the optimum value of k (see text).
Sz 120 is not considered in this analysis because it is a binary HAeBe star (Correia et al. 2006).
Acknowledgments
It is a pleasure to thank Johny Setiawan and Maren Mohler for their assistance with the FEROS Data Reduction Software and cross-correlation techniques used for radial velocity computations. We thank Kengo Tachihara for providing the image for Figs. 3, 6, 14, and 15, and the PASJ editorial office for granting permission to reproduce this material. We also acknowledge the referee for the careful report. This research was funded by a doctoral/post-doctoral fellowship from FAPESP and made use of the SIMBAD database operated at the CDS, Strasbourg, France.
References
- Appenzeller, I., & Mundt, R. 1989, A&ARv, 1, 291 [NASA ADS] [CrossRef] [Google Scholar]
- Ball, N. M., Brunner, R. J., Myers, A. D., et al. 2007, ApJ, 663, 774 [NASA ADS] [CrossRef] [Google Scholar]
- Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barbier-Brossat, M., & Figon, P. 2000, A&AS, 142, 217 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barbier-Brossat, M., Petit, M., & Figon, P. 1994, A&AS, 108, 603 [NASA ADS] [Google Scholar]
- Basri, G., Martin, E. L., & Bertout, C. 1991, A&A, 252, 625 [NASA ADS] [Google Scholar]
- Bertout, C. 1987, in Circumstellar Matter, eds. I. Appenzeller, & C. Jordan, IAU Symp., 122, 23 [Google Scholar]
- Bertout, C., & Genova, F. 2006, A&A, 460, 499 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bertout, C., Basri, G., & Bouvier, J. 1988, ApJ, 330, 350 [NASA ADS] [CrossRef] [Google Scholar]
- Bertout, C., Robichon, N., & Arenou, F. 1999, A&A, 352, 574 [NASA ADS] [Google Scholar]
- Blaauw, A. 1964, ARA&A, 2, 213 [NASA ADS] [CrossRef] [Google Scholar]
- Bobylev, V. V. 2006, Astron. Lett., 32, 816 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, C. H., Mamajek, E. E., Bitner, M. A., et al. 2011, ApJ, 738, 122 [NASA ADS] [CrossRef] [Google Scholar]
- Comerón, F. 2008, in The Lupus Clouds, ed. B., Reipurth, 295 [Google Scholar]
- Correia, S., Zinnecker, H., Ratzka, T., & Sterzik, M. F. 2006, A&A, 459, 909 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Bruijne, J. H. J. 1999, MNRAS, 310, 585 [NASA ADS] [CrossRef] [Google Scholar]
- de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, A. G. A., & Blaauw, A. 1999, AJ, 117, 354 [NASA ADS] [CrossRef] [Google Scholar]
- Dubath, P., Reipurth, B., & Mayor, M. 1996, A&A, 308, 107 [NASA ADS] [Google Scholar]
- Ducourant, C., Teixeira, R., Périé, J. P., et al. 2005, A&A, 438, 769 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Duflot, M., Figon, P., & Meyssonnier, N. 1995, A&AS, 114, 269 [NASA ADS] [Google Scholar]
- Dzib, S., Loinard, L., Mioduszewski, A. J., et al. 2010, ApJ, 718, 610 [NASA ADS] [CrossRef] [Google Scholar]
- Dzib, S., Loinard, L., Rodríguez, L. F., Mioduszewski, A. J., & Torres, R. M. 2011, ApJ, 733, 71 [NASA ADS] [CrossRef] [Google Scholar]
- ESA 1997, VizieR Online Data Catalog: I/239 [Google Scholar]
- Feigelson, E. D. 1996, ApJ, 468, 306 [NASA ADS] [CrossRef] [Google Scholar]
- Fix, E., & Hodges, J. 1951, Discriminatory analysis. Nonparametric discrimination: Consistency properties, Technical Report 4, Project Number 21-49-004, USAF School of Aviation Medicine, Randolph Field, Texas [Google Scholar]
- Galli, P. A. B., Teixeira, R., Ducourant, C., Bertout, C., & Benevides-Soares, P. 2012, A&A, 538, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ghez, A. M., McCarthy, D. W., Patience, J. L., & Beck, T. L. 1997, ApJ, 481, 378 [NASA ADS] [CrossRef] [Google Scholar]
- Girard, T. M., van Altena, W. F., Zacharias, N., et al. 2011, AJ, 142, 15 [NASA ADS] [CrossRef] [Google Scholar]
- Gontcharov, G. A. 2006, Astron. Lett., 32, 759 [NASA ADS] [CrossRef] [Google Scholar]
- Gregorio-Hetem, J., Lepine, J. R. D., Quast, G. R., Torres, C. A. O., & de La Reza, R. 1992, AJ, 103, 549 [NASA ADS] [CrossRef] [Google Scholar]
- Grenier, S., Burnage, R., Faraggiana, R., et al. 1999, A&AS, 135, 503 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Guenther, E. W., Esposito, M., Mundt, R., et al. 2007, A&A, 467, 1147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Herbig, G. H. 1962, Adv. Astron. Astrophys., 1, 47 [NASA ADS] [CrossRef] [Google Scholar]
- Herbig, G. H., & Bell, K. R. 1988, in Third Catalog of Emission-Line Stars of the Orion Population, 3 [Google Scholar]
- Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [NASA ADS] [Google Scholar]
- Hughes, J., Hartigan, P., & Clampitt, L. 1993, AJ, 105, 571 [NASA ADS] [CrossRef] [Google Scholar]
- James, D. J., Melo, C., Santos, N. C., & Bouvier, J. 2006, A&A, 446, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [NASA ADS] [CrossRef] [Google Scholar]
- Jones, B. F., & Herbig, G. H. 1979, AJ, 84, 1872 [NASA ADS] [CrossRef] [Google Scholar]
- Joy, A. H. 1945, ApJ, 102, 168 [NASA ADS] [CrossRef] [Google Scholar]
- Kaufer, A., Stahl, O., Tubbesing, S., et al. 1999, The Messenger, 95, 8 [Google Scholar]
- Kharchenko, N. V., Scholz, R.-D., Piskunov, A. E., Röser, S., & Schilbach, E. 2007, Astron. Nachr., 328, 889 [NASA ADS] [CrossRef] [Google Scholar]
- Knude, J., & Hog, E. 1998, A&A, 338, 897 [NASA ADS] [Google Scholar]
- Knude, J., & Nielsen, A. S. 2001, A&A, 373, 714 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Krautter, J., Wichmann, R., Schmitt, J. H. M. M., et al. 1997, A&AS, 123, 329 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Loinard, L., Torres, R. M., Mioduszewski, A. J., et al. 2007, ApJ, 671, 546 [NASA ADS] [CrossRef] [Google Scholar]
- Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 480, 785 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- López Martí, B., Jiménez-Esteban, F., & Solano, E. 2011, A&A, 529, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Luhman, K. L., Mamajek, E. E., Allen, P. R., & Cruz, K. L. 2009, ApJ, 703, 399 [NASA ADS] [CrossRef] [Google Scholar]
- Madsen, S., Dravins, D., & Lindegren, L. 2002, A&A, 381, 446 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Makarov, V. V. 2007, ApJ, 658, 480 [NASA ADS] [CrossRef] [Google Scholar]
- Malaroda, S., Levato, H., & Galliani, S. 2006, VizieR Online Data Catalog: III/249 [Google Scholar]
- Mamajek, E. E. 2005, ApJ, 634, 1385 [NASA ADS] [CrossRef] [Google Scholar]
- Mamajek, E. E. 2008, Astron. Nachr., 329, 10 [Google Scholar]
- Mamajek, E. E., Meyer, M. R., & Liebert, J. 2002, AJ, 124, 1670 [NASA ADS] [CrossRef] [Google Scholar]
- Marengo, M., & Sanchez, M. C. 2009, AJ, 138, 63 [NASA ADS] [CrossRef] [Google Scholar]
- Mathieu, R. D. 1986, Highlights of Astronomy, 7, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Melo, C. H. F. 2003, A&A, 410, 269 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Merín, B., Jørgensen, J., Spezzi, L., et al. 2008, ApJS, 177, 551 [NASA ADS] [CrossRef] [Google Scholar]
- Nordström, B., Mayor, M., Andersen, J., et al. 2004, A&A, 418, 989 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Preibisch, T., & Mamajek, E. 2008, in The Nearest OB Association: Scorpius-Centaurus (Sco OB2), ed. B. Reipurth, 235 [Google Scholar]
- Racine, R. 1968, AJ, 73, 233 [CrossRef] [Google Scholar]
- Reipurth, B., Pedrosa, A., & Lago, M. T. V. T. 1996, A&AS, 120, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roeser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440 [NASA ADS] [CrossRef] [Google Scholar]
- Schwartz, R. D. 1977, ApJS, 35, 161 [Google Scholar]
- Setiawan, J., Pasquini, L., da Silva, L., von der Lühe, O., & Hatzes, A. 2003, A&A, 397, 1151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sterzik, M. F., & Durisen, R. H. 1995, A&A, 304, L9 [NASA ADS] [Google Scholar]
- Tachihara, K., Toyoda, S., Onishi, T., et al. 2001, PASJ, 53, 1081 [NASA ADS] [Google Scholar]
- Torres, C. A. O., Quast, G. R., da Silva, L., et al. 2006, A&A, 460, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Torres, R. M., Loinard, L., Mioduszewski, A. J., & Rodríguez, L. F. 2007, ApJ, 671, 1813 [NASA ADS] [CrossRef] [Google Scholar]
- Torres, R. M., Loinard, L., Mioduszewski, A. J., & Rodríguez, L. F. 2009, ApJ, 698, 242 [NASA ADS] [CrossRef] [Google Scholar]
- Torres, R. M., Loinard, L., Mioduszewski, A. J., et al. 2012, ApJ, 747, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Udry, S., Mayor, M., & Queloz, D. 1999, in IAU Colloq. 170, Precise Stellar Radial Velocities, eds. J. B. Hearnshaw, & C. D. Scarfe, ASP Conf. Ser., 185, 367 [Google Scholar]
- van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Venables, W. N., & Ripley, B. D. 2002, Modern Applied Statistics, ed. S. Fourth (Springer Verlag) [Google Scholar]
- Walter, F. M. 1986, ApJ, 306, 573 [NASA ADS] [CrossRef] [Google Scholar]
- Weise, P. 2010, Ph.D. Thesis, University of Heidelberg [Google Scholar]
- Weise, P., Launhardt, R., Setiawan, J., & Henning, T. 2010, A&A, 517, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- White, R. J., Gabor, J. M., & Hillenbrand, L. A. 2007, AJ, 133, 2524 [NASA ADS] [CrossRef] [Google Scholar]
- Wichmann, R., Krautter, J., Covino, E., et al. 1997a, A&A, 320, 185 [NASA ADS] [Google Scholar]
- Wichmann, R., Sterzik, M., Krautter, J., Metanomski, A., & Voges, W. 1997b, A&A, 326, 211 [NASA ADS] [Google Scholar]
- Wichmann, R., Covino, E., Alcalá, J. M., et al. 1999, MNRAS, 307, 909 [NASA ADS] [CrossRef] [Google Scholar]
- Zacharias, N., Urban, S. E., Zacharias, M. I., et al. 2004, AJ, 127, 3043 [NASA ADS] [CrossRef] [Google Scholar]
- Zacharias, N., Finch, C. T., Girard, T. M., et al. 2012, VizieR Online Data Catalog: I/322 [Google Scholar]
All Tables
Median positions, proper motions, radial velocities, and the number of stars for the various subgroups in Lupus.
All Figures
![]() |
Fig. 1 Proper motion and associated errors for the 217 stars considered in this work. |
In the text |
![]() |
Fig. 2 Proper motion vectors of the 217 Lupus candidate stars. |
In the text |
![]() |
Fig. 3 Location of the 217 Lupus PMS stars overlaid on the 12CO intensity map from Tachihara et al. (2001). Stars from US and UCL subgroups of the Sco-Cen association as given in de Zeeuw et al. (1999) and Preibisch & Mamajek (2008) are also included in this figure. Black filled circles denote the Lupus candidate stars and open symbols indicate US and UCL members. Stars marked with green filled squares were classified as both Lupus and UCL members in the literature. The blue dashed box encloses the Lupus star-forming clouds. |
In the text |
![]() |
Fig. 4 Histogram of RVs found in the literature for Lupus stars. |
In the text |
![]() |
Fig. 5 Comparison of performance between the classic and weighted k-NN methods as a function of k derived from our simulations. Each point represents an average value of 1000 iterations. |
In the text |
![]() |
Fig. 6 Location of the 217 PMS stars overlaid on the 12CO intensity map of Tachihara et al. (2001). Different symbols and colors indicate the membership classification (Lupus or UCL) that results from our k-NN analysis. |
In the text |
![]() |
Fig. 7 Realizations of CP coordinates for 1000 Monte Carlo simulations (red dots) overlaid on the X2 contours (solid lines) for the CP solution derived in Sect. 5.3. The blue square denotes the CP coordinates for the Lupus moving group derived in that section, and the green triangle denotes the centroid distribution of simulated CPs (see Sect. 5.4). The dashed lines indicate the 1σ, 2σ, and 3σ contour levels of our CP solution. |
In the text |
![]() |
Fig. 8 Histograms of the Galactic velocity components for members of the Lupus moving group with known RVs. The red dashed line denotes the average values given in Sect. 6.1. |
In the text |
![]() |
Fig. 9 Spatial velocities of the 19 group members of the Lupus moving group with known RVs projected on the XY, YZ and ZX planes. |
In the text |
![]() |
Fig. 10 Comparison of parallaxes computed with RVs (πind) and the spatial velocity (πapp) for the Lupus core moving group of 19 stars. The red solid line indicates perfect correlation. The mean difference between parallaxes computed with both strategies is 0.2 mas, and the rms is 1.4 mas. |
In the text |
![]() |
Fig. 11 Parallaxes derived in this paper compared with the trigonometric parallaxes given in HIP97 (upper panel) and HIP07 (lower panel). The red solid line represents the expected distribution for equal results. |
In the text |
![]() |
Fig. 12 Comparison between the parallaxes derived in this work and by Makarov (2007). Parallaxes are calculated using the proper motion catalogs included in this work (see Sect. 3.1) and UCAC2 following Makarov (2007). The red solid line indicates perfect correlation. |
In the text |
![]() |
Fig. 13 Histogram of the individual distances of the 109 Lupus moving group members to the center of the association. |
In the text |
![]() |
Fig. 14 Location of the moving group members overlaid on the 12CO intensity map from Tachihara et al. (2001). Different symbols and colors mark the various YSO subclasses. |
In the text |
![]() |
Fig. 15 Location of the moving group members overlaid on the 12CO intensity map from Tachihara et al. (2001). Different symbols and colors mark the various clouds of the Lupus complex. |
In the text |
![]() |
Fig. 16 Histogram of parallaxes for the various subgroups of the Lupus complex. |
In the text |
![]() |
Fig. 17 Blaauw’s linear expansion model applied to the Lupus association of young stars. The red dashed line indicates the weighted least-squares fit to the data. |
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.