One-Boson Exchange Models


Scott Bogner, National Superconducting Cyclotron Laboratory and Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA
Morten Hjorth-Jensen, National Superconducting Cyclotron Laboratory and Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA

Date: 2017

A typical form of the nuclear force

Here we display a typical way to parametrize (non-relativistic expression) the nuclear two-body force in terms of some operators, the central part, the spin-spin part and the central force.

0

< < < ! ! M A T H _ B L O C K

$$ \left. + C_{SL} \left( {1\over m_\alpha r} + {1\over \left( m_\alpha r\right)^2} \right) \mathbf{L}\cdot \mathbf{S} \right\} \frac{e^{-m_\alpha r}}{m_\alpha r} $$

How do we derive such terms? (Note: no isospin dependence and that the above is an approximation)

Nuclear forces

To derive the above famous form of the nuclear force using field theoretical concepts, we will need some elements from relativistic quantum mechanics. These derivations will be given below. The material here gives some background to this. I know that many of you have not taken a course in quantum field theory. I hope however that you can see the basic ideas leading to the famous non-relativistic expressions for the nuclear force.

Furthermore, when we analyze nuclear data, we will actually try to explain properties like spectra, single-particle energies etc in terms of the various terms of the nuclear force. Moreover, many of you will hear about these terms at various talks, workshops, seminars etc. Then, it is good to have an idea of what people actually mean!!

Dramatis Personae

Baryons Mass (MeV) Mesons Mass (MeV)
$p,n$ 938.926 $\pi$ 138.03
$\Lambda$ 1116.0 $\eta$ 548.8
$\Sigma$ 1197.3 $\sigma$ $\approx 550.0$
$\Delta$ 1232.0 $\rho$ 770
$\omega$ 782.6
$\delta$ 983.0
$K$ 495.8
$K^{\star}$ 895.0

Components of the force and quantum numbers

But before we proceed, we will look into specific quantum numbers of the relative system and study expectation vaues of the various terms of

2

< < < ! ! M A T H _ B L O C K

$$ \left. + C_{SL} \left( {1\over m_\alpha r} + {1\over \left( m_\alpha r\right)^2} \right) \mathbf{L}\cdot \mathbf{S} \right\} \frac{e^{-m_\alpha r}}{m_\alpha r} $$

Relative and CoM system, quantum numbers

When solving the scattering equation or solving the two-nucleon problem, it is convenient to rewrite the Schroedinger equation, due to the spherical symmetry of the Hamiltonian, in relative and center-of-mass coordinates. This will also define the quantum numbers of the relative and center-of-mass system and will aid us later in solving the so-called Lippman-Schwinger equation for the scattering problem.

We define the center-of-mass (CoM) momentum as

$$ \mathbf{K}=\sum_{i=1}^A\mathbf{k}_i, $$

with $\hbar=c=1$ the wave number $k_i=p_i$, with $p_i$ the pertinent momentum of a single-particle state. We have also the relative momentum

$$ \mathbf{k}_{ij}=\frac{1}{2}(\mathbf{k}_i-\mathbf{k}_j). $$

We will below skip the indices $ij$ and simply write $\mathbf{k}$

Relative and CoM system, quantum numbers

In a similar fashion we can define the CoM coordinate

$$ \mathbf{R}=\frac{1}{A}\sum_{i=1}^{A}\mathbf{r}_i, $$

and the relative distance

$$ \mathbf{r}_{ij}=(\mathbf{r}_i-\mathbf{r}_j). $$

Relative and CoM system, quantum numbers

With the definitions

$$ \mathbf{K}=\sum_{i=1}^A\mathbf{k}_i, $$

and

$$ \mathbf{k}_{ij}=\frac{1}{2}(\mathbf{k}_i-\mathbf{k}_j). $$

