After Deuterium Ignition

One gram of interstellar matter generates energy by the reaction of $^2{\rm H}(p,\gamma)^3{\rm He}$ as

\begin{displaymath}
\delta\equiv \left(\frac{1}{m_p\mu_H}\right) \left(\frac{n_...
... D}}{n_{\rm H}}\right)
\epsilon_{\rm D\rightarrow He} \simeq
\end{displaymath} (4.119)

where we used an energy per reaction of $^2{\rm H}(p,\gamma)^3{\rm He}$ as $\epsilon_{\rm D\rightarrow He}= 5.5 {\rm MeV}$ and a composition of Deuterium of ${n_{\rm D}}/{n_{\rm H}}=2.5\times 10^{-5}$. This gives an estimate of the Deuterium-produced luminosity as
\begin{displaymath}
L_{\rm D}=\dot{M}\delta
=15L_\odot \left(\frac{\dot{M}}{10^{-5}M_\odot{\rm yr}^{-1}}\right),
\end{displaymath} (4.120)

for a star with the accretion rate of $\dot{M}\simeq 10^{-5}M_\odot{\rm yr}^{-1}$.

Entropy inversion $\partial s/\partial r<0$ is unstable for the convection (see $\S$2.7 in p.[*]). Convection is generated in the region of $\partial s/\partial M<0$. Convective motion transfers the energy and the entropy. The specific entropy peak generated by the Deuterium burning is smoothed out by the convection. The convection induced with the Deuterium burning brings a fresh gas containing Deuterium into the hot Deuterium burning region. Thus, although the consent of Deuterium is small and Deuterium burning could not continue long, Deuterium burning continues actually as long as the convection can mix the outer Deuterium rich fresh gas into the burning region due to the convection.

Figure: Luminosities in protostars accreting at $\dot{M}_*=10^{-5}M_\odot{\rm yr}^{-1}$ (left) and $\dot{M}_*=10^{-4}M_\odot{\rm yr}^{-1}$ (right) (Palla & Staler 1993). Accretion luminosity $L_{\rm acc}$, Deuterium burning luminosity $L_{\rm D}$, maximum radiatively transferable luminosity $L_{\rm rad}$ and luminosity at the surface $L_{\rm surf}$ are plotted.
\begin{figure}\begin{center}
\centering\leavevmode
\epsfxsize =.45\columnwidth...
...=.45\columnwidth \epsfbox{eps/Palla_Stahler_93_L-4.ps}
\end{center}
\end{figure}

In Figure 4.22, luminosities in protostars accreting at $\dot{M}_*=10^{-5}M_\odot{\rm yr}^{-1}$ (left) and $\dot{M}_*=10^{-4}M_\odot{\rm yr}^{-1}$ (right) are plotted (Palla & Stahler 1993). $L_{\rm rad}$ represents the maximum luminosity transferable radiatively for a star with mass $M_*$ and radius $R_*$ as

\begin{displaymath}
L_{\rm rad}\sim 1 L_\odot \left(\frac{M_*}{1M_\odot}\right)^{11/2}
\left(\frac{R_*}{R_\odot}\right)^{-1/2},
\end{displaymath} (4.121)

where the Krammers opacity $\kappa \propto \rho T^{-3.5}$ is applied for stars with $M\mbox{\raisebox{0.3ex}{$<$}\hspace{-1.1em}
\raisebox{-0.7ex}{$\sim$}} 5M_\odot$. Actual luminosity at the surface of the star
\begin{displaymath}
L_{\rm surf}=4\pi R_*^2 T_{\rm eff}^4,
\end{displaymath} (4.122)

