Open access peer-reviewed chapter

Exergy: Mechanical Nuclear Physics Measures Pressure, Viscosity and X-Ray Resonance in K-Shell in a Classical Way

Written By

Edward Henry Jimenez

Submitted: 09 September 2020 Reviewed: 07 December 2020 Published: 15 January 2021

DOI: 10.5772/intechopen.95405

From the Edited Volume

Nuclear Materials

Edited by Pavel V. Tsvetkov

Chapter metrics overview

589 Chapter Downloads

View Full Metrics


First, the liquid drop model assumes a priori; to the atomic nucleus composed of protons and neutrons, as an incompressible nuclear fluid that should comply with the Navier–Stokes 3D equations (N-S3D). Conjecture, not yet proven, however, this model has successfully predicted the binding energy of the nuclei. Second, the calculation of nuclear pressure p0∈1.42,1.94]1032Pa and average viscosity η=1.71×1024fm2/s, as a function of the nuclear decay constant k=p02η=1T1/2, not only complements the information from the National Nuclear Data Center, but also presents an analytical solution of (N- S3D). Third, the solution of (N-S3D) is a Fermi Dirac generalized probability function Pxyzt=11+ep02ηt−μx2+y2+z21/2, Fourth, the parameter μ has a correspondence with the Yukawa potential coefficient μ=αm=1/r, Fifth, using low energy X-rays we visualize and measure parameters of the nuclear surface (proton radio) giving rise to the femtoscope. Finally, we obtain that the pressure of the proton is 8.14 times greater than the pressure of the neutron, and 1000 times greater than the pressure of the atomic nucleus. Analyzed data were isotopes: 9≤Z≤92 and 9≤N≤200.


  • femtoscope
  • Navier Stokes 3D
  • nuclear viscosity
  • minimum entropy

1. Introduction

Neutron stars are among the densest known objects in the universe, withstanding pressures of the order of 1034Pa. However, it turns out that protons [1], the fundamental particles that make up most of the visible matter in the universe, contain pressures 10 times greater, [2, 3] 1035Pa. This has been verified from two perspectives at the Jefferson Laboratory, MIT [1, 2, 3, 4, 5]. High-energy physics continue to guide the study of the mechanical properties of the subatomic world.

Viscosity is a characteristic physical property of all fluids, which emerges from collisions between fluid particles moving at different speeds, causing resistance to their movement (Figure 1). When a fluid is forced to move by a closed surface, similar to the atomic nucleus, the particles that make up the fluid move slower in the center and faster on the walls of the sphere. Therefore, a shear stress (such as a pressure difference) is necessary to overcome the friction resistance between the layers of the nuclear fluid. For the same radial velocity profile, the required tension is proportional to the viscosity of the nuclear fluid or its composition given by ZN [6, 7, 8].

Figure 1.

Obtaining nuclear viscosity and nuclear pressure from the speed of the neutron particles in the disintegration of chemical element. Figure 1a. indicates that the BE/A ratio is proportional to the nuclear pressure, p0 represented in Figure 1b. Figure 1c is the graph of the viscosity in equilibrium and out of equilibrium, the viscosity in equilibrium is greater than the viscosity at the moment of nuclear decay. Figure 1d is the average half-time of each isotope, 9Z92 and 9N200.

Radioactive decay is a stochastic process, at the level of individual atoms. According to quantum theory, it is impossible to predict when a particular atom will decay, regardless of how long the atom has existed. However, for a collection of atoms of the same type [8], the expected decay rate is characterized in terms of its decay constant k=1T1/2. The half-lives of radioactive atoms do not have a known upper limit, since it covers a time range of more than 55 orders of magnitude, from almost instantaneous to much longer than the age of the universe [8, 9].

Other characteristics of the proton such as its size have been studied in many institutes such as Max Plank, where it has been measured with high precision ranges rp=0.8418467fm, providing new research methods [10, 11].

The atomic nucleus is an incompressible fluid, justified by the formula of the nuclear radius, R=1.2A13, where it is evident that the volume of the atomic nucleus changes linearly with A=Z+N, giving a density constant [11]. All incompressible fluid and especially the atomic nucleus comply with the Navier Stokes equations. We present a rigorous demonstration on the incomprehensibility of the atomic nucleus, which allows to write explicitly the form of the nuclear force FN=gμ28πA1P1Pr, which facilitates the understanding of nuclear decay.

The Navier Stokes equations are a problem of the millennium [12, 13], that has not been resolved yet in a generalized manner. We present a solution that logically meets all the requirements established by the Clay Foundation [14, 15]. This solution coherently explains the incompressible nuclear fluid and allows calculations of the nuclear viscosity and nuclear pressure [1, 2].

The alpha particle is one of the most stable. Therefore it is believed that it can exist as such in the heavy core structure. The kinetic energy typical of the alpha particles resulting from the decay is in the order of 5 MeV.

For our demonstrations, we will use strictly the scheme presented by Fefferman in [13, 14], where six demonstrations are required to accept as valid a solution to the Navier–Stokes 3D Equation [16, 17, 18]. An understanding of the mechanics of the atomic nucleus cannot do without fluid equations.


2. Model

The velocity defined as u=2νPP, with a radius noted as r=x2+y2+z21/2 where Pxyzt is the logistic probability function Pxyzt=11+ektμr, and the expected value Err0<C exist. The term P is defined in xyzR3t0, where constants k>0,μ>0 and Pxyxt is the general solution of the Navier–Stokes 3D equation, which has to satisfy the conditions (1) and (2), allowing us to analyze the dynamics of an incompressible fluid [12, 13, 14].


With, uR3 an known velocity vector, ρ0 constant density of fluid, η dynamic viscosity, ν cinematic viscosity, and pressure p=p0P in xyzR3t0.

Where velocity and pressure are depending of r and t. We will write the condition of incompressibility.


The initial conditions of fluid movement u0xyz, are determined for t=0. Where speed u0 must be C divergence-free vector.


For physically reasonable solutions, we make sure uxyzt does not grow large as r. We will restrict attention to initial conditions u0 that satisfy.


The Clay Institute accepts a physically reasonable solution of (1), (2) and (3), only if it satisfies:


and the finite energy condition [14, 15, 16].

R3uxyzt2dxdydzCforallt0bounded energy.E6

The problems of Mathematical Physics are solved by the Nature, guiding the understanding, the scope, the limitations and the complementary theories. These guidelines of this research were: the probabilistic elements of Quantum Mechanics, the De Broglie equation and the Heisenberg Uncertainty principle.

2.1 Definitions

Nuclear reaction velocity coefficient.

We will use an equation analogous to concentration equation of Physical Chemistry C=C0ekt, where k=p02η, is velocity coefficient, p0 is the initial pressure of our fluid, η the dynamic viscosity and C0 the initial concentration of energetic fluid molecules.

It is evident that, in equilibrium state we can write μr=kt, however, the Navier–Stokes equation precisely measures the behavior of the fluids out of equilibrium, so that: μrkt.

Fortunately, there is a single solution for out-of-equilibrium fluids, using the fixed-point theorem for implicit functions, 11+ektμr=2μr, the proof is proved in Theorem 1.

Attenuation coefficient.

We will use the known attenuation formula of an incident flux I0, for which I=I0eμr. Where, I0 initial flux and μ attenuation coefficient of energetic molecules that enter into interaction and/or resonance with the target molecules, transmitting or capturing the maximum amount of energy [5].

Dimensional analysis and fluid elements.

We will define the respective dimensional units of each one of variables and physical constants that appear in the solution of the Navier–Stokes 3D equation [12, 13, 14, 16].

Nuclear decay Nt=N0ekt. Where k=1T1/2 [1/s], the velocity coefficient, and T1/2 ground state half-life.

Kinematic viscosity ν=ηρ0,m2s.

Dynamic viscosity η, [pa.s], where pa represents pascal pressure unit.

Initial Pressure of out of equilibrium. p0,[pa].

Fluid density ρ0,kgm3, where kg is kilogram and m3 cubic meters.

Logistic probability function, Pxyzt=11+ektμr, it is a real number 0P1.

Equilibrium condition, r=kμt=p02ρ0νμt=uet, [m].

Fluid velocity in equilibrium, ue,m/s. Protons, neutrons and alpha particles are the elements of the fluid.

Fluid field velocity out of equilibrium, u=2νμ1Pr. [m/s]. All nuclear decay is a process out of nuclear equilibrium.

Position, r=x2+y2+z21/2, [m].

