Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

A publishing partnership

The following article is Open access

The IGRINS YSO Survey. III. Stellar Parameters of Pre-main-sequence Stars in Ophiuchus and Upper Scorpius

, , , , , , , , , , , and

Published 2023 January 24 • © 2023. The Author(s). Published by the American Astronomical Society.
, , Citation Ricardo López-Valdivia et al 2023 ApJ 943 49 DOI 10.3847/1538-4357/acab04

Download Article PDF
DownloadArticle ePub

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

As featured in:
Youtube
opens in new tab

Article metrics

868 Total downloads

Share this article

Author e-mails

rlopezv@astro.unam.mx

Author affiliations

1 The University of Texas at Austin, Department of Astronomy, 2515 Speedway, Stop C1400, Austin, TX 78712-1205, USA; rlopezv@astro.unam.mx

2 Universidad Nacional Autónoma de México, Instituto de Astronomía, AP 106, Ensenada 22800, BC, México

3 Lowell Observatory, 1400 W. Mars Hill Road, Flagstaff. AZ 86001, USA

4 Physics & Astronomy Department, Rice University, 6100 Main Street, Houston, TX 77005, USA

5 Korea Astronomy and Space Science Institute, 776 Daedeok-daero, Yuseong-gu, Daejeon 34055, Republic of Korea

ORCID iDs

Dates

  1. Received 2022 September 2
  2. Revised 2022 December 8
  3. Accepted 2022 December 9
  4. Published

Check for updates using Crossmark

Unified Astronomy Thesaurus concepts

Fundamental parameters of stars; Infrared sources; Pre-main sequence stars; High resolution spectroscopy

Journal RSS
Create or edit your corridor alerts

What are corridors?opens in new tab

0004-637X/943/1/49

Abstract

We used the Immersion GRating Infrared Spectrometer (IGRINS) to determine fundamental parameters for 61 K- and M-type young stellar objects (YSOs) located in the Ophiuchus and Upper Scorpius star-forming regions. We employed synthetic spectra and a Markov chain Monte Carlo approach to fit specific K-band spectral regions and determine the photospheric temperature (T), surface gravity ($\mathrm{log}g$), magnetic field strength (B), projected rotational velocity ($v\sin i$), and K-band veiling (rK). We determined B for ∼46% of our sample. Stellar parameters were compared to the results from Taurus-Auriga and the TW Hydrae association presented in Paper I of this series. We classified all the YSOs in the IGRINS survey with infrared spectral indices from Two Micron All Sky Survey and Wide-field Infrared Survey Explorer photometry between 2 and 24 μm. We found that Class II YSOs typically have lower $\mathrm{log}g$ and $v\sin i$, similar B, and higher K-band veiling than their Class III counterparts. Additionally, we determined the stellar parameters for a sample of K and M field stars also observed with IGRINS. We have identified intrinsic similarities and differences at different evolutionary stages with our homogeneous determination of stellar parameters in the IGRINS YSO survey. Considering $\mathrm{log}g$ as a proxy for age, we found that the Ophiuchus and Taurus samples have a similar age. We also find that Upper Scorpius and TWA YSOs have similar ages, and are more evolved than Ophiuchus/Taurus YSOs.

Export citation and abstract BibTeX RIS

Previous article in issue
Next article in issue

1. Introduction

Young stellar objects (YSOs) are stars at an early stage of evolution and their properties have a dominant impact on the environments of planet development. Stellar parameters like temperature (T) and surface gravity ($\mathrm{log}g$) permit comparisons with evolutionary models (Baraffe et al. 2015; Simon et al. 2019). With these fundamental parameters, the differences between models and observations for stars with masses <1.5 M⊙ and age <10 Myr (e.g., Simon et al. 2013) can be better understood. The accurate and precise determination of YSO parameters is limited by the impacts of interstellar reddening, continuum veiling, and magnetic fields. These processes, if not considered, can result in inaccurate YSO parameters and thus lead to wrong conclusions about star formation and evolution.

For example, the veiling is a nonstellar continuum emission that reduces the depth of the photospheric lines in YSO spectra (Joy 1949). Such reduction of the line depth could lead to erroneous T or $\mathrm{log}g$ values as these parameters have similar effects on certain atomic lines. On the other hand, the presence of strong magnetic fields (B) changes the absorption line profiles through Zeeman broadening (e.g., Johns-Krull et al. 1999; Yang et al. 2005; Lavail et al. 2017; Sokal et al. 2018; Lavail et al. 2019; Sokal et al. 2020).

The Immersion GRating INfrared Spectrometer (IGRINS; Yuk et al. 2010; Park et al. 2014; Mace et al. 2016) was designed to minimize the problems associated with determining YSO parameters. IGRINS employs a silicon immersion grating as the primary disperser (Jaffe et al. 1998; Gully-Santiago et al. 2010; Wang et al. 2010) and volume phase holographic gratings to cross disperse the H- and K-band echellograms onto Teledyne Hawaii-2RG arrays. This setup provides a compact design with high sensitivity and a significant single-exposure spectral grasp at a spectral resolution of R ∼ 45,000 (Yuk et al. 2010; Park et al. 2014). IGRINS has a fixed spectral format and no moving optics, so the science products are consistent over time baselines of several years. IGRINS has increased its scientific value by traveling among the McDonald Observatory, the Lowell Discovery Telescope (LDT), and the Gemini South telescope (Mace et al. 2018).

The IGRINS YSO survey is uniquely suited to determining physical parameters because (i) the simultaneous coverage of the H and K bands (1.45–2.5 μm) of IGRINS includes separable photospheric and disk contributions to the YSO spectrum. (ii) The fixed spectral format of IGRINS provides similar spectral products for each object at each epoch. (iii) Multiple-epoch observations can provide a means to characterize or average over variability. (iv) The velocity resolution of ∼7 km s−1 (R ∼ 45,000) is smaller than the typical young star rotational velocity and can resolve magnetic fields ≳1 kG.

In López-Valdivia et al. (2021, hereafter Paper I), we presented the first results of the IGRINS YSO survey. This first study consisted of the determination of photospheric temperature (T), surface gravity ($\mathrm{log}g$), magnetic field (B), projected rotational velocity ($v\sin i$), and the K-band veiling (rK ) for 110 YSOs located in the Taurus-Auriga star-forming region and 19 young stars in the TW Hydrae association (TWA). The IGRINS YSO survey analysis continued with the determination of low-resolution veiling spectra for 144 Taurus-Auriga members in Kidder et al. (2021, hereafter Paper II). In this paper, the third of the series, we extend our analysis to K- and M-type YSOs in the Ophiuchus and Upper Scorpius star-forming regions (Elias 1978; Lada & Wilking 1984; Wilking et al. 1989). Combining the results of Paper I with those obtained here improved our understanding of YSO stellar parameters as a function of age.

2. Observations and Sample

Among the IGRINS YSO survey, we included bright objects (K < 11 mag) with spectral types between K0 and M5 (Wilking et al. 2005; Torres et al. 2006; Ricci et al. 2010; Esplin & Luhman 2020; Luhman & Esplin 2020), which were classified as members of Ophiuchus and Upper Scorpius by Hsieh & Lai (2013), Rebull et al. (2018), and Esplin & Luhman (2020).

We observed our sample with IGRINS at the McDonald Observatory 2.7 m telescope, the LDT, and the Gemini South Telescope between 2014 and 2019. We followed the observing and data reduction methods as outlined in Paper I. In brief, we observed the YSOs and A0V telluric standards 6 by nodding the stars along the slit in patterns made up of AB or BA pairs, where A and B are two different positions on the slit. We aimed to obtain a minimum signal-to-noise ratio (S/N) of 70 in the K band by adjusting the total exposure time based on the IGRINS exposure time calculator. 7 An S/N above 70 is sufficient to use our methods and to produce reliable stellar parameters.

We employed the IGRINS pipeline (Lee et al. 2017) 8 to reduce all the spectroscopic data. The pipeline produces a telluric corrected spectrum with a wavelength solution derived from OH night sky emission lines and telluric absorption lines. The wavelength solution was corrected for the barycenter velocity determined with zbarycorr (Wright & Eastman 2014).

About a third of the YSOs in our sample have more than one epoch available. For the determination of average parameters (this work) a weighted-average spectrum is produced from multi-epoch observations using the S/N at each data point as the weight. The standard deviation of the mean gives the final uncertainties per data point. Future studies of the IGRINS YSO survey will look at multi-epoch variability.

We identified binaries (or multiples) in the YSO sample with a separation of less than 2''. When the separation is greater than 2'', the objects can be observed individually by IGRINS. If the separation is less than 2'', the resultant IGRINS spectrum will contain some flux from both components. A complete determination of the binary nature of these stars is beyond the scope of this work. Binary searches for the YSOs in our sample have identified some binaries and provided limits on detection methods (Simon et al. 1995; Haisch et al. 2002; Prato et al. 2003; Ratzka et al. 2005; Prato 2007; Kraus et al. 2008; Duchêne 2010; Cieza et al. 2010; Barenfeld et al. 2019). We marked known binaries and left as singles those objects that appear as limits in binary classification studies or meet our separation criterion. Those objects that lack information in the literature have a question mark. We have excluded double-lined binaries from our sample through careful inspection of the YSO survey spectra. The binary status and basic information of our sample are available in Table 1.

Table 1. Basic Information and Results for Our Ophiuchus and Upper Scorpius Sample

2MASSName N K SpTrefclusterRefBinRef AV ${\alpha }_{{K}_{s}-{\rm{W}}4}$ Flag T $\mathrm{log}g$ B $v\sin i$ rK
   mag      mag  KdexkGkm s−1  
