Understanding the origin of quark confinement in hadrons remains one of the most challenging problems in modern physics. Recently, the pressure distribution inside the proton was measured via deeply virtual Compton scattering. Surprisingly, strong repulsive pressure up to 1035 pascals, the highest so far measured in our universe, was obtained near the center of the proton up to 0.6 fm, combined with strong binding energy at larger distances. We show here that this profile can be derived semiquantitatively without any adjustable parameters using the rotating lepton model of composite particles (RLM), i.e. a proton structure comprising a ring of three gravitationally attracting rotating ultrarelativistic quarks. The RLM synthesizes Newton’s gravitational law, Einstein’s special relativity, and the de Broglie’s wavelength expression, thereby conforming with quantum mechanics, and also yields a simple analytical formula for the proton radius and for the maximum measured pressure which are in excellent agreement with the experimental values.
The proton consists of fundamental particles called quarks and gluons. Gluons carry the force that binds quarks together. Quarks are always confined in the composite particles in which they are located. The origin of quark confinement is still a subject of intensive study in modern physics. Recently, for the first time, the pressure distribution inside the proton was measured via deeply virtual Compton scattering [1, 2, 3, 4]. Strong repulsive pressure up to 1035 pascals, the highest so far measured in our universe, was obtained near the center of the proton up to 0.6 fm, combined with strong binding energy at larger distances. This recent pioneering experimental [1, 4] and theoretical [1, 2, 3] work on the deeply virtual Compton scattering has opened a new area of research on the fundamental gravitational properties of protons, neutrons and nuclei and has underlined that gravity plays an important, perhaps dominant, role inside protons and other hadrons [1–3].
The dominant role of gravity inside hadrons has also emerged in recent years from the rotating lepton model of composite particles (RLM) [5, 6], in which gravity causes confinement of highly energetic neutrinos in bound rotational states and thus leads to formation of quarks, hadrons and bosons [5–8]. This model follows exactly the steps of the Bohr treatment of the H atom and contains no adjustable parameters (Figure 1). It synthesizes Newton’s gravitational law, Einstein’s special relativity , and de Broglie’s wavelength expression, thereby conforming with quantum mechanics. The model does not require any new theory. Furthermore, it fits extremely well with current experimental evidence for the masses and other properties of hadrons. The simple structure of the RLM for the proton is also suggested by the earlier proposed  bagel-shaped proton geometry (Figure 1). Without any adjustable parameters, the RLM predicts the masses of hadrons [5, 6], but also of bosons [7, 8] with an accuracy of one percent [5, 6, 7, 8, 11, 12].
Two important suggestions have emerged from the RLM analysis [5, 6]: First, quarks and hadrons consist primarily of rotating neutrinos, and their mass is due to the kinetic energy of the rotating neutrinos. Second, the Strong Force can be viewed as the relativistic gravitational force. Here the RLM is used to describe quantitatively the measured pressure distribution in protons and to derive simple analytical formulae for the proton radius and for the maximum measured pressure of 1035 pascals. We also show that a proton structure comprising a ring of three rotating ultrarelativistic quarks with radius 0.63 fm describes the proton’s measured pressure profile semi-quantitatively without any adjustable parameters.
Similar to the Bohr model of the H atom, but utilizing gravitational instead of electrostatic forces, the RLM for the proton comprises only two equations, i.e.
where re is the radius of the circle defined by the centers of the three rotating particles, which form an equilateral triangle, mo is the neutrino rest mass, and mg is the neutrino gravitational mass, together with the corresponding de Broglie wavelength equation
which is the historical basis of quantum mechanics and introduces quantization via the integer number n. Accounting for the equivalence principle, mg equals the inertial mass, mi, which as Einstein showed in his pioneering special relativity paper in 1905 , is equal, for linear motion, to the longitudinal mass, [9, 13]. Originally derived for linear motion , this useful result has been shown [5, 6], via the use of instantaneous reference frames , to remain valid for arbitrary motion including circular motion. Thus it follows
where is the Planck mass and the subscript “q” denotes quark.
Setting mp=938.272 MeV/c2, the proton mass, in equation (4), one computes mo = 0.0437 eV/c2, which remarkably is within the experimental limits of the mass of the heaviest electron neutrinos (0.048±0.01 eV/c2) [14, 15]. Consequently, the rest mass of quarks appears to be that of electron neutrinos. Also the radius, re, of the proton computed as the quark de Broglie wavelength from equation (2) for n = 1 and v ≈ c, is given by
in very good agreement with the experimental value [1, 16]. It is worth noting that the first equation (4), together with (3), dictate that the gravitational mass, of the relativistic rotational quarks is very close to the Planck mass, i.e.
Consequently the relativistic gravitational force and potential energy computed from eq (1) is
which are the values anticipated for the strong force between two quarks . These values are a factor of α−1 = 137.035 stronger than the Coulombic attraction and potential energy of an e+e− pair at the same distance.
The rotating quark ring structure of the proton according to the RLM allows for the quantitative description of the potential energy and pressure distribution inside a proton and their comparison with those measured in the pioneering reference . We first show in Figure 1, bottom, that the rotating quark ring of the RLM creates two types of attractive strong relativistic gravitational forces. First, those pointing inward and confining the three neutrinos to their rotational orbit at fm (Figure 1, bottom) and second, those pointing outwards and giving, due to their outward direction, the impression of repulsive forces (and pressure ) but in reality being equal and symmetric with the confining ones on the two sides of the rotating quark ring.
Thus, considering the distance of the quarks A and B fixed, and quark C at a variable location r on the perpendicular bisector of AB (Figure 2, bottom) one obtains, as shown in the methods section, that the total gravitational force, F, exerted by quarks A and B to quark C in the x direction is
which, in dimensionless form, is given by
Equation (8) is plotted in Figure 3, as F(r)/4π vs r and z, since this quantity is equal to r2p(r) (expressed in 10−2 GeV/fm) which has been used in the y axis of ref.  as shown in Figures 2 and 3. In this way, direct comparison of equation (8) with reference  is possible. The shaded part of Figure 3 is the one compared in Figure 2 with the experimental data of . There is remarkable semiquantitative agreement between reference  and equations (8) and (9).
It is also worth noting that the force expressions (8) and (9) describe in a simple way the two key characteristics of the Strong Force, i.e. asymptotic freedom (F = 0 at z+1 = 0, i.e. at fm) and confinement, i.e. dF/dr > 0 up to the force minimum at thus r=0.071 fm, in excellent agreement with the experimental results , as shown in Figures 2 and 3.
Utilizing F = −dU/dr, where U is the potential energy of quark C due to its interaction with A and B one obtains
which is plotted in Figure 4. The minimum U value occurs at z = −1, i.e. at and equals
while the corresponding value at the rotational center z = 0 and r = 0 fm is
for quark C and
for the entire proton, as shown in Figure 4. This is consistent with the virial theorem, i.e.
as shown in Figure 4, since, according to the RLM, the kinetic energy, Tp, of the relativistic neutrinos is the rest energy of the proton.
and the pressure, p, computed from p = F/4πr2 is of the order of
in very good agreement with reference . Indeed, as shown in Figures 2, 3 and 4, these values provide a good measure of the interquark forces and pressure inside the proton, in very good agreement with the experimental values measured in .
In summary, there is at least semiquantitative agreement between the RLM and the recent pioneering experimental results for force and pressure in the proton . This agreement provides strong support for the simple rotating quark model and demonstrates the excellent quality of the experimental data of references  and  and of their theoretical treatment . It also underlines the importance of gravitational forces inside the proton and suggests that, as previously proposed , the Strong Force can be viewed as relativistic gravity. This suggests that the simple methodology used here to model the proton structure and to compute the mass and basic properties of protons can be also applied to other hadrons as well.
Indeed the use of RLM has been explored recently to compute the masses of several hadrons [5, 6], but also bosons [7, 8], with an accuracy of one percent [5, 6, 7, 8, 11, 12], using similar simple rotating lepton structures, the components of which are chosen on the basis of the decay products of the corresponding composite particles. It thus appears that the present methodology may be of broader usefulness for the study of the structure and mass of composite particles.
Figure 5 shows the model geometry. Particles A, B and C, of gravitational mass mg each, lie on a circle of center M and radius R. The distance of points A and B is fixed and equal to Point C can move on the perpendicular bisector, x, of AB. When point C coincides with point C*, then the isoscele triangle ABC becomes equilateral, i.e. (AC*)=(BC*)=(AB) and point M coincides with point O. The distances OA, OB are always equal to Point C′ is the antidiametric of point C on the circle of radius, R, thus (CC′)=2R.
One notes in Figure 5 that
When C reaches C*, then and φ = 30∘, thus cos2 and equation (21) gives confirming its validity.
We denote FA and FR the gravitational forces exerted by particles A and B on particle C. It is
and we consider the centripetal force, F, acting on particle C. It is
Noting from (18) that
we can write equation (23) as
From Pythagoras’s theorem it follows
From (24) and (26) it follows
The force component in Figure 5 is given by
and its minimum occurs when
which leads to
it follows from (32) that
Consequently it follows
which implies that the minimum occurs at
and the maximum at
The recently measured for the first time internal pressure profile in a proton  is described semiquantitatively by the Rotating Lepton Model (RLM) of composite particles. There is excellent agreement between model and experiment regarding the radius (0.63 fm) and the maximum pressure (~2.5·1035 Pa).
This agreement provides very strong support for the validity of the rotating lepton model (RLM) and its great advantage over the standard model (SM) which neglects gravitational forces and does not allow for any such pressure computations except perhaps via inclusion of more adjustable parameters.
The RLM also explains that the strong repulsive pressure reported in  near the center of the proton is in reality due to the strong gravitational attraction by the rotating neutrino ring, which attracts matter equally strongly inside and outside its radius.
We are very thankful to Professor Donna Greenberg from Harvard University for numerous helpful discussions and style corrections.
 Burkert V.D., Elouadrhiri L., Girod F.X., The pressure distribution inside the proton, Nature, 2018, 557, 396. Search in Google Scholar
 Ji X.D., Deeply virtual Compton scattering, Phys. Rev. D, 1997, 55, 7114–7125. Search in Google Scholar
 Teryaev O.V., Gravitational form factors and nucleon spin structure, Front. Phys. 2016, 11, 111207. Search in Google Scholar
 Jo H.S. et al., Cross Sections for the Exclusive Photon Electro-production on the Proton and Generalized Parton Distributions CLAS Collaboration. Phys. Rev. Lett., 2015, 115, 212003. Search in Google Scholar
 Vayenas C.G., Souentie S., Gravity, special relativity and the strong force: A Bohr-Einstein-de-Broglie model for the formation of hadrons, (Springer, New York, 2012) Search in Google Scholar
 Vayenas C.G., Souentie S., & Fokas A., A Bohr-type model of a composite particle using gravity as the attractive force, Physica A, 2014, 405, 360. Search in Google Scholar
 Vayenas C.G., Fokas A.S., Grigoriou D., On the structure, masses and thermodynamics of the W± bosons, Physica A, 2016, 450, 37. Search in Google Scholar
 Fokas A.S., Vayenas C.G., Grigoriou D.P., On the mass and thermodynamics of the Higgs boson, Physica A, 2018, 492, 737. Search in Google Scholar
 Einstein A., Zür Elektrodynamik bewegter Körper. Ann. der Physik., 1905, 17, 891. Search in Google Scholar
 Miller G.A., Shapes of the proton, Phys. Rev. C, 2003, 68, 022201(R) Search in Google Scholar
 Vayenas C.G., Thermodynamics and catalysis of the generation of mass, Proc. Acad. of Athens, 2018, 93A, 97. Search in Google Scholar
 Vayenas C.G., Grigoriou D.P., Hadronization via gravitational confinement in Proc. of the 18th Lomonosov conference, in Particle Physics at the Year of 25th Anniversary of the Lomonosov Conferences, p. 517 (2019). Search in Google Scholar
 French A.P., Special relativity. (W.W. Norton and Co., New York, 1968). Search in Google Scholar
 Mohapatra R.N. et al., Theory of Neutrinos: A White paper, Rep. Prog. Phys, 2007, 70, 1757. Search in Google Scholar
 Aartsen M.G. et al., Measurement of Atmospheric Neutrino Oscillations at 6–56 GeV with IceCube DeepCore, Phys. Rev. Lett., 2018, 120, 071801. Search in Google Scholar
 Griflths D., Introduction to Elementary Particles. (2nd ed. Wiley-VCH Verlag GmbH & Co. KgaA, Weinheim, 2008). Search in Google Scholar
© 2019 C. G. Vayenas et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 Public License.