## 1. Introduction

One of the crises of contemporary mathematics belongs in part to the subject of the infinite and infinitesimals [1]. It originates from the barest necessity to develop a rigorous language for description of observable physical phenomena. It was a time when foundations of integral and differential calculi were developed. A theoretical foundation for facilitation of understanding of classical mechanics is provided by the concepts of absolute time and space originally formulated by Sir Isaac Newton [2]. Space is distinct from body. And time passes uniformly without regard to whether anything happens in the world. For this reason Newton spoke of absolute space and absolute time as of a "container" for all possible objects and events. The space-time container is absolutely empty until prescribed metric and a reference frame are introduced. Infinitesimals are main tools of differential calculus [3, 4] within chosen reference frames.

Infinitesimal increment being a cornerstone of theoretical physics has one receptee default belief, that increment *δV* tending to zero contains a lot of events to be under consideration. Probability of detection of a particle within this infinitesimal volume *δV* tends to zero we lose information about amount of the events. What is more, the information becomes uncertain. It means the infinitesimal increment being applied in physics faces with a conflict of depth of understanding physical processes on such minuscule scales. This trouble is avoided in quantum mechanics by proclamation that infinitesimal increments are operators, whereas observables are averaged on an ensemble.

In light of classical views Newton maintained the theory that light was made up of tiny particles, corpuscles. They spread through space in accordance with law of classical mechanics. Christian Huygens (a contemporary of Newton), believed that light was made up of waves vibrating up and down perpendicularly to direction of the light propagation. It comes into contradiction with Newtonian idea about corpuscular nature of light. Huygens was a proclaimer of wave mechanics as opposite to classical mechanics [5].

We abstain from allusion to physical vacuum but expand Huygensian idea to its logical completion. Let us imagine that all Newtonian absolute space is not empty but is populated everywhere densely by Huygensian vibrators. The vibrators are silent at absence of wave propagating through. But as soon as a wave front reaches some surface all vibrators on this surface begin to radiate at a frequency resonant with that of incident wave. From here it follows, that in each point of the space there are vibrators with different frequencies ranging from zero frequency up to infinite. All are silent in absence of an external wave perturbance. Thus, the infinitesimal volume *δV* is populated by infinite amount of the vibrators with frequencies ranging from zero to infinity. They are virtual vibrators facilitating propagation of waves through space.

Let us return to our days. One believes that besides matter and physical fields there is nothing more in the universe. Elementary particles are only a building material of "eternal and indestructible" substance of the cosmos. However we should avow that all observed matter and physical fields, are not the basis of the world, but they are only a small part of the total quantum reality. Physical vacuum in this picture constitutes a basic part of this reality. In particular, modern conception of the physical vacuum covers Huygens's idea perfectly. All space is fully populated by virtual particle-antiparticle pairs situated on a ground level. Such a particle-antiparticle pair has zero mass, zero charge, and zero magnetic moment. Famous Dirac’ sea (Dirac postulated that the entire universe is entirely filled by particles with negative energy) is boundless space of electron-positron pairs populated everywhere densely - each quantum state possessing a positive energy is accompanied by a corresponding state with negative energy. Electron has positive mass and positron has the same mass but negative; electron has negative charge and positron has the same positive charge; when electron and positron dance in pair theirs magnetic moments have opposite orientations, so magnetic moment of the pair is zero.

Let electron and positron of a virtual pair rotate about mass center of this pair. Rotation of the pairs happens on a Bohr orbit. Energy level of the first Bohr orbit, for example, *mv*^{ 2}/2, is about 14 eV. Here *m* is electron mass and *v* is its velocity (on the first Bohr orbit the velocity is about 2.188 10^{6} m/s). Energy of the pair remains zero since positron has the same energy but with negative sign. Quantum fluctuations around this zero energy are as zero oscillations of electromagnetic field. Observe that, energy releasing of electron and positron from vacuum occurs at *mc*^{2} = 1.022 MeV. So we see that there are about 7.5 10^{4} Bohr orbits lying below this energy. That is, there is a vast scope for occupation of different Bohr orbits by the virtual electron-positron pairs.

A short outline given above is a basis for understanding of interference effects to be described below.

