In the proton-proton chain reaction in the Sun, the first step is $$p + p \rightarrow \; ^2_2{\rm He} .$$ After this, the most likely thing to happen next is that the reverse reaction occurs: the $^2_2{\rm He}$ splits up and we are back where we started. Sometimes, very much more rarely, there is a weak beta+ decay process which converts a proton to a neutron and then we get a deuteron and further reactions can occur. I know how to find the rate for the latter process, but I would like to find the rate for the first process.

My question is: In the solar core, how many times does the fuse-fission sequence $p+p \rightarrow {\rm He} \rightarrow p+p$ occur, on average, for any given proton, before we get the weak decay and a deuteron? Even an order-of-magnitude estimate from rough cross section formula would be helpful. (I have tried to find it with no luck yet.)

Edit added 19 Feb 2018: My research so far has yielded numbers between $10^{15}$ and $10^{28}$. So anyone who can pin it down better than that range (13 orders of magnitude) is making progress!

  • 2
    $\begingroup$ Lots and lots. :) I've never seen a raw number, or probability, just the estimated half-life of 1 billion years for a solar core proton to be converted to He-4. Wikipedia says that the deuteron converts to He-3 in about 4 seconds, and He-3 to He-4 can proceed in various ways, taking about 400 years. $\endgroup$
    – PM 2Ring
    Commented Feb 12, 2019 at 11:54

1 Answer 1


There is a nice analysis by Gillian Knapp here. She argues that it's natural to define the cross-section for collision as the cross-section for an S-wave collision, and this comes out to be four orders of magnitude bigger than the geometrical cross-section. Using this along with the number density of protons in the core, she gets a collision rate of $\sim 10^{12}\ \text{s}^{-1}$. This is the rate of assaults on the Coulomb barrier, not the rate of close collisions in which we actually penetrate the barrier. The observed fusion rate, based on how much of the sun's hydrogen remains unburned, is $5\times10^{-18}\ \text{s}^{-1}$. The WKB barrier penetration probability for S waves is $\sim e^{-16}\sim10^{-7}$. So we have

$\frac{\text{rate of barrier-penetrating collisions}}{\text{rate of fusion}} = \frac{(\text{rate of S-wave collisions})(\text{WKB probability})}{\text{rate of fusion}}\sim 10^{23}.$

IIRC this kind of thing is tricky to estimate because you're looking at the overlap of two thin tails: the upper tail of the Maxwellian distribution and the lower tail of the cross-section as a function of energy. There is some energy at which the product of these two functions has a pronounced, narrow peak, and it's this energy at which you should be estimating the WKB probability, not the typical thermal energy.

  • 1
    $\begingroup$ Thanks again. At the moment I think your number is in the right ball-park for a general elastic nuclear process between the protons, and may be somewhat of an overestimate for a process where the protons linger close together, but I am not yet very confident. I will wait a little longer on the off-chance of any more answers before ticking "accept". $\endgroup$ Commented Feb 21, 2019 at 19:14

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.