Recoil Proton Momentum Spectrum in Neutron Decay

Waleed Khalid
Messages
2
Reaction score
0
I wish to draw the proton momentum spectrum by transforming the energy spectrum of recoil protons. I have calculated the energy spectrum using Nachtmann's spectrum: wp=g1[T]+a*g2[T]
Where:
g1[T]=(1 - x2/σ[T])2 * Sqrt[1 - σ[T]] * (4*(1 + x2/σ[T]) - (4/3*(σ[T] - x2)/σ[T])*(1 - σ[T]));
g2[T]=(1 - x2/σ[T])2 * Sqrt[1 - σ[T]] * (4*(1 + x2/σ[T] - 2*σ[T]) - 4/3*(σ[T] - x2)/σ[T]*(1 - σ[T]));
and σ[T]=1 - 2 * T * mn/(mn-mp)2
and a is the electron neutrino correlation.

To get the momentum spectrum, I am transforming the functions (non relativistically):
TofP[p]=p2/(2*mn)

wmom=wp[TofP[p]]

However this doesn't yield the correct spectrum for the momentum of the recoiled protons, as far as I have gotten it I have to multiply it with p and TofP[p] to get the shape of the correct spectrum (wmom=wp[TofP[p]]*p*TofP[p]). Which doesn't make sense to me, so if anyone can explain I would be highly thankful.
 
Physics news on Phys.org
Nevermind, I was being stupid, the answer was simple, since I was using a non relativistic conversion I had to multiply with a factor of p/mn

Waleed Khalid said:
I wish to draw the proton momentum spectrum by transforming the energy spectrum of recoil protons. I have calculated the energy spectrum using Nachtmann's spectrum: wp=g1[T]+a*g2[T]
Where:
g1[T]=(1 - x2/σ[T])2 * Sqrt[1 - σ[T]] * (4*(1 + x2/σ[T]) - (4/3*(σ[T] - x2)/σ[T])*(1 - σ[T]));
g2[T]=(1 - x2/σ[T])2 * Sqrt[1 - σ[T]] * (4*(1 + x2/σ[T] - 2*σ[T]) - 4/3*(σ[T] - x2)/σ[T]*(1 - σ[T]));
and σ[T]=1 - 2 * T * mn/(mn-mp)2
and a is the electron neutrino correlation.

To get the momentum spectrum, I am transforming the functions (non relativistically):
TofP[p]=p2/(2*mn)

wmom=wp[TofP[p]]

However this doesn't yield the correct spectrum for the momentum of the recoiled protons, as far as I have gotten it I have to multiply it with p and TofP[p] to get the shape of the correct spectrum (wmom=wp[TofP[p]]*p*TofP[p]). Which doesn't make sense to me, so if anyone can explain I would be highly thankful.
 
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top