Ones suppose that random fluctuations of electron-positron pairs take place always. What is more, these fluctuations are induced by other pairs and by particles traveling through this random conglomerate. Edward Nelson has described mathematical models [6, 7] representing the above random fluctuations as Brownian motions of particles that are subjected by random impacts from particles populating aether (Nelson's title of a lower environment). The model is viewed as a Markov process

loaded by a Wiener term *w* with diffusion coefficient equal to *ћ*/2*m*, where *m* is a mass of the particle and *ћ* is the reduced Planck constant. In this perspective Nelson has considered two Markov processes complementary to each other. One is described by forward-difference operator in the time, here *b*(*x*(*t*),*t*) is a velocity calculated forward. And other equation is described by backward-difference operator in the time with the velocity *b*^{†}(*x*(*t*),*t*) calculated backward. In general *b*(*x*(*t*),*t*) ≠ *b*^{†}(*x*(*t*),*t*) [6, 7]. The two complementary processes, by means of transition to two new variables, real and imaginary, finally lead to emergence of the Schrödinger equation.

Nelson' vision that aether fluctuations look as Brownian movements of subparticles with *ћ*/2*m* being the diffusion coefficient of the movements, correlates with Feynman's ideas about quantum fluctuations of virtual particles in vacuum [8]. The Feynman path integral is akin to the Einstein-Smoluchowski integral equation [9]. The latter computes transition probability density. We shall deal with the modified Feynman path integral loaded by a temperature multiplied by the Boltzmann constant. At that, probability amplitude stays as a fundamental mathematical object at all stages of computations. As an example we shall consider emergence of interference patterns at scattering heavy particles on gratings. The particles are heavy in the sense that they adjoin to both realms, quantum and classical. They are nanoparticles. Such nanoparticles have masses about 100 amu and more, as, for example, fullerene molecules [10] shown in Fig. 1. It is remarkable that there are many experiments with such molecules showing interference patterns in the near field [11-18]. On the other hand these molecules are so large, that they behave themselves as classical particles at ordinary conditions.

The article consists of five sections. Sec. 2 introduces a general conception of the path integral that describes transitions along paths both of classical and quantum particles. Here we fulfill expansions in the Taylor series of terms presented in the path integral. Depending on type of presented parameters we disclose either the Schrödinger equation or diffusion-drift equation containing extra term, osmotic diffusion. In the end of the section we compute passing nanoparticles through *N*-slit grating. Sec. 3 deals with interference patterns from the *N*-slit grating. Specifically, we study blurring of the Talbot carpet (an interference pattern emergent under special conditions imposed on the grating [19, 20]) arising under decoherence of incident on the grating nanoparticles. In Sec. 4 we find equations for computing the Bohmian trajectories. Also we compute variance of momenta along the trajectories. These computations lead to emergence of the uncertainty conditions. In concluding Sec. 5 we summarize results. For confirmation of existence of the Bohmian trajectories here we mention interference experiments with single silicon oil droplet [21].

## 2. Generalized path integral

Let many classical particles occupy a volume *V* and they move with different velocities in different directions. Let us imagine that there is a predominant orientation along which ensemble of the particles drifts. As a rule, one chooses a small volume *δV* in order to evaluate such a drift, Fig. 2. Learning of statistical mechanics begins with assumption that the volume should contain many particles.

The problem is to find transition probabilities that describe transition of the particle ensemble from one statistical state to another. These transient probabilities can be found through solution of the integral Einstein-Smoluchowski equation [9]. This equation in mathematical physics is known as Chapman-Kolmogorov [22-24] equation.. This equation looks as follows

This equation describes a Markovian process without memory. That is, only previous state bears information for the next state. Integration here is fulfilled over all a working scene R, encompassing finite space volume. Infinitesimal volume *δV* in the integral tends to zero. It should be noted, however, that this volume should contain as many particles as possible for getting a satisfactory statistical pattern. One can see that this claim enters in conflict with the assumption *δV* → 0.

Next we shall slightly modify approach to this problem. Essential difference from the classical probability theory is that instead of the probabilities we shall deal with probability amplitudes. The transition amplitudes can contain also imaginary terms. They bear information about phase shifts accumulated along paths. In that way, a transition from an initial state ^{3} (see Fig. 2) is represented by the following path integral

Here function *ψ*(...) is a probability amplitude. Probability density *p*(...), in turn, is represented by square of modulo of the probability amplitude, namely, *p*(...)=|*ψ*(...)|^{2}. Integral kernel

where denominator Γ under exponent is a complex-valued quantity, i.e., Γ = β+**i***ћ*. The both parameters, β and the reduced Planck constant *ћ*, have dimensionality of energy multiplied by time. From here it follows that β = 2*k*_{B}*Tδt*. Here *k*_{B} is Boltzmann constant and *T* is temperature. So, we can write

Factor 2 at the first term is conditioned by the fact that the kernel *K* relates to transitions of the probability amplitude *ψ*, not the probability *p*. Observe that a fullerene molecular beam for interference experiment prepared in a Knudsen cell at T=1070 K [27] spreads further in a vacuum. That is, fullerene molecules keep this temperature. From here it follows, that thermal fluctuations of carbon atoms from equilibrium occur at that temperature as the fullerene molecules propagate further within the vacuum chamber. One can see that at *T*=1000 K the term 2*k*_{B}*Tδt* may be about *ћ* if *δt* is about 10^{-14} s.

Next let us imagine that particles pass through a path length one by one. That is, they do not collide with each other along the path length. The particles are complex objects, however. They are nanoparticles. Fullerene molecule, for example, contains 60 carbon atoms, Fig. 1. Conditionally we can think that the atoms are connected with each other by springs simulating elastic vibrations. In this view the Lagrangian

Here *N* is amount of atoms, constituent complex molecule, and *m* is mass of a single atom. By supposing that there are no quantum permutations between atoms we may expand the Lagrangian further

(7) |

Here we admit that *m*_{N} = *Nm* is a mass of the molecule referred to its center of mass. Small deviations *k*th atoms with respect to the center of mass.

Let us consider the second row in Eq. (7). First of all we note that the term

The path integral (3) contains functions depending only on coordinates *ћ* the parameter 2*k*_{B}*Tδt* in Eq. (5) can be different from zero as well.

### 2.1. Expansion of the path integral

The next step is to expand terms, ingoing into the integral (3), into Taylor series. The wave function written on the left is expanded up to the first term

As for the terms under the integral, here we preliminarily make some transformations. We define a small increment

The Lagrangian (7) is rewritten, in such a case, in the following view

Here *K* contains the Lagrangian given from Eq. (10). The under integral function

The potential energy

Taking into account the expressions (8)-(11) and substituting theirs into Eq. (3) we get

(12) |

First, we consider terms enveloped by braces (a) and (b): (a) here displacement *δt* is a velocity *th* atom at its deviation from a steady position relative to the center of mass. Summation through all deviations of atoms divided by *N* gives averaged velocity, *T*_{N}, of the atoms oscillating about the center of mass. Observe that energy of thermal fluctuations, *T*_{N}, is proportional to *k*_{B}*T*. Because of its presence in the propagator intensity of an interference pattern diminishes in general. Further we shall add this term into the potential energy as some constant component.

In the light of the above observation we rewrite Eq. (12) as follows

(13) |

Here we have expanded preliminarily exponents to the Taylor series up to the first term of the expansion. Exception relates to the term exp{ -*m*_{N}
*ξ*^{ 2}/2Γ*δt*} which remains in its original form. This exponent integrated over all space R ^{3} results in

To derive outcomes of integration of terms containing factors

In the light of this observation let us now solve integral (13) accurate to terms containing *δt* not higher the first order:

We note that the term

The parameter

is seen to be as a complex-valued diffusion coefficient consisting of real and imaginary parts.

#### 2.1.1. Temperature is zero

Let *k*_{B}*Tδt* = 0. It means that Γ = i*ћ*. Also *T*_{N} = 0. One can see that Eq. (17) is reduced to

It is the Schrödinger equation.

#### 2.1.2. Temperature is not zero

Let *k*_{B}*Tδt* >> *ћ*. We can suppose that Γ = *k*_{B}*Tδt*. Eq. (17) takes a form

Here

is the diffusion coefficient. The coefficient has dimensionality of [length^{2} time^{-1}]. It is a factor of proportionality representing amount of substance diffusing across a unit area in a unit time - concentration gradient in unit time.

We can see that Eq. (20) deals with the amplitude function *ψ*, not a concentration. However, a measurable function is *ρ*=|*ψ*|^{2} - concentration having dimensionality of [(amount of substance) length^{-3}]. Let us multiply Eq. (20) from the left by 2*ψ*. First we note that the combination

Extra term

is an osmotic velocity required of the particle to counteract osmotic effects [6]. Namely, imagine a suspension of many Brownian particles within a physical volume acted on by an external, virtual in general, force. This force is balanced by an osmotic pressure force of the suspension [6]:

From here it is seen, that the osmotic pressure force arises always when density difference exists and especially when the density tends to zero. And vice-versa, the force disappears in extra-dense media with spatially homogeneous distribution of particles. As states the second law of thermodynamics spontaneous processes happen with increasing entropy. The osmosis evolves spontaneously because it leads to increase of disorder, i.e., with increase of entropy. When the entropy gradient becomes zero the system achieves equilibrium, osmotic forces vanish.

Due to appearance of the term

Reduction to PDEs, Eq. (19) and Eq. (22), was done with aim to show that the both quantum and classical realms adjoin with each other much more closely, than it could seem with the first glance. Further we shall return to the integral path paradigm [8] and calculate patterns arising after passing particles through gratings. We shall combine quantum and classical realms by introducing the complex-valued parameter Γ = 2*k*_{B}*Tδt* + i*ћ*.

### 2.2. Paths through *N*-slit grating

Computation of a passing particle through a grating is based on the path integral technique [8]. We begin with writing the path integral that describes passing the particle through a slit made in an opaque screen that is situated perpendicularly to axis *z*, Fig. 3. For this reason we need to describe a movement of the particle from a source to the screen and its possible deflection at passing through the slit, see Fig. 4. At that we need to evaluate all possible deflections.

We believe, that before the screen and after it, the particle (fullerene molecule, for example) moves as a free particle. Its Lagrangian, rewritten from Eq. (17), describes its deflection from a straight path in the following form

The first term relates to movement of the center of mass of the molecule. So that *m*_{N} is mass of the molecule and *z*. The second term is conditioned by collective fluctuations of atoms constituent this molecule. This term is nonzero when atoms have predominant fluctuations along axis *x*. For the sake of simplicity we admit that *p*_{z}, is much greater than its transverse component [16, 17, 29] and we believe it is constant also. By translating a particle's position on a small distance *δx* = (*x*_{b} - *x*_{a})<< 1 for a small time *δt* = (*t*_{b} - *t*_{a})<< 1 we find that a weight factor of such a translation has the following form

The particles flying to the grating slit along a ray α, Fig. 4, pass through the slit within a range from *x*_{0}-*b*_{0} to *x*_{0}+*b*_{0} with high probability. The path integral in that case reads

Integral kernel (propagator) for the particle freely flying is as follows [8, 26]

Here *v*_{(a,b)} = (*x*_{b}
*– x*_{a})/(*t*_{b}
*– t*_{b}) is a velocity of the molecule on a segment from *x*_{a} to *x*_{b}. And *v* is an average velocity of collective deflection of atoms constituent the molecule. It was defined in Eq. (12). We shall believe that the ratio *v*/*v*_{ (a,b)} is small enough. We can define a new renormalized mass {*m*}_{N} = *m*_{N} (1 - *v*/*v*_{(a,b)}) and further we shall deal with this mass.

#### 2.2.1. Passing of a particle through slit

By substituting the kernel (28) into the integral (27) we obtain the following detailed form

(29) |

The integral is computed within a finite interval [-*b*_{0},+*b*_{0}]. Observe, that the integrating can be broadened from -∞ to +∞. But in this case we need to load the integral by the step function equal to unit within the finite interval [-*b*_{0},+*b*_{0}] and it vanishes outside of the interval. The step function, that simulate a single slit, can be approximated by the following a set of the Gaussian functions [30]

Here parameter *b* is a half-width of the slit, real *η >* 0 is a tuning parameter, and *K* takes integer values. At *K* = 1 this form-factor degenerates to a single Gaussian function. And at *K* → ∞ this function tends to an infinite collection of the Kronecker deltas which fill everywhere densely the step function. We rewrite Eq. (29) with inserting this form-factor

We do not write parameters *η* and *K* in the Gaussian form-factor and for the sake of simplicity further we shall consider they equal to 1. That is, for simulating the slit we select a single Gaussian function.

#### 2.2.2. Definition of new working parameters

First we replace the flight times *τ*_{0} and *τ*_{1} by flight distances (*z*_{0}-*z*_{s}) and (*z*_{1}-*z*_{0}), see Fig. 4. This replacement reads

where *v*_{z} is a particle velocity along the axis *z*.

There is, however, one more parameter of time that is represented in definition of the coefficient Γ= 2*k*_{B}*Tδt* + **i***ћ*. It is a small time increment *δt*. The parameter *δt* first appears in the path integral (3) as the time increment along a path. In accordance with the uncertainty principle *δt* should be more or equal to the ratio of *ћ*, Planck constant, to energy of occurring events. In a case of a flying particle through vacuum it can be minimal energy of vacuum fluctuations (it is about energy of the first Bohr orbit of electron-positron pair that is about 14 eV). Evaluation gives *δt* ~ 2.8 10^{-16} s. From here it follows, that 2*k*_{B}*Tδt* is less than *ћ* on about one order at *T* = 1000 K (almost temperature of fullerene evaporation from the Knudsen cell [27]).

Emergence of the term 2*k*_{B}*Tδt* can be induced by existence of quantum drag [31] owing to different conditions for quantum fluctuations both inside of the fullerene molecule and outside what can induce weak Casimir forces. Because of the weak Casimir force the quantum drag does not lead to decoherence at least in the near zone. However further we shall see that a weak washing out of the Talbot interference pattern is due to existence of this term.

Let us divide the parameter Γ by {*m*}_{N}
*v*_{z}

Here *p*_{z} = {*m*}_{N}
*v*_{z} is a particle momentum along axis *z*. We can define the de Broglie wavelength *λ*_{dB} = *h*/*p*_{z} where *h* = 2*πћ* is the Planck constant. Let us also define a length *δ*_{T} = 4*πk*_{B}*Tδt*/({*m*}_{N}
*v*_{z}). In this view we can rewrite Eq. (33) as follows

