A New Electron Charge Distribution Function of Electron ()

Teruo Kurai^{}

Unoki, Ota-ku, Tokyo, Japan.

**DOI: **10.4236/jmp.2023.142008
PDF
HTML XML
94
Downloads
483
Views
Citations

Unoki, Ota-ku, Tokyo, Japan.

Here we derive a new charge distribution function for an electron by using as an equation of motion a segment of charge whose self energy interaction is due to electric field potential. Our method is based on the consideration that a charged distribution function should be represented as an eigenfunction of electron mass energy. We compare our electron charge distribution function to that of Weinberg’s *η*(*r*) and our charged electron radius to that obtained by Kim.

Keywords

Share and Cite:

Kurai, T. (2023) A New Electron Charge Distribution Function of Electron. *Journal of Modern Physics*, **14**, 111-126. doi: 10.4236/jmp.2023.142008.

1. Introduction

The theory that particle mass is generated by Higgs mechanism is widely accepted. Especially after the discovery of the Higgs particle [1], the Higgs mechanism has been convinced as the true origin of Bosons. Encouraged by this discovery, many physicists have been trying to show that the Higgs mechanism is also true for fermions, specifically by showing a mass hierarchy for quarks and lepton [2]. However, some physicists are still not convinced that the Higgs mechanism is true for fermions, primarily because of lepton consideration. Quantum Electrodynamics (QED) has been used to describe the interaction between leptons, especially electrons, and it is well known that it describes such electromagnetic phenomena quite well. The only ambiguity appears in the determination of self-energy, and physicists disagree on how to consider this ambiguity. Some physicists argue that to describe leptons, especially electron, QED is sufficient even though this ambiguity exists.

For instance, Weinberg shows an estimated charged radius in his famous book Quantum Theory of Field I [3]. He recently also published the paper [4] in which he opposed the consideration in some leading approximations that only quarks and leptons with nonzero masses are third generation with the other lepton and quark masses arising from some sort of radiative correction. Weinberg points out that this kind of theory was firstly proposed by S.M. Barr and A. Zee [5] and recently has been revived [6]. The opposition by Weinberg seems to show that he still favors the same scheme mentioned in his book. The motivation of this paper is also based on the work of Weinberg. We investigate an electron charged radius by considering that a charge distribution function is described by a segment of charge function. By constructing a segment of charge function as an eigenfunction of electron mass energy, we find that the electron charge distribution function appears simultaneously when the electron gains mass energy by the absorption of an infinite number of photons, as described by Weinberg.

2. Formalism

In order to obtain a charge distribution function, we introduce a charge matrix state as

${\rho}_{\lambda \eta}\left(t,\stackrel{\to}{r}\right)=\left(-e\right)\langle 0|{q}_{\lambda}\left(t,\stackrel{\to}{r}\right){q}_{\eta}^{\u2020}\left(t,\stackrel{\to}{r}\right)|\text{chargestate}\rangle =\left(-e\right){\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(t,\stackrel{\to}{r}\right)$ (1)

where $\lambda ,\eta $ are Dirac indices. This kind of charge matrix (only operator part) is also used by Karnieli [7].

Since electron charge is invariable over time, time derivative of this quantity must be zero.

Dirac equation with mass less field is represented as

${\gamma}^{\mu}{D}_{\mu}q=0$ (2)

where ${D}_{\mu}={\partial}_{\mu}-ie{A}_{\mu}$.

We employ $\gamma $ matrices following Weinberg way [8] as follows

${\gamma}^{0}=-i\left(\begin{array}{cc}0& {\sigma}_{0}\\ {\sigma}_{0}& 0\end{array}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\gamma}^{k}=-i\left(\begin{array}{cc}0& {\sigma}_{k}\\ -{\sigma}_{k}& 0\end{array}\right)$

Here, ${\sigma}_{0}$ is a unit matrix of 2 × 2 matrix and ${\sigma}_{k}$ are the 2 × 2 Pauli matrices ( $k=1,2,3$ ), ${\alpha}^{\mu}={\gamma}^{0}{\gamma}^{\mu}$ ( $\mu =0,1,2,3$ ) and $\beta =i{\gamma}^{0}$.

From Equation (2) we can obtain following equation as

$i\frac{\partial q}{\partial t}=-i{\alpha}^{k}{\partial}_{k}q-e{\alpha}^{\mu}{A}_{\mu}q$ (3)

$i\frac{\partial {q}^{\u2020}}{\partial t}=-i{\alpha}^{k}{\partial}_{k}{q}^{\u2020}+e{\alpha}^{\mu}{A}_{\mu}{q}^{\u2020}$ (4)

where $\mu =0,1,2,3$.

Note that we use the same metric system that Weinberg employs in book [8]

${\eta}^{00}=-1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\eta}^{11}={\eta}^{22}={\eta}^{33}=1$

