Accretion Rates for T Tauri Stars Using Nearly Simultaneous Ultraviolet and Optical Spectra

Accretion Rates for T Tauri Stars Using Nearly Simultaneous Ultraviolet and Optical Spectra


We analyze the accretion properties of 21 low mass T Tauri stars using a dataset of contemporaneous near ultraviolet (NUV) through optical observations obtained with the Hubble Space Telescope Imaging Spectrograph (STIS) and the ground based Small and Medium Aperture Research Telescope System (SMARTS), a unique dataset because of the nearly simultaneous broad wavelength coverage. Our dataset includes accreting T Tauri stars (CTTS) in Taurus, Chamaeleon I, Chamaeleon and the TW Hydra Association. For each source we calculate the accretion rate () by fitting the NUV and optical excesses above the photosphere, produced in the accretion shock, introducing multiple accretion components characterized by a range in energy flux (or density) for the first time. This treatment is motivated by models of the magnetospheric geometry and accretion footprints, which predict that high density, low filling factor accretion spots co-exist with low density, high filling factor spots. By fitting the UV and optical spectra with multiple accretion components, we can explain excesses which have been observed in the near infrared. Comparing our estimates of to previous estimates, we find some discrepancies; however, they may be accounted for when considering assumptions for the amount of extinction and variability in optical spectra. Therefore, we confirm many previous estimates of the accretion rate. Finally, we measure emission line luminosities from the same spectra used for the estimates, to produce correlations between accretion indicators (H, Ca II K, C II] and Mg II) and accretion properties obtained simultaneously.

Accretion, accretion disks, Stars: Pre Main Sequence, Ultraviolet: Stars, Stars: Chromospheres

1 Introduction

Classical T Tauri stars (CTTS) are pre-main sequence objects that are accreting gas from their circumstellar disks. The currently accepted paradigm is magnetospheric accretion, where the circumstellar disk is truncated at a few stellar radii by the stellar magnetosphere (Hartmann et al., 1994; Bouvier et al., 1995; Johns-Krull et al., 1999; Johns-Krull & Gafford, 2002; Muzerolle et al., 1998, 2001), with strong magnetic field strengths of a few kG (Johns-Krull et al., 2000). The gas at the truncation radius is channeled along the magnetic field lines at nearly free fall velocities, until it impacts the stationary photosphere, producing an accretion shock. Accreting sources are identified by broad emission lines (like H, Ca II and FUV lines of C IV and He II), tracing the fast moving material in the accretion streams (Edwards et al., 1994; Muzerolle et al., 1998; Ardila et al., in press), or by an ultraviolet (UV) and band excess over the stellar photosphere, produced by hot gas in the accretion shock (Calvet & Gullbring, 1998, hereafter CG98).

CTTS are known to vary in brightness over short timescales, due to changes in the accretion luminosity and the modulation of spots on the stellar surface (Herbst et al., 1994; Alencar et al., 2012). In addition, a spread in the accretion rate is observed for sources of the same age (Hartmann et al., 1998; Calvet et al., 2005; Sicilia-Aguilar et al., 2010). Therefore, the accretion properties of young stars are best characterized by studying a large sample, including non-accreting T Tauri stars (WTTS) for comparison. WTTS are no longer accreting, but they retain strong magnetic fields, due to magnetic dynamo effects, which heat the chromosphere. The ramped up activity in the chromosphere is seen as a UV excess over a dwarf star photosphere (Houdebine et al., 1996) and emission in lines of hydrogen and calcium, which are the same diagnostics of magnetospheric accretion; however, the chromospheric contribution is much weaker (Ingleby et al., 2011). Due to the similarity in tracers, it is essential to use a WTTS template when estimating the accretion rate onto the star (), especially for low objects where the accretion emission is comparable to that from the active chromosphere. Until now, high signal to noise UV spectra of WTTS were not available so dwarf spectra, with lower chromospheric activity levels, took their place as stellar templates, resulting in overestimated UV excesses. Accurate s are vital for the understanding of disk physics. The accretion rate provides information about the surface density of the inner circumstellar disk; for the lowest accretors, the accretion properties reveal the characteristics of the final stages of the inner disk, shortly before it is dissipated (D’Alessio et al., 1999; Ingleby et al., 2011).

Initial attempts to fit the UV excess produced in the accretion shock assumed accretion spots on the stellar surface were characterized by either a single temperature and density in treatments of the shock as a slab (Valenti et al., 1993; Gullbring et al., 1998; Herczeg & Hillenbrand, 2008; Rigliaco et al., 2012), or a single energy flux in the accretion column for accretion shock models (Calvet & Gullbring, 1998, see Section 3.4). These models indicated that the accretion spots cover a small fraction of the stellar surface, tenths to a few percent. Recent models of the magnetosphere, calculated to reproduce spectro-polarimetric observations, revealed that the magnetic field geometry of accreting stars is complex, including high order and tilted fields (Donati et al., 2008; Gregory & Donati, 2011; Gregory et al., 2012). The accretion footprints on the stellar surface which result from field lines of varying strengths and geometries are themselves complex, ranging in size and density of the shocked material, with filling factors which may exceed 10% of the stellar surface (Long et al., 2011). These large filling factors have yet to be reproduced when fitting the accretion excess without significantly over-estimating the UV fluxes.

In addition to the large filling factors, near infrared (IR) veiling observations are not explained by current accretion shock models. Veiling occurs when the excess emission produced in the shock fills in, or “veils”, photospheric absorption lines, causing them to appear shallow when compared to a WTTS template (Hartigan et al., 1991; Valenti et al., 1993; Johns-Krull et al., 1999; Dodin & Lamzin, 2012). Veiling at progressively redder wavelengths can be found in the literature, out to 1 m (Basri & Batalha, 1990; Hartigan et al., 1991; White & Hillenbrand, 2004; Edwards et al., 2006). More recent results on near IR veiling were discussed in Fischer et al. (2011) and McClure et al. (in prep), who showed that the amount of veiling becomes nearly constant from 0.8 to 1.4 m. For an accretion shock spectrum which peaks in the UV and decreases toward long wavelengths, veiling should be negligible in the near-IR. The spectrum of the disk, assuming it is produced by circumstellar dust at the dust sublimation temperature around 1400 K, also has little contribution to the veiling near 1 m and does not begin to contribute significantly until 2 m (Fischer et al., 2011). Long wavelength veiling may originate in a cool accretion component, cooler than that which describes the UV excess (Calvet & Gullbring, 1998; White & Hillenbrand, 2004); however, Fischer et al. (2011) found that in some cases the surface area of the cool accretion column required to explain the veiling was larger than the stellar surface, so they suggested that the veiling came from hot gas in the inner circumstellar disk. On the other hand, by varying the sublimation temperature of the dust over a range consistent with different materials, McClure et al. (in prep) showed that the veiling between 0.8 and 2.32 microns could be explained by the combined emission of dust at the sublimation radius and a cool accretion column with a reasonable filling factor.

In this paper, we assume that 1 m veiling is produced by cool accretion components and use the accretion shock models described in CG98 to fit UV and optical observations of a large sample of CTTS. The sample is a part of the large Hubble Space Telescope (HST) program, Disks, Accretion and Outflows of T Tau stars (DAO; PI G. Herczeg) which compiled a dataset for each source covering a long range in wavelength with observations as close to simultaneous as possible. Observations included HST far UV (FUV), near UV (NUV) and optical spectra obtained with the Cosmic Origins Spectrograph (COS) and the Space Telescope Imaging Spectrograph (STIS). We focus on the NUV and optical observations which are ideal for measuring accretion excesses, whereas the FUV spectrum is complicated by molecular emission, both in lines and the continuum (Bergin et al., 2004; Ingleby et al., 2011; France et al., 2011). Attempts were made to contemporaneously (within a few nights of the HST observations) obtain additional optical spectra in order to observe H, which was not covered in the STIS spectra; these attempts were successful for over half of the sample (Section 2.3).

Here, we improve upon treatments of the accretion shock as a single spot on the stellar surface by including multiple accretion components, a scenario more physically accurate given the geometry of the magnetosphere. With the added components, we attempt to explain long wavelength veiling by including cooler accretion columns, as suggested by White & Hillenbrand (2004) and CG98, from which the emission peaks toward redder wavelengths. In Section 2 we discuss the data used in this paper, both new observations and some previously published. Section 3 describes the process we use for calculating accretion rates and our results are presented in Section 4. Finally, in Section 4.1 we explore the range of accretion rates possible, assuming that some accretion luminosity is undetectable above the star and in Section 5 we show how our new results affect correlations between accretion indicators and , which are commonly used when UV observations are not feasible.

2 Sample and Observations

2.1 Description of Sample

The sample consists of 21 low mass CTTS and 4 WTTS, primarily from the large HST program GO 11616 (PI: G. Herczeg), including three CTTS observed earlier with STIS, BP Tau and TW Hya (GO program 9081; PI. Calvet) and LkCa 15 (GO program 9374; PI. Bergin). The total DAO sample is larger and will be presented in Herczeg (in prep), but here we select those sources which are low mass (spectral types late G to early M). The high mass sources in the sample lack accurate templates which are necessary for our analysis. The sources included in this paper are in the 1–2 Myr old Taurus Molecular Cloud (Kenyon & Hartmann, 1995), the 2–3 Myr old Chamaeleon I star forming region (Gauvin & Strom, 1992), the 9 Myr Chamaeleon region (Lawson et al., 2004), the 10 Myr TW Hydra Association (Webb et al., 1999), and the 16 Myr Lower Centaurus Crux subgroup (Pecaut et al., 2012). Each of these regions is characterized as having isolated low mass star formation, with relatively low extinction, so the sources are not affected by remaining molecular cloud material or nearby high mass stars. The WTTS and CTTS studied in the paper are listed in Tables 1 and 2.

