Gamow factor explained

The Gamow factor, Sommerfeld factor or Gamow–Sommerfeld factor,[1] named after its discoverer George Gamow or after Arnold Sommerfeld, is a probability factor for two nuclear particles' chance of overcoming the Coulomb barrier in order to undergo nuclear reactions, for example in nuclear fusion. By classical physics, there is almost no possibility for protons to fuse by crossing each other's Coulomb barrier at temperatures commonly observed to cause fusion, such as those found in the Sun. When George Gamow instead applied quantum mechanics to the problem, he found that there was a significant chance for the fusion due to tunneling.

The probability of two nuclear particles overcoming their electrostatic barriers is given by the following equation:[2]

PG(E)=

-\sqrt{{EG
e

/{E}}}

where

EG

is the Gamow energy,

EG\equiv2mrc2(\pi\alphaZa

2
Z
b)

Here,

mr=

mamb
ma+mb
is the reduced mass of the two particles. The constant

\alpha

is the fine-structure constant,

c

is the speed of light, and

Za

and

Zb

are the respective atomic numbers of each particle.

While the probability of overcoming the Coulomb barrier increases rapidly with increasing particle energy, for a given temperature, the probability of a particle having such an energy falls off very fast, as described by the Maxwell–Boltzmann distribution. Gamow found that, taken together, these effects mean that for any given temperature, the particles that fuse are mostly in a temperature-dependent narrow range of energies known as the Gamow window.

Derivation

Gamow[3] first solved the one-dimensional case of quantum tunneling using the WKB approximation. Considering a wave function of a particle of mass m, we take area 1 to be where a wave is emitted, area 2 the potential barrier which has height V and width l (at

0<x<l

), and area 3 its other side, where the wave is arriving, partly transmitted and partly reflected. For a wave number k and energy E we get:

\Psi1=Aei(kx+\alpha)e-i{Et/{\hbar}}

\Psi2=B1e-k'x+B2ek'x