Since we consider the equation of motion for a segment of charge function ${\stackrel{\xaf}{\rho}}_{\lambda \eta}$, we consider the case that ${\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(t,\stackrel{\to}{r}\right)=\mathrm{exp}\left(-i{P}_{0}t\right){\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(\stackrel{\to}{r}\right)$.

Then from Equation (3) and Equation (4) we obtain

${P}_{0}{\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(\stackrel{\to}{r}\right)=-i{\alpha}^{k}{\partial}_{k}{\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(\stackrel{\to}{r}\right)$ (5)

Here we use the following notation,

$-i{\alpha}^{k}{\partial}_{k}{\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(\stackrel{\to}{r}\right)=\langle 0|-i\left({\alpha}^{k}{\partial}_{k}{q}_{\lambda}\left(\stackrel{\to}{r}\right)\right){q}_{\eta}^{\u2020}\left(\stackrel{\to}{r}\right)+{q}_{\lambda}\left(\stackrel{\to}{r}\right)\left({q}_{\eta}^{\u2020}{\alpha}^{k}{\stackrel{\leftarrow}{\partial}}_{k}\left(-i\right)\right)|\text{chargestate}\rangle $

Here $|0\rangle $ denotes vacuum state.

We interpolate that Equation (5) implies free motion of a segment of charge. There for to obtain an appropriate equation of motion, we have to also consider self-interaction energy. First though, we decompose ${\stackrel{\xaf}{\rho}}_{\lambda \eta}$ to Lorentz invariant form and consider that a segment of charge is depend on only the radius $r=\left|\stackrel{\to}{r}\right|$, then ${\stackrel{\xaf}{\rho}}_{\alpha \beta}$ is described as

${\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(r\right)=I{\stackrel{\xaf}{\rho}}_{0}\left(r\right)+\left(-i\stackrel{\to}{\alpha}\cdot \stackrel{^}{r}\right){\stackrel{\xaf}{\rho}}_{1}\left(r\right)+\beta {\stackrel{\xaf}{\rho}}_{2}\left(r\right)+\beta \left(i\stackrel{\to}{\alpha}\cdot \stackrel{^}{r}\right){\stackrel{\xaf}{\rho}}_{3}\left(r\right)$ (6)

Recalling that actual charge is a scalar quantity, a segment of charge function which is affected by self-interaction energy should be defined as

$\text{actualcharge}=\frac{1}{4}Tr{\rho}_{\lambda \eta}\left(r\right)=\left(-e\right)\frac{1}{4}Tr{\stackrel{\xaf}{\rho}}_{\lambda \eta}\left(r\right)=\left(-e\right){\stackrel{\xaf}{\rho}}_{0}\left(r\right)$ (7)

where *Tr* denotes Trace of matrix.

To consider self-interaction energy, we use Gauss’s Law as

$\text{div}\stackrel{\to}{E}=\frac{{\rho}_{0}}{{\u03f5}_{0}}$ (8)

Since we consider only *r* dependent case, thus volume integral using Green’s theorem gives

${E}_{r}\left(r\right)=\frac{1}{4\pi {\u03f5}_{0}}\frac{4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\rho}_{0}\left({r}^{\prime}\right)}{{r}^{2}}=\frac{1}{4\pi {\u03f5}_{0}}\frac{\left(-e\right)4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)}{{r}^{2}}$ (9)

Then we obtain an electric field potential as

${V}_{r}\left(r\right)=-\frac{-e}{4\pi {\u03f5}_{0}{a}_{0}}+\frac{-e}{4\pi {\u03f5}_{0}}\frac{4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)}{{r}^{\prime}}+\frac{-e}{4\pi {\u03f5}_{0}}4\pi {\displaystyle {\int}_{0}^{{a}_{0}}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{r}^{\prime}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)$ (10)

Note that to obtain Equation (8) we use the setting consideration that charge vanishes at and beyond $r={a}_{0}$.

This means that

$\begin{array}{l}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(r>{a}_{0}\right)\\ {\stackrel{\xaf}{\rho}}_{0}\left(r\right)\ne 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(r\le {a}_{0}\right)\end{array}$ (11)

Then we have to set the following two conditions.

The first condition is

$4\pi {\displaystyle {\int}_{0}^{{a}_{0}}\text{d}r}\text{\hspace{0.05em}}{r}^{2}{\rho}_{0}\left(r\right)=-e$ (12)

Rearrangement gives

${\int}_{0}^{{a}_{0}}\text{d}r}\text{\hspace{0.05em}}{r}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)=\frac{1}{4\pi$ (13)

The second condition is

${\rho}_{0}\left({a}_{0}\right)={\stackrel{\xaf}{\rho}}_{0}\left({a}_{0}\right)=0$ (14)

From Equation (10), the self-interaction energy is described as

$\begin{array}{l}{V}_{e}=\left(-e\right){\stackrel{\xaf}{\rho}}_{0}\left(r\right){V}_{r}\left(r\right)\\ =-\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)+\frac{{e}^{2}}{4\pi {\u03f5}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\left[\frac{4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)}{{r}^{\prime}}+4\pi {\displaystyle {\int}_{r}^{{a}_{0}}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{r}^{\prime}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)\right]\end{array}$ (15)

Then the equation of motion of a segment of charge function becomes

$\begin{array}{c}{P}_{0}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)=-\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{1}\left(r\right)-\frac{2}{r}{\stackrel{\xaf}{\rho}}_{1}\left(r\right)+[-\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\begin{array}{c}\underset{}{\overset{}{\text{\hspace{0.05em}}}}\\ \text{\hspace{0.05em}}\end{array}\\ \text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{{e}^{2}}{4\pi {\u03f5}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\left[\frac{4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)}{{r}^{\prime}}+4\pi {\displaystyle {\int}_{r}^{{a}_{0}}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{r}^{\prime}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)\right]]\end{array}$ (16)

${P}_{0}{\stackrel{\xaf}{\rho}}_{1}\left(r\right)=\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)$ (17)

${P}_{0}{\stackrel{\xaf}{\rho}}_{2}\left(r\right)=-\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{3}\left(r\right)-\frac{2}{r}{\stackrel{\xaf}{\rho}}_{3}\left(r\right)$ (18)

${P}_{0}{\stackrel{\xaf}{\rho}}_{3}\left(r\right)=\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{2}\left(r\right)$ (19)

Derivation of kinetic term of each equation above is shown in Appendix.

Note that because we investigate electron charge distribution function by considering the equation of motion of its segment of charge with self-interaction energy, Equations (16)-(19) represent an electron in the rest frame. Also note that these equations are described in units of $\hslash =c=1$.

To obtain a proper unit description of the charge distribution of an electron, we rewrite Equations (16)-(19) as follows.

$\begin{array}{l}{P}_{0}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)=-\hslash c\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{1}\left(r\right)-\hslash c\frac{2}{r}{\stackrel{\xaf}{\rho}}_{1}\left(r\right)+[-\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\begin{array}{c}\underset{}{\overset{}{\text{\hspace{0.05em}}}}\\ \text{\hspace{0.05em}}\end{array}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{{e}^{2}}{4\pi {\u03f5}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\left[\frac{4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)}{{r}^{\prime}}+4\pi {\displaystyle {\int}_{r}^{{a}_{0}}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{r}^{\prime}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)\right]]\end{array}$ (20)

${P}_{0}{\stackrel{\xaf}{\rho}}_{1}\left(r\right)=\hslash c\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)$ (21)

${P}_{0}{\stackrel{\xaf}{\rho}}_{2}\left(r\right)=-\hslash c\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{3}\left(r\right)-\hslash c\frac{2}{r}{\stackrel{\xaf}{\rho}}_{3}\left(r\right)$ (22)

${P}_{0}{\stackrel{\xaf}{\rho}}_{3}\left(r\right)=\hslash c\frac{\partial}{\partial r}{\stackrel{\xaf}{\rho}}_{2}\left(r\right)$ (23)

Because we are interested in only charge distribution, we need only to obtain the solution of ${\stackrel{\xaf}{\rho}}_{0}\left(r\right)$, however, this solution must correspond to a certain eigenvalue that should be the rest mass of an electron.

To do this, we need to consider only Equation (20) and Equation (21). By inserting the description of ${\stackrel{\xaf}{\rho}}_{1}\left(r\right)$ from Equation (21) into Equation (20), we can obtain the equation for ${\stackrel{\xaf}{\rho}}_{0}\left(r\right)$ as follows.

