Exercise 26
In Lyttleton and Bondi suggested that the expansion of the Universe could be explained if matter carried a net charge. Suppose that the Universe is made up of hydrogen atoms with a number density , which is maintained a constant. Let the charge on the proton be: where is the electronic charge.
(a)
Find the critical value of y such that expansion may start.
(b)
Show that the velocity of expansion is proportional to the distance from the centre.
NCERT’s answer
(a)
Let the Universe have a radius \(\displaystyle R\). Assume that the hydrogen atoms are uniformly distributed. The charge on each hydrogen atom is \(\displaystyle e_H=-(1+y) e+e=-y e=|y e|\)
The mass of each hydrogen atom is \(\displaystyle \sim m_p\) (mass of proton). Expansion starts if the Coulumb repulsion on a hydrogen atom, at \(\displaystyle R\), is larger than the gravitational attraction.
Let the Electric Field at \(\displaystyle R\) be \(\displaystyle \mathbf{E}\). Then
\[\begin{aligned}
& 4 \pi R^2 E=\frac{4}{3 \varepsilon_{\mathrm{o}}} \pi \mathrm{R}^3 \mathrm{~N}|\mathrm{ye}| \text { (Gauss's law) } \\
& \mathbf{E}(\mathrm{R})=\frac{1}{3} \frac{\mathrm{~N}|\mathrm{ye}|}{\varepsilon_{\mathrm{o}}} \mathrm{R} \hat{\mathbf{r}}
\end{aligned}
\]
Let the gravitational field at \(\displaystyle R\) be \(\displaystyle G_{\mathrm{R}}\). Then
\[\begin{aligned}
& -4 \pi R^2 \mathrm{G}_R=4 \pi G \mathrm{~m}_{\mathrm{p}} \frac{4}{3} \pi \mathrm{R}^3 \mathrm{~N} \\
& \mathrm{G}_{\mathrm{R}}=-\frac{4}{3} \pi \mathrm{Gm}_\rho \mathrm{NR} \\
& \mathbf{G}_{\mathrm{R}}(\mathbf{R})=-\frac{4}{3} \pi \mathrm{Gm}_\rho N R \hat{\mathbf{r}}
\end{aligned}
\]
Thus the Coulombic force on a hydrogen atom at \(\displaystyle R\) is
\[\mathrm{ye} \boldsymbol{E}(\mathrm{R})=\frac{1}{3} \frac{\mathrm{Ny}^2 \mathrm{e}^2}{\varepsilon_{\mathrm{o}}} \mathrm{R} \hat{\mathbf{r}}
\]
The gravitional force on this atom is
\[\mathrm{m}_{\mathrm{p}} \mathbf{G}_{\mathrm{R}}(\mathrm{R})=-\frac{4 \pi}{3} \mathrm{GNm}_{\mathrm{p}}^2 \mathrm{R} \hat{\mathbf{r}}
\]
The net force on the atom is
\[\mathbf{F}=\left(\frac{1}{3} \frac{\mathrm{Ny}^2 \mathrm{e}^2}{\varepsilon_{\mathrm{o}}} \mathrm{R}-\frac{4 \pi}{3} \mathrm{GNm}_{\mathrm{p}}^2 \mathrm{R}\right) \hat{\mathbf{r}}
\]
The critical value is when
\[\begin{aligned}
& \frac{1}{3} \frac{N y_{\mathrm{C}}^2 \mathrm{e}^2}{\varepsilon_{\mathrm{o}}} R=\frac{4 \pi}{3} G N m_{\mathrm{p}}^2 R \\
& \Rightarrow y_{\mathrm{c}}^2=4 \pi \varepsilon_{\mathrm{o}} \mathrm{G} \frac{\mathrm{~m}_{\mathrm{p}}^2}{\mathrm{e}^2} \\
& \simeq \frac{7 \times 10^{-11} \times 1.8^2 \times 10^6 \times 81 \times 10^{-62}}{9 \times 10^9 \times 1.6^2 \times 10^{-38}} \\
& \sim 63 \times 10^{-38} \\
& \therefore y_{\mathrm{C}} \sim 8 \times 10^{-19} \sim 10^{-18}
\end{aligned}
\]
(b)
Because of the net force, the hydrogen atom experiences an acceleration such that
\[m_p \frac{d^2 R}{d t^2}=\left(\frac{1}{3} \frac{N y^2 e^2}{e_o} R-\frac{4 p}{3} G N m_p^2 R\right)
\]
Or, \(\displaystyle \frac{d^2 R}{d t^2}=a^2 R\) where \(\displaystyle \alpha^2=\frac{1}{m_p}\left(\frac{1}{3} \frac{N y^2 e^2}{e_o}-\frac{4 p}{3} G N m_p^2\right)\)
This has a solution \(\displaystyle R=A e^{a t}+B e^{-a t}\)
As we are seeking an expansion, \(\displaystyle B=0\).
\[\begin{aligned}
& \therefore R=A e^{\alpha t} \\
& \Rightarrow \dot{R}=\alpha A e^{\alpha t}=\alpha R
\end{aligned}
\]
Thus, the velocity is proportional to the distance from the centre.
About the centre, an atom at distance \(\displaystyle r\) feels only the \(\displaystyle \tfrac43\pi r^3N\) atoms inside radius \(\displaystyle r\) (Gauss's law, and its gravitational analogue).
Atom mass \(\displaystyle \approx m_p\); with \(\displaystyle e\) the electron's charge, the net charge per atom is
\[q_H=e_p+e=-ye \]
\[F_e=\frac{1}{4\pi\varepsilon_0}\,\frac{\left(\tfrac43\pi r^3N\,ye\right)(ye)}{r^2}=\frac{4\pi}{3}Nr\,\frac{y^2e^2}{4\pi\varepsilon_0} \]
\[F_g=\frac{G\left(\tfrac43\pi r^3N\,m_p\right)m_p}{r^2}=\frac{4\pi}{3}Nr\,Gm_p^2 \]
(a) Expansion starts when \(\displaystyle F_e=F_g\); \(\displaystyle r\) cancels:
\[\frac{y^2e^2}{4\pi\varepsilon_0}=Gm_p^2\;\Rightarrow\;y_c=\frac{m_p}{e}\sqrt{4\pi\varepsilon_0G} \]
\[y_c=\frac{1.67\times10^{-27}\,\mathrm{kg}}{1.6\times10^{-19}\,\mathrm{C}}\sqrt{\frac{6.67\times10^{-11}}{9\times10^{9}}}\ \mathrm{C\,kg^{-1}}\approx0.9\times10^{-18} \]
(b) For \(\displaystyle y>y_c\) the net outward force on the atom at distance \(\displaystyle r\) gives
\[m_p\ddot r=\frac{4\pi}{3}Nr\left[\frac{y^2e^2}{4\pi\varepsilon_0}-Gm_p^2\right]\;\Rightarrow\;\ddot r=\kappa^2r,\qquad \kappa^2=\frac{4\pi N}{3m_p}\left[\frac{y^2e^2}{4\pi\varepsilon_0}-Gm_p^2\right] \]
\(\displaystyle N\) constant \(\displaystyle \Rightarrow\) \(\displaystyle \kappa\) constant, and the general solution is
\[r=Ae^{\kappa t}+Be^{-\kappa t} \]
Expansion means the decaying term must go, \(\displaystyle B=0\):
\[r=Ae^{\kappa t}\;\Rightarrow\;v=\dot r=\kappa Ae^{\kappa t}=\kappa r \]
Answer: \(\displaystyle y_c=\dfrac{m_p}{e}\sqrt{4\pi\varepsilon_0G}\approx10^{-18}\); \(\displaystyle v=\kappa r\), so \(\displaystyle v\propto r\).