we can rewrite the two-particle kinetic energy (note that we use $\hbar=c=1$ as

$$ \frac{\mathbf{k}_1^2}{2m_n}+\frac{\mathbf{k}_2^2}{2m_n}=\frac{\mathbf{k}^2}{m_n}+\frac{\mathbf{K}^2}{4m_n}, $$

where $m_n$ is the average of the proton and the neutron masses.

Relative and CoM system, quantum numbers

Since the two-nucleon interaction depends only on the relative distance, this means that we can separate Schroedinger's equation in an equation for the center-of-mass motion and one for the relative motion.

With an equation for the relative motion only and a separate one for the center-of-mass motion we need to redefine the two-body quantum numbers.

Previously we had a two-body state vector defined as $|(j_1j_2)JM_J\rangle$ in a coupled basis. We will now define the quantum numbers for the relative motion. Here we need to define new orbital momenta (since these are the quantum numbers which change). We define

$$ \hat{l}_1+\hat{l}_2=\hat{\lambda}=\hat{l}+\hat{L}, $$

where $\hat{l}$ is the orbital momentum associated with the relative motion and $\hat{L}$ the corresponding one linked with the CoM. The total spin $S$ is unchanged since it acts in a different space. We have thus that

$$ \hat{J}=\hat{l}+\hat{L}+\hat{S}, $$

which allows us to define the angular momentum of the relative motion

$$ { \cal J} = \hat{l}+\hat{S}, $$

where ${ \cal J}$ is the total angular momentum of the relative motion.

Models for nuclear forces and derivation of non-relativistic expressions

To describe the interaction between the various baryons and mesons of the previous table we choose the following phenomenological lagrangians for spin $1/2$ baryons

1 4

< < < ! ! M A T H _ B L O C K

$$ {\cal L}_{s} =g^{s}\overline{\Psi}\Psi\phi^{(s)}, $$

and

$$ {\cal L}_{v} =g^{v}\overline{\Psi}\gamma_{\mu}\Psi\phi_{\mu}^{(v)} +g^{t}\overline{\Psi}\sigma^{\mu\nu}\Psi\left (\partial_{\mu}\phi_{\nu}^{(v)} -\partial_{\nu}\phi_{\mu}^{(v)}\right), $$

for pseudoscalar (ps), scalar (s) and vector (v) coupling, respectively. The factors $g^{v}$ and $g^{t}$ are the vector and tensor coupling constants, respectively.

Models for nuclear forces and derivation of non-relativistic expressions

For spin $1/2$ baryons, the fields $\Psi$ are expanded in terms of the Dirac spinors (positive energy solution shown here with $\overline{u}u=1$)

$$ u(k\sigma)=\sqrt{\frac{E(k)+m}{2m}} \left(\begin{array}{c} \chi\\ \\ \frac{\mathbf{\sigma}\mathbf{k}}{E(k)+m}\chi \end{array}\right), $$

with $\chi$ the familiar Pauli spinor and $E(k) =\sqrt{m^2 +|\mathbf{k}|^2}$. The positive energy part of the field $\Psi$ reads

$$ \Psi (x)={\displaystyle \frac{1}{(2\pi )^{3/2}} \sum_{\mathbf{k}\mathbf{\sigma}}u(k\mathbf{\sigma})\exp{-(ikx)}a_{\mathbf{k}\mathbf{\sigma}}}, $$

with $a$ being a fermion annihilation operator.

Models for nuclear forces and derivation of non-relativistic expressions

Expanding the free Dirac spinors in terms of $1/m$ ($m$ is here the mass of the relevant baryon) results, to lowest order, in the familiar non-relativistic expressions for baryon-baryon potentials. The configuration space version of the interaction can be approximated as

1 9

< < < ! ! M A T H _ B L O C K

$$ + C_{SL}\left. \left( {1\over m_\alpha r} + {1\over \left( m_\alpha r\right)^2} \right) \mathbf{L}\cdot \mathbf{S} \right\} \frac{\exp{-(m_\alpha r)}}{m_\alpha r}, $$

where $m_{\alpha}$ is the mass of the relevant meson and $S_{12}$ is the familiar tensor term.

Models for nuclear forces and derivation of non-relativistic expressions

We derive now the non-relativistic one-pion exchange interaction.

Here $p_{1}$, $p_{1}'$, $p_{2}$, $p_{2}'$ and $k=p_{1}-p_{1}'$ denote four-momenta.
The vertices are given by the pseudovector Lagrangian

$$ {\cal L}_{pv}=\frac{f_{\pi}}{m_{\pi}}\overline{\psi}\gamma_{5}\gamma_{\mu} \psi\partial^{\mu}\phi_{\pi}. $$

From the Feynman diagram rules we can write the two-body interaction as