$\begin{array}{l}{P}_{0}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)=-\frac{{\left(\hslash c\right)}^{2}}{{P}_{0}}\frac{{\partial}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)}{\partial {r}^{2}}-\frac{{\left(\hslash c\right)}^{2}}{{P}_{0}}\frac{2}{r}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(r\right)}{\partial r}+[-\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\begin{array}{c}\underset{}{\overset{}{\text{\hspace{0.05em}}}}\\ \text{\hspace{0.05em}}\end{array}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{{e}^{2}}{4\pi {\u03f5}_{0}}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\left[\frac{4\pi {\displaystyle {\int}_{0}^{r}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{{r}^{\prime}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)}{{r}^{\prime}}+4\pi {\displaystyle {\int}_{r}^{{a}_{0}}\text{d}{r}^{\prime}}\text{\hspace{0.05em}}{r}^{\prime}{\stackrel{\xaf}{\rho}}_{0}\left({r}^{\prime}\right)\right]]\end{array}$ (24)

Because we did not find out exact treatment of the integral part of Equation (24), we adopt a Tayler expansion around $r=0$. This requires us to make a variable dimensionless.

Thus, we change a variable as $\stackrel{\xaf}{r}=\frac{r}{{a}_{0}}$. Then Equation (24) becomes

$\begin{array}{c}{\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=-\frac{{\partial}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial {\stackrel{\xaf}{r}}^{2}}-\frac{2}{\stackrel{\xaf}{r}}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}\\ \text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{{P}_{0}}{{\left(\frac{\hslash c}{{a}_{0}}\right)}^{2}}[-\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}+\frac{{e}^{2}}{4\pi {\u03f5}_{0}}\frac{4\pi {a}_{0}^{2}{\displaystyle {\int}_{0}^{\stackrel{\xaf}{r}}\text{d}\stackrel{\xaf}{{r}^{\prime}}}\text{\hspace{0.05em}}{\stackrel{\xaf}{{r}^{\prime}}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{{r}^{\prime}}\right)}{\stackrel{\xaf}{r}}\\ \begin{array}{c}\underset{}{\overset{}{\text{\hspace{0.05em}}}}\\ \text{\hspace{0.05em}}\end{array}+\frac{{e}^{2}}{4\pi {\u03f5}_{0}}\left[4\pi {a}_{0}^{2}{\displaystyle {\int}_{0}^{1}\text{d}\stackrel{\xaf}{{r}^{\prime}}}\text{\hspace{0.05em}}{\stackrel{\xaf}{{r}^{\prime}}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{{r}^{\prime}}\right)+4\pi {a}_{0}^{2}{\displaystyle {\int}_{0}^{\stackrel{\xaf}{r}}\text{d}\stackrel{\xaf}{{r}^{\prime}}}\text{\hspace{0.05em}}\stackrel{\xaf}{{r}^{\prime}}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{{r}^{\prime}}\right)\right]]\end{array}$ (25)

Note that the condition Equation (13) becomes

$\int}_{0}^{1}\text{d}\stackrel{\xaf}{r}}\text{\hspace{0.05em}}{\stackrel{\xaf}{r}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=\frac{1}{4\pi {a}_{0}^{3$ (26)

To make dimension correct, we need to set one more condition as

$\int}_{0}^{1}\text{d}\stackrel{\xaf}{r}}\text{\hspace{0.05em}}\stackrel{\xaf}{r}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=\frac{{C}_{1}}{4\pi {a}_{0}^{2$ (27)

Here ${C}_{1}$ is constant and is determined after obtaining the corresponding solution ${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)$ to the eigenvalue.

To obtain the Tayler expansion, we define these integral as

${F}_{1}\left(\stackrel{\xaf}{r}\right)={\displaystyle {\int}_{0}^{\stackrel{\xaf}{r}}\text{d}\stackrel{\xaf}{{r}^{\prime}}}\text{\hspace{0.05em}}{\stackrel{\xaf}{{r}^{\prime}}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{{r}^{\prime}}\right)$ (28)

${F}_{2}\left(\stackrel{\xaf}{r}\right)=-{\displaystyle {\int}_{0}^{\stackrel{\xaf}{r}}\text{d}\stackrel{\xaf}{{r}^{\prime}}}\text{\hspace{0.05em}}\stackrel{\xaf}{{r}^{\prime}}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{{r}^{\prime}}\right)$ (29)

Then Tayler expansion around $\stackrel{\xaf}{r}=0$ of ${F}_{1}\left(\stackrel{\xaf}{r}\right)$ and ${F}_{2}\left(\stackrel{\xaf}{r}\right)$ become

${F}_{1}\left(\stackrel{\xaf}{r}\right)=\stackrel{\xaf}{r}\left({\stackrel{\xaf}{r}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\right)+\frac{1}{2}{\stackrel{\xaf}{r}}^{2}\left(2\stackrel{\xaf}{r}{\rho}_{0}\left(\stackrel{\xaf}{r}\right)+{\stackrel{\xaf}{r}}^{2}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}\right)$ (30)

${F}_{2}\left(\stackrel{\xaf}{r}\right)=-\stackrel{\xaf}{r}\left(\stackrel{\xaf}{r}{\stackrel{\xaf}{\rho}}_{0}\left(r\right)\right)-\frac{1}{2}{\stackrel{\xaf}{r}}^{2}\left({\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)+\stackrel{\xaf}{r}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}\right)$ (31)

Note that usual Tayler expansion around
$\stackrel{\xaf}{r}=0$ involves using the description of the value of *n*th derivative at
$\stackrel{\xaf}{r}=0$ times
$\frac{{\stackrel{\xaf}{r}}^{n}}{n!}$. However, we use only the nth derivative term itself because we are considering the case that
$\stackrel{\xaf}{r}$ itself is very close to zero. And note that this consideration restricts the order of singularity at most 1, which means that if
${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)$ has
${\stackrel{\xaf}{r}}^{-n}$ term largest *n* is 1.

For the time being, we neglect derivative expression part in Equation (30) and Equation (31). Then we obtain

$\frac{1}{\stackrel{\xaf}{r}}{F}_{1}\left(\stackrel{\xaf}{r}\right)+{F}_{2}\left(\stackrel{\xaf}{r}\right)=\frac{1}{2}{\stackrel{\xaf}{r}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)$ (32)

We will demonstrate that this neglecting process does not affect the result by checking the corresponding eigenvalue solution ${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)$.

Then Equation (25) becomes

$\begin{array}{c}{\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=-\frac{{\partial}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial {\stackrel{\xaf}{r}}^{2}}-\frac{2}{\stackrel{\xaf}{r}}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}-\frac{{P}_{0}}{{\left(\frac{\hslash c}{{a}_{0}}\right)}^{2}}\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}\left(1-{a}_{0}{C}_{1}\right){\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{{P}_{0}}{{\left(\frac{\hslash c}{{a}_{0}}\right)}^{2}}\frac{{e}^{2}}{4\pi {\u03f5}_{0}{a}_{0}}\frac{1}{2}{\stackrel{\xaf}{r}}^{2}{\left({\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\right)}^{2}\end{array}$ (33)

Equation (33) is a homogeneous nonlinear second order differential equation. Because we have not discovered an exact treatment for this equation, we set the following condition to linearize it.

${\left({\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\right)}^{2}={\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)+f\left(\stackrel{\xaf}{r}\right)$ (34)

Then Equation (33) becomes the following nonhomogeneous linear second order differential equation.

$\begin{array}{l}\frac{{\partial}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial {\stackrel{\xaf}{r}}^{2}}+\frac{2}{\stackrel{\xaf}{r}}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}+\left({\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}+\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)\left(1-{a}_{0}{C}_{1}\right)\right){\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)-\frac{1}{2}\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right){\stackrel{\xaf}{r}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\\ =\frac{1}{2}\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)f\left(\stackrel{\xaf}{r}\right)\end{array}$ (35)

