Appendix A — Mathematical Foundations for Nuclear Physics

This appendix collects the mathematical results used throughout the textbook. It is a reference, not a derivation — proofs and physical motivation appear in the chapters where each tool is first introduced (especially Chapter 5). When you are working a problem and need a Clebsch-Gordan coefficient, a spherical harmonic, or the asymptotic form of a Bessel function, this is where you look.

A.1 Angular Momentum Algebra

A.1.1 Orbital Angular Momentum

The orbital angular momentum operator $\hat{\mathbf{L}}$ satisfies the commutation relations

$$[\hat{L}_i, \hat{L}_j] = i\hbar\, \epsilon_{ijk}\, \hat{L}_k$$

with simultaneous eigenstates $|l, m_l\rangle$:

$$\hat{\mathbf{L}}^2 |l, m_l\rangle = \hbar^2 l(l+1)|l, m_l\rangle, \qquad \hat{L}_z |l, m_l\rangle = \hbar m_l |l, m_l\rangle$$

where $l = 0, 1, 2, \ldots$ and $m_l = -l, -l+1, \ldots, l$. The spectroscopic labeling is $s, p, d, f, g, h, i, \ldots$ for $l = 0, 1, 2, 3, 4, 5, 6, \ldots$

The raising and lowering operators $\hat{L}_\pm = \hat{L}_x \pm i\hat{L}_y$ act as:

$$\hat{L}_\pm |l, m_l\rangle = \hbar\sqrt{l(l+1) - m_l(m_l \pm 1)}\; |l, m_l \pm 1\rangle$$

The parity of an orbital angular momentum state is $(-1)^l$.

A.1.2 Spin Angular Momentum

Nucleons are spin-$\frac{1}{2}$ fermions. The spin operator $\hat{\mathbf{S}}$ obeys the same algebra:

$$[\hat{S}_i, \hat{S}_j] = i\hbar\, \epsilon_{ijk}\, \hat{S}_k$$

with $\hat{\mathbf{S}}^2 |s, m_s\rangle = \hbar^2 s(s+1)|s, m_s\rangle$ and $\hat{S}_z |s, m_s\rangle = \hbar m_s |s, m_s\rangle$, where $s = \frac{1}{2}$ and $m_s = \pm\frac{1}{2}$ for a single nucleon. In matrix form (Pauli matrices):

$$\hat{S}_x = \frac{\hbar}{2}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \qquad \hat{S}_y = \frac{\hbar}{2}\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}, \qquad \hat{S}_z = \frac{\hbar}{2}\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$$

Useful identities for Pauli matrices $\boldsymbol{\sigma} = (\sigma_x, \sigma_y, \sigma_z)$:

• $\sigma_i \sigma_j = \delta_{ij} \mathbb{I} + i\epsilon_{ijk}\sigma_k$
• $(\boldsymbol{\sigma} \cdot \mathbf{A})(\boldsymbol{\sigma} \cdot \mathbf{B}) = \mathbf{A}\cdot\mathbf{B}\;\mathbb{I} + i\boldsymbol{\sigma}\cdot(\mathbf{A}\times\mathbf{B})$
• $\text{Tr}(\sigma_i) = 0$, $\text{Tr}(\sigma_i \sigma_j) = 2\delta_{ij}$

A.1.3 Total Angular Momentum and Coupling Rules