The majority of the sources are single stars; however there are a few wide binaries. DK Tau A, HN Tau A and RW Aur A all have companions at large separations of 1.4″or 200 AU (Kraus et al., 2011; White & Ghez, 2001), resolved by STIS. CS Cha is a spectroscopic binary (Guenther et al., 2007; Nguyen et al., 2012), surrounded by a circumbinary disk (Takami et al., 2003; Nagel et al., 2012). Several of the sources have disks which show evidence for gaps or holes in their infrared spectral energy distributions, including CS Cha, DM Tau, GM Aur and LkCa 15, but all of these sources are still accreting (Strom et al., 1989; Espaillat et al., 2007a, b, 2010).

2.2 Hst Observations

Observations were obtained with STIS between 2009 and 2012. STIS NUV observations used the MAMA detector and the G230L grating providing spectral coverage from 1570 to 3180 Å with R500-1000. Optical observations were completed during the same orbit as the NUV using the G430L grating which covers 2900-5700 Å with R530–1040, resulting in almost simultaneous NUV to optical coverage with STIS. The bright source, CV Cha, was observed with the echelle grating in the NUV (E230M), which we use in place of the G230L spectrum.

The low-resolution STIS spectra were calibrated with custom written IDL routines following the procedures described in the STIS data handbook. The wavelengths were calibrated from the location of identified emission lines within the spectrum, and fluxes were calibrated from spectra of WD 1337+705 in the NUV and HIP 45880 in the optical. The flux calibration also includes a wavelength-dependent aperture correction. Two of the STIS spectra (those of RECX 1 and RECX 11) were previously discussed in Ingleby et al. (2011)

2.3 Ground Based Observations

We also obtained low dispersion spectra of much of the sample, covering the H line, using the Small and Medium Aperture Research Telescope System (SMARTS) 1.5m telescope at the Cerro Tololo Inter-American Observatory (CTIO) with spectral coverage between 5600 – 7000 Å. Sources were observed with SMARTS for a few days before and after the HST observations when possible, providing contemporaneous observations of H. The sample of CTTS with SMARTS coverage includes AA Tau, CS Cha, CV Cha, DE Tau, DM Tau, DN Tau, FM Tau, PDS 66, RECX 11, RECX 15, TWA 3a and V836 Tau. For the following analysis, we used the epoch of SMARTS spectra which was nearest in time to the HST observations. The reduction of SMARTS spectra was described in Ingleby et al. (2011). The flux calibration of the SMARTS spectra was not completely accurate due to flux losses in the slit and changes in the observing conditions. To best estimate the absolute fluxes, we scaled the SMARTS spectra to the HST optical spectra in the region where the two overlap, assuming that the shape of the spectrum is accurate. We assessed the degree of variability in the shape of our SMARTS spectra by computing the standard deviation of the fluxes in the different epochs and dividing it by the median of the observed spectra. With spectral shape variability on the order of 10–25% for our sample (whether due to observing conditions or intrinsic variability), there is little error introduced by fitting only one SMARTS optical spectrum per source.

Due to the low resolution of the SMARTS data, when available, we supplemented the observations with non-simultaneous high resolution spectra obtained with MIKE (Magellan Inamori Kyocera Echelle) on the Magellan-Clay telescope at Las Campanas Observatory in Chile. MIKE has a spectral coverage of 4800–9000 Å and resolution of R35,000. The data were reduced using the Image Reduction and Analysis Facility (IRAF) tasks CCDPROC, APFLATTEN, and DOECSLIT (Tody, 1993). MIKE spectra were used to calculate V band veiling, which would not be accurate at the low resolution of the SMARTS data (Section 3.2). We use MIKE spectra to obtain band veiling for CS Cha, DM Tau, LkCa 15, RECX 11, RECX 15 and TWA 3a. With its location in Chile, Taurus observations are difficult and therefore much of our sample was not covered with MIKE. Two sources were observed with the echelle spectrograph on the SMARTS 1.5m telescope, CV Cha and PDS 66, allowing for veiling estimates.

Additional observations were obtained with CRIRES, the CRyogenic high-resolution InfraRed Echelle Spectrograph on the Very Large Telescope (VLT). CRIRES has a spectral resolution up to in the range of 1 to 5 m. These spectra are used in this paper to calculate veiling for a small number of sources in our sample with unknown long wavelength veiling measurements (see Section 3.2). A description of the CRIRES data reduction may be found in Ingleby et al. (2011). A log of all HST and ground based observations is given in Table 3.

3 Calculating Accretion Rates

3.1 Stellar Template: The Active Chromosphere

Enhanced chromospheric activity is characteristic of stars which have not yet reached the main sequence (Bertout, 1989; Guinan et al., 2003). A higher level of activity is inferred from lines like H and the Ca II infrared triplet, which are in emission and also variable in young stars (Gálvez et al., 2009). WTTS also show excess emission at UV through blue optical wavelengths compared to dwarf stars (Houdebine et al., 1996). Since H, Ca II emission lines and the UV excess are also typical tracers of accretion, emission from the active chromosphere must be taken into account when the accretion luminosity is estimated from any of the aforementioned indicators. While chromospheric emission is not a large contaminant for strong accretors, it is extremely significant when is low. In Ingleby et al. (2011) we demonstrated how the emission from an active chromosphere can mask all evidence of an accretion shock excess for the lowest accretors, though line profiles observed at high resolution revealed complex absorption features produced by the accretion flows.

Previous estimates of have mainly relied on dwarf photospheres as templates against which to measure the U band or UV excess (Romaniello et al., 2004). A better estimate of the excess could be found using a WTTS template (Valenti et al., 1993) but, until now, WTTS with good signal in the UV continuum were not available. As part of the DAO sample, WTTS were included to act as templates, covering the range of spectral types of the accreting sample (Table 1). Figure 1 shows the WTTS templates compared to main sequence dwarf stars with the same spectral type; dwarf spectra were scaled to the WTTS at 5500 Å. STIS observations of main sequence stars were taken from the HST Next Generation Spectral Library (Heap & Lindler, 2007). We compared the luminosity of the WTTS to that of the dwarf between 2000 and 3000 Å and found that the WTTS had NUV luminosities 3 higher than the dwarf stars. Findeisen et al. (2011) showed that WTTS of a given age and spectral type exhibit a range in UV fluxes, so the example WTTS shown in Figure 1 may not exactly represent the chromosphere for each individual source, but they are the best templates to date. In the following analysis, we use the WTTS with the closest spectral type match to each CTTS as the stellar template.

Figure 1: Comparison of WTTS and dwarf stars. In each panel the black line is the WTTS observed as part of the DAO sample and the magenta line is a dwarf standard of the same spectral type, taken from the STIS Next Generation Spectral Library. The WTTS excess in the NUV is produced by an active chromosphere in young stars.

3.2 Veiling By Shock Emission

The excess emission produced in the shock can contribute significantly to the total luminosity of the CTTS, making it difficult to determine the luminosity of the star itself. However, veiled absorption lines provide a diagnostic of the relative contributions from the star and shock (Gullbring et al., 1998). Veiling occurs when excess emission is added to the spectrum of the star, filling in the photospheric absorption lines (Hartigan et al., 1989, 1991), so by comparing the depth of absorption lines in the CTTS to a WTTS, we have an estimate of the veiling continuum. There are some uncertainties in this method; in particular Gahm et al. (2008) and Dodin & Lamzin (2012) show that emission lines can also fill in photospheric absorption lines, so treating the veiling emission as a continuum may not be correct; however, this effect will likely only be important for sources with the strongest emission line spectra, or highest veiling (Petrov et al., 2011).

The continuum flux responsible for the veiling of absorption lines is added to the flux of the photosphere, here taken to be a WTTS, to produce the observed spectrum. So, veiling at one wavelength allows us to estimate the intrinsic photospheric flux at that wavelength from the observed spectrum. If we have an accurate photospheric template, this scaling provides the stellar component flux over the entire spectrum. When available, we used the veiling at , published in Edwards et al. (2006), which provides a compilation from the literature. For a small sample of the sources with unpublished veiling (CS Cha, CV Cha, DM Tau, LkCa 15, PDS 66, RECX 11, RECX 15 and TWA 3a) we calculated from high resolution optical spectra obtained with MIKE or the SMARTS echelle. For templates, we used echelle spectra of WTTS with the same spectral type as each of our sources. We then added a continuum excess to the WTTS spectrum until the depth of the absorption lines matched those in the CTTS, giving us the veiling continuum. We show three examples of different degrees of veiling observed in our MIKE spectra in Figure 2.

Figure 2: Veiling in MIKE spectra. The panels show three CTTS (solid, black) observed with MIKE which have different degrees of veiling compared to a WTTS of the same spectral type (dashed, red). The CTTS are BP Tau (left) with spectral type of K7, CS Cha (middle) with spectral type K6 and RECX11 (right) with spectral type of K5. BP Tau has significantly more veiling than RECX 11, observed as shallower absorption lines. Typical errors on the veiling are 0.1.

Based on the veiling measurements at band, we scaled our WTTS templates to each CTTS using,


where and are the continuum fluxes of the WTTS and CTTS at , respectively, and , is the veiling at V, where is the excess continuum emission added to the photospheric spectrum at . It is important to note that none of our band veilings are simultaneous to the STIS UV and optical spectra. This introduces uncertainties in our estimates of because CTTS accretion properties are known to be variable and as (and the excess continuum emission) change, so does the veiling (Alencar et al., 2012). Ideally, s should be measured using simultaneous long wavelength veiling information to accurately assess the flux contribution from the star and measure the accretion excess.

According to accretion shock models, veiling should drop to 0 in the near IR as the shock spectrum is peaked in the UV (CG98); however, non-zero veiling at long wavelengths has been measured in some CTTS, which we attempt to explain in this paper. Veilings at 1 m were found in the literature for a subset of 10 CTTS in this sample (Edwards et al., 2006). For another five sources, we measured veiling near 1 m using our CRIRES spectra with either the WTTS RECX 1 as a template (for CS Cha, LkCa 15 and RECX 11) or using a Near-Infrared Spectrograph (NIRSPEC) spectrum of a WTTS (Edwards et al., 2006). The NIRSPEC template was used to calculate 1 m veiling for DM Tau, FM Tau and IP Tau, after convolving the CRIRES observations to the lower resolution of NIRSPEC. For CV Cha, PDS 66, RECX 15 and TWA3a, we did not have an appropriate template to calculate the veiling, so we used the relation, (Fischer et al., 2011) to obtain . We also used this relation to estimate for IP Tau, which was not in the literature and for which we did we have high resolution optical spectra; however we did have CRIRES spectra to calculate . All values of and are listed in Table 2

3.3 Extinction and Stellar Parameters

A large source of error in our estimates of comes from the assumed amount of extinction, or . Our STIS dataset did not provide the resolution necessary to measure veiling in the optical and determine based on the veiling, as in Gullbring et al. (1998), so we use extinction estimates from the literature. Table 4 shows that values for a given source vary widely in the literature. Some of the discrepancies in may be due to true variability in the extinction, perhaps by inhomogeneities in the intervening molecular cloud or warps in the circumstellar disk. Variable accretion hot spots can also affect the colors of the star, affecting the determinations (Carpenter et al., 2001, 2002). Gullbring et al. (1998) showed that calculating by comparing the colors of reddened stars to standard photospheric colors results in different values depending on the colors used. In particular, they found that colors are sensitive to the spectral type classification.

A recently developed method to estimate uses IR veiling estimates to determine the shape of the reddened photosphere below the excess and correct accordingly; however, this method has only been applied to a few sources in our sample (see Fischer et al. 2011 and McClure et al. 2013). Both Fischer et al. (2011) and McClure et al. (2012) note that their estimates are higher than the s from frequently quoted sources, like Kenyon & Hartmann (1995). There are still discrepancies among estimates of for the same source, using the IR veiling method (Table 4). For consistency, we use the estimates from Furlan et al. (2009, 2011) when correcting Taurus and Chameleon I sources for extinction. Furlan et al. (2009, 2011), computed extinction values by comparing observed , , or colors to expected photospheric colors in Kenyon & Hartmann (1995), using the reddest colors available to minimize the impact from the shock excess which is typically strongest in the blue. For the remaining sources (TW Hya, TWA 3a, RECX 11 and RECX 15) we assume =0 (Webb et al., 1999; Luhman & Steeghs, 2004). Errors of .5 in can lead to an uncertainty up to 1 order of magnitude in making extinction estimates the largest source of error in our calculations.

The STIS and SMARTS spectra were de-reddened using the reddening law of Whittet et al. (2004). UV emission is extremely sensitive to reddening assumptions, with near 2500 Å and corrections are complicated by uncertain UV reddening laws. Calvet et al. (2004) compared UV reddening laws and found that of Whittet et al. (2004) was appropriate for sources in environments like the Taurus Molecular Cloud. We calculated the stellar luminosity from the flux in the band of the photosphere using the bolometric correction of Kenyon & Hartmann (1995). The photospheric band flux was obtained by scaling de-reddened 2MASS magnitudes by the veiling at 1 m. Although there is an excess above the photosphere at , with known veiling we separate the flux of the star from that of any continuum excess and use the flux from the star alone when calculating the luminosity. We then used the Siess et al. (2000) evolutionary tracks to determine the masses of the sources in the sample and estimate radii from the luminosities. The stellar properties for our WTTS and CTTS are listed in Tables 1 and 2, respectively.

3.4 Accretion Shock Model

A full description of the model used to characterize the accretion shock can be found in CG98 but here we review the main points. The current picture of material accretion in the inner disk of CTTS is magnetospheric accretion. Columns of accreting material fall onto the star along the magnetic field lines traveling at the free fall velocity, , hit the stationary photosphere and create a shock. The velocity of the material is given by


where and are the stellar mass and radius, respectively, and (assumed to be 5) is the radius at which the magnetosphere truncates the disk (CG98).

The model is simplified by assuming that the accretion column has a plane parallel geometry and is perpendicular to the stellar photosphere. The shock formed at the base of the accretion column reduces the velocity of the infalling material in order for it to join the star at the photosphere, converting the kinetic energy into thermal energy, and causing the temperature to increase sharply. In the strong shock approximation, which is used because the material is traveling at high velocities, the temperature immediately after the shock is described by


(CG98). At these high temperatures, the shocked material emits soft X-rays into the pre-shock, post-shock and the photosphere below the shock. This radiation heats the material in these regions causing it to emit the observed excess continuum emission.

In this treatment, the emission from the accretion column is characterized by two parameters, and , which are the total energy flux in the accretion column (, where is the density of the material in the accretion column) and the filling factor (or fraction of the surface of the star which is covered by the accretion hot spot), respectively. Changes in the filling factor, , cause the shock spectrum to increase or decrease independent of wavelength, in essence scaling the luminosity of the emission. acts to change the wavelength of the peak of the accretion shock emission. Figure 3 shows how the shock spectrum, normalized by the maximum flux of each spectrum, changes for different values of , assuming typical values for and of 0.5 and 2.0 , respectively. As increases, more energy is deposited on the stellar surface, increasing the temperature of the photosphere below the shock. This region behaves as the photosphere of a star with an earlier spectral type than that of the undisturbed photosphere. High columns have temperatures up to 9000 K, whereas low columns are cooler, around 5000–6000 K, producing the wavelength shift between columns with different s. For the lowest value of , the emission peaks around 6000–7000 Å whereas for the model with the highest , the emission peaks between 2000 and 3000 Å.

Figure 3: Fluxes of accretion shock models with varying energy flux (). Each line shows the emission from an accretion column with a different value of , normalized to unity. The energy flux of the column is listed next to each spectrum, spanning the range from . The peak of the emission shifts to longer wavelengths as decreases.

We use this property of the accretion column emission to explain excesses at both short and long wavelengths when attempting to fit the DAO observations. In the following text we refer to accretion columns with energy fluxes of as “low ” columns and those characterized by as “high ” columns. depends on the density and velocity of the accreting material; we assume that the velocity is constant, as set by the geometry (see Equation 2), so the range of represents high and low density columns. Note that as decreases, larger values of are needed to produce the same amount of flux as a higher accretion column. can be calculated for each accretion column, with and known, using


(Gullbring et al., 2000). The total is a sum of the contributions from each column and from we determine using,


This treatment does not distinguish between multiple accretion spots, each with a distinct density, or a single accretion spot with a range of densities. The latter scenario is suggested by models produced using the Zeeman Doppler Imaging technique (ZDI) (Gregory & Donati, 2011; Donati et al., 2008).

In this paper, the excess which veils the photospheric absorption lines is primarily a continuum. The main opacity sources included in the calculation are H bound-free and free-free, H- bound-free and free-free, C, Si and Mg, plus additional sources described in Calvet et al. (1991). Line blanketing is also included, using the line list from Kurucz & Bell (1995). A more thorough treatment of the spectrum of the excess continuum including spectral lines has been presented in Dodin & Lamzin (2012) where it was found that emission lines are important to consider when estimating veiling. Dodin & Lamzin showed that veiling may be overestimated when lines are neglected in the veiling calculation. An important difference between the models of CG98 and Dodin & Lamzin (2012) is the structure of the region where the heated photosphere joins the shock above it, which in both treatments is taken at the ram pressure of the shock. In CG98, when the heated photosphere is joined to the postshock region the temperature rises quickly to K because this region has small physical dimensions. Metals would become quickly ionized and few photospheric emission cores would be expected, especially since the location of this transition region is close to the temperature minimum. In contrast, since Dodin & Lamzin (2012) do not have a rapid temperature rise, strong emission lines form in this region. Additional work is needed to sort out the differences between the two models. For now, we acknowledge that emission lines may affect our results, especially for the highly veiled sources like DR Tau and RW Aur A (Gahm et al., 2008; Petrov et al., 2011), and accretion spot sizes for those sources, in particular, may be overestimated due to the omission of emission lines in the veiling spectrum.

4 Results

In this analysis we include, for the first time, multiple accretion columns covering a range in energy flux; specifically, we add low accretion columns to high columns to fit both the UV and optical excesses. Evidence for multiple accretion spots is seen in maps of the magnetosphere (Donati et al., 2008) and also in fits of the broad He and Ca emission, which require inhomogeneity in the hot spots (Dodin et al., 2013). We also compare the amount of 1 m veiling predicted by our models with observed values. Our spectra do not extend into the near IR, so to estimate the 1 m veiling produced by our model we first find of our WTTS template. We scaled by the veiling, , and then using colors for the given spectral type, find the magnitude of the WTTS which reflects the veiling. We assume that , which is valid because the veiling continuum is constant between 0.8 and 1.4 m (Fischer et al., 2011; McClure et al., in prep). We estimate of the model by calculating the flux that would be measured assuming the transmission curve for the Two Micron All Sky Survey (Skrutskie et al., 2006, 2MASS).

To find the best fit of the accretion shock models to the STIS and SMARTS spectra for each source, we calculated the emission from accretion columns spanning erg s cm, using the masses and radii from Table 2. We allowed the filling factor, , to vary independently for each accretion column and then summed the flux from each column with that of the WTTS, scaled to the CTTS by the veiling at (which tells us the location of the photosphere below the excess). The spectrum of the final model is given by;


where is the flux in the model, is the flux from an accretion column with energy flux and is the filling factor of a column with energy flux , where = erg s cm. The values of we have chosen represent the typical range of values for CTTS (CG98).

When determining the best fit model, we isolated continuum regions because the accretion shock models do not attempt to reproduce line emission produced in the shock. We calculated the , as


where is the number of continuum wavelengths which contribute to the fit and is the error in the observed fluxes. The filling factor of each contributing column in the best fit model is given in Table 5; columns with do not contribute to the model fluxes. Finally, we calculated using Equations 2, 4 and,


Estimated s are listed in Table 5, along with the sum of the filling factors of all columns.