To obtain the Equation (35), we used the fact that $\frac{{e}^{2}}{4\pi {\u03f5}_{0}}=\alpha \hslash c$.

To solve Equation (35), we need to solve the following homogeneous equation.

$\begin{array}{l}\frac{{\partial}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial {\stackrel{\xaf}{r}}^{2}}+\frac{2}{\stackrel{\xaf}{r}}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}+\left({\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}+\left(\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)\left(1-{a}_{0}{C}_{1}\right)\right)\right){\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\\ -\frac{1}{2}\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right){\stackrel{\xaf}{r}}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=0\end{array}$ (36)

For the simplicity we use the following notation

$T={\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}+\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)\left(1-{a}_{0}{C}_{1}\right)$ (37)

$S=\frac{1}{2}\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)$ (38)

Then Equation (36) is expressed as

$\frac{{\partial}^{2}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial {\stackrel{\xaf}{r}}^{2}}+\frac{2}{\stackrel{\xaf}{r}}\frac{\partial {\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial \stackrel{\xaf}{r}}+\left(T-S{\stackrel{\xaf}{r}}^{2}\right){\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=0$ (39)

To solve Equation (39), we use ${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=\frac{1}{\stackrel{\xaf}{r}}{\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}\left(\stackrel{\xaf}{r}\right)$ and rearrange Equation (39) to become the following after factor out $\frac{1}{\stackrel{\xaf}{r}}$

$\frac{{\partial}^{2}{\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}\left(\stackrel{\xaf}{r}\right)}{\partial {\stackrel{\xaf}{r}}^{2}}+\left(T-S{\stackrel{\xaf}{r}}^{2}\right){\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}\left(\stackrel{\xaf}{r}\right)=0$ (40)

To determine the specific eigenvalue, we use $z={\stackrel{\xaf}{r}}^{2}$.

Then Equation (40) becomes

$\frac{{\partial}^{2}{\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}}{\partial {z}^{2}}+\frac{1}{2z}\frac{\partial {\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}}{\partial z}+\left[\frac{T}{4z}-\frac{S}{4}\right]{\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}=0$ (41)

By setting as ${\stackrel{\xaf}{\stackrel{\xaf}{\rho}}}_{0}\left(z\right)={z}^{-\frac{1}{4}}{\stackrel{^}{\rho}}_{0}\left(z\right)$ and factor out ${z}^{-\frac{1}{4}}$, Equation (41) becomes

$\frac{{\partial}^{2}{\stackrel{^}{\rho}}_{0}}{\partial {z}^{2}}+\left(-\frac{S}{4}+\frac{T}{4z}+\frac{\frac{3}{16}}{{z}^{2}}\right){\stackrel{^}{\rho}}_{0}=0$ (42)

To obtain the standard form of Whittaker equation, we use $z=\xi \stackrel{\xaf}{z}$.

Then Equation (42) becomes

$\frac{{\partial}^{2}{\stackrel{^}{\rho}}_{0}}{\partial {\stackrel{\xaf}{z}}^{2}}+\left(-\frac{S{\xi}^{2}}{4}+\frac{\frac{1}{4}T\xi}{\stackrel{\xaf}{z}}+\frac{\frac{3}{16}}{{\stackrel{\xaf}{z}}^{2}}\right){\stackrel{^}{\rho}}_{0}=0$ (43)

By setting the value of $\xi $ at $S{\xi}^{2}=1$, Equation (43) becomes

$\frac{{\partial}^{2}{\stackrel{^}{\rho}}_{0}}{\partial {\stackrel{\xaf}{z}}^{2}}+\left(-\frac{1}{4}+\frac{\frac{1}{4}T{S}^{-\frac{1}{2}}}{\stackrel{\xaf}{z}}+\frac{\frac{3}{16}}{{\stackrel{\xaf}{z}}^{2}}\right){\stackrel{^}{\rho}}_{0}=0$ (44)

The standard form of the Whittaker equation is presented in Ref. [9] as

$\frac{{\partial}^{2}W}{\partial {\stackrel{\xaf}{z}}^{2}}+\left(-\frac{1}{4}+\frac{\kappa}{\stackrel{\xaf}{z}}+\frac{-{\mu}^{2}+\frac{1}{4}}{{\stackrel{\xaf}{z}}^{2}}\right)W=0$ (45)

By comparison to the standard form of Whittaker equation, we obtain

$\begin{array}{l}\mu =\pm \frac{1}{4}\\ \kappa =\frac{1}{4}T{S}^{-\frac{1}{2}}\end{array}$

Basic solutions of Equation (45) are

${M}_{\kappa ,\mu}\left(\stackrel{\xaf}{z}\right)$ and ${M}_{\kappa ,-\mu}(\; z\; \xaf\; )$

The definition of ${M}_{\kappa ,\mu}\left(\stackrel{\xaf}{z}\right)$ is given in Ref [9] as

$\begin{array}{c}{M}_{\kappa ,\mu}\left(\stackrel{\xaf}{z}\right)={\stackrel{\xaf}{z}}^{\mu +\frac{1}{2}}\mathrm{exp}\left(-\frac{\stackrel{\xaf}{z}}{2}\right)F\left(\mu -\kappa +\frac{1}{2},2\mu +1,\stackrel{\xaf}{z}\right)\\ ={\stackrel{\xaf}{z}}^{\mu +\frac{1}{2}}\mathrm{exp}\left(-\frac{\stackrel{\xaf}{z}}{2}\right)\underset{n=0}{\overset{\infty}{{\displaystyle \sum}}}\frac{\Gamma \left(2\mu +1\right)\Gamma \left(\mu -\kappa +n+\frac{1}{2}\right)}{\Gamma \left(2\mu +n+1\right)\Gamma \left(\mu -\kappa +\frac{1}{2}\right)}\frac{{\stackrel{\xaf}{z}}^{n}}{n!}\end{array}$ (46)

Because we do not want continuous eigenvalue, the series must be terminated. This condition gives

$\mu -\kappa +n+\frac{1}{2}=0$ (47)

The Lowest value of corresponds to lowest eigenvalue comes from the case of *n *= 0 and
$\mu =-\frac{1}{4}$.

Then we obtain the determining equation for eigenvalue as

$\kappa =\frac{1}{4}=\frac{1}{4}T{S}^{-\frac{1}{2}}=\frac{1}{4}\left({\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}+\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)\left(1-{a}_{0}{C}_{1}\right)\right){\left(\frac{2}{\alpha}\right)}^{\frac{1}{2}}{\left(\frac{\frac{\hslash c}{{a}_{0}}}{{P}_{0}}\right)}^{\frac{1}{2}}$ (48)

Rearrangement yields the following for the eigenvalue determining equation.

$\left({\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{2}+\alpha \left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)\left(1-{a}_{0}{C}_{1}\right)\right){\left(\frac{2}{\alpha}\right)}^{\frac{1}{2}}{\left(\frac{\frac{\hslash c}{{a}_{0}}}{{P}_{0}}\right)}^{\frac{1}{2}}=1$ (49)

Recalling that $\frac{1}{\stackrel{\xaf}{r}}$ and ${\stackrel{\xaf}{z}}^{-\frac{1}{4}}$, as well as $\stackrel{\xaf}{z}=\frac{z}{\xi}$ and $z={\stackrel{\xaf}{r}}^{2}$, can be factored out, the corresponding eigenfunction becomes

${\stackrel{\xaf}{\rho}}_{0}^{\left(1\right)}\left(\stackrel{\xaf}{r}\right)=const\frac{1}{\stackrel{\xaf}{r}}\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)$ (50)

This solution has maximum, actually it is a singularity, at $\stackrel{\xaf}{r}=0$, so it has an appropriate solution from Tayler expansion around $\stackrel{\xaf}{r}=0$. The condition for the order of singularity that we mentioned before is also satisfied. And we can

obtain the same factor of $\frac{1}{2}$ when calculating the second order Tayler series with the derivative term that was neglected before. Thus, this solution is very acceptable.

Another basic solution is

${\stackrel{\xaf}{\rho}}_{0}^{\left(2\right)}=const\frac{{\xi}^{\frac{1}{2}}}{{\stackrel{\xaf}{r}}^{\frac{3}{2}}}{M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)$ (51)

Note that ${\stackrel{\xaf}{\rho}}_{0}^{\left(2\right)}$ does not have a singularity at $\stackrel{\xaf}{r}=0$. We need this solution also because in Equation (35) we are dealing with nonhomogeneous equation and we must construct a particular solution to complete it.

We adopt
$\mu =-\frac{1}{4}$ solution because in the case the solution corresponds to eigenvalue. If we take
$\mu =+\frac{1}{4}$ in Equation (47), the determining equation is
$\kappa =\frac{3}{4}$ (*n *= 0) instead
$\kappa =\frac{1}{4}$, which has no singularity at
$\stackrel{\xaf}{r}=0$. Therefore the type of solution in Equation (50) is more appropriate for the consideration besides the determined eigenvalue is larger than
$\mu =-\frac{1}{4}$ case.

To obtain $f\left(\stackrel{\xaf}{r}\right)$, we describe ${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)$ as a function of $f\left(\stackrel{\xaf}{r}\right)$ by solving the following equation

${\left({\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)\right)}^{2}={\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)+f\left(\stackrel{\xaf}{r}\right)$ (52)

We insert the solution of Equation (52), that is a function of $f\left(\stackrel{\xaf}{r}\right)$, into Equation (33) and obtain the following equation.

$\begin{array}{l}\frac{{\partial}^{2}f}{\partial {\stackrel{\xaf}{r}}^{2}}-\frac{2{\left(\frac{\partial f}{\partial \stackrel{\xaf}{r}}\right)}^{2}}{1+4f}+\frac{2}{\stackrel{\xaf}{r}}\frac{\partial f}{\partial \stackrel{\xaf}{r}}+T\left(\frac{1\pm \sqrt{1+4f}}{2}\right)\frac{\sqrt{1+4f}}{2}\\ \text{\hspace{0.05em}}-S{\stackrel{\xaf}{r}}^{2}\left(\frac{1\pm \sqrt{1+4f}}{2}+f\right)\frac{\sqrt{1+4f}}{2}=0\end{array}$ (53)

We could not find the exact solution of Equation (54), so instead we seek a plausible solution by recalling that the equation we are dealing with is under the condition that $\stackrel{\xaf}{r}$ is very small (near 0). Then as a trial solution, we set the following form.

$1+4f\left(\stackrel{\xaf}{r}\right)=a{\stackrel{\xaf}{r}}^{2}+b\stackrel{\xaf}{r}+c$ (54)

Recalling that $\stackrel{\xaf}{r}$ is small, and noticing the fact that

$\frac{2}{\stackrel{\xaf}{r}}\frac{\partial f}{\partial \stackrel{\xaf}{r}}=2a+\frac{\frac{b}{4}}{\stackrel{\xaf}{r}}$

the $\frac{1}{\stackrel{\xaf}{r}}$ term does not appear separately from the other terms, thus must be 0.

Then, we can use the following approximations to describe Equation (53) second term of $1+4f=c$

$\sqrt{1+4f}=\sqrt{c}{\left(1+\frac{a{\stackrel{\xaf}{r}}^{2}}{\sqrt{c}}\right)}^{\frac{1}{2}}=\sqrt{c}\left(1+\frac{a{\stackrel{\xaf}{r}}^{2}}{2\sqrt{c}}\right)$

Comparing each ${\stackrel{\xaf}{r}}^{0}$ (order o term), ${\stackrel{\xaf}{r}}^{2}$ term, we obtain the following equations.

${\stackrel{\xaf}{r}}^{0}$ term $0=\pm \left(\frac{3a}{2}\right)+\frac{\sqrt{c}\pm c}{2}T$ (55)

${\stackrel{\xaf}{r}}^{2}$ term $0=\pm \left(-\frac{{a}^{2}}{2c}\right)+\frac{1}{2}\left(\frac{1}{\sqrt{c}}\pm 1\right){a}^{2}T-S\left(\frac{c\pm c\sqrt{c}}{2}+\frac{c\sqrt{c}-\sqrt{c}}{4}\right)$ (56)

The remaining terms are only *S* multiplies of order
${\stackrel{\xaf}{r}}^{4}$ and
${\stackrel{\xaf}{r}}^{6}$ terms, and because we are dealing with Equation (35) under the condition that
$\stackrel{\xaf}{r}$ is small, we neglect them.

For the + sign, we obtain the determining equation for *c* as

$\begin{array}{l}0=\left(\frac{3}{4}S+\frac{T}{2}\right){\left(\sqrt{c}\right)}^{4}+\left(\frac{S}{2}+\frac{3}{2}T\right){\left(\sqrt{c}\right)}^{3}\\ \text{}+\left(\frac{S}{4}-\frac{{T}^{2}}{9}+T\right){\left(\sqrt{c}\right)}^{2}-\frac{{T}^{2}}{18}\sqrt{c}-\frac{{T}^{2}}{18}\end{array}$ (57)

Once the value of *c* is obtained, a is given from Equation (55) as

$a=-\frac{\sqrt{c}+c}{3}$ (58)

Here, *T* and *S* are defined by Equation (37) and Equation (38).

Note that the obtained $f\left(r\right)$ is only an approximation and includes uncertainty. However, this uncertainty does not affect the eigenvalue because it is only involved in the ambiguity for ${a}_{0}$, as shown in later.

Then our $f\left(\stackrel{\xaf}{r}\right)$ is expressed as

$f\left(\stackrel{\xaf}{r}\right)=\stackrel{\xaf}{a}{\stackrel{\xaf}{r}}^{2}+\stackrel{\xaf}{c}$ (59)

where $\stackrel{\xaf}{a},\stackrel{\xaf}{c}$ denote the values by Equation (58) and Equation (57).

To obtain a particular solution of Equation (35), we use the Wronskian method. Recalling that Equation (36) is homogeneous and that its solutions are obtained in Equation (50) and Equation (51), the Wronskian becomes

$W\left(\stackrel{\xaf}{r}\right)=\left|\begin{array}{cc}\frac{1}{\stackrel{\xaf}{r}}{\text{e}}^{-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}}& \frac{\partial}{\partial \stackrel{\xaf}{r}}\left(\frac{1}{\stackrel{\xaf}{r}}{\text{e}}^{-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}}\right)\\ \frac{{\xi}^{\frac{1}{2}}}{{\stackrel{\xaf}{r}}^{\frac{3}{2}}}{M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)& \frac{\partial}{\partial \stackrel{\xaf}{r}}\left(\frac{{\xi}^{\frac{1}{2}}}{{\stackrel{\xaf}{r}}^{\frac{3}{2}}}{M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)\right)\end{array}\right|$ (60)

