DV Calculations of 4-momentum components

From New IAC Wiki
Jump to navigation Jump to search

Calculations of 4-momentum components

Initial Conditions

Lab Frame

Lab.png
Figure 1: Definition of variables in the Lab Frame


Begining with the assumption that the incoming electron, p1, has momentum of 11000 MeV in the positive z direction.


[math]\vec p_{1(z)}\equiv \vec p_{1}=11000 MeV \widehat {z}[/math]


We can also assume the Moller electron, p2, is initially at rest


[math]\vec p_{2}\equiv 0[/math]


This gives the total energy in this frame as


[math]E\equiv \sqrt{p^2+m^2}[/math]


[math]E\equiv \sqrt{(p_{1}+p_{2})^2+(m_{1}+m_{2})^2}[/math]


[math]\Longrightarrow E\equiv \sqrt{(11000 MeV)^2+(.511 MeV+.511 MeV)^2}\approx 11000 MeV[/math]

Center of Mass Frame

CM.png
Figure 2: Definition of variables in the Center of Mass Frame


Starting with the definition for the total relativistic energy:


[math]E^2\equiv p^2c^2+m^2c^4[/math]


[math]\Longrightarrow \frac {E^2}{c^2}-p^2=m^2 c^2[/math]


We can use 4-momenta vectors, i.e. [math]{\mathbf P}\equiv \left(\begin{matrix} E/c\\ p_x \\ p_y \\ p_z \end{matrix} \right)[/math] ,with c=1, to describe the variables in the CM Frame.


Using the fact that the scalar product of a 4-momenta with itself,


[math]{\mathbf P_1}\cdot {\mathbf P_1}=E_1E_1-\vec p_1\cdot \vec p_1 =m_{1}^2=s[/math]


is invariant. Since we can assume that the frame of reference is an inertial frame, it moves at a constant velocity, the mass should remain constant.

[math]\frac {d\vec p}{dt}=0\Rightarrow \frac{d(m\vec v)}{dt}=\frac{c\ dm}{dt}\Rightarrow \frac{dm}{dt}=0[/math]


<center[math] \therefore m=const[/math]


Using this notation, the sum of two 4-momenta forms a 4-vector as well

[math]{\mathbf P_1}+ {\mathbf P_2}= \left( \begin{matrix}E_1+E_2\\\vec p_1 +\vec p_2 \end{matrix} \right)= {\mathbf P}[/math]

The length of this four-vector is an invariant as well

[math]{\mathbf P^2}=({\mathbf P_1}+{\mathbf P_2})^2=(E_1+E_2)^2-(\vec p_1 +\vec p_2 )^2=(m_1+m_2)^2=s[/math]



[math]{\mathbf P^2}=({\mathbf P_1}+{\mathbf P_2})^2=(E_1+E_2)^2-(\vec p_1 +\vec p_2 )^2=s[/math]


We can perform a Lorentz transformation to the Center of Mass frame, with zero total momentum


[math]\left( \begin{matrix}E^*_{1}+E^*_{2}\\ 0 \\ 0 \\ 0\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}E_{1}+E_{2}\\ 0 \\ 0 \\ p_{1(z)}+p_{2(z)}\end{matrix} \right)[/math]

Without knowing the values for gamma or beta, we can show that lengths of the two 4-momenta are invariant

[math]s={\mathbf P^*}^2=(E^*_{1}+E^*_{2})^2-(\vec p\ ^*_{1}+\vec p\ ^*_{2})^2=(m_{1}^*+m_{2}^*)^2[/math]


[math]s={\mathbf P}^2=(E_{1}+E_{2})^2-(\vec p_{1}+\vec p_{2})^2=(m_{1}+m_{2})^2[/math]


[math]\Longrightarrow (m_{1}^*+m_{2}^*)^2=(m_{1}+m_{2})^2[/math]


Using the fact that

[math]\begin{cases} m_{1}=m_{2} \\ m_{1}^*=m_{2}^* \end{cases}[/math]


[math](m_{1}^*+m_{1}^*)^2=(m_{1}+m_{1})^2[/math]


[math]2m_{1}^*=2m_{1}[/math]


[math]m_{1}^*=m_{1}[/math]


[math]\Longrightarrow m_{1}=m^*_{1}\ ; m_{2}=m^*_{2}[/math]




For incoming electrons moving only in the z-direction, we can write


[math]{\mathbf P_1}+ {\mathbf P_2}= \left( \begin{matrix}E_1+E_2\\ 0 \\ 0 \\ p_{1(z)}+p_{2(z)}\end{matrix} \right)={\mathbf P}[/math]