The best fit of the accretion shock models to the de-reddened NUV and optical spectra, corresponding to the stellar and accretion properties listed in Tables 2 and 5, are shown in Figures 48. Overall, we found good fits between the models and the observed spectra, with a few exceptions. For the later spectral types, especially the M stars, the models did not reproduce a rise in the observed spectra between 2000 and 3000 Å (see, for example, RECX 15 and DE Tau in Figure 8). This spectral region is populated with Fe emission lines, which are unresolved at the STIS resolution. The Fe lines may be produced in accretion related processes (Herczeg et al., 2005; Petrov et al., 2011); however, they are also observed in the WTTS templates, though weaker, they appear to have a chromospheric component. Fe lines become more apparent at later spectral types, as the photospheric emission in the UV decreases. Finally, we do not attempt to fit the FUV spectrum (1700 Å) which has contributions from H in the disk (Ingleby et al., 2009) that are not included in the accretion shock model.

Figures 48 also include non-simultaneous photometry from the literature, de-reddened using the s from Table 2. When available, we used the range of optical photometry from Herbst et al. (1994) and this range is shown as the green error bars. The photometry for the remaining CTTS came from the following sources; CV Cha (Lawson et al., 1996), GM Aur, IP Tau, LkCa 15, TW Hya and V836 Tau (Kenyon & Hartmann, 1995), PDS 66 (Batalha et al., 1998; Cortes et al., 2009), RECX 11 and RECX 15 (Sicilia-Aguilar et al., 2009) and TWA 3a (Gregorio-Hetem et al., 1992). The models shown do not attempt to fit the photometry because it is not simultaneous, but we use the photometry to look for evidence of high amplitude variability. Most of the photometry agrees with the STIS and SMARTS spectra and is therefore well fit by the models. For a number of the sources, RW Aur A, FM Tau, IP Tau, V836 Tau, RECX 15 and TWA 3A, the STIS spectra are slightly lower than expected from the photometry. Of these sources with mis-matched photometry and spectra, with the exceptions of RW Aur A and FM Tau, photometry for only one epoch is available, so the difference may be due to variability. The DAO NUV and optical spectra of HN Tau A are significantly higher than the range of photometry from Herbst et al. (1994); however, Grankin et al. (2007) observed that the brightness of HN Tau A was steadily decreasing during the 10 years between 1985 and 1995. It may be that Herbst et al. (1994) observed HN Tau A during a dim period and the source has since brightened; such changes in brightness are observed in other CTTS, so it is not unexpected (Grankin et al., 2007).

In Table 5 we also give the values of predicted by our models. Comparing the model to those from the literature or our observations (Table 2), the model values for all but four sources are within 0.1 of the observed , indicating that near IR veiling naturally occurs from a multi-column accretion shock model. Sources with large are fit by models with large filling factors for the low columns. The four sources which do not agree are AA Tau, DR Tau, DK Tau A and HN Tau A. For DK Tau A, we were unable to reproduce the 1 m veiling, likely because of variability. In Edwards et al. (2006), the veiling at is equal to the veiling at 1 m, requiring an abnormally red accretion spectrum, which is not supported by our NUV and optical observations. The and 1 m veiling measurements were not simultaneous, so it is possible that was higher when the IR observations were obtained. HN Tau A also appears to be a remarkably variable object, as mentioned above, so it is possible that variability is responsible for not being able to fit the non-simultaneous optical and near IR veiling data. Veiling is very difficult to measure for DR Tau because the excess continuum is so strong that many absorption lines are completely filled in; therefore, we expect large error bars on the near IR veiling.

Figure 4: Spectra of late G and early K spectral type CTTS in the DAO sample: CV Cha, PDS 66 and RW Aur A. We fit the STIS NUV and optical spectra with emission from accretion columns plus a WTTS spectrum. In each panel the black spectrum is the CTTS and the blue spectrum is the WTTS, LkCa 19. The broken black lines represent accretion shock models with different values, defined as in Figure 3. The red line is the best model fit to the data, adding the emission from the different shock models to the WTTS spectrum. The green error bars indicate non-simultaneous photometry, representing the range observed in multi-epoch observations when available. Although PDS 66 does not have an NUV excess, examination of the H emission line reveals that it is broad, and therefore accretion is likely ongoing, though at low levels.
Figure 5: Spectra of mid K CTTS in the DAO sample: CS Cha, DR Tau, HN Tau A, LkCa 15 and RECX 11. Lines are defined as in Figure 4. The adopted WTTS in each panel is RECX 1. RECX 11 does not have a significant NUV excess, but it was shown in Ingleby et al. (2011) that it is still accreting.
Figure 6: Spectra of late K and early M CTTS in the DAO sample: AA Tau, BP Tau, DK Tau A, DN Tau. Lines are defined as in Figure 4. The adopted WTTS in each panel is HBC 427.
Figure 7: Spectra of late K and early M CTTS in the DAO sample continued: FM Tau, GM Aur, IP Tau, TW Hya and V836 Tau. Lines are defined as in Figure 4. The adopted WTTS in each panel is HBC 427.
Figure 8: Spectra of mid M CTTS in the DAO sample: DE Tau, DM Tau, RECX 15 and TWA 3a. Lines are defined as in Figure 4. The adopted WTTS in each panel is TWA 7.

The models in Figures 48 represent the best fit to the UV and optical data; however there is some degeneracy in and . There is a range in and which will produce a shock spectrum that results in a comparable fit to the data. Figure 9 shows this range for CS Cha and FM Tau, where the colored points represent values of for models with different total filling factors and characteristic values. The characteristic value is the average of the contributing columns, weighted by the filling factor of each column. Figure 9 reveals that there is a range of roughly one order of magnitude in both and for which models will result in fits with the best fit model. While interesting for the shock geometry, this degeneracy has little affect on the ’s because each model is accurately measuring the total excess from which the accretion luminosity is measured.

Figure 9: Degeneracy in characteristic and total filling factor f. The color of each point represents the value of calculated from Equation 7, assuming a model with a given total filling factor and characteristic (the average of the contributing columns weighted by the filling factor) for CS Cha (left) and FM Tau (right). The plot range represents a typical parameter space over which the accretion shock models were calculated. The red box shows the point where is at a minimum.

An important result is that, by including the low columns, we find significantly higher filling factors than would be needed to fit the UV excess alone (Table 5). CG98 fit optical spectra of Taurus CTTS using a single accretion column model. For seven out of nine of our sources which overlap with CG98, we find higher values. The characteristic of the contributing accretion columns is typically lower than that found by CG98, because of the addition of low accretion columns which have large filling factors. However, both DN Tau and AA Tau have larger filling factors and higher values in CG98 than in our analysis. DN Tau was fit by an accretion column characterized by and in CG98; however, their fit to the observed spectrum cannot explain the blue excess very well. Our model includes a high column along with the low column, better fitting both the blue and red excesses and skewing the characteristic to a higher value. We mentioned above that our model for AA Tau did not reproduce near IR veiling measurements. There may be variability in the cool accretion component, where it was less prevalent during the DAO observations. Were we to increase the contribution from low columns to fit the 1 m veiling, it would bring the characteristic to a lower value and the to a larger value than found in CG98. Given that models of the magnetosphere predict large (Long et al., 2011), multi-component accretion columns are a better representation of the physical geometry of the system and now reproduce the flux distribution as well.

We compare the accretion properties calculated here with those from two previous studies of young stars in Taurus in Figure 10 and Table 6. Valenti et al. (1993), hereafter V93, fit optical spectra with coverage from 3400 to 5000 Å of Taurus CTTS using a hydrogen slab model and a WTTS template. The slab was characterized by the temperature, number density, thickness and coverage on the stellar surface. Gullbring et al. (1998), hereafter G98, used similar assumptions to those of V93 and therefore, their accretion rates are in good agreement. G98 notes that the major differences between their s and the V93 s come from choices of and the evolutionary tracks used to determine . Different s change the stellar luminosity, and therefore radius, so the value of will be affected by the extinction estimate. To avoid the additional variables of mass and radius in the analysis, we compare the accretion luminosities from each study instead of .

In Figure 10, left, we plot our values of compared to those of G98 and find that our estimates tend to be higher. The values of we calculate here (shown by the asterisks) have less than a 30% difference from those of the G98 for half of the sample. After correcting our accretion luminosities by the difference in , shown by the arrows, we find that all but one are in good agreement. The outlier is HN Tau A, where our is over one magnitude larger than G98. This is likely due to real variability; we observed that the STIS fluxes are higher than the range in optical photometry observed by Herbst et al. (1994). It appears that two factors in our analysis, the excess flux decrement in the UV due to the WTTS template and the additional excess flux in the red due to the low column, cancel each other out. Therefore, we find similar values of , after accounting for differences in , to those of G98. We chose to use the extinction estimates of Furlan et al. (2009, 2011) who used IR colors to estimate the amount of extinction, avoiding spectral regions where the shock contribution is highest. Furlan et al. (2009, 2011) also covered all of our sources in the Taurus and Chamaeleon I regions, allowing us to use ’s which were derived consistently for each source.

Figure 10: DAO accretion luminosities versus values from the literature. We compare our estimates of to those from Gullbring et al. (1998) on the left and Valenti et al. (1993) on the right. The solid line shows where sources would fall if the accretion luminosities agreed. The dashed lines represent a 30% difference in the values of . Asterisks represent accretion luminosities calculated assuming the value of listed in Table 2 while the arrows show how would change if we use the s from Gullbring et al. (1998) or Valenti et al. (1993).

When comparing the DAO accretion luminosities to those of V93, our estimates are high for three sources; DK Tau A, HN Tau A and DN Tau and low for two; FM Tau and RW Aur A. V93 obtained optical spectra from the UV Schmidt spectrometer at Lick Observatory and provided an atlas of the observations making direct comparison of the optical fluxes possible. For DK Tau A, HN Tau A and DN Tau, the optical fluxes are lower in V93 than in our observations while for RW Aur A and FM Tau they are higher, consistent with the discrepancies in . These variations in the observed fluxes may be intrinsic variability but the difficulty of flux calibrating ground based slit spectra may also contribute. In particular for FM Tau, V93 notes that their slit loss correction was greater than 25%, indicating that it was observed at a high zenith angle or in poor seeing, making flux calibration more uncertain.