For a single nucleon, the total angular momentum is $\hat{\mathbf{j}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$. The coupled states $|j, m_j\rangle$ are labeled by $j = l \pm \frac{1}{2}$ (for $l \geq 1$; for $l = 0$, only $j = \frac{1}{2}$). In nuclear spectroscopic notation, an orbit is labeled $nl_j$, e.g., $1d_{5/2}$ means $n = 1$, $l = 2$, $j = \frac{5}{2}$.

General coupling of two angular momenta: Given $\hat{\mathbf{J}} = \hat{\mathbf{j}}_1 + \hat{\mathbf{j}}_2$, the coupled states $|j_1, j_2; J, M\rangle$ have quantum numbers:

$$J = |j_1 - j_2|, \; |j_1 - j_2| + 1, \; \ldots, \; j_1 + j_2$$

$$M = m_1 + m_2 = -J, -J+1, \ldots, J$$

The number of coupled states equals the number of uncoupled states: $(2j_1+1)(2j_2+1)$.

L-S coupling (Russell-Saunders): First couple orbital angular momenta to $\mathbf{L} = \sum_i \mathbf{l}_i$ and spins to $\mathbf{S} = \sum_i \mathbf{s}_i$, then couple $\mathbf{J} = \mathbf{L} + \mathbf{S}$. Used when the residual interaction is stronger than spin-orbit coupling (light nuclei, $A \lesssim 20$).

j-j coupling: First couple each nucleon's $\mathbf{l}_i + \mathbf{s}_i = \mathbf{j}_i$, then couple individual $\mathbf{j}_i$ to total $\mathbf{J} = \sum_i \mathbf{j}_i$. Dominant in nuclear physics because the spin-orbit interaction is strong.

Isospin coupling: Proton ($t_z = -\frac{1}{2}$) and neutron ($t_z = +\frac{1}{2}$) form isospin doublet. Coupling follows identical angular momentum algebra: $\hat{\mathbf{T}} = \sum_i \hat{\mathbf{t}}_i$, with $T = |T_z|, |T_z|+1, \ldots$ and $T_z = \frac{1}{2}(N - Z)$.

A.2 Clebsch-Gordan Coefficients

A.2.1 Definition

The Clebsch-Gordan (CG) coefficient $\langle j_1, m_1; j_2, m_2 | J, M\rangle$ relates the uncoupled and coupled angular momentum bases:

$$|j_1, j_2; J, M\rangle = \sum_{m_1, m_2} \langle j_1, m_1; j_2, m_2 | J, M\rangle\; |j_1, m_1\rangle |j_2, m_2\rangle$$

The inverse relation is:

$$|j_1, m_1\rangle |j_2, m_2\rangle = \sum_{J} \langle j_1, m_1; j_2, m_2 | J, M\rangle\; |j_1, j_2; J, M\rangle$$

The CG coefficient vanishes unless: 1. $M = m_1 + m_2$ 2. $|j_1 - j_2| \leq J \leq j_1 + j_2$ (triangle inequality) 3. $j_1 + j_2 + J$ is an integer

A.2.2 Symmetry Properties

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{j_1 + j_2 - J}\langle j_2, m_2; j_1, m_1 | J, M\rangle$$

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{j_1 + j_2 - J}\langle j_1, -m_1; j_2, -m_2 | J, -M\rangle$$

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{j_1 - m_1}\sqrt{\frac{2J+1}{2j_2+1}}\;\langle j_1, m_1; J, -M | j_2, -m_2\rangle$$

Orthogonality relations:

$$\sum_{m_1, m_2} \langle j_1, m_1; j_2, m_2 | J, M\rangle \langle j_1, m_1; j_2, m_2 | J', M'\rangle = \delta_{JJ'}\delta_{MM'}$$

$$\sum_{J, M} \langle j_1, m_1; j_2, m_2 | J, M\rangle \langle j_1, m_1'; j_2, m_2' | J, M\rangle = \delta_{m_1 m_1'}\delta_{m_2 m_2'}$$

A.2.3 Tables of Clebsch-Gordan Coefficients

Table A.1: $j_1 = \frac{1}{2} \otimes j_2 = \frac{1}{2}$ (proton-neutron isospin, two-nucleon spin coupling)

The spin-triplet ($J = 1$) states are symmetric; the spin-singlet ($J = 0$) state is antisymmetric under particle exchange.

Table A.2: $j_1 \otimes j_2 = \frac{1}{2}$ (coupling single-particle $j$ with nucleon spin — used in $l \otimes s$ coupling)

For $j_1 = l$ (integer) and $j_2 = \frac{1}{2}$:

These are used extensively in the shell model (Chapter 6) and electromagnetic transitions (Chapter 9).

Table A.3: $j_1 = 1 \otimes j_2 = 1$ (vector coupling, relevant for $d$-wave, quadrupole operators)

Only $M = m_1 + m_2$ entries are shown (others vanish).

A.3 Wigner 3j and 6j Symbols

A.3.1 Wigner 3j Symbol

The Wigner 3j symbol is a symmetrized form of the Clebsch-Gordan coefficient:

$$\begin{pmatrix} j_1 & j_2 & J \\ m_1 & m_2 & M \end{pmatrix} = \frac{(-1)^{j_1 - j_2 + M}}{\sqrt{2J+1}}\;\langle j_1, m_1; j_2, m_2 | J, -M\rangle$$

The inverse relation is:

$$\langle j_1, m_1; j_2, m_2 | J, M\rangle = (-1)^{-j_1 + j_2 - M}\sqrt{2J+1}\begin{pmatrix} j_1 & j_2 & J \\ m_1 & m_2 & -M \end{pmatrix}$$

Selection rules (the 3j symbol vanishes unless): 1. $m_1 + m_2 + M = 0$ 2. Triangle condition: $|j_1 - j_2| \leq J \leq j_1 + j_2$ 3. $j_1 + j_2 + J$ is a non-negative integer

Symmetry properties: - Even permutation of columns: unchanged - Odd permutation of columns: factor $(-1)^{j_1 + j_2 + J}$ - Sign reversal of all $m$ values: factor $(-1)^{j_1 + j_2 + J}$

$$\begin{pmatrix} j_1 & j_2 & J \\ m_1 & m_2 & M \end{pmatrix} = \begin{pmatrix} j_2 & J & j_1 \\ m_2 & M & m_1 \end{pmatrix} = \begin{pmatrix} J & j_1 & j_2 \\ M & m_1 & m_2 \end{pmatrix}$$

$$= (-1)^{j_1+j_2+J}\begin{pmatrix} j_2 & j_1 & J \\ m_2 & m_1 & M \end{pmatrix} = (-1)^{j_1+j_2+J}\begin{pmatrix} j_1 & j_2 & J \\ -m_1 & -m_2 & -M \end{pmatrix}$$

Orthogonality:

$$\sum_{m_1, m_2}\begin{pmatrix} j_1 & j_2 & J \\ m_1 & m_2 & M \end{pmatrix}\begin{pmatrix} j_1 & j_2 & J' \\ m_1 & m_2 & M' \end{pmatrix} = \frac{\delta_{JJ'}\delta_{MM'}}{2J+1}$$

A.3.2 Wigner 6j Symbol

The Wigner 6j symbol arises in the recoupling of three angular momenta. If we can couple $(j_1, j_2)$ to $j_{12}$, then couple $j_{12}$ with $j_3$ to $J$, or alternatively couple $(j_2, j_3)$ to $j_{23}$, then couple $j_1$ with $j_{23}$ to $J$, the unitary transformation between these two coupling schemes involves the 6j symbol:

$$\begin{Bmatrix} j_1 & j_2 & j_{12} \\ j_3 & J & j_{23} \end{Bmatrix}$$

Relation to CG coefficients (Racah formula):

$$\begin{Bmatrix} j_1 & j_2 & j_{12} \\ j_3 & J & j_{23} \end{Bmatrix} = \sum_{\text{all }m} (-1)^{j_3+J+j_{23}+m_{23}} \begin{pmatrix} j_1 & j_2 & j_{12} \\ m_1 & m_2 & -m_{12} \end{pmatrix} \begin{pmatrix} j_1 & j_{23} & J \\ m_1 & m_{23} & -M \end{pmatrix} \begin{pmatrix} j_2 & j_3 & j_{23} \\ m_2 & m_3 & -m_{23} \end{pmatrix} \begin{pmatrix} j_{12} & j_3 & J \\ m_{12} & m_3 & -M \end{pmatrix}$$

Selection rules: The 6j symbol vanishes unless all four triangle conditions are simultaneously satisfied: $\Delta(j_1, j_2, j_{12})$, $\Delta(j_1, j_{23}, J)$, $\Delta(j_2, j_3, j_{23})$, and $\Delta(j_{12}, j_3, J)$.

Important special cases:

When one argument is zero:

$$\begin{Bmatrix} j_1 & j_2 & J \\ 0 & J & j_2 \end{Bmatrix} = \frac{(-1)^{j_1+j_2+J}}{\sqrt{(2j_2+1)(2J+1)}}$$

The 6j symbol appears in nuclear physics in: - Particle-hole conjugation (Chapter 7) - Reduced matrix elements of tensor products (Chapters 9, 15) - Two-body matrix elements in the shell model (Chapter 7) - Angular correlation functions (Chapter 9)

A.4 Spherical Harmonics

A.4.1 Definition and Properties

The spherical harmonics $Y_l^m(\theta, \phi)$ are the angular part of the solutions to Laplace's equation in spherical coordinates, and the eigenfunctions of $\hat{\mathbf{L}}^2$ and $\hat{L}_z$:

$$Y_l^m(\theta, \phi) = (-1)^m \sqrt{\frac{(2l+1)}{4\pi}\frac{(l-m)!}{(l+m)!}}\; P_l^m(\cos\theta)\; e^{im\phi}$$

where $P_l^m$ is the associated Legendre function (using the Condon-Shortley phase convention).

Orthonormality:

$$\int_0^{2\pi}\int_0^{\pi} Y_l^m(\theta,\phi)^* Y_{l'}^{m'}(\theta,\phi)\;\sin\theta\;d\theta\;d\phi = \delta_{ll'}\delta_{mm'}$$

Completeness:

$$\sum_{l=0}^{\infty}\sum_{m=-l}^{l} Y_l^m(\theta,\phi)^* Y_l^m(\theta',\phi') = \delta(\cos\theta - \cos\theta')\delta(\phi - \phi')$$

Complex conjugation:

$$Y_l^m(\theta,\phi)^* = (-1)^m Y_l^{-m}(\theta,\phi)$$

Parity:

$$Y_l^m(\pi - \theta, \phi + \pi) = (-1)^l Y_l^m(\theta, \phi)$$

Addition theorem:

$$P_l(\cos\gamma) = \frac{4\pi}{2l+1}\sum_{m=-l}^{l} Y_l^m(\theta_1,\phi_1)^* Y_l^m(\theta_2,\phi_2)$$

where $\gamma$ is the angle between directions $(\theta_1,\phi_1)$ and $(\theta_2,\phi_2)$.

A.4.2 Explicit Forms

Table A.4: Spherical harmonics through $l = 3$

A.4.3 Product and Integration Rules

The integral of three spherical harmonics (Gaunt integral) is:

$$\int Y_{l_1}^{m_1*}(\Omega)\; Y_{l_2}^{m_2}(\Omega)\; Y_{l_3}^{m_3}(\Omega)\; d\Omega = \sqrt{\frac{(2l_1+1)(2l_2+1)(2l_3+1)}{4\pi}} \begin{pmatrix} l_1 & l_2 & l_3 \\ 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} l_1 & l_2 & l_3 \\ -m_1 & m_2 & m_3 \end{pmatrix}$$