is larger than $L_{\rm rad}$ for low mass stars ( $M_* < M_1= 2.4 M_\odot$ for $\dot{M}_*=10^{-5}M_\odot{\rm yr}^{-1}$ and $M_1= 3.4 M_\odot$ for $\dot{M}_*=10^{-4}M_\odot{\rm yr}^{-1}$). Excess energy is transferred by the convection, $L_{\rm surf}=L_{\rm rad}+L_{\rm con}$. This figure also shows that except for extremely low-mass object the accretion luminosity is the major energy source over Deuterium burning $L_{\rm acc} > L_{\rm D}$. Since massive objects with $M_* > M_3$ satisfies $L_{\rm rad}\simeq L_{\rm surf} > L_{\rm acc}$, interior of such stars are fully radiative (energy is transferred by the radiation). $M_3=3.9 M_\odot$ for $\dot{M}_*=10^{-5}M_\odot{\rm yr}^{-1}$ and $M_3=5.9 M_\odot$ for $\dot{M}_*=10^{-5}M_\odot{\rm yr}^{-1}$.

Figure 4.23 (left) illustrates how the Deuterium burning in protostars evolves. Accumulating infalling gas and increasing $M_*$, the structure of protostars is changed. (a) in low-mass stars, convection feeds the Deuterium of the accreting matter to the burning center. (b) as increasing the mass and thus the temperature, the opacity of the stellar interior decreases. Finally a radiative zone appears and it prevents from mixing Deuterium. (c) as a consequence, central Deuterium is depleted and energy generation rate is reduced. Thus the central region becomes radiative. (d) For higher mass stars, Deuterium ignites as a shell source just outside the radiative central region. Due to this Deuterium shell burning, which is similar to the Hydrogen shell burning in the red giants, radius of the protostar expands. Figure 4.23 (right) plots the mass-radius relation of a protostar with accretion rate of $\dot{M}_*=10^{-5}M_\odot{\rm yr}^{-1}$. The left open circle represents the beginning of central Deuterium burning (state a). On the second circle, radiative barrier appears (state b) and Deuterium shell burning begins. Beyond the mass, the radius swells dramatically. Further increase in the mass increases the gravity and the radius shrinks. At the third open circle, Hydrogen begins to burn at the center of the protostar.

In the HR diagram, the mass-radius relation for the protostars shown in Figure 4.23(right) corresponds to the dotted line (Fig.4.24). Observed premain-sequence stars: T Tauri stars ( $M_* \mbox{\raisebox{0.3ex}{$<$}\hspace{-1.1em}
\raisebox{-0.7ex}{$\sim$}} 2M_\odot$) and Herbig Ae/Be stars ( $2M_\odot\mbox{\raisebox{0.3ex}{$<$}\hspace{-1.1em}
\raisebox{-0.7ex}{$\sim$}} ...
...\mbox{\raisebox{0.3ex}{$<$}\hspace{-1.1em}
\raisebox{-0.7ex}{$\sim$}} 5M_\odot$) distribute below the dotted line. Thus, this line is often called the birth line for stars. In Figure 4.24, evolutionary tracks without mass accretion $\dot{M}=0$ are plotted from the birth line. Modulation of the birth line is related to the swell and shrink of radius shown in Figure 4.23. That is, increase of luminosity around $M\simeq 1M_\odot$ corresponds to the first swell after Deuterium burning begins (between the first and second open circles). Another increase around $M\simeq 2.5M_\odot$ is related to the Deuterium shell burning (after the second open circle).

Figure 4.23: Schematic view of Deuterium burning in protostars (left). Evolution of the radius of a protostar is plotted against its mass (right). Both figures are taken from Palla & Stahler (1990).
\begin{figure}\begin{center}
\centering\leavevmode
\epsfxsize =.45\columnwidth...
...=.45\columnwidth \epsfbox{eps/Palla_Stahler_90_M-R.ps}
\end{center}
\end{figure}

Figure 4.24: Premain-sequence evolutionary path in the H-R diagram. (Palla & Staler 1993). Evolution of stars with different masses from $0.6 M_\odot $ to $6 M_\odot $ are calculated from the birth line (dotted line).
\begin{figure}\begin{center}
\centering\leavevmode
\epsfxsize =.45\columnwidth \epsfbox{eps/Palla_Stahler_93_HR.ps}
\end{center}
\end{figure}

Kohji Tomisaka 2009-12-10