Soft pattern of Rutherford scattering from heavy target mass expansion

Figures(1)

Get Citation
Yu Jia and Jia-Yue Zhang. Soft pattern of Rutherford scattering from heavy target mass expansion[J]. Chinese Physics C. doi: 10.1088/1674-1137/acdb56
Yu Jia and Jia-Yue Zhang. Soft pattern of Rutherford scattering from heavy target mass expansion[J]. Chinese Physics C.  doi: 10.1088/1674-1137/acdb56 shu
Milestone
Received: 2023-06-02
Article Metric

Article Views(1706)
PDF Downloads(76)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Soft pattern of Rutherford scattering from heavy target mass expansion

  • 1. Institute of High Energy Physics and Theoretical Physics Center for Science Facilities, Chinese Academy of Sciences, Beijing 100049, China
  • 2. School of Physics, University of Chinese Academy of Sciences, Beijing 100049, China

Abstract: We investigate the soft behavior of the tree-level Rutherford scattering process. We consider two types of Rutherford scattering processes: One in which a low-energy massless point-like projectile (say, a spin-${1/ 2}$ or spin-$ 0 $ electron) hits a static massive composite target particle carrying various spins (up to spin-$ 2 $), and one where a slowly-moving light projectile hits a heavy static composite target. For the first type, the unpolarized cross sections in the laboratory frame are found to exhibit universal forms in the first two orders of $ 1/M $ expansion yet differ at the next-to-next-to-leading order (though some terms at this order still remain universal or depend on the target spin in a definite manner). For the second type, at the lowest order in electron velocity expansion, through all orders in $ 1/M $, the unpolarized cross section is universal (also not sensitive to the projectile spin). The universality partially breaks down at relative order-$ v^2/M^2 $, though some terms at this order are still universal or depend on the target spin in a specific manner. We also employ the effective field theory approach to reproduce the soft behavior of the differential cross sections for when the target particle is a composite spin-${1/ 2}$ fermion.

    HTML

    I.   INTRODUCTION
    • Rutherford scattering is one of the most classic experiments in the history of physics. Originally, Geiger and Marsden bombarded a gold foil with a nonrelativistic α particle beam in $ 1909 $. Shortly after, in $ 1911 $, Rutherford introduced the revolutionary concept of the atomic nucleus and successfully explained the unexpected large-angle rebouncing events by simply exploiting classical mechanics. Without exaggeration, the Rutherford scattering experiment heralded the advent of nuclear physics and quantum mechanics.

      In the late $ 1950 $s, a new form of the Rutherford scattering experiment conducted by McAllister and Hofstadter, that is, bombarding proton and α targets with a $ 188 $ MeV electron beam, played a pivotal role in unraveling the nuclear structure. Through $ ep $ elastic scattering experiments, the electromagnetic form factors (FFs) of the proton have been measured over a large range of $ Q^2 $(for comprehensive review see [1, 2]). From their profiles at the lower $ Q^2 $ end, one can infer the proton's gross features, such as the charge radius and magnetic dipole [3-5]. It is interesting to note the decade-long puzzle regarding the proton's charge radius [6-8], that is, the five standard discrepancy between the value extracted from $ ep $ elastic scattering/ordinary hydrogen spectra and from muonic hydrogen Lamb shift measurement [5, 9].

      To infer the gross features of a composite nucleus from elastic Rutherford scattering, the exchanged photon necessarily carries a long wavelength and hence bears low resolution. To this purpose, it is helpful to concentrate on the low-energy Rutherford scattering process exemplified by the pioneering Hofstadter's experiment. It is worth mentioning that another basic QED process, Compton scattering, in which a photon beam shines on a composite spinning target particle, can also be used to probe the internal structure of the atomic nucleus [10]. The angular distribution of Compton scattering in the laboratory frame in the soft photon limit has been thoroughly studied by Gell-Mann and Low in the $ 1960 $, which turns out to possess some simple and universal structures. Based on the intuitive multipole expansion picture, one naturally anticipates that the soft limit of Rutherford scattering may also exhibit some universal and recognizable patterns.

      The goal of this study is to comprehensively investigate the soft behavior of two typical categories of Rutherford scattering, that is, processes in which a low-energy massless or a slowly-moving light projectile hits a heavy, static, composite spinning target. For simplicity, we assume the projectile to be a structureless point particle, a spin-$ {1\over 2} $ or spin-$ 0 $ electron. For concreteness, we choose the spin of the composite target particle to range from $ 0 $ to $ 2 $. We find in both cases that the differential cross sections of the Rutherford scattering processes possess the universal form in the first two terms upon heavy target mass expansion yet differ at the next-to-next-to-leading order (NNLO) (depending on the target spin). We conjecture that this pattern may persist for a heavy target particle with arbitrary spin.

      The rest of the paper is structured as follows: In Section II, we present the expression of the tree-level Rutherford scattering amplitude involving a heavy composite spinning target particle and specify the parameterization of the electromagnetic FFs of various target particles. Section III is the main body of the paper, where we present the soft behavior of two types of unpolarized Rutherford scattering cross sections up to NNLO in heavy target mass expansion, assuming the projectile to be a point spin-${1/2}$ electron. For the first type of Rutherford scattering, the first two terms in the differential cross section upon heavy target mass expansion are universal, and the NNLO term starts to exhibit spin-dependence, but with an interesting pattern. For the second type of Rutherford scattering, we observe that the differential cross section is universal at the lowest order in the velocity of the incident slowly-moving projectile, but to all orders in $ 1/M $expansion. At the next-to-leading order (NLO) in projectile velocity, we demonstrate that the first two terms in heavy target mass expansion remain universal, and the NNLO term ($ {\cal O}(v^2/M^2) $ term) exhibits some interesting patterns of target spin dependence. In Section IV, for both types of Rutherford scattering, we attempt to apply heavy particle effective theory (HPET) as well as nonrelativistic QED (NRQED) to reproduce the observed soft behaviors, assuming the target to be a composite spin-${1/2}$ fermion. We summarize the study in Section V. In the appendix, we demonstrate that the observed soft pattern of Rutherford scattering still holds once the projectile is replaced by a point-like spinless electron.

    II.   AMPLITUDE OF RUTHERFORD SCATTERING INVOLVING A HEAVY COMPOSITE TARGET PARTICLE
    • To be specific, let us consider the Rutherford scattering process $ e(k) N(p)\to e(k^\prime) N(p^\prime) $, where N represents a heavy target particle. At the tree-level, Rutherford scattering is induced by a single photon t-channel exchange, as depicted in Fig. 1. The scattering amplitude can be written as

      Figure 1.  Tree-level Feynman diagram of the Rutherford scattering process $ e N\to e N $. The double line represents the heavy target particle, and the black dot denotes the electromagnetic vertex given in (2).

      $ \mathcal{M}=\frac{e^2 g_{\mu\nu}}{q^2}\langle e^-\left(k'\right) \vert J^\mu\vert e^-\left(k\right)\rangle\langle N\left(p',\lambda'\right)|J^\nu|N\left(p,\lambda\right)\rangle, $

      (1)

      where $ J^\mu $ denotes the electromagnetic current, $ q=k-k' $ represents the momentum exchange due to the virtual photon, and $ \lambda,\lambda' $ denote the polarization indices for the massive spinning target particle. For simplicity, we suppress the spin index of the electron.

      The electromagnetic transition matrix element involving the nucleus in (1) is generally a nonperturbative object because the heavy target N is assumed to be any massive composite particle. However, this matrix element can be generally decomposed into a linear combination of independent electromagnetic FFs according to Lorentz group representation [9]:

      $ \begin{aligned}[b]\\[-5pt] \langle N\left(p',\lambda'\right)|J^\mu|N\left(p,\lambda\right)\rangle_{s=0}= 2P^\mu F_{1,0}\left( \frac {q^2}{M^2}\right),\end{aligned} $

      (2a)

      $ \begin{aligned}[b] \langle N\left(p',\lambda'\right)|J^\mu|N\left(p,\lambda\right)\rangle_{s=\frac{1}{2}}= \bar{u}(p',\lambda')\left[2P^\mu F_{1,0}\left( \frac {q^2}{M^2}\right)+ {\rm i} \sigma^{\mu\nu}q_\nu F_{2,0}\left( \frac {q^2}{M^2}\right)\right]u(p,\lambda), \end{aligned} $

      (2b)

      $ \begin{aligned}[b] \langle N\left(p',\lambda'\right)|J^\mu|N\left(p,\lambda\right)\rangle_{s=1}=& -\varepsilon^*_{\alpha'}(p',\lambda') \bigg\{2P^ \mu\left[g^{\alpha'\alpha}F_{1,0}\left( \frac {q^2}{M^2}\right)- \frac {q^{\alpha'}q^\alpha}{{2M^2}}F_{1,1}\left( \frac {q^2}{M^2}\right)\right]\\ & -\left(g^{\mu\alpha'}q^{\alpha}-g^{\mu\alpha}q^{\alpha'}\right)F_{2,0}\left( \frac {q^2}{M^2}\right)\bigg\}\varepsilon_{\alpha}(p,\lambda), \end{aligned} $

      (2c)

      $ \begin{aligned}[b] \langle N\left(p',\lambda'\right)|J^\mu|N\left(p,\lambda\right)\rangle_{s=\frac{3}{2}}=& -\bar{u}_{\alpha'}(p',\lambda') \bigg\{2P^ \mu\left[g^{\alpha'\alpha}F_{1,0}\left( \frac {q^2}{M^2}\right)- \frac {q^{\alpha'}q^\alpha}{{2M^2}}F_{1,1}\left( \frac {q^2}{M^2}\right)\right] \\ & + {\rm i}\sigma^{\mu\nu}q_\nu\left[g^{\alpha'\alpha}F_{2,0}\left( \frac {q^2}{M^2}\right)- \frac {q^{\alpha'}q^\alpha}{{2M^2}}F_{2,1}\left( \frac {q^2}{M^2}\right)\right]\bigg\}u_{\alpha}(p,\lambda), \end{aligned} $

      (2d)

      $ \begin{aligned}[b] \langle N\left(p',\lambda'\right)|J^\mu|N\left(p,\lambda\right)\rangle_{s=2}=& \varepsilon^*_{\alpha'_1\alpha'_2}(p',\lambda')\bigg\{2P^\mu\bigg[g^{\alpha'_1\alpha_1}g^{\alpha'_2\alpha_2}F_{1,0}\left( \frac {q^2}{M^2}\right)- \frac {q^{\alpha'_1}q^{\alpha_1}}{2M^2}g^{\alpha'_2\alpha_2}F_{1,1}\left( \frac {q^2}{M^2}\right) \\ & + \frac {q^{\alpha'_1}q^{\alpha_1}}{2M^2} \frac {q^{\alpha'_2}q^{\alpha_2}}{2M^2}F_{1,2}\left( \frac {q^2}{M^2}\right)\bigg]-\left(g^{\mu\alpha'_2}q^{\alpha_2}-g^{\mu\alpha_2}q^{\alpha'_2}\right) \\ & \times\bigg[g^{\alpha'_1\alpha_1}F_{2,0}\left( \frac {q^2}{M^2}\right)- \frac {q^{\alpha'_1}q^{\alpha_1}}{2M^2}F_{2,1} \left( \frac {q^2}{M^2}\right)\bigg]\bigg\}\varepsilon_{\alpha_1\alpha_2}(p,\lambda). \end{aligned} $

      (2e)

      The various electromagnetic FFs are normalized to be dimensionless. $ P=(p+p')/2 $ is the average momentum of the target particle, and M is the mass of the target particle. u, $ \varepsilon^\mu $, $ u^\mu $, and $ \varepsilon^{\alpha\beta} $ denote the wave function for the spin-${1/ 2}$, $ 1 $, ${3/2}$, and $ 2 $ particles, respectively. Only keeping the Lorentz structures that obey the current conservation, we find that the number of independent electromagnetic FFs is $ 2s+1 $ for the target particle with spin s. Note that the decomposition of the electromagnetic transition matrix element involving a charged particle carrying various spins has been widely studied [11].

      The electromagnetic FFs in (2) encode the internal structure of the composite target particle. In principle, they can be extracted from experiments or computed using nonperturbative theoretical tools. Although the concrete profiles of various FFs depend on the specific target particle, their values near the zero momentum transfer characterize the electromagnetic multipole moments of the composite target particle. For example, $ F_{1,0}(0)=Z $ denotes the total electric charge of the target particle in units of e. $ F_{1,0}(0)+F_{1,1}(0) $, $ F_{2,0}(0) $, and $ F_{2,0}(0)+F_{2,1}(0) $ are the electric quadrupole moment, magnetic dipole moment, and magnetic octupole moment of the composite target particle, in units of $ \dfrac{e}{2M} $, $ \dfrac{e}{ M^2} $, and $ \dfrac{e}{2M^3} $, respectively. The charge radius of a proton can also be expressed as $\langle r_p^2\rangle= \dfrac{3}{2 M^2} [-{F_{1,0}}(0)+{4F'_{1,0}}(0)+F_{2,0}(0)]$ 1. It is worth noting that the commonly used electromagnetic FFs of the proton are the Sachs FFs[12, 13], which are defined as

      $ G_E\left( \frac {q^2}{M^2}\right)=\left(1- \frac {q^2}{4M^2}\right)F_{1,0}\left( \frac {q^2}{M^2}\right)+ \frac {q^2}{4M^2}F_{2,0}\left( \frac {q^2}{M^2}\right), $

      (3a)

      $ G_M\left( \frac {q^2}{M^2}\right)=F_{2,0}\left( \frac {q^2}{M^2}\right). $

      (3b)

      The proton charge radius is defined by $\langle r_p^2\rangle= \dfrac{6}{M^2} G^\prime_E(0)$.

    III.   LOW-ENERGY RUTHERFORD SCATTERING IN HEAVY TARGET MASS EXPANSION
    • Squaring amplitude (1), averaging over spins in the initial state, and summing over the polarizations in the final states, we can straightforwardly obtain the unpolarized differential cross sections of Rutherford scattering in the laboratory frame for various target particle species. In deriving the unpolarized cross sections, the following spin sum relations are useful:

      $ \sum\limits_{\lambda}u(p,\lambda)\bar{u}(p,\lambda)= \frac {\not p+M}{2M}, $

      (4a)

      $ \sum\limits_{\lambda}\varepsilon_\alpha(p,\lambda)\varepsilon^*_{\alpha'}(p,\lambda)=\eta_{\alpha\alpha'}, $

      (4b)

      $ \sum\limits_{\lambda}u_\alpha(p,\lambda)\bar{u}_{\alpha'}(p,\lambda)$

      $ \begin{aligned}[b] =- \frac {\not p+M}{2M}\Bigg({g_{\alpha\alpha'}}{-\frac{1}{3}\gamma_\alpha\gamma_{\alpha'}-\frac{2p_\alpha p_{\alpha'}}{3M^2}+\frac{\gamma_{\alpha'}p_\alpha-\gamma_{\alpha'}p_{\alpha}}{3M}}\Bigg), \end{aligned} $

      (4c)

      $ \begin{aligned}[b]& \sum\limits_{\lambda}\varepsilon_{\alpha_1\alpha_2}(p,\lambda)\varepsilon^*_{\alpha'_1\alpha'_2}(p,\lambda)\\=&\eta_{\alpha_1\alpha'_1}\eta_{\alpha_2\alpha'_2}+\eta_{\alpha_1\alpha'_2}\eta_{\alpha_2\alpha'_1}- \frac {2}{3}\eta_{\alpha_1\alpha_2}\eta_{\alpha'_1\alpha'_2}, \end{aligned} $

      (4d)

      with $\eta_{\alpha\beta} \equiv -g_{\alpha\beta}+ \dfrac {p_\alpha p_\beta}{M^2}$. Note the Dirac spinor wave function is normalized as $ \bar{u}(p,r) u(p,s)=\delta^{rs} $.

    • A.   Massless spin-1/2 projectile

    • We first consider the $ ep $ elastic scattering experiment, with the famous Hofstadter's experiment as the prototype [3]. In such cases, the incident electron is treated as massless, and we are concerned with the low-energy limit $ |{\bf k}| \ll M $.

      We focus on Rutherford scattering in the laboratory frame, with the four-momentum of the target particle in the initial state signified by $ p^\mu=(M,{\bf 0}) $. The corresponding differential unpolarized cross section is defined by

      $ \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}= \frac {1}{2 |{\bf k}|}\cdot \frac {1}{2M}\cdot \frac {{\bf k'}^2}{8\pi|{\bf k}| M}\left( \frac {1}{2} \frac {1}{2s+1}\sum\limits_\text{spins}\left|\mathcal{M}\right|^2\right), $

      (5)

      where θ denotes the polar angle between the incident and reflected electron, and $ |{{\bf k}^\prime}| $ is a function of $ |{\bf k}| $, $ \cos\theta $, and M:

      $ |{\bf k}'|= {|{\bf k}| \over 1+ \dfrac{|{\bf k}|}{M}\left ( 1-\cos\theta\right)}. $

      (6)

      The full expressions of the unpolarized cross sections are generally lengthy and appear cumbersome, and therefore it is difficult to recognize any clear pattern regarding the dependence on the heavy target particle spin. Hopefully, once heavy target mass expansion is conducted, the soft behavior of Rutherford scattering will become transparent and we may readily identify a simple pattern.

      After expanding both the squared amplitude and phase space measure (the factor $ {{\bf k}'}^2/{\bf k}^2 $) in the (5) powers of $ 1/M $, the differential Rutherford scattering cross sections become considerably simpler. We find that the first two orders in heavy target expansion are universal, for example, independent of the heavy target spin,

      $ \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta} = \dfrac {\pi\alpha^2Z^2\cos^2\dfrac{\theta}{2}}{2 {\bf k}^2\sin^4\left(\dfrac{\theta}{2}\right)}-\dfrac{\pi\alpha^2Z^2\cos^2\dfrac{\theta}{2}}{M |{\bf k}|\sin ^2\left(\dfrac{\theta}{2}\right)}+\mathcal{O}\left(\dfrac{1}{M^2}\right). $

      (7)

      For clarity, we substitute $ F_{1,0}=Z $. The result is intuitively clear: In the soft limit, the long wavelength photon can only feel the total charge of the composite target particle and is insensitive to any further details about its internal structure.

      In contrast, the NNLO terms in heavy target mass expansion vary with different heavy target particles:

      $ \begin{aligned}[b]\\[-8pt] \left( \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}\right)_\text{NNLO}^{s}= & - \frac {4\pi\alpha^2}{M^2\sin^2\dfrac{\theta}{2}} \bigg\{ \Theta\left(s- \frac {1}{2}\right)\frac{s+1}{48s}F^2_{2,0}\left(\cos\theta-3\right) + \frac {1}{4}\cos^2 \frac {\theta}{2} \bigg[4F'_{1,0}Z -\Theta\left(s-1\right) \frac {2}{3}F_{1,1}Z \\ & +\Theta\left(s- \frac {1}{2}\right) \frac {5-(-1)^{2s}}{6}F_{2,0}Z+Z^2\cos\theta - \frac {1}{3}\left(\lceil s\rceil +s+3\right)Z^2 \bigg]\bigg\}, \qquad \left(s=0, \frac {1}{2},1, \frac {3}{2},2 \right) \end{aligned} $

      (8)

      where $ \Theta(s) $ is the Heaviside step function with $ \Theta(0)=1 $, and $ \lceil s\rceil $ denotes the ceiling function mapping s to the least integer greater than or equal to s. For notational brevity, we neglect the argument $ 0 $ in various FFs. We observe that the $ F'_{1,0}Z $, $ Z^2\cos\theta $, and $ F_{1,1}Z $ terms with a prefactor $ \cos^2(\theta/2) $ are still universal, that is, independent of the target spin. In fact, the $ F'_{1,0}Z $ and $ Z^2\cos\theta $ terms have the same origin of the leading order (LO) and NLO cross sections, which correspond to different terms in the Taylor expansion of $ F^2_{1,0}(q^2/M^2) $ in the squared LO amplitude and phase space measure. The coefficient of the $ F_{2,0}Z $ term seems to reflect the spin-statistic characteristic of the target particle. For fermions, the coefficient is $ 1 $, whereas for bosons, it is $ 2/3 $.

      Although the coefficients of $ F^2_{2,0}(\cos\theta -3) $ inside the bracket depend on the target particle spin s, they seem to fit into the expression $\dfrac{1+s}{ 48s}$ (for $ s=1/2,1,3/2,2 $). It will be interesting to see whether this pattern still persists for higher target spin.

      Owing to tremendous phenomenological interest, it is interesting to specialize in the spin-${1/ 2}$ target. Upon heavy target mass expansion, the differential $ ep $ elastic scattering cross section can be expressed in terms of the Sachs FFs:

      $ \begin{aligned}[b] \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta} = \dfrac {\pi\alpha^2\cos^2\dfrac{\theta}{2}}{2 {\bf k}^2\sin^4\left(\dfrac{\theta}{2}\right)}-\dfrac{\pi\alpha^2\cos^2\dfrac{\theta}{2}}{M |{\bf k}|\sin ^2\left(\dfrac{\theta}{2}\right)}- \dfrac {4\pi\alpha^2}{M^2\sin^2\dfrac{\theta}{2}} \end{aligned} $

      $ \begin{aligned}[b] \times\bigg[\dfrac{1}{16}G_M^2\left(\cos\theta-3\right)+ \dfrac {1}{4} \cos^2 \dfrac {\theta}{2}\left( 4G'_E+\cos\theta- \dfrac {1}{2}\right)\bigg], \end{aligned} $

      (9)

      or equivalently, expressed in terms of proton charge radius and magnetic dipole moment $ \mu_p $ by making the replacements $ G_M\to \mu_p $ and $G^\prime_E\to \dfrac{M^2}{ 6} \langle r^2_p \rangle$. The contributions of $ \langle r^2_p\rangle $ and $ \mu_p $ begin at NNLO in heavy target mass expansion, whose values may be extracted by carefully studying the angular distributions.

    • B.   Light non-relativistic spin-1/2 projectile

    • Next, we turn to the soft limit of the original prototype of the Rutherford scattering process, that is, a slowly moving light particle hits a heavy static target. We again assume that the projectile is a Dirac fermion, whose mass and momentum are denoted by m and $ \bf k $.

      The differential cross section for this type of Rutherford scattering in the laboratory frame is defined by

      $ \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}= \frac {1}{32\pi M}\left[{p'^0+k'^0\left(1- \frac {|\mathbf{k}|}{|\mathbf{k}'|}\cos\theta\right)}\right]^{-1} \frac {|\mathbf{k}'|} {|\mathbf{k}|}\left|\mathcal{M}\right|^2. $

      (10)

      The resulting expressions are rather lengthy. Fortunately, we are only interested in its soft behavior. Since there are three widely separated scales in this process, which obey $ |{\bf k}|\ll m \ll M $, the appropriate way of extracting the soft behavior is to expand the differential cross sections in powers of $ v=|{\bf k}|/m $ (velocity of the projectile) and $ 1/M $ simultaneously. The necessity of performing double expansion renders this case somewhat more complicated than the preceding case discussed in Section IIIA.

      Interestingly, at the lowest order in velocity yet to all orders in $ 1/M $, the differential cross sections scale as $ 1/{\bf k}^4 $, which takes a uniform form:

      $ \begin{aligned}[b] \left( \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}\right)^{s}_{(v^0)} =& {2 \pi Z^2 \alpha^2 \over {\bf k}^4} { m^2 (M+m)^2 \left(\sqrt{M^2-m^2\sin^2\theta}+m \cos\theta\right)^2 \over M \sqrt{M^2-m^2\sin^2\theta}\left(M-\cos\theta \sqrt{M^2-m^2\sin^2\theta}+m \sin^2\theta\right)^2} \\ =& \frac{8 \pi Z^2 \alpha ^2 m^2}{\mathbf{k}^4 \sin^4\dfrac{\theta}{2}}-\frac{\pi Z^2 \alpha ^2 m^4}{M^2\mathbf{k}^4}+ \mathcal{O}\left({m^6\over M^4 {\bf k}^4}\right). \end{aligned} $

      (11)

      The LO term is the extremely well-known expression obtained by Rutherford using classical mechanics, dubbed the Rutherford formula.

      At NLO in velocity expansion, the differential cross sections scale as $ 1/{\bf k}^2 $, whose explicit expressions are still rather complicated and vary with different target species. Nevertheless, once heavy target mass expansion is conducted, a clear pattern emerges:

      $ \begin{aligned}[b] \left( \dfrac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}\right)^{s}_{(v^2)} = & \dfrac {\pi\alpha^2}{\mathbf{k}^2\sin^2\dfrac{\theta}{2}} \Bigg[\dfrac{Z^2 \cos^2\dfrac{\theta}{2}}{2\sin^2\dfrac{\theta}{2}} -\dfrac{Z^2 m\cos^2\dfrac{\theta}{2}}{M} \\&-\dfrac{Z m^2}{4M^2 } f^s_{\rm NNLO} + \mathcal{O}\left( \dfrac {1}{M^3}\right)\Bigg], \end{aligned} $

      (12)

      where

      $ \begin{aligned}[b] f_{\rm NNLO}^{s} =& 16F'_{1,0}+Z\cos\theta -\Theta\left(s-1\right) \frac {8}{3}F_{1,1}\\&+\Theta\left(s- \frac {1}{2}\right) \frac {2\left[5-(-1)^{2s}\right]}{3}F_{2,0} \\ & - \frac {4}{3}\left(\lceil s\rceil +s+ \frac {3}{4}\right)Z. \qquad \left( s=0, \frac {1}{2},1, \frac {3}{2},2\right) \end{aligned} $

      (13)

      The LO and NLO terms in$ 1/M $ expansion are universal. The NNLO terms begin to exhibit target spin dependence. However, even at $ {\cal O}(v^2/M^2) $, the $ F^\prime_{1,0} $, $ F_{1,1} $, and $ Z\cos\theta $ terms still seem to be universal, that is, independent of the target particle spin. The coefficient of $ F_{2,0} $ seems to reflect the spin-statistic characteristic of the target particle. For fermions, the coefficient is $ 4 $, whereas for bosons, it is ${8/ 3}$.

    IV.   REPRODUCING THE SOFT BEHAVIOR FROM EFFECTIVE FIELD THEORY
    • The low-energy limit of Rutherford scattering is largely dictated by a heavy target particle interacting with a soft photon. Therefore, it is natural to expect that the soft behavior can be reproduced by an effective field theory (EFT) analogous to heavy quark effective theory (HQET), which automatically incorporates heavy target mass expansion. In this section, we specialize in the case of a spin-$ 1/2 $ composite target particle.

      Originally, HQET was designed to describe a structureless heavy quark interacting with soft gluons [13]. Owing to the asymptotic freedom property of QCD, Wilson coefficients can be computed in perturbation theory through a perturbative matching procedure.

      The key idea of HQET can be readily transplanted to the case of a heavy composite particle interacting with a soft photon, as long as the photon wavelength is too long to deeply probe the internal structure of the composite target. Consequently, one is generally unable to calculate various Wilson coefficients from the top-down perspective. The internal structure of the composite heavy target particle is encoded in various Wilson coefficients, which essentially represent various multipole moments. They can, in principle, be evaluated by nonperturbative means or determined using the bottom-up approach, for example, extracted from low-energy Rutherford scattering experiments.

      Analogous to HQET, we build up an EFT, dubbed HPET, describing a static heavy composite fermionic target particle interacting with a soft photon:

      $ \mathcal{L}_{\rm HPET}=\bar{h}_v\Bigg(i D_0+c_2 \frac {\mathbf{D}^2}{2M}+c_Fe \frac {\boldsymbol{\sigma}\cdot\mathbf{B}}{2M}+c_De \frac {\left[\boldsymbol{\nabla}\cdot\mathbf{E}\right]}{8M^2}+ ic_Se \frac {\boldsymbol{\sigma}\cdot\left(\mathbf{D}\times\mathbf{E}-\mathbf{E}\times\mathbf{D}\right)}{8M^2}\Bigg) h_v + {\mathcal O}(1/M^3),$

      (14)

      where we truncate the effective Lagrangian through order $ 1/M^2 $. $ h_v $ represents the heavy target HPET field, with the label velocity $ v^\mu=(1,{\bf 0}) $, $D^\mu=\partial^\mu+ {\rm i} ZeA^\mu$ signifies the covariant derivative, $ \mathbf{E} $ and $ \mathbf{B} $ denote the electric and magnetic field, the coefficient $ c_2=1 $ is a rigorous consequence of Lorentz symmetry, and the $ c_F $, $ c_D $, and $ c_S $-related terms are often referred to as Fermi, Darwin, and spin-orbital terms. The organization of the HPET Lagrangian is governed by powers of $ |{\bf q}|/M $, with $ {\bf q} $ signifying the photon momentum.

    • A.   HPET description of a massless spin-1/2 projectile hitting a static spin-${1/ 2}$ target

    • In contrast with Fig. 1, up to order $ 1/M^2 $, there are five Feynman diagrams in the context of HPET for the tree-level process $ e(k) N(p)\to e(k^\prime) N(p^\prime) $. The corresponding amplitude reads

      $ \begin{aligned}[b] \mathcal{M}_{\rm HPET} = & -\sqrt{1+c_2 \frac {\mathbf{p}'^2}{2M^2}} \frac {e^2}{q^2} \bigg\{-Z\bar{u}_\text{NR}u_\text{NR}\bar{u}(k')\gamma^0u(k) + \frac {c_2Z}{2M}\bar{u}_\text{NR}u_\text{NR}\bar{u}(k')\mathbf{p}'\cdot\boldsymbol{\gamma} u(k) \\ & - \frac {c_F}{4M}\bar{u}_\text{NR}\left[\not q,\gamma^\mu\right]u_\text{NR}\bar{u}(k')\gamma_\mu u(k) - \frac {c_Dq^2}{8M^2}\bar{u}_\text{NR}u_\text{NR}\bar{u}(k')\gamma^0 u(k)\bigg\} \\ = & \frac {e^2}{q^2} \bigg\{ Z\bar{u}_\text{NR}u_\text{NR}\bar{u}(k')\gamma^0u(k) + \frac {c_F}{4M}\bar{u}_\text{NR}\left[\not q,\gamma^\mu\right]u_\text{NR}\bar{u}(k')\gamma_\mu u(k) + \frac {c_Dq^2}{8M^2}\bar{u}_\text{NR}u_\text{NR}\bar{u}(k')\gamma^0 u(k) \bigg\}, \end{aligned} $

      (15)

      where $ u_\text{NR} $ denotes the HPET spinor wave function that satisfies $\not v u_\text{NR}(v,s)=u_\text{NR}(v,s)$ and $ \bar{u}_\text{NR}(v,s)u_\text{NR}(v,s)=2v^0 $. Notice that the final expression does not depend on $ c_2 $. Moreover, note that the contribution from the $ c_S $ term is proportional to $ q^0\sim \mathbf{k}^2/M $; hence, it does not need to be considered at the prescribed accuracy of $ 1/M^2 $.

      Squaring the amplitude in (15) and summing/averaging over polarizations, we find the differential unpolarized cross section obtained from the HPET side:

      $ \begin{aligned}[b] \dfrac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}\bigg|_{\rm EFT} = & \dfrac {\pi\alpha^2Z^2\cos^2\dfrac{\theta}{2}} {2 {\bf k}^2 \sin^4\dfrac{\theta}{2}}-\dfrac{\pi\alpha^2Z^2\cos^2\dfrac{\theta}{2}}{M |{\bf k}| \sin^2\dfrac{\theta}{2}} \\& - \dfrac {\pi\alpha^2}{8M^2\sin^2\dfrac{\theta}{2}}\Big[Z^2\left(\cos 2\theta-1\right) \\&+c_D Z\left(\cos\theta+1\right)+c_F^2\left(\cos\theta-3\right)\Big]. \end{aligned} $

      (16)

      Note that the LO and NLO terms are only sensitive to the total charge $ Ze $ of the target particle and do not depend on any nontrivial Wilson coefficients (hence they are insensitive to the target's internal structure). This indicates that these two terms are solely dictated by the leading HPET Lagrangian. Since the leading HPET Lagrangian possesses heavy particle spin symmetry, the spin degree of freedom is completely decoupled at the lowest order in $ 1/M $ expansion. From a technical perspective, the first two terms arise solely from expanding the squared LO amplitude in (15) as well as expanding the factor $ |{\bf k}^\prime|/|{\bf k}| $ in the phase space measure in Eq. (10).

      At first sight, one may worry that the interference between the $ c_F $ term and the LO amplitude would nominally generate a $ {\cal O}(1/M) $ correction, thus breaking the universality at NLO. A closer examination reveals that this contribution actually vanishes after summing over polarizations. This cancellation is anticipated to persist for other species of spinning target particles.

      The first two terms in (16) are indeed identical to the universal behaviors revealed in (7). The EFT approach helps us to better understand why they are independent of the species of the composite target particle, though our HPET Lagrangian only specializes in the spin-$ 1/2 $ target.

      At NNLO in $ 1/M $ expansion, three terms emerge, which stem from different sources. The first term clearly stems from expanding the squared LO amplitude in combination with phase space expansion, which is of the same origin as the LO and NLO contributions and is anticipated to be universal. The second term comes from the interference between the $ c_D $ term and the LO amplitude, and the last term stems from the square of the $ c_F $ term. Concretely, the last two terms depend on the composite target particle's charge radius and magnetic dipole. Interestingly, the $ c_F $ term can be identified with the $ F^2_{2,0}(\cos\theta-3) $ term in (8). As discussed in the paragraph after (8), the coefficient of this term may depend on the target spin in a specific manner.

      To verify that the EFT amplitude does capture the correct soft behavior, we can perform heavy target mass expansion from the full QED amplitude in Eq. (1):

      $ \begin{aligned}[b] {\cal M}_{\rm QED} = & \frac{e^2}{q^2}\bar{u}(k^\prime)\gamma^\mu u(k) \bar{u}(p',\lambda')\left[2P^\mu F_{1,0}\left( \frac {q^2}{M^2}\right)+i\sigma^{\mu\nu}q_\nu F_{2,0}\left( \frac {q^2}{M^2}\right)\right]u(p,\lambda) \\ =& \frac{e^2}{q^2}\bar{u}\gamma^\mu u\bar{u}_\text{NR}^{\lambda'}\sqrt{ \frac {p'^0}{M}}\left(1- \frac {\mathbf{p}'\cdot\boldsymbol{\gamma}}{2M}- \frac {\mathbf{p}'^2}{8M^2}\right)\left[2P^\mu \left(F_{1,0}+F'_{1,0} \frac {q^2}{M^2}\right)+i\sigma^{\mu\nu}q_\nu F_{2,0}\right]u_\text{NR}^\lambda \\ =& \frac{2Me^2}{q^2}\left[ Z\bar{u}\gamma^0 u\bar{u}_\text{NR}^{\lambda'} u_\text{NR}^\lambda + \frac {F_{2,0}}{4M}\bar{u}\gamma_\mu u\bar{u}_\text{NR}^{\lambda'}\left[\not q, \gamma^\mu\right] u_\text{NR}^\lambda + \frac {q^2\left(8F'_{1,0}-F_{1,0} + 2F_{2,0}\right)}{8M^2}\bar{u}\gamma^0 u\bar{u}_\text{NR}^{\lambda'}u_\text{NR}^\lambda\right], \end{aligned} $

      (17)

      where we not only expand the FF $ F_{1,0} $ to the first order in $ q^2/M^2 $, but also expand the Dirac spinor using $u(p')=\sqrt{\dfrac{p^{\prime 0}}{M}}\left(1- \dfrac {\mathbf{p}'\cdot\boldsymbol{\gamma}}{2M}- \dfrac {\mathbf{p}'^2}{8M^2}\right) u_\text{NR}+\mathcal{O}(1/M^3)$.

      Note that the HPET amplitude assumes nonrelativistic normalization for the target particle; therefore, we must include an overall factor $ 2M $ prior to comparing it with the full QED amplitude. By equating (15) and (17), we are able to identify the relations between the Wilson coefficients in HPET and the electromagnetic FFs near the zero-momentum transfer:

      $ c_F= F_{2,0}, $

      (18a)

      $ c_D= 2F_{2,0}+8F'_{1,0}-F_{1,0}, $

      (18b)

      which are identical to the relations obtained for the structureless quark in HQET .

      Substituting relation (18) into (16), we fully reproduce the NNLO contribution for a heavy spin-$ 1/2 $ target, as recorded in (8).

    • B.   NRQED+HPET description of a slowly-moving spin-1/2 projectile hitting a static spin-${1/ 2}$ target

    • Next, we turn to the EFT approach to understand the second type of Rutherford scattering, where a light non-relativistic particle hits a static heavy composite target. To be specific, we specialize in a spin-$ 1/2 $ structureless projectile and a spin-$ 1/2 $ target particle. The treatment of the static composite fermionic target is identical to that in Section IV.A. It is natural to apply NRQED to describe the incident slowly-moving electron.

      Up to the relative order $ v^2 $, the electron sector of the NRQED Lagrangian reads as

      $ \begin{aligned}[b] \mathcal{L}_\text{NRQED}=&\psi^\dagger\bigg[{\rm i} D^0 + d_2 \frac {\mathbf{D}^2}{2 m}+d_4 \frac {\mathbf{D}^4}{8 m^3}+d_Fe \frac {\boldsymbol{\sigma}\cdot\mathbf{B}}{2 m}\\&+d_De \frac {\left[\boldsymbol{\nabla}\cdot\mathbf{E}\right]}{8 m^2} + {\rm i}d_Se \frac {\boldsymbol{\sigma}\cdot\left(\mathbf{D}\times\mathbf{E}-\mathbf{E}\times\mathbf{D}\right)}{8 m^2}\bigg]\psi, \end{aligned} $

      (19)

      where ψ denotes a Pauli spinor field that annihilates a nonrelativistic electron, $ d_2=d_4=1 $ is a rigorous consequence of Lorentz invariance, and the $ d_4 $ term, together with the $ d_F $, $ d_D $, and $ d_S $ terms (referred to as the Fermi, Darwin, and spin-orbital terms), represent the $ O(v^2) $ corrections to the NRQED Lagrangian. At the tree level, the Wilson coefficients $ d_F = d_D= d_S = 1 $.

      Our starting point is the HPET Lagrangian (14) and the NRQED Lagrangian (19). It is convenient to work in Coulomb gauge. Up to $ {\cal O}(v^2/M^2) $, the relevant tree-level EFT amplitude for $ e N\to e N $ reads as

      $ \begin{aligned}[b]\\[-8pt] \mathcal{M}_{\rm EFT}= & \frac {e^2}{\mathbf{q}^2}\xi^\dagger \left[1+ \frac {d_2}{4m^2}(\mathbf{k}^2+\mathbf{k}'^2)- \frac {d_D}{8m^2} |\mathbf{k}'-\mathbf{k}|^2- \frac {{\rm i} d_S}{4m^2}\boldsymbol{\sigma}\cdot(\mathbf{k}\times \mathbf{k}')\right] \xi \, \bar{u}_\text{NR}^{\lambda'} \Big[-Z+ \frac {\left(c_D-2c_2Z\right)\mathbf{p}'^2}{8M^2} \Big]u_\text{NR}^\lambda \\& - \frac {1}{q^2}\left(\delta^{ij}- \frac {q^iq^j}{\mathbf{q}^2}\right)\xi^\dagger \bigg\{ (k^i+k^{'i}) \left[ \frac {d_2}{2m}+ \frac {d_2^2-d_4}{8m^3}(\mathbf{k}^2+\mathbf{k}'^2)\right] + \frac {{\rm i} d_F}{2m}\left(1+ d_2 {\mathbf{k}^2+\mathbf{k}'^2 \over 4 m^2}\right)\left[\boldsymbol{\sigma}\times(\mathbf{k}'-\mathbf{k})\right]^i \end{aligned} $

      (20)

      $ \begin{aligned}[b] - \frac {d_D}{16m^3}(k'^i-{k}^i)(\mathbf{k}'^2-\mathbf{k}^2) - \frac {{\rm i} d_S}{16m^3}(\mathbf{k}'^2-\mathbf{k}^2)\left[\boldsymbol{\sigma}\times(\mathbf{k}+\mathbf{k}')\right]^i\bigg\}\xi \, \bar{u}_\text{NR}^{\lambda'} \Big[- \frac {c_2Z}{2M}{p}'^j- {\rm i} \frac {c_F}{2M}\sigma^{jl}{p}'^l \Big]u_\text{NR}^\lambda, \end{aligned} $

      (20)

      where the first line represents temporal photon exchange, and the remaining lines represent transverse photon exchange. ξ denotes the two-component spinor wave function.

      After some simplification, (20) reduces to

      $ \begin{aligned}[b] \mathcal{M}_{\rm EFT} = & \frac {e^2}{\mathbf{q}^2}\Bigg[-Z+ \frac {\left(c_D-2c_2Z\right)\mathbf{p}'^2}{8M^2} \Bigg] \xi^\dagger\left[1+ \frac {d_2}{4m^2}\left(\mathbf{k}^2+\mathbf{k}'^2\right)- \frac {d_D}{8m^2}|\mathbf{k}'-\mathbf{k}|^2- \frac {{\rm i} d_S}{4m^2}\boldsymbol{\sigma}\cdot(\mathbf{k}\times \mathbf{k}')\right]\xi\bar{u}_\text{NR}^{\lambda'} u_\text{NR}^\lambda - \frac {c_Fe^2}{4M\mathbf{q}^2}\xi^\dagger\\& \times\bigg\{ (k^i+{k}'^i) \left[ \frac {d_2}{2m}+ \frac {d_2^2-d_4}{8m^3}(\mathbf{k}^2+\mathbf{k}'^2)\right] + \frac {{\rm i} d_F}{2m}\left[\boldsymbol{\sigma}\times(\mathbf{k}'-\mathbf{k})\right]^i- \frac {{\rm i} d_S}{16m^3}(\mathbf{k}'^2-\mathbf{k}^2)\left[\boldsymbol{\sigma}\times(\mathbf{k}'-\mathbf{k})\right]^i\bigg\}\xi \times \bar{u}_\text{NR}^{\lambda'} \Big[\gamma^i,\boldsymbol{\gamma}\cdot \mathbf{q} \Big]u_\text{NR}^\lambda. \end{aligned} $

      (21)

      Squaring the amplitude in (21), and summing/averaging over various spins, we obtain the differential unpolarized Rutherford scattering cross section in the context of EFT:

      $ \begin{aligned}[b] \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta} \bigg|_{\rm EFT} =& \: \frac{\pi \alpha ^2 m^2Z^2}{2\mathbf{k}^4 \sin^4\frac{\theta}{2}} -\frac{\pi \alpha ^2 m^4Z^2}{M^2\mathbf{k}^4} + \frac {\pi\alpha^2Z}{2\mathbf{k}^{2}\sin^2\frac{\theta}{2}} \bigg\{ \frac {Z\left(d_D\cos\theta-d_D+2\right)}{2\sin^2\frac{\theta}{2}} \\& - {m\over M} (d_D\cos\theta-d_D+2) -{m^2\over 2 M^2} \left[ Z(2+d_D-4c_2)+Z\cos\theta (2-d_D) + 2c_D\right]\bigg\}. \end{aligned} $

      (22)

      The $ d_S $ term in (21) does not contribute to the squared amplitude because its interference with the LO amplitude in velocity expansion only contains a single Pauli matrix and hence vanishes upon summing over polarization.

      Substituting $ c_2=d_D=1 $ in (22), and utilizing the relations given in (18), we exactly reproduce (11), which encodes the LO and NLO terms in heavy target expansion, and (13), which encapsulates the NNLO term. (13) indicates that the $ Z\cos\theta $ and $ F^\prime_{1,0} $ terms in NNLO correction are universal, for example, independent of the target spin. This may indicate that structures such as $ Z\cos\theta (2-d_D) + 2c_D $ arise ubiquitously in an EFT calculation for heavy targets other than spin-$ 1/2 $ fermions.

      To verify that the EFT amplitude indeed reproduces the correct soft behavior, we conduct both nonrelativistic and heavy target mass expansion from the full QED amplitude in (1).

      Working again in the Coulomb gauge, and employing the following relation between the relativistic electron spinor and nonrelativistic electron spinor:

      $u(k)= \frac {1}{\sqrt{k^0+m}} \begin{pmatrix} \left(k^0+m\right)\xi \nonumber\\ \mathbf{k}\cdot\boldsymbol{\sigma}\xi \end{pmatrix} ,\quad \bar{u}(k)= \frac {1}{\sqrt{k^0+m}} \begin{pmatrix} \left(k^0+m\right)\xi^\dagger -\xi^\dagger\mathbf{k}\cdot\boldsymbol{\sigma} \end{pmatrix},$

      we expand the full QED amplitude through $ {\cal O}(v^2/M^2) $:

      $ \begin{aligned}[b] \mathcal{M}_{\rm QED} = & -\frac{e^2}{\mathbf{q}^2}\bar{u}(k^\prime)\gamma^0 u(k) \bar{u}_\text{NR}^{\lambda'}\left[ 2ZP^0 + \frac {P^0}{4M^2}\left(8F'_{1,0}q^2+Z{\mathbf{q}^2}\right) +i \frac {\mathbf{p}'\cdot\boldsymbol{\gamma}}{2M}\sigma^{0i}q_i F_{2,0}\right]u_\text{NR}^\lambda \\ & -\frac{e^2}{\mathbf{q}^2}\left(\delta_{ij}- \frac {q^iq^j}{\mathbf{q}^2}\right)\bar{u}\gamma^i u\bar{u}_\text{NR}^{\lambda'}\left[ 2ZP^j + \frac {1}{2}\left[\not q, \gamma^j\right] F_{2,0} + \frac {P^j}{4M^2}\left(8F'_{1,0}q^2+Z{\mathbf{q}^2}\right) \right]u_\text{NR}^\lambda \\ =& \frac{2Me^2}{\mathbf{q}^2}2m\bigg[ -Z + \frac {\mathbf{p}'^2}{8M^2}\bigg(2F_{2,0}+8F'_{1,0}-3Z\bigg) \bigg]\xi^\dagger\left[ 1 + \frac {\left|\mathbf{k}+\mathbf{k}'\right|^2}{8m^2} - \frac {i}{4m^2} \boldsymbol{\sigma}\cdot\left(\mathbf{k}\times\mathbf{k}'\right)\right]\xi \bar{u}_\text{NR}^{\lambda'}u_\text{NR}^\lambda \\ & -\frac{F_{2,0}e^2}{2\mathbf{q}^2}2m\xi^\dagger\bigg\{ \frac {1}{2m}\left({k}^i+{k}'^i\right)+ \frac {i}{2m}\left[\boldsymbol{\sigma}\times\left(\mathbf{k}'-\mathbf{k}\right)\right]^i- \frac {i}{16m^3}\left(\mathbf{k}'^2-\mathbf{k}^2\right)\left[\boldsymbol{\sigma}\times\left(\mathbf{k}'-\mathbf{k}\right)\right]^i \bigg\}\xi \bar{u}_\text{NR}^{\lambda'}\left[\gamma^i,\boldsymbol{\gamma}\cdot\mathbf{q} \right]u_\text{NR}^\lambda. \end{aligned} $

      (23)

      After including the normalization factor $ (2M)(2 m) $, employing the relations for the heavy target Wilson coefficients in (18), and taking $ d_2=d_4=d_F=d_D=d_S=1 $, we find that the EFT amplitude (21) exactly agrees with the full QED amplitude (24).

    V.   SUMMARY
    • In this study, we conduct a comprehensive investigation of the soft behavior of the tree-level Rutherford scattering process. We consider two classes of Rutherford scattering experiments: One in which a low-energy point-like massless projectile (for example, a spin-${1/ 2}$ or spin-$ 0 $ electron) bombards a static massive composite spinning target particle (for example, an atomic nucleus), and one where a slowly-moving light structureless projectile hits a static heavy composite spinning target. We consider various composite target particles with spin up to 2.

      The soft limit of the differential Rutherford cross sections in the laboratory frame in both cases exhibits some universal pattern. For the first type of Rutherford scattering process, given a specific projectile, the first two terms in the differential cross section are universal upon heavy target mass expansion, whereas the universality starts to break down at NNLO. Nevertheless, many terms at NNLO still remain spin-independent or have some definite spin-dependence pattern. For the second type, we must perform both nonrelativistic and heavy target mass expansion to infer the correct soft limit. At the lowest order in projectile velocity expansion, but to all orders in $ 1/M $expansion, the differential cross section has a universal form (insensitive to the projectile spin). At NLO in velocity expansion, the first two terms in the differential cross section in$ 1/M $ expansion are still universal. The $ {\cal O}(v^2/M^2) $ piece starts to partially violate the universality. Despite this, some terms at this order still remain target spin independent.

      It is of special interest that the $ F_{2,0} $ term (magnetic dipole of the target particle) at $ {\cal O}(1/M^2) $ (or $ {\cal O}(v^2/M^2) $ for the second type of Rutherford scattering) seems to reflect a peculiar spin-statistics feature. Its coefficient keeps one constant for the fermionic target and another constant for the bosonic target. It is interesting to verify this observation by investigating target particles with even higher spin.

      We also attempt to apply the EFT approach to understand the soft pattern of the Rutherford scattering cross sections, taking the target particle as a composite Dirac fermion for concreteness. Some useful insight is gained from the EFT perspective. However, since we exclusively consider the spin-${1/ 2}$ target particle, the EFT approach itself is unable to account for the specific pattern of target spin dependence at NNLO in heavy target mass expansion observed in this study. Further exploration is needed to understand this issue.

    ACKNOWLEDGMENTS
    • We are grateful for useful discussions with Zhewen Mo.

    APPENDIX A: RUTHERFORD SCATTERING WITH MASSLESS SPINLESS PROJECTILE
    • We can repeat our investigation in Section III.A by replacing the projectile with a massless spin-0 electron, which is described by scalar QED. The electromagnetic vertex involving a scalar electron is simply given by

      $ \langle e (k')\vert J^\mu\vert e(k)\rangle=-\left(k^\mu+k'^\mu\right). $

      (A1)

      Upon heavy target mass expansion, we again observe that the unpolarized cross sections exhibit some universal feature. Concretely, the LO and NLO pieces in$ 1/M $ expansion are independent of the target particle spin:

      $ \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}= \frac {\pi\alpha^2Z^2}{2{\bf k}^2\sin^4\left(\dfrac{\theta}{2}\right)} -\frac{\pi\alpha^2Z^2}{M |{\bf k}|\sin ^2\left(\dfrac{\theta}{2}\right)}+{\cal O}\left(\dfrac{1}{M^2}\right). $

      (A2)

      The universality becomes partially violated at NNLO. For various target particles, the NNLO contributions are

      $ \begin{aligned}[b]\\[-11pt] \left( \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}\right)_\text{NNLO}^{s}= & - \frac {4\pi\alpha^2}{M^2\sin^2\dfrac{\theta}{2}} \bigg[F'_{1,0}Z -\Theta\left(s- \frac {1}{2}\right)\frac{s+1}{48s}F^2_{2,0}\left(\cos\theta+1\right) +\Theta\left(s- \frac {1}{2}\right) \frac {5-(-1)^{2s}}{24}F_{2,0}Z \\ & +\frac{5}{16}Z^2\cos\theta -\Theta\left(s-1\right) \frac {1}{6}F_{1,1}Z - \frac {1}{12}\left(\lceil s\rceil +s+ \frac {15}{4}\right)Z^2 \bigg]. \qquad \left( s=0, \frac {1}{2},1, \frac {3}{2},2\right) \end{aligned} $

      (A3)

      Similar to the pattern indicated in (8) for a massless spin-${1/ 2}$ projectile, we observe that the $ F'_{1,0}Z $, $ Z^2\cos\theta $, and $ F_{1,1}Z $ terms are independent of the target spin. The $ F'_{1,0}Z $ and $ Z^2\cos\theta $ terms actually have the same origin of the LO and NLO cross sections, which correspond to different terms in the Taylor expansion of $ F^2_{1,0}(q^2/M^2) $ in the squared LO amplitude and phase space measure. The coefficient of the $ F_{2,0}Z $ term seems to reflect the spin-statistic characteristic of the target particle. For fermions, the coefficient is $ 1/4 $, whereas for bosons, it is $ 1/6 $.

      Although the coefficients of $ F^2_{2,0}(\cos\theta+1) $ inside the square brackets explicitly depend on the target spin s, they seem to be expressed as $-\dfrac{1+s}{48s}$, at least for $ s=1/2 $, 1, 3/2, 2. It will be interesting to see whether this parameterization persists for an arbitrary s.

      Analogous to what is done in Section III.A, for a spin-$ 1/2 $ composite target particle, the HPET-based calculation yields the following unpolarized cross section:

      $ \begin{aligned}[b] \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}=& \frac {\pi\alpha^2Z^2}{2{\bf k}^2\sin^4\dfrac{\theta}{2}} -\frac{\pi\alpha^2Z^2}{M |{\bf k}|\sin^2\dfrac{\theta}{2}}+ \frac {\pi\alpha^2}{4M^2\sin^2\dfrac{\theta}{2}} \\&\times \Big[-2c_DZ+c_F^2(\cos\theta+1)+5Z^2(1-\cos\theta)\Big]. \end{aligned} $

      (A4)

      Reassuringly, this EFT result exactly agrees with that obtained from (A3).

    APPENDIX B: RUTHERFORD SCATTERING WITH NONRELATIVISTIC SPINLESS PROJECTILE
    • We can repeat our investigation in Section III.B by replacing the projectile with a light slowly-moving spinless electron. At the lowest order in electron velocity, but to all orders in $ 1/M $, the resulting unpolarized cross section is identical to (11), which is obtained for a spin-${1/2}$ projectile. This is well anticipated because the spin degree of freedom decouples in the nonrelativistic limit.

      At relative order-$ v^2 $, after heavy target mass expansion, the differential unpolarized cross section becomes particularly simple:

      $ \begin{aligned}[b] \left( \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}\right)_{(v^2)}= & \frac {\pi\alpha^2}{2\mathbf{k}^2\sin^2\dfrac{\theta}{2}} \Bigg[\frac{Z^2}{\sin^2\dfrac{\theta}{2}}-\frac{2mZ^2}{M}\\&- \frac {4m^2Z}{M^2}\tilde{f}^s_\text{NNLO}+\mathcal{O}\left( \frac {1}{M^3}\right)\Bigg], \end{aligned} $

      (B1)

      where

      $ \begin{aligned}[b] \tilde{f}^{s}_\text{NNLO}=& 2F'_{1,0}+\Theta\left(s- \frac {1}{2}\right) \frac {5-(-1)^{2s}}{12}F_{2,0} -\Theta\left(s-1\right) \frac {1}{3}F_{1,1}\\&+ \frac {Z}{4}\cos\theta - \frac {1}{6}\left(\lceil s\rceil +s+ \frac {2}{3}\right)Z. \\& \left( s=0, \frac {1}{2},1, \frac {3}{2},2\right) \end{aligned} $

      (B2)

      Clearly, the $ {\cal O}(v^2/M^n) $ ($ n=0,1 $) terms remain universal. At $ {\cal O}(v^2/M^2) $, the universality becomes partially violated. However, the $ F^\prime_{1,0} $, $ F_{1,1} $, and $ Z\cos\theta $ terms still do not depend on the target particle spin. The coefficient of $ F_{2,0} $ seems to reflect the spin-statistic characteristic of the target particle. For fermions, the coefficient is $ 1/2 $, whereas for bosons, it is $ 1/3 $.

      Similar to Section IVB, we can combine NRQED and HPET to study the soft behavior of this type of Rutherford scattering. Since the incident electron is assumed to be spinless, it is natural to work with scalar NRQED plus HPET. Up to the relative order-$ v^2 $, the scalar NRQED Lagrangian reads as

      $ \mathcal{L}_\text{sNRQED}= Q^\dagger\left({\rm i}D^0+d_2 \frac {\mathbf{D}^2}{2m}+d_4 \frac {\mathbf{D}^4}{8m^3}\right)Q, $

      (B3)

      with Q signifying the field that annihilates a nonrelativistic scalar electron. Again, $ d_2=d_4=1 $ is a rigorous consequence of Lorentz symmetry.

      Based on scalar NRQED and HPET, we are able to obtain the following unpolarized Rutherford cross section, which is accurate to the relative order-$ v^2/M^2 $:

      $ \begin{aligned}[b] \frac {\mathrm{d}\sigma}{\mathrm{d}\cos\theta}=& \frac{\pi \alpha^2 m^2Z^2}{2\mathbf{k}^4 \sin^4\dfrac{\theta}{2}} -\frac{\pi \alpha ^2 m^4Z^2}{M^2\mathbf{k}^4} + \frac {Z\pi\alpha^2}{2\mathbf{k}^2\sin^2\dfrac{\theta}{2}}\\&\times\left[ \frac {Z}{\sin^2\dfrac{\theta}{2}} - \frac {2mZ}{M} -\frac{m^2\left(-2c_2Z+c_D+Z\cos\theta+Z\right)}{M^2}\right]. \end{aligned} $

      (B4)

      Reassuringly, this EFT result exactly reproduces the soft limit obtained from full QED, (B2).

      We can further verify that the EFT amplitude reproduces the correct soft behavior, which is deduced by conducting both nonrelativistic and heavy target mass expansion from the full QED amplitude in (1). Working again in Coulomb gauge, and using the following electromagnetic matrix element involving a spinless electron:

      $ \langle k'\vert J^0\vert k\rangle= -1, $

      (B5a)

      $ \langle k'\vert \mathbf{J}\vert k\rangle= (\mathbf{k}+\mathbf{k}')\left[ \frac {d_2}{2m}- \frac {d_4}{8m^3}(\mathbf{k}^2+|\mathbf{k}'|^2)\right], $

      (B5b)

      we can readily obtain the expanded Rutherford amplitude through order-$ v^2/M^2 $, which is indeed compatible with the EFT amplitude.

Reference (13)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return