This integral vanishes unless $m_1 = m_2 + m_3$, the triangle condition holds, and $l_1 + l_2 + l_3$ is even. The Gaunt integral is central to computing electromagnetic transition matrix elements (Chapter 9) and multipole expansion coefficients.

A.5 Legendre Polynomials and Associated Legendre Functions

A.5.1 Legendre Polynomials

The Legendre polynomials $P_l(x)$ are defined by the Rodrigues formula:

$$P_l(x) = \frac{1}{2^l l!}\frac{d^l}{dx^l}(x^2 - 1)^l$$

First several:

Orthogonality:

$$\int_{-1}^{1} P_l(x) P_{l'}(x)\; dx = \frac{2}{2l+1}\delta_{ll'}$$

Recursion relation:

$$(l+1)P_{l+1}(x) = (2l+1)x P_l(x) - l P_{l-1}(x)$$

Special values: $P_l(1) = 1$, $P_l(-1) = (-1)^l$, $P_l(0) = 0$ for odd $l$, $P_{2n}(0) = (-1)^n \binom{2n}{n}/2^{2n}$.

Nuclear physics application: The differential cross section for Rutherford scattering and nuclear reactions is expanded in Legendre polynomials:

$$\frac{d\sigma}{d\Omega}(\theta) = \sum_l A_l P_l(\cos\theta)$$

