# Main sequence

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

### Main Sequence

\We understand now how to calculate the luminosity produced by fusion, in addition to understanding how luminous a star needs to be from energy transport. We have

$L_{\rm {fusion}}=\int dM_{r}\epsilon (\rho ,T)\,\!$ and

$L\propto M^{3}{\rm {\;\;(Thomson)}}\,\!$ $L\propto M^{5.5}R^{-0.5}{\rm {\;\;(free-free)}}\,\!$ $L\propto M^{4/7}R^{2}{\rm {\;\;(Convection)}}\,\!$ When the fusion luminosity is equal to the luminosity that can leak out, then the star is in steady state. This condition, plus hydrostatic equilibrium, give the properties of stars on the main sequence. In other words, given the mass, you can solve for $L,T_{c},T_{\rm {eff}},R$ . \\We can start with an object that is not fusing, but is radiating energy away. That is, it is undergoing Kelvin Helmholtz contraction. As the object contracts, its central temperature increases, which in turn increases the luminosity due to fusion (once it is hot enough for at least some fusion to occur). This fusion luminosity will eventually increase to the point that it is equal to the luminosity being radiated away, at which point the object ceases contracting. Once this happens, the star ceases contracting, the central temperature stays fixed, the luminosity is fixed, and the star is on the main sequence. Well, almost fixed. The composition of the star will change as it fuses, which will cause very slow changes in these properties even while it is on the main sequence. This phase of hydrogen fusing lasts longer than any other stage of a star’s lifetime. \\For the Sun, we showed that radiative diffusion, Thomson scattering, and the virial theorem imply that $L\sim L_{\odot }$ , independent of fusion. Then, knowing that fusion of hydrogen in to helium is the energy source, we determined that the central temperature of the Sun is $T_{c}\sim 10^{7}$ K. We then used the virial theorem to determine the radius of the star to be $R\sim 10^{11}$ cm. The way to think about this is that the central temperature is not set by the virial theorem, but fusion sets the central temperature, and this and the virial theorem set the radius. Once we know the radius and luminosity, we can use the Stefan-Boltzmann equation to determine the effective temperature to be $T_{\rm {eff}}\sim 6000$ K. \\What we want to do now is to do this for other stars, and explain the main sequence portion of the Hertzsprung-Russell diagram. In principle, this is easy, as we can just follow the exact same procedure for other stars as we did for the Sun. What complicates this is that for other stars, the opacity, energy generation, energy transport, and dominant pressure may be different. This is because all of these properties depend on density and temperature, both of which depend on the mass of the star. Thus, in order to do this correctly, we need to know what the dominant mechanism for each of these things is for other stars. \\We showed on homework that at the center of the Sun, free-free opacity and Thomson opacity are roughly of equal importance. Thomson opacity is constant with density and temperature, but $\kappa _{f-f}\propto \rho T^{-7/2}$ . For stars more massive than the Sun, the central temperature is larger and the density is lower, meaning Thomson scattering dominates for stars more massive than the Sun, and free-free dominates for stars less massive than the Sun. \\The energy produced by the CNO cycle is about 1% of the energy produced by the proton-proton chain in the Sun. Since the CNO cycle scales as $T^{20}$ whereas the p-p chain scales as $T^{4}$ , the CNO cycle is dominant for stars somewhat more massive than the Sun, and the proton-proton chain is dominant for stars the Sun’s mass and lower. \\Gas pressure dominantes for all but the most massive stars. Stars less massive than the Sun have convective envelopes and radiative cores. Stars more massive than the Sun have convective cores and radiative envelopes. Knowing all of this, we can look at the main sequence of mass ranges of stars. \\First, we’ll look at stars less massive than the Sun down to about one third of a solar mass. These have free-free opacity, photons transport the energy, proton-proton chain, and are gas pressure dominated. We know that

$L_{\rm {rad}}\propto M^{5.5}R^{-0.5},\,\!$ $L_{p-p}=\int \epsilon _{p-p}dM_{r}=L_{\rm {rad}},\,\!$ $L_{p-p}\sim \epsilon _{c}M,\,\!$ $L_{p-p}\sim \rho T^{4.5}M,\,\!$ $T_{c}\propto {\frac {M}{R}},\,\!$ $L_{p-p}\sim {\frac {M^{6.5}}{R^{7.5}}}.\,\!$ Equating the radiative losses to the energy generation,

$L_{p-p}\sim L_{\rm {rad}},\,\!$ ${\frac {M^{6.5}}{R^{7.5}}}\propto {\frac {M^{5.5}}{R^{0.5}}},\,\!$ $R\propto M^{1/7},\,\!$ $R=R_{\odot }\left({\frac {M}{M_{\odot }}}\right)^{1/7}\,\!$ The last step is possible because we are able to normalize to the Sun. Likewise for central temperature, we find

