IPhO 2009 · Theory problem 3

Why Are Stars So Large?

The stars are spheres of hot gas. Most of them shine because they are fusing hydrogen into helium in their central parts. In this problem we use concepts of both classical and quantum mechanics, as well as of electrostatics and thermodynamics, to understand why stars have to be big enough to achieve this fusion process and also derive what would be the mass and radius of the smallest star that can fuse hydrogen.

Figure 1

Figure 1. Our Sun, as most stars, shines as a result of thermonuclear fusion of hydrogen into helium in its central parts.

USEFUL CONSTANTS
Gravitational constant = \(G=6.7\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}}\) (misprint in the official paper: \(\mathrm{s^2}\) instead of \(\mathrm{s^{-2}}\))
Boltzmann’s constant = \(k=1.4\times10^{-23}\ \mathrm{J\,K^{-1}}\)
Planck’s constant = \(h=6.6\times10^{-34}\ \mathrm{m^2\,kg\,s^{-1}}\)
Mass of the proton = \(m_p=1.7\times10^{-27}\ \mathrm{kg}\)
Mass of the electron = \(m_e=9.1\times10^{-31}\ \mathrm{kg}\)
Unit of electric charge = \(q=1.6\times10^{-19}\ \mathrm{C}\)
Electric constant (vacuum permittivity) = \(\varepsilon_0=8.9\times10^{-12}\ \mathrm{C^2\,N^{-1}\,m^{-2}}\)
Radius of the Sun = \(R_S=7.0\times10^8\ \mathrm{m}\)
Mass of the Sun = \(M_S=2.0\times10^{30}\ \mathrm{kg}\)

1. A classical estimate of the temperature at the center of the stars.

Assume that the gas that forms the star is pure ionized hydrogen (electrons and protons in equal amounts), and that it behaves like an ideal gas. From the classical point of view, to fuse two protons, they need to get as close as \(10^{-15}\ \mathrm{m}\) for the short range strong nuclear force, which is attractive, to become dominant. However, to bring them together they have to overcome first the repulsive action of Coulomb’s force. Assume classically that the two protons (taken to be point sources) are moving in an antiparallel way, each with velocity \(v_{rms}\), the root-mean-square (rms) velocity of the protons, in a one-dimensional frontal collision.

1a What has to be the temperature of the gas, \(T_c\), so that the distance of closest approach of the protons, \(d_c\), equals \(10^{-15}\ \mathrm{m}\)? Give this and all numerical values in this problem up to two significant figures. 1.5

2. Finding that the previous temperature estimate is wrong.

To check if the previous temperature estimate is reasonable, one needs an independent way of estimating the central temperature of a star. The structure of the stars is very complicated, but we can gain significant understanding making some assumptions. Stars are in equilibrium, that is, they do not expand or contract because the inward force of gravity is balanced by the outward force of pressure (see Figure 2). For a slab of gas the equation of hydrostatic equilibrium at a given distance \(r\) from the center of the star, is given by

\[ \frac{\Delta P}{\Delta r}=-\frac{GM_r\rho_r}{r^2}, \]

where \(P\) is the pressure of the gas, \(G\) the gravitational constant, \(M_r\) the mass of the star within a sphere of radius \(r\), and \(\rho_r\) is the density of the gas in the slab.

Figure 2

Figure 2. The stars are in hydrostatic equilibrium, with the pressure difference balancing gravity.

An order of magnitude estimate of the central temperature of the star can be obtained with values of the parameters at the center and at the surface of the star, making the following approximations:

\[ \Delta P\approx P_o-P_c, \]

where \(P_c\) and \(P_o\) are the pressures at the center and surface of the star, respectively. Since \(P_c\gg P_o\), we can assume that

\[ \Delta P\approx-P_c. \]

Within the same approximation, we can write

\[ \Delta r\approx R, \]

where \(R\) is the total radius of the star, and

\[ M_r\approx M_R=M, \]

with \(M\) the total mass of the star.

The density may be approximated by its value at the center,

\[ \rho_r\approx\rho_c. \]

You can assume that the pressure is that of an ideal gas.

2a Find an equation for the temperature at the center of the star, \(T_c\), in terms of the radius and mass of the star and of physical constants only. 0.5

We can use now the following prediction of this model as a criterion for its validity:

2b Using the equation found in (2a) write down the ratio \(M/R\) expected for a star in terms of physical constants and \(T_c\) only. 0.5
2c Use the value of \(T_c\) derived in section (1a) and find the numerical value of the ratio \(M/R\) expected for a star. 0.5
2d Now, calculate the ratio \(M(Sun)/R(Sun)\), and verify that this value is much smaller than the one found in (2c). 0.5

3. A quantum mechanical estimate of the temperature at the center of the stars

The large discrepancy found in (2d) suggests that the classical estimate for \(T_c\) obtained in (1a) is not correct. The solution to this discrepancy is found when we consider quantum mechanical effects, that tell us that the protons behave as waves and that a single proton is smeared on a size of the order of \(\lambda_p\), the de Broglie wavelength. This implies that if \(d_c\), the distance of closest approach of the protons is of the order of \(\lambda_p\), the protons in a quantum mechanical sense overlap and can fuse.

3a Assuming that \(d_c=\dfrac{\lambda_p}{2^{1/2}}\) is the condition that allows fusion, for a proton with velocity \(v_{rms}\), find an equation for \(T_c\) in terms of physical constants only. 1.0
3b Evaluate numerically the value of \(T_c\) obtained in (3a). 0.5
3c Use the value of \(T_c\) derived in (3b) to find the numerical value of the ratio \(M/R\) expected for a star, using the formula derived in (2b). Verify that this value is quite similar to the ratio \(M(Sun)/R(Sun)\) observed. 0.5

Indeed, stars in the so-called main sequence (fusing hydrogen) approximately do follow this ratio for a large range of masses.

4. The mass/radius ratio of the stars.

The previous agreement suggests that the quantum mechanical approach for estimating the temperature at the center of the Sun is correct.

4a Use the previous results to demonstrate that for any star fusing hydrogen, the ratio of mass \(M\) to radius \(R\) is the same and depends only on physical constants. Find the equation for the ratio \(M/R\) for stars fusing hydrogen. 0.5

5. The mass and radius of the smallest star.

The result found in (4a) suggests that there could be stars of any mass as long as such a relationship is fulfilled; however, this is not true.

The gas inside normal stars fusing hydrogen is known to behave approximately as an ideal gas. This means that \(d_e\), the typical separation between electrons is on the average larger that \(\lambda_e\), their typical de Broglie wavelength. If closer, the electrons would be in a so-called degenerate state and the stars would behave differently. Note the distinction in the ways we treat protons and electrons inside the star. For protons, their de Broglie waves should overlap closely as they collide in order to fuse, whereas for electrons their de Broglie waves should not overlap in order to remain as an ideal gas.

The density in the stars increases with decreasing radius. Nevertheless, for this order-of-magnitude estimate assume they are of uniform density. You may further use that \(m_p\gg m_e\).

5a Find an equation for \(n_e\), the average electron number density inside the star. 0.5
5b Find an equation for \(d_e\), the typical separation between electrons inside the star. 0.5
5c Use the \(d_e\geq\dfrac{\lambda_e}{2^{1/2}}\) condition to write down an equation for the radius of the smallest normal star possible. Take the temperature at the center of the star as typical for all the stellar interior. 1.5
5d Find the numerical value of the radius of the smallest normal star possible, both in meters and in units of solar radius. 0.5
5e Find the numerical value of the mass of the smallest normal star possible, both in kg and in units of solar masses. 0.5

6. Fusing helium nuclei in older stars.

As stars get older they will have fused most of the hydrogen in their cores into helium (He), so they are forced to start fusing helium into heavier elements in order to continue shining. A helium nucleus has two protons and two neutrons, so it has twice the charge and approximately four times the mass of a proton. We saw before that \(d_c=\dfrac{\lambda_p}{2^{1/2}}\) is the condition for the protons to fuse.

6a Set the equivalent condition for helium nuclei and find \(v_{rms}(He)\), the rms velocity of the helium nuclei and \(T(He)\), the temperature needed for helium fusion. 0.5

Official English text of the problem.

Answers