Attenuation coefficient, μ,[1/m].

Growth coefficient, k=p02ρ0ν=p02η, [1/s].

Concentration C=C01PP.

Theorem 1The velocity of the fluid is given byu=2νPP,wherePxyztis the logistic probability functionPxyxt=11+ektμr,andppressure such thatp=p0P,both defined onxyzR3t0.The functionPis the general solution of the Navier Stokes equations, which satisfies conditions(1)and(2).

Proof. To verify condition (2), .u=0, we must calculate the gradients and laplacians of the radius. r=xryrzr, and 2r=.r=y2+z2+x2+z2+x2+y2x2+y2+z23/2=2r.


Replacing the respective values for the terms: 2r and r2 in the Eq. (7).


Where the gradient modulus of P=μPP2r, has the form P2=μ2PP22r2=μ2PP22.


Simplifying for 1P0, we obtain the main result of this paper, which represents a fixed point of an implicit function ftr where ftr=P2μr=0. In Nuclear Physics, r0<r<1.2A1/3.


Eq. (10) has a solution according to the fixed-point theorem of an implicit function, and it is a solution to the Navier Stokes stationary equations, which are summarized in: 2P=2μ21r=0. Furthermore, it is the typical solution of the Laplace equation for the pressure of the fluid 2p=p02P=0. Kerson Huang (1987).

To this point, we need to verify that Eq. (10) is also a solution of requirement (1), ut+u.u=ν2upρ0. We will do the equivalence u=θ after we replace in Eq. (1). Taking into account that θ=2νlnP, and that θ is irrotational, ×θ=0, we have: u.u=θ.θ=12θ.θθ××θ=12θ.θ, and 2u=.u××u=.θ××θ=2θ. Simplifying terms in order to replace these results in Eq. (1) we obtain


The explicit form of velocity is u=2μν1Pr. Next, we need the partial derivative ut


After replacing the last four results u.u,2u,ut and pρ0 in Eq. (1) we obtain (11).


The Eq. (11) is equivalent to Eq. (1). After obtaining the term P2P2 from the incompressibility equation 2θ=2νP2P2+2PP=0 and replacing in Eq. (11).


Eq. (10) simultaneously fulfills requirements (1) expressed by Eq. (12) and requirement (2) expressed by Eq. (7), for a constant k=p02ρ0ν=p02η. Moreover, according to Eq. (10), the probability P=2μr which allows the Laplace equation to be satisfied: 2P=2μ21r=0. In other words, the Navier–Stokes 3D equation system is solved.▪

Implicit Function.

An implicit function defined as (10), ftr=11+ektμr2μr=0 has a fixed point tr of R=tr0<atb0<r<+, where m and M are constants, such as: mM. Knowing that the partial derivative exists: rftr=νP1P+2μr2 we can assume that: 0<mrftrM. If, in addition, for each continuous function φ in ab the composite function gt=ftφt is continuous in ab, then there is one and only one function: r=φt continuous in ab, such that ftφt=0 for all t in ab.

Theorem 2An implicit function defined as(10)ftr=11+ektμr2μr=0has a fixed pointtrofR=tr0<atb0<r<+. In this way, the requirements (1) and (2) are fulfilled.

Proof. Let C be the linear space of continuous functions in ab, and define an operator T:CC by the equation:


Then we prove that T is a contraction operator, so it has a unique fixed point r=φt in C. Let us construct the following distance.


Using the mean value theorem for derivation, we have


Where ϕt is situated between φt and ψt. Therefore, the distance equation can be written as:


Using the hypothesis 0<mrftrM we arrive at the following result:


with which we can write the following inequality:


Where α=1mM. Since 0<mM, we have 0α<1. The above inequality is valid for all t of ab. Where T is a contraction operator and the proof is complete, since for every contraction operator T:CC there exists one and only one continuous function φ in C, such that Tφ=φ. Using Eq. (10), which represents the fundamental solution of the Navier–Stokes 3D equation, we verify Eq. (2), which represents the second of the six requirements of an acceptable solution.▪

Proposition 3Requirement (3). The initial velocity can be obtained from:uxyz0=2νPP,where each of the componentsux,uyanduzare infinitely derivable.


Proof. Taking the partial derivatives of xnxr,ynxr and znxr.


Recalling the derivatives of special functions (Legendre), it is verified that there exists the derivative C.