$\,\!$ $T_{c}\propto {\frac {M}{R}}\propto M^{6/7},\,\!$ $T_{c}\propto (1.5\times 10^{7}K)\left({\frac {M}{M_{\odot }}}\right)^{6/7}.\,\!$ The radial change with mass is tiny, so we can just say

$L=L_{\odot }\left({\frac {M}{M_{\odot }}}\right)^{5.5}.\,\!$ We can also solve for $R$ in terms of $L$ , and find

$R\propto L^{1/40}.\,\!$ This is tiny, so we can just say that

$T_{\rm {eff}}=5800K\left({\frac {L}{L_{\odot }}}\right)^{1/4}.\,\!$ This is everything we would want to know for a star that is between roughly 0.3 and 1 solar mass, but outside this mass range, these results do not hold. Let’s do this again then for stars that are more massive than the Sun, ranging from roughly 1–50 solar masses. Here, CNO is dominant, gas pressure is dominant, and Thomson opacity is dominant. We showed in lecture that

$L_{\rm {rad}}=L_{\odot }\left({\frac {M}{M_{\odot }}}\right)^{3}.\,\!$ The energy generated by fusion is

$L_{\rm {fusion}}=\int \epsilon dM_{r}\sim \epsilon _{c}M,\,\!$ $\epsilon _{c}\propto \rho T^{20},\,\!$ $\epsilon _{c}\propto {\frac {M^{21}}{R^{23}}}.\,\!$ $L_{\rm {fusion}}\propto {\frac {M^{22}}{R^{23}}}.\,\!$ Equating the radiative luminosity to the fusion luminosity gives

${\frac {M^{22}}{R^{23}}}\propto M^{3},\,\!$ $R\propto M^{19/23}\propto M^{0.8},\,\!$ $R=R_{\odot }\left({\frac {M}{M_{\odot }}}\right)^{0.8}.\,\!$ Using the virial theorem gives the central temperature dependence of

$T_{c}=(1.5\times 10^{7}K)\left({\frac {M}{M_{\odot }}}\right)^{0.2}.\,\!$ For stars on the CNO cycle, the temperature only needs to change slightly to keep up with the extra luminosity. Finally, Stefan-Boltzmann gives

$T_{\rm {eff}}=5800K\left({\frac {M}{M_{\odot }}}\right)^{3/8}.\,\!$ Or, in terms of luminosity, useful for plotting the HR diagram,

$T_{\rm {eff}}=5800K\left({\frac {L}{L_{\odot }}}\right)^{1/8}.\,\!$ Thus there is a large range of luminosities for stars on the main sequence, despite a relatively small range in effective temperatures. \\We could do this again, this time for stars more massive than 50 solar masses. These only change in that the gas pressure no longer dominates, with radiation pressure taking over. This gives $L\propto M$ . We will not actually do this. Likewise, we could do fully convective stars, but we’ll leave that for homework. \\Instead, let’s explain why massive stars have convective cores. Convection sets in when there is a negative entropy gradient. A more useful condition in practice is that entropy sets in when

${\frac {d\ln T}{d\ln P}}>{\frac {\gamma -1}{\gamma }}.\,\!$ Remember that $\gamma$ , the adiabatic index, is for the pressure providing particles. Then, for gas pressure dominated stars, we have $\gamma =5/3$ . We showed that the left hand side can be written

${\frac {d\ln T}{d\ln P}}={\frac {1}{4}}{\frac {P}{P_{\rm {rad}}}}{\frac {L}{L_{edd}}}{\frac {L_{r}/L}{M_{r}/M}}.\,\!$ Thus convection sets in if

${\frac {P}{P_{\rm {rad}}}}{\frac {L}{L_{edd}}}{\frac {L_{r}/L}{M_{r}/M}}>{\frac {8}{5}}.\,\!$ Rearranging this in terms of the enclosed mass, convection sets in if

${\frac {M_{r}}{M}}<{\frac {5}{8}}{\frac {P}{P_{\rm {rad}}}}{\frac {L}{L_{edd}}}{\frac {L_{r}}{L}}.\,\!$ For the CNO cycle, $L_{r}\approx L$ even for very small radii. We can get a rough sense of when this is important by looking at the Sun (even though it does not strictly apply for the Sun). We have $P/P_{\rm {rad}}\sim 3000$ , and $L/L_{edd}\sim 4\times 10^{-5}$ . This gives convection for

${\frac {M_{r}}{M}}<0.1.\,\!$ This is why stars more massive than the Sun have convective cores. The Sun only barely avoids this, because it is dominated by the proton-proton chain, which is slightly less temperature sensitive. The convective core sets in for stars about 1.3 $M_{\odot }$ and above.