All results use the approximations prescribed in the problem. Numerical results are rounded to two significant figures, with unrounded values retained in intermediate calculations.

  • 1a. \(\displaystyle T_c=\frac{q^2}{12\pi\varepsilon_0kd_c}=5.4\times10^9\ \mathrm{K}\).
  • 2a. \(\displaystyle T_c\simeq\frac{GMm_p}{2kR}\).
  • 2b. \(\displaystyle \frac{M}{R}\simeq\frac{2kT_c}{Gm_p}\).
  • 2c. \(\displaystyle \frac{M}{R}=1.3\times10^{24}\ \mathrm{kg\,m^{-1}}\).
  • 2d. \(\displaystyle \frac{M_S}{R_S}=2.9\times10^{21}\ \mathrm{kg\,m^{-1}}\). The classical prediction is \(4.7\times10^2\) times larger.
  • 3a. \(\displaystyle T_c=\frac{m_pq^4}{24\pi^2\varepsilon_0^2kh^2}\).
  • 3b. \(\displaystyle T_c=9.7\times10^6\ \mathrm{K}\).
  • 3c. \(\displaystyle \frac{M}{R}=2.4\times10^{21}\ \mathrm{kg\,m^{-1}}\simeq0.84\,\frac{M_S}{R_S}\).
  • 4a. \(\displaystyle \frac{M}{R}=\frac{q^4}{12\pi^2\varepsilon_0^2Gh^2}\).
  • 5a. \(\displaystyle n_e=\frac{3M}{4\pi m_pR^3}\).
  • 5b. \(\displaystyle d_e\simeq n_e^{-1/3}=\left(\frac{4\pi m_pR^3}{3M}\right)^{1/3}\).
  • 5c. \(\displaystyle R_{\min}=\left(\frac{\varepsilon_0h^4}{2Gq^2m_p^{5/2}m_e^{3/2}}\right)^{1/2}\).
  • 5d. \(\displaystyle R_{\min}=6.9\times10^7\ \mathrm{m}=0.099\,R_S\).
  • 5e. \(\displaystyle M_{\min}=1.7\times10^{29}\ \mathrm{kg}=0.083\,M_S\).
  • 6a. The condition is \(\displaystyle d_c(He)=\frac{\lambda_{He}}{\sqrt2}=\frac{h}{4\sqrt2\,m_pv_{rms}(He)}\), giving \[ v_{rms}(He)=\frac{\sqrt2\,q^2}{\pi\varepsilon_0h} =2.0\times10^6\ \mathrm{m\,s^{-1}}, \qquad T(He)=\frac{8m_pq^4}{3\pi^2\varepsilon_0^2kh^2} =6.2\times10^8\ \mathrm{K}. \]

Solution

1a. Classical temperature required for fusion

For a nonrelativistic ideal gas, equipartition gives the proton rms speed through

\[ \frac12m_pv_{rms}^2=\frac32kT_c. \]

The collision is head-on, but \(v_{rms}\) is the rms speed in the three-dimensional gas.

The two protons have equal and opposite momenta, so their center of mass is stationary. At closest approach, both momentarily stop. Taking their initial electrostatic potential energy to be negligible, conservation of energy gives

\[ 2\left(\frac12m_pv_{rms}^2\right) =\frac{q^2}{4\pi\varepsilon_0d_c}. \]

Thus

\[ 3kT_c=\frac{q^2}{4\pi\varepsilon_0d_c}, \qquad \boxed{T_c=\frac{q^2}{12\pi\varepsilon_0kd_c}}. \]

For \(d_c=10^{-15}\ \mathrm{m}\),

\[ T_c= \frac{(1.6\times10^{-19})^2} {12\pi(8.9\times10^{-12})(1.4\times10^{-23})(10^{-15})}, \]

so

\[ \boxed{T_c=5.4\times10^9\ \mathrm{K}}. \]

2a. Central temperature from hydrostatic equilibrium

Using the prescribed approximations, including a characteristic radial scale \(r\simeq R\),

\[ -\frac{P_c}{R}\simeq-\frac{GM\rho_c}{R^2}, \]

hence

\[ P_c\simeq\frac{GM\rho_c}{R}. \]

In fully ionized hydrogen, the proton and electron number densities are equal. Neglecting electron mass in the mass density,

\[ n_{p,c}=n_{e,c}\simeq\frac{\rho_c}{m_p}. \]

Both species contribute to the ideal-gas pressure:

\[ P_c=(n_{p,c}+n_{e,c})kT_c =\frac{2\rho_ckT_c}{m_p}. \]

Equating the pressure expressions and cancelling \(\rho_c\),

\[ \boxed{T_c\simeq\frac{GMm_p}{2kR}}. \]

2b. Mass-to-radius ratio

Rearranging the result of (2a),

\[ \boxed{\frac{M}{R}\simeq\frac{2kT_c}{Gm_p}}. \]