\Psi3=(C1e-i(kx+\beta)+C2ei(kx+\beta'))e-i{Et/{\hbar}}

where

k=\sqrt{2mE}

and k' = \sqrt.This is solved for given A and α by taking the boundary conditions at the both barrier edges, at

x=0

and

x=l

, where both

\Psi

and its derivative must be equal on both sides.For

k'l\gg1

, this is easily solved by ignoring the time exponential and considering the real part alone (the imaginary part has the same behavior). We get, up to factors depending on the phases which are typically of order 1, and up to factors of the order of /=\sqrt (assumed not very large, since V is greater than E not marginally):

B1,B2A

C1,C2

1A
2
k'
k

ek'l

Next Gamow modeled the alpha decay as a symmetric one-dimensional problem, with a standing wave between two symmetric potential barriers at

q0<x<q0+l

and

-(q0+l)<x<-q0

, and emitting waves at both outer sides of the barriers.Solving this can in principle be done by taking the solution of the first problem, translating it by

q0

and gluing it to an identical solution reflected around

x=0

.

Due to the symmetry of the problem, the emitting waves on both sides must have equal amplitudes (A), but their phases (α) may be different. This gives a single extra parameter; however, gluing the two solutions at

x=0

requires two boundary conditions (for both the wave function and its derivative), so in general there is no solution. In particular, re-writing

\Psi3

(after translation by

q0

) as a sum of a cosine and a sine of

kx

, each having a different factor that depends on k and α, the factor of the sine must vanish, so that the solution can be glued symmetrically to its reflection. Since the factor is in general complex (hence its vanishing imposes two constraints, representing the two boundary conditions), this can in general be solved by adding an imaginary part of k, which gives the extra parameter needed. Thus E will have an imaginary part as well.

The physical meaning of this is that the standing wave in the middle decays; the emitted waves newly emitted have therefore smaller amplitudes, so that their amplitude decays in time but grows with distance. The decay constant, denoted &lambda;, is assumed small compared to

E/\hbar

.

&lambda; can be estimated without solving explicitly, by noting its effect on the probability current conservation law. Since the probability flows from the middle to the sides, we have:

\partial
\partialt
(q0+l)
\int
-(q0+l)

\Psi*\Psidx=2 ⋅

\hbar
2mi
*
\left(\Psi
1
\partial\Psi1
\partialx

-\Psi1

\partial
*
\Psi
1
\partialx

\right),

Note the factor of 2 is due to having two emitted waves.

Taking

\Psi\sime

, this gives:

λ ⋅

1
4

2(q0+l)A2

k'2
k2

e2k'l ≈ 2

\hbar
m

A2k,

Since the quadratic dependence in

k'l

is negligible relative to its exponential dependence, we may write:

λ

\hbark
m(q0+l)
k2
k'2

e-2k'l

Remembering the imaginary part added to k is much smaller than the real part, we may now neglect it and get:

λ

\hbark
m2(q0+l)

⋅ 8

E
V-E

e-2\sqrt{2m(V-E)l/\hbar}

Note that
\hbark
m
is the particle velocity, so the first factor is the classical rate by which the particle trapped between the barriers hits them.

Finally, moving to the three-dimensional problem, the spherically symmetric Schrödinger equation reads (expanding the wave function

\psi(r,\theta,\phi)=\chi(r)u(\theta,\phi)

in spherical harmonics and looking at the n-th term):
\hbar2\left(
2m
d2\chi
dr2

+

2
r
d\chi
dr

\right)=\left(V(r)+

\hbar2
2m
n(n+1)
r2

-E\right)\chi

Since

n>0

amounts to enlarging the potential, and therefore substantially reducing the decay rate (given its exponential dependence on

\sqrt{V-E}

), we focus on

n=0

, and get a very similar problem to the previous one with

\chi(r)=\Psi(r)/r

, except that now the potential as a function of r is not a step function.

The main effect of this on the amplitudes is that we must replace the argument in the exponent, taking an integral of 2\sqrt/\hbar over the distance where

V(r)>E

rather than multiplying by l. We take the Coulomb potential:

V(r)=

z(Z-z)e2
4\pi\varepsilon0r
where

\varepsilon0

is the vacuum electric permittivity, e the electron charge, z = 2 is the charge number of the alpha particle and Z the charge number of the nucleus (Z-z after emitting the particle). The integration limits are then

r2=

z(Z-z)e2
4\pi\varepsilon0E
, where we assume the nuclear potential energy is still relatively small, and

r1

, which is where the nuclear negative potential energy is large enough so that the overall potential is smaller than E. Thus, the argument of the exponent in &lambda; is:
2\sqrt{2mE
} \int_^ \sqrt \, dr = 2\frac \int_^ \sqrt \,dr This can be solved by substituting

t=\sqrt{r/r2}

and then

t=cos(\theta)

and solving for θ, giving:

2 ⋅

r
2\sqrt{2mE
} \cdot(\cos^(\sqrt) - \sqrt\sqrt) = 2\frac \cdot(\cos^(\sqrt) - \sqrt\sqrt)where

x=r1/r2

.Since x is small, the x-dependent factor is of order 1.

Gamow assumed

x\ll1

, thus replacing the x-dependent factor by

\pi/2

, giving:

λ\sim

-\sqrt{{Eg
e

/{E}}}

with:

Eg=

2\pi2m\left[z(Z-z)e2\right]2
4\pi\varepsilon0\hbar2
which is the same as the formula given in the beginning of the article with

Za=z

,

Zb=Z-z

and the fine-structure constant

\alpha=

e2
4\pi\varepsilon0\hbarc
.

For a radium alpha decay, Z = 88, z = 2 and m = 4mp, EG is approximately 50 GeV. Gamow calculated the slope of

log(λ)

with respect to E at an energy of 5 MeV to be ~ 1014 J−1, compared to the experimental value of .

External links

Notes and References

  1. Yoon, Jin-Hee . Wong, Cheuk-Yin . Relativistic Modification of the Gamow Factor . February 9, 2008 . nucl-th/9908079 . 10.1103/PhysRevC.61.044905 . 61 . Physical Review C . 2000PhRvC..61d4905Y.
  2. Web site: Dept. Physics & Astronomy University College London. Nuclear reactions in stars. 2014-11-12. 2017-01-15. https://web.archive.org/web/20170115214447/http://zuserver2.star.ucl.ac.uk/~idh/PHAS2112/Lectures/Current/Part7.pdf. dead.
  3. http://web.ihep.su/dbserv/compas/src/gamow28/eng.pdf Quantum Theory of the Atomic Nucleus, G. Gamow