Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.mso.anu.edu.au/~wood/lectures/2007/assignment_2.pdf Дата изменения: Thu Mar 15 02:33:22 2007 Дата индексирования: Tue Oct 2 08:15:00 2012 Кодировка: Поисковые слова: п п п п п п п п п п п п п п п

Assignment 2

Semester 1, 2007

1. The Lane-Emden equation is

1d d 2 2 d d
1

= -

n

for a polytrope of index n in which P = K1+ n and = c n , c being the central density. Show that the total mass of the polytrope is given by (n + 1)Pc M = 4 1 3 (4 ) 3 Gc where s is the value of at the stellar surface. 2. As a crude approximation, we can replace an equation such as dP GMr =- 2 dr r Pc GM M 2 3, R RR where Pc is the central pressure and R and M are the stellar radius and mass, respectively. Use approximations of this type to show that for a sequence of stars in hydrostatic equilibrium, consisting of a perfect gas, with a similar average opacity throughout, and in which energy is transported by radiation Tc M , R and by a dimensional relation

3 2

-

2

d d

,
=
s

L M 3. Here, Tc is the central temperature. (Note. The equation L M 3 is not too bad an approximation for main-sequence stars with M > M : they actually follow a relation L M 3.8 .) The above approximations do not use any information about the source of the luminosity (nuclear fusion)! Given that most of the stellar luminosity in a main-sequence star is produced by nuclear fusion near the centre of the star, and that the rate of nuclear energy production increases very rapidly with central temperature Tc , use the results above to describe how a newly-formed star comes to a state of thermal and nuclear-burning equilibrium after it condenses out from a molecular cloud. 3. This example explores the evolution of the sun. The attached figure shows the evolution of the sun in the theoretical Hertzsprung-Russell diagram, as well as a plot of the hydrogen abundance XH against mass fraction Mr /M in the sun. (a) Use the plot to make a rough estimate the mass of hydrogen that has been converted to helium in the sun during its lifetime, from birth to the present day. Given that the mass of proton is 1.00870 atomic mass units and the mass of a helium nucleus is 4.0026 atomic mass


units, estimate the sun's average rate of energy production by nuclear reactions (in solar luminosities) over its lifetime. Take the age of the sun as 4.5 в 109 yr. Does nuclear energy production seem to successfully account for the energy output of the sun? (b) Comment on the changes in the abundance profiles as the sun ages. In particular, for each profile, where is hydrogen burning occurring in each model, why is there a discontinuity in the abundance profile near mass-fraction 0.31 when the sun is near the base of the red giant branch, and why is the surface hydrogen abundance slightly reduced in the giant branch model. Some constants you may need are: M = 1.989 в 1033 g; L = 3.9 в 10 light c = 3 в 1010 cm s-1 ; atomic mass unit = 1.66 в 10-24 g.
33

erg s

-1

; speed of


The CNO bi-cycles The first way to process H into 4 He is by the PP chains, while a second way is through a series of nuclear reactions involving different isotopes of carbon, nitrogen and oxygen. The reactions are as follows: The CN cycle
12

C + p -
13

13 13 14 15 15 12

N+ C+e + N+ O+ N + e+ + C + 4 He
e + e

(1) (2) (3) (4) (5) (6)

N -

13 14

C + p - N + p -
15

O -

15

N + p -

or the ON cycle, which starts with
14

N(p, )15 O( + )15 N

(4, 5 )

then follows
15 16

N + p - O + p -
17

16 17 17 14

O+ F+ O + e+ + N + 4 He.
18 e

(7) (8) (9) (10)
19

F -

17

O + p -

Note that we are ignoring the extra cycles that involve

O and

F.

Part A: By summing the reacting and the product particles for the CN and ON cycle, show that the C, N and O nuclei only play the role of catalysts. Compute the Gamow peak and width for the 12 C(p, )13 N reaction using
22 E0 = 0.122 Za ZX AT 2 9 1/3

(MeV)

(11)

and
22 = 0.237 Za ZX AT 5 9

1/6

(MeV)

(12)

where Za and Z

X

are the element number, A the reduced mass A= Aa AX Aa + A X (13)

where Aa , AX are the atomic masses. Compute E0 and for i) 15 million K, appropriate for the centre of the Sun, and ii) for 25 million K, characteristic of CNO burning on the upper main sequence. Comment on your answers, in particular compared to the lowest experimental energy measurements available. Hint: Look at values in the recent NACRE compilation article, Angulo, C., 1999, Nuclear Physics A, 656, p3-187.


Part B: We can define lifetime of a nuclear species in a given environment a (X ) = (
aX

Na )

-1

,

(14)

where a (X ) is the lifetime of species X against reactions with species a, and aX =< v >, the reaction rate cross section, and Na the number density of species a, given by Na = (X (a))/Amu , where X (a) is the mass fraction of a, is the density, A the atomic mass and mu the atomic mass unit, given by 1.66 в10-24 grams. If we assume that the radioactive nuclides 13 N, 15 O, 17 F beta decay essentially instantaneously to 13 C, 15 N, 17 O, then we can write a simple set of differential equations that govern the abundances of the C, N and O nuclei as a function of time. These equations are: d12 C dt d13 C dt 14 N d dt 15 N d dt d16 O dt 17 O d dt
12

=-
12

C

15

= = =

12 13 C 13 14 N



12 C - -

+
13

N



(15) (16)
17

15

C + O

13 14 N 14 15 N







(17) (18) (19) (20)

17



-

14 15



= =

N

15 16

15 16 16 O 17 O - , 16 17

-

O

where, in the short-hand notation used here, 12 refers to the lifetime of 12 C against reactions with protons, and 12 C refers to the number density of 12 C atoms in the plasma. The constants and are branching ratios from the CN cycle into the ON cycle. Compute the equilibrium values (12 C/13 C)e , (12 C/14 N)e , and (16 O/14 N)e , assuming a temperature of 25 million K (MK) and XH = 50. Assume that the 14 N/15 N ratio has reached equilibrium. Comment on your answers. Hint: Use Table 5­3 from Clayton's textbook, where he provides log( XH ) 100 (21)

as a function of T6 , the temperature in units of 106 K, and where the lifetimes are expressed in years. The branching ratio has a weak temperature dependence, and multiplication factors at different temperatures and densities are also provided in Table 5­3. Table 5­3 is attached to the end of this assignment sheet.


Part C: The variation of the abundance ni of a species i at a given time t, dni (t)/dt, can be approximated by the exponential law, dni = ke dt
-t/

(22)

where is the lifetime, defined in Part B, and k = -1/ . This differential equation has the solution n(t) = C e-t/ , where C is some constant. Using the lifetimes computed previously in Part B for 25MK and XH = 50, how long does it take the initial abundances of 12 C to be reduced by a factor of 10? What about 16 O? Comment on the difference between the two times? What about at 75 million K and X H = 50? [Note: The CN cycle reaches equilibrium in a time of the order of 12 , which is much faster than any significant interchange of nuclei between the CN cycle and the ON cycle. Those two interchange times are characteristically 103 14 and 16 .]