This is the partial wave expansion in angle (Chapter 17).

A.5.2 Associated Legendre Functions

$$P_l^m(x) = (-1)^m (1 - x^2)^{m/2}\frac{d^m}{dx^m}P_l(x), \qquad 0 \leq m \leq l$$

For negative $m$: $P_l^{-m}(x) = (-1)^m \frac{(l-m)!}{(l+m)!}P_l^m(x)$.

A.6 Spherical Bessel Functions

A.6.1 Definitions

The spherical Bessel functions of the first and second kind are:

$$j_l(x) = \sqrt{\frac{\pi}{2x}}\; J_{l+1/2}(x), \qquad n_l(x) = (-1)^{l+1}\sqrt{\frac{\pi}{2x}}\; J_{-l-1/2}(x)$$

where $J_\nu$ is the ordinary Bessel function. They arise as the radial solutions of the free-particle Schrodinger equation in spherical coordinates (or equivalently, the Helmholtz equation).

Low-order explicit forms:

$$j_0(x) = \frac{\sin x}{x}$$

$$j_1(x) = \frac{\sin x}{x^2} - \frac{\cos x}{x}$$

$$j_2(x) = \left(\frac{3}{x^3} - \frac{1}{x}\right)\sin x - \frac{3}{x^2}\cos x$$

$$n_0(x) = -\frac{\cos x}{x}$$

$$n_1(x) = -\frac{\cos x}{x^2} - \frac{\sin x}{x}$$

$$n_2(x) = -\left(\frac{3}{x^3} - \frac{1}{x}\right)\cos x - \frac{3}{x^2}\sin x$$

A.6.2 Asymptotic Forms

Small argument ($x \to 0$):

$$j_l(x) \to \frac{x^l}{(2l+1)!!}, \qquad n_l(x) \to -\frac{(2l-1)!!}{x^{l+1}}$$

where $(2l+1)!! = 1 \cdot 3 \cdot 5 \cdots (2l+1)$.

Large argument ($x \to \infty$):

$$j_l(x) \to \frac{1}{x}\sin\left(x - \frac{l\pi}{2}\right), \qquad n_l(x) \to -\frac{1}{x}\cos\left(x - \frac{l\pi}{2}\right)$$

A.6.3 Spherical Hankel Functions

The spherical Hankel functions (outgoing and incoming spherical waves) are:

$$h_l^{(1)}(x) = j_l(x) + i n_l(x), \qquad h_l^{(2)}(x) = j_l(x) - i n_l(x)$$

For large $x$: $h_l^{(1)}(x) \to (-i)^{l+1} e^{ix}/x$ (outgoing wave) and $h_l^{(2)}(x) \to (i)^{l+1} e^{-ix}/x$ (incoming wave).

Nuclear physics applications: Spherical Bessel functions appear in: - Partial wave decomposition of scattering states (Chapter 17) - Free-particle radial wavefunctions inside and outside the nuclear potential (Chapters 3, 6) - Born approximation matrix elements (Chapter 19) - Nuclear form factors from electron scattering (Chapter 2)

A.7 Laguerre Polynomials and Associated Laguerre Polynomials

A.7.1 Definition

The Laguerre polynomials $L_n(x)$ are defined by the Rodrigues formula:

$$L_n(x) = \frac{e^x}{n!}\frac{d^n}{dx^n}(x^n e^{-x})$$

The associated Laguerre polynomials are:

$$L_n^k(x) = \frac{d^k}{dx^k}L_{n+k}(x) = (-1)^k \frac{(n+k)!}{n!\; k!}\; {}_1F_1(-n; k+1; x)$$

First several (standard Laguerre, $k = 0$):

Orthogonality (associated):

$$\int_0^{\infty} x^k e^{-x} L_n^k(x) L_m^k(x)\; dx = \frac{(n+k)!}{n!}\;\delta_{nm}$$

Nuclear physics application: The radial wavefunctions of the three-dimensional harmonic oscillator potential (used as a starting approximation for the nuclear mean field) are:

$$R_{nl}(r) = N_{nl}\; \left(\frac{r}{b}\right)^l\; L_{n}^{l+1/2}\left(\frac{r^2}{b^2}\right)\; e^{-r^2/(2b^2)}$$

where $b = \sqrt{\hbar/(m\omega)}$ is the oscillator length parameter and $n$ is the number of radial nodes (Chapter 6).

A.8 Gamma Function and Error Function

A.8.1 Gamma Function

The gamma function generalizes the factorial to complex arguments:

$$\Gamma(z) = \int_0^{\infty} t^{z-1} e^{-t}\; dt, \qquad \text{Re}(z) > 0$$

Key properties: - $\Gamma(n+1) = n!$ for non-negative integers $n$ - $\Gamma(z+1) = z\,\Gamma(z)$ (recursion) - $\Gamma(\frac{1}{2}) = \sqrt{\pi}$ - $\Gamma(n + \frac{1}{2}) = \frac{(2n)!}{4^n n!}\sqrt{\pi}$ - Reflection formula: $\Gamma(z)\Gamma(1-z) = \frac{\pi}{\sin(\pi z)}$ - Duplication formula: $\Gamma(z)\Gamma(z + \frac{1}{2}) = \frac{\sqrt{\pi}}{2^{2z-1}}\Gamma(2z)$ - Double factorial: $(2n-1)!! = \frac{2^n}{\sqrt{\pi}}\Gamma(n + \frac{1}{2})$

Table A.5: Useful values

A.8.2 Error Function

$$\text{erf}(x) = \frac{2}{\sqrt{\pi}}\int_0^x e^{-t^2}\; dt$$

Properties: - $\text{erf}(0) = 0$, $\text{erf}(\infty) = 1$ - $\text{erf}(-x) = -\text{erf}(x)$ (odd function) - Complementary error function: $\text{erfc}(x) = 1 - \text{erf}(x) = \frac{2}{\sqrt{\pi}}\int_x^{\infty} e^{-t^2}\; dt$ - Small-$x$ expansion: $\text{erf}(x) \approx \frac{2x}{\sqrt{\pi}}\left(1 - \frac{x^2}{3} + \frac{x^4}{10} - \cdots\right)$ - Large-$x$ asymptotic: $\text{erfc}(x) \approx \frac{e^{-x^2}}{x\sqrt{\pi}}\left(1 - \frac{1}{2x^2} + \cdots\right)$

The error function appears in Gaussian beam profiles, thermal averaging of reaction rates (Gamow peak width, Chapter 21), and stopping power calculations (Chapter 16).

A.9 The WKB Approximation

A.9.1 General WKB Wavefunctions

The Wentzel-Kramers-Brillouin (WKB) approximation provides semiclassical solutions to the one-dimensional Schrodinger equation $-\frac{\hbar^2}{2m}\psi''(x) + V(x)\psi(x) = E\psi(x)$.

Classically allowed region ($E > V(x)$, local wavenumber $k(x) = \sqrt{2m(E - V(x))}/\hbar$):

