Slide 1

Slide 1 text

Physical properties of starspots and their effects on young stellar ages Michael Gully-Santiago Kepler/K2 Guest Observer Office Thursday, June 22, 2017 at NASA Ames Kepler/K2 Science Conference Collaborators: Greg Herczeg (KIAA), Ian Czekala (CfA/KICAP), Garrett Somers (OSU/Vanderbilt), J.F. Donati (CNRS), Konstantin Grankin (CrAO), Kevin Covey (WWU), G. Mace (UTexas), MATYSSE team, ASASSN team, ++

Slide 2

Slide 2 text

Disk fraction as a function of age gives us the disk lifetime. Haisch et al. 2001

Slide 3

Slide 3 text

Disk fraction as a function of age gives us the disk lifetime. The disk lifetime limits the planet formation timescale. Haisch et al. 2001, Alexander & Armitage 2009, Kraus et al. 2012

Slide 4

Slide 4 text

The disk lifetime limits the planet formation timescale. Haisch et al. 2001, Alexander & Armitage 2009, Kraus et al. 2012 How accurate are cluster ages?

Slide 5

Slide 5 text

How accurate are cluster ages? ? ? ? ?

Slide 6

Slide 6 text

We want to know the ages and masses of young stars so that we can understand how star and planet formation proceed in time.

Slide 7

Slide 7 text

We want to know the ages and masses of young stars so that we can understand how star and planet formation proceed in time. Age is not a direct observable. Mass is sometimes a direct observable, but only in rare eclipsing binary systems or resolved gas disk systems.

Slide 8

Slide 8 text

We want to know the ages and masses of young stars so that we can understand how star and planet formation proceed in time. Age is not a direct observable. Mass is sometimes a direct observable, but only in rare eclipsing binary systems or resolved gas disk systems.

Slide 9

Slide 9 text

Ages are estimated by placing a young star on a pre-main sequence HR diagram.

Slide 10

Slide 10 text

Ages are estimated by placing a young star on a pre-main sequence HR diagram.

Slide 11

Slide 11 text

Ages are estimated by placing a young star on a pre-main sequence HR diagram.

Slide 12

Slide 12 text

Ages are estimated by placing a young star on a pre-main sequence HR diagram.

Slide 13

Slide 13 text

Ages are estimated by placing a young star on a pre-main sequence HR diagram.

Slide 14

Slide 14 text

Ages are estimated by placing a young star on a pre-main sequence HR diagram.

Slide 15

Slide 15 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why?

Slide 16

Slide 16 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots

Slide 17

Slide 17 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots Herczeg & Hillenbrand 2014 Spectral types map imprecisely to Teff

Slide 18

Slide 18 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots APOGEE spectra of thousands of young stars show large disagreement Cottaar et al. 2014

Slide 19

Slide 19 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots

Slide 20

Slide 20 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots 290 300 310 320 330 340 350 360 l 10 5 0 5 10 15 20 25 30 b ⇢ Oph Lup I Lup II Lup III Lup IV F/G/K/M Pre-MS Accretors Turnoff B-type Other B1-B5 24 Myr 20 Myr 16 Myr 13 Myr 11 Myr 9 Myr 7 Myr Age Pecaut & Mamajek 2016

Slide 21

Slide 21 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots

Slide 22

Slide 22 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots

Slide 23

Slide 23 text

Observed star clusters (~1 Myr) show large age spreads of 1-10 Myr. Why? - Observational uncertainties - True age spreads - Episodic accretion - Physics beyond the standard evolutionary models - Magnetic fields - Starspots

Slide 24

Slide 24 text

No content

Slide 25

Slide 25 text

No content

Slide 26

Slide 26 text

Conspicuous starspot transits have been relatively rare. Little/no enhanced in-transit variance has been observed. - B. Morris & E. Krause from discussions at this conference

Slide 27

Slide 27 text

Starspots imbue periodic flux diminution as they enter-and-exit the unresolved stellar disk. Spots have been taken for granted for measuring rotation periods. Spot amplitudes encode information about longitudinally asymmetric spots. Notsu et al. 2013

Slide 28

Slide 28 text

Starspots have both theoretical and observational consequences.

Slide 29

Slide 29 text