4.1 Hidden Accretion Emission

Not all sources require a low energy column to fit the UV and optical excesses in our analysis; however, there could be emission from cold columns hidden by the photospheric flux. If these columns were present, they would increase our estimates of the mass accretion rates. Here, we determine how much hidden flux may be present for sources with no detectable red excess. In Figure 11 we show two model fits to the spectrum of V836 Tau, which has no observed veiling at 1 m (Edwards et al., 2006). In the left panel, the best fit model is produced with only a high column, no additional accretion columns are needed (Table 5). In the right panel, we assume that a low column exists, but the emission from the column is not detectable above the stellar component. For V836 Tau we find that the contribution to the accretion emission that may be hidden below the stellar emission could be equal to that in the high column, doubling the estimated accretion rate. This new model meets the constraint of having limited veiling at 1 m, with 0.1, within the errors of veiling estimates.

We perform this analysis for all sources which had for the low columns in the initial fit, with the constraints that the NUV and optical fluxes are not overestimated by the new models, and that the veiling at 1 m remains within 0.1 of the observed 1 m veiling. Using this method, we calculate upper limits on and the filling factor, listed in Table 7, where the differences in the two filling factors comes from increasing ) from its value in Table 5. These additional low accretion columns are important because they increase the expected along with the filling factors. For the example shown in Figure 11, the accretion rate would double if we include the maximum shock emission that could be hidden by the photosphere, but the filling factor may be up to 100 times higher. The increase of the filling factor is much higher than because there is less mass per unit area in the low columns than in the high columns. This would increase the filling factor from 0.1% to a few percent, better in line with the models of the accretion footprints produced by complex magnetic field geometries (Long et al., 2011) and the distribution of excess accretion related emission regions now being found from magnetic mapping studies (Donati et al., 2011). The addition of a hidden accretion column is increasingly significant for sources with a lower accretion rate, where the flux that may be hidden by the star becomes comparable to the excess observed in the UV.

Figure 11: Range of possible accretion shock models for V836 Tau. In each panel the lines are defined as in Figure 4. The left panel shows the best fit model of the accretion shock to the optical spectra. The right panel shows the fit if we assume that some accretion flux is hidden by the stellar photosphere. The upper limit on the luminosity of the low accretion column (dotted line), constrained by the lack of veiling at 1 m, is equal to that of the column which fits the NUV excess (long dashed line), doubling .

5 Correlations with Accretion Indicators

Measuring the UV excess is ideal for calculating all but the lowest s; however, UV observations are difficult to obtain. For this reason it is common to use an indicator of accretion easily accessible in the optical or near-IR which has been calibrated with s calculated for the few sources with UV observations. We assume that broad optical line emission, primarily in lines of hydrogen and calcium, originates in the material free-falling onto the star in the accretion flows supported by the fact that the emission lines are well reproduced by models which assume this geometry (Muzerolle et al., 1998, 2001). Recently, Dupree et al. (2012) suggested that the lines are instead formed in the turbulent region directly below the shock. However, Calvet (1981) found that even a chromosphere covering the entire stellar surface did not have enough emitting volume to explain the lower Balmer lines in CTTS; the shock region also has a small emitting volume and so it is unlikely to explain the observed emission in these lines.

Both the equivalent width of H, EW(H), and the width of H at 10% of the maximum flux are often used as accretion proxies (White & Basri, 2003; Natta et al., 2004). The width of H at 10% is expected to be a better indicator of accretion, since the wings of the line trace the fast moving material in the accretion flows (Muzerolle et al., 2001). Also, the EW(H) saturates for the highest accretors, or sources with the largest veiling (Muzerolle et al., 1998). Correlations between and H 10% width have been widely used; however there is considerable scatter because H and the UV excess have not been measured simultaneously. The dataset presented here, with UV and H observations separated by less than one to a few days, would be ideal for this analysis; however, the resolution of the SMARTS optical spectra is too low to accurately measure the 10% width of H. Even observing the wings of H may not be the best method for determining , as Costigan et al. (2012) showed that variability in the wings exceeded variability observed in other tracers.

In Figure 12, and EW(H) are plotted for the sample of our sources which had nearly simultaneous SMARTS optical and UV observations. There is no clear relation between EW(H) and . Two sources, DM Tau and RECX 15, have relatively low but significant EW(H), showing that it is not a good tracer of . However, Herczeg & Hillenbrand (2008) and Manara et al. (2012) found correlations between the luminosity of H and or . Herczeg & Hillenbrand (2008) determine using optical spectra with coverage down to the atmospheric limit in the blue while Manara et al. (2012) observed the excess (to estimate ) and H with the HST Wide-Field Planetary Camera 2 (WFPC2). Our analysis is the next step in defining these correlations because we have long wavelength spectral coverage, in particular the STIS spectra provide the shape of the excess in the UV. We can measure the flux of the H line in the SMARTS data and is then calculated using distances from Table 2. In Figure 13 we see a correlation between H luminosity and (or ) with a Pearson correlation coefficient of 0.9. The least-square fit to the data in Figure 13 yields


For Equations 916, is in units of and all luminosities are in units of . Equation 10 is in good agreement with the fits found in Herczeg & Hillenbrand (2008) and differs only slightly from that of Manara et al. (2012), who found a y-intercept of 2.60.1 as opposed to our value of 1.30.7.

Figure 12: versus EW(H). We compare our calculated s given in Table 5 to the EW(H) measured from the nearly simultaneous SMARTS spectra. The asterisks show the initial fits to the UV and optical excesses. Arrows show the range of for sources which were assumed to have a hidden cool accretion component, with the point of the arrow representing the upper limit on . We see no correlation between the EW(H) and .
Figure 13: and versus H line luminosity. The luminosity of H, measured in the low resolution SMARTS spectra, is correlated with the accretion properties. Asterisks and arrows are defined as in Figure 12. Least square fits to the data are given in Equations 9 and 10.

H is less often used to trace accretion than H but has been shown to correlate with (Muzerolle et al., 2001; Herczeg & Hillenbrand, 2008; Fang et al., 2009). We compare the luminosities in H () measured from the simultaneous STIS optical spectra to both and in Figure 14. We find strong correlations between the line luminosity and the accretion rates and luminosities, both with correlation coefficients 0.9. The lines shown in Figure 14 which describe the trends are given by the following equations;


Again, our relation between and agrees with that found by Herczeg & Hillenbrand (2008).

Figure 14: and versus H line luminosity. Asterisks and arrows are defined as in Figure 12. We measured the luminosity of H from the simultaneous STIS spectra. H clearly traces the accretion properties of the source. Fits to the data are given in Equations 11 and 12.

Another commonly used tracer of accretion is the Ca II infrared triplet line emission (Mohanty et al., 2005; Rigliaco et al., 2011). Our SMARTS spectra did not cover the infrared triplet lines of Ca II, but in our STIS optical data, we observe the Ca II H and K lines. The H line at 3969 Å is blended with the 3970 Å H line at the resolution of STIS, so we only consider the K line which is free of contamination. Both the H and K lines are also used as tracers of chromospheric activity (Mamajek & Hillenbrand, 2008, and references therein). We measured the luminosity of the K line in our sample of WTTS and found that the maximum value of log . For CTTS with similar Ca II K luminosities, there will be some contamination in the line due to chromospheric activity. Similar to the H and H line luminosities, the Ca II K line at 3934 Å appears to have an origin in accretion related processes, as there is a clear trend between the Ca II K luminosity () and or (with correlation coefficients of 0.8 in each case). We show the trends in Figure 15 and give the relations in the following;


Equation 14 is once again in good agreement with the relation between and produced in Herczeg & Hillenbrand (2008).

Figure 15: and versus the Ca II K line luminosity. Asterisks and arrows are defined as in Figure 12. We measured the luminosity of the Ca II K line from the simultaneous STIS spectra. Like the infrared triplet lines of Ca II, the K line luminosity correlates with the accretion properties. Fits to the data are given in Equations 13 and 14.

Included in our NUV spectra are two emission lines which have been shown to correlate with accretion properties. Calvet et al. (2004) analyzed NUV spectra for a sample of intermediate mass T Tauri stars and found that luminosities of C II] 2325 Å and Mg II 2800 Å increased with . We expected a similar trend due to simultaneous observations of the lines and accretion excess and Figure 16 shows a clear correlation. Mg II is also a tracer of chromospheric activity, and is detected in each of the WTTS templates, but the chromospheric contribution to the luminosity is likely saturated at the age of our sample (Cardini & Cassatella, 2007) and therefore not contributing to the spread in luminosity. In the CTTS, Mg II can have strong absorption features due to the presence of outflows, which was seen in a sample of NUV spectra of T Tauri stars obtained with the International Ultraviolet Explorer (Ardila et al., 2002). C II] is not observed in WTTS; however, it is easily identified even in low objects and correlates with . Therefore, C II] appears to be a very sensitive tracer of accretion or accretion related outflows (Calvet et al., 2004; Gómez de Castro & Ferro-Fontán, 2005), observable even when the NUV continuum excess is not. Relations between line fluxes and are given by;

Figure 16: Accretion luminosity versus luminosity of NUV emission lines. We compare the luminosity of C II] 2325 Å (left panel) and Mg II 2800 Å (right panel) to the range of accretion luminosities calculated from the s in Table 5. Asterisks and arrows are defined as in Figure 12. We find strong correlations between both C II] and Mg II and . Fits to the data are given in Equations 15 and 16.

In conclusion, there are several secondary tracers of accretion which are easily accessible in optical spectra and the good agreement between the correlations found here and in the literature, using different techniques, show that these correlations are robust. With well calibrated accretion indicators, it is not necessary to go to the UV to obtain accretion rate estimates. These results will facilitate the study of accretion properties for large samples of sources using optical spectra, which are significantly easier to obtain than UV data, and may offer more opportunities to measure simultaneously to data tracing other T Tauri phenomena, like outflows or circumstellar disk dynamics.

6 Summary and Conclusions