$$ V^{pv}=\frac{f_{\pi}^{2}}{m_{\pi}^{2}}\frac{\overline{u}(p_{1}')\gamma_{5} \gamma_{\mu}(p_{1}-p_{1}')^{\mu}u(p_{1})\overline{u}(p_{2}')\gamma_{5} \gamma_{\nu}(p_{2}'-p_{2})^{\nu}u(p_{2})}{(p_{1}-p_{1}')^{2}-m_{\pi}^{2}}. $$

Models for nuclear forces and derivation of non-relativistic expressions

The factors $p_{1}-p_{1}'=p_{2}'-p_{2}$ are both the four-momentum of the exchanged meson and come from the derivative of the meson field in the interaction Lagrangian. The Dirac spinors obey

$$ \begin{eqnarray*} \gamma_{\mu}p^{\mu}u(p)&=&mu(p) \nonumber \\ \overline{u}(p)\gamma_{\mu}p^{\mu}&=&m\overline{u}(p). \nonumber \end{eqnarray*} $$

Models for nuclear forces and derivation of non-relativistic expressions

Using these relations, together with $\{\gamma_{5},\gamma_{\mu}\}=0$, we find

$$ \begin{eqnarray*} \overline{u}(p_{1}')\gamma_{5}\gamma_{\mu}(p_{1}-p_{1}')^{\mu}u(p_{1}) &=&m\overline{u}(p_{1}')\gamma_{5}u(p_{1})+\overline{u}(p_{1}')\gamma_{\mu} p_{1}'^{\mu}\gamma_{5}u(p_{1}) \nonumber \\ &=&2m\overline{u}(p_{1}')\gamma_{5}u(p_{1}) \nonumber \end{eqnarray*} $$

and

$$ \overline{u}(p_{2}')\gamma_{5}\gamma_{\mu}(p_{2}'-p_{2})^{\mu}= -2m\overline{u}(p_{2}')\gamma_{5}u(p_{1}). $$

Models for nuclear forces and derivation of non-relativistic expressions

We get

$$ V^{pv}=-\frac{f_{\pi}^{2}}{m_{\pi}^{2}}4m^{2}\frac{\overline{u}(p_{1}') \gamma_{5}u(p_{1})\overline{u}(p_{2}')\gamma_{5}u(p_{2})}{(p_{1}-p_{1}') ^{2}-m_{\pi}^{2}}. $$

By inserting expressions for the Dirac spinors, we find

$$ \begin{eqnarray*} \overline{u}(p_{1}')\gamma_{5}u(p_{1})&=&\sqrt{\frac{(E_{1}'+m)(E_{1}+m)} {4m^{2}}}\left(\begin{array}{cc}\chi^{\dagger}&-\frac{\sigma_{1}\cdot{ \bf p_{1}}}{E_{1}' +m}\chi^{\dagger}\end{array}\right)\left(\begin{array}{cc}0&1\\1&0\end{array} \right)\nonumber \\ &&\times \left(\begin{array}{c}\chi\\ \frac{\sigma_{1}\cdot\mathbf{p_{1}}}{E_{1}+m}\chi \end{array}\right) \nonumber \\ &=&\sqrt{\frac{(E_{1}'+m)(E_{1}+m)}{4m^{2}}}\left(\frac{\sigma_{1}\cdot \mathbf{p_{1}}}{E_{1}+m}-\frac{\sigma_{1}\cdot\mathbf{p_{1}'}}{E_{1}'+m}\right) \nonumber \end{eqnarray*} $$

Models for nuclear forces and derivation of non-relativistic expressions

Similarly

$$ \overline{u}(p_{2}')\gamma_{5}u(p_{2})=\sqrt{\frac{(E_{2}'+m)(E_{2}+m)} {4m^{2}}}\left(\frac{\sigma_{2}\cdot \mathbf{p}_{2}}{E_{2}+m}- \frac{\sigma_{2}\cdot\mathbf{p'}_{2}}{E_{2}'+m}\right). $$

In the CM system we have $\mathbf{p}_{2}=-\mathbf{p}_{1}$, $\mathbf{p'}_{2}= -\mathbf{p'}_{1}$ and so $E_{2}=E_{1}$, $E_{2}'=E_{1}'$.
We can then write down the relativistic contribution to the NN potential in the CM system:

$$ \begin{eqnarray} V^{pv}&=&-\frac{f_{\pi}^{2}}{m_{\pi}^{2}}4m^{2}\frac{1}{(p_{1}-p_{1}')^{2}- m_{\pi}^{2}}\frac{(E_{1}+m)(E_{1}'+m)}{4m^{2}} \nonumber \\ &\times&\left(\frac{\sigma_{1}\cdot\mathbf{p}_{1}}{E_{1}+m}-\frac{\sigma_{1} \cdot\mathbf{p'}_{1}}{E_{1}'+m}\right)\left(\frac{\sigma_{2}\cdot\mathbf{p}_{1}} {E_{1}+m}-\frac{\sigma_{2}\cdot\mathbf{p'}_{1}}{E_{1}'+m}\right). \nonumber \end{eqnarray} $$

Models for nuclear forces and derivation of non-relativistic expressions

In the non-relativistic limit we have to lowest order

$$ E_{1}=\sqrt{\mathbf{p}_{1}^{2}+m^{2}}\approx m \approx E_{1}' $$

and then $(p_{1}-p_{1}')^{2}=-\mathbf{k}^{2}$, so we get for the contribution to the NN potential

$$ \begin{eqnarray} V^{pv}&=&-\frac{f_{\pi}^{2}}{m_{\pi}^{2}}4m^{2}\frac{1}{\mathbf{k}^{2}+m^{2}} \frac{2m\cdot 2m}{4m^{2}}\frac{\sigma_{1}}{2m}\cdot(\mathbf{p}_{1}-\mathbf{p'}_{1}) \frac{\sigma_{2}}{2m}\cdot (\mathbf{p}_{1}-\mathbf{p'}_{1}) \nonumber \\ &=&-\frac{f_{\pi}^{2}}{m_{\pi}^{2}} \frac{(\sigma_{1}\cdot\mathbf{k})(\sigma_{2}\cdot\mathbf{k})}{\mathbf{k}^{2}+m_{\pi}^{2}}. \nonumber \end{eqnarray} $$

We have omitted exchange terms and the isospin term $\mathbf{\tau}_1\cdot\mathbf{\tau}_2$.

Models for nuclear forces and derivation of non-relativistic expressions

We have

$$ V^{pv}(k)=-\frac{f_{\pi}^{2}}{m_{\pi}^{2}} \frac{(\sigma_{1}\cdot\mathbf{k})(\sigma_{2}\cdot\mathbf{k})}{\mathbf{k}^{2}+m_{\pi}^{2}}. $$

In coordinate space we have

$$ V^{pv}(r)=\int\frac{d^3k}{(2\pi)^3}e^{i\mathbf{kr}}V^{pv}(k) $$

resulting in

$$ V^{pv}(r)=-\frac{f_{\pi}^{2}}{m_{\pi}^{2}} \sigma_{1}\cdot{\nabla}\sigma_{2}\cdot{\nabla} \int\frac{d^3k}{(2\pi)^3}e^{i\mathbf{kr}}\frac{1}{\mathbf{k}^{2}+m_{\pi}^{2}}. $$

Models for nuclear forces and derivation of non-relativistic expressions

We obtain

$$ V^{pv}(r)=-\frac{f_{\pi}^{2}}{m_{\pi}^{2}}\sigma_{1}\cdot{\nabla}\sigma_{2}\cdot{\nabla}\frac{e^{-m_{\pi}r}}{r}. $$

Carrying out the differentation of

$$ V^{pv}(r)=-\frac{f_{\pi}^{2}}{m_{\pi}^{2}}\sigma_{1}\cdot{\nabla}\sigma_{2}\cdot{\nabla}\frac{e^{-m_{\pi}r}}{r}. $$

we arrive at the famous one-pion exchange potential with central and tensor parts

$$ V(\mathbf{r})= -\frac{f_{\pi}^{2}}{m_{\pi}^{2}}\left\{C_{\sigma}\mathbf{\sigma}_1\cdot\mathbf{\sigma}_2+ C_T \left( 1 + \frac{3}{m_\alpha r} + \frac{3}{\left(m_\alpha r\right)^2}\right) S_{12}(\hat r)\right\}\frac{\exp{-m_\pi r}}{m_\pi r}. $$

For the full potential add the exchange part and the $\mathbf{\tau}_1\cdot\mathbf{\tau}_2$ term as well. (Subtle point: there is a divergence which gets cancelled by using cutoffs) This leads to coefficients $C_{\sigma}$ and $C_T$ which are fitted to data.

Models for nuclear forces and derivation of non-relativistic expressions

When we perform similar non-relativistic expansions for scalar and vector mesons we obtain for the $\sigma$ meson

$$ V^{\sigma}= g_{\sigma NN}^{2}\frac{1}{\mathbf{k}^{2}+m_{\sigma}^{2}}\left (-1+\frac{\mathbf{q}^{2}}{2M_N^2} -\frac{\mathbf{k}^{2}}{8M_N^2}-\frac{\mathbf{LS}}{2M_N^2}\right). $$

We note an attractive central force and spin-orbit force. This term has an intermediate range. We have defined $1/2(p_{1}+p_{1}')=\mathbf{q}$. For the full potential add the exchange part and the isospin dependence as well.

Models for nuclear forces and derivation of non-relativistic expressions

We obtain for the $\omega$ meson

$$ V^{\omega}= g_{\omega NN}^{2}\frac{1}{\mathbf{k}^{2}+m_{\omega}^{2}}\left (1-3\frac{\mathbf{LS}}{2M_N^2}\right). $$

We note a repulsive central force and an attractive spin-orbit force. This term has short range. For the full potential add the exchange part and the isospin dependence as well.

Models for nuclear forces and derivation of non-relativistic expressions

Finally for the $\rho$ meson

$$ V^{\rho}= g_{\rho NN}^{2}\frac{\mathbf{k}^{2}}{\mathbf{k}^{2}+m_{\rho}^{2}}\left ( -2\sigma_{1}\sigma_{2}+S_{12}(\hat{k})\right)\tau_{1}\tau_{2}. $$

We note a tensor force with sign opposite to that of the pion. This term has short range. For the full potential add the exchange part and the isospin dependence as well.