Then a particular solution of Equation (35) is constructed as

${\stackrel{\xaf}{\rho}}_{0}^{P}\left(\stackrel{\xaf}{r}\right)=\frac{{\xi}^{\frac{1}{2}}}{{\stackrel{\xaf}{r}}^{\frac{3}{2}}}{M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right){\displaystyle {\int}_{0}^{\stackrel{\xaf}{r}}\text{d}{r}^{\prime}}\frac{f\left({r}^{\prime}\right)\frac{1}{{r}^{\prime}}{\text{e}}^{-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}}}{W\left({r}^{\prime}\right)}-\frac{1}{\stackrel{\xaf}{r}}{\text{e}}^{-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}}{\displaystyle {\int}_{0}^{\stackrel{\xaf}{r}}\text{d}{r}^{\prime}}\frac{f\left({r}^{\prime}\right)\frac{{\xi}^{\frac{1}{2}}}{{{r}^{\prime}}^{\frac{3}{2}}}{M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)}{W\left({r}^{\prime}\right)}$ (61)

This formula is presented as an example in book of Ince [10].

Thus, the general solution of Equation (35) is written as

${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)={\stackrel{\xaf}{\rho}}_{0}^{P}\left(\stackrel{\xaf}{r}\right)+A\frac{1}{\stackrel{\xaf}{r}}\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)+B\left(\frac{{\xi}^{\frac{1}{2}}}{{\stackrel{\xaf}{r}}^{\frac{3}{2}}}\right){M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)$ (62)