We used a large sample of T Tauri stars with nearly simultaneous spectral coverage in the FUV, NUV and optical to analyze the accretion properties of young stars. The main results are summarized here.

  • The NUV fluxes of WTTS are 3 stronger than dwarf stars of the same spectral type, due to enhanced chromospheric activity. This emission can introduce uncertainties in the estimate of the NUV excess produced by accretion if the correct templates are not used. By using WTTS templates with good signal to noise in the UV, which first became available in the DAO sample, we are able to distinguish chromospheric excesses from accretion excess when calculating accretion rates.

  • Accretion shock models which explain the NUV excess cannot describe both large filling factors predicted by models of the magnetic field distribution and non-zero 1 m veiling. Here, we modeled the shock emission as arising in columns carrying a range of energy fluxes (or densities) and covering different fractions of the stellar surface, to simulate the highly complex geometry of the accretion region. The shock emission is assumed to be a combination of high energy flux () columns which peak in the UV and low accretion columns which peak at redder, optical wavelengths. By including a variety of accretion columns, optical excesses and veiling at 1m can be explained. Also, the large filling factors of low energy columns are in better agreement with analysis based on magnetic field geometries.

  • Comparing our estimates of the accretion properties to those calculated in previous analyses Valenti et al. (1993) and Gullbring et al. (1998), we found that our assumptions for and intrinsic variability in the optical fluxes accounted for the majority of the differences in our or estimates. Two factors in our technique appear to offset each other; 1) the UV flux is decreased by using a WTTS template and 2) fitting 1 m veiling estimates increases the red shock excess.

  • Low accretion columns which can explain 1 m veiling in high accretors, may be present even when no red excess is observed. The intrinsic stellar emission, mainly the photosphere in the red and the active chromosphere in the blue, may hide this accretion component at low . We measured the maximum emission of a cool column that can be hidden by the stellar emission and determined the filling factor of the cool column. We found large surface coverage, even in sources without veiling at 1 m. In some cases the upper limit on the red excess could be equal to the excess flux in the NUV, doubling .

  • We used simultaneous measurements of emission line luminosities and accretion luminosities to calibrate correlations between the two. We found clear trends between the luminosities of H, Ca II K, Mg II and CII]. No correlation is found between EW(H) and , even with quasi-simultaneous observations, however the luminosity of H is correlated with the accretion properties. C II] may be a sensitive tracer of accretion at low , as it is observed in all the CTTS in our sample, even those with no NUV accretion excess, yet absent from the spectra of WTTS.

7 Acknowledgments

We would like to thank Melissa McClure for valuable discussions regarding extinction estimates and veiling in the near IR and providing results prior to publication. We would also like to thank Ted Bergin, Lee Hartmann, Jon Miller and Fred Adams for providing comments on an early version of this paper as part of L. Ingleby’s thesis. This work was supported by NASA grants for Guest Observer program 11616 to the University of Michigan, Caltech and the University of Colorado. Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Space Telescope Science Institute data archive. STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555. C.E. was supported by a Sagan Exoplanet Fellowship from NASA and administered by the NASA Exoplanet Science Institute. RDA acknowledges support from the Science & Technology Facilities Council (STFC) through an Advanced Fellowship (ST/G00711X/1). SGG acknowledges support from STFC via a Ernest Rutherford Fellowship (ST/J003255/1).

Object Luminosity Radius Mass SpT
mag () ()
HBC 427 0.0 0.8 1.9 0.8 K7
LkCa 19 0.0 1.7 1.6 1.3 K0
RECX 1 0.0 1.0 1.8 0.9 K5
TWA 7 0.0 0.5 1.8 0.5 M1

Note. – and SpT references: HBC 427 and LkCa 19 (Kenyon & Hartmann, 1995); RECX 1 (Luhman & Steeghs, 2004); TWA 7 (Webb et al., 1999).

Table 1: WTTS Sample and Properties
Object SpT Luminosity Radius Mass Distance
() () () (pc)
AA Tau K7 1.9 1.0 2.1 0.8 140 0.3 0.2
BP Tau K7 1.1 1.0 2.1 0.8 140 0.7 0.3
CS Cha K6 0.3 1.9 2.7 0.9 160 0.2 0.0
CV Cha G9 1.5 3.1 2.0 1.5 160 1.1 0.6
DE Tau M2 0.9 0.8 2.4 0.4 140 0.6 0.2
DK Tau A K7 1.3 1.6 2.6 0.7 140 0.5 0.5
DM Tau M1 0.7 0.2 1.1 0.5 140 0.7 0.2
DN Tau M0 0.9 1.5 2.8 0.6 140 0.1 0.0
DR Tau K5 1.4 0.4 1.1 0.9 140 8.1 2.0
FM Tau M0 0.7 0.1 0.7 0.6 140 2 0.8
GM Aur K7 0.6 1.2 2.3 0.8 140 0.2 0.0
HN Tau A K5 1.1 0.7 1.5 1.1 140 0.8 0.5
IP Tau M0 1.7 0.7 1.9 0.6 140 0.2 0.1
LkCa 15 K5 1.1 0.8 1.6 1.1 140 0.2 0.0
PDS 66 K1 0.2 0.9 1.3 1.1 86 0.2 0.1
RECX 11 K5 0.0 0.6 1.4 1.0 97 0.0 0.0
RECX 15 M3 0.0 0.1 0.9 0.3 97 0.8 0.4
RW Aur A K3 0.5 0.5 1.1 0.9 140 2.0 0.9
TWA 3a M3 0.0 0.4 1.8 0.3 50 0.0 0.0
TW Hya K7 0.0 0.3 1.1 0.8 56 0.3 0.0
V836 Tau K7 1.5 1.0 2.1 0.8 140 0.0 0.0

Note. – and SpT references: AA Tau, BP Tau, DE Tau, DN Tau, DR Tau, FM Tau, GM Aur, HN Tau A, IP Tau, LkCa 15, RW Aur A and V836 Tau (Furlan et al., 2011); DK Tau A, DM Tau, CS Cha and CV Cha (Furlan et al., 2009); RECX 11 and RECX 15 (Luhman & Steeghs, 2004); TWA 3a and TW Hya (Webb et al., 1999); PDS 66 (Mamajek et al., 2002). Distance references: Taurus Molecular Cloud (Kenyon et al., 1994); Chamaeleon I (Whittet et al., 1997); Chamaeleon (Mamajek et al., 1999); TWA (Webb et al., 1999); PDS 66 (Mamajek et al., 2002). and values were taken from Edwards et al. (2006) and Hartmann & Kenyon (1990) except for those sources identified in Section 3.2. and from the literature are in regular font, whereas those determined from new observations or found assuming a relation between and are in italics.

