ASTR 498N                   Lecture 16             Stars That Can Be

 

(online at www.astro.umd.edu/~drabin/)

               

 Be all that you can be.

               

 

Minimum mass of a main sequence star

 

Back in Lecture 5, we made rough order-of-magnitude estimates of the central pressure and temperature of a star, based on hydrostatic support and the ideal gas law, :

 

                                                               

 

                     

 

In Homework 4, you derived an upper bound on Pc related to the central density rather than the average density:

 

                     

 

Let us set the central pressure to this upper limit.  For actual stars, which are centrally concentrated, this approximation is better than expressions pegged to the average density; realistic models yield central pressures that are roughly half the maximum value. 

 

Consider the initial contraction of a protostellar cloud.  Equating our approximation to the central pressure to the perfect gas pressure, we find

 

                     

 

which shows how the central temperature rises with the central density.

 

What happens as the contraction continues?  (We largely follow the presentation in The Physics of Stars by A. C. Phillips).  Set aside for the moment the possibility of nuclear burning.  When the central density becomes high enough, the pressure of degenerate, nonrelativistic electrons becomes significant.  The total pressure becomes

 

                              where  

 

Assume for simplicity a composition of pure hydrogen, so that

 

                     

 

Again equating the two expressions for central pressure,

 

                     

 

Clearly this expression has a maximum, shown by differentiation to have the value  at a density of .  Substituting for A and B, we find

 

                     

 

Suppose that proton-proton fusion “ignites” at a characteristic temperature Tign.  Then the mass of a star that barely reaches this central temperature is, from the previous equation,

 

                     

 

The ignition temperature depends on the environment of the contracting material.  When the nuclear power produced in a particular region slightly exceeds the power that is lost, the region heats up and burning takes hold.  If we suppose that the ignition temperature of hydrogen is about 1.5 × 106 K (one tenth the central temperature of the Sun), the expression above yields

 

Divertissement

 

Taking the ignition temperature to be 0.1Tc is unsatisfyingly arbitrary.  Could we instead relate the ignition temperature to some intrinsic characteristic of the nuclear reaction? 

 

You showed as a homework problem that the fusion rate may be expressed approximately as

           

 

If we ignore the temperature dependence of S(E0), this has the form

 

                     

 

Differentiating this expression yields a maximum at x = 0.296 or kT = 0.84 EG

 

The first step of the proton-proton chain, , has Gamow energy EG = 0.493 MeV.  Thus, the p-p chain generates maximum energy at a temperature 0.84(0.493) MeV or 4.9 × 109 K (!), whereas the actual burning temperature at the center of the Sun is 300 times cooler.  Gravitational contraction stops when the central temperature is high enough for nuclear energy generation to supply the needed luminosity, not when the nuclear efficiency reaches some particular fraction of its maximum.

 

 

Maximum mass of a main sequence star

 

In Lecture 5 we derived two forms of the virial theorem for a self-gravitating system with internal pressure:

 

                    2K + U = 0   (nonrelativistic)

                    K + U = 0     (ultrarelativistic)

 

where K is the kinetic energy of the system,   U is the potential energy, and the total energy is always E = K + U.  [Recall that the difference between the two forms arises because pressure is equal to 2/3 of energy density for an ideal gas but only 1/3 for a completely relativistic gas such as radiation.]  The total energy is negative for a gravitationally bound system. 

 

Thus, if the internal pressure of a gravitating system is dominated by radiation pressure, the system will be on the margin between bound and unbound.  Note also that contraction, which causes a nonrelativistic system to become more bound, does not help in this case:  the total energy stays near zero.  We expect such a system to be highly prone to disruption, particularly in a region of dynamic star formation where there are a variety of potentially significant external perturbations.

 

As in Part 2/Problem 2 of the mid-term exam, characterize the relative importance of ideal gas pressure and radiation pressure by the parameter β:

 

                     

 

where  and  .  From these equations it is easy to derive

 

                     

 

and, by equating these two expressions for T,

 

                     

 

This is a general expression, valid for any combination of radiation and ideal gas pressure (except β = 0 or 1).  If we equate this expression at the center of the star to the same estimate of the central pressure we used to estimate the minimum mass above, we get

 

                       

or

 

                        (μ = 0.5)

The graph of this function decreases monotonically with β, or, equivalently, increases monotonically with 1 β (fraction of pressure due to radiation).  In the absence of an explicit stability analysis, we can make a rough guess that a star will be prone to instability if more than half the pressure is due to radiation, 1 β > 0.5.  The equation gives

 

                     

 

From observation, stars more massive than 50 solar masses are quite rare.

 

 

A fundamental unit for stellar masses

 

We have estimated that the mass range of main sequence stars spans something less than four orders of magnitude (0.01100 M), with one solar mass squarely in the middle of the range.  This prompts us to ask whether there is a dimensionless number that sets the mass scale. 

 

Recall that the fine structure constant characterizes the strength of the electromagnetic force by comparing the electrostatic potential energy between fundamental units of charge, separated by a fundamental distance, with a fundamental energy, the rest-mass energy:

 

                                          electrostatic potential energy

                                               reduced Compton wavelength

                                fine structure constant

 

Let’s apply this logic to the gravitational interaction between nucleons:

 

                               gravitational potential energy

                                             reduced Compton wavelength

                                  dimensionless constant = 5.9 × 1039

 

Another way to write αG  is  where  is a fundamental length known as the Planck length. 

 

Our expression for the minimum stellar mass can be written in terms of αG  as

 

                     

 

if we again use the estimate Tign = 1.5 × 106 K.  Similar, our expression for the maximum mass becomes

                   

                     

 

for β = 0.5 and μ = 0.5.

 

The commonality of the expressions for Mmin and Mmax  prompts us to introduce

 

                     

 

as a fundamental unit for stellar masses.   Note that the number of nucleons in a fundamental stellar mass is determined solely by αG: