1. Introduction: Shortcomings of Boltzmann Physical Kinetics
In 1872 L Boltzmann published his famous kinetic equation for the one-particle distribution function (DF)  . He expressed the equation in the form
where is the collision (“stoß”) integral, and
is the substantial (particle) derivative, and being the velocity and radius-vector of the particle, respectively. Transport processes in open dissipative systems are considered in physical kinetics. Therefore, the kinetic description is inevitably related to the system diagnostics. Such an element of diagnostics in the case of theoretical description in physical kinetics is the concept of the physically infinitely small volume ( ). The correlation between theoretical description and system diagnostics is well-known in physics. Suffice it to recall the part played by test charge in electrostatics or by test circuit in the physics of magnetic phenomena.
Let us consider the hydrodynamic description in more detail from this point of view. Assume that we have two neighboring physically infinitely small volumes and in a non-equilibrium system. The one-particle distribution function (DF) corresponds to the volume , and the function to the volume . It is assumed in a first approximation that does not vary within , same as does not vary within the neighboring volume . This assumption of locality is implicitly contained in the Boltzmann equation (BE)  - . However, the assumption is too crude; PhSV is an open thermodynamic system.
The Boltzmann equation (BE) fully ignores non-local effects and contains only the local collision integral . But these nonlocal effects are irrelevant only in equilibrium systems, where the kinetic approach goes into methods of statistical mechanics. As a result, the difficulties of classical Boltzmann physical kinetics arise. Conclusion:
1) Kinetic theory should be non-local.
2) The effect is of the order of Knudsen number.
3) The effect is due to the reduced description and is not related to the specific division of the physical system by the PhSV grid.
4) Accurate derivation of the kinetic equation with respect to the one-particle DF should lead to corrections of the order of the Knudsen number before uncoupling the Bogolyubov chain.
5) This means that in the Boltzmann equation, the terms of the order of the Knudsen number are lost, significant for both large and small Knudsen numbers.
6) The Boltzmann equation does not even belong to the class of minimal models, being only a “plausible” equation.
7) Boltzmann equation in this sense is the wrong equation.
2. Generalization of the Boltzmann Kinetic Equation. Nonlocal Physical Kinetics
A rigorous approach to the derivation of the kinetic equation for the distribution function (DF) ( ) is based on the hierarchy of the Bogolyubov- Born-Green-Kirkwood-Yvon (BBGKY) equations     . The structure of the is generally as follows
where is the non-local integral term incorporating in particular the time delay effect. The generalized Boltzmann physical kinetics, in essence, involves a local approximation
for the second collision integral, here being proportional to the mean time between the particle collisions. All of the known methods  -  of deriving kinetic equation relative to one-particle DF lead to approximation (2.2), including the method of many scales, the method of correlation functions, and the iteration method. We can draw here an analogy with the Bhatnagar-Gross-Krook (BGK,  ) approximation for ,
which popularity as a means to represent the Boltzmann collision integral is due to the huge simplifications it offers.
In other words, the local Boltzmann collision integral admits approximation via the BGK algebraic expression (2.3), but more complicated non-local integral can be expressed as differential form (2.2). The ratio of the second to the first term on the right-hand side of Equation (2.1) is given to an order of magnitude as and at large Knudsen numbers (defining as ratio of mean free path of particles to the character hydrodynamic length) these terms become of the same order of magnitude. It would seem that at small Knudsen numbers answering to hydrodynamic description the contribution from the second term on the right-hand side of Equation (2.1) is negligible.
This is not the case,however.When one goes over to the hydrodynamic approximation (by multiplying the kinetic equation by collision invariants and then integrating over velocities), the Boltzmann integral part vanishes, and the second term on the right-hand side of Equation (2.1) gives a single-order contribution in the generalized Navier-Stokes description. Mathematically, we cannot neglect a term with a small parameter in front of the higher derivative. Physically, the appearing additional terms are due to viscosity and they correspond to the small-scale Kolmogorov turbulence .
The integral term turns out to be important both at small and large Knudsen numbers in the theory of transport processes.
Thus, is the distribution function fluctuation, and we find
Writing Equation (2.4) without taking into account Equation (2.5) makes the BE non-closed. From viewpoint of the fluctuation theory, Boltzmann employed the simplest possible closure procedure . Then, the additional GBE terms (as compared to the BE) are significant for any Kn, and the order of magnitude of the difference between the BE and GBE solutions is impossible to tell beforehand. For GBE the generalized H-theorem is proven .
Madelung’s quantum hydrodynamics is equivalent to the Schrödinger equation (SE) and furnishes the description of the quantum particle evolution in the form of Euler equation and continuity equation. Madelung’s interpretation of SE (connected with wave function ) leads to the probability density
and velocity . Schrödinger-Madelungquantum hydrodynamics does not lead to the energy equation in principle; then the Schrödinger-Madelungquantum mechanics cannot be used for the description of the dissipative processes.
The dependent variable p in the energy equation of the generalized quantum hydrodynamics can be titled as the rest quantum pressure or simply quantum pressure.
Generalized Boltzmann physical kinetics brings the strict approximation of non-local effects in space and time and after transfer to the local approximation leads to parameter , which on the quantum level corresponds to the uncertainty principle “time-energy”.
The appearance of the nonlocal parameter is consistent with the Heisenberg uncertainty relation. But in principle generalized nonlocal kinetic equations (and therefore GHE) needn’t in using of the “time-energy” uncertainty relation for estimation of the value of the non-locality parameter . Moreover, the “time-energy” uncertainty relation does not produce the exact relations and from position of non-local physics is only the simplest estimation of the non-local effects. Really, let us consider two neighboring physically infinitely small volumes and in a non-equilibrium system. Obviously the time should tends to diminish with increasing of the velocities u of particles invading in the nearest neighboring physically infinitely small volume ( or ):
But the value cannot depend on the velocity direction and naturally to tie with the particle kinetic energy, then
where is a coefficient of proportionality, which reflects the state of physical system. In the simplest case is equal to Plank constant and relation (2.7) became compatible with the Heisenberg relation. The non-locality parameter plays the same role as the transport coefficients in local hydrodynamics. The different models can be introduced for the definition, but the corresponding results not much different like in local kinetic theory for different models of the particles interaction.
In the general case, the parameter is the non-locality parameter; in quantum hydrodynamics, its magnitude is correlated with the Heisenberg “time-energy” uncertainty relation   and with the Ehrenfest adiabatic theorem .
Now we can turn our attention to the quantum hydrodynamic description of individual particles. The abstract of the classical Madelung’s paper  contains only one phrase: “It is shown that the Schrödinger equation for one-electron problems can be transformed into the form of hydrodynamic equations”. The following conclusion of principal significance can be done from the previous consideration:
Madelung’s quantum hydrodynamics is equivalent to the Schrödinger equation (SE) and leads to the description of the quantum particle evolution in the form of Euler equation and continuity equation. Generalized Boltzmann physical kinetics leads to the strict approximation of non-local effects in space and time and in the local limit leads to parameter , which on the quantum level corresponds to the uncertainty principle “time-energy”. Generalized hydrodynamic equations (GHE) lead to SE as a deep particular case of the generalized Boltzmann physical kinetics and therefore of non-local hydrodynamics.
Finally of Item 2 we can state that introduction of control volume by the reduced description for ensemble of particles of finite diameters leads to fluctuations (proportional to Knudsen number) of velocity moments in the volume. This fact leads to the significant reconstruction of the theory of transport processes. The violation of Bell’s inequalities  is found for local statistical theories, and the transition to the non-local description is inevitable.
3. Basic Quantum Nonlocal Hydrodynamic Equations. Superconducting Soliton Motion in the Two Component Physical System
In general case the strict consideration leads to the following system of the non-local quantum hydrodynamic equations written in the Generalized Hydrodynamic Form (GHE) for multi-component species (see also  -  ):
Continuity equation for species :
Continuity equation for mixture:
Momentum equation for species :
Momentum equation for mixture
Energy equation for species
Energy equation for mixture:
Here are the forces (acting on the mass unit) of the non-magnetic origin, —magnetic induction, —unit tensor, —charge of the -component particle, —static pressure for -component, —internal energy for the particles of -component, —hydrodynamic velocity for mixture, —non-local parameter.
In the following we intend to obtain the soliton’s type of solution of the generalized hydrodynamic equations (GHE). The non-stationary 1D model will be used with taking into account the energy equation, external forces and non-locality parameter defined by the “time-energy” uncertainty relation of Heisenberg. Then GHE contain Poisson equation (reflected fluctuations of charges and flux of the charges density), two continuity equations for positive (lattice ions) and negative (electrons) species, momentum equation and two energy equations for positive and negative species. This system of six non-stationary 1D equations is written as:
(Generalized Poisson equation):
(Continuity equation for ions):
(Continuity equation for electrons):
(Energy equation for ions):
(Energy equation for electrons):
where u is velocity of the directed motion of combined quantum object (phonon-electron), and —numerical density of the charged species, and —forces (of potential and non-potential origin), acting on the mass unit of the charged particles. The right hand sides of the energy equations are written in the relaxation forms following from BGK kinetic approximation.
For acting potential forces of the electrical origin the relations are valid
where —scalar potential.
Let us introduce approximations for and using (2.7)
, . (3.15)
For electron-phonon non-local parameter the following relation is applicable
Formula (3.17) is obvious consequence of uncertainty relation for combined particle which mass is . Energy equation of the generalized quantum hydrodynamics contains pressures , which can be named as the quantum pressure of the non-local origin. In the definite sense these pressures can be considered as analog of the Bose condensate pressure.
4. Combined Quantum Solitons in the Self-Consistent Electric Field
Let us formulate the problem in detail. The non-stationary 1D motion of the combined phonon-electron soliton is considered under influence of the self- consistent electric forces of the potential and non-potential origin. It should be shown that mentioned soliton can exists without a chemical bond formation. We introduce the coordinate system moving along the positive direction of the x axis in 1D space with the velocity , which C is equal to the phase velocity of this quantum object.
Taking into account de Broglie relation we should wait that the group velocity is equal to . Really the energy of a relativistic particle is
and c is the light velocity, is the group velocity, —the particle rest mass. Relation (4.2) can be written as
is the particle impulse. In the non-relativistic approach the relation (4.4) takes the form
Using the dualism principle in the de Broglie interpretation we have for the particle energy
where is the circular frequency, —the phase velocity, is the wave number and is the wave length. Correspondingly the particle impulse p is
and using (4.8),
Then in the non-relativistic case we have
From (4.9) and (4.10) for the non-relativistic case one obtains
Then we should wait that the indestructible soliton has the velocity in the coordinate system moving with the phase velocity . If we pass on the moving coordinate system, all dependent hydrodynamic values will be functions of . But we investigate the possibility of the creation of the combined quantum object of the soliton type. For this case the explicit time dependence of solutions does not exist in mentioned coordinate system moving with the phase velocity .
Write down the system of Equations (3.7)-(3.12) for the two component mixture of charged particles (without taking into account the component’s internal energy) in the dimensionless form, where dimensionless symbols are marked by
tildes. We introduce the scales for velocity , for coordinate x , for the potential scale and for the density scale (e is absolute electron charge).
Generalized Poisson Equation (3.7) takes the form
Scaling forces are
Analogical transformations should be applied to other equations of the system (3.7)-(3.12). As result one obtains the six non-linear dimensionless ordinary differential equations
Some comments to Equations (4.12, 4.15-4.20):
1) Equations. (4.12, 4.15-4.20) contain 6 dependent variables: . Every equation from the system (4.12, 4.15-4.20) is of the second order and needs two conditions. The problem belongs to the class of Cauchy problems.
2) In comparison with the Schrödinger theory connected with behavior of the wave function, no special conditions are applied for dependent variables including the domain of the solution existing. This domain is defined automatically in the process of the numerical solution of the concrete variant of calculations.
3) From the introduced scales
, , , , (4.21)
only one parameter is independent – the phase velocity of the combined quantum object. The value can be titled as quantum kinematic viscosity . From this point of view the obtained solutions which will be discussed below have the universal character defined only by Cauchy conditions.
5. Maple Program Realizing the Solution of the System of Equations (4.15)-(4.20)
The system of generalized quantum hydrodynamic Equations (4.12), (4.15)-(4.20) have the great possibilities of mathematical modeling as result of changing of twelve Cauchy conditions describing the character features of initial perturbations which lead to the soliton formation.
On this step of investigation we intend to demonstrate the influence of difference conditions on the soliton formation. The following figures reflect some results of calculations realized according to the system of Equations (4.12), (4.15)-(4.20) with the help of Maple. The following notations on figures are used: r—density , s—density , u—velocity , p—pressure , q—pressure and v—self consistent potential . Explanations placed under all following figures, Maple program contains Maple’s notations, for example the expression
means in usual notations , independent variablet responds to .
There is the problem of principle significance—is it possible after a perturbation (defined by Cauchy conditions) to obtain the quantum object of the soliton’s kind as result of the self-organization of ionized matter? In the case of the positive answer, what is the origin of existence of this stable object?
Appendixes 1 and 2 contain the corresponding Maple programs. The programs are ready for using.
6. Super-Conductivity as Movement of Solitons without Destruction
Let us demonstrate some examples of application of non-local GHE with the help of the Maple program in the Item 5. Significant remarks:
1) We investigate the wave movement of two component ionized matter in the form of moving solitons. It is known that the Schrödinger-Madelung quantum mechanics leads to the soliton destruction.
Extremely important that the really working super-conductivity theory should lead to the soliton conservation in the frame of Cauchy problem,but not as a boundary problem.It means that we should observe the self-organization of matter.
2) The program contains the following Cauchy dimensionless parameters:
—self-consistent potential and , —velocity and , density of negative particles and , density of positive particles and , pressure of negative particles and , pressure of positive particles and , and the mass ratio of hard (H) and light (L) particles. The program (Appendix 1) is ready for calculations.
3) Let us demonstrate some calculations changing the pressure parameters by ten orders of magnitude. We introduce the dimensionless soliton energetic temperature using the definition
Analogically we have
Figures 1-6 reflect the calculations for the Cauchy conditions
Figures 7-12 reflect the calculations for the Cauchy conditions
Table 1 contains the Cauchy conditions for pressure, the boundaries (designated as lim1 and lim2) of the soliton, dimensionless energetic temperature for positive particles and dimensionless concentration of the captured electrons.
From our results obtained in the frame of nonlocal physics follow:
1) Nonlocal quantum hydrodynamics leads to the following physical picture of super-conductivity. So called “quasi-particles” are in reality the solitons which cannot be discovered in Schrödinger-Madelung theory in principal.
Figure 1. Dimensionless energetic temperature R for positive particles and dimensionless concentration S of the captured electrons.
Figure 2. q— pressure of negative particles, p— pressure of positive particles.
Figure 3. u—velocity , r—density .
Figure 4. Velocity u— .
Figure 5. s—density .
Figure 6. Self-consistent potential v— .
Figure 7. Dimensionless energetic temperature R for positive particles and dimensionless concentration S of the captured electrons.
Figure 8. q— pressure of negative particles, p— pressure of positive particles.
Figure 9. u—velocity .
Figure 10. r—density .
Figure 11. s—density .
Figure 12. Self-consistent potential v— .
Table 1. Cauchy conditions for pressure, the boundaries (designated as lim1 and lim2) of the soliton, dimensionless energetic temperature for positive particles and dimensionless concentration of the captured electrons.
2) The stability of the quantum object is result of the self organization of the ionized matter. In other words—it is self-consistent influence of electric forces and quantum pressures. Solitons defines the extremely “stiff” structure. Look at Figure 4 and Figure 9—all solitons parts are moving with the same velocity.
3) Stability can be achieved if soliton has negative “shell” and positive “nuclei” and , or if soliton has positive “shell” and negative “nuclei” and . The moving plane waves remind the structure of “flat atoms” or paired electrons in BCS theory. The previous figures display the typical quantum objects placed in the bounded region of 1D space, all parts of these objects are moving with the same velocity. Namely from calculations follow that in coordinate system moving with the phase velocity, indestructible soliton has the velocity equal to the phase velocity. Moreover the attempt to impose to soliton to move with another group velocity leads to the soliton destruction.
4) Table 1 demonstrates practically the linear dependence between energetic temperature and quantity of captured electrons. All quantum fluctuations in quantum GHE (see Equations 3.1-3.6, 3.7-3.12 and relations 3.15-3.17) are proportional to the Planck constant .
7. The Charge Carrier’s Behavior near the Critical Point
The step function for profiles in the vicinity of the critical temperature can be obtained analytically from the nonlocal GHE equations. Really let us consider the generalized continuum Equations (4.16) and (4.17). Equation
can be immediately integrated
But for all parts of solitons velocity . Then
In the vicinity of we have and Equation (7.3) is written as follows
After (7.5) integration we have practically step function.
Analogically for equation
Some remarks concerning the scales relations. Let us remind Homes-Zaanen relation    
or taking into account that in the soliton’s boundaries
we find (see (3.15))
where is the phase velocity of the soliton’s movement. An attempt to impose a velocity different from the phase velocity using Cauchy conditions leads to the destruction of the soliton, and hence to the destruction of superconductivity. Then
and we reach the scaling law (7.9) in the frame of nonlocal physics
This time is very short. By analogy with the theory of gravity, it was called the “Planck time” for dissipation. The attempt to tie the superconductivity effects with the cosmological problems seems unacceptable only from the first glance. The origin of this analogy is connected with nonlocal physical kinetics which applicable from the Planck time in the Big Bang theory to the Universe evolution, .
It means that we find the phase velocity
or (the value can be titled as quantum kinematic viscosity ,
Planck constant , Boltzmann constant )
, where is the temperature measured in Kelvin.
The time scale we used
Let us introduce and calculate the corresponding Reynolds number
In other words we consider the flow of nonlocal quantum liquid by the condition .
8. The Generalized Non-Local London’s Theory. Longitudinal Electro-Magnetic Waves. The Basic Equations
We use generalized nonlocal Maxwell equations to prove the existence of longitudinal electromagnetic waves. To this end let’s differentiate in time both parts of the Maxwell equation.
where , is fluctuation of the current density. We find
Obviously the second term in the square bracket can be omitted with the good accuracy, we find
It is known that the double vector product is
Then the Equation (8.4) is written as follows
or if the electric permeability does not depend on time we have
We write (8.10) in the form
where (see   )
Then the nonlocal Equation (8.12) takes the form
where is the coefficient of conductivity. In the simplest Drude model
where n is numerical electron density, is effective mass, is the relaxation time, the dimension of the conductivity is s−1. As we see from the definition of , the introduction of an “effective” mass translates the relaxation timeto the level of the fitting parameter. If nonlocal parameter and the coefficient of conductivity , we have
Equation (8.19) is the basement relation for following investigation.
9. The Nonlocal Theory of Longitudinal Electromagnetic Waves. The 1D Non-Stationary Case
Let us investigate the 1D non-stationary case. For the 1D case (x—direction), we find from (8.19)
After integration we find
where the Boltzmann constant .
Let us consider the main features of the solution of Equation (9.2). With this aim we suppose that , and the last term of the equation can be omitted. The penultimate term in Equation (9.2) can be omitted if the pressure gradient can compensate the influence of the gravitation force acting on the volume unit and . We have
The solution of this equation is a damping longitudinal E-wave. Really
1) The validity of the solution (9.4) could be verified by the direct substitution of this relation (9.4) into Equation (9.3).
2) If the non-locality parameter is equal to zero, the electric intensity for the longitudinal waves (EMLW) also turns into zero. It means that EMLW cannot exist in the frame of the classic Maxwell electrodynamics.
Equation (9.2) can be written in the form
The last term in square bracket takes into account the influence of the pressure gradient and the gravitation force. If nonlocal parameter , then
and we have in the local electrodynamics the exponential attenuation (if in (9.2)) without the wave creation
The existence of longitudinal electromagnetic waves does not contradict (nonlocal) electrodynamics if longitudinal waves are actually detected and the medium in which they are generated is known. Then in this case the problem of substantiating such waves in electrodynamics is reduced to the search for material equations characterizing the response of the medium to the influence of the field and the joint solution of these material equations and Maxwell’s equations, with appropriate boundary conditions.
Low permittivity (low epsilon) materials are now attracting wide attention due to potential novel applications in optics and radio communications. Surface plasmon polaritons (SPPs) are electromagnetic waves that travel along a metal-dielectric or metal-air interface, practically in the infrared or visible-frequency. The term “surface plasmon polariton” explains that the wave involves both charge motion in the metal (“surface plasmon”) and electromagnetic waves in the air or dielectric (“polariton”). Application of SPPs enables subwavelength optics in microscopy and lithography beyond the diffraction limit. It also enables the first steady-state micro-mechanical measurement of a fundamental property of light itself: the momentum of a photon in a dielectric medium.
Longitudinal electromagnetic waves play important role in plasma, in surface plasmon polaritons in anisotropic materials, in space-charge waves in semiconductor materials. Then for so called that in Epsilon Near Zero (ENZ) materials the velocity (the usual light velocity). For the following details of the wave process one needs the explicit form of the function and boundary and initial conditions. As we see from (9.10)
if . (9.11)
Obviously to overcome this velocity limit is possible only in the systems with negative conductivity where formally
if . (9.12)
The possibility of the appearance of negative conductivity in a non-equilibrium electron system, i.e., a situation in which the current flows opposite to the electric field, was apparently indicated for the first time by Krömer in the late 1950s . The mechanism of absolute negative conductivity (ANC) in a two-dimensional electron system placed into magnetic and ac electric fields, which is associated with two-dimensional electron scattering by impurities, accompanied by ac field photon absorption, was proposed in  . The state with negative conductivity is unstable, the system decays into domains, and the measured macroscopic resistance becomes zero. The existence of this effect was experimentally confirmed in 2002 .
10. Non-Local Magnetic Field Evolution in Plasma
We consider the 3D non-stationary magnetic field evolution in plasma. With this aim let us transform Equation (8.19) written as follows
Let us apply the vector product operator to the left and right sides of (10.1)
Let us remind that in the coordinate notations we have
We use e as absolute charge of electron, as a mean velocity of the electron motion and the operator relation
The corresponding current density is
We introduce also the charge number density
and the character relaxation time
We transform (10.2):
1) Keeping only the terms proportional to the first power of the non-local parameter ,
2) Using the relations (for simplicity) , ; the following transformation should be taken into account
The last line in (10.8) corresponds the physical system in which the charge current is realized by electrons. Extremely important that the last term in (10.8) contains the hydrodynamic velocity , which can be defined only as a result of the whole hydrodynamic problem.
Then if we want to separate the electro-dynamic and hydro-dynamic problems we should neglect the changes in the space of the hydrodynamic velocity fluctuations. Therefore we use the relation
Using also the Maxwell relation
we write (10.2) as follows (omitting the gravitation influence)
Keeping as before only the terms proportional to the first power of the non-local parameter , we can use the Maxwell relation
Let us use the identity
and the Maxwell relation
or, if the reverse relaxation time is
or after integration on time we reach
11. Quantization in the Theory of the Magnetic Field Penetration. Comparison of the Non-Local and London’s Theories
Let us consider the 1D non-stationary electron evolution in plasma. We suppose that
After using (11.1) we have from (10.23)
In the following and are character time and length relaxation. Let us remind the London’s equation written here as
where is the London penetration depth.
This London equation predicts that the magnetic field in a superconductor decays exponentially decays from whatever value it possesses at the surface.
We choose the function in the form
is telegraph equation which contains now and as calculations parameters.
Solution of Equation (11.8) can be written in the form
where is the character scale of the relaxation length.
Solution (11.9) responds to natural initial and boundary conditions, namely
Solution (11.9) can satisfy the Equation (11.8) if a condition should be satisfied. This condition can be found after the direct substitution of (11.9) in (11.8). We have
The relation (11.14) is fulfilled if the expression in the curly bracket is equal to zero.
As we see in the stationary case we have known equation
which has a solution
Finally, we have non-stationary solution as a wave with the time attenuation.
and approximately we reach
1) Relation (11.16) defines penetration of the longitudinal magnetic field in space.
2) This non-stationary attenuation has the character of wave damping.
3) The stationary case corresponds to the London’s regime.
4) It is interesting to notice that the regime of the non-stationary wave attenuation exists for the evolution of the spherical object without fields (see Appendix 3).
12. Nonlocal Pointing-Umov Theorem
Let us obtain now the generalized nonlocal formulation of the Pointing-Umov theorem. In other words we intend to obtain the law of the energy conservation for electro-magnetic processes in the frame of non-local physics. We use the Equation (8.19) for following transformations.
We use the vector identity
Let us transform now the first term of the right-hand-side of (12.3)
After integration on time, we reach
where a function is defined by the initial and boundary conditions. Scalar multiplication by of the both parts of Equation (12.5) gives the relation
Let us use the vector identity
or (for simplicity , and )
It is the nonlocal formulation of the Umov-Pointing theorem which takes into account transport processes in physical system. In the local case we find
As usual a function is originated by integration on time and can be taken as zero; we obtain the classical formulation of the Umov-Pointing theorem.
This equation is well known local energy equation, where the left hand side of this equation corresponds to the Joule heating.
13. The Soliton Movement in One Species Physical System about the Heat Transfer in Graphene
Particular attention of researchers has been recently attracted to a new carbon material, i.e., graphene, consisting of a single layer of carbon atoms and having a planar hexagonal structure. Electromagnetic waves propagating in carbon structures become highly nonlinear even at relatively weak fields, which results in possible propagation of electromagnetic solitary waves (which are soliton analogs, or even solitons) in carbon nanotubes and graphene (see review in    ). The discussed properties of carbon structures have generated both increased theoretical interest and attempts at application in nonlinear optical devices.
The fundamental result also consists in discovering of the extremely high thermal conductivity in two-dimensional crystals including graphene. This effect (which takes place even in room-temperature) can be explained as a result of soliton movements without destruction. This effect was forecasted by me many years ago. Really, the transport processes in graphene as the effect of the soliton movement are investigated by me many years ago. The fundamental monograph  contains Chapter 6 (Quantum Solitons in Solid Matter), Item 6.2 (Application of non-local quantum hydrodynamics to the description of the charge density waves in the graphene crystal lattice) with the fundamental conclusion (p. 178): “Important conclusion: high temperature superconductors demonstrate new type of electronic order and modulation of atomic positions. The above mentioned graphene properties can be explained only in the frame of the self-consistent non-local quantum theory which leads to the appearance of the soliton waves moving in graphene.”
Then it is impossible to talk about individual particles as heating transfer carriers. In the definite sense the high temperature super conductivity and super heat conductivity have the same origin—the appearance of moving solitons without destruction.
Let us demonstrate the example of the electron soliton movement in 1D physical system.
We transform nonlocal Maxwell equation
where is the electron charge density, or
where is electrical potential and is the electron number density. Equation (13.2) can be written in terms of wave parameters ( ). Namely
where is the mass density, or introducing the absolute electron charge and potential we find
We use the following system of scales:
, , , , , (13.5)
The single independent scale is .
Taking into account , , we have
, , . (13.7)
The considered physical system works in Meissner regime, and then we needn’t to use the influence of magnetic field. We reach the system of dimensionless equations
Appendix 2 contains the corresponding Maple program.
The following Figure 13 reflects the result of calculations for the Cauchy conditions:
As we see from Figure 13: the pressure distribution (leading to the heat transport) and the density distribution (leading to the charge transport) have the character of moving solitons.
Figure 13. Density r ( )—solid line, velocity u ( )—dashed line, pressure p ( )—dash-dot line, self-consistent potential v ( )—dotted line.
14. Discussion, Principal Derivations, Conclusion and Proposals
Research of superconductors is curried out very actively. But in spite of obvious success, the following conclusion could be established:
1) Contemporary theories of superconductivity based on the Schrödinger equation, practically exhaust their arguments and have no possibility to explain effects of the high temperature superconductivity.
2) Contemporary theories of superconductivity (including BCS) based on the Schrödinger equation, can’t propose the principles of search and creation of superconducting materials.
The most impressive demonstration of these difficulties consists in the fundamental distinction between a strange metal and a conventional metal; in other words in the absence of well-defined quasi-particles in the frame of local physics. This is manifested in transport properties which defy conventional theory, the most famous of which is a T-linear resistivity that persists from nearly 0 K to high temperatures above the proposed Mott-Ioe-Regel (MIR) limit, beyond which Boltzmann theory ceases to be valid.
Without exception, all existing before proposals fail this test. Most of these theories depart from the assumption that the electrical currents are carried by one or the other system of quasi-particles. This is fundamental: it is impossible to identify the simplicity principle dealing with particle physics. The transport is assumed to be due to thermally excited quasi-particles behaving like classical balls being scattered in various ways, dumping eventually their momentum in the lattice. However, it is a matter of principle that the physics of such quasi-particles in real solids is never simple. These interact with phonons which are very efficient sources of momentum dissipation which should be strongly temperature dependent for elementary reasons.
The conclusion is that quasi-particles are quantum solitons which are moving without destruction.
From position of the quantum non-local hydrodynamics, the problem of search and creation of superconductive materials come to the search of materials which lattices ensure the soliton movement without destruction. In my opinion, the mentioned materials can be created artificially using the technology of the special introduction of quantum dots in matrices on the basement of proposed quantum hydrodynamics. It is known that technology of material creation with special quantum dots exists now in other applications.
Maple program for modeling of soliton motion in two species physical system.
(the program is ready for application)
The ratio L/H is ratio of masses of the light and heavy particles; for example L= 1, H= 1838.
Maple program for modeling of soliton motion in one species physical system.
(the program is ready for application)
The non-stationary wave attenuation of the spherical object without fields.
Let us consider the nonlocal space—time evolution of the spherical object. We suppose:
1) The nonlocal parameter .
2) Radial velocity of the object .
3) The thermal velocity , then
The nonlocal system of hydrodynamic equation has the following form—the continuity equation is 
The system consists of the continuity equation and the energy equation ( —density, p—quantum pressure). The equation of motion in a spherical coordinate system is absent if the radial component of the velocity is absent in the case of radial symmetry. Using (A3.1) we find in the dimensionless form; continuity equation
The following scales are used: , . Nonlinear parabolic Equations (A3.4), (A3.5) can be solved by Fourier method. Really, for (A3.4) we separate the unknown variables using
Relations (A3.7) and (A3.8) lead to two ordinary differential equations
, . (A3.11)
Equation (A3.10) can be solved by numerical methods. For example
The solution for the quantum pressure can be found by the analogical way. It should be noticed that the solution of the homogeneous Equations (A3.4) and (A3.5) for unknown variables p and can be found up to an additive arbitrary constant (Figure 14).
Figure 14. Radial density damping.
 Bogolyubov, N.N. (1946) Problemy Dinamicheskoi Teorii v Statisticheskoi Fizike (Dynamic Theory Problems in Statistical Physics). Gostekhizdat, Moscow Leningrad. [Translated into English The Dynamical Theory in Statistical Physics (Hindustan Publ. Corp., Delhi, 1965)]