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
SMART_READER_LITE
LIVE PREVIEW

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


slide-1
SLIDE 1

IGRINS Survey of T Tauri Stars in Tau-Aur

Status Report

Kyoung Hee Kim (KHU)

2017 IGRINS Users Workshop 2017.07.28

slide-2
SLIDE 2

Star Formation & Young Stellar Objects

Overview of star formation (Greene, 2001)

slide-3
SLIDE 3

Dullemond+2010

slide-4
SLIDE 4

Borrowed from Elisabetta Rigliaco’s homepage (orginally from Hartmann (2009))

Schematic picture of accretion in T Tauri stars

slide-5
SLIDE 5

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.

slide-6
SLIDE 6

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)

slide-7
SLIDE 7

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

slide-8
SLIDE 8

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)

slide-9
SLIDE 9

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

slide-10
SLIDE 10

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

slide-11
SLIDE 11

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

slide-12
SLIDE 12

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

slide-13
SLIDE 13

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

slide-14
SLIDE 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)

slide-15
SLIDE 15

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

slide-16
SLIDE 16
slide-17
SLIDE 17

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.