J16081474-1908327EPIC 20515254828.4K218USco14Y160.6 ± 0.1−2.79 ± 0.0614652 ± 1724.55 ± 0.341.05 ± 0.708.4 ± 3.80.31 ± 0.15
J16082324-1930009EPIC 20508061619.5K918USco14N80.8 ± 0.5−1.09 ± 0.2213782 ± 2774.16 ± 0.442.03 ± 0.9114.3 ± 4.10.15 ± 0.11
J16090075-1908526EPIC 20515138729.2M118USco14N80.7 ± 0.5−0.70 ± 0.2703636 ± 1884.22 ± 0.321.88 ± 0.458.8 ± 2.90.26 ± 0.09
J16093030-2104589RX J1609.5-210518.9M018USco15N160.6 ± 0.6−2.75 ± 0.0404021 ± 2783.98 ± 0.471.44 ± 0.9010.7 ± 3.90.20 ± 0.11
J16110890-1904468ScoPMS 4447.7K418USco15N161.4 ± 0.2−2.78 ± 0.0714220 ± 1604.50 ± 0.33<2.6627.1 ± 3.40.44 ± 0.09
J16113134-1838259V* V866 Sco105.8K518USco14Y33.4 ± 0.5−0.49 ± 0.0513919 ± 2303.49 ± 0.39<1.5216.5 ± 3.67.45 ± 1.00
J16142029-1906481EPIC 20515823917.8M018USco14N132.8 ± 0.8−0.71 ± 0.0713911 ± 3914.13 ± 0.60<1.9724.5 ± 8.31.64 ± 0.53
J16153456-2242421V* VV Sco17.9M018USco15Y81.2 ± 0.7−0.92 ± 0.2213685 ± 3564.17 ± 0.562.47 ± 1.0510.2 ± 5.91.16 ± 0.31
J16211848-2254578EPIC 204290918110.2M218USco14N112.1 ± 0.7−0.84 ± 0.1603414 ± 2294.09 ± 0.45<1.5119.4 ± 3.40.67 ± 0.16
J16220961-1953005EPIC 20500067618.9M3.717USco14Y161.6 ± 0.7−2.08 ± 0.0913371 ± 2283.80 ± 0.40<1.2416.1 ± 2.90.19 ± 0.11
J16232454-1717270EPIC 20548325839.7M2.518USco15?⋯⋯⋯3⋯⋯⋯⋯⋯

Note. The "Flag" columns listing 0, 1, 2, and 3 indicate the good, acceptable, at the edge of the grid, and bad determinations of atmospheric parameters. The K magnitude is from Cutri et al. (2003). The full version of this table is available in the online version of this paper.

References. (1) Simon et al. (1995) (2) Haisch et al. (2002) (3) Prato et al. (2003) (4) Ratzka et al. (2005) (5) Wilking et al. (2005) (6) Torres et al. (2006) (7) Prato (2007) (8) Kraus et al. (2008) (9) Ricci et al. (2010) (10) Duchêne (2010) (11) Cieza et al. (2010) (12) Hsieh & Lai (2013) (13) Lafreniére et al. (2014) (14) Rebull et al. (2018) (15) Luhman et al. (2018) (16) Barenfeld et al. (2019) (17) Esplin & Luhman (2020) (18) Luhman & Esplin (2020).

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

3. Stellar Parameter Determination

We followed the same methods as in Paper I to obtain T, $\mathrm{log}g$, B, $v\sin i$, and rK . Our method used a Markov chain Monte Carlo (MCMC) analysis as implemented in the code emcee (Foreman-Mackey et al. 2013). We computed a four-dimensional (T, $\mathrm{log}g$, B, and $v\sin i$) grid of synthetic spectra using the moogstokes code (Deen 2013).