where *A* and *B* are arbitrary constants determined by condition Equation (26) and Equation (27).

The corresponding eigenvalue is given in Equation (49). Equation (49) is rewritten as

${\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{\frac{3}{2}}{\left(\frac{2}{\alpha}\right)}^{\frac{1}{2}}+{\left(2\alpha \right)}^{\frac{1}{2}}\left(1-{a}_{0}{C}_{1}\right){\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{\frac{1}{2}}=1$ (63)

To see whether this scheme works or not, we use the following form for ${M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)$.

${M}_{\frac{1}{4},\frac{1}{4}}\left(\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}\right)={\left(\frac{\stackrel{\xaf}{r}}{\xi}\right)}^{\frac{3}{2}}\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)$ (64)

Then the second basic solution becomes

${\stackrel{\xaf}{\rho}}_{0}^{\left(2\right)}\left(\stackrel{\xaf}{r}\right)=\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)$ (65)

Then Wronskian becomes $W=\frac{1}{{\stackrel{\xaf}{r}}^{2}}{\text{e}}^{-\frac{{\stackrel{\xaf}{r}}^{2}}{\xi}}$. Thus we obtain the particular solution

${\stackrel{\xaf}{\rho}}_{P}\left(\stackrel{\xaf}{r}\right)=\left[c+\frac{5}{2}\left|a\right|\sqrt{\frac{2}{\alpha}}\right]{\stackrel{\xaf}{r}}^{2}$ (66)

Then the general solution of Equation (35) is expressed as

${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)={\stackrel{\xaf}{f}}_{P}+A\left(\frac{1}{\stackrel{\xaf}{r}}\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)\right)+B\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)$ (67)

Inserting Equation (67) to the condition equations Equation (26) and Equation (27), we obtain arbitrary constants *A* and *B* as

$A=6\left(\frac{1}{4\pi {a}_{0}^{3}}-\frac{\eta}{5}+\frac{\eta}{3}{\text{e}}^{\frac{1}{2\xi}}\right)$ (68)

$B=-3\eta {\text{e}}^{\frac{1}{2\xi}}+\frac{6\eta}{5}-\frac{6}{4\pi {a}_{0}^{3}}$ (69)

where $\eta =c+\frac{5}{2}\left|a\right|\sqrt{\frac{2}{\alpha}}$. Then, the solution of this system is

$\begin{array}{c}{\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)=\eta {\stackrel{\xaf}{r}}^{2}+\frac{6}{4\pi {a}_{0}^{3}}\left(\frac{1}{\stackrel{\xaf}{r}}-1\right)\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)+\left(-\frac{6\eta}{5}+2\eta {\text{e}}^{\frac{1}{2\xi}}\right)\frac{1}{\stackrel{\xaf}{r}}\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left(-3\eta {\text{e}}^{\frac{1}{2\xi}}+\frac{6\eta}{5}\right)\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(\stackrel{\xaf}{r}\in \left[0,1\right]\right)\end{array}$ (70)

To calculate ${C}_{1}$, we use only the second term because it is proportional to $\frac{1}{{a}_{0}^{3}}$ and reminding the fact that ${a}_{0}$ is very small quantity.

Then we obtain ${a}_{0}{C}_{1}=3$. Note that the second term is not dependent on $\eta $, so we can calculate anal most exact eigenvalue without knowing its value.

Then Equation (49) (the determining equation for eigenvalue) becomes as

${\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{\frac{3}{2}}-2\alpha {\left(\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}\right)}^{\frac{1}{2}}=\sqrt{\frac{\alpha}{2}}$ (71)

The solution of Equation (71) is

$\frac{{P}_{0}}{\frac{\hslash c}{{a}_{0}}}={4.74}^{2}\alpha $ (72)

Thus ${P}_{0}$ is expressed as

${P}_{0}=\frac{\alpha \hslash c}{{a}_{c}}\frac{{4.74}^{2}{a}_{c}}{{a}_{0}}$

where ${a}_{c}$ denotes classical electron radius.

Thus if ${a}_{0}={4.74}^{2}{a}_{c}$, then must be ${P}_{0}={m}_{e}{c}^{2}$ (electron rest mass energy).

Using Equation (72) we can determine $\sqrt{c}$ to be 0.02284. With Equation (37), we can find $c=0.000522$ and with Equation (38), we can find $a=-0.00779$. From these values $\eta $ becomes as $\eta =0.322822$.