Table 2: CTTS Sample and Properties
Object RA DEC Telescope/ Instrument Date of Obs
(J2000) (J2000)
AA Tau 04 34 55.42 +24 28 52.8 HST/ STIS G230L/G430L 2011-01-07
CTIO/ SMARTS RC Spectrograph 2010-12-31
CTIO/ SMARTS RC Spectrograph 2011-01-02
CTIO/ SMARTS RC Spectrograph 2011-01-04
CTIO/ SMARTS RC Spectrograph 2011-01-05
BP Tau 04 19 15.86 +29 06 27.2 HST/ STIS G230L 2002-01-12
CS Cha 11 02 25.20 -77 33 36.3 HST/ STIS G230L/G430L 2011-06-01
CTIO/ SMARTS RC Spectrograph 2011-05-30
CTIO/ SMARTS RC Spectrograph 2011-05-31
CTIO/ SMARTS RC Spectrograph 2011-06-02
Magellan/ MIKE 2012-02-16
VLT/ CRIRES 2010-12-23
CV Cha 11 12 27.65 -76 44 22.1 HST/STIS E230M/G430L 2011-04-13
CTIO/ SMARTS RC Spectrograph 2011-04-08
CTIO/ SMARTS RC Spectrograph 2011-04-13
CTIO/ SMARTS RC Spectrograph 2011-04-15
CTIO/ SMARTS RC Spectrograph 2011-04-17
CTIO/ Bench-Mounted Echelle Spectrograph 2011-04-12
Magellan/ MIKE 2010-03-10
DE Tau 04 21 55.69 +27 55 06.1 HST/ STIS G230L/G430L 2010-08-20
CTIO/ SMARTS RC Spectrograph 2010-08-15
CTIO/ SMARTS RC Spectrograph 2010-08-18
CTIO/ SMARTS RC Spectrograph 2010-08-20
CTIO/ SMARTS RC Spectrograph 2010-08-21
DK Tau A 04 30 44.25 +26 01 24.5 HST/ STIS G230L/G430L 2010-02-04
DM Tau 04 33 48.74 +18 10 09.7 HST/ STIS G230L/G430L 2010-08-22
CTIO/ SMARTS RC Spectrograph 2010-08-15
CTIO/ SMARTS RC Spectrograph 2010-08-18
CTIO/ SMARTS RC Spectrograph 2010-08-20
CTIO/ SMARTS RC Spectrograph 2010-08-21
Magellan/ MIKE 2011-01-04
VLT/ CRIRES 2011-08-11
DN Tau 04 35 27.44 +24 14 59.1 HST/ STIS G230L/G430L 2011-09-10
CTIO/ SMARTS RC Spectrograph 2011-09-07
CTIO/ SMARTS RC Spectrograph 2011-09-11
DR Tau 04 47 06.22 +16 58 42.6 HST/ STIS G230L/G430L 2010-02-15
FM Tau 04 14 13.56 +28 12 48.8 HST/ STIS G230L/G430L 2011-09-21
CTIO/ SMARTS RC Spectrograph 2011-09-21
CTIO/ SMARTS RC Spectrograph 2011-09-27
VLT/ CRIRES 2010-12-02
GM Aur 04 55 10.98 +30 21 59.1 HST/ STIS G230L/G430L 2010-08-19
HBC 427 04 56 02.02 +30 21 03.2 HST/ STIS G230L/G430L 2011-03-30
HN Tau A 04 33 39.37 +17 51 52.1 HST/ STIS G230L/G430L 2010-02-10
IP Tau 04 24 57.14 +27 11 56.4 HST/ STIS G230L/G430L 2011-03-21
VLT/ CRIRES 2010-11-13
LkCa 15 04 39 17.73 +22 21 03.8 HST/ STIS G230L 2003-02-13
Magellan/ MIKE 2011-01-04
VLT/ CRIRES 2010-02-03
VLT/ CRIRES 2010-02-05
LkCa 19 04 55 36.97 +30 17 55.0 HST/ STIS G230L/G430L 2011-03-31
PDS 66 13 22 07.45 -69 38 12.6 HST/ STIS G230L/G430L 2011-05-23
CTIO/ SMARTS RC Spectrograph 2011-05-16
CTIO/ SMARTS RC Spectrograph 2011-05-20
CTIO/ Bench-Mounted Echelle Spectrograph 2011-05-22
RECX 1 08 36 56.12 -78 56 45.3 HST/ STIS G230L/G430L 2010-01-22
VLT/ CRIRES 2011-05-30
RECX 11 08 47 01.28 -78 59 34.1 HST/ STIS G230L/G430L 2009-12-12
CTIO/ SMARTS RC Spectrograph 2009-11-27
CTIO/ SMARTS RC Spectrograph 2009-12-15
CTIO/ SMARTS RC Spectrograph 2009-12-19
Magellan/ MIKE 2010-03-10
VLT/ CRIRES 2011-05-30
RECX 15 08 43 18.43 -79 05 17.7 HST/ STIS G230L/G430L 2010-02-05
CTIO/ SMARTS RC Spectrograph 2010-02-01
CTIO/ SMARTS RC Spectrograph 2010-02-04
CTIO/ SMARTS RC Spectrograph 2010-02-06
Magellan/ MIKE 2010-03-11
RW Aur A 05 07 49.51 +30 24 04.8 HST/ STIS G230L/G430L 2011-03-25
TWA 3a 11 10 27.80 -37 31 51.2 HST/ STIS G230L/G430L 2011-03-26
CTIO/ SMARTS RC Spectrograph 2011-03-18
CTIO/ SMARTS RC Spectrograph 2011-03-22
CTIO/ SMARTS RC Spectrograph 2011-03-25
CTIO/ SMARTS RC Spectrograph 2011-03-29
Magellan/ MIKE 2010-03-10
TWA 7 10 42 29.94 -33 40 16.7 HST/ STIS G230L/G430L 2011-05-05
TW Hya 11 01 51.95 -34 42 17.7 HST/ STIS G230 2002-05-10
V836 Tau 05 03 06.62 +25 23 19.6 HST/ STIS G230L/G430L 2011-02-05
CTIO/ SMARTS RC Spectrograph 2011-01-28
CTIO/ SMARTS RC Spectrograph 2011-02-01
CTIO/ SMARTS RC Spectrograph 2011-02-03
CTIO/ SMARTS RC Spectrograph 2011-02-08
Table 3: Log of Observations
Object F09/11 KH95 G98 V93 F11 M13
AA Tau 1.9 0.5 0.7 1.3 1.3
BP Tau 1.1 0.5 0.5 0.9 1.8 0.6
DE Tau 0.9 0.6 0.6 1.7 0.9
DK Tau A 1.3 0.8 1.4 1.2 1.8
DM Tau 0.7 0.0 0.1
DN Tau 0.9 0.5 0.3 0.5
DR Tau 1.4 1.0 1.8 2.0
FM Tau 0.7 0.7 0.8
GM Aur 0.6 0.1 0.3 0.5
HN Tau A 1.1 0.5 0.7 1.0 3.1
IP Tau 1.7 0.2 0.3
LkCa 15 1.1 0.6
RW Aur A 0.5 1.2
V836 Tau 1.5 0.6 1.4

Note. – F09/11 (Furlan et al., 2009, 2011), KH95 (Kenyon & Hartmann, 1995), G98 (Gullbring et al., 1998), V93 (Valenti et al., 1993), F11 (Fischer et al., 2011), M13 (McClure et al., in prep)

Table 4: Literature s for Taurus Sources
AA Tau 0 0 0 0 0.002 0.002 0.0 1.5
BP Tau 0 0 0.02 0 0.002 0.022 0.2 2.9
CS Cha 0.008 0 0.003 0 0 0.011 0.0 5.3
CV Cha 0 0.4 0 0.02 0 0.42 0.6 5.9
DE Tau 0.1 0 0.0008 0.0007 0 0.1 0.3 2.8
DK Tau A 0.01 0 0.005 0.005 0 0.02 0.1 3.4
DM Tau 0.08 0 0 0.001 0.0006 0.082 0.3 2.9
DN Tau 0 0 0 0.002 0 0.002 0.0 1.0
DR Tau 0 0 0.3 0.07 0.006 0.37 1.6 5.2
FM Tau 0 0.07 0 0 0.001 0.071 0.8 1.2
GM Aur 0 0 0.001 0.003 0 0.004 0.0 9.6
HN Tau A 0 0 0 0.01 0.004 0.014 0.1 1.4
IP Tau 0 0 0 0 0.001 0.001 0.0 7.2
LkCa 15 0 0 0.01 0.0007 0.0001 0.011 0.1 3.1
PDS 66 0 0 0 0 0.0001 0.0001 0.0 1.3
RECX 11 0 0 0.001 0 0 0.001 0.0 1.7
RECX 15 0.02 0.001 0 0.0007 8 0.022 0.3 8.0
RW Aur A 0 0.03 0.2 0 0 0.23 0.8 2.0
TWA 3a 0 0 0 0 8 8 0.0 9.8
TW Hya 0 0 0 0.0013 0.0013 0.0026 0.0 1.8
V836 Tau 0 0 0 0 0.0002 0.0002 0.0 1.1

Note. – The values in parentheses represents the energy flux () of each column in erg s cm.

Table 5: Results from Multi-Component Model Fits to UV and Optical Spectra
Object (G98) (V93)
() ()
AA Tau 3.3 7.1
BP Tau 2.9 2.4
DE Tau 2.6 1.8
DK Tau A 3.8 6.1
DM Tau 2.9
DN Tau 3.5 1.3
FM Tau 7.7
GM Aur 9.6 7.4
HN Tau A 1.3 3.9
IP Tau 8.0
RW Aur A 3.3

Note. –
Gullbring et al. (1998)
Valenti et al. (1993)

Table 6: Accretion Rates from Literature for Taurus Objects
AA Tau 1.5 0.002 1.6 0.08
BP Tau 2.9 0.022 3.3 0.2
DN Tau 1.0 0.002 1.7 0.06
DR Tau 5.2 0.37 5.6 0.8
GM Aur 9.6 0.004 1.3 0.04
HN Tau A 1.4 0.014 1.5 0.06
IP Tau 7.2 0.001 9.4 0.03
LkCa 15 3.1 0.011 3.6 0.03
PDS 66 1.3 0.0001 3.8 0.02
RECX 11 1.7 0.001 2.6 0.006
TWA 3a 9.8 3.4 0.002
TW Hya 1.8 0.0026 2.3 0.05
V836 Tau 1.1 0.0002 2.6 0.02

Note. –
Values which give the best fit of the model to the UV and optical data (see Table 5).
Values calculated when allowing for a low column which produces emission not detectable above the intrinsic stellar emission.