moogstokes synthesizes the emergent spectrum of a star taking into account the Zeeman broadening produced by the presence of a photospheric magnetic field. In this work, we used the MARCS atmospheric models (Gustafsson et al. 2008) with solar metallicity (suitable for YSOs; D'Orazi et al. 2011) and the astrophysical-inferred modifications to the Vienna Atomic Line Database (Ryabchikova et al. 2015) line transition data ($\mathrm{log}{gf}$ and van der Waals constant values) presented by Flores et al. (2019). We used a microturbulence of 1 km s−1 as low-mass stars have microturbulence values between 0 and 2 km s−1 (e.g., Gray 2005; Bean et al. 2006)

Our grid of synthetic spectra matched the IGRINS spectral resolution (R ∼ 45,000) and covers the parameter space as follows: from 3000–5000 K in T (steps of 100 K up to 4000 K, and 250 K above 4000 K), from 3.0–5.0 dex in $\mathrm{log}g$ (steps of 0.5 dex), from 0–4 kG in B (steps of 0.5 kG), and values from 2–50 km s−1 in $v\sin i$ (steps of 2 km s−1). The grid steps 9 are spaced enough to differentiate the effects of the parameters on the spectra.

The synthetic spectral grid covers four spectral intervals in the K band, namely, Na, Ti, Ca, and CO. The intervals include neutral atomic (Fe, Ca, Na, Al, and Ti) and molecular lines (CO) that are sensitive to changes in the different stellar parameters (see Figure 2 of Paper I). The four spectral intervals are shown in Figure 1. Each spectral region was normalized using an interactive Python script. The script fits a polynomial of order n to a custom number of flux bins, usually between 10 and 20 bins, and a polynomial of order 1–4. After normalizing the continuum, we carried out the MCMC trial comparing observed and synthetic spectra in the K-band spectral regions by allowing T, $\mathrm{log}g$, B, $v\sin i$, and veiling to vary along with small continuum (<6%) and wavelength (<1.0 Å) offsets. In each MCMC trial, we linearly interpolated within the four-dimensional (T, $\mathrm{log}g$, B, and $v\sin i$) synthetic spectral grid to obtain the corresponding spectrum with the sampled set of parameters. The interpolated synthetic spectrum was then artificially veiled by the veiling parameter and re-normalized for each region. A single veiling value was used for all the K-band wavelength regions.

Figure 1.

Figure 1. The spectral regions used in this work for the star DoAr 24 (solid black line, 2MASS J16261706-2420216). We included the best-fit synthetic spectrum (red-dashed line), which has T = 4092 K, $\mathrm{log}g$ = 4.11 dex, B = 2.12 kG, $v\sin i$ = 10.6 km s−1, and a K-band veiling of 0.79. The bottom panels show the residuals between the observed and the synthetic spectra. The level of agreement between the best-fit synthetic spectrum and the DoAr 24 observed spectrum was considered a good determination. We also indicate some atomic and molecular lines present in the spectral regions.

Download figure:

Standard image High-resolution image

Finally, from the posterior probability distributions of the MCMC, we took the 50th percentile as the most likely value for each stellar parameter. We computed the total uncertainty for each parameter as the quadrature sum of the formal fit errors (the larger of the 16th and 84th percentiles of the posterior probability distribution) and the systematic errors. The systematic errors were determined in Paper I and we assume the same values for this work. They are 75 K, 0.13 dex, 0.26 kG, and 1.7 km s−1 for T, $\mathrm{log}g$, B, and $v\sin i$, respectively.

4. Results

Once we determined the stellar parameters of our sample (Table 1), we categorized our determinations, based on a visual inspection and χ2 statistics, with a numerical quality flag equal to 0, 1, or 3 if they are good, acceptable, or poor determinations. We identified objects with lower/upper limits outside or at the edge of our grid with a flag of 2. We used each stellar parameter's total uncertainty to compute its lower and upper limits. In Figure 1, we show an illustrative example of how the synthetic spectra reproduce the observations.

We identified 31, 25, 5, and 9 stars whose parameters are good, acceptable, at the edge of the grid or poor determinations, respectively. We excluded poor determinations in our further analysis and referred to the 61 remaining stars as our sample.

4.1. Photospheric Temperature

The derived temperatures for nonbinary Ophiuchus and Upper Scorpius YSOs follow the same trend with spectral type as was found for single Taurus YSOs in Paper I. To improve our previous IGRINS temperature scale, we combined our previous T values of Taurus and TWA (Paper I) with the T determinations made in this work. The combined sample (Paper I and this work) comprises 190 YSOs.

We computed the mean T and the standard deviation of the mean in bins of ±0.5 spectral type subclasses to construct the IGRINS temperature scale. We included in our IGRINS temperature scale spectral type bins with more than one object.

In Figure 2, we compare the updated IGRINS temperature scale, which is in Table 2, to the published scales of Pecaut & Mamajek (2013), Herczeg & Hillenbrand (2014) and Luhman et al. (2003). In general, the shape of our temperature scale is similar to the published ones. The T values for M0–M3 stars agree with the three temperature scales. However, our scale agrees better with the Luhman et al. (2003) for late spectral types. For the K stars, we found that our temperature values are cooler than the temperature scales found in Pecaut & Mamajek (2013) and Herczeg & Hillenbrand (2014).

Figure 2.

Figure 2. Temperature as a function of the literature spectral type. The dots represent the IGRINS temperature scale, determined using the temperatures obtained for the nonbinary YSOs of Ophiuchus, Upper Scorpius, TWA, and Taurus. The different lines represent the temperature scales of Herczeg & Hillenbrand (2014), Luhman et al. (2003), and Pecaut & Mamajek (2013). Our temperature scale agrees with the published temperature scales for the M0–M3 stars but follows the Luhman et al. (2003) scale for late types. However, our temperatures are cooler than the published scales in K stars. The error bar in our temperature is the standard deviation of the mean within the ±0.5 spectral type subclasses.

Download figure:

Standard image High-resolution image

Table 2. Temperature Scale Determined Using IGRINS Data Along with Those of Luhman et al. (2003, L03), (Pecaut & Mamajek (2013, PM13), and (Herczeg & Hillenbrand (2014, HH14)

SpTL03PM13HH14IGRINSN
 (K)(K)(K)(K) 
K0⋯50304870⋯⋯
K1⋯4920⋯⋯⋯
K2⋯476047104424 ± 182
K3⋯4550⋯⋯⋯
K4⋯4330⋯4158 ± 432
K5⋯414042104025 ± 495
K6⋯4020⋯4016 ± 3014
K7⋯397040203861 ± 374
K8⋯3940⋯3815 ± 295
K9⋯3880⋯3795 ± 92
M0⋯377039003856 ± 4315
M13705363037203663 ± 3316
M23560349035603525 ± 2016
M33415336034103462 ± 1819
M43270316031903332 ± 2916
M53125288029803237 ± 2611

Note. This IGRINS temperature scale is an update of that presented in Paper I. The last column lists the number of objects for which we computed the mean T value and the standard error of the mean within ±0.5 spectral type subclasses. We only included spectral type bins with more than one determination.

Download table as:  ASCIITypeset image

The differences between our T values and those in the literature could be an effect of having a small sample of K stars. Also, it can be explained as a cumulative result of using different atmospheric models and line lists, different methodologies (spectroscopy versus photometry), different wavelength intervals (optical versus infrared), the effect of the B field on T, and even differing SpT classifications.

For example, typical SpT uncertainties are about one subclass (∼100 K); however, in many cases, the SpT determined by different groups for a certain YSO disagrees by more than two subclasses, even more than five in extreme cases. Additionally, omitting the B field in the T determination produces temperatures, on average, ∼40–70 K cooler than if the B is considered (Paper I).

One way that we reduce SpT uncertainties in this work is to use homogeneous SpT classifications. The SpT of the Taurus sample comes from Luhman et al. (2017), while for Ophiuchus and Upper Scorpius members, the SpT comes mostly from Esplin & Luhman (2020) and Luhman & Esplin (2020). These three studies measured SpT by comparing absorption bands (TiO, VO, Na i, K i, and H2O) between observed and standard optical and infrared spectra. Their reported uncertainties in SpT are ± 0.25 and ± 0.5 subclasses for optical and infrared types.

4.2. Surface Gravity

The $\mathrm{log}g$ distributions of Ophiuchus and Upper Scorpius are differentiated in Figure 3. Their Kolmogorov–Smirnov (K-S) probability 10 of 0.1% also shows that they are statistically different. The Ophiuchus distribution has a mean $\mathrm{log}g$ (and error on the mean) of 3.83 ± 0.03 dex, while Upper Scorpius has 4.13 ± 0.08 dex.

Figure 3.

Figure 3. Probability density of $\mathrm{log}g$ for 48, 13, 110, and 19, YSOs members of Ophiuchus, Taurus, Upper Scorpius, and TWA. We also show the $\mathrm{log}g$ determinations obtained for 133 K and M field stars.

Download figure:

Standard image High-resolution image

Considering the $\mathrm{log}g$ as a proxy for age, we can speculate that such differences in $\mathrm{log}g$ between samples might be due to age. To test this, we used the combined sample, and we also included the $\mathrm{log}g$ values determined for 133 field stars (see Appendix A.2) as a more evolved counterpart.

From Figure 3, we can see that the $\mathrm{log}g$ distributions of Taurus and Ophiuchus are very similar as its K-S probability of 78% also confirms. We also found that the $\mathrm{log}g$ distribution of Upper Scorpius is comparable to the TWA distribution (K-S probability of 35%).

We can also see in Figure 3 that the $\mathrm{log}g$ distribution of field stars is distinguishable from the YSOs, as is expected for a more evolved sample. There is some overlapping between the low $\mathrm{log}g$ values of the field stars, which is about 4.0 dex, and the higher $\mathrm{log}g$ determinations of Taurus, Ophiuchus, Upper Scorpius, and TWA.

Based on our $\mathrm{log}g$ determinations and using them as an age proxy, we can conclude that Ophiuchus and Taurus have a similar age, which is younger than the indistinguishable Upper Scorpius and TWA samples.

In Figure 4, we produce the Kiel (T versus $\mathrm{log}g$) diagram for the YSOs in Taurus, Ophiuchus, Upper Scorpius, and TWA. We also included in Figure 4 the magnetic evolutionary models of Feiden (2016) as a theoretical comparison. The Taurus and Ophiuchus objects seem to populate similar parts of the Kiel diagram, and there is no clear difference in $\mathrm{log}g$ between them. These objects also have a larger $\mathrm{log}g$ dispersion than their counterparts in Upper Scorpius or TWA.

Figure 4.

Figure 4. Spectroscopic Hertzprung–Russell or Kiel diagram for 48 Ophiuchus (circles), 110 Taurus-Auriga (plus), 13 Upper Scorpius (stars), and 19 TWA (squares) YSOs. The Taurus and Ophiuchus samples seem to populate similar parts of this diagram, and it is no clear difference in $\mathrm{log}g$. Some objects of Upper Scorpius and Taurus have a $\mathrm{log}g$ value relatively high compared with most of our determinations. We included magnetic evolutionary tracks of Feiden (2016). Finally, the error bar in the right-bottom corner represents the median uncertainties of our combined sample.

Download figure:

Standard image High-resolution image

4.3. Spectral Indices and Classification

Based on the shape of their de-reddened spectral energy distributions (SEDs), Lada & Wilking (1984) divided the YSOs into three different morphological classes (I, II, and III), by means of the spectral index $\alpha =d\mathrm{log}(\lambda {F}_{\lambda })/d\mathrm{log}\lambda $.

The class scheme represents an evolutionary sequence of low-mass stars, and it is helpful to understand the star formation process. Although most of our YSOs have a previous classification, we computed de-reddened α indices for our YSOs to investigate trends with stellar parameters.

We first constructed the infrared SED from the Two Micron All Sky Survey (2MASS; J, H, and Ks ) and Wide-field Infrared Survey Explorer (WISE; W1, W2, W3, and W4) magnitudes reported in the ALLWISE catalog (Cutri et al. 2021). We included those magnitudes with a flux signal-to-noise ratio greater than 2 and ignored the upper limits.

We then corrected the observed SEDs by reddening. We used a visual extinction (AV) estimation that we converted into a extinction in the Ks band (${A}_{{K}_{{\rm{s}}}}$) by means of ${A}_{{K}_{{\rm{s}}}}$/AV = 0.112 (Rieke & Lebofsky 1985; Großschedl et al. 2019). We then employed the zero magnitude flux densities, and extinction laws reported in Indebetouw et al. (2005) and Großschedl et al. (2019) to compute the de-reddened J, H, K, W1, W2, W3, and W4 fluxes. Using the de-reddened fluxes and the central wavelengths of each filter, we created the extinction-corrected infrared SED for those YSOs in our combined sample with AV estimations. We obtained AV values with the code MassAge (J. Hernández et al. 2023, in preparation). MassAge input includes T, Gaia EDR3 (Gp, Rp, and Bp), and 2MASS (J and H) photometry to compute AV by minimizing the differences between the observed and expected intrinsic colors of Luhman & Esplin (2020) affected by reddening. The reddening AV is changed until the best comparison is found using the minimal χ2 method.

The uncertainties in AV values are obtained using the Monte Carlo method of error propagation (Anderson 1976), assuming Gaussian distributions for the uncertainties in the input parameters. We found photometric data for only 171 objects, for which we computed the AV (see Table 1).

Continuum veiling (rK ) and AV affect the stellar light similarly, i.e., they redden it. Since we determine T by taking into account rk , we ensure that massage is correctly determining the AV values. Moreover, if we compare the rk values as a function of their MassAge AV we do not find any trend. Furthermore, these data have a correlation coefficient of 0.18, suggesting a weak correlation between our rk and AV values.

The presence, structure, accretion status, and evolutionary state of a circumstellar disk impact the photometric colors, producing an incorrect AV value. As some of our YSOs might still possess a disk, we ran MassAge a second time. This time we did not include the Gaia Bp magnitude, which should be the photometric band most affected by the accretion, to assess the impact of accretion in our sample. We found that the AV obtained using all the photometric bands are higher, on average 0.14 mag, than those values determined ignoring Bp magnitude. This difference is lower than the mean uncertainty determined for AV, which is about 0.4 mag.

Finally, to obtain the infrared spectral index (${\alpha }_{{K}_{s}-{\rm{W}}4}$) we perform a linear fit between the Ks and W4 de-reddened fluxes. Our infrared indices are reported in Tables 1 and 3.

As a quality check, in Figure 5, we compared our spectral indices with those determined by Dunham et al. (2015) and Luhman et al. (2010) for YSOs in the Ophiuchus and Taurus star-forming regions, respectively. Dunham et al. (2015) determined the spectral index α through a linear least-squares fit of all available 2MASS and Spitzer photometry between 2 and 24 μm. On the other hand, Luhman et al. (2010) computed the spectral index between four pairs of photometric bands, including the pair Ks and Spitzer 24 μm. These spectral indices are compatible with those determined here, as we obtained them from the same wavelength range (2–24 μm). The only difference between the works of Dunham et al. (2015), Luhman et al. (2010), and ours is that we used WISE photometry instead of Spitzer.

Figure 5.

Figure 5. Comparison of the de-reddened α indices determined in this work with those determined in Dunham et al. (2015) and Luhman et al. (2010). The gray squares represent the area where they agree on the class status. We used the classification thresholds of Großschedl et al. (2019) for Class III, Class II, Flat spectrum, and Class 0+I, respectively. The error bar on our determinations is the uncertainty in the best-fit slope. Our indices are in good agreement for Taurus, but they present a systematic offset toward higher values (on average, +0.25) for Ophichus YSOs. The reason for this offset is unclear but might be related to using different photometric data, zero-points, and extinction laws.

Download figure:

Standard image High-resolution image

Our classifications are consistent with the literature. For Taurus, we found that ${\alpha }_{{K}_{s}-{\rm{W}}4}$ values are around the one-to-one line for the entire interval. On the other hand, our Ophiuchus α indices are higher than those determined by Dunham et al. (2015), on average, ∼0.25. The source of this discrepancy is unclear but is likely related to the difference between the WISE and Spitzer photometry and the use of different photometric zero-points and extinction laws.

4.4. Spectral Indices versus Stellar Parameters

Once we homogeneously computed ${\alpha }_{{K}_{s}-{\rm{W}}4}$ values, we looked for trends between them and the determined stellar parameters. As we mentioned before, the α index informs us of the evolutionary status of the YSO.

We did not find any correlation between ${\alpha }_{{K}_{s}-{\rm{W}}4}$ and T. The dispersion of the T values in the four populations is very similar, and also between Class II and Class III YSOs. 11 Regarding $\mathrm{log}g$, we found that, on average, $\mathrm{log}g$ values of Class III are higher than their Class II counterparts. We found a mean and standard deviation of the mean for Class II YSOs of 3.86 ± 0.02 dex, while for Class III we found 4.02 ± 0.04 dex.

Similarly, we found that $v\sin i$ values of Class III are slightly higher than those of Class II. We found a mean value of 15.0 ± 0.7 km s−1 and 18.2 ± 1.4 km s−1 for the Class II and Class III YSOs, respectively. This result can be related to the disk locking scenario in which the presence of an accretion disk regulates the stellar rotation (Serna et al. 2021).

We consider 65 B determinations from Taurus, Ophiuchus, Upper Scorpius, and TWA. We did not find any trend with the ${\alpha }_{{K}_{s}-{\rm{W}}4}$ values and the B field. Both Class II and Class III YSOs have B fields very similar between them. We found a mean value of 1.96 ± 0.06 kG for Class II and 2.04 ± 0.15 kG for Class III YSOs, respectively.

The strongest and expected trend of ${\alpha }_{{K}_{s}-{\rm{W}}4}$ is seen for the K-band veiling values (see Figure 6). More evolved objects have lower K-band veiling values (rK ≲0.6). Interestingly, the K-band veiling values show a large dispersion for Class II YSOs. Such a distribution could be related to the variable nature of the Class II YSOs or intrinsic differences between the YSOs.

Figure 6.

Figure 6. Comparison between the K-band veiling and T (top), $\mathrm{log}g$ (middle), and ${\alpha }_{{K}_{s}-{\rm{W}}4}$ (bottom) for the Ophiuchus (black circles), Upper Scorpius (green stars), Taurus (red crosses), and TWA (blue squares). The veiling values of Taurus, Ophiuchus, and Upper Scorpius are more spread out than those of TWA. We also included in each panel the mean error bar.

Download figure:

Standard image High-resolution image

We also compared T and $\mathrm{log}g$ with the K-band veiling in Figure 6. We can see that the TWA objects have a K-band veiling less than ∼0.8. The veiling values of Taurus, Ophiuchus, and Upper Scorpius are more spread out than those of TWA. Most are between 0 and 3.2, while a few YSOs have K-band veiling greater than 3.2 and up to ∼7.3.

Although of low significance, there seems to be a trend toward higher rK as the $\mathrm{log}g$ decreases. The latter suggestion implies that less evolved objects have high veiling values., which aligns with the idea that the veiling is related to the presence of a circumstellar disk. Also, this behavior might be due in part to a degeneracy between rK and $\mathrm{log}g$, as their effects on the spectral lines are similar.

In terms of T, between ∼3600 and ∼4100 K, the rK could be as high as ∼7.3. This differs from other temperature ranges, where the highest K-band veiling is around 2.

4.5. Magnetic Field

The ability to determine the B field from the Zeeman broadening is related directly to the spectral resolution and the $v\sin i$ of the star. Based on the results of Hussaini et al. (2020), the minimum B field we can detect from IGRINS spectra should be greater than 1.0 kG, and this scales with the $v\sin i$ as B ≥ $v\sin i$/8 (kG/km s−1). We consider B values that do not meet this criterion to be a non-detection. We found that seven and 20 B-field values meet the detection threshold in the Upper Scorpius and Ophiuchus star-forming regions, respectively. The Kolmogorov–Smirnov probability computed between the two samples of B-field determinations is 50%, which means that both samples are not statistically different at that level. In Table 1 we report both determinations and limits for the B field of our targets.

The combined sample (Paper I and this work) has 78 B-field determinations (20 in Ophiuchus, seven in Upper Scorpius, 41 in Taurus, and 10 in TWA) and 112 B-field limits. We have compared stellar parameters for the YSOs with B-field detections and non-detections looking for possible trends. In terms of T, we found that the mean T for the B-field determinations and B-field limits is very similar, 3675 and 3655 K, respectively. This result suggests that the two groups are indistinguishable in terms of T, as it is also confirmed by their K-S test probability of 35%. In terms of K-band veiling, we found a K-S probability of 5% with a mean K-band veiling for the B-field determinations of 0.88 and the B-field non-detections of 0.92. We then focused on the B-field determinations and split our 78 determinations into two B-field strength bins: low B field (1 ≤ B <2.0 kG) and high B field (B ≥2.0 kG). According to the K-S probability, the low and high B-field distributions are statistically the same for T and K-band veiling at 40% and 14%, respectively.

Interestingly, we found a mean $\mathrm{log}g$ of 3.88 and 4.08 dex for the low B- and high B-field bins, respectively. We also found a K-S probability of about 1%, meaning that both strength bins are different in terms of $\mathrm{log}g$. This is the expected behavior if most of the B field comes from primordial field conservation during radius contraction of pre-main-sequence stars (Moss 2003).

4.6. Temperature Bins

The homogeneous stellar parameter determination that we perform for the YSOs in Ophiuchus, Upper Scorpius, Taurus, and TWA, along with the field stars, allowed us to analyze the similarities and differences of these samples from an evolutionary point of view.

As we see in the previous sections, there are differences between the YSO populations in their $\mathrm{log}g$ distributions, which we attribute to differences in age. However, the stellar mass plays a crucial parameter in the evolution of the stars. T is a rough proxy for stellar mass in most evolutionary models, as constant mass tracks are primarily vertical in the Kiel diagram (e.g., Figure 4).

For that reason, and to fully understand if the differences found in $\mathrm{log}g$ are due to age, we separated our star-forming regions into three different T bins. Our YSO T values span a range of ∼1500 K, between ∼3100 and ∼4600 K, so our bins comprise stars within 500 K. We name our bins as low (3100 ≤ T < 3600 K), med (3600 ≤ T < 4100 K), and high (T ≥ 4100 K) temperature, and they contain 82, 92, and 16 objects, respectively. In Figure 7, we show how the $\mathrm{log}g$, $v\sin i$, and K-band veiling behave in each mass bin and each star-forming region. We also include field stars in the $\mathrm{log}g$ and $v\sin i$ analysis.

Figure 7.

Figure 7. Mean $\mathrm{log}g$, $v\sin i$, and K-band veiling as a function of T bins and by star-forming region. The YSOs of Taurus, Ophiuchus, Upper Scorpius, and TWA are divided into low (circles, 3100 ≤T < 3600 K), med (squares, 3600 ≤T < 4100 K), and high (triangles, T ≥ 4100 K) temperature bins. The error bar is the standard error on the mean ($\sigma /\sqrt{n}$) of each T bin. We also included in the left and middle panels the field stars.

Download figure:

Standard image High-resolution image

In the four star-forming regions, the behavior with $\mathrm{log}g$ (left panel Figure 7) is the same. Objects in the low and med T bins have lower mean $\mathrm{log}g$ values than those in the high T bin. However, in the field stars, the high T bin has low mean $\mathrm{log}g$ values. According to the evolutionary tracks of Marigo et al. (2017), for an age of ∼500 Myr, stars with T between 3000 and 4000 K (roughly our low and med T bins) have $\mathrm{log}g$ values ranging from ∼5.0 to 4.7 dex, while stars with K spectral types (high T bin) have $\mathrm{log}g$ in the range of ∼4.66–4.63 dex. Our $\mathrm{log}g$ values agree with the theoretical predictions of Marigo et al. (2017) for the M stars. However, our results for the K stars are lower by ∼0.3 dex than such predictions. Considering that $\mathrm{log}g$ is a challenging parameter to determine, we need a more detailed analysis of the field star sample to address the source of this $\mathrm{log}g$ discrepancy in the K stars, which is beyond the scope of this paper. Our $\mathrm{log}g$ values follow the trend of the theoretical predictions where K stars have higher $\mathrm{log}g$ than M stars, and they are helpful in the relative sense we have used them throughout the paper. We see from the middle panel of Figure 7 that the mean $v\sin i$ values create a descending ladder of the mean $v\sin i$ values in the following order: Taurus, Ophiuchus, Upper Scorpius, TWA, and field stars. Finally, the mean values of the K-band veiling (right panel Figure 7) are very similar between T bins and between star-forming regions.

5. Summary and Conclusions

Using high signal-to-noise K-band spectra obtained with IGRINS and an MCMC approach to spectral fitting, we determined stellar parameters for a sample of 61 YSOs in the Ophiuchus and Upper Scorpius complexes. We also determined the stellar parameters for some K and M field stars reported in López-Valdivia et al. (2019). We simultaneously fit four observed spectral regions with a synthetic spectral grid computed with solar metallicity MARCS models and the spectral synthesis code moogstokes to obtain T, $\mathrm{log}g$, $v\sin i$, rK , and B-field strength. Our B-field determination relies on the Zeeman broadening, whose detection is limited by the $v\sin i$ of the star and the spectral resolution. For this reason, we determined B just for 27 YSOs and provided limits on the B field for the remaining stars in the sample.

We found that our temperatures for M0–M3 stars agree with three published temperature scales (Luhman et al. 2003; Pecaut & Mamajek 2013; Herczeg & Hillenbrand 2014). For late M stars, our T values agree more with those of Luhman et al. (2003) while for K stars, we found that our temperatures are colder than the temperature scales of Pecaut & Mamajek (2013) and Herczeg & Hillenbrand (2014). The differences in the T can be explained as a cumulative result of small number statistics, differences in the employed methodology, use of different atmospheric models, line lists, and wavelength intervals, the effect of the B field on T, and inconsistent spectral type classifications.

We found that the mean $\mathrm{log}g$ value of Ophiuchus (3.83 ± 0.03 dex) is similar to that of Taurus (3.87 ± 0.03 dex), while the mean $\mathrm{log}g$ of Upper Scorpius (4.13 ± 0.08 dex) is similar to the value found for TWA (4.22 ± 0.03 dex). If we consider $\mathrm{log}g$ as a proxy for age, our results suggest that the Ophiuchus and Taurus samples have a similar age. We also find that Upper Scorpius and TWA YSOs have similar ages, and are more evolved than Ophiuchus/Taurus YSOs. We derived visual extinction and gathered from literature 2MASS and WISE photometry to construct de-reddened infrared SEDs for our sample. We compute spectral indices and classify our YSO samples into morphological classes with these SEDs. In the combined sample of YSOs, which includes objects from Ophiuchus, Taurus, Upper Scorpius, and TWA, we found that Class II YSOs have, on average, lower $\mathrm{log}g$ and $v\sin i$, similar B field, and higher K-veiling values than their Class III counterparts.

We appreciate the detailed comments from the anonymous referee, which improved the science and text of this manuscript. This work used IGRINS, which was developed under a collaboration between the University of Texas at Austin and the Korea Astronomy and Space Science Institute (KASI) with the financial support of the US National Science Foundation under grant Nos. AST-1229522 and AST-1702267, of the University of Texas at Austin, and of the Korean GMT Project of KASI. R.L.V. acknowledges support from CONACYT through a postdoctoral fellowship within the program Estancias Posdoctorales por México. This paper includes data taken at The McDonald Observatory of The University of Texas at Austin. These results made use of the Lowell Discovery Telescope. Lowell is a private, nonprofit institution dedicated to astrophysical research and public appreciation of astronomy and operates the LDT in partnership with Boston University, the University of Maryland, the University of Toledo, Northern Arizona University, and Yale University. Based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). This material is based upon work supported by the National Science Foundation under grant No. AST-1908892 to G.M.

Software: IGRINS pipeline package (Lee et al. 2017), MoogStokes (Deen 2013), MARCS models (Gustafsson et al. 2008), zbarycorr (Wright & Eastman 2014), MassAge (J. Hernández et al. 2023, in preparation).

Appendix: Stellar Parameters of Taurus, TWA, and Field Stars

A.1. Taurus and TWA

We included in this work the Taurus and TWA YSOs analyzed in Paper I. We compare their stellar parameters to establish similarities and differences between these four YSO populations.

In Table 3, we gathered the stellar parameters reported in Paper I for the Taurus and TWA YSOs. We also added the values of AV and ${\alpha }_{{K}_{s}-{\rm{W}}4}$ computed in this work for many of them.

Table 3. Stellar Parameters for Taurus and TWA YSOs

2MASSSpTCluster AV ${\alpha }_{{K}_{s}-{\rm{W}}4}$ Flag T $\mathrm{log}g$ B $v\sin i$ rK
   mag   K dexkGkm s−1  
J04043984+2158215M3.2Tau0.7 ± 0.6−2.6 ± 0.013461 ± 2164.54 ± 0.352.75 ± 0.6017.7 ± 3.40.09 ± 0.06
J04053087+2151106M2.7Tau0.6 ± 0.6−2.6 ± 0.013496 ± 1834.61 ± 0.243.21 ± 0.395.0 ± 2.40.04 ± 0.03
J04131414+2819108M3.6Tau0.4 ± 0.4−2.8 ± 0.103393 ± 1543.40 ± 0.24<1.4331.2 ± 2.30.23 ± 0.06
J04132722+2816247M0.5Tau2.6 ± 0.3−2.8 ± 0.113875 ± 1813.95 ± 0.34<2.4927.3 ± 4.00.38 ± 0.11
J04141358+2812492M4.5Tau0.9 ± 0.5−0.6 ± 0.203318 ± 2204.19 ± 0.422.67 ± 0.609.7 ± 3.21.62 ± 0.21
J04141458+2827580M3.5Tau0.8 ± 0.5−0.4 ± 0.203251 ± 1573.34 ± 0.271.14 ± 0.334.6 ± 2.30.37 ± 0.08
J04141760+2806096M5.0Tau1.9 ± 0.6−0.3 ± 0.103281 ± 1963.88 ± 0.40<1.099.9 ± 2.91.10 ± 0.21
J04143054+2805147M2.2Tau5.1 ± 0.60.2 ± 0.103620 ± 1863.64 ± 0.34<2.3622.7 ± 3.90.73 ± 0.12
J04144730+2646264M2.6Tau0.7 ± 0.8−1.2 ± 0.203530 ± 2573.72 ± 0.42<1.7231.0 ± 4.20.44 ± 0.15
J04144786+2648110M2.5Tau0.5 ± 0.5−0.7 ± 0.203520 ± 1803.61 ± 0.30<1.3920.5 ± 2.50.28 ± 0.08
J04144928+2812305M3.5Tau1.6 ± 0.3−0.9 ± 0.103246 ± 1043.47 ± 0.21<1.0311.1 ± 1.90.52 ± 0.06
J04154278+2909597M1.0Tau2.2 ± 0.3−1.4 ± 0.503723 ± 1253.87 ± 0.221.36 ± 0.328.2 ± 2.00.12 ± 0.04
J04162810+2807358M2.0Tau1.1 ± 0.3−2.7 ± 0.013713 ± 1383.85 ± 0.21<2.2627.3 ± 2.30.23 ± 0.04
J04173372+2820468M2.3Tau0.1 ± 0.5−1.3 ± 0.103445 ± 1553.73 ± 0.28<1.3110.8 ± 2.30.81 ± 0.13
J04173893+2833005M2.2Tau0.9 ± 0.3−2.6 ± 0.013658 ± 1364.06 ± 0.25<2.5334.4 ± 2.60.15 ± 0.05
J04174965+2829362M3.7Tau1.4 ± 0.3−0.8 ± 0.103315 ± 1293.37 ± 0.251.22 ± 0.339.7 ± 2.00.28 ± 0.06
J04181078+2519574M1.0Tau−0.7 ± 0.3−0.7 ± 0.103721 ± 1493.67 ± 0.23<1.1215.6 ± 2.20.32 ± 0.07
J04182909+2826191M3.5Tau⋯⋯13618 ± 3274.05 ± 0.59<2.3230.7 ± 7.40.39 ± 0.20
J04183112+2816290M3.5Tau2.0 ± 0.6−0.5 ± 0.103250 ± 2073.47 ± 0.37<1.4816.9 ± 3.11.82 ± 0.26
J04183158+2816585M4.2Tau1.3 ± 0.70.2 ± 0.213385 ± 2103.83 ± 0.35<1.7031.4 ± 3.10.64 ± 0.12
J04183444+2830302M0.0Tau⋯⋯04230 ± 1424.53 ± 0.271.84 ± 0.4910.9 ± 3.00.44 ± 0.08
J04184703+2820073K8.0Tau1.8 ± 0.1−2.6 ± 0.013806 ± 813.69 ± 0.142.63 ± 0.2716.8 ± 1.70.26 ± 0.01
J04191281+2829330M3.5Tau1.4 ± 0.4−1.2 ± 0.213296 ± 1703.93 ± 0.30<1.4318.9 ± 2.40.63 ± 0.10
J04191583+2906269M0.5Tau0.2 ± 0.3−0.9 ± 0.103719 ± 1304.01 ± 0.252.19 ± 0.379.9 ± 2.31.36 ± 0.11
J04192625+2826142K8.0Tau1.1 ± 0.2−2.2 ± 0.203830 ± 1243.91 ± 0.212.22 ± 0.309.8 ± 2.00.20 ± 0.03
J04194127+2749484M0.0Tau0.9 ± 0.2−2.7 ± 0.003853 ± 1534.08 ± 0.26<1.9415.6 ± 2.50.19 ± 0.05
J04202606+2804089M3.5Tau0.4 ± 0.6−0.5 ± 0.303407 ± 1974.09 ± 0.351.37 ± 0.6310.5 ± 2.80.22 ± 0.10
J04214323+1934133M2.4Tau3.5 ± 0.6−0.5 ± 0.123543 ± 1683.27 ± 0.271.04 ± 0.357.9 ± 2.20.70 ± 0.09
J04215563+2755060M2.3Tau0.3 ± 0.5−0.9 ± 0.103463 ± 1573.45 ± 0.27<1.019.4 ± 2.21.11 ± 0.11
J04215943+1932063K0.0Tau1.1 ± 0.40.0 ± 0.104065 ± 2323.51 ± 0.40<2.3421.5 ± 4.24.90 ± 0.85
J04220217+2657304M0.0Tau2.6 ± 0.3−0.3 ± 0.213810 ± 1864.12 ± 0.36<2.3420.8 ± 3.71.38 ± 0.19
J04220313+2825389M2.5Tau1.2 ± 0.4−2.4 ± 0.103659 ± 2264.17 ± 0.36<2.6743.7 ± 3.80.13 ± 0.09
J04221675+2654570M1.5Tau3.7 ± 0.6−0.8 ± 0.103586 ± 2363.83 ± 0.442.20 ± 0.7716.4 ± 3.51.74 ± 0.26
J04244457+2610141M2.8Tau4.1 ± 0.7−0.8 ± 0.013624 ± 3564.01 ± 0.641.55 ± 0.9811.1 ± 5.52.47 ± 0.63
J04245708+2711565M0.6Tau0.7 ± 0.4−1.0 ± 0.103686 ± 1834.16 ± 0.282.46 ± 0.4311.5 ± 2.51.23 ± 0.15
J04265352+2606543M0.0Tau4.7 ± 0.4⋯04034 ± 1823.74 ± 0.35<2.3637.5 ± 3.20.95 ± 0.20
J04265440+2606510M2.5Tau3.5 ± 0.5−1.1 ± 0.103456 ± 1363.53 ± 0.23<2.0138.6 ± 2.40.33 ± 0.06
J04270469+2606163K7.0Tau2.1 ± 0.4−0.0 ± 0.123969 ± 2163.20 ± 0.28<1.7826.9 ± 4.34.82 ± 0.94
J04293606+2435556M3.0Tau3.1 ± 0.6−1.3 ± 0.213557 ± 1723.74 ± 0.28<2.2824.6 ± 2.70.38 ± 0.08
J04294155+2632582M2.3Tau0.3 ± 0.3−1.0 ± 0.303477 ± 1253.89 ± 0.222.21 ± 0.328.4 ± 2.11.24 ± 0.10
J04294247+2632493M0.7Tau0.8 ± 0.2−1.7 ± 0.403675 ± 1023.61 ± 0.17<1.0512.2 ± 1.80.23 ± 0.03
J04295156+2606448M1.1Tau2.6 ± 0.5−1.0 ± 0.103612 ± 1763.72 ± 0.311.81 ± 0.5314.3 ± 2.51.42 ± 0.15
J04300357+1813494M2.0Tau0.6 ± 0.4−2.7 ± 0.003527 ± 1274.07 ± 0.192.41 ± 0.3211.7 ± 2.10.02 ± 0.03
J04300399+1813493K0.0Tau1.9 ± 0.2−1.1 ± 0.404606 ± 2334.43 ± 0.42<1.7026.0 ± 4.51.81 ± 0.44
J04302961+2426450M2.2Tau2.2 ± 0.3−1.1 ± 0.203587 ± 1143.82 ± 0.21<1.028.3 ± 2.00.61 ± 0.07
J04304425+2601244K8.5Tau3.2 ± 0.3−0.9 ± 0.013809 ± 1834.00 ± 0.352.55 ± 0.6417.6 ± 3.33.16 ± 0.29
J04305137+2442222M4.3Tau0.9 ± 0.4−1.2 ± 0.113277 ± 1443.61 ± 0.27<1.5218.4 ± 2.40.41 ± 0.08
J04311444+2710179K8.0Tau0.3 ± 0.2−2.2 ± 0.203802 ± 1394.06 ± 0.262.01 ± 0.4113.9 ± 2.40.12 ± 0.05
J04312382+2410529M4.5Tau0.9 ± 0.6−2.6 ± 0.113203 ± 1714.14 ± 0.37<1.5516.6 ± 3.10.92 ± 0.14
J04314007+1813571M2.0Tau−0.1 ± 0.9−0.3 ± 0.223606 ± 2623.24 ± 0.35<1.8717.1 ± 3.63.85 ± 0.56
J04315056+2424180M1.0Tau2.3 ± 0.6−0.6 ± 0.303598 ± 2033.64 ± 0.35<1.5822.9 ± 2.90.93 ± 0.14
J04315779+1821350M3.3Tau2.6 ± 0.4⋯13492 ± 1613.74 ± 0.29<1.8423.0 ± 2.90.36 ± 0.09
J04315779+1821380M1.7Tau2.2 ± 0.5⋯03518 ± 1623.65 ± 0.26<1.9318.2 ± 2.50.30 ± 0.07
J04320926+1757227K6.0Tau0.4 ± 0.1−2.7 ± 0.004122 ± 964.20 ± 0.18<2.6230.9 ± 2.10.16 ± 0.04
J04321456+1820147M2.0Tau0.7 ± 0.2−2.7 ± 0.013610 ± 913.91 ± 0.15<2.4220.8 ± 1.80.13 ± 0.02
J04321885+2422271M0.8Tau2.0 ± 0.4−2.8 ± 0.103684 ± 1393.65 ± 0.22<1.8832.1 ± 2.30.15 ± 0.05
J04323034+1731406K7.5Tau0.8 ± 0.3−0.8 ± 0.103711 ± 1213.68 ± 0.20<1.1010.4 ± 1.90.98 ± 0.07
J04323058+2419572M0.1Tau3.2 ± 0.5−1.2 ± 0.203925 ± 2223.90 ± 0.371.71 ± 0.7613.0 ± 3.10.78 ± 0.14
J04323176+2420029M0.5Tau3.0 ± 0.8−0.9 ± 0.123751 ± 3183.37 ± 0.481.20 ± 0.888.2 ± 4.65.53 ± 1.05
J04324282+2552314M2.0Tau0.6 ± 0.5⋯03426 ± 1793.79 ± 0.33<1.4718.6 ± 2.50.24 ± 0.10
J04324303+2552311M1.9Tau1.6 ± 1.8−0.8 ± 0.113670 ± 4373.90 ± 0.64<2.2331.0 ± 9.73.06 ± 0.97
J04324373+1802563K6.0Tau0.4 ± 0.2−2.7 ± 0.103912 ± 1264.11 ± 0.231.40 ± 0.357.5 ± 2.20.12 ± 0.04
J04324911+2253027K5.5Tau3.6 ± 0.6−0.8 ± 0.204105 ± 3944.15 ± 0.64<2.4038.2 ± 8.21.58 ± 0.59
J04324938+2253082M4.5Tau3.9 ± 0.9−1.2 ± 0.123380 ± 2473.42 ± 0.42<1.1818.1 ± 3.40.37 ± 0.14
J04325323+1735337M2.0Tau0.5 ± 0.5−2.0 ± 0.203543 ± 1323.65 ± 0.22<1.1011.4 ± 2.00.09 ± 0.04
J04330622+2409339M2.0Tau1.0 ± 0.5−1.1 ± 0.203553 ± 1603.64 ± 0.26<1.5928.1 ± 2.50.57 ± 0.09
J04330664+2409549K7.0Tau1.0 ± 0.3−1.4 ± 0.213880 ± 1774.20 ± 0.352.50 ± 0.7119.9 ± 4.51.16 ± 0.15
J04331003+2433433K7.5Tau0.5 ± 0.1−2.8 ± 0.103878 ± 803.89 ± 0.14<2.4031.6 ± 1.80.16 ± 0.02
J04333405+2421170M0.4Tau1.9 ± 0.4−0.7 ± 0.103689 ± 1853.79 ± 0.321.92 ± 0.4412.1 ± 2.61.42 ± 0.15
J04333456+2421058K6.5Tau1.5 ± 0.3−0.7 ± 0.103842 ± 2033.83 ± 0.33<2.0521.2 ± 3.01.57 ± 0.18
J04333678+2609492M0.0Tau1.9 ± 0.4−0.9 ± 0.103581 ± 1523.98 ± 0.252.28 ± 0.3711.8 ± 2.40.71 ± 0.09
J04333906+2520382K5.5Tau0.9 ± 0.5−0.8 ± 0.103930 ± 2223.80 ± 0.431.95 ± 0.7011.0 ± 3.73.23 ± 0.42
J04334871+1810099M3.0Tau0.1 ± 0.3−0.6 ± 0.503449 ± 1034.09 ± 0.181.63 ± 0.295.7 ± 2.00.08 ± 0.03
J04335200+2250301K5.5Tau1.3 ± 0.2−0.9 ± 0.103951 ± 943.77 ± 0.171.95 ± 0.3112.5 ± 1.92.34 ± 0.09
J04335470+2613275K5.0Tau2.6 ± 0.4−1.2 ± 0.103977 ± 1954.09 ± 0.35<2.8336.4 ± 3.80.71 ± 0.14
J04341099+2251445M1.5Tau1.8 ± 0.5−2.7 ± 0.003589 ± 1543.97 ± 0.262.16 ± 0.4013.3 ± 2.30.09 ± 0.06
J04345542+2428531M0.6Tau4.9 ± 0.3−0.8 ± 0.203751 ± 1713.87 ± 0.262.16 ± 0.4112.5 ± 2.31.56 ± 0.12
J04352020+2232146M3.2Tau0.5 ± 0.7−1.1 ± 0.103493 ± 2404.00 ± 0.48<1.6220.2 ± 3.81.04 ± 0.22
J04352089+2254242K8.0Tau1.9 ± 0.4−2.8 ± 0.103856 ± 1953.87 ± 0.30<0.906.2 ± 2.40.17 ± 0.06
J04352450+1751429M2.6Tau0.9 ± 0.4−2.6 ± 0.013494 ± 1464.28 ± 0.212.44 ± 0.335.0 ± 2.30.04 ± 0.04
J04352737+2414589M0.3Tau0.1 ± 0.4−0.9 ± 0.203648 ± 1443.74 ± 0.221.54 ± 0.3710.9 ± 2.10.43 ± 0.05
J04354093+2411087M0.5Tau4.6 ± 0.4−1.0 ± 0.003665 ± 1543.85 ± 0.262.33 ± 0.348.2 ± 2.10.77 ± 0.08
J04354733+2250216K2.0Tau2.4 ± 0.3−0.7 ± 0.104265 ± 2554.22 ± 0.44<2.8543.0 ± 5.01.18 ± 0.32
J04355277+2254231K4.0Tau3.0 ± 0.5−0.6 ± 0.214096 ± 3804.07 ± 0.62<1.7918.6 ± 8.82.23 ± 0.90
J04355349+2254089M0.6Tau2.3 ± 0.3⋯03760 ± 1643.84 ± 0.26<2.6543.3 ± 2.90.25 ± 0.07
J04355684+2254360M2.0Tau2.0 ± 0.6−0.9 ± 0.103406 ± 2073.99 ± 0.43<1.2212.5 ± 3.10.88 ± 0.20
J04355892+2238353K8.0Tau0.9 ± 0.3−2.7 ± 0.103869 ± 1694.10 ± 0.29<2.5429.1 ± 3.50.16 ± 0.07
J04361909+2542589K5.0Tau0.3 ± 0.2−2.9 ± 0.104134 ± 1264.07 ± 0.24<1.7222.8 ± 2.50.15 ± 0.07
J04382858+2610494M0.3Tau1.3 ± 0.6−0.5 ± 0.123704 ± 3503.44 ± 0.501.78 ± 0.8512.7 ± 5.45.44 ± 1.04
J04390163+2336029M4.9Tau0.1 ± 0.5−1.3 ± 0.113221 ± 1883.86 ± 0.36<1.2717.8 ± 2.70.52 ± 0.12
J04391779+2221034K5.5Tau0.8 ± 0.2−1.1 ± 0.204156 ± 1234.11 ± 0.23<1.8315.4 ± 2.31.08 ± 0.10
J04392090+2545021M2.5Tau3.9 ± 0.4−1.0 ± 0.103536 ± 1523.76 ± 0.28<1.5813.8 ± 2.41.47 ± 0.15
J04400800+2605253 Tau7.2 ± 0.5−0.6 ± 0.103749 ± 2733.84 ± 0.492.32 ± 0.9716.6 ± 4.31.52 ± 0.28
J04410470+2451062M0.9Tau0.7 ± 0.2−2.7 ± 0.003663 ± 1154.07 ± 0.192.28 ± 0.298.3 ± 2.00.10 ± 0.04
J04411681+2840000M1.1Tau2.4 ± 0.3−0.6 ± 0.603706 ± 1213.70 ± 0.21<1.9026.3 ± 2.20.22 ± 0.04
J04413882+2556267M0.0Tau3.3 ± 0.3−0.1 ± 0.103738 ± 1613.74 ± 0.26<1.8024.2 ± 2.60.59 ± 0.09
J04420548+2522562M0.5Tau3.0 ± 0.3−2.2 ± 0.103779 ± 1463.89 ± 0.24<2.2221.7 ± 2.60.25 ± 0.05
J04420777+2523118K7.0Tau3.0 ± 0.3⋯03763 ± 1603.93 ± 0.281.59 ± 0.4911.0 ± 2.51.40 ± 0.12
J04423769+2515374M0.8Tau0.7 ± 0.5−0.3 ± 0.103758 ± 2823.91 ± 0.43<1.8919.9 ± 4.11.60 ± 0.25
J04430309+2520187M2.3Tau1.0 ± 0.6−0.9 ± 0.203487 ± 1813.89 ± 0.31<1.3012.3 ± 2.50.49 ± 0.11
J04465305+1700001M0.6Tau1.3 ± 0.8−0.7 ± 0.113704 ± 3193.88 ± 0.57<1.8617.7 ± 4.90.59 ± 0.23
J04465897+1702381K8.0Tau2.4 ± 0.4−1.1 ± 0.103903 ± 2433.53 ± 0.35<1.3012.8 ± 3.12.36 ± 0.36
J04474859+2925112M0.4Tau0.9 ± 0.2−1.2 ± 0.103879 ± 1383.89 ± 0.252.09 ± 0.4113.4 ± 2.31.79 ± 0.13
J04514737+3047134K7.0Tau0.8 ± 0.6−0.2 ± 0.003780 ± 2793.79 ± 0.49<1.8021.4 ± 4.33.03 ± 0.45
J04551098+3021595K6.0Tau0.9 ± 0.3−0.6 ± 0.504020 ± 994.02 ± 0.192.11 ± 0.3414.4 ± 2.00.48 ± 0.05
J04553695+3017553K2.0Tau−0.1 ± 0.1−2.5 ± 0.104450 ± 1224.26 ± 0.25<1.9521.8 ± 2.80.24 ± 0.09
J04560201+3021037K6.0Tau0.6 ± 0.3−2.4 ± 0.104034 ± 1003.99 ± 0.19<1.6114.2 ± 2.00.18 ± 0.04
J05030659+2523197M0.8Tau1.1 ± 0.2−1.0 ± 0.213760 ± 1584.01 ± 0.262.67 ± 0.4011.8 ± 2.40.68 ± 0.08
J05071206+2437163K6.0Tau0.5 ± 0.3−2.6 ± 0.003987 ± 1314.11 ± 0.25<2.0021.6 ± 2.50.09 ± 0.05
J05074953+3024050K2.0Tau3.4 ± 0.7−0.8 ± 0.103885 ± 2154.10 ± 0.39<2.0216.6 ± 3.92.09 ± 0.29
J10120908-3124451M4.0TWA1.1 ± 0.2−1.6 ± 0.113315 ± 874.06 ± 0.18<1.3418.0 ± 2.00.28 ± 0.03
J10423011-3340162M3.2TWA⋯⋯03327 ± 804.27 ± 0.132.19 ± 0.267.4 ± 1.70.10 ± 0.00
J11015191-3442170M0.5TWA⋯⋯03783 ± 1074.35 ± 0.182.75 ± 0.308.4 ± 2.00.53 ± 0.05
J11091380-3001398M2.2TWA⋯⋯03557 ± 924.07 ± 0.16<1.3815.9 ± 1.90.06 ± 0.02
J11102788-3731520M4.1TWA⋯⋯13284 ± 754.00 ± 0.13<0.9711.9 ± 1.70.29 ± 0.00
J11210549-3845163M2.75TWA⋯⋯03534 ± 1044.19 ± 0.172.67 ± 0.3219.7 ± 2.00.02 ± 0.02
J11211723-3446454M1.1TWA⋯⋯03638 ± 954.14 ± 0.162.12 ± 0.2914.2 ± 1.90.04 ± 0.02
J11211745-3446497M1.0TWA⋯⋯03672 ± 934.13 ± 0.16<1.7113.7 ± 1.90.04 ± 0.02
J11220530-2446393K6.0TWA⋯⋯04257 ± 1024.51 ± 0.21<1.1010.2 ± 2.10.31 ± 0.04
J11324124-2651559M2.9TWA⋯⋯03397 ± 754.30 ± 0.132.86 ± 0.267.5 ± 1.70.13 ± 0.00
J11482373-3728485M3.4TWA0.2 ± 0.2−2.7 ± 0.103351 ± 754.30 ± 0.14<1.3010.9 ± 1.70.09 ± 0.00
J11482422-3728491K6.0TWA0.1 ± 0.3−2.7 ± 0.004043 ± 944.32 ± 0.172.29 ± 0.3012.0 ± 2.00.17 ± 0.03
J12072738-3247002M3.5TWA⋯⋯13404 ± 1014.15 ± 0.18<1.4219.1 ± 1.90.16 ± 0.03
J12153072-3948426M0.5TWA⋯⋯03707 ± 994.14 ± 0.172.21 ± 0.2915.0 ± 1.90.06 ± 0.03
J12345629-4538075M3.0TWA0.0 ± 0.4−2.7 ± 0.003445 ± 904.15 ± 0.161.85 ± 0.2812.4 ± 1.90.07 ± 0.02
J12350424-4136385M3.0TWA⋯⋯13357 ± 754.23 ± 0.132.54 ± 0.269.5 ± 1.70.10 ± 0.00
J12354893-3950245M4.5TWA0.9 ± 0.4−2.5 ± 0.013391 ± 1224.45 ± 0.22<1.9822.8 ± 2.30.36 ± 0.04
J12360055-3952156M2.5TWA⋯⋯03531 ± 1154.22 ± 0.182.23 ± 0.3414.8 ± 2.10.02 ± 0.02
TWA 3BM4.0TWA⋯⋯13355 ± 754.20 ± 0.13<1.7815.8 ± 1.70.16 ± 0.00

Note. Columns 1–3 provide target names, spectral types, and membership information. Columns 4 and 5 are the values of the visual extinction and the infrared index computed in this work. Finally, Columns 6–11 present the quality flag, T, $\mathrm{log}g$, B field, $v\sin i$, and K-band veiling determined in Paper I.

A machine-readable version of the table is available.

Download table as:  DataTypeset images: 1 2 3

A.2. Field Stars

Between 2014 and 2018, as part of various scientific projects, IGRINS observed 254 K and M field stars at the McDonald Observatory 2.7 m telescope, Lowell Discovery Telescope, and Gemini South Telescope. López-Valdivia et al. (2019) presented an accurate determination of their temperatures using absorption line-depths of Fe I, OH, and Al I, located at the H-band and BT-Settl models (Allard et al. 2013).

These K and M stars are an excellent sample of more evolved counterparts for our YSO sample. The field stars allow us to complete the age ladder we are creating with the IGRINS YSO survey, establish the differences or similarities between the populations at different evolutionary stages, and identify the role of the stellar parameters at every step.

We started selecting from López-Valdivia et al. (2019) just the stars with spectral type between K0 and M5 and with a minimum mean signal-to-noise ratio of ∼50 in the K band, finishing with a sample of 207 field stars. Then, we determined the stellar parameters in the field stars following the MCMC approach as in our YSO sample, but with minor variations.

The main differences between the MCMC applied to the field stars and the YSOs are that (i) due to the more evolved status of the field stars, there is no need to include any K-band veiling contribution, so we removed this free parameter from the MCMC analysis. (ii) We extended the synthetic spectral grid to include $\mathrm{log}g$ values up to 5.5 dex, as M field stars might have $\mathrm{log}g$ values close to or even higher than 5.0 dex (e.g., Ségransan et al. 2003). The MARCS library has models with $\mathrm{log}g$ equal to 5.5 dex for T lower than 3900 K. We extended the synthetic grid for T colder than 3900 K.

As with the YSOs, we classified the stellar parameter determination of the field stars with the same quality flag as the YSOs. We assigned a flag equal to 0, 1, 2, or 3 if the determination is good, acceptable, at the edge of our grid, or poor. We found that almost 36% of our field stars (74 of 207) have poor determinations. The reason for this high number is not clear. We identified some common issues in the stellar spectra that prevented our MCMC analysis from converging into a good fit. We identify some stars with very high $v\sin i$, narrow atomic or molecular lines, and some double lines. The remaining 133 stars, reported in Table 4, are considered our field sample, and we used them throughout the paper.

Table 4. Stellar Parameters of the K and M Field Stars

NameSpTFlag T $\mathrm{log}g$ B $v\sin i$
[RSP2011] 315⋯03558 ± 1864.84 ± 0.31<0.424.7 ± 2.9
G 43-43⋯03655 ± 1444.90 ± 0.27<0.324.0 ± 2.5
LP 611-70⋯03494 ± 1494.24 ± 0.20<0.423.4 ± 2.0
UCAC4 545-148763⋯03793 ± 1404.59 ± 0.221.45 ± 0.395.5 ± 2.7
G 194-18⋯23748 ± 1265.32 ± 0.25<0.404.9 ± 3.5
HD 285482K0.004744 ± 914.63 ± 0.17<0.343.9 ± 2.0
HD 285690K0.004843 ± 964.61 ± 0.19<0.334.3 ± 2.2
HD 285876K0.004190 ± 1114.58 ± 0.19<0.263.4 ± 2.1
HD 286363K0.004606 ± 934.54 ± 0.20<0.364.5 ± 2.1
BD+45 598K0.024924 ± 1114.12 ± 0.37<0.6322.5 ± 3.6
HD 182488K0.024983 ± 784.09 ± 0.22<0.244.6 ± 2.1
* 54 PscK0.524993 ± 754.27 ± 0.15<0.093.4 ± 1.8
* 107 PscK1.024996 ± 754.25 ± 0.15<0.055.0 ± 1.8
HD 125455K1.024990 ± 764.10 ± 0.18<0.134.8 ± 1.9
HD 285348K2.004637 ± 904.56 ± 0.18<0.373.3 ± 2.1
HD 3765K2.004815 ± 804.28 ± 0.15<0.092.3 ± 1.7
BD+20 2720K2.024980 ± 793.31 ± 0.24<0.153.9 ± 2.0
HD 21845K2.024957 ± 994.36 ± 0.28<0.5218.4 ± 2.9
HD 220339K2.024936 ± 1283.98 ± 0.46<0.598.4 ± 4.4
HD 88925K2.024992 ± 753.14 ± 0.19<0.065.4 ± 1.8
HD 122064K3.004554 ± 774.24 ± 0.14<0.082.5 ± 1.7
HD 219134K3.004687 ± 904.34 ± 0.18<0.183.0 ± 1.9
HD 52919K4.024616 ± 914.48 ± 0.21<0.212.9 ± 2.0
BD+49 2125K4/514535 ± 1304.36 ± 0.30<0.8212.0 ± 2.6
BD-09 2926BK514584 ± 1064.61 ± 0.25<0.284.0 ± 2.3
BD+20 1790K5.004352 ± 3354.77 ± 0.222.61 ± 0.4910.8 ± 2.3
HD 286554K5.004067 ± 1174.51 ± 0.18<0.283.3 ± 2.0
BD-09 4191K5.013719 ± 1505.16 ± 0.30<0.397.8 ± 3.6
HD 98800K5.014340 ± 904.20 ± 0.18<0.5910.3 ± 2.0
HD 122120K5.024369 ± 864.32 ± 0.16<0.182.6 ± 1.8
BD+05 378K6.004144 ± 2904.67 ± 0.412.27 ± 0.7212.1 ± 3.3
HD 88230K6.003843 ± 1044.31 ± 0.16<0.263.2 ± 2.0
HD 283869K7.004369 ± 834.28 ± 0.15<0.152.7 ± 1.8
HD 6440BK7.003854 ± 974.77 ± 0.16<0.173.3 ± 1.9
LP 533-57K7.003873 ± 834.36 ± 0.15<0.294.5 ± 1.9
HD 21845BK7.013758 ± 1984.53 ± 0.35<2.7226.8 ± 4.6
HD 157881K7.023877 ± 804.37 ± 0.14<0.292.5 ± 1.8
HD 97101K7.023874 ± 804.25 ± 0.14<0.132.5 ± 1.8
BD+06 2986K8.013846 ± 885.09 ± 0.16<0.123.1 ± 1.9
BD-13 6424M0.003682 ± 1764.49 ± 0.261.96 ± 0.419.1 ± 2.8
BD+02 1729M0.003743 ± 914.44 ± 0.15<0.283.6 ± 1.9
BD+13 2618M0.003631 ± 1024.37 ± 0.162.21 ± 0.2912.5 ± 1.9
BD+23 2063BM0.003749 ± 1424.23 ± 0.21<0.433.7 ± 2.2
BD+33 1505M0.003692 ± 1174.58 ± 0.18<0.373.5 ± 2.0
BD+36 2322M0.003732 ± 1234.73 ± 0.212.22 ± 0.3612.8 ± 2.4
HD 19305M0.003860 ± 804.54 ± 0.14<0.152.6 ± 1.8
HD 79211M0.003759 ± 1284.66 ± 0.19<0.233.1 ± 2.0
HIP 17248M0.003754 ± 1944.54 ± 0.252.24 ± 0.435.4 ± 2.8
Ross 1050M0.003793 ± 1255.10 ± 0.24<0.284.7 ± 2.7
StKM 1-82M0.003816 ± 1264.76 ± 0.21<0.495.1 ± 2.4
TYC 5174-242-1M0.003749 ± 1374.69 ± 0.24<0.689.9 ± 2.6
BD+30 397BM0.013563 ± 1474.74 ± 0.242.84 ± 0.378.1 ± 2.3
CCDM J10095+4126ABM0.014768 ± 1154.71 ± 0.24<0.899.6 ± 2.6
GSC 04686-00596M0.014001 ± 2104.53 ± 0.31<2.8737.9 ± 3.8
HD 59438CM0.013821 ± 1205.26 ± 0.26<0.2610.8 ± 3.1
StKM 1-516M0.013826 ± 1025.16 ± 0.18<0.133.6 ± 2.1
HD 209290M0.503662 ± 1004.55 ± 0.16<0.322.8 ± 1.9
HD 28343M0.503821 ± 964.19 ± 0.15<0.313.0 ± 1.9
GJ 3653M0.513794 ± 1454.77 ± 0.29<2.1834.6 ± 3.8
MCC 124M0.713847 ± 1045.16 ± 0.201.60 ± 0.539.2 ± 3.0
2MASS J14040922+2044314M1.003655 ± 1414.78 ± 0.24<0.785.3 ± 2.5
BD+61 1678CM1.003627 ± 1234.62 ± 0.19<0.343.6 ± 2.1
HD 199305M1.003600 ± 884.60 ± 0.15<0.212.5 ± 1.8
HD 42581M1.003589 ± 1474.50 ± 0.21<0.403.5 ± 2.2
LP 353-51M1.003772 ± 1154.79 ± 0.181.59 ± 0.324.0 ± 2.2
BD+02 348M1.013882 ± 775.46 ± 0.14<0.073.3 ± 1.9
HD 263175BM1.013845 ± 1025.35 ± 0.20<0.214.0 ± 2.4
BD+44 2051M1.023872 ± 815.47 ± 0.14<0.103.4 ± 2.0
2MASS J00393579-3816584M1.413535 ± 1164.55 ± 0.172.54 ± 0.329.0 ± 1.9
* tet Per BM1.503584 ± 964.76 ± 0.16<0.252.8 ± 1.9
BD+18 3421M1.503856 ± 1085.38 ± 0.18<0.075.5 ± 1.9
HD 31867BM1.503589 ± 1824.71 ± 0.30<0.364.5 ± 2.8
LP 298-53M1.503622 ± 1594.92 ± 0.261.77 ± 0.396.1 ± 2.6
NLTT 19115M1.503703 ± 1305.01 ± 0.21<0.243.7 ± 2.2
Wolf 58M1.503648 ± 1534.69 ± 0.221.82 ± 0.356.8 ± 2.4
HD 154363BM1.513771 ± 1085.13 ± 0.21<0.123.6 ± 2.1
G 202-48M1.523852 ± 945.46 ± 0.14<0.155.1 ± 2.3
LP 493-31M1.523591 ± 1295.12 ± 0.21<0.183.2 ± 2.2
BD+01 2447M2.003684 ± 955.11 ± 0.17<0.123.6 ± 1.9
BD+25 3173M2.003599 ± 1084.73 ± 0.18<0.213.2 ± 1.9
BD+44 3567M2.003716 ± 1095.14 ± 0.19<0.163.5 ± 2.1
HD 119850M2.003693 ± 875.06 ± 0.16<0.063.3 ± 1.8
HD 217987M2.003742 ± 1215.13 ± 0.21<0.213.8 ± 2.2
V* AN SexM2.003550 ± 1244.67 ± 0.20<0.323.3 ± 2.1
Ross 799M2.013476 ± 1564.61 ± 0.26<0.364.0 ± 2.4
HD 349726M2.023876 ± 805.46 ± 0.14<0.093.9 ± 2.1
Ross 730M2.023875 ± 805.47 ± 0.14<0.083.1 ± 1.9
* 18 Pup BM2.503585 ± 1564.87 ± 0.28<0.325.2 ± 2.6
* tau Cyg CM2.503696 ± 1445.16 ± 0.26<0.354.7 ± 2.8
G 145-11M2.503686 ± 1664.98 ± 0.31<0.384.7 ± 3.1
HD 285968M2.503477 ± 1284.64 ± 0.20<0.363.5 ± 2.2
HD 38529BM2.503404 ± 1784.47 ± 0.27<0.393.9 ± 2.4
HD 50281BM2.503633 ± 1095.00 ± 0.18<0.263.5 ± 2.1
Ross 265M2.503512 ± 1124.67 ± 0.21<0.418.3 ± 2.0
UCAC4 315-070111M2.503634 ± 1165.15 ± 0.212.14 ± 0.389.6 ± 2.4
G 155-29M2.513633 ± 1294.99 ± 0.23<0.176.2 ± 2.2
LP 222-50M2.523776 ± 1155.35 ± 0.21<0.223.8 ± 2.3
LP 636-19M2.913554 ± 1244.94 ± 0.26<0.7115.0 ± 2.5
G 121-42M3.003681 ± 1065.12 ± 0.19<0.103.8 ± 2.0
GJ 752M3.003492 ± 994.68 ± 0.18<0.213.3 ± 1.9
GJ569M3.003651 ± 1294.99 ± 0.22<0.995.1 ± 2.5
Ross 905M3.003569 ± 1074.91 ± 0.18<0.143.2 ± 1.9
BD+68 946M3.013600 ± 994.99 ± 0.17<0.103.3 ± 2.0
G 80-21M3.013603 ± 2104.92 ± 0.362.51 ± 0.799.2 ± 5.6
HD 173739M3.013878 ± 795.48 ± 0.13<0.062.8 ± 1.8
UCAC4 595-047332M3.013675 ± 1495.26 ± 0.261.67 ± 0.557.0 ± 3.3
G 83-44M3.4503494 ± 1164.77 ± 0.21<0.193.5 ± 2.1
* ksi Peg BM3.503749 ± 1395.23 ± 0.26<0.315.0 ± 2.9
BD+05 1668M3.513746 ± 1095.32 ± 0.19<0.063.5 ± 1.9
CD-44 11909M3.513433 ± 1074.94 ± 0.22<0.157.1 ± 2.0
HD 127871BM3.513713 ± 1545.24 ± 0.26<0.244.6 ± 2.7
BD-15 6290M3.523467 ± 1015.06 ± 0.18<0.162.7 ± 1.9
G 106-36M3.523680 ± 1525.28 ± 0.29<0.316.7 ± 2.9
HD 173740M3.523882 ± 785.49 ± 0.13<0.064.5 ± 2.1
V* CW UMaM3.523812 ± 1245.39 ± 0.193.60 ± 0.495.2 ± 3.5
Wolf 1062M3.523883 ± 785.48 ± 0.13<0.046.9 ± 1.8
Ross 104M4.003648 ± 975.10 ± 0.18<0.173.7 ± 2.0
BD+20 2465M4.013548 ± 915.11 ± 0.172.55 ± 0.332.9 ± 1.9
G 96-10M4.013603 ± 1285.25 ± 0.24<0.353.6 ± 2.3
HD 18143CM4.013515 ± 1884.93 ± 0.36<0.484.9 ± 3.2
LP 642-48M4.013614 ± 1375.21 ± 0.243.44 ± 0.3614.2 ± 2.3
Ross 689M4.013570 ± 1255.05 ± 0.21<0.183.6 ± 2.2
UCAC4 468-040412M4.013575 ± 1525.21 ± 0.24<0.334.0 ± 2.4
V* WW PsAM4.013578 ± 1314.89 ± 0.262.39 ± 0.5316.1 ± 2.6
V* YZ CMiM4.013606 ± 1115.36 ± 0.183.81 ± 0.293.3 ± 2.0
Wolf 358M4.013463 ± 1115.05 ± 0.20<0.163.4 ± 2.1
Wolf 414M4.013493 ± 1064.95 ± 0.22<0.138.0 ± 2.0
UCAC4 536-150368M4.023798 ± 1305.43 ± 0.222.22 ± 2.0715.6 ± 5.5
2MASS J23263962+4521141M4.503766 ± 1494.95 ± 0.32<0.629.5 ± 4.1
* rho01 Cnc BM4.523391 ± 1084.89 ± 0.21<0.142.8 ± 1.9
2MASS J04324938+2253082M4.613415 ± 2714.07 ± 0.45<0.4220.6 ± 3.6
V* TX PsAM5.013720 ± 1875.04 ± 0.36<0.3327.2 ± 2.7
LP 214-42M5.023696 ± 2295.13 ± 0.49<1.0745.7 ± 4.7

Note. The stellar parameters are determined using the same MCMC approach, but ignoring K-band veiling. The spectral type references are those indicated in López-Valdivia et al. (2019).

A machine-readable version of the table is available.

Download table as:  DataTypeset images: 1 2 3

In Figure 8, we compared the T values determined in here with those values reported by López-Valdivia et al. (2019), color coded by our $\mathrm{log}g$ values. There is an acceptable agreement between the two studies, but we found a trend with our $\mathrm{log}g$ values. The line-depth temperatures are hotter than our values for objects with $\mathrm{log}g$≲4.8 dex, while the opposite behavior occurs if the $\mathrm{log}g$ that we determined is greater than 4.8 dex. Such trend might be the result of two main things: (i) López-Valdivia et al. (2019) used a different set of synthetic spectra (BT-Settl models), and (ii) they fixed the $\mathrm{log}g$ value to 4.5 dex in their determination method for all the sample. Although this, the majority of our determinations (97 of 133) are within ±300 K; however, the differences could be as high as 600 K, especially in those objects that have $\mathrm{log}g$ values significantly different from the fixed 4.5 dex used in López-Valdivia et al. (2019). We also see that the line-depth temperatures piled up about ∼4400 K, where the method loses its effectiveness in determining temperature.

Figure 8.

Figure 8. Comparison between our T values and the line-depth temperatures of López-Valdivia et al. (2019) for 133 field stars, color coded by our $\mathrm{log}g$ determinations. The dashed lines denote ± 300 K from the one-to-one relation (solid line). About 73% of the sample is within ± 300 K, but the differences could be as high as 600 K. The error bar in the right corner represents the mean errors of each method, which are of ∼120 and ∼170 K, for this work and the line-depth methods, respectively.

Download figure:

Standard image High-resolution image

It is important to mention that we selected and fine-tuned the method, synthetic spectra, line list, and the spectral regions used in this paper to determine the stellar parameters of YSOs, which may not work properly for more evolved stars. Despite this, our stellar parameters follow the theoretical predictions for K and M stars.

Footnotes

  • 6  

    An A0V telluric standard star was observed at a similar airmass within 2 hr before or after the science target.

  • 7  

    https://wikis.utexas.edu/display/IGRINS/SNR+Estimates+and+Guidelines

  • 8  

    https://github.com/igrins/plp/tree/v2.1-alpha.3

  • 9  

    The grid steps of T and $\mathrm{log}g$ are defined by the model atmosphere, while we selected the steps of B and $v\sin i$.

  • 10  

    The K-S test computes the probability that two populations come from the same parent distribution.

  • 11  

    We used the classification thresholds of Großschedl et al. (2019) and our ${\alpha }_{{K}_{s}-{\rm{W}}4}$ values. Class II: −1.6 <α < −0.3; Class III: α ≤ −1.6 .

Please wait… references are loading.
10.3847/1538-4357/acab04