We describe this scheme in terms of choosing only the first term of ${M}_{\frac{1}{4},\frac{1}{4}}$, but we can also obtain ${P}_{0}={m}_{e}{c}^{2}$ when using the exact form of ${M}_{\frac{1}{4},\frac{1}{4}}$. This is because we can absorb the difference caused by including infinite series terms of ${M}_{\frac{1}{4},\frac{1}{4}}$ into ${a}_{0}$ used to obtain the electron mass energy eigenvalue. Our viewpoint is that the result using exact form of ${M}_{\frac{1}{4},\frac{1}{4}}$ is not very different from the result shown in this paper because the maximum value of the nth term of ${M}_{\frac{1}{4},\frac{1}{4}}$ is

$\frac{1}{\left(2n+1\right)n!}{\left(\frac{\alpha}{2}\right)}^{\frac{n}{2}}\approx \frac{1}{\left(2n+1\right)n!}{\left(\frac{1}{17}\right)}^{n}$

3. Results

Using the equation of motion for a segment of charge, we obtain the charge wave function for an electron as Equation (63). We call ${\stackrel{\xaf}{\rho}}_{0}\left(\stackrel{\xaf}{r}\right)$ as a wave function because it appears as the eigenfunction corresponding to eigenvalue. An important point is that this is interpreted as a charge distribution function. To make this clearer, we performed a calculation example by choosing only the first term of the infinite series of ${M}_{\frac{1}{4},\frac{1}{4}}$. This result in an obtained charge distribution function that comes from an eigenfunction of electron mass energy. Also, its dominant function is $\frac{1}{\stackrel{\xaf}{r}}\mathrm{exp}\left(-\frac{{\stackrel{\xaf}{r}}^{2}}{2\xi}\right)$, which can be favorably compared to Weinberg’s charge distribution function relating $\eta \left(\stackrel{\to}{r}\right)$ to $-\frac{const}{{r}^{3}}{\displaystyle {\int}_{0}^{1}\text{d}x}\text{\hspace{0.05em}}x\left(1-x\right)\left(1+\frac{mr}{\sqrt{x\left(1-x\right)}}\right)\mathrm{exp}\left(-\frac{mr}{\sqrt{x\left(1-x\right)}}\right)$. We can say that the characteristic form of this expression is $-\frac{1}{{r}^{3}}\mathrm{exp}\left(-mr\right)$, which has a minus sign that, as pointed out by Weinberg, becomes plus sign after volume integration. This is exactly the same as our ${\stackrel{\xaf}{\rho}}_{0}\left(r\right)$ expression. Weinberg adds the $\delta \left(\stackrel{\to}{r}\right)$ term to satisfy this condition. This modification seems to be inconsistent because it creates dipole moment, which is unlikely for reasons described later. However, apart from this and sign issue, we think that the behavior of Weinberg expression as a charge distribution function carries real physical meaning. When we compare our ${\stackrel{\xaf}{\rho}}_{0}\left(r\right)$ to the characteristic form of Weinberg’s $\eta \left(\stackrel{\to}{r}\right)$, both functions have a singularity at the origin, although the order of them is different and the corresponding term of the exponential part is Gaussian for our case. Thus we can say that the overall behavior is similar. The most interesting point is that our derivation is totally different from that of Weinberg even though both results are similar. Thus, we can say that our result may reflect Weinberg’s comment that by absorbing an infinite number of photons, an electron has mass while charge is simultaneously spread out to finite volume [3] because our result is obtained as an eigenfunction of electron mass energy.

Our result has another interesting property. It shows a charged electron radius, denoted by
${a}_{0}$, that is almost 23 times as much as the classical electron radius. This result is consistent with the result obtained by Kim *et al*. [11], who reformulate an electron equation of motion of electron by using nonrelativistic quantum mechanics. Anequation of motion of electron was first founded by Abraham [12] and Lorentz [13] more than century ago considering an electron as a point-like which caused run-aways solutions and no causal behavior. LaterCaldiola [14] extended it, still in classical theory, showed that no run-aways and causality solution occur when the electron charged radius is more than
$\frac{2}{3}$ times classical electron radius. Kim *et al*. used nonrelativistic quantum mechanics consideration based on neglecting magnetic term to this problem and obtained the result that are alistic electron with no run-aways and no causality must satisfy the condition that the charged radius is smaller than 1.75 times Compton wave lengths. However, it has also pointed out that nonrelativistic quantum mechanics cannot strictly valid in point-particle limit because of vacuum polarization. In addition, QED calculation in this matter has not been reported. This is presumably because no physicists know how to regularize soft photons. Even Weinberg showed his electron charged radius had some remaining divergence [3]. The ACME collaboration report that they could not find a dipole moment inside electron electric field suggests that regularization based on the idea that creation and annihilation of tiny plus charge and minus charge in electric field is unlikely [15]. Therefore, Weinberg’s expression of an electron charged radius is still meaningful. This suggests that the lower bound of an electron charged radius according to quantum theory remains unclear. Thus, at the present time, we can say safely that the realistic electron satisfy the condition that its charged radius is more than
$\frac{2}{3}$ times the classical electron radius and less than 1.75 times the Compton wave lengths. As we mentioned in sec. 2, our obtained
${a}_{0}$ is only approximation, but even using full expression of
${M}_{\frac{1}{4},\frac{1}{4}}$, would not yield a very different result. Thus, with respect to the electron charged radius, our result is also plausible.

4. Discussion

To this point, there has not been a decisive experiment result to determine the size of a charged electron, the few results we do have illustrate why this problem is difficult. Let us consider the example of
${e}^{+}-{e}^{-}$ collision. Approaches to using these data can be divided into two groups. One group says that an electron is a point-like because no resolution data show up even though collision energy is very high, very much larger than the mass energy of electron. The other group uses data for extrapolating to low energy, (
${q}^{2}\approx 0$ ), then claims that the electron charged radius is a finite and not point-like. The consideration that the size of a charged electron is obtained at low energy, (almost
${q}^{2}\approx 0$ ) during a
${e}^{+}-{e}^{-}$ collision experiment is based on Cabbibo’s paper [16]. Also, this consideration corresponds to the determination method of pion charged radius [17] together with the determination equation of
${\left|{F}_{\pi}\right|}^{2}$ from cross section of pion-electron collision data given by NA7 Collaboration (CERN) [18]. However, there is a paper that might suggest an answer to this problem. Karnieli *et al*. suggest that the temporal duration of Cherenkov radiation, envisioned for almost a century as a shock wave, is limited by underlying entanglement between particle and light defined current operator as
$\stackrel{^}{J}\left(\stackrel{\to}{r},t\right)=ec{\psi}^{\u2020}\left(\stackrel{\to}{r},t\right)\stackrel{\to}{\alpha}\psi \left(\stackrel{\to}{r},t\right)$ where
$\psi \left(\stackrel{\to}{r},t\right)$ is a delocalized wave function instead of being limited by
$\stackrel{\to}{J}\left(\stackrel{\to}{r},t\right)=e\stackrel{\to}{v}\delta \left(\stackrel{\to}{r}-\stackrel{\to}{v}t\right)$ which uses current density and emits shock wave [7]. This means that a free charged electron is not a point-like but spread out in space-time. However, a complete and accurate description of this property depends importantly on whether or not the Higgs mechanism generates the masses of leptons, or at least electrons.

Appendix