$$\psi_{\text{WKB}}(x) = \frac{A}{\sqrt{k(x)}}\;\exp\left(\pm i\int^x k(x')\; dx'\right)$$

Classically forbidden region ($E < V(x)$, local decay constant $\kappa(x) = \sqrt{2m(V(x) - E)}/\hbar$):

$$\psi_{\text{WKB}}(x) = \frac{B}{\sqrt{\kappa(x)}}\;\exp\left(\pm\int^x \kappa(x')\; dx'\right)$$

The WKB approximation is valid when the potential varies slowly on the scale of the local wavelength: $|dk/dx| \ll k^2$, equivalently $|dV/dx| \ll 2m(E-V)^{3/2}/\hbar^2$.

A.9.2 Connection Formulas at Turning Points

At a classical turning point $x_0$ where $E = V(x_0)$, the WKB wavefunctions diverge. The connection formulas, derived from the Airy function solutions near the turning point, are:

Right-hand turning point (classically allowed for $x < x_0$, forbidden for $x > x_0$):

$$\frac{2}{\sqrt{k(x)}}\cos\left(\int_x^{x_0} k\; dx' - \frac{\pi}{4}\right) \longleftrightarrow \frac{1}{\sqrt{\kappa(x)}}\exp\left(-\int_{x_0}^x \kappa\; dx'\right)$$

Left-hand turning point (forbidden for $x < x_0$, allowed for $x > x_0$):

$$\frac{1}{\sqrt{\kappa(x)}}\exp\left(-\int_x^{x_0} \kappa\; dx'\right) \longleftrightarrow \frac{2}{\sqrt{k(x)}}\cos\left(\int_{x_0}^x k\; dx' - \frac{\pi}{4}\right)$$

A.9.3 Tunneling Probability

For a barrier between turning points $r_1$ (inner) and $r_2$ (outer) where $V(r) > E$:

$$T_{\text{WKB}} = \exp\left(-\frac{2}{\hbar}\int_{r_1}^{r_2}\sqrt{2m(V(r) - E)}\; dr\right)$$

This is the Gamow factor. For alpha decay through the Coulomb barrier (Chapter 13), $V(r)$ is the Coulomb potential beyond the nuclear surface:

$$V(r) = \frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 r} \qquad \text{for } r > R$$

with turning points $r_1 = R$ (nuclear radius) and $r_2 = Z_1 Z_2 e^2 / (4\pi\epsilon_0 E)$ (classical turning point). The resulting Gamow factor is:

$$G = \frac{2}{\hbar}\int_R^{r_2}\sqrt{2\mu\left(\frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 r} - E\right)}\; dr = \frac{2Z_1 Z_2 e^2}{\hbar v}\left[\arccos\sqrt{\rho} - \sqrt{\rho(1-\rho)}\right]$$

where $\mu$ is the reduced mass, $v = \sqrt{2E/\mu}$ is the asymptotic velocity, and $\rho = R/r_2$. For $\rho \ll 1$ (the typical case for alpha decay):

$$G \approx \frac{2\pi Z_1 Z_2 e^2}{\hbar v} - 4\sqrt{\frac{2\mu Z_1 Z_2 e^2 R}{4\pi\epsilon_0 \hbar^2}}$$

The first term is the Sommerfeld parameter $2\pi\eta$; the second is a finite-size correction.

A.9.4 Bound-State Quantization (Bohr-Sommerfeld)

For a bound state in a potential well between turning points $x_1$ and $x_2$:

$$\int_{x_1}^{x_2} k(x)\; dx = \left(n + \frac{1}{2}\right)\pi, \qquad n = 0, 1, 2, \ldots$$

This determines the allowed energies semiclassically. Applied to the harmonic oscillator, it gives the exact result $E_n = \hbar\omega(n + \frac{1}{2})$. Applied to nuclear potentials, it provides useful estimates of single-particle energies and helps understand the shell structure (Chapter 6).

A.10 Fourier Transforms

A.10.1 Convention

This textbook uses the symmetric Fourier transform convention:

$$\tilde{f}(k) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} f(x)\; e^{-ikx}\; dx$$

$$f(x) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} \tilde{f}(k)\; e^{ikx}\; dk$$

In three dimensions:

$$\tilde{f}(\mathbf{k}) = \frac{1}{(2\pi)^{3/2}}\int f(\mathbf{r})\; e^{-i\mathbf{k}\cdot\mathbf{r}}\; d^3r$$

A.10.2 Important Transform Pairs

A.10.3 Nuclear Form Factor

The nuclear charge form factor measured in electron scattering (Chapter 2) is the Fourier transform of the charge density:

$$F(q) = \frac{1}{Ze}\int \rho_{\text{ch}}(\mathbf{r})\; e^{i\mathbf{q}\cdot\mathbf{r}}\; d^3r$$

For a spherically symmetric charge distribution, this reduces to:

$$F(q) = \frac{4\pi}{Ze}\int_0^{\infty} \rho_{\text{ch}}(r)\; j_0(qr)\; r^2\; dr = \frac{4\pi}{Ze}\int_0^{\infty} \rho_{\text{ch}}(r)\; \frac{\sin(qr)}{qr}\; r^2\; dr$$

where $q = |\mathbf{q}|$ is the momentum transfer.