The length *δ*_{T} tends to zero as *T* → 0. At *T* = 1000 K and at adopted *δt* = 2.8 10^{-16} s we have *δ*_{T} ≈ 0.4 pm. On the other hand, the de Broglie wavelength, *λ*_{dB}, evaluated for the fullerene molecule moving with the velocity *v*_{z} =100 m/s is about 5 pm [18]. So, we can see that the length *δ*_{T} is less of the de Broglie wavelength on about one order and smaller. A signification of the length *δ*_{T} is that it determines decoherence of a particle beam. Decoherence of flying particles occurs the quickly, the larger *δ*_{T}. Observe that the length *δ*_{T} has a close relation with the coherence width - a main parameter in the generalized Gaussian Schell-model [32, 33].

## 3. Wave function behind the grating

Wave function from one slit after integration over *ξ*_{0} from -∞ to +∞ has the following view [28]

Here argument of ψ-function contains apart *x* also *z* in order to emphasize that we carry out observation in the point (*z*,*x*), see Fig. 4. The factor (2/*π*)^{1/2} comes from (30). Parameters Ξ_{0} and Σ_{0} read

In order to simplify records here we omit subscript 1 at *x* and *z* - an observation point that is situated after the slit.

Let us consider that an opaque screen contains *N*_{0} slits spaced through equal distance, *d*, from each other. Numeration of the slits is given as it is shown in Fig. 3, *n*_{0} = 0,1,2, …, *N*_{0} – 1. Sum of all wave functions (35), each of which calculates outcome from an individual slit, gives a total effect in the point (*z*,*x*) where a detector is placed:

Probability density in the vicinity of the observation point (*x,z*) reads

Calculation of the wave function (37) is fulfilled for the grating containing *N*_{0} = 32 slits. Distance between slits is *d* = 10^{5}*λ*_{dB}. So at *λ*_{dB} = 5 pm the distance is equal to 500 nm. Requirement *λ*_{dB} << *d* and *N*_{0} tending to infinity together with a condition that the particle beam is paraxial, that is, *x*_{s} = 0 and *z*_{s} → -∞, provides emergence in the near-field of an interference pattern, named Talbot carpet [19,20]. Here a spacing along interference patterns is measured in the Talbot length

which is a convenient natural length at representation of interference patterns. Since we restrict themselves by finite *N*_{0} we have a defective carpet, which progressively collapses as a spacing from the slit increases. Fig. 5 shows the Talbot carpet, being perfect in the vicinity of the grating slit; it is destroyed progressively with increasing *z*_{T}. As for the Talbot carpet we have the following observation. We see that in a cross-section *z* = *z*_{T}/2 image reproduces radiation of the slits but phase-shifted by half period between them. At *z* = *z*_{T} radiation of the slits is reproduced again on the same positions where the slits are placed. And so forth.

Evaluation of sizes of the interference pattern is given by ratio of the Talbot length to a length of the slit grating. In our case the Talbot length is *z*_{T} = 0.1 m. And length of the slit grating is about *N*_{0}*d* = 1.6 10^{-5} m. From here we find that the ratio is 6250. It means that the interference pattern shown in Fig. 5 represents itself a very narrow strip.

Zigzag curves, drawn in the upper part of Fig. 5 by dark blue color, show Bohmian trajectories that start from the slit No. 15. One can see that particles prefer to move between nodes having positive interference effect and avoid empty lacunas. However the above we noted, that the ratio of the Talbot length to the length of the grating is about 6250 >> 1. It means that really the Bohmian trajectories look almost as straight lines slightly divergent apart. Zigzag-like behavior of the trajectories is almost invisible. Such an almost feebly marked zigzag-like behavior may be induced by fluctuations of virtual particles escorting the real particle.

As soon as we add the term *k*_{B}*Tδt* different from zero (*T* > 0 K) we observe blurring the interference pattern. The blurring is the stronger, the larger *k*_{B}*Tδt*. For comparison see Figs. 6 and 7. Here instead of *k*_{B}*Tδt* we write a more convenient parameter, the coherence length *δ*_{T}. This length characterizes a dispersed divergence from initially tuned the de Broglie wavelength. Such a disperse medium can be due to quantum drag on the vacuum fluctuations. Here the Bohmian trajectories are not shown, since because of the thermal term *k*_{B}*Tδt* > 0 a Brownian like scattering of the trajectories arises. This scattering we shall discuss later on.

We may think that the technical vacuum can be not perfect. It causes additional scattering of particles on residual gases. Because of this additional scattering the interference pattern can be destroyed entirely, as shown in Fig. 8.

Now let us draw dependence of the probability density *p*(*x, z*) as a function of *x* at fixed *z*. In other words, we calculate interference fringes in a cross-section of the interference patterns at *z* = *z*_{T}/2 for different values of the length *δ*_{T}. Such a cross-section is chosen because a self image of the slit grating appears phase-shifted by half period of the grating. For that reason we should see the interference fringes spaced between the slit sources of radiation.

Fig. 9 shows three characteristic patterns of the interference fringes. In Fig. 9(a) almost ideal interference fringes are shown obtained at *δ*_{T} < *λ*_{dB}. Fig. 9(b) shows interference fringes obtained at *δ*_{T} ~ *λ*_{dB}. It is instructive to compare these interference fringes with those that have been measured in experiments [14, 34]. And Fig. 9(c) shows disappearance of interference fringes because of strong scattering of the particles on residual gases in vacuum, *δ*_{T} > *λ*_{dB}.

Disappearance of interference fringes is numerically evaluated by calculating a characteristic called visibility [14, 34]. The fringe visibility [27] is represented as a ratio of difference between maximal and minimal intensities of the fringes to their sum:

Evaluation of *P*_{max} and *P*_{min} is shown in Fig. 9(b). As follows from the figure, the evaluations are fulfilled in a central region of the grating. That is, edges of the grating have to be left far off from the measured zone. The visibility *V* as a function of the parameter *δ*_{T} is shown in Fig. 10. One can see that crossover from almost perfect interference fringes, *V* = 1, up to their absence, *V* = 0, begins near *δ*_{T}. ~ *λ*_{dB}. Transition from almost coherent particle beam to incoherent is a cause of such a crossover [30].

## 4. Bohmian trajectories and variance of momenta and positions along paths

Here we repeat computations of David Bohm [35] which lead to the Hamiltoton-Jacobi equation loaded by the quantum potential and, as consequence, to finding Bohmian trajectories. But instead of the Schrödinger equation we choose Eq. (17) that contains complex-valued parameter Γ=β+i*ћ*:

Here *β* = 2*k*_{B}*Tδt* (in particular, the diffusion coefficient reads *D* = *β*/2*m*_{N}) and

