SciELO - Scientific Electronic Library Online

 
vol.64 issue2Integrals of the motion and Green function for dual damped oscillators and coupled harmonic oscillatorsNewton’s law of cooling with fractional conformable derivative author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista mexicana de física

Print version ISSN 0035-001X

Rev. mex. fis. vol.64 n.2 México Mar./Apr. 2018

https://doi.org/10.31349/revmexfis.64.158 

Research

A relativistic formulation of the de la Peña-Cetto stochastic quantum mechanics

A. González Lezcanoa 

A. Cabo Montes de Ocab 

a Departamento de Física, Universidad de Pinar del Río, Pinar del Río, Cuba.

b Departamento de Física Teórica, Instituto de Cibernética, Matemática y Física, Calle E, No. 309, e/ 13 y 15, Vedado, La Habana, Cuba.


Abstract

A covariant generalization of a non-relativistic stochastic quantum mechanics introduced by de la Peña and Cetto is formulated. The analysis is done in space-time and avoids the use of a non-covariant time evolution parameter in order to search for Lorentz invariance. The covariant form of the set of iterative equations for the joint coordinate and momentum distribution function Q(x, p) is derived and expanded in power series of the coupling of the particle with the stochastic forces. Then, particular solutions of the zeroth order in the charge of the iterative equations for Q(x, p) are considered. For them, it follows that the space-time probability density ρ(x) and the function S(x) which gradient defines the mean value of the momentum at the space time point x, define a complex function ψ(x) which exactly satisfies the Klein-Gordon (KG) equation. These results for the zeroth order solution reproduce the ones formerly and independently derived in the literature. It is also argued that when the KG solution is either of positive or negative energy, the total number of particles conserves in the random motion. Other solutions for the joint distribution function in lowest order, satisfying the positive condition are also presented here. The are consistent with the assumed lack of stochastic forces implied by the zeroth order equations. It is also argued that such joint distributions, after considering the action of the stochastic forces, might furnish an explanation of the quantum mechanical properties, as associated to ensembles of particles in which the vacuum makes such particles behave in a similar way as Couder’s droplets moving over oscillating liquid surfaces. Some remarks on the solutions of the positive joint distribution problem proposed in the Olavos’s analysis are also presented.

Keywords: Stochastic QED; Couder’s experiments

PACS: 02.50.-r; 03.50.De; 03.65.-w

1. Introduction

The search for stochastic descriptions of quantum mechanics and quantum fields theory has a large history. By example, in the works 1-6, it was considered that the random forces which determine the phase space density of particle, were given in a stochastic and relativistic invariant way. This defined the so called Stochastic Electrodynamics (SQED). Specifically, these forces were generated by an electromagnetic field configuration, obtained by exciting all the oscillation modes with one half a quanta of energy. In addition, the phases of the photon modes were assumed to be stochastically defined within the interval (0.2π) in a uniform way. This force had been argued to be invariant under Lorentz transformations 2,5. In addition, in reference 1, it was argued that the stochastic motion guided by those forces, when taken in the non-relativistic limit, leads to the Schrodinger equation. Further, in the series of works 8-13, it was proposed a solution for the a central open issue of the theory: the apparent lack of positive definiteness of the coordinate-momentum joint distribution function emerging in the previous discussions. More recently, and with the purpose of start considering a generalization of the results in Ref. 1, in Ref. 18 it was suggested that random forces, showing the same statistical distribution in all Lorentz frames, can be expected to imply the satisfaction of the Klein-Gordon equation, for which the non-relativistic limit is the Schrodinger one.

The present work is devoted to present a derivation of a relativistic version of the stochastic electrodynamics. That is, we intend to relativistically generalize the discussion done in Ref. 1. For this purpose we start from the relativistic formulation of the kinetic equations given in Ref. 7. The implementation of the stochastic force is assumed to coincide with the one employed in the stochastic electrodynamics 5. Then, a formula for the equation satisfied by the fluctuation independent space-momenta joint distribution function is derived. It directly generalizes the expression obtained in 1 for the non-relativistic limit. A formula for the joint distribution is derived.

Then, the solution of the equations for the joint particle distribution is searched as expanded in a power series in the squared particle charge. Further, it is shown the existence of particular solutions for the joint distribution in the zeroth order in the charge. This implies that the density of particles at a given space time point is defined as a square of complex function satisfying the Klein-Gordon equation, in this considered interaction free approximation. This indicates that the Klein Gordon equation (or its non-relativistic limit the Schrodinger one) can be described as satisfying some of the equations of SQED in the first approximations. It can be cited that recently, in Ref. 15 a derivation was presented of the Klein Gordon equation, from a modified classical Hamilton-Jacobi equation for a particle interacting with random background forces. In addition an alternative derivation of the KG equations was also obtained in Ref. 13. The present discussion independently generalizes the kinetic discussion of the non relativistic analysis done in 1. The presentation here also clarifies the role of the positive and negative energy solutions, by showing that both of them separately imply the conservation of the total number of particles in the stochastic motions, assumed that external electromagnetic fields are absent.

By the side, it can be observed here that a special circumstance could help to overcome the known lack of strict positiveness of the ansatz for the joint distribution adopted in 1 and here. The KG equation, when seen as theory of particles lacks a standard definition for the position operator having eigen-functions like the Dirac’s Delta function. The position operator for this theory shows Gaussian like spatial behavior with non vanishing values within an spatial neighborhood of the size of the Compton wavelength of the particle (17). Then, it looks reasonable that a clear interpretation of the Q~0(x,p) as describing particles with a well defined position x, can run in troubles. Thus, one can imagine that a proper modification of the kinetic equations to take account of an extensive nature of the particles could lead to a consistent hidden variable interpretation of SQED in describing quantum mechanics. Assumed that above mentioned difficulty can be surmounted, an interesting extension of the work could seem feasible. For this, after including an external electromagnetic field, it seem possible to develop a picture in which both types of particles move randomly: one kind of them guided by the positive energy solutions and the other one (with opposite charges) moves as driven by the negative energy waves. The development of such a picture is an interesting envisioned extension of the work.

Further in the work, we also present some solutions for the zeroth order joint distribution function Q(x, p) which have positive values in all the phase space points. They have this property independently of the assumption done in the Olavo’s works 8 about the infinitesimal character of the of the Fourier conjugate variable of the particle momenta in the stochastic motions. It also can be stressed that the zeroth order equations for joint distribution function corresponds to the limit in which no stochastic action of the particles are effected. Therefore, these solutions could be more reasonable to be adopted in the zeroth order, since they describe uniform motions of free particles. Two kinds of localized solutions are found: One constructed as point-like localized spatial dependence which moves with the four-velocity associated to the 4momentum of the particle. The other kind is attained by employing Yukawa like localized solutions of the Klein-Gordon equation in their construction. Both types of joint coordinate momentum distributions in turn strongly suggest the possibility of describing extended particles, showing the surprising experimental properties exhibited by droplets moving over oscillating liquid surfaces 19,20. This idea comes from the suspicion about that after the action of stochastic forces (in higher orders in the coupling) both sorts of solutions might transform in extended wave packets surrounding a stochastic mean position of the particle moving with constant velocity υ=p2+m2. Such outcome is suggested in reference 21, in which the existence of such configurations is argued from a given proof of a stochastic Noether theorem. The investigation of the scattering properties of such solutions on two slits screens and potential walls, by example is expected to be considered elsewhere.

In Sec. 2 we introduce the basic notions of the relativistic kinetic theory. Next, in Sec. 3, the relativistically invariant equations for the mean value of the distribution and its random fluctuations are written. Further, in Sec. 4, the momentum Fourier transform of the mean joint distribution is introduced and the equations for it, are written. Section 5 considers the equations following in the first order zeroth approximation in the coupling with the stochastic forces. It is exposed how solutions of the Klein-Gordon equation define particular solutions of the relativistic kinetic equations, determining a possible joint distribution function in the assumed zeroth order in the charge. Next, Sec. 6 discusses how these special solutions determine particle distributions which conserve the total number of particles when the KG waves are assumed to be alternatively positive or negative energy modes. Further, in Sec. 7, we present other special solutions for the equation for the joint distribution function, satisfying the positiveness condition and being consistent with the lack of action of the stochastic forces in zeroth order of their coupling. Also, their possible links with the Couder’s experiments are identified. Finally, we advance some remarks linked with the argue presented in Olavo’s analysis about the positiveness of the joint distribution function 9,12.

The results are reviewed and commented in the Summary section.

2. The equation for the joint distribution function

Let us start by writing the relativistic invariant equation for the density of points in phase space R(x, p) for an ensemble of massive particles all evolving under the action of a stochastic 4-force Fμ (x, p) which was derived in Ref. 7

pμμRx,p+mFμx,pμRx,p=0. (1)

The 4-coordinates xμ will be considered in the metric

gμν=10000-1000000-100-1, (2)

where the natural system of coordinates will be also employed, in which the light velocity c = 1 and the time is the x0 coordinate and the Planck constant ħ = 1. The four momentum as usual, is given in terms of the 3-velocity as

pμ=m(1,υ)1-υ2 . (3)

In order to simplify the discussion, we will firstly consider that the external force vanishes. The particle density n(x,t) and the particle flow j(x,t) in this relativistically invariant case have the form

Nμx=nx,t,jx, t=dpp0pμ  Rx,p, (4)

in which the integration is over all the 3-momenta. In general the conventions defined in Ref. 7 will be employed. It will be assumed that the momenta values are defined on the mass-shell

p2-m2=0 .

In the present work, the force will assumed to be stochastically defined as in Ref. 5. That force had been argued to be invariant under Lorentz transformations 2-5. Therefore, the Eq. (1) also becomes relativistically invariant in form. As remarked before in Ref. 1, it was argued that the stochastic motion guided by that force, when taken in the non-relativistic limit, leads to the Schrodinger equation in the first steps of an iterative process of solution of the equations for the non relativistic distribution function. Therefore, as it was argued in Ref. 18, it can be suspected that the relativistic invariant motions determined by (1) could be related with the satisfaction of the Klein-Gordon equation, for which the non-relativistic limit is the Schrodinger one. This work is devoted to investigate this possibility. In order to make the discussion clearer let us argue in the next section that in the non-relativistic limit, the stochastic equations reproduce the ones employed in Ref. 1.

2.1. The non-relativistic limit of the equation

In this case since υ2/c21 the momentum and the external force can be approximately given by

pμ=m1,υ. (5)

Then, after considering υ=p/m the Eq. (1) reduces to

x0Rx,p+pmxR(x,p)+Fx,ppRx,p=x0R(x,p)+xpmRx,p+Fx,ppRx,p=0 . (6)

Further, assuming that Fx,p=Fx, in other words that the force is independent of the momentum, leads to

x0Rx,p+xpmRx,p+pFxRx,p=0. (7)

which is the same starting formula employed in 1. However, in order to arrive to this expression it was assumed that the force is not momentum dependent. But, the stochastic electric force term employed in 1 is momentum independent, an thus it makes the non-relativistic equation employed in Ref. 1 and the one employed here, equivalent in the non relativistic limit.

2.2. The adopted SQED relativistic random vacuum forces

Let us give a precise definition of the relativistic stochastic process under consideration. Note first that we had omitted a time t argument in the distribution in order to avoid the use of the non-relativistic invariant definition of the time. Therefore, the stochastic character of the process will be implemented by defining a large ensemble of particle trajectories in the phase space (x, p). Each of these trajectories will be defined by a solution of the Eq. (1) for a force given by a random realization of the relativistic invariant Lorentz force employed in SQED 5

Fμx,p=qmFvμxpv=qfμx,p, (8)

fμx,p=Fvμxpvm, (9)

in which the stochastic space-time dependent field intensity F γβ (x) is given by

Fγβx=γAβx-βAγx, (11)

Aβx=A0x,A(x), (11)

Ax=λ=12d k1wkϵ k,λhk,λ×sinkx-wkx0+θk,λ, (12)

where wk=k, ϵ(k,λ) are two unit polarization vectors associated to the wave vector k and satisfying

ϵk,λϵk,λ'=δλλ',   kϵk,λ=0, (13)

and the number h is defined as satisfying

π2h2=12wk. (14)

Finally, the phases θ(k, λ) are defined as independent random (one for each value of (k, λ)) and uniformly distributed in the interval (0, 2π) 5. In what follows, in place of the force expression, we will prefer to work with the defined above force per unit of charge

fμx,p=1qFμx,p=Fνμxpνm. (15)

2.3. Definitions for the operators and their kernels in joint coordinate-momentum space

We will consider in what follows linear kernels in the space of functions PS = {g(x, p)}of the particle coordinates x and momenta p, which explicitly written will make the expression to appear as cumbersome. Then, for any of such kernels K, which action on functions of the space g is defined as

g'x,p=dx' dp' Kx,p;x',p'gx',p', (16)

its compact operator expression will be defined according to the following equivalence rules

gx,pg, (17)

dx' dp' Kx,p;x',p' gx', p'K^ g. (18)

The special Delta function kernel δ8(x − x′, p − p′) = δ4(x−x′) δ(4)(p−p') will be simply defined as the identity Î, which will mean for the kernel associated to the inverse of K, the relation

K^-1K^=K^K^-1=I^

The local operators, like L = pμ (∂/∂xμ) and f = f μ (x,p)(∂/∂pμ) are also considered as kernels in the usual way

pμxμpμxμδ(8)x-x',p-p'L^, (19)

fμx,ppμfμx,p×pμδ(8)x-x',p-p'f^. (20)

3. The equations for the joint distribution function

Now, we will apply the method of smoothing (See 16 and 1) in order to reduce the Eq. (1) to a non random one for the coordinate-momenta joint distribution function over the defined ensemble of trajectories. The ensemble is generated by samples of the stochastic force, generated by the random phases of the electromagnetic modes θ(k,λ) taken for all the values of momenta and polarization (k,λ). For the further analysis, the distribution function will be decomposed in its average coordinate-momenta joint distribution Q(x,p) and its random fluctuations δQ(x,p) as

Rx,p=Qx,p+δQx,p, (21)

with

Qx,p=P^Rx,p, (22)

δQx,p=1-P^Rx,p, (23)

where the P^ is a projection operator satisfying P^2=1. After substituting these expression in equation (1) and applying alternatively P^ or (1-P^), the following two equations follow

L  Qx,p+qP^fμx,ppμδQx,p=0,L  δQx,p+q1-P^fμx,ppμQx,p=0, (24)

L=pμxμ. (25)

Consider now the retarded Green function in the space PS of the differential operator L = pμ (∂/∂xμ), satisfying

pμxμGx,p;x',p'=δ(4)x,x'δ(4)p-p'. (26)

In terms of its Fourier transform in the two arguments, the Green function G takes the form

Gx,p;x',p'=dq2π4dz2π41-ipμqμ×exp(-ix-x'μqμ-ip-p'μzμ)=dq2π41-ipμqμ×exp(-ix-x'μqμ)δ(p-p')=Gx-x'p'δp-p'L^-1. (27)

It is important to note here, that this expression for G indicates that the derivatives ∂/∂pμ do not commute with the operator G, since

pμGx,p;x',p'=G(x-x'|p')pμδ(p-p')=-G(x-x'|p')pμδ(p-p')=G(x-x'|p')δ(p-p')p'μ+p'μGx-x'p'δ(p-p')Gx-x'p'δp-p'p'μ. (28)

Therefore, this non commutativity of the momentum derivative with the propagator G made difficulty to show in this relativistic case an important property derived in Ref. 1: the quadratic dependence in the Fourier transform variable z of the momentum p, of some relevant quantities in the discussion. This lack of commutativity led us, further ahead in this work, to consider the expansion in the coupling in the equations, in place of the expansion in powers of z. Now, acting with the product of G and pμ (∂/∂xμ) on an arbitrary function g(x, p) it follows

dx'dp'G(x,p;x'p')p'μx'μg(x'p')=dx'dp' dq2πμδp-p'-ipμqμexp-ix-x'μqμ×p'μx'μgx',p'=gx,p, (29)

which himplies

dx'dp'pμxμGx,p,x',p'δx-x'δp-p'. (30)

As in reference we will now define a compact notation in order to eliminate the cumbersome appearance determined by the kernel structure of the Green function. This notation is described by expressing the above relation in the form

dx'dp'p'μx'μGx,p;x'p'p'μx'μL^L^-1=I^δx-x'δp-p'.

The use of these relations written above, after acting with the kernel G at the left of the relations (24) gives for the average of the distribution and its random part, the expressions

Q=-qL^-1P^f^δQ, (31)

δQ=-qL^-1I^+qL^-1I-P^f^-1P^f^Q. (32)

Then, the substitution of these relations in the first of the Eqs. (24) leads to the following equations for the joint distribution function Q (x, p)

L^Q=q2P^f^L^-1I^+qL^-1I^-P^f^-1f^Q,=q2P^f^L^-1n=0-1nqL^-1(I^-P^)f^nf^Q,=q2P^f^L^-1m=0q2mL^-1I^-P^f^2mf^Q. (33)

Therefore, the equation for the joint distribution function can be written in a compact form, which after expanded in powers of q2, takes the form

L^Q=q2J^q2Q, (34)

J^q2=m=0q2mP^f^L^-1L^-1I^-P^f^2mf^=m=0q2mJ^mq2. (35)

In this relation it has been employed that the mean value of an odd number of the random force functions 1,5 vanishes.

We recall that in that relation f^ is the operator corresponding to the kernel

f^fμx,ppμ. (36)

4. The momentum Fourier transformed joint distribution function

Let us perform now the Fourier transformation of the joint distribution over the momentum variable as follows

Qx,p=dzQ~x,zexp-ipμzμ, (37)

Q~x,z=dp2π4Qx,pexpi pμzμ. (38)

It can be mentioned that this concept had been defined and employed, by example in Refs. 1 and 8. It is also named as the Characteristic Function. Then, after Fourier transforming the Eq. (34) the following equation for Q(x, p) can be written

pμxμQx,p=q2m=0q2m×dx'dp'Jmq2x,p;x',p'Qx',p', (39)

ixμzμQ~x,z=q2m=0q2mdx'd z'×Jmq2x,iz;x',iz'Q~x',z', (40)

with the operator

Jmq2x,iz;x'iz'

operating in the space of function of the variables (x, z) is defined by

Jmq2x,iz;x'iz'=P^f^L^-1L^-1I^-P^f^2mf^×x,p;x'p'p8iz,p'iz'. (41)

Now, the mean value of a function of the coordinates and momenta at a specific space-time position x can be written in two forms as

Ax,px=1ρt(x)dp Ax,pQx,p=1ρt(x)Ax,izQ~(x,z)z=0, (42)

where the distribution function in the 3D-space points x and a given time x0 is given by

ρtx=dp Qx,p=m=0q2mρm(x)=m=0q2mdpQm(x,p)=m=0q2mQ~mx,0, (43)

in which the general expressions for the distribution function have been expanded in series of the squared charge as follow

Qx,p=m=0q2mQmx,p, (44)

Q~mx,z=m=0q2mQ~mx,z. (45)

Therefore, the general equation (40) can be written in the form

i zμxμQ~0x,z=0, (46)

i zμxμQ~mx,z=q2n=0q2dx'd z'Jnq2×x,iz;x'iz'Q~m-n-1x',z',m1. (47)

4.1. The interaction free limit q2 → 0

Let us consider now the satisfaction of the first of the iterative equations in which the random movement had been decomposed. In this case all the Qm (x,p) for m ≥ 1 will vanish and thus Q(x, p) = Q0 (x, p). It should be remarked here, that zeroth order equation for Q0 (x, p) exactly coincides with one obtained in Ref. 13 for the relativistic characteristic function in the absence of interaction with the random force. Then, the total density reduces to

ρtx=dp Q0(x,p)=ρ0x=dp Q0x,p=Q~0(x,0). (48)

The mean values in this limit have the expression

A(x,p)x=1ρ0(x)dp Ax,pQ0x,p=1ρ0(x)Ax,izQ~0(x,z)z=0, (49)

in which the distribution function has been expanded in series of the squared charge.

Let us define for what follow

ρx=ρ0x. (50)

Then, the equation for Q~0x,z.

i zμ xμQ~0x,z=0, (51)

can be derived after expanding the exponential in powers of z. A similar equation had been obtained also in Refs. 13 where relativistic expressions for the characteristic function, not coming from a stochastic discussion was independently considered.

It is helpful to write the Fourier transforms in the zero order in q2

Q0x,p=dz Q~0x,zexpi pμ zμ, (52)

Q~0x,z=dp2π4Q0x,pexp-i pμzμ. (53)

The mean value formula reduces to

A(x,p)x=1ρ(x)dp Ax,p)Q0x,p=1ρxA x,izQ~0(x,z)z=0, (54)

which allows to write the lowest order equation in the form

izμxμQ~0x,z=-xμdp2π4Q0x,ppμexp(-ipz)=xμ(ρxpμexp-i p zx=0. (55)

We will now introduce new variables z+ and z in substitution of the variables x and z. The change is defined as

zμ+=xμ+β zμ,  zμ-=xμ-βzμ, xμ=12zμ++βzμ+,  zμ=12βzμ+-βzμ+,  (56)

and for the derivatives

zμ++μ=12xμ+β zμ,  zμ--μ=12xμ-β zμ, (57)

zμ=zμ++zμ-,  zμ=βzμ+-zμ-. (58)

The above kind of transformations had been suggested by reference 1 and were also employed in reference 13, 14 in deriving the Klein-Gordon equation. These relations allow to derive the identities

2zμzν=β22xμxν-2β2-μ+v++μ-v (59)

2zμxμ=β2zμ+zμ+-2zμ-zμ-. (60)

This last equation permits to write the zeroth order equation for Q~0 in the form

izμxμQ~0x,z=βi×2zμ+zμ+-2zμ-zμ-Q~0x,z=0. (61)

4.2. Satisfying the higher order equations

In order to directly satisfy the set of equations for n ≥ 3 let us assume Q~0(x,z) in the form

Q~0x,zQ~0z+,z-=Ψ*z+Ψz-. (62)

It should be noted that this form is suggested by the analysis done in Ref. 1. This assumption directly led to the character of solutions of the Schrodinger equation to the entering function Ψ. Then we will follow this assumption in searching for a covariant generalization of the discussion in 1. However, it can be noted that the satisfaction of the quantum equations as implied by the stochastic theory being constructed could be more naturally expected to appear after including the stochastic effects by considering higher than zeroth orders in the coupling with the stochastic forces. Therefore, in the last part of this work we also search for zeroth order joint distributions being positive definite and also consistent with a free motion of a localized particles when the random forces are absent in the zeroth order in the coupling expansion.

Then, substituting the above commented assumed form in the zeroth equation leads to

1Ψ*(z+)2zμ+zμ+Ψ*z+-1Ψz-2zμ-zμ-Ψz-=0. (63)

But this relation is directly satisfied if Ψ obey the linear equation for any argument u and fixed value of the parameter M

2uμuμΨu-M2Ψu=0. (64)

Then, the satisfaction of this equation implies

ixμzμzμ1zμ2zμnQ~0x,z=0,  m=1,2,, (65)

a condition which will be helpful in the further discussion.

Thus, we had been able to find a solution of the equation describing the free approximation of the random process.

The equation for Q~0 can be also written as

ixμzμQ~0x,z=-xμ×dp2π4Q0x,ppμexp-ipμzμ=-xμρxpμexp-ipμzμx=0, (66)

which after expanding the exponential in powers of p.z gives the following set of equations

xμρxpμx=0. (67)

xμρ(x)pμpvx=0, (68)

xμρ(x)pμpμ1pμ2pμnx=0, n=1,2,. (69)

The last of these relations is directly implied by Eq. (65). Thus, let us study in what follows the satisfaction of the first two equations after the adopted ansatz

Q~0x,z=Ψ*zμ+Ψzμ-.

We will use now

pμpvx=-1Q~0x,0 zμzvQ~0x,z, (70)

and the general relations

zμ1Q~0(x,z)zμQ~0x,z=-1Q~0x,z2×zvQ~0x,zzμQ~0(x,z)+1Q~0x,zzvzμQ~0x,z, (71)

which after evaluated in z = 0 permits to write

pμpvx=pμxpvx-2zμzvlnQ~0x,z)z=0=pμxpvx-β22xμxvlnQ~0x,zz=0+σμv, (72)

σμv=2β2-μ+v++μ-vlnQ~0z=0, (73)

where it was used relation (59).

However, the assumed form of the zeroth order distribution allows also to find

σμv=2β2-μ+v++μ-vlnQ~0z=0=2β2-μ+v++μ-vlnΨ*z++lnΨz-z=0=0. (74)

Thus, the resting two equations which remaining to be verified in their compatibility with the ansatz can be written. as

xμρxpμx=0,xμρx(pμxpvx-β22xμxvlnQ~0(x,z))z=0=0. (75)

The derivation of these equations was directly suggested by their non relativistic counterparts in 1. They also coincide with the ones derived in Ref. 13. In this work an alternative relativist discussion of a similar ansatz for the characteristic function, not coming from an stochastic formulation, was earlier and independently presented.

5. The Klein-Gordon equation in the noninteracting limit

The q2 = 0, will be called the ”non interacting or free approximation”. In this section we will study the compatibility of the two Eqs. (75) with the expression assumed for the zeroth order joint distribution function. After checking this, it will follow that in this free limit, the equations admit solutions for joint distribution function which are defined by waves solving the Klein-Gordon equations. This result directly generalizes the derivation of the Schrodinger equation in Ref. 1. Analogous derivation of the KG equation was formerly and independently given in Ref. 13. Firstly, let us search for solutions of the set of the two equations

xμρ(x)pμx=0, (76)

xμρxpμxpvx-β2ρ(x)2xμxvlnρ(x)=0, (77)

where it has been substituted ρx=Q~(x,0). Further, let search for solutions in which the mean momentum value at a given space-time point x is defined by

pμx=xμSx. (78)

After substituting in (76) and (77), it follows

xμρ(x)xμS(x)=0, (79)

xμρ(x)xμS(x)xvS(x)-β2ρx2xμxvlnρ(x)=0. (80)

The second equation can be transformed as follows

xμSxxμxvSx-β2ρ(x)xμρx2xμxvlnρx=0, (81)

xμ12xμS(x)xμS(x)-β2ρxxμρx2xμxvlnρx=0. (82)

Further, the density dependent term can be expressed as the divergence of a vector (as the first term also is) as follows

1ρxxμρx2xμxvlnρx=xv12xμlnρ(x)×xμlnρx+ρ(x)2xμxμlnρ(x)=xv-12×xμlnρxxμlnρx+1ρx2xμxμρ(x). (83)

Therefore, the following relation arises

xv12xμSxxμSx+β22xμlnρ(x)xμlnρx-β21ρ(x)2xμxμρ(x)=0. (84)

This equation implies,

12xμSxxμSx+β22xμlnρ(x)xμlnρx-β21ρx2xμxμρx=ctc. (85)

Now, if we fix the constant to be a positive value, given by ctc = m2, the searched joint distribution function satisfies the two equations

xμρ(x)xμS(x)=0, (86)

12xμSxxμSx+β22xμlnρ(x)xμlnρx-β21ρ(x)2xμxμρx=m2. (87)

5.1. The satisfaction of the Klein-Gordon equation

Consider now expressing the KG equation for a complex scalar field ϕ(x) in terms of the phase function S(x) and a positive density ρ(x) defined as follows

ϕx=ρ(x)12exp(i S(x)). (88)

Therefore, substituting in the KG equation it follow

2+m2ϕx=0,

2xμxμ+m2ρ(x)12expi S(x)=0.

which after separating the real and imaginary parts and equalizing both of them to zero, leads to the two equations

xμSxxμlnρx+1ρ(x)xμxμSx=0, (89)

12xμSxxμSx+14xμlnρ(x)xμlnρ(x)-121ρx2xμxμρxm2. (90)

It can be noted that Eqs. (89) and (90) become fully equivalent to (86) and (87) after assuming that the constant β 2 takes the value

β2=12 . (91)

Observe that this value coincides with the one before derived for the non-relativistic situation in 1. Therefore, under the defined non-interacting approximation, it followed that the space time distribution

ρx=dp Qx,p=Q~(x,0)