Physically, this solution is valid for the initial velocity, indicated by Eq. (4), where the components of the initial velocity are infinitely differentiable, and make it possible to guarantee that the velocity of the fluid is zero when r [6, 7, 8, 9].

Proposition 4Requirement (4). Using the initial velocity of a moving fluid given byuxyz0=u0xyz=2νμ1P0xryrzr, it is evident that

xαu0CαK1+rKonR3foranyα andK

Proof. Using the initial velocity of a moving fluid given by u0xyz=2νμ1P0xryrzr, we can find each of the components: xαux0,yαuy0 and zαuz0.


For the three components x,y,z the results of the partial derivatives are as follows:


Replacing Eq. (17) with the explanatory form of the Legendre polynomials, for the following terms xα11r and xα1r.


Also, knowing that for each α0, the maximum value of Pα1=1. We can write the following inequality


Grouping terms for xαxr2,yαyr2 and zαzr2 we have the next expressions.


The module of xαu0 is given by xαu0=xαxr2+yαyr2+zαzr21/2. Simplifying and placing the terms of Eq. (20) we have


Taking into consideration that xr1,yr1,zr1 the last term xαu0 can be easily written that.


It is verified that there exists Cα=10α!2 such that if r0, then xαu00. Thus, we proved requirement (4).▪

According to Mathematics, and giving an integral physical structure to the study, we need to prove that there are the spatial and temporal derivatives of the velocity and pressure components, satisfying the requirement (5).

Proposition 5Requirement (5). The velocity can be obtained from:uxyzt=2νPPand each of the componentsux,uyanduzare infinitely derivable.


Proof. Taking partial derivatives for xnxr2,ynxr2 and znxr2.


Recalling the derivatives of special functions, it is verified that the derivative C exists. These derivatives appear as a function of the Legendre polynomials Pn..


There are the spatial derivatives n and the time derivative which is similar to Eq. (25).▪

Proposition 6Requirement (5). The pressure is totally defined by the equivalencepxyzt=p0Pxyztand is infinitely differentiable in each of its components.


Proof. Taking partial derivatives for xn1r,yn1r and zn1r, recalling the derivatives of special functions of Eq. (16), it is shown that the derivative C. We only have to find the time derivatives: tnp0P=p0tnP. Using Eq. (21) for P, we have.


It is always possible to find the derivative tnP as a function of the previous derivative, since the resulting polynomial of each derivative n1 is of degree n.▪

Proposition 7Requirement (6). The energy must be limited in a defined volume and fundamentally it must converge at any time, such thatt0.

R3uxyzt2dxdydzCforallt0bounded energy.

Proof. We will use the explicit form of velocity given in Eq. (21)uxyzt=2νμ1Pr, to obtain the vector module: u2=4ν2μ21P2. Rewriting Eq. (21), and applying a change of variable in: dxdydx=4πr2dr.


Making another change of variable dP=μP1Pdr. Using (10), replacing r2=2μP2 we have


Where radius r, when t0, we have limrP=limr11+expktexpμr=P=1. Moreover, physically if rr00 then t0 we have limr0P=limr011+expktexpμr=P0=12. Here, a probability 12 represents maximum entropy.


In this way the value of the constant C is C=32πν2μ. Verifying the proposition (6) completely. In general, Eq. (10) can be written ftr+r0=11+ektμr+r02μμr+r0=0 and in this way discontinuities are avoided when r0, but this problem does not occur since in the atomic nucleus r0<r<1.2A1/3 is satisfied.▪

Lemma 8The irrotational field represented by the logistic probability functionPxyztassociated with the velocityu=2νPP, can produce vortices, due to the stochastic behavior of the physical variablesp0,η,μ.These stochastic variations are in orders lower than the minimum experimental value.

Proof. The implicit function representing the solution of the Navier–Stokes 3D equation, 11+eμrkμt=2μr depends on the values of initial pressure p0, viscosity ν and attenuation coefficient μ. Due to Heisenberg uncertainty principle, these parameters have a variation when we measure and use them, as is the case of the estimate of ξ=rϕtkμ. The function ϕtkμ=kμt, expressly incorporates these results, when <ξ<+. The physical and mathematical realities are mutually conditioned and allow for these surprising results. For a definite t there exist infinities xyz that hold the relationship r=x2+y2+z21/2. Moreover, for a definite r there are infinities t that respect the fixed-point theorem and create spherical trajectories. When the physical variables k,μ vary, even at levels of 1/100 or 1/1000, they remain below the minimum variation of the experimental value. We could try to avoid the existence of trajectories on the spherical surface, for which we must assume that the fluid is at rest or it is stationary, which contradicts the Navier–Stokes 3D equation, where all fluid is in accelerated motion ut0. In short, if there are trajectories in the sphere as long as it is probabilistically possible, this is reduced to showing that the expected value of the radius Err0 exists and it is finite.

Derivation of Err0.

The logistic density function for ξ when Eξ=0 and Varξ=σ2 is defined by.

hξ=μexpξ1+expξ2, where 1μ=σ3/π is a scale parameter. Given that r=ϕtp0ημ+ξ function for r is then fr=exprϕ/ττ1+exp(rϕ/τ2, to facilitate the calculations we put ϕ=ϕtkμ=kμt. By definition, the truncated density for r when r0 is given by frr0=frPr0 for r0. Given that the cumulative distribution function for r is given by Fr=11+expktμr, it follows that Pr0=1F0=expϕ)1+expϕ)=11+expϕ). The derivation of Err0 then proceeds as follows:


We replaced in Eq. (29)dP=μP1Pdr and r2=2μP2 of this manner we obtain


where we have used the fact that


where the last equality follows from an application of the L’Hopital’s rule Pr0=limtexpkt1+expkt=1.▪


3. Results

The main results of applying the Navier Stokes equation to the atomic nucleus, which behaves like an incompressible nuclear fluid, are:

  • The nuclear force and the Navier Stokes force are related.

  • The Cross Sections in Low energy X ray can explain the Golden Ratio σ1σ2 that appears in the femtoscope.

  • Navier Stokes Equation and Cross Section in Nuclear Physics.

  • The principle of the femtoscope explains that low energy X-rays produce resonance in layer K.

3.1 The nuclear force and the Navier Stokes force are related

Firstly, we will use the concepts of Classic Mechanics and the formulation of the Yukawa potential, Φr=g4πrA1eμr to find the nuclear force exerted on each nucleon at interior of the atomic core FN=Φr. Also, replace the terms of the potential eμr=1PP and 1r=μ2P by the respective terms already obtained in Eq. (10).


The general form of the Eq. (32), is a function of xyzt.


Secondly, we will obtain the Navier Stokes force equation given by:


Theorem 9The Nuclear Force and Navier Stokes Force are proportional inside the atomic nucleusFN=CFNS.

Proof.Eq. (34) rigorously demonstrated by theorems and propositions 1 through 8, represent the acceleration of a particle within the atomic nucleus. According to Classical Mechanics the force of Navier Stokes applied to a particle of mass m, would have the form:


Proof. The nuclear force on its part would be calculated as follows FN=Φr.


Replacing the term P=μPP2r of Eq. (7), we obtain


It is possible to write nuclear force as a function of speed.


Finally, we can show that the nuclear force and force of Navier Stokes differ at most in a constant C. Equating (35) and (37), we find the value g as a function of the parameters nuclear viscosity ν, attenuation μ and growth coefficient of the nuclear reaction k, nucleon mass m and C1.


3.2 Cross section and Golden ratio σ1σ2 are important elements of the femtoscope

According to NIST and GEANT4 [17], current tabulations of μρ rely heavily on theoretical values for the total cross section per atom, σtot, which is related to μρ by the following equation:


In (Eq. 39), u=1.6605402×1024gr is the atomic mass unit (1/12 of the mass of an atom of the nuclide 12C)4.

The attenuation coefficient, photon interaction cross sections and related quantities are functions of the photon energy. The total cross section can be written as the sum over contributions from the principal photon interactions


Where σpe is the atomic photo effect cross section, σcoh and σincoh are the coherent (Rayleigh) and the incoherent (Compton) scattering cross sections, respectively, σpair and σtrip are the cross sections for electron-positron production in the fields of the nucleus and of the atomic electrons, respectively, and σph.n. is the photonuclear cross section 3,4.

We use data of NIST and simulations with GEANT4 for elements Z=11 to Z=92 and photon energies 1.0721×103 MeV to 1.16×101 MeV, and have been calculated according to:


The attenuation coefficient μ of a low energy electron beam 10,100eV will essentially have the elastic and inelastic components. It despises Bremsstrahlung emission and Positron annihilation.


3.3 The principle of the femtoscope explains that low energy X-rays produce resonance in K-shell

A resonance region is created in a natural way at the K-shell between the nucleus and the electrons at S-level. The condition for the photons to enter in the resonance region is given by rarn+λ. This resonance region gives us a new way to understand the photoelectric effect. There is experimental evidence of the existence of resonance at K-level due to photoelectric effect, represented by the resonance cross section provided by NIST and calculated with GEANT4 for each atom. In the present work we focus on the resonance effects but not on the origin of resonance region.

The resonance cross section is responsible for large and/or abnormal variations in the absorbed radiation I2I1.


Theorem 10Resonance region. The resonance cross section is produced by interference between the atomic nucleus and the incoming X-rays inside the resonance region, where the boundaries are the surface of the atomic nucleus and K-shell.

The cross section of the atomic nucleus is given by:


The photon cross section at K-shell depends on the wave length and the shape of the atomic nucleus:


Subtracting the cross sections (13) and (14) we have:


The resonance is produced by interactions between the X-rays, the K-shell electrons and the atomic nucleus. The cross sections corresponding to the nucleus is weighted by probability pn and should have a simple dependence of an interference term. This last depends on the proton radius rp or the difference between the nucleus and proton radius (rnrp) according to the following relation


We note that left hand side of Eqs. (46) and (47) should have a factor larger than one due to resonance. The unique factor that holds this requirement is σ1σ2b Where a,b constants.


After performing some simulations it is shown that the thermal a represents the dimensionless Rydberg constant a=R=1.0973731568539107.


We use σ1Z,σ2Z for represent cross section in resonance, R is the generalized Rydberg constant for all elements of periodic table, and Z atomic number.

Last equation with a R2=0.9935 was demonstrated and constructed using the elements of the solution of the Navier Stokes equations.


This equation was obtained with a R2=0.9939, and indicates that the ratio of the effective sections fully explain each element of the periodic table.


This equation obtained for a R2=0.9996, complements the system of equations that allow to know the simulation values as a function of Z for σ1Z,σ2Z and EZ, where σ1Z<σ2Z. The Femtoscope equations further demonstrate that energy E=minE is minimum and Shannon entropy S=maxS is maximum in resonance, because in equilibrium σ1Zi=σ2Zi.

The radius of the neutron can be obtained using Eq. (49) in the following way.


3.4 Navier Stokes equation and cross section in nuclear physics

The speed needs to be defined as u=2νPP, where Pxyzt is the logistic probability function Pxyxt=11+ektμr,r=x2+y2+z21/2 defined in xyzR3t0 This P is the general solution of the Navier Stokes 3D equations, which satisfies the conditions (50) and (51), allowing to analyze the dynamics of an incompressible fluid.


Where, uR3 an known velocity vector, ρ0 constant density of fluid and pressure p=p0PR.

With speed and pressure dependent on r and t. We will write the condition of incompressibility as follows.


Theorem 11The velocity of the fluid given by:u=2νPP,wherePxyztis6the logistic probability functionPxyxt=11+ektμx2+y2+z21/2,defined inxyzR3t0is the general solution of the Navier Stokes equations, which satisfies conditions(50)and(51).

Proof. Firstly, we will make the equivalence u=θ and replace it in Eq. (50). Taking into account that θ is irrotational, ×θ=0, we have.


We can write,


It is equivalent to,


where Δp is the difference between the actual pressure p and certain reference pressure p0. Now, replacing θ=2νlnP, Navier Stokes equation becomes.


The external force is zero, so that there is only a constant force F due to the variation of the pressure on a cross section σ. Where σ is the total cross section of all events that occurs in the nuclear surface including: scattering, absorption, or transformation to another species.


putting (52) in (53) we have


In order to verify Eq. (51), .u=0, we need to obtain r=xryrzr,2r=.r=y2+z2+x2+z2+x2+y2x2+y2+z23/2=2r.


Replacing the respective values for the terms: 2P and P2 of Eq. (55). The Laplacian of P can be written as follows.


Using gradient P=μPP2r, modulus P2=μ2PP22r2 and 2P in (56).


Replacing Eqs. (55) and (56) in (57) we obtain the main result of the Navier Stokes equations, the solution represents a fixed point of an implicit function ftr where ftr=P2μr=0.


An important result of the Navier Stokes 3D equation, applied to the nuclear fluid of an atom, allows us to advance our understanding of nuclear dynamics and nuclear force.

Corollary 12The nuclear decay constantkis determined by the nuclear pressurep0,and the dynamic nuclear viscosityη,as follows:k=p02η.

Proof. A nuclear decay is a reaction of degree 1, which is explained by the exponential law Nt=N0ekt. Rewriting Eq. (58) as a function of the initial nuclear pressure p0 and the nuclear viscosity η we can prove that they are related to the nuclear decay constant, in an intrinsic way and allow to explain dynamic nuclear phenomena.


3.5 Model of nuclear pressure

The work in the nuclear fluid is given by the variation of the energy necessary to form that atomic nucleus, that is, by the excess mass required in the process. Using Eq. (59) we can find the explicit form of nuclear pressure p=p0P=p02μr.


Applying the mean value theorem of integrals, we know that there is a mean value of the nuclear pressure p¯ and the volume of the atomic nucleus V¯ of the integral (60)


that according to Quantum Mechanics, pζ,Vr0 are the observable values. By this restriction we can equalize Eq. (60), (61) and find the value of the initial nuclear pressure as a function of the experimentally measured nuclear pressure: pζ4π3r03=4πp0r02μ with r0=1.2A1/2fm, from where:


For Yukawa’s potential, it is often assumed μ1r0.

3.6 Model of nuclear viscosity

Using the fundamental expression k=p02η=1T1/2 obtained from the resolution of the Navier Stokes 3D equations, Eqs. (58) and (59), it is possible to find the average or most probable values of the variables involved pressure, nuclear viscosity and nuclear decay constant as follows: k¯=p¯2η¯.

Explicitly we find the nuclear viscosity as a function of the nuclear decay constant k=k¯ and the average value of the nuclear pressure p¯ as follows:


There is another experimental way of determining nuclear viscosity, through the fuzziness of alpha particles, protons or neutrons ejected in a nuclear decay using the fluid velocity equation u=2νμ1Pr, in which modulo u=2νμ1Pr, replacing the dynamic viscosity ν=ηρ0m2s, we obtain: u=2ηρ0μ1P.

So the second way to find the nuclear viscosity has the form:


Eq. (64) has a more complicated form and depends on the speed of the fluid particles and the radius from which these particles leave.

3.7 Calculation of nuclear pressure and viscosity

The nuclear decay constant k is determined by the nuclear pressure p0, and the dynamic nuclear viscosity η, as follows: k=p02η.


4. Discussion of results

4.1 Calculation of the pressure relation of proton and neutron

The fine-structure constant, α, has several physical interpretations, we use the most known.

The fine-structure constant, α=1137.03599917435 is the ratio of two energies: the energy needed to overcome the electrostatic repulsion between two protons a distance of 2r apart, and (ii) the energy of a single photon of wavelength λ=2πr (or of angular wavelength r; according to Planck relation).


We can find the relationship of energies between the proton and the atomic nucleus. Knowing that the two occupy the same nuclear volume V. This relationship is identical to 1α', since the atomic nucleus interacts with the proton through the electromagnetic field and the nuclear force.


Thus, we already know the pressure relation between the proton and the atomic nucleus pppN=2α.

Now we find the relation of pressures between the neutron and the atomic nucleus, which are under the action of the same nuclear force, F.


For the chemical element with maximum nuclear pressure, 2862Ni we have pnpN=1.2621/320.841842=31.830.

If we divide Eq. (67) for Eq. (66) we have.


4.2 An action on the nuclear surface produces a reaction in the nuclear volume and vice versa

It is created by the friction between the layers of nucleons.

Theorem 13An action on the nuclear surface produces a reaction in the nuclear volume and vice versa.

Proof. The volume and the nuclear surface are connected through the Gaussian divergence theorem and the Navier Stokes equations.

For an incompressible fluid, whose velocity field uxyz is given, .u=0 is fulfilled.

Logically, the integral of this term remains zero, that is:


Writing the Divergence theorem. u.ndS=.udxdydz=0, the first term must be equal to zero, that is:


The only possible trajectory is circular, because in this case the vector n is perpendicular to the surface of the sphere. In this way the equation of the outer sphere corresponding to the surface is: x2+y2+z2=1.2A1/3.

Within the nuclear fluid there are layers of nucleons that move in spherical trajectories.▪


5. Conclusions

High energy physics is the guide to low energy physics, because, in certain processes such as in the measurement of the internal pressure of protons and neutrons. However, we demonstrate that there is a trend compatibility of the two in the characterization of the atomic nucleus.

The Femtomathematics corresponds to the tools and the logic that allows to calculate the parameters of the order of 1015m. In an indirect way, it is in full correspondence with Femtophysics.

The nuclear viscosity at equilibrium is much larger than the nuclear viscosity at the moment of nuclear decay, which is totally logical, because at equilibrium the atomic nucleus is totally compact in pressure and density, while at the moment protons, neutrons, alpha particles come out of the atomic nucleus from a nuclear decay, indicating that the nuclear viscosity decreased.



  1. To the participants of the congress: The 4th International Conference on Materials Sciences and Nanomaterials (ICMSN 2020), Cambridge, United Kingdom during July 08-10, 2020.

  2. To the research directorate (DI) of the Central University of Ecuador, as well as the group of the

GIIP process of the Faculty Chemical Engineering, UCE.


Conflicts of interest

The author declare no conflict of interest.


  1. 1. Burkert VD. Elouadrhiri, L. & Girod, F.X (2018). The pressure distribution inside the proton. Nature. 2018;557:396-399. DOI:
  2. 2. P. E. Shanahan and W. Detmold (2019). Pressure Distribution and Shear Forces inside the Proton. Phys. Rev. Lett. 122, 072003 – Published 22 February 2019
  3. 3. Ozel F, Freire P. Masses, Radii, and Equation of State of Neutron Stars. Ann. Rev. Astron. Astrophys. 2016;401-440(2016):54
  4. 4. Pagels, H (1966). Energy-Momentum Structure Form Factors of Particles", Phys. Rev. 144, 1250–1260 (1966).
  5. 5. Ji, X. D (1997). Deeply virtual Compton scattering", Phys. Rev. D 55, 7114–7125 (1997).
  6. 6. Teryaev, O. V (2016). Gravitational form factors and nucleon spin structure", Front. Phys. (Beijing) 11, 111207-111207-8 (2016).
  7. 7. Yeverechyahu A. Nuclear viscosity and widths of giant resonances. Annals of Physics. Volume 95, Issue 1. November. 1975:1975
  8. 8. Wong. Introductory to Nuclear Physics. WILEY-VCH Verlag Gmbl I & Co. Weinheim: KGaA; 2004. p. 2004
  9. 9. Loveland W, Morrissey D, Seaborg GT. Modern Nuclear Chemistry. Wiley-Interscience. p. In: ISBN 978–0–471-11532-8. 2006
  10. 10. Pohl R et al. The size of the proton. Nature. 2010;2010(466):213-216
  11. 11. Jiménez, Recalde, Jimenez Chacon (2017). Extraction of the Proton and Electron Radii from Characteristic Atomic Lines and Entropy Principles. Entropy 2017.
  13. 13. Kulish VV. On the Relationship between Fluid Velocity and de Broglie’s Wave Function and the Implications to the Navier – Stokes Equation. International Journal of Fluid Mechanics Research. 2002;29(1):2002
  14. 14. Leray J, l’espace S l m d’u l v e. Acta Math. J. 1934:63
  15. 15. Caffarelli L, Kohn R, Nirenberg L. Partial regularity of suitable weak solutions of the Navier–Stokes equations. Comm. Pure & Appl. Math. 1982;35:771-831
  16. 16. Kerson Huang (1987). Statistical Mechanics, 2nd Edition, John Willey & Sons, 1987.
  17. 17. Geant 4. (2016). Geant4 User’s Guide for Application Developers Version: geant4 10.3. Publication date 9 December 2016.
  18. 18. Oganessian and Utyonkov. Superheavy nuclei from 48Ca-induced reactions. Nuclear Physics A. 2015;944(2015)

Written By

Edward Henry Jimenez

Submitted: 09 September 2020 Reviewed: 07 December 2020 Published: 15 January 2021