Further we apply polar representation of the wave function, *ψ* = *R* exp{i*S*/*ћ*}. It leads to obtaining two equations for real and imaginary parts that deal with real-valued functions *R* and *S*. The function *R* is the amplitude of the wave function and *S*/*ћ* is its phase. After series of computations, aim of which is to put together real and imaginary terms, we obtain the following equations

Firstly, we can see that at *β* = 0 Eq. (42) reduces to the modified Hamilton-Jacobi equation due to loaded the quantum potential that is enveloped here by brace (b). And Eq. (43) reduces to the continuity equation. These equations read

Terms enveloped by braces (a), (b), and (c) are the kinetic energy of the particle, the quantum potential, and the right part is a kernel of the continuity equation (45), respectively. In particular, the term

Momentum of the particle reads

where *p* to the wavelength *λ*_{dB} = *h/p*, where *h =* 2*πћ* is the Planck constant. Now, as soon as we found the current velocity

positions of the particle in each current time beginning from the grating' slits up to a detector is calculated by the following formula

Here *t* is a current time that starts from *t*=0 on a slit source and *δt* is an arbitrarily small increment of time. Some calculated trajectories of particles, the Bohmian trajectories, are shown in upper part of Fig. 5. It should be noted that the Bohmian trajectories follow from exact solutions of Eqs. (44)-(45). These equations give a rule for finding geodesic trajectories and secants of equal phases, *S/ћ*, at given boundary conditions. The geodesic trajectories point to tendency of the particle migration along paths. And the secant surfaces describe a coherence of all the passing particles created on a single source.

In case of *β* >> *ћ* we have the following two equations

Here we take into consideration that in the first equation we may replace *β*/2*m*_{N} =*D* is the diffusion coefficient. Now we may rewrite Eqs. (49)-(50) as follows

Here

### 4.1. Dispersion of trajectories and the uncertainty principle

As for the Bohmian trajectories there is a problem concerning their possible existence. As follows from Eqs. (46) and (48) in each moment of time there are definite values of the momentum and the coordinate of a particle moving along the Bohmian trajectory. This statement enters in conflict with the uncertainty principle.

Here we try to retrace emergence of the uncertainty principle stemming from standard probability-theoretical computations of expectation value and variance of a particle momentum. We adopt a wave function in the polar representation

where *R* = (*ρ*)^{1/2} is the amplitude of the wave function (*ρ* =*R*^{2} = <*Ψ*|*Ψ*> is the probability density) and *S/ћ* is its phase. Momentum operator *m* is mass of the particle. Expectation value of the velocity operator reads

The velocity V _{g} is seen to be complex-valued. Here *ћ/*2*m*) is the quantum diffusion coefficient [6, 7]. Therefore its imaginary part is a quantum osmotic velocity

It is instructive to compare this velocity with the classical osmotic velocity given in Eq. (23). As can see the osmotic velocity stems from gradient of entropy that evaluates degree of order and disorder on a quantum level, likely of vacuum fluctuations.

Real part of Eq. (54) gives the current velocity

Let us now calculate variance of the velocity V _{g}. This computation reads

Terms over bracket (d) kill each other as follows from Eq. (54). It is reasonable in the perspective to multiply Var(V _{g}) by *m*/2

So, this expression has a dimensionality of energy. The first term to be computed represents the following result

Here the term enveloped by bracket (a) is a kinetic energy of the particle, the term enveloped by bracket (b) with negative sign added is the quantum potential *Q*, and the term enveloped by bracket (c) comes from the continuity equation. See for comparison Eqs. (44) and (45). We rewrite the quantum potential as follows

As for the second term in Eq. (57) we have

One can see that the variance consists of real and imaginary parts. Observe that the right side is represented through square of gradient of the complexified action [28], namely

The first term in this expression represents kinetic energy, *E*, of the particle. The second term, stemming from the quantum potential, contains under braces a term having dimensionality of inverse time, that is, of frequency

This frequency multiplied by *ћ* represents an energy binding a particle with vacuum fluctuations. This energy, as follows from Eq. (61), is equal to the particle mass multiplied by squared the osmotic velocity (55) and divided by 2. It is an osmotic kinetic energy. In the light of the above said we rewrite Eq. (61) in the following way

Let we have two Bohmian trajectories. Along one trajectory we have *E*_{1} - *ћω*_{Q,1}, and along other trajectory we have a perturbed value *E*_{2} - *ћω*_{Q,2}. Subtracting one from other we have

One can suppose that emergence of the second trajectory was conditioned by a perturbation of the particle moving along the first trajectory. If it is so, then emergence of the second trajectory stems from an operation of measurement of some parameters of the particle. One can think that duration of the measurement is about *δt* = 1/*δω*_{Q}. From here we find

Now let us return to Eq. (48) and rewrite it in the following view

The initial Bohmian trajectory is marked here by subscript 1. Observe that *δE* into Eq. (66) we obtain finally

Here we take into account

## 5. Concluding remarks

Each nanoparticle incident on a slit grating passes only through a single slit. Its path runs along a Bohmian trajectory which is represented as an optimal path for the nanoparticle migrating from a source to a detector. Unfortunately, the Bohmian trajectory can not be observable since a serious obstacle for the observation comes from the uncertainty principle. In other words, an attempt to measure any attribute of the nanoparticle, be it position or orientation, i.e., the particle momentum, leads to destroying information relating to future history of the nanoparticle. What is more, any collision of the nanoparticle with a foreign particle destroys the Bohmian trajectory which could give a real contribution to the interference pattern. It relates closely with quality of vacuum. In the case of a bad vacuum such collisions will occur frequently. They lead to destruction of the Bohmian trajectories. Actually, they degenerate to Brownian trajectories.

Excellent article [21] of Couder & Fort with droplets gives, however, a clear hint of what happens when the nanoparticle passes through a single grating slit. In the light of this hint we may admit that the particle “bouncing at moving through vacuum” generates a wave at each bounce. So, a holistic quantum mechanical object is the particle + wave. Here the wave to be generated by the particle plays a role of the pilot-wave first formulated by Lui de Broglie and later developed by Bohm [37]. It is interesting to note in this context, that the pilot-waves have many common with Huygens waves [5].

A particle passing through vacuum generates waves with wavelength that is inversely proportional to its momentum (it follows from the de Broglie formula, *λ* = *h/p*, where *h* is Planck's constant). One can guess that a role of the vacuum in the experiment of Couder & Fort takes upon itself a silicon oil surface with subcritical Faraday ripples activated on it [21]. Observe that pattern of the ripples is changed in the vicinity of extraneous bodies immersed in the oil which simulate grating slits. Interference of the ripples with waves generated by the bouncing droplets provides optimal paths for the droplets traveling through the slits and further. As a result we may observe an interference pattern emergent depending on amount of slits in the grating and distance between them.

Now we may suppose that the subcritical Faraday ripples on the silicon oil surface simulate vacuum fluctuations. Consequently, the vacuum fluctuations change their own pattern near the slit grating depending on amount of slits and distance between them. We may imagine that the particle passing through vacuum (bouncing through, Fig. 11) initiates waves which interfere with the vacuum fluctuations. As a result of such an interference the particle moves along an optimal path - along the Bohmian trajectory. Mathematically the bounce is imitated by an exponential term exp{i*S/ћ*}, where the angle *S/ћ* parametrizes the group of rotation given on a circle of unit radius. So, the path along which the particle moves is scaled by this unitary group, U(1), due to the exponential mapping of the phase *S/ћ* on the circle.

In conclusion it would be like to remember remarkable reflection of Paul Dirac. In 1933 Paul Dirac drew attention to a special role of the action S in quantum mechanics [38] - it can exhibit itself in expressions through exp{iS/ћ}. In 1945 he emphasized once again, that the classical and quantum mechanics have many general points of crossing [39]. In particular, he had written in this article: "We can use the formal probability to set up a quantum picture rather close to the classical picture in which the coordinates q of a dynamical system have definite values at any time. We take a number of times t_{1}, t_{2}, t_{3}, … following closely one after another and set up the formal probability for the q's at each of these times lying within specified small ranges, this being permissible since the q's at any time all commutate. We then get a formal probability for the trajectory of the system in quantum mechanics lying within certain limits. This enables us to speak of some trajectories being improbable and others being likely."