and the also space time function S (which determines the mean value of the omentum at a given space time point through

pμx=xμS(x)

both define a complex function

ϕx=ρ(x)12expi S(x)

satisfying the KG equation.

In order that the searched solution of the distribution

Q~0x,0=Ψ*z+Ψz-, (92)

can be compatible with the solution for the KG waves generating the density ρ(x) and the phase function S, the relation

Q~0x,0=Ψ*z+Ψz-|z0=Ψ*xΨx=ρx, (93)

implies that the constant M defining the ansatz, and m defining the solution of the Hamilton-Jacobi equation, should coincide.

Up to now we have discussed the generalization of the de la Peña-Cetto derivation of the set of equations which should be satisfied by the joint distribution function in the non-relativistic limit of the SQED. We also derived a particular solution for the joint distribution function which implies the satisfaction of the Klein-Gordon equation.

6. About the role of positive energy solutions of the KG equations

Let us discuss an important physical question related with the obtained solutions. Since the particles which are assumed to undergo the random forces, are relativistic classical particles with rest mass m, it is natural to assume that the KG equation solution describing the stochastic motion should be expected to be a positive energy one. This circumstance is in certain form confirmed by a special property of the positive energy solutions (and also the negative energy ones): their total probability conserves in time, that is

x0dxρx=x0dx ϕ*xϕx=0. (94)

This property can be easily derived after considering the following equations also satisfied by the positive or negative energy solutions

ix0ϕx=±m2-2ϕx=±m=0-1m2m2mϕx.

For proving the condition, consider that the spatial integral of the density for a positive energy (or a negative energy one) solution can be written and transformed in the way

x0dx ρx=x0dxϕ*(x)ϕ(x)=dxx0ϕ*xϕx-ϕ*(x)x0ϕ(x)=1idxm=0-1m2m2mϕ*xϕ(x)-ϕ*xm=0-1m2m2mϕx=0, (95)

in which there had been performed integration by parts over the all the derivatives forming the appearing series. Thus, either the positive or negative energy solutions define a stochastic motions conserving the total number of particles.

This result suggests the possibility of simultaneously consider two stochastic motions: both of them associated with positive energy solutions: but having opposite values of the charges. Let us indicate the positive energy solutions by ϕ + (x) and the negative energy ones by ϕ − (x). But we can also define positive energy functions starting from the negative waves by defining φ+x=φ+x0,x=ϕ-(x0,x). In this way all the solutions of the KG equations could participate in defining a combined stochastic process in which two kinds of particles participate: one kind with positive charges and the another with negative ones. We suspect that this construction can be extended to a full covariant stochastic theory of complex scalar particles which could appear to be a hidden variable theory for the quantum field theory of such particles. This question will be explored in extensions of this study. However, before considering this problem, the approach should made consistent, by finding a positive definite joint distribution function. The possibilities for this will be discussed in the next section.

In conclusion, for defining the found solution of the zeroth order equation for the joint distribution as a well defined hidden variable theory, it rest only (in the relativistic as a well as in the non-relativistic cases) to check wether or not, the derived joint distribution can obey the important positivity property, which is required by its character as a coordinate and momenta distribution of classical particles of mass m, all showing a relativistic momenta obeying the mass shell condition p2m2 = 0. In the coming sections we will remark on these questions.

7. Positiveness of the joint distribution and Olavo’s analysis

Finally, in this section we want to discuss the question about the required positive character of the joint coordinate-momentum distribution, needed if at the end, the SQED approach can furnish a consistent hidden variable approach to quantum mechanics. It is known that the Fourier transformed joint distribution of the form Q~0 x,z=Ψ*(z+)Ψ(z-), does not give a solution to this difficulty, since in general, the joint coordinate-momentum distributions following from the inverse Fourier transform of Q~0(x,z) (for all values of z) are not positive definite in general. Thus, the before discussed here solution for the joint distribution has a formal value, but for every solution of the KG equation does not furnish a physical positive joint distribution.

A step in the solution of this relevant interpretation problem, was given in Refs. 8 to 13. In these works it has been emphasized the idea about that by retaining only the second order in the expansion in the conjugate variable of the momenta in the Fourier transform of the joint distribution, a proper positive distribution is obtained 9,12. Thus, the infinitesimal values of these Fourier conjugate variable of the momentum was elevated to a central assumption. The consistency of this analysis had been argued by few alternative derivations 8-13.

It should recalled that we consider that the derivation of the Schrodinger (or the Klein-Gordon one) from the equations for Q(x, p) in the absence of the vacuum stochastic force, is in some sense inconsistent with the spirit of the SQED, in which the stochastic action of the force is expected to define the quantum properties. What could be more natural is to start the iterative construction of higher order in the stochastic coupling, by employing also a zeroth order solutions for the joint distribution, but reflecting a free classical motion of the particle. This approach directly suggests that that SQED could be valid as describing quantum mechanical effects, but in a form compatible with the recent experimental results of Ives Couder. Those surprising findings show that mechanical systems, like liquid droplets moving on oscillating liquid surfaces, can exhibit quantum mechanical properties, as tunnel effects and double slit interference 19,20.

Finally, in a last subsection we present some remarks on some special issues in the discussions given in 9,12.

7.1. Positive zeroth order joint distributions and Couder’s findings

In this section we will present a solution of the equation for the joint distribution Q(x, p) being positive definite in phase space. These solutions strongly suggests a possible link with the recent studies on quantum mechanical properties in the movements of droplets 19,20. The picture could be as follows. After considering the action of stochastic forces on the free solutions to be presented below, the stochastic movements of the particle could form a localized standing wave which center of mass could be in uniform-motions. The possible existence of such standing waves, showing counterparts moving with different velocities under a Lorentz transformation, had been suggested in a recent work 21. In it a form of a stochastic Noether theorem had been introduced. The extender nature of these Lorentz invariant stochastic solutions can be imagined to show properties being similar to the Couder’s droplets moving over liquid surfaces. If such is the case, the quantum particles of SQED could be imagined to describe the quantum properties in nature through the Couder’s mechanism. We expect to study this possibility elsewhere.

The construction of the positive joint distributions is based in the Yukawa potential like solution of the Klein Gordon equation in the presence of sources

ψx=14πexp-m rr, (96)

r=x12+x22+x32 , (97)

satisfying the KG equation

2-m2ψx=-δ(3)x (98)

where δ3 is the three dimensional Dirac delta function. After performing a Lorentz transformation to a frame moving with velocity υ along, let say, the x1 axis, followed by a shift in the position of the origin of coordinates (at vanishing time t) to an arbitrary point xk, this Yukawa like solution becomes a “moving” one, of the form

ψxk,x,p=14πexp-mp2+m2m2x1-p1p2+m2t-x1k2+x2-x2k22+x3-x3kp2+m2m2x1-p1p2+m2t-x1k2+x2-x2k22+x3-x3k (99)

x=t, x,    p= p0,p=m1-υ2,m υ1-υ2,    xk=0,xk. (100)

In this change it had been used the relations

υ=p1p2+m2,    1-υ2=mp2+m2 . (101)

This function also can be expressed in a rotational invariant form as

ψxk,x,p=14πexp-mp2+m2m2x-pp2+m2t-xkpp2+x-xkI-ppp2x-xkp2+m2m2x-pp2+m2t xkpp2+x-xkI-ppp2x-xk (102)

where

I-ppp2

is the projection tensor on the plane orthogonal to the velocity and p/p is a unit vector in the direction of the velocity.

The function ψxk(x,p) is positive definite in the phase space (x, p). Also, in the rest frame p=(m,0), it directly satisfies the equation for the joint momenta-coordinate distributions, since it is time independent and the three velocity vanishes. Then

pμxμψxk,x,p=0. (103)

Having this equation a covariant form, it should be also valid after performing any Lorentz transformation of the coordinates and momenta. Now, we can define a set of N points xk , k = 1, 2, ...., N. Then, by superposing the functions of the type (102) for all the values of k, more general solutions can be constructed. They will describe a set of N localized solutions of the Klein-Gordon equation with sources. Also, the new distributions also will be positive quantities in all the mass-shell phase space. These solutions of the equations pμ (∂/∂xμ)Q(x, p) = 0, showing a positive joint distribution have the forms

Qx,p=kψxk,x,p, (104)

where for each k value, the velocity for the distribution with index k, points in an arbitrary direction.

It can be noticed that the proposed solution is not rigorously point like. It has distributions showing a Compton wavelength extension. However, strictly point-like solutions for the joint distribution function are directly obtained in the form

Qx,p=kΦxk,x,p, (105)

pμxμΦxk,x,p=0, (106)

in which

Φxk,x,p=14πδ(1)p2+m2m2x1-p1p2+m2t-x1k2+x2-x2k22+x3-x3k) (107)

x=t,x,  p=p0,p=m1-υ2,m υ1-υ2,  xk=0,xk. (108)

These joint distributions are positive in the whole phase space and solves the zeroth order equation. They represent the uniform motions of localized particles in the absence of stochastic perturbations, being more compatible with considered zeroth order equation under consideration. Below, we remark on the possible connections of these distributions with the Couder’s findings of quantum mechanical properties droplets moving over oscillating liquid surfaces.

7.1.1. Possible links with the Couder’s systems