Starspots inhibit convective efficiency. Somers & Pinnsoneault 2015 Less efficient energy transport means stars get larger (R increases) but cooler (T decreases) Postdoc at Vanderbilt

Slide 30

Slide 30 text

Starspots confound measurements of L and Teff. Tamb=Teff 0th order assumption No starspots Tamb

Slide 31

Slide 31 text

Starspots confound measurements of L and Teff. Tamb=Teff 0th order assumption No starspots 1st order correction Non-emitting starspots Tamb Tspot = 0 K fspot Tamb

Slide 32

Slide 32 text

Starspots confound measurements of L and Teff. Tamb=Teff 0th order assumption No starspots 1st order correction Non-emitting starspots Tamb Tspot > 0 K fspot Tamb Tspot = 0 K fspot Tamb 2nd order correction Emitting starspots

Slide 33

Slide 33 text

Starspots confound measurements of L and Teff. Tamb=Teff 0th order assumption No starspots 1st order correction Non-emitting starspots Tamb Tspot > 0 K fspot Tamb Tspot = 0 K fspot Tamb 2nd order correction Emitting starspots Teff < Tamb

Slide 34

Slide 34 text

Photometric modulation probes longitudinally asymmetric spots. ΔV LkCa 4 ΔV 2015: 0.5 2004: 0.8 1986: 0.2 LkCa 4 light curve

Slide 35

Slide 35 text

You can find a minimum coverage of starspots for LkCa 4 ΔV LkCa 4 ΔV 2015: 0.5 2004: 0.8 1986: 0.2 Spot-free Spotted, non-emitting spots

Slide 36

Slide 36 text

You can find a minimum coverage of starspots for LkCa 4 ΔV LkCa 4 ΔV 2015: 0.5 2004: 0.8 1986: 0.2 Spot-free Spotted, emitting spots Tspot

Slide 37

Slide 37 text

Expected range of starspot properties

Slide 38

Slide 38 text

"any  photometric  modula2ons  due  to  starspots  is  the  asymmetric   component  of  the  starspot  coverage"     Neff  et  al.  1995,  Harrison  et  al  2011 Smith  1994 Uniformly distributed spots; pole-on star spots circumpolar spots Differential rotation confounding periods Geometrical effects hinder interpretation of lightcurve amplitudes.

Slide 39

Slide 39 text

K2 Campaign 2 Observed the ~few Myr Oph and Upper Sco. What can be learned from periodic lightcurve amplitudes in K2?

Slide 40

Slide 40 text

Young stars have larger amplitude of periodic variability than typical stars (we already knew that). K2  Cycle  2  light  curves  for     1658  candidate  or  confirmed  young   stars  towards  Oph/Sco.   compared  to     everything  else  in  that  Cycle. Uncorrected Vanderburg and Johnson lightcurves

Slide 41

Slide 41 text

Caveats: non-members non-periodic structures Vanderburg and Johnson lightcurves Periods < 10 days K2 C2 lightcurve amplitudes: ~0.5-5% peak to valley

Slide 42

Slide 42 text

How to resolve geometric degeneracies to assess total starspot coverage fraction?

Slide 43

Slide 43 text

Directly detect the spectrum arising from starspots. Tamb Tspot fspot 1-fspot My Research - This has historically been challenging due to low flux ratio of spots. - My strategy leverages new high resolution near-IR spectrographs.

Slide 44

Slide 44 text

Starspot emission Tamb Tspot fspot 1-fspot Tspot = 2800 K Tamb = 4100 K fspot = 0.5 (!) Example Key insight: - In the visible, starspot flux is 5-20x weaker than the ambient photosphere. - In the near-IR, starspot flux is only 2.5-4x weaker than ambient

Slide 45

Slide 45 text

Starspot emission Tamb Tspot fspot 1-fspot Tspot = 2800 K Tamb = 4100 K fspot = 0.5 (!) Example Key insight: - In the visible, starspot flux is 5-20x weaker than the ambient photosphere. - In the near-IR, starspot flux is only 2.5-4x weaker than ambient IGRINS ESPaDOnS

Slide 46

Slide 46 text

How to figure out which lines are attributable to starspots or ambient photosphere? portion of LkCa 4 IGRINS spectrum from November 2015

Slide 47

Slide 47 text

We forward  model the IGRINS spectra. Starfish is an open source spectral inference framework for stellar spectra. Czekala et al. 2015 github.com/iancze/Starfish Starfish parameters: 1. Tamb 2. logg 3. [Fe/H] 4. v sini 5. vz 6. Ω 7-9. c0, c1, c2... 10. GP scale 11. GP amplitude 12. σ scale 13. Tspot 14. fspot Intrinsic Starspots = + Composite Ambient Starspot Tspot = 2800 K Tamb = 4100 K **Lots  of  assump,ons  embedded  here

Slide 48

Slide 48 text

The spectrum has features from both ambient photosphere and starspots. λ (Angstrom) The constraint on filling factor comes from the range of flux ratios.

Slide 49

Slide 49 text

No content

Slide 50

Slide 50 text

Expected range of starspot properties

Slide 51

Slide 51 text

No content

Slide 52

Slide 52 text

LkCa 4 looks younger and less massive in *unspotted* tracks.

Slide 53

Slide 53 text

Unspotted tracks from G. Somers

Slide 54

Slide 54 text

Spotted tracks from G. Somers

Slide 55

Slide 55 text

Take aways 1. Stars are probably more spotted than we previously thought 2. Polar spots would have evaded most conventional methods of detecting and characterizing starspots, since they induce zero photometric modulation. 3. If LkCa 4 is representative of other young stars, the masses and ages of all young stars are considerably biased. 4. The stellar age biases change timescale available for planet formation. 5. What matters is the starspot coverage history, which is generally unobservable. 6. Teff measurement is hindered for highly inclined young spotted stars.

Slide 56

Slide 56 text

1. Stars are probably more spotted than we previously thought 2. Polar spots would have evaded most conventional methods of detecting and characterizing starspots, since they induce zero photometric modulation. 3. If LkCa 4 is representative of other young stars, the masses and ages of all young stars are considerably biased. 4. The stellar age biases change timescale available for planet formation. 5. What matters is the starspot coverage history, which is generally unobservable. 6. Teff measurement is hindered for highly inclined young spotted stars. Take aways

Slide 57

Slide 57 text

1. Stars are probably more spotted than we previously thought 2. Polar spots would have evaded most conventional methods of detecting and characterizing starspots, since they induce zero photometric modulation. 3. If LkCa 4 is representative of other young stars, the masses and ages of all young stars are considerably biased. 4. The stellar age biases change timescale available for planet formation. 5. What matters is the starspot coverage history, which is generally unobservable. 6. Teff measurement is hindered for highly inclined young spotted stars. Take aways

Slide 58

Slide 58 text

LkCa 4 has photometric monitoring going back 31 years. May 6, 1985 June 22, 2017

Slide 59

Slide 59 text

1. Stars are probably more spotted than we previously thought 2. Polar spots would have evaded most conventional methods of detecting and characterizing starspots, since they induce zero photometric modulation. 3. If LkCa 4 is representative of other young stars, the masses and ages of all young stars are considerably biased. 4. The stellar age biases change timescale available for planet formation. 5. What matters is the starspot coverage history, which is generally unobservable. 6. Need spectroscopy to assess full spot coverage fraction. Take aways

Slide 60

Slide 60 text

1. Stars are probably more spotted than we previously thought 2. Polar spots would have evaded most conventional methods of detecting and characterizing starspots, since they induce zero photometric modulation. 3. If LkCa 4 is representative of other young stars, the masses and ages of all young stars are considerably biased. 4. The stellar age biases change timescale available for planet formation. 5. What matters is the starspot coverage history, which is generally unobservable. 6. Need spectroscopy to assess full spot coverage fraction. Take aways

Slide 61

Slide 61 text

Recent evidence for large spot coverage in Pleiades 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 2600 3000 3500 4000 4500 5000 5500 6000 6500 TiO2n Teff (K) Inactive dwarfs PHOENIX(4.5) PHOENIX(5.0) Cubic splines fits Estimate of Tamb, Tspot, fspot in 304 LAMOST spectra 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 3000 3500 3800 4000 4500 5000 5500 6000 6500 fs Teff (K) Pleiades? Pleiades

Slide 62

Slide 62 text

Thanks!