In-Class Assignment 8

Contents

In-Class Assignment 8#

Exploring Opacities in Stellar Models#

import numpy as np
import matplotlib.pyplot as plt
import pandas as pd

Your favorite model files are again:

\(1 M_\odot\): M1_default_profile8.data; M1_default_profile218.data;
\(15 M_\odot\): M15_aprox21_profile8.data; M15_aprox21_profile19.data;

\(\kappa_{\rm{rad}}\)#

Using the 1\(M_\odot\) and 15\(M_\odot\) MESA models at their early evolutionary phase, plot the opacity. Consider using the mass scale in Figure 4.6 of HKT to emphasis the outer layers.

Hint: this requires you to normalize by the total mass of the star!

On the same plot, plot an estimate of the \(\kappa_{\rm{rad}}\) of your choosing assuming a uniform mixture. An example of this is on HKT 216. Do the opacity profiles match your estimated value? Which opacity source did you choose and why?

## a results here

\(\kappa_{\rm{cond}}\)#

Perform the same exercise here for the later evolutionary profiles. This time, compare with an estimate for the \(\kappa_{\rm{cond}}\) in HKT 4.72. Do the profiles match? Why?

## b results here