Finally, let us here very roughly argue about the possible links of these special joint distributions with the Couder’s experimental results 19,20. For this purpose, let us qualitatively discuss a situation in which one of the solution ψxk(x,p) with a given velocity, perpendicularly approaches a wall having two slits holes. Let us first note that the singularity of the solution is similar to the one in the Coulomb potential field and its classical energy is infinite. Thus, since we will assume that the particle has a finite mass m, some negative contribution to the bound energy (which is concentrated in the singularity point) should cancel the infinite and positive contribution of the classical energy outside the point. Then, let us consider a sphere centered in the singularity at any instant, such that the field energy outside the sphere coincides with the total mass of the particle. Therefore, the contribution to the total energy inside the defined sphere, should vanish (the negative cohesive term should cancel the infinite positive energy laying outside the point, but inside the sphere). However, having not net mass, the system inside the sphere, might be suspected to weakly contribute to the free dynamics of the particle. Then, we have that outside the small sphere the system will satisfy the Klein-Gordon equation and also will move as a whole with constant velocity, Thus, the solution ψxk(x,p) might perhaps also be approximately represented by a wave packet solving KG equation, but with a momenta distribution showing non vanishing values only in a small neighborhood of the momenta component p1=mυ/1-υ2, in order that the packet shows a constant velocity. If this idea is valid (in spite of its very rough nature) then

ψxk,x,pdpfpexp-ip+m2t+px idp1dp2dp3 fpexp-i pp+m2pt-x1+i p2x2+i p3x3dp1dp2dp3 fp×exp-i pυ t-x1+i  p2x2+i p3x3. (109)

However, it should be noted that the singular solution is localized within region of the size of a Compton associated to the free mass parameter m. This is a small quantity for usual particles as, by example electrons. Thus, in order that the representation (109) could be valid, the momentum bandwidth of the integral in (109) should be larger than one over the spatial width of the singular solution δx ≃ (1/m), that is, larger than m. In the relativistic limit p1 ≫ m, this condition can be satisfied. However, for mass parameters m larger than the electron’s one, and in the non relativistic limit, this rule can not be imposed. Thus, we expect that the parameter m should have a small value in order to allow implementing the representation (109) in the non relativistic limit. Since the Couder’s experiments are done for movements of the droplets over surfaces having massless propagating modes, even the vanishing mass parameter m could be allowed. As positive factor it can be noticed that after considering the effect of the stochastic forces, the wave packet bound to the particle will be constituted by real waves stochastically excited by the random moving of the particle. Therefore, opening possibilities for showing Couder’s like effects.

To end the argue, note that (109) is a wave corresponding to particles of mass m and wavelength λ = 2π/p. Therefore, assumed that the two slits to which the wave approaches, have a separation of a similar size to the wavelength λ, these modes should tend to be scattered by the action of the slits. Clearly, if the waves were completely free ones, this scattering should work. However, since the singularity is expected to maintain the structure of the particle when the scattering process occurs (as it happens in the Couder’s experiments) the whole effect of such ”dispersion forces” could be suspected to be reduced to control the movement of the singularity (if it passes through the slits) to be pointing in the directions of the usual interference maxima. This argue also indicates (as noted above) that the results of Couder could perhaps be described by the found positive solutions of the equation for the joint distribution function after incorporating the action of the vacuum stochastic forces. The study of the two slits scattering on such configurations is expected to be considered elsewhere. The clear limitations of the above discussion for pure free solutions could be solved by the noted consideration of the stochastic forces. They have the chance of converting, the uniformly moving point like distribution functions in sorts of solitonic stochastic wave packets showing similar properties as the Couder’s droplets.

7.2. Remarks on the Olavo analysis

In the discussion given in 9,12 it was obtained the general relation

-2ρx,t2lnZx,δx,tδx2δq=0=-2ρ(x,t)2lnRx;t24 x2 (110)

=dp  p-xSx,t2Fx,p;t. (111)

between the characteristic function Z and its inverse Fourier transform, the joint coordinate momentum distribution function F.

We intend to underline here that this exact equation (only assuming the product structure for Z(x, δx, t) = Ψ* (x + (δx/2), x − (δx/2)) ) implies a strong link between the sign of the second derivative (−∂2 ln R(x; t)2/ ∂x2) and the strict positivity of the joint distribution F(x, p; t). This connection determines that when (∂2 ln R(x; t)2/ ∂x2) in some spacial points turns to be negative, then F(x, p, t) can not be strictly positive. But, it can be noted that it is possible to give initial conditions at some time t o to defining arbitrary values for the spatial dependence of the wave functions. This property allows fixing the spatial dependence of the coordinates x for the density function in an arbitrary form. Therefore, the relation

-2lnRx;t2x2=1ρ(x,t)-2ρ(x,t)x2+1ρ(x,t)p(x,t)xρ(x,t)x (112)

can be always made negative, at a some time to and some point xo, by fixing a minimum of the density at to and xo. This conclusion means that solutions will always exist showing a non positive joint distribution function. Therefore, since this result is determined by the characteristic function up to second order in δx, that is in the limit considered in the Olavo’s discussion, it is concluded that this approach still present a difficulty in consistently predict an always positive definite joint distribution function.

The above conclusion seems contradictory with the fact the in references 9,12, it was argued the positivity of the joint distribution function. Let us remark on this point. In the notation of these works (in which the discussion was non relativistic), the characteristic function for infinitesimal values of the conjugate variable to the momenta δx was firstly obtained up to second order in the expansion over δx, in the form

ZQx,δx;t=Rx;t2+i δxRx;t2s(x;t)x (113)

+δx22Rx;t242lnRx;t2x2-Rx;t22S(x;t)x2, (114)

where ρ(x, t) = R(x; t)2 is the particle density of the derived Schrodinger equation expressed as the absolute value of the wave function R(x, t) and p-x,t=S(x;t)/x is the gradient of the phase of the wave function S (x, t). In a central step in the discussion in Ref. 9, after employing the assumed differential character of δx, the above expression was substituted by the following one.

ZQx,δx;t=Rx;t2exp+i δxs(x;t)x+δx282lnRx;t2x2. (115)

This is a valid transformation up to the assumed second order in δx. In other words, a series of terms “completing” the appearing exponential function were added. Note that all such terms have at least a cubic dependence in δx. If the further use of the form of the characteristic function is reduced to evaluate its derivatives up to second order in δx, or integrals over δx, for small intervals in which the second order approximation gives a good approximation, the use of (115) is allowed. However, the obtaining of the expression for the joint distribution function in 9 rests in evaluating the Fourier inverse transformation

Fx,p;t=Rx;t2-dδxexp-i dxp-sx;tx+δx282lnRx;t2x2=Rx;t2πexpp-s(x;t)x22lnR(x;t)22x2-2lnR(x;t)28 x2, (116)

