1. Introduction
The proton radius puzzle is one of the most fundamental challenges of modern physics. Before the year 2010, the proton charge radius r
p was determined by the spectroscopic method, relying on the electron energy levels in hydrogen atoms, and by the elastic scattering of electrons on protons. The mean value of the proton charge radius, recommended by CODATA (Committee on Data of the International Science Council), was r
p = (0.8775 ± 0.0051) × 10
−13 cm—see, e.g., the reviews by Pohl et al. [
1] and by Gao and Vanderhaenghen [
2], as well as references therein.
In 2010, Pohl et al. [
3], and then in 2013, Antognini et al. [
4] determined r
p from the experiment on muonic hydrogen atoms. Because the ratio of the muon mass m
μ to the electron mass m
e is m
μ/m
e ≈ 207, the average muon–proton distance in muonic hydrogen atoms is about 200 smaller than the electron–proton distance in electronic hydrogen atoms. Therefore, the shift in the energy of an S-state, caused by the finite proton size, for muonic hydrogen atoms is about 8 million times greater than for electronic hydrogen atoms. Consequently, muonic measurements should be much more sensitive to r
p than the corresponding electronic measurements. The resulting proton charge radius was claimed to be r
p = (0.84087 ± 0.00039) × 10
−13 cm, e.g., about 4% (or 5 standard deviations) smaller than the above CODATA value. This result prompted calls for a new physics model beyond the standard model.
In 2019, Bezginov et al. [
5] remeasured the n = 2 Lamb shift for electronic hydrogen atoms. They deduced the value of r
p consistent with the muonic measurements from papers [
3,
4]. In the same year, Xiong et al. [
6] remeasured r
p in the electron scattering experiment and found it to be consistent with the muonic measurements from [
3,
4].
The results from [
5,
6] favor the smaller charge radius of the proton. However, they do not explain why the experimental values of r
p, found before the year 2010, yielded the larger value. Besides, Fleurbaey et al. [
7] reported the larger value of r
p = (0.877 ± 0.013) × 10
−13 cm, obtained from the two-photon measurements in the electronic hydrogen (they measured the 1S–3S two-photon transition frequency of hydrogen using a continuous-wave excitation laser at 205 nm).
So, the puzzle is not considered to be resolved yet—see, e.g., the conclusions of Karr-Marchand of 2019 [
8] and of Gao-Vanderhaenghen’s review of 2022 [
2].
There are many theoretical factors contributing to the shift in S-states of muonic hydrogen atoms—see, e.g., reviews by Pohl et al. [
1] and by Karshenboim et al. [
9]. In the present paper, we bring to the attention of the research community, dealing with the proton radius puzzle, the contributing factor never taken into account in any previous calculations. This factor has to do with the hydrogen atoms of the second flavor, whose existence is confirmed in four different types of atomic experiments.
There are two analytical solutions of the Dirac equation for hydrogen atoms (two coupled differential equations for the components of the Dirac bispinor have two solutions). One solution is only weakly singular at small r, while the other solution is more strongly singular at small r. The second solution is rightly rejected for the model where the proton is point-like, as well as for the models where the charge distribution inside the proton is a uniform spherical shell or a uniformly charged sphere. However, well-known experiments on the elastic scattering of electrons on protons, performed in the previous century, revealed that the actual charge distribution has the maximum at r = 0, thus being significantly different from the above models (see, e.g., Simon et al. (1980) [
10] and Perkins (1987) [
11]).
In [
12,
13], the following was shown analytically. After taking into account the actual charge distribution inside the proton, the second solution outside the proton can be tailored with the regular solution inside the proton for any S-state. In other words, the second solution outside the proton is legitimate for all S-states. This second type of hydrogen atom possessing only the S-states (the energies of the S-states being the same as for the usual first solution) was later named the second flavor of hydrogen atoms (SFHA)—using an analogy with the quantum chromodynamics where up and down quarks are named two flavors [
14].
Outside the proton, for the S-states at small r, the radial wave function R(r) for the first solution scales is ~1/r
β/2, where
where α is the fine structure constant (α = e
2/(ħc) ≈ 0.007297), while for the SFHA, R(r) scales as ~ 1/r
2−β/2. Consequently, for relatively large values of the linear momentum p >> p
0 = me
2/ħ (where m is the mass of the atomic lepton, whether it is electron or muon), the corresponding wave function in the momentum representation φ(p) for the SFHA falls off much slower than for the hydrogen atoms of the first (usual) flavor. This is because φ(p) and R(r) are interconnected by the Fourier transform, so that, for the SFHA, the more rapid increase in R(r) as r decreases translates into the slower decrease in φ(p) as p increases in the range of p >> p
0.
By now, the existence of the SFHA is proven in four various types of atomic experiments, as follows.
For p
0 << p << mc, i.e., in the non-relativistic part of the tail of the distribution (we note that p
0/mc = α ≈ 0.007297), the experimental result, deduced by Gryzinski [
15] from the analysis of atomic experiments, was F
exper(p) ~ (mc/p)
4, while the corresponding theoretical result by Fock [
16] was F
theor(p) ~ (mc/p)
6. Here, F(p)dp is the probability of finding the linear momentum in the interval (p, p + dp). This means that, for the ratio F
theor(p)/F
exper(p) = (mc/p)
2, for the values of p ~ 10p
0, the discrepancy F
theor(p)/F
exper(p) between the experimental and theoretical results was ~200 times (!).
In [
12], it was shown that, with the allowance for the SFHA, this huge discrepancy was completely eliminated. No alternative explanation of this huge discrepancy was ever offered.
- B.
Experiments on the electron impact excitation of electronic hydrogen molecules
There was a discrepancy by at least a factor of two between the experimental and theoretical cross sections of the excitation to the lowest triplet states, as pointed out in [
17]. In the same paper, it was shown that this large discrepancy can be eliminated if the SFHA was present in the experimental gas. Again, no alternative explanation of this significant discrepancy was ever offered.
- C.
Experiments on the electron impact excitation of electronic hydrogen atoms
The theoretical ratio of the cross section for the excitation of the state 2s to the cross section for the excitation of the state 2p was 20% higher than the corresponding experimental ratio—well beyond the experimental error margin of 9%, as pointed out in [
18]. In the same paper, it was shown that this significant discrepancy can be eliminated if the SFHA was present in the experimental gas. Again, no alternative explanation of this significant discrepancy was ever offered.
- D.
Experiments on the charge exchange between electronic hydrogen atoms and protons
There was a noticeable discrepancy between the experimental and theoretical cross sections, as pointed out in [
19]. In the same paper, it was shown that this noticeable discrepancy can be eliminated if the SFHA was present in the experimental gas. Again, no alternative explanation of this significant discrepancy was ever provided.
The present paper has two central points. The first one is that muonic hydrogen atoms should also have two flavors—because all analytical results from [
12] for the ground state and their generalization in [
13] for any S-state are valid for muonic hydrogen atoms after replacing m
e in those calculations by m
μ. So, there should exist the second flavor of muonic hydrogen atoms (SFMHA).
The second central point of the present paper is that, because, for the SFMHA, the radial wave function R(r) in the vicinity of the proton—and consequently inside the proton (because both the outside and inside parts of R(r) match at the proton boundary)— is significantly different compared to the usual muonic hydrogen atoms, even a relatively small admixture of the SFMHA to the usual muonic hydrogen atoms in the experimental gas can affect the shift of the S-states, and thus modify the determination of the proton charge radius from the experimental Lamb shift of muonic hydrogen atoms.
We present a simple model illustrating that even a relatively small admixture of the SFMHA to the usual muonic hydrogen atoms in the experimental gas can lead to the false conclusion that the proton charge radius is about 4% smaller than its actual value.
2. Model
For the ground state of muonic hydrogen atoms, outside the proton, the radial part of the Dirac bispinor, based on Equation (17) from [
12], can be represented in the following form:
In Equation (2), ε is the relatively small share of the SFMHA in the experimental muonic hydrogen gas (ε << 1), R
p is the proton radius in units of the muonic Bohr radius a
0μ = ħ
2/(m
μe
2), and r is the distance from the origin in units of the muonic Bohr radius a
0μ. Equation (2) was simplified compared with Equation (17) from [
12], using the fact that β = α
2 << 1. We also note that, in Equation (17) from [
12], the second term in f(r) and g(r) was proportional to the quantity:
which is the shift (with the minus sign) in the ground state energy due to the finite proton size, with the shift being in units of m
μc
2. Because, in our Equation (2), the second term in f(r) and g(r) is assumed to be a relatively small correction to the first term (because ε << 1), while deriving Equation (2), we used for the shift the following approximate textbook expression (see, e.g., Flügge textbook [
20]):
The squared absolute value of the wave function of the ground state is
From Equation (2), it is seen that f
2(r)/g
2(r) ~ α
2 << 1, so that
The shift of the ground state energy δE due to the proton finite size is (in analogy to Equation (3) from [
1] or to Equation (66) from [
2])
where b is a constant of no importance for the purpose of the present paper. We would like to find out whether there exists a value of ε << 1, such that
so that, while disregarding a relatively small admixture of the SFMHA to the experimental muonic hydrogen gas, one would deduce—from the experimental shift—the value of R
p that would be 4% smaller than the actual value of R
p.
Equation (8) is quadratic with respect to ε—so, it has the following two solutions:
The numerical value of the proton charge radius r
p (defined as the root-mean-square radius of the proton charge distribution) in units of the muonic Bohr radius a
0μ is 0.00343. The proton “sphere” radius R
p would be a factor of (5/3)
1/2 greater than r
p (it would be equal to 0.00443) if the proton would be a uniformly charged sphere (which the proton is not). The actual value of R
p should be between 0.00343 and 0.00443. For further numerical estimates of ε
1 and ε
2, we adopt the value R
p ≈ 0.004, so that
Physically, the share of the SFMHA ε
2 = 0.13 seems to be slightly more preferable (compared with ε
1 = 0.003). This is because it is of the same order of magnitude as the share of the SFHA in the experimental gas of the electronic hydrogen molecules, which (the share) was required for eliminating the large discrepancy (by at least of a factor of two) between the theoretical and experimental cross sections of the excitation by the electron impact [
17].
As the proton charge radius rp is proportional to Rp, the above result about the determination of Rp from the energy shift is also true for rp. Namely, indeed, even a relatively small admixture of the SFMHA to the usual muonic hydrogen atoms in the experimental gas can lead to the false conclusion that the proton charge radius rp is about 4% smaller than its actual value.