Considering polar coordinate each ${x}_{1},{x}_{2}$ and ${x}_{3}$ are represented as

${x}_{1}=r\mathrm{sin}\theta \mathrm{cos}\theta ,{x}_{2}=r\mathrm{sin}\theta \mathrm{sin}\theta ,{x}_{3}=\mathrm{cos}\theta $

$\begin{array}{c}-i{\alpha}^{k}{\partial}_{k}\left(1{q}_{0}\left(r\right)\right)=i\left({\alpha}^{1}\frac{\partial}{\partial {x}_{1}}+{\alpha}^{2}\frac{\partial}{\partial {x}_{2}}+{\alpha}^{3}\frac{\partial}{\partial {x}_{3}}\right)\left(1{q}_{0}\left(r\right)\right)\\ =-i\left({\alpha}^{1}\mathrm{sin}\theta \mathrm{cos}\theta \frac{\partial}{\partial r}+{\alpha}^{2}\mathrm{sin}\theta \mathrm{sin}\theta \frac{\partial}{\partial r}+{\alpha}^{3}\mathrm{cos}\theta \frac{\partial}{\partial r}\right)\left(1{q}_{0}\left(r\right)\right)\\ =-i\left({\alpha}^{1}\frac{{x}_{1}}{r}+{\alpha}^{2}\frac{{x}_{2}}{r}+{\alpha}^{3}\frac{{x}_{3}}{r}\right)\frac{\partial}{\partial r}\left(1{q}_{0}\left(r\right)\right)\\ =-i\left(\stackrel{\to}{\alpha}\cdot \stackrel{^}{r}\right)\frac{\partial}{\partial r}{q}_{0}\left(r\right)1\end{array}$

where 1 denote unit matrix.

Note that because
${q}_{\mu}$ is a function only of *r*, the
$\frac{\partial}{\partial \theta},\frac{\partial}{\partial \phi}$ terms become 0.
$-i{\alpha}^{k}{\partial}_{k}\left(-i\stackrel{\to}{\alpha}\cdot \stackrel{^}{r}\right){q}_{1}\left(r\right)$ case

Recalling that ${\left({\alpha}^{1}\right)}^{2}={\left({\alpha}^{2}\right)}^{2}={\left({\alpha}^{3}\right)}^{2}=1$ (case 1) ${\alpha}^{k}{\alpha}^{l}=-{\alpha}^{l}{\alpha}^{k}\left(k\ne l\right)$ (case 2), we obtain the following.

From case 1 term, multiplying by factor ${\left(-i\right)}^{2}=-1$ yields

$\begin{array}{l}\left(\frac{3}{r}-\frac{{x}_{1}^{2}+{x}_{2}^{2}+{x}_{3}^{2}}{{r}^{3}}\right){q}_{1}+\left({\left(\mathrm{sin}\theta \mathrm{cos}\phi \right)}^{2}+{\left(\mathrm{sin}\theta \mathrm{sin}\phi \right)}^{2}+{\left(\mathrm{cos}\theta \right)}^{2}\right)\frac{\partial}{\partial r}{q}_{1}\\ =\frac{2}{r}{q}_{1}+\frac{\partial}{\partial r}{q}_{1}\end{array}$

From case 2 term,

We show only the ${\alpha}^{1}{\alpha}^{2}+{\alpha}^{2}{\alpha}^{1}$ part. Recalling that ${\alpha}^{2}{\alpha}^{1}=-{\alpha}^{1}{\alpha}^{2}$

$\begin{array}{l}\frac{\partial}{\partial {x}_{1}}\left(\frac{{x}_{2}}{r}\right)+\frac{{x}_{2}}{r}\frac{\partial}{\partial {x}_{1}}-\frac{\partial}{\partial {x}_{2}}\left(\frac{{x}_{1}}{r}\right)-\frac{{x}_{1}}{r}\frac{\partial}{\partial {x}_{2}}\\ =-\frac{{x}_{2}{x}_{1}}{{r}^{3}}+\frac{{x}_{1}{x}_{2}}{{r}^{3}}+{\left(\mathrm{sin}\theta \right)}^{2}\mathrm{sin}\phi \mathrm{cos}\phi \frac{\partial}{\partial r}-{\left(\mathrm{sin}\theta \right)}^{2}\mathrm{cos}\phi \mathrm{sin}\phi \frac{\partial}{\partial r}=0\end{array}$

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

[1] |
Cho, A. (2012) Science, 338, 1524-1525. https://doi.org/10.1126/science.338.6114.1524 |

[2] |
Gonzalez, E., Kane, G., Nguyen, K.D. and Perry, M. (2021) Physical Review D, 103, Article ID: 126027. https://doi.org/10.1103/PhysRevD.103.126027 |

[3] | Weinberg, S. (1995) The Quantum Theory of Field I. Cambridge University Press, Cambridge. |

[4] |
Weinberg, S. (2020) Physical Review D, 101, Article ID: 035020. https://doi.org/10.1103/PhysRevD.101.035020 |

[5] |
Barr, S.M. and Zee, A. (1978) Physical Review D, 17, 1854. https://doi.org/10.1103/PhysRevD.17.1854 |

[6] |
Dobrescu, B.A. and Fox, P. (2008) Journal of High Energy Physics, 8, 100. https://doi.org/10.1088/1126-6708/2008/08/100 |

[7] |
Karnieli, A., Rivera, N. and Kaminer, I. (2021) Science Advances, 7. https://doi.org/10.1126/sciadv.abf8096 |

[8] | Weinberg, S. (2000) The Quantum Theory of Field III. Cambridge University Press, Cambridge. |

[9] | Moriguchi, S., Udagawa, K. and Hitotsumatsu, S. (1975) Mathematics Formula III. Iwanami. |

[10] | Ince, E.L. (1956) Ordinary Differential Equations. Dover Publication, Mineola. |

[11] | Kim, K.J. and Sessler, A.M. (1998) 8th Workshop on Advanced Acceleration Concepts July 5-11, 1998. |

[12] | Abraham, M. (1905) “Theorie der Elektrizitat” Vol. II: Elektromagnetische Theorie der Strahlung. Teubner, Leipzig. |

[13] | Lorentz, H.A. (1909) The Theory of Electrons. Teubner, Leipzig. |

[14] |
Caldirola, P. (1956) Nuovo Cimento, 3, 297. https://doi.org/10.1007/BF02743686 |

[15] |
Andreev, V., Ang, D.G., DeMille, D., et al. (2018) Nature, 562, 355-360. https://doi.org/10.1038/s41586-018-0599-8 |

[16] |
Cabbibo (1961) Physical Review, 124, 1577. https://doi.org/10.1103/PhysRev.124.1577 |

[17] |
Cui, Z.F. and Schmidt, S.M. (2021) Physics Letters B, 882, Article ID: 136631. https://doi.org/10.1016/j.physletb.2021.136631 |

[18] | Amendolia, S.R., et al. (1986) Nuclear Physics B, 277, 168-196. |

Journals Menu

Contact us

+1 323-425-8868 | |

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2024 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.