A.11 Useful Integrals

A.11.1 Gaussian Integrals

$$\int_{-\infty}^{\infty} e^{-ax^2}\; dx = \sqrt{\frac{\pi}{a}}, \qquad a > 0$$

$$\int_{-\infty}^{\infty} x^{2n} e^{-ax^2}\; dx = \frac{(2n-1)!!}{(2a)^n}\sqrt{\frac{\pi}{a}}$$

$$\int_0^{\infty} x^{2n+1} e^{-ax^2}\; dx = \frac{n!}{2a^{n+1}}$$

$$\int_0^{\infty} x^n e^{-ax}\; dx = \frac{n!}{a^{n+1}} = \frac{\Gamma(n+1)}{a^{n+1}}$$

A.11.2 Angular Momentum Integrals

$$\int Y_{l_1}^{m_1}(\theta,\phi)^* Y_{l_2}^{m_2}(\theta,\phi)\; d\Omega = \delta_{l_1 l_2}\delta_{m_1 m_2}$$

$$\int Y_{l_1}^{m_1*}\; Y_{l_2}^{m_2}\; Y_{l_3}^{m_3}\; d\Omega = \text{(Gaunt integral — see Section A.4.3)}$$

$$\int_0^{\pi} P_l(\cos\theta) P_{l'}(\cos\theta)\;\sin\theta\; d\theta = \frac{2}{2l+1}\delta_{ll'}$$

A.11.3 Radial Integrals

For the harmonic oscillator with $\psi \propto r^l e^{-r^2/(2b^2)}$:

$$\int_0^{\infty} r^{l+2} e^{-r^2/b^2}\; r^2 dr = \frac{b^{l+5}}{2}\;\Gamma\left(\frac{l+5}{2}\right)$$

Coulomb integrals (used in alpha decay and fusion calculations):

$$\int_R^{b}\sqrt{\frac{b}{r} - 1}\; dr = b\left[\arccos\sqrt{\frac{R}{b}} - \sqrt{\frac{R}{b}\left(1 - \frac{R}{b}\right)}\right]$$

where $b = Z_1 Z_2 e^2/(4\pi\epsilon_0 E)$ is the classical turning point.

A.11.4 Miscellaneous Useful Results

Dirac delta representations:

$$\delta(x) = \frac{1}{2\pi}\int_{-\infty}^{\infty} e^{ikx}\; dk = \lim_{\epsilon\to 0}\frac{1}{\pi}\frac{\epsilon}{x^2 + \epsilon^2} = \lim_{N\to\infty}\frac{\sin(Nx)}{\pi x}$$

In spherical coordinates:

$$\delta^3(\mathbf{r} - \mathbf{r}') = \frac{\delta(r-r')}{r^2}\delta(\cos\theta - \cos\theta')\delta(\phi - \phi') = \frac{\delta(r-r')}{r^2}\sum_{l,m}Y_l^{m*}(\Omega')Y_l^m(\Omega)$$

Stirling's approximation:

$$\ln(n!) \approx n\ln n - n + \frac{1}{2}\ln(2\pi n), \qquad n \gg 1$$

Used in nuclear level density calculations (Chapter 18).

Sommerfeld parameter:

$$\eta = \frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 \hbar v} = \frac{Z_1 Z_2 \alpha}{v/c}\;\frac{m c}{\hbar}\cdot\frac{\hbar}{mc} = \frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 \hbar v}$$

This dimensionless parameter characterizes the importance of Coulomb effects in nuclear reactions (Chapters 17, 21). When $\eta \gg 1$, the Coulomb barrier dominates (low-energy reactions between heavy nuclei). When $\eta \ll 1$, the interaction is essentially nuclear (high-energy or light-ion reactions).

A.12 Summary of Notation

For quick reference, the following notation is used throughout this textbook:

This appendix is a reference. The derivations, physical context, and worked examples are in the main text — particularly Chapter 5 (Quantum Mechanics Review), Chapter 6 (Shell Model), Chapter 9 (Electromagnetic Transitions), Chapter 13 (Alpha Decay), and Chapter 17 (Reaction Fundamentals).