Table 7: Maximum and Assuming Hidden Accretion Emission


  1. affiliation: Department of Astronomy, University of Michigan, 830 Dennison Building, 500 Church Street, Ann Arbor, MI 48109, USA;,,
  2. affiliation: Department of Astronomy, University of Michigan, 830 Dennison Building, 500 Church Street, Ann Arbor, MI 48109, USA;,,
  3. affiliation: The Kavli Institute for Astronomy and Astrophysics, Peking University, Yi He Yuan Lu 5, Hai Dian Qu, 100871 Beijing, P. R. China
  4. affiliation: Department of Astronomy, University of Michigan, 830 Dennison Building, 500 Church Street, Ann Arbor, MI 48109, USA;,,
  5. affiliation: Stony Brook University, Stony Brook NY 11794-3800, USA
  6. affiliation: NASA Herschel Science Center, California Institute of Technology, Mail Code 100-22, Pasadena, CA 91125, USA
  7. affiliation: Department of Physics & Astronomy, University of Leicester, University Road, Leicester, LE1 7RH, UK
  8. affiliation: Department of Astronomy, Smith College, Northampton, MA 01063, USA
  9. affiliation: NASA Sagan Postdoctoral Fellow. Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-78, Cambridge, MA 02138, USA
  10. affiliation: California Institute of Technology, Department of Astrophysics, MC 249-17, Pasadena, CA 91125, USA
  11. affiliationmark:
  12. affiliation: School of Physics & Astronomy, University of St. Andrews, St. Andrews, KY16 9SS, UK
  13. affiliation: California Institute of Technology, Department of Astrophysics, MC 249-17, Pasadena, CA 91125, USA
  14. affiliation: Center for Astrophysics and Space Astronomy, University of Colorado, Boulder, CO 80309-0389, USA


  1. Alencar, S. H. P., Bouvier, J., Walter, F. M., et al. 2012, A&A, 541, A116
  2. Ardila, D., Herczeg, G., Gregory, S., et al. in press, ApJ
  3. Ardila, D. R., Basri, G., Walter, F. M., Valenti, J. A., & Johns-Krull, C. M. 2002, ApJ, 566, 1100
  4. Basri, G., & Batalha, C. 1990, ApJ, 363, 654
  5. Batalha, C. C., Quast, G. R., Torres, C. A. O., et al. 1998, A&AS, 128, 561
  6. Bergin, E., Calvet, N., Sitko, M. L., et al. 2004, ApJ, 614, L133
  7. Bertout, C. 1989, ARA&A, 27, 351
  8. Bouvier, J., Covino, E., Kovo, O., et al. 1995, A&A, 299, 89
  9. Calvet, N. 1981, Rev. Mexicana Astron. Astrofis., 6, 227
  10. Calvet, N., Briceño, C., Hernández, J., et al. 2005, AJ, 129, 935
  11. Calvet, N., & Gullbring, E. 1998, ApJ, 509, 802
  12. Calvet, N., Muzerolle, J., Briceño, C., et al. 2004, AJ, 128, 1294
  13. Calvet, N., Patino, A., Magris, G. C., & D’Alessio, P. 1991, ApJ, 380, 617
  14. Cardini, D., & Cassatella, A. 2007, ApJ, 666, 393
  15. Carpenter, J. M., Hillenbrand, L. A., & Skrutskie, M. F. 2001, AJ, 121, 3160
  16. Carpenter, J. M., Hillenbrand, L. A., Skrutskie, M. F., & Meyer, M. R. 2002, AJ, 124, 1001
  17. Cortes, S. R., Meyer, M. R., Carpenter, J. M., et al. 2009, ApJ, 697, 1305
  18. Costigan, G., Scholz, A., Stelzer, B., et al. 2012, ArXiv e-prints
  19. D’Alessio, P., Calvet, N., Hartmann, L., Lizano, S., & Cantó, J. 1999, ApJ, 527, 893
  20. Dodin, A. V., & Lamzin, S. A. 2012, Astronomy Letters, 38, 649
  21. Dodin, A. V., Lamzin, S. A., & Sitnova, T. M. 2013, arXiv:1302.1825
  22. Donati, J.-F., Jardine, M. M., Gregory, S. G., et al. 2008, MNRAS, 386, 1234
  23. Donati, J.-F., Bouvier, J., Walter, F. M., et al. 2011, MNRAS, 412, 2454
  24. Dupree, A. K., Brickhouse, N. S., Cranmer, S. R., et al. 2012, ApJ, 750, 73
  25. Edwards, S., Fischer, W., Hillenbrand, L., & Kwan, J. 2006, ApJ, 646, 319
  26. Edwards, S., Hartigan, P., Ghandour, L., & Andrulis, C. 1994, AJ, 108, 1056
  27. Espaillat, C., Calvet, N., D’Alessio, P., et al. 2007a, ApJ, 670, L135
  28. —. 2007b, ApJ, 664, L111
  29. Espaillat, C., D’Alessio, P., Hernández, J., et al. 2010, ApJ, 717, 441
  30. Fang, M., van Boekel, R., Wang, W., et al. 2009, A&A, 504, 461
  31. Findeisen, K., Hillenbrand, L., & Soderblom, D. 2011, AJ, 142, 23
  32. Fischer, W., Edwards, S., Hillenbrand, L., & Kwan, J. 2011, ApJ, 730, 73
  33. France, K., Schindhelm, E., Burgh, E. B., et al. 2011, ApJ, 734, 31
  34. Furlan, E., Watson, D. M., McClure, M. K., et al. 2009, ApJ, 703, 1964
  35. Furlan, E., Luhman, K. L., Espaillat, C., et al. 2011, ApJS, 195, 3
  36. Gahm, G. F., Walter, F. M., Stempels, H. C., Petrov, P. P., & Herczeg, G. J. 2008, A&A, 482, L35
  37. Gálvez, M. C., Montes, D., Fernández-Figueroa, M. J., De Castro, E., & Cornide, M. 2009, AJ, 137, 3965
  38. Gauvin, L. S., & Strom, K. M. 1992, ApJ, 385, 217
  39. Gómez de Castro, A. I., & Ferro-Fontán, C. 2005, MNRAS, 362, 569
  40. Grankin, K. N., Melnikov, S. Y., Bouvier, J., Herbst, W., & Shevchenko, V. S. 2007, A&A, 461, 183
  41. Gregorio-Hetem, J., Lepine, J. R. D., Quast, G. R., Torres, C. A. O., & de La Reza, R. 1992, AJ, 103, 549
  42. Gregory, S. G., & Donati, J.-F. 2011, Astronomische Nachrichten, 332, 1027
  43. Gregory, S. G., Donati, J.-F., Morin, J., et al. 2012, ApJ, 755, 97
  44. Guenther, E. W., Esposito, M., Mundt, R., et al. 2007, A&A, 467, 1147
  45. Guinan, E. F., Ribas, I., & Harper, G. M. 2003, ApJ, 594, 561
  46. Gullbring, E., Calvet, N., Muzerolle, J., & Hartmann, L. 2000, ApJ, 544, 927
  47. Gullbring, E., Hartmann, L., Briceno, C., & Calvet, N. 1998, ApJ, 492, 323
  48. Hartigan, P., Hartmann, L., Kenyon, S., Hewett, R., & Stauffer, J. 1989, ApJS, 70, 899
  49. Hartigan, P., Kenyon, S. J., Hartmann, L., et al. 1991, ApJ, 382, 617
  50. Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385
  51. Hartmann, L., Hewett, R., & Calvet, N. 1994, ApJ, 426, 669
  52. Hartmann, L. W., & Kenyon, S. J. 1990, ApJ, 349, 190
  53. Heap, S. R., & Lindler, D. J. 2007, in Astronomical Society of the Pacific Conference Series, Vol. 374, From Stars to Galaxies: Building the Pieces to Build Up the Universe, ed. A. Vallenari, R. Tantalo, L. Portinari, & A. Moretti, 409
  54. Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 1994, AJ, 108, 1906
  55. Herczeg, G. in prep
  56. Herczeg, G. J., & Hillenbrand, L. A. 2008, ApJ, 681, 594
  57. Herczeg, G. J., Walter, F. M., Linsky, J. L., et al. 2005, AJ, 129, 2777
  58. Houdebine, E. R., Mathioudakis, M., Doyle, J. G., & Foing, B. H. 1996, A&A, 305, 209
  59. Ingleby, L., Calvet, N., Bergin, E., et al. 2009, ApJ, 703, L137
  60. —. 2011, ApJ, 743, 105
  61. Johns-Krull, C. M., & Gafford, A. D. 2002, ApJ, 573, 685
  62. Johns-Krull, C. M., Valenti, J. A., Hatzes, A. P., & Kanaan, A. 1999, ApJ, 510, L41
  63. Johns-Krull, C. M., Valenti, J. A., & Linsky, J. L. 2000, ApJ, 539, 815
  64. Kenyon, S. J., Dobrzycka, D., & Hartmann, L. 1994, AJ, 108, 1872
  65. Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117
  66. Kraus, A. L., Ireland, M. J., Martinache, F., & Hillenbrand, L. A. 2011, ApJ, 731, 8
  67. Kurucz, R. L., & Bell, B. 1995, Atomic line list
  68. Lawson, W. A., Feigelson, E. D., & Huenemoerder, D. P. 1996, MNRAS, 280, 1071
  69. Lawson, W. A., Lyo, A.-R., & Muzerolle, J. 2004, MNRAS, 351, L39
  70. Long, M., Romanova, M. M., Kulkarni, A. K., & Donati, J.-F. 2011, MNRAS, 413, 1061
  71. Luhman, K. L., & Steeghs, D. 2004, ApJ, 609, 917
  72. Mamajek, E. E., & Hillenbrand, L. A. 2008, ApJ, 687, 1264
  73. Mamajek, E. E., Lawson, W. A., & Feigelson, E. D. 1999, ApJ, 516, L77
  74. Mamajek, E. E., Meyer, M. R., & Liebert, J. 2002, AJ, 124, 1670
  75. Manara, C. F., Robberto, M., Da Rio, N., et al. 2012, ApJ, 755, 154
  76. McClure, M., Calvet, N., Espaillat, C., et al. in prep
  77. Mohanty, S., Jayawardhana, R., & Basri, G. 2005, ApJ, 626, 498
  78. Muzerolle, J., Calvet, N., & Hartmann, L. 1998, ApJ, 492, 743
  79. —. 2001, ApJ, 550, 944
  80. Nagel, E., Espaillat, C., D’Alessio, P., & Calvet, N. 2012, ApJ, 747, 139
  81. Natta, A., Testi, L., Muzerolle, J., et al. 2004, A&A, 424, 603
  82. Nguyen, D. C., Brandeker, A., van Kerkwijk, M. H., & Jayawardhana, R. 2012, ApJ, 745, 119
  83. Pecaut, M. J., Mamajek, E. E., & Bubar, E. J. 2012, ApJ, 746, 154
  84. Petrov, P. P., Gahm, G. F., Stempels, H. C., Walter, F. M., & Artemenko, S. A. 2011, A&A, 535, A6
  85. Rigliaco, E., Natta, A., Randich, S., Testi, L., & Biazzo, K. 2011, A&A, 525, A47
  86. Rigliaco, E., Natta, A., Testi, L., et al. 2012, ArXiv e-prints
  87. Romaniello, M., Robberto, M., & Panagia, N. 2004, ApJ, 608, 220
  88. Sicilia-Aguilar, A., Henning, T., & Hartmann, L. W. 2010, ApJ, 710, 597
  89. Sicilia-Aguilar, A., Bouwman, J., Juhász, A., et al. 2009, ApJ, 701, 1188
  90. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593
  91. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
  92. Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 97, 1451
  93. Takami, M., Bailey, J., & Chrysostomou, A. 2003, A&A, 397, 675
  94. Tody, D. 1993, in Astronomical Society of the Pacific Conference Series, Vol. 52, Astronomical Data Analysis Software and Systems II, ed. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, 173
  95. Valenti, J. A., Basri, G., & Johns, C. M. 1993, AJ, 106, 2024
  96. Webb, R. A., Zuckerman, B., Platais, I., et al. 1999, ApJ, 512, L63
  97. White, R. J., & Basri, G. 2003, ApJ, 582, 1109
  98. White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265
  99. White, R. J., & Hillenbrand, L. A. 2004, ApJ, 616, 998
  100. Whittet, D. C. B., Prusti, T., Franco, G. A. P., et al. 1997, A&A, 327, 1194
  101. Whittet, D. C. B., Shenoy, S. S., Clayton, G. C., & Gordon, K. D. 2004, ApJ, 602, 291
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description