which again is a positive definite expression only for the special class of wave functions having −(∂2 ln R(x; t)2/8 ∂x2) > 0, as argued before.

That is, even after assuming the modified expression for the characteristic function (115) as valid for finite values of δx, (as required to evaluate the inverse Fourier transform) the resulting expression can result as a non positive definite quantity for existing special wave functions.

Summary

We have presented a generalization of the non-relativistic stochastic quantum mechanics introduced by de la Peña and Cetto 1. The discussion starts form the description of the random motions of a particle under the action of a relativistically invariant stochastic force defined in reference 5. It is checked that in the non relativistic limit the starting equations reduce to the ones employed in 1. Then, the set of equations for the joint distribution is expanded as a series of the particle charge. The free approximation, that is, the equation following in the zero order of the expansion in the charge is solved by considering the ansatz for the distribution adopted in 1. After this, it is argued that the spacetime probability distribution of the stochastic process ρ(x) and a phase function S(x) which gradient determines the momentum mean value at a definite space-time point, define a complex scalar function satisfying the KG equation through ϕx=ρ(x) exp(iS(x)). It is argued that the total number of particles N=dx ρ(x) determined by the spacetime distribution conserves in time, if the KG solutions considered for determining the distributions are either positive or negative energy modes.

It is also addressed the fact that above mentioned ansatz solution, is associated to a joint distribution function possibly giving negative values. In this sense, it is remarked that there are special circumstances that could help to overcome this limitations. One of them is the following, the KG equation, when seen as theory of particles shows the interesting effect that it lacks a standard definition for the position operator eigenfunction as the Dirac’s Delta function. In place of it, the appropriate position operator for this theory has a Gaussian like appearance showing non vanishing values in spatial regions of the size of the Compton wavelength of the particles 17. Therefore, it looks reasonable that the proper interpretation of the Q~0(x,p) as describing particles with a well defined position x can show difficulties. This leads to the idea that a proper modification of the kinetic equations to take account of this extensive nature of the particles occasionally could still lead to a consistent equivalence of the SQED with QM and QFT.

The work is also presenting solutions for the joint distribution functions which obey the positive conditions required by a proper classical distribution. Such solutions in a sense, seems that can be more natural (in the zeroth order in the coupling with the stochastic forces) than the ansatz considered in 1. It is also argued that these solutions could be related with systems of particles showing the Couder’s experimental results, if the localized solutions employed for their construction can be approximately represented as Couder’s like wave packets after including the effects of the stochastic forces 19. The investigation of this particular possibility is expected to be considered elsewhere.

They are exposed some remarks in connection the Olavo’s discussion in reference 8-13 directed to define a positive definite joint distribution function. They identify some estimated limitations of the approach in connection with the central question of defining a positive result for this quantity.

Finally, assumed that the identified difficulties with positive condition of the joint distribution (the one adopted in 1) can be properly solved, the discussion opens possibilities for the extension of the SQED analysis. Of particular interest is: The possibility of generalize the discussion to describe the stochastic evolution of two sets of particles: one described by the positive energy solutions and the other by the negative energy ones. After to also including the presence of an external electromagnetic fields, this construction seem to offer opportunities for describing the creation and annihilation particles by the electromagnetic field. The interaction terms, could result to be sources of the variation in the total numbers of positive or negative charged particles generated by the annihilation or creation of particles due to the action of the electromagnetic field. The search for the connection of this construction with the quantum field theory of the complex scalar field is a further question of interest to explore.

Acknowledgments

The support granted by the N-35 OEA Network of the ICTP is greatly appreciated.

References

1. L. de la Peña and A.M. Cetto, J. Math. Phys. 18 (1977) 1612. [ Links ]

2. T.W. Marshall, Proc. R. Soc. A 276 (1963) 475. [ Links ]

3. T.W. Marshall, Proc. Camb. Philos. Soc. 61(1965)537. [ Links ]

4. T.W. Marshall, Nuovo Cimento 38 (1965) 206. [ Links ]

5. T.H. Boyer, Phys. Rev. D 11 (1975) 790. [ Links ]

6. R. T. Cahill and S. M. Gunner, Fizika B 7 (1998) 171. [ Links ]

7. S.R. de Groot, W.A. van Leewen and Ch. G. van Weert, Relativistic kinetic theory, North Holland, Amsterdam (1980). [ Links ]

8. L.S.F. Olavo, Physica A 262 (1999) 197. [ Links ]

9. L.S.F. Olavo, Physica A 271 (1999) 260. [ Links ]

10. L.S.F. Olavo, Phys. Rev. A, 61 (2000) 052109. [ Links ]

11. L.S.F. Olavo, Foundations of Physics 34 (2004) 891. [ Links ]

12. L.S.F. Olavo, Annals of Physics 327 (2012) 1391. [ Links ]

13. L. S. F. Olavo, Foundations of Quantum Mechanics: Special and General Relativistic Extensions, in Contemporary Research in Quantum Systems, Zoheir Ezziane, Ed. Nova Science Publishers, Inc. (2014). [ Links ]

14. L.S.F. Olavo, Foundations of Quantum Mechanics: Special and General Relativistic Extensions in Contemporary Research in Quantum Systems, Zoheir Ezziane, Ed. Nova Science Publishers, Inc. (2014). [ Links ]

15. G. Grössing, Derivation of the Schrodinger equation and the Klein-Gordon equation from first principles, quant-ph/0311109 (2003). [ Links ]

16. U. Frisch, in Probabilistic methods in applied physics, Vol 1, (Academic Press, New York 1968). [ Links ]

17. S. S. Schweber, An introduction to relativistic quantum field theory, (Dover Publications, Mineola, New York 2005). [ Links ]

18. A. Cabo-Bizet and A. Cabo, Phys. Lett. A 359 (2006) 265-271. [ Links ]

19. Y. Couder, S. Protiére, E. Fort and A. Boudaoud, Nature 437 (2005) 208. [ Links ]

20. A. Eddi, E. Fort, F. Moisy and Y. Couder, Phys. Rev. Lett. 102 (2009) 240401. [ Links ]

21. A. González Lezcano and A. Cabo Montes de Oca, A stochastic version of the Noether Theorem, arXiv:1709.09295v1 [condmat.stat-mech] 27 Sep 2017. [ Links ]

Received: September 14, 2017; Accepted: November 27, 2017

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License