Setting these equal to each other, we can use this for the collision of two particles of mass m1 and m2. Since the total momentum is zero in the Center of Mass frame, we can express total energy in the center of mass frame as

[math](E^*_{1}+E^*_{2})^2-(\vec p\ ^*_{1}+\vec p\ ^*_{2})^2=s=(E_{1}+E_{2})^2-(\vec{p_{1}}+\vec p_{2})^2[/math]


[math](E^*)^2-(\vec p\ ^*)^2=(E_{1}+E_{2})^2-(\vec{p_{1}}+\vec p_{2})^2[/math]


[math](E^*)^2=(E_{1}+E_{2})^2-(\vec{p_{1}}+\vec p_{2})^2[/math]


[math]E^*=[(E_{1}+E_{2})^2-(\vec{p_{1}}+\vec p_{2})^2]^{1/2}[/math]


[math]E^*=[E_{1}^2+2E_{1}E_{2}+E_{2}^2-\vec p_{1} . \vec p_{2} -\vec p_{1} . \vec p_{1} -\vec p_{2} . \vec p_{1} -\vec p_{2} . \vec p_{2} ]^{1/2}[/math]


[math]E^*=[(E_{1}^2- p_{1}^2 )+(E_{2}^2-p_{2}^2 )+2E_{1}E_{2}-\vec p_{1} . \vec p_{2} -\vec p_{2} . \vec p_{1} ]^{1/2}[/math]


[math]E^*=[(E_{1}^2- p_{1}^2 )+(E_{2}^2-p_{2}^2 )+2E_{1}E_{2}-p_{1} p_{2}\cos(\theta) - p_{2} p_{1}\cos(\theta) ]^{1/2}[/math]


[math]E^*=[m_{1}^2+m_{2}^2+2E_{1}E_{2}-p_{1} p_{2}\cos(\theta) - p_{2} p_{1}\cos(\theta) ]^{1/2}[/math]


[math]E^*=[m_{1}^2+m_{2}^2+2E_{1}E_{2}-2p_{1} p_{2}\cos(\theta) ]^{1/2}[/math]


[math]E^*=[m_{1}^2+m_{2}^2+2E_{1}E_{2}-2p_{1} p_{2}\cos(\theta) ]^{1/2}[/math]


Using the relations [math]\beta\equiv \vec p/E\Longrightarrow \vec p=\beta E[/math]


[math]E^*=[m_{1}^2+m_{2}^2+2E_{1}E_{2}(1-\beta_{1}\beta_{2}\cos(\theta))]^{1/2}[/math]


where [math] \theta [/math] is the angle between the particles in the Lab frame.



In the frame where one particle (m2) is at rest


[math]\Longrightarrow \beta_{2}=0[/math]


[math]\Longrightarrow p_{2}=0[/math]


which implies,


[math] E_{2}=[p_{2}^2+m_{2}^2]^{1/2}=m_{2}[/math]



[math]E^*=(m_{1}^2+m_{2}^2+2E_{1} m_{2})^{1/2}=(.511MeV^2+.511MeV^2+2(\sqrt{(11000 MeV)^2+(.511 MeV)^2})(.511 MeV))^{1/2}\approx 106.030760886 MeV[/math]

where [math]E_{1}=\sqrt{p_{1}^2+m_{1}^2}\approx 11000 MeV[/math]


Inspecting the Lorentz transformation to the Center of Mass frame:


[math]\left( \begin{matrix}E^*_{1}+E^*_{2}\\ 0 \\ 0 \\ 0\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}E_{1}+E_{2}\\ 0 \\ 0 \\ p_{1(z)}+p_{2(z)}\end{matrix} \right)[/math]


For the case of a stationary electron, this simplifies to:

[math]\left( \begin{matrix} E^* \\ p^*_{x} \\ p^*_{y} \\ p^*_{z}\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^*\gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}E_{1}+m_{2}\\ 0 \\ 0 \\ p_{1(z)}+0\end{matrix} \right)[/math]


which gives,


[math]\Longrightarrow\begin{cases} E^*=\gamma^* (E_{1}+m_{2})-\beta^* \gamma^* p_{1(z)} \\ p^*_{z}=-\beta^* \gamma^*(E_{1}+m_{2})+\gamma^* p^*_{1(z)} \end{cases}[/math]


Solving for [math]\beta^*[/math], with [math]p^*_{z}=0[/math]

[math]\Longrightarrow \beta^* \gamma^*(E_{1}+m_{2})=\gamma^* p_{1(z)}[/math]


[math]\Longrightarrow \beta^*=\frac{p_{1}}{(E_{1}+m_{2})}[/math]


Similarly, solving for [math]\gamma^*[/math] by substituting in [math]\beta^*[/math]


[math]E^*=\gamma^* (E_{1}+m_{2})-\frac{p_{1}}{(E_{1}+m_{2})} \gamma^* p_{1(z)}[/math]


[math]E^*=\gamma^* \frac{(E_{1}+m_2)^2}{(E_{1}+m_{2})}-\gamma^*\frac{(p_{1(z)})^2}{(E_{1}+m_{2})}[/math]


Using the fact that [math]E^*=[(E_{1}+E_{2})^2-(\vec p_{1}+\vec p_{2})^2]^{1/2}[/math]


[math]E^*=\gamma^* \frac{E^*\ ^2}{(E_{1}+m_{2})}[/math]
[math]\Longrightarrow \gamma^*=\frac{(E_1+m_2)} {E^*}[/math]



Using the relation

[math]\left( \begin{matrix} E^*_{1}+E^*_{2} \\ p^*_{1(x)}+p^*_{2(x)} \\ p^*_{1(y)}+p^*_{2(y)} \\ p^*_{1(z)}+p^*_{2(z)}\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}E_{1}+m_{2}\\ 0 \\ 0 \\ p_{1(z)}+0\end{matrix} \right)[/math]


[math]\Longrightarrow \left( \begin{matrix} E^*_{2} \\ p^*_{2(x)} \\ p^*_{2(y)} \\ p^*_{2(z)}\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}m\\ 0 \\ 0 \\ 0\end{matrix} \right)[/math]


[math]\Longrightarrow\begin{cases} E^*_{2}=\gamma^* (m_{2}) \\ p^*_{2(z)}=-\beta^* \gamma^* (m_{2}) \end{cases}[/math]


[math]\Longrightarrow\begin{cases} E^*_{2}=\frac{(E_{1}+m_{2})}{E^*} (m_2) \\ p^*_{2(z)}=-\frac{p_{1}}{(E_{1}+m_2)} \frac{(E_{1}+m_{2})}{E^*} (m_{2}) \end{cases}[/math]


[math]\Longrightarrow\begin{cases} E^*_{2}=\frac{(11000 MeV+.511 MeV)}{106.031 MeV} (.511 MeV) \approx 53.0129177 MeV\\ p^*_{2(z)}=-\frac{11000 MeV}{106.031 MeV} (.511 MeV) \approx -53.013 MeV \end{cases}[/math]


[math]\Longrightarrow \left( \begin{matrix} E^*_{1}\\ p^*_{1(x)} \\ p^*_{1(y)} \\ p^*_{1(z)}\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}E_{1}\\ 0 \\ 0 \\ p_{1(z)}\end{matrix} \right)[/math]


[math]\Longrightarrow\begin{cases} E^*_{1}=\gamma^* (E_{1})-\beta^* \gamma^* p_{1(z)} \\ p^*_{1(z)}=-\beta^* \gamma^*(E_{1})+\gamma^* p_{1(z)} \end{cases}[/math]


[math]\Longrightarrow\begin{cases} E^*_{1}=\frac{(E_{1}+m_{2})}{E^*} (E_{1})-\frac{p_{1(z)}}{(E_{1}+m_2)} \frac{(E_{1}+m_{2})}{E^*} p_{1(z)} \\ p^*_{1(z)}=-\frac{p_{1}}{(E_{1}+m_2)} \frac{(E_{1}+m_{2})}{E^*}(E_{1})+\frac{(E_{1}+m_{2})}{E^*} p_{1(z)} \end{cases}[/math]


[math]\Longrightarrow\begin{cases} E^*_{1}=\frac{(11000 MeV+.511 MeV)}{106.031 MeV} (11000 MeV)-\frac{11000 MeV}{106.031 MeV} 11000 MeV \approx 53.013 MeV\\ p^*_{1(z)}=-\frac{11000 MeV}{106.031 MeV}(11000 MeV)+\frac{(11000 MeV+.511 MeV)}{106.031 MeV} 11000 MeV \approx 53.013 MeV \end{cases}[/math]


[math]p^*_{1} =\sqrt {(p^*_{1(x)})^2+(p^*_{1(y)})^2+(p^*_{1(z)})^2} \Longrightarrow p^*_{1}=p^*_{1(z)}[/math]


[math]p^*_{2} =\sqrt {(p^*_{2(x)})^2+(p^*_{2(y)})^2+(p^*_{2(z)})^2} \Longrightarrow p^*_{2}=p^*_{2(z)}[/math]


This gives the momenta of the particles in the center of mass to have equal magnitude, but opposite directions.

[math] p^*_{1}=p^*_{2}\approx 53.013 MeV[/math]



[math]E^*_{1}= E^*_{2}\approx 53.013 MeV[/math]

Final Conditions

Moller electron Lab Frame

Finding the correct kinematic values starting from knowing the momentum of the Moller electron, [math]p^'_{2}[/math] , in the Lab frame,

Xz lab.png
Figure 3: Definition of Moller electron variables in the Lab Frame in the x-z plane.
Using [math]\theta '_2=\arccos \left(\frac{p^'_{2(z)}}{p^'_{2}}\right)[/math]


[math]\Longrightarrow {p^'_{2(z)}=p^'_{2}\cos(\theta '_2)}[/math]

^'


Checking on the sign resulting from the cosine function, we are limited to:

[math]0^\circ \le \theta '_2 \le 60^\circ \equiv 0 \le \theta '_2 \le 1.046\ Radians[/math]

Since,

[math]\frac{p^'_{2(z)}}{p^'_{2}}=cos(\theta '_2)[/math]


[math]\Longrightarrow p^'_{2(z)}\ should\ always\ be\ positive[/math]
Xy lab.png
Figure 4: Definition of Moller electron variables in the Lab Frame in the x-y plane.
Similarly, [math]\phi '_2=\arccos \left( \frac{p^'_{2(x) Lab}}{p^'_{2(xy)}} \right)[/math]


where [math]p_{2(xy)}^'=\sqrt{(p_{2(x)}^')^2+(p^'_{2(y)})^2}[/math]


[math](p^'_{2(xy)})^2=(p^'_{2(x)})^2+(p^'_{2(y)})^2[/math]


and using [math]p^2=p_{(x)}^2+p_{(y)}^2+p_{(z)}^2[/math]


this gives [math](p^'_{2})^2=(p^'_{2(xy)})^2+(p^'_{2(z)})^2[/math]


[math]\Longrightarrow (p'_{2})^2-(p'_{2(z)})^2 = (p'_{2(xy)})^2[/math]


[math]\Longrightarrow p_{2(xy)}^'=\sqrt{(p^'_{2})^2-(p^'_{2(z)})^2}[/math]


which gives[math]\phi '_2 = \arccos \left( \frac{p_{2(x)}'}{\sqrt{p_{2}^{'\ 2}-p_{2(z)}^{'\ 2}}}\right)[/math]
[math]\Longrightarrow p_{2(x)}'=\sqrt{p_{2}^{'\ 2}-p_{2(z)}^{'\ 2}} \cos(\phi)[/math]


Similarly, using [math]p_{2}^2=p_{2(x)}^2+p_{2(y)}^2+p_{2(z)}^2[/math]


[math]\Longrightarrow p_{2}^{'\ 2}-p_{2(x)}^{'\ 2}-p_{2(z)}^{'\ 2}=p_{2(y)}^{'\ 2}[/math]
[math]p_{2(y)}'=\sqrt{p_{2}^{'\ 2}-p_{2(x)}^{'\ 2}-p_{2(z)}^{'\ 2}}[/math]


Checking on the sign from the cosine results for [math]\phi '_2[/math]

[math]-\pi \le \phi '_2 \le \pi\ Radians[/math]


[math]For\ 0 \le \phi '_2 \le \frac{-\pi}{2}\ Radians[/math]
x=POSITIVE
y=NEGATIVE
[math]For\ 0 \le \phi '_2 \le \frac{\pi}{2}\ Radians[/math]
x=POSITIVE
y=POSITIVE
[math]For\ \frac{-\pi}{2} \le \phi '_2 \le -\pi\ Radians[/math]
x=NEGATIVE
y=NEGATIVE
[math]For\ \frac{\pi}{2} \le \phi '_2 \le \pi\ Radians[/math]
x=NEGATIVE
y=POSITIVE

Moller electron Center of Mass Frame

Relativistically, the x and y components remain the same in the conversion from the Lab frame to the Center of Mass frame, since the direction of motion is only in the z direction.


[math]p^*_{2(x)}\Leftrightarrow p_{2(x)}'[/math]


[math]p^*_{2(y)}\Leftrightarrow p_{2(y)}'[/math]


[math]p^*_{2(z)}=\sqrt {(p^*_2)^2-(p^*_{2(x)})^2-(p^*_{2(y)})^2}[/math]



[math]\Longrightarrow E^*_{2}=\approx 53.013 MeV[/math]

Electron Center of Mass Frame

Relativistically, the x, y, and z components have the same magnitude, but opposite direction, in the conversion from the Moller electron's Center of Mass frame to the electron's Center of Mass frame.


[math]p^*_{2(x)}= -p^*_{1(x)}[/math]


[math]p^*_{2(y)}= -p^*_{1(y)}[/math]


[math]p^*_{2(z)}=-p^*_{1(z)}[/math]

where previously it was shown

[math]p^*_{1}\approx 53.013 MeV[/math]



[math]E^*_{1}\approx 53.013 MeV[/math]

Electron Lab Frame

We can perform a Lorentz transformation from the Center of Mass frame, with zero total momentum, to the Lab frame.


[math]\left( \begin{matrix}E'_{1}+E'_{2}\\ p'_{1(x)}+p'_{2(x)} \\p'_{1(y)}+ p'_{2(y)} \\ p'_{1(z)}+p'_{2(z)}\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & \beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ \beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}E^*_{1}+E^*_{2}\\ p^*_{1(x)}+p^*_{2(x)} \\ p^*_{(1(y)}+p^*_{2(y)} \\ p^*_{1(z)}+p^*_{2(z)}\end{matrix} \right)[/math]


[math]p^*_{1(x)}= p'_{1(x)}[/math]


[math]p^*_{1(y)}= p'_{1(y)}[/math]


[math]\Longrightarrow\begin{cases} E'=\gamma E^*+\beta \gamma p^*_{z} \\ p'_{z}=\beta \gamma E^*+ \gamma p^*_{z} \end{cases}[/math]


[math]\Longrightarrow\begin{cases} \gamma = \frac{E'}{E^*} \\ \beta =\frac{p'_{z}}{E'} \end{cases}[/math]



Without knowing the values for gamma or beta, we can use the fact that lengths of the two 4-momenta are invariant

[math]s={\mathbf P^*} ^2=(E^*_{1}+E^*_{2})^2-(\vec p\ ^*_{1}+\vec p\ ^*_{2})^2[/math]


[math]s={\mathbf P}^2=(E'_{1}+E'_{2})^2-(\vec p\ '_{1}+\vec p\ '_{2})^2[/math]


Setting these equal to each other, we can use this for the collision of two particles of mass m1 and m2. Since the total momentum is zero in the Center of Mass frame, we can express total energy in the center of mass frame as

[math](E^*_{1}+E^*_{2})^2-(\vec p\ ^*_{1}+\vec p\ ^*_{2})^2=s=(E'_{1}+E'_{2})^2-(\vec p\ '_{1}+\vec p\ '_{2})^2[/math]


[math](E^*)^2-(\vec p\ ^*)^2=(E'_{1}+E'_{2})^2-(\vec p\ ')^2[/math]


[math](E^*)^2=(E')^2-(\vec p\ ')^2[/math]


[math](E^*)^2+(\vec{p'})^2=(E')^2[/math]


[math]E'=\sqrt{(E^*)^2+(\vec p\ ')^2}[/math]

Since momentum is conserved, the initial momentum from the incident electron and stationary electron is still the same in the Lab frame, therefore [math]p\ '_{Lab}=11000 MeV[/math]


[math]E'=\sqrt{(106.031 MeV)^2+(11000 MeV)^2}\approx 11000.511 MeV[/math]


This is the same as the initial total energy in the Lab frame, which should be expected since scattering is considered to be an elastic collision.


Since,

[math]E'\equiv E'_{1}+E'_{2}[/math]


[math]\Longrightarrow E'_{1}=E'-E'_{2}[/math]

Using the relation

[math]E'^2\equiv p^2+m^2[/math]


[math]\Longrightarrow p_{1}^{'\ 2}=E_{1}^{'\ 2}-m_{1}^2=E_{1}^{'\ 2}-(.511 MeV)^2[/math]


Using

[math]p^2 \equiv p_x^2+p_y^2+p_z^2[/math]


[math]\Longrightarrow p_{1(z)}'=\sqrt {p_{1}^{'\ 2}-p_{1(x)}^{'\ 2}-p_{1(y)}^{'\ 2}}[/math]


DV_RunGroupC_Moller#Momentum distributions in the Center of Mass Frame

DV_MollerTrackRecon#Calculations_of_4-momentum_components