IGRINS Survey of T Tauri Stars in Tau-Aur Status Report Kyoung Hee - - PowerPoint PPT Presentation
IGRINS Survey of T Tauri Stars in Tau-Aur Status Report Kyoung Hee - - PowerPoint PPT Presentation
IGRINS Survey of T Tauri Stars in Tau-Aur Status Report Kyoung Hee Kim (KHU) 2017 IGRINS Users Workshop 2017.07.28 Star Formation & Young Stellar Objects Overview of star formation (Greene, 2001) Dullemond+2010 Schematic picture of
Star Formation & Young Stellar Objects
Overview of star formation (Greene, 2001)
Dullemond+2010
Borrowed from Elisabetta Rigliaco’s homepage (orginally from Hartmann (2009))
Schematic picture of accretion in T Tauri stars
Purpose of the program
- Goal: to understand the disk evolution and the effect of
planet formation in the accretion system via systematic studies using a large sample
- Method: We measure properties related to mass accretion
(e.g., veiling, mass accretion rate, radial velocity, magnetic field strength, other accretion tracers).
(1) Accretion rates from Brackett lines. (2) Any trends of the disk properties (e.g., spectral index). (3) Line profiles of Br gamma. (4) Search correlations among other lines. (5) Investigate methods to figure out mass accretion activity and magnetic fields (use Ti I lines). (6) Eventually test magnetospheric accretion models.
Observation Status
- Proposed: 126 Class II YSOs selected (Kmag < 11) + α
- Observations: from Oct. 2014 to Jan. 2015 (UT/KASI/KHU)
- 86 objects observed out of 126
– Single epoch: 56 – 2 epochs: 23 – 3 epochs: 3 – 3+ : 4 (CI Tau, DS Tau, DH Tau A, DO Tau)
(Updated since IGRINS 2015 user meeting)
Reduction process (2015)
- Spectrum Extraction with PLP_v2.1
- IRAF: scombine, splot, telluric
- Extract K band 12th order (2.155~2.185µm)
- Correct Br γ absorption of A0 Standard star’s
spectrum
- Telluric correction of target spectrum (adjust scale
factor and wavelength shift)
- Normalize the telluric corrected target spectrum
- Measure equivalent width (EW) of Br γ emission
feature
From Br γ to dM/dt
- EW = Fline/Fcont
- Adopted 2MASS K band magnitude as Fcont
- Distance to Taurus, d=140 pc
- LBrγ = 4πd2Fline
(Muzerolle+1998)
Reduction Status
- Reduction done by PLP 2.1-alpha4 (2015) and PLP 2.2-alpha1 (2017)
- Keep Brγ analysis with the products after IRAF tasks (plp 2.1-alpha4)
- 64: available dM/dt right now
– Adopted averaged dM/dt for multi epochs data for the current purpose – Median(dM/dt): 10-7.7 Msun/yr (2.1 x 10-8 Msun/yr)
- others: too faint, very weak emission or very noisy, can’t fit; no
available info of R*, M*, etc by now
Mstar, Lstar, Rstar, Teff, and Av from Andrews+2013 (many from Furlan+2011) Mdisk from Andrews+(2005) (31/64 available)
Reproduction of the basic correlations
correlated correlated no correlation! Disk masses to be updated w. recent
- bservations
(ALMA)
N=64 s=1.3+-0.3 corr=0.4+-0.1 N=64 s=0.9+-0.1 corr=0.8+-0.1 N=64 s=-0.1+-0.2 corr=-0.1+-0.2
n(K-6), n(13-31), and EW(10) from Spitzer/IRS data
correlated tendency no corr no corr Correlation No noticeable trends between mass accretion rates and properties of
- uter disks
N=64 s=0.14+-0.03 corr=0.5+-0.1 N=58 s=0.4+-0.2 corr=0.4+-0.1 N=58 s=-0.1+-0.1 corr=0.1+-0.1
2MASS J, H, and K band magnitude (extinction corrected)
correlated correlated correlated The strength of the correlation increased from J to K. Why? Is this related to disk structures (gas and dust)? Could RADMC-3 continuum modeling help to understand it?
N=63 s=-0.5+-0.1 corr=-0.83+-0.1 N=63 s=-0.6+-0.1 corr=-0.94+-0.06 N=63 s=-0.6+-0.1 corr=-0.97+-0.1
spectral index of 2MASS bands (extinction corrected)
It looks scattered, but focus on distribution between -1 and 1. Stellar contribution to J? Disk contribution to K (and H)?
N=63 s=0.1+-0.1 corr=0.12+-0.13 N=63 s=0.2+-0.1 corr=0.2+-0.13 N=63 s=0.2+-0.1 corr=0.24+-0.14
Next: Examine Spectra
H2 2.119 2.120 2.121 2.122 2.123 2.124 0.90 0.95 1.00 1.05 1.10 Brc 2.162 2.164 2.166 2.168 2.170 0.9 1.0 1.1 1.2 Ti I 2.220 2.222 2.224 2.226 0.85 0.90 0.95 1.00 1.05 1.10 1.15 2.300 2.305 2.310 2.315 0.8 0.9 1.0 1.1 1.2 2.360 2.365 2.370 2.375 2.380 0.7 0.8 0.9 1.0 1.1 1.2 1.3
dhtau (20141130)
H2 2.119 2.120 2.121 2.122 2.123 2.124 0.90 0.95 1.00 1.05 1.10 Brc 2.162 2.164 2.166 2.168 2.170 0.8 1.0 1.2 1.4 1.6 Ti I 2.220 2.222 2.224 2.226 0.85 0.90 0.95 1.00 1.05 1.10 1.15 2.300 2.305 2.310 2.315 0.9 1.0 1.1 1.2 2.360 2.365 2.370 2.375 2.380 0.7 0.8 0.9 1.0 1.1 1.2 1.3
fztau (20141124)
ex) To figure out how to analyze dM/dt and other lines
H2 2.119 2.120 2.121 2.122 2.123 2.124 0.9 1.0 1.1 1.2 Brc 2.162 2.164 2.166 2.168 2.170 0.6 0.8 1.0 1.2 1.4 Ti I 2.220 2.222 2.224 2.226 0.8 0.9 1.0 1.1 2.300 2.305 2.310 2.315 0.8 0.9 1.0 1.1 2.360 2.365 2.370 2.375 2.380 0.0 0.5 1.0 1.5 2.0 2.5 3.0
gntau (20150126)
the magnetospheric accretions in a stable (left-hand panel) and an unstable (right-hand panel) regime
Kurosawa & Romanova 2013
Future Goal: Test Magnetospheric accretion models
Summary
- There has been no progress for last one and half year.
- But, we have resumed IGRINS data analysis this summer.
- It is early to make conclusions based on the current correlation analysis.
- We need to refine correlations by dealing the variability and multiple system
effects.
- Utilizing other lines and modeling to trace and confine mass accretion
features is necessary.
- Radial velocity analysis will be followed soon.
- We may develop an idea for a Gemini proposal.