2c. Ratio predicted using the classical fusion temperature

Substitution of the temperature from (1a) gives

\[ \frac{M}{R} =\frac{2k}{Gm_p} \frac{q^2}{12\pi\varepsilon_0kd_c} =\frac{q^2}{6\pi\varepsilon_0Gm_pd_c}. \]

Numerically,

\[ \boxed{\frac{M}{R}=1.3\times10^{24}\ \mathrm{kg\,m^{-1}}}. \]

2d. Comparison with the Sun

The supplied solar mass and radius give

\[ \frac{M_S}{R_S} =\frac{2.0\times10^{30}}{7.0\times10^8} =\boxed{2.9\times10^{21}\ \mathrm{kg\,m^{-1}}}. \]

Using unrounded values, the comparison is

\[ \frac{(M/R)_{\mathrm{classical}}}{M_S/R_S} =\boxed{4.7\times10^2}. \]

The Sun’s ratio is therefore much smaller than the classical prediction. Within the hydrostatic model, the classical fusion temperature is far too high.

3a. Quantum estimate of the fusion temperature

For a proton moving at \(v_{rms}\), the de Broglie wavelength is

\[ \lambda_p=\frac{h}{m_pv_{rms}}. \]

The collision-energy relation from (1a) gives

\[ d_c=\frac{q^2}{4\pi\varepsilon_0m_pv_{rms}^2}. \]

Imposing the problem’s fusion condition,

\[ d_c=\frac{\lambda_p}{\sqrt2}, \]

we obtain

\[ \frac{q^2}{4\pi\varepsilon_0m_pv_{rms}^2} =\frac{h}{\sqrt2\,m_pv_{rms}}, \]

and consequently

\[ v_{rms}=\frac{\sqrt2\,q^2}{4\pi\varepsilon_0h}. \]

Using \(m_pv_{rms}^2=3kT_c\),

\[ T_c =\frac{m_p}{3k} \left(\frac{\sqrt2\,q^2}{4\pi\varepsilon_0h}\right)^2. \]

Therefore,

\[ \boxed{T_c=\frac{m_pq^4}{24\pi^2\varepsilon_0^2kh^2}}. \]

3b. Numerical quantum temperature

Substituting the given constants,

\[ T_c= \frac{(1.7\times10^{-27})(1.6\times10^{-19})^4} {24\pi^2(8.9\times10^{-12})^2 (1.4\times10^{-23})(6.6\times10^{-34})^2}, \]

which gives

\[ \boxed{T_c=9.7\times10^6\ \mathrm{K}}. \]

The corresponding proton rms speed is \(4.9\times10^5\ \mathrm{m\,s^{-1}}\).

3c. Mass-to-radius ratio using the quantum temperature

Using (2b) with the temperature from (3b),

\[ \frac{M}{R}=\frac{2kT_c}{Gm_p} =\boxed{2.4\times10^{21}\ \mathrm{kg\,m^{-1}}}. \]

The ratio to the solar value is

\[ \frac{M/R}{M_S/R_S}=\boxed{0.84}. \]

Thus the prediction is about \(16\%\) below the solar value, a close agreement for the order-of-magnitude approximations used.

4a. Common mass-to-radius ratio

Within this model, hydrogen fusion has the same characteristic temperature in every star. Combining (2b) and (3a),

\[ \frac{M}{R} =\frac{2k}{Gm_p} \frac{m_pq^4}{24\pi^2\varepsilon_0^2kh^2}. \]

The factors \(m_p\) and \(k\) cancel:

\[ \boxed{\frac{M}{R} =\frac{q^4}{12\pi^2\varepsilon_0^2Gh^2}}. \]

This depends only on physical constants. For the following parts, denote this common ratio by

\[ A\equiv\frac{q^4}{12\pi^2\varepsilon_0^2Gh^2}, \qquad M=AR. \]

5a. Mean electron number density

Each hydrogen proton supplies one electron. Neglecting electron mass,

\[ N_e=N_p\simeq\frac{M}{m_p}. \]

For a uniform spherical star of volume \(V=4\pi R^3/3\),

\[ \boxed{n_e=\frac{N_e}{V} =\frac{3M}{4\pi m_pR^3}}. \]

5b. Typical electron separation

The typical volume per electron is \(1/n_e\). Taking its cube root as the characteristic spacing,

\[ \boxed{d_e\simeq n_e^{-1/3} =\left(\frac{4\pi m_pR^3}{3M}\right)^{1/3}}. \]

5c. Radius of the smallest normal star

Equipartition for electrons gives

\[ m_ev_{e,rms}^2=3kT_c. \]

Their characteristic de Broglie wavelength is therefore

\[ \lambda_e=\frac{h}{m_ev_{e,rms}} =\frac{h}{\sqrt{3m_ekT_c}}. \]

The prescribed condition for avoiding electron degeneracy becomes

\[ d_e\geq\frac{\lambda_e}{\sqrt2} =\frac{h}{\sqrt{6m_ekT_c}}. \]

Using \(M=AR\) in the spacing formula,

\[ d_e^3=\frac{4\pi m_pR^2}{3A}. \]

Thus, along the model’s mass-radius relation, the spacing decreases as \(R^{2/3}\). The minimum permitted radius occurs at equality:

\[ \frac{4\pi m_pR_{\min}^2}{3A} =\left(\frac{h}{\sqrt{6m_ekT_c}}\right)^3. \]

Hence

\[ R_{\min}^2 =\frac{3Ah^3}{4\pi m_p(6m_ekT_c)^{3/2}}. \]

First substitute \(A=2kT_c/(Gm_p)\):

\[ R_{\min}^2 =\frac{h^3} {4\pi Gm_p^2m_e^{3/2}\sqrt{6kT_c}}. \]

From (3a),

\[ \sqrt{6kT_c} =\frac{q^2\sqrt{m_p}}{2\pi\varepsilon_0h}. \]

Therefore,

\[ R_{\min}^2 =\frac{\varepsilon_0h^4} {2Gq^2m_p^{5/2}m_e^{3/2}}, \]

or

\[ \boxed{ R_{\min} =\left( \frac{\varepsilon_0h^4} {2Gq^2m_p^{5/2}m_e^{3/2}} \right)^{1/2} }. \]

For a smaller radius along \(M=AR\), the electron waves would overlap according to the stated criterion, and the assumed ideal-gas description would cease to apply.

5d. Numerical minimum radius

Substitution gives

\[ \boxed{R_{\min}=6.9\times10^7\ \mathrm{m}}. \]

Dividing the unrounded result by the supplied solar radius,

\[ \boxed{\frac{R_{\min}}{R_S}=0.099}. \]

Thus the smallest normal star in this model has a radius of approximately one-tenth of the solar radius.

5e. Numerical minimum mass

Because \(M=AR\) with fixed positive \(A\), the minimum radius also gives the minimum mass:

\[ M_{\min} =AR_{\min} =\frac{q^4R_{\min}}{12\pi^2\varepsilon_0^2Gh^2}. \]

Numerically,

\[ \boxed{M_{\min}=1.7\times10^{29}\ \mathrm{kg}}, \qquad \boxed{\frac{M_{\min}}{M_S}=0.083}. \]

6a. Helium fusion

A helium nucleus has mass \(m_{He}\simeq4m_p\) and charge \(Q_{He}=2q\). Its de Broglie wavelength is

\[ \lambda_{He}=\frac{h}{4m_pv_{rms}(He)}. \]

The corresponding fusion condition is

\[ \boxed{ d_c(He)=\frac{\lambda_{He}}{\sqrt2} =\frac{h}{4\sqrt2\,m_pv_{rms}(He)} }. \]

For two helium nuclei approaching with equal and opposite rms velocities, energy conservation gives

\[ 4m_pv_{rms}(He)^2 =\frac{(2q)^2}{4\pi\varepsilon_0d_c(He)}. \]

Substituting the fusion condition and simplifying,

\[ \boxed{ v_{rms}(He)=\frac{\sqrt2\,q^2}{\pi\varepsilon_0h} }. \]

This is four times the proton speed obtained in (3a), so

\[ \boxed{v_{rms}(He)=2.0\times10^6\ \mathrm{m\,s^{-1}}}. \]

Equipartition now gives

\[ T(He)=\frac{4m_pv_{rms}(He)^2}{3k}. \]

Since the mass is four times greater and the speed is four times greater,

\[ T(He)=4(4^2)T_c=64T_c, \]

where \(T_c\) is the hydrogen fusion temperature from (3b). Therefore,

\[ \boxed{ T(He)=\frac{8m_pq^4}{3\pi^2\varepsilon_0^2kh^2} =6.2\times10^8\ \mathrm{K} }. \]

The solution was checked against the official solution.