MESA/GYRE & gravity-mode oscillations of massive stars MESA: thank you Bill Paxton GYRE: thank you Rich Townsend
MESA/GYRE & gravity-mode oscillations of massive stars Lecturer: Conny Aerts TAs: Timothy Van Reeth & Evan Bauer
Setting the scene Stellar evolution mainly tested from surface properties while life directed by stellar interior This MESA/GYRE class focuses on young stars with a well-developed convective core AIMS: probe physics in near-core regions from gravity- mode oscillations
Looking inside stars: quakes
The beauty of gravity modes Quakes = solution of perturbed SSE in terms of periodic eigenfunctions Oscillation modes described in terms of spherical harmonics: Depending on dominance of restoring force: 1. pressure modes (acoustic waves) 2. gravity modes (horizontal waves) Gravity waves propagate in the radiative zone,from surface to the near-core region
The beauty of gravity modes PhD Thesis Valentina Schmid Leuven, 2016 Green: l=0 Red: l=1, n=3 Blue: l=1, n=-20
MESA/GYRE computations of quakes Compute stellar models with MESA & their pulsations with
Topic 1: MESA predictions of gravity modes
Theoretical Predictions No rotation, no magnetic field, convective core, radiative envelope
Theoretical Predictions No rotation, no magnetic field, convective core, radiative envelope MESA
Minilab 1: compute the asymptotic period spacing for dipole g-modes of a massive star from ZAMS to TAMS
Topic 2: GYRE predictions of gravity modes
A real B star: KIC 10526294 Pápics et al. (2014): Teff=13000 K log g=4.1 vsini < 10 km/s M=3.25M
2 real F stars: KIC 10080943 Schmid et al. (2015)
2 real F stars: KIC 10080943 Schmid & Aerts (2016)
Minilab 2: compute the real period spacing for dipole zonal g-modes of a massive star from ZAMS to TAMS for n=-5,,-75; compare with the asymptotic value: how good is the latter? (How) does its quality change during evolution? P = Pn+1,l -Pn,l
Mode trapping due to chemical composition gradient As pioneered by Miglio et al. (2008) prior to space data! M=6M, Xc=0.50 l=1 modes n=18,20,22 Xc=0.50 chemical composition gradient left behind by receding convective core
Topic 3: turning on rotation in GYRE (not in MESA )
Effect of Rotation: tilt Rotation lifts frequency degeneracy:
Rotation & g-modes: TA rigid rotation Solve Laplace tidal equation, as function of spin parameter Convention in GYRE: m>0: prograde m<0: retrograde
An example: KIC11721304 Teff=7160 K, log g=4.13, vsini=28 km/s (Van Reeth et al. 2015)
Prograde & Retrograde Modes KIC 8375138: Teff=7070 K, log g=4.05, vsini=140 km/s (Van Reeth et al. 2015)
Apply TA to derive interior rotation Sample of 40 F-type g-mode pulsators (Van Reeth et al. 2016)
Minilab 3: compute GYRE frequencies of prograde dipole g-modes of massive stars in TA from ZAMS to TAMS for rigid rotation at 5,10,15,20% critical; compare with non-rotating case
Topic 4: turning on mixing & core overshooting in MESA & rotation in GYRE
Effect of Chemical Gradient: Dips Van Reeth et al. (2015) Stars with shrinking core during# core-h burning: chemical gradient in near-core region
A real star:kic 7760680 From Moravveji et al. (2015) Convective core overshooting fov or Dov Radial diffusive mixing in radiative envelope log Dmix (in cm^2/s)
Effect of Chemical Gradient: Dips KIC 5708550: Teff=7000 K, log g=3.97, vsini=68 km/s (Van Reeth et al. 2015)
A real star:kic 7760680 Pápics et al. (2014)
Dips & chemical mixing Various types of mixing Mixing disconnected with rotation Rotationally induced mixing Schmid & Aerts (2016)
Dips & chemical mixing Various types of mixing in MESA Mixing disconnected with rotation Rotationally induced mixing log Dmix Schmid & Aerts (2016)
A real star:kic 7760680 KIC7760680: M=3.25M, Xc=0.50, frot=0.48/d, fov=0.024 Hp log Dmix=0.75±0.25 (Moravveji et al. 2016)
2 real F stars: KIC 10080943 Maxilab: primary needs more complex mixing: rotational or other origin?
Maxilab: compute GYRE period spacings of prograde dipole gmodes of massive stars in TA from ZAMS to TAMS for rigid rotation at 5,10,15,20% critical and for diffusive mixing with D =10, 50,100 cm^2/s. Then repeat for different core overshooting values fov=0.010, 0.015, 0.020. Plot the period spacing patterns. mix For the ultrafast or supermotivated: compare with data in Van Reeth et al. (2016)