# Noise perturbations in the Brownian motion and quantum dynamics

## Abstract

The third Newton law for mean velocity fields is utilised to generate anomalous (enhanced) or non-dispersive diffusion-type processes which, in particular, can be interpreted as a probabilistic counterpart of the Schrödinger picture quantum dynamics.

If we consider a fluid in thermal equilibrium as the noise carrier, a kinetic theory viewpoint amounts to visualizing the constituent molecules that collide not only with each other but also with the tagged (colloidal) particle, so enforcing and maintaining its observed erratic motion. The Smoluchowski approximation takes us away from those kinetic theory intuitions by projecting the phase-space theory of random motions into its configuration space image which is a spatial Markovian diffusion process, whose formal infinitesimal encoding reads:

$$d\stackrel{}{X}(t)=\frac{\stackrel{}{F}}{m\beta }dt+\sqrt{2D}d\stackrel{}{W}(t).$$
$`(1)`$   
In the above $`m`$ stands for the mass of a diffusing particle, $`\beta `$ is a friction parameter, D is a diffusion constant and $`\stackrel{}{W}(t)`$ is a normalised Wiener process. The Smoluchowski forward drift can be traced back to a presumed selective action of the external force $`\stackrel{}{F}=\stackrel{}{}V`$ on the Brownian particle that has a negligible effect on the thermal bath but in view of frictional resistance imparts to a particle the mean velocity $`\stackrel{}{F}/m\beta `$ on the $`\beta ^1`$ time scale, . The noise carrier (fluid in the present considerations) statistically remains in the state of rest, with no intrinsic mean flows.

An implicit phase-space scenario of the Brownian motion refers to minute acceleration/deceleration events which modify (infinitely gently, , at a generic rate of $`10^{21}`$ times per second) velocities of realistic particles. Clearly, the microscopic energy-momentum conservation laws need to be respected in each separate collision event. In contrast to derivations based on the Boltzmann equation, this feature is completely alien to the Brownian motion theory. That energy-momentum deficit is one of ”forgotten” or tacitly ”overlooked” problems in the standard theory of the Brownian motion, see e.g. .

An issue of origins of general isothermal flows in the noisy environment is far from being settled, . Namely, the standard theory does not account (or refers to phenomena for which that is not necessary) for generic perturbations of the random medium as a reaction to the enforced (exclusively due to the impacts coming from the medium) particle motion. An ability of the medium to perform work i. e. to give a kinetic energy to Brownian particles (strictly speaking, in the local mean, hence on the ensemble average) appears to be a universal feature of the thermal bath even in the absence of any external force.

Then, there arise problems with the thermal equilibrium, . The tagged particle always propagates ”at the expense” of the bath, which nonetheless remains ”close” to its thermal equilibrium. In view of the involved local cooling and heating phanomena, which are necessary to keep in conformity with the second law of thermodynamics, the bath should develop local flows and thus in turn actively react back to what is being happenning to the particle in the course of its random propagation. Such environmental reaction, while interpreted as a feedback effect, surely can be neglected in an individual particle propagation sample. However, an obvious Brownian particle energy-momentum ”deficit” on the ensemble average (kinetic energy and momentum is carried locally by diffusion currents), indicates that the medium reaction may give observable contributions on the ensemble average, as a statistically accumulated feature. For a statistical ensemble of weakly-out-of-equilibrium systems, an effective isothermal scenario can be used again, but then mean flows in the noise carrier are ”born” as a consequence of the averaging, .

It is well known that a spatial diffusion (Smoluchowski) approximation of the phase-space process, allows to reduce the number of independent local conservation laws (cf. -) to two only. Therefore the Fokker-Planck equation can always be supplemented by another (independent) partial differential equation to form a closed system.

If we assign a probability density $`\rho _0(\stackrel{}{x})`$ with which the initial data $`\stackrel{}{x}_0=\stackrel{}{X}(0)`$ for Eq. (1) are distributed, then the emergent Fick law would reveal a statistical tendency of particles to flow away from higher probability residence areas. This feature is encoded in the corresponding Fokker-Planck equation:

$$_t\rho =\stackrel{}{}(\stackrel{}{v}\rho )=\stackrel{}{}[(\frac{\stackrel{}{F}}{m\beta }D\frac{\stackrel{}{}\rho }{\rho })\rho ]$$
$`(2)`$   
where a diffusion current velocity is $`\stackrel{}{v}(\stackrel{}{x},t)=\stackrel{}{b}(\stackrel{}{x},t)D\frac{\stackrel{}{}\rho (\stackrel{}{x},t)}{\rho (\stackrel{}{x},t)}`$ while the forward drift reads $`\stackrel{}{b}(\stackrel{}{x},t)=\frac{\stackrel{}{F}}{m\beta }`$. Clearly, the local diffusion current (a local flow that might be experimentally observed for a cloud of suspended particles in a liquid) $`\stackrel{}{j}=\stackrel{}{v}\rho `$ gives rise to a non-negligible matter transport on the ensemble average, .

It is interesting to notice that the local velocity field $`\stackrel{}{v}(\stackrel{}{x},t)`$ obeys the natural (local) momentum conservation law which directly originates from the rules of the Itô calculus for Markovian diffusion processes, , and from the first moment equation in the diffusion approximation (!) of the Kramers theory, :

$$_t\stackrel{}{v}+(\stackrel{}{v}\stackrel{}{})\stackrel{}{v}=\stackrel{}{}(\mathrm{\Omega }Q).$$
$`(3)`$   
An effective potential function $`\mathrm{\Omega }(\stackrel{}{x})`$ can be expressed in terms of the forward drift $`\stackrel{}{b}(\stackrel{}{x})=\frac{\stackrel{}{F}(\stackrel{}{x})}{m\beta }`$ as follows: $`\mathrm{\Omega }=\frac{\stackrel{}{F}^2}{2m^2\beta ^2}+\frac{D}{m\beta }\stackrel{}{}\stackrel{}{F}`$.

Let us emphasize that it is the diffusion (Smoluchowski) approximation which makes the right-hand-side of Eq. (3) substantially different from the usual moment equations appropriate for the Brownian motion. In particular, the force $`\stackrel{}{F}`$ presumed to act upon an individual particle, does not give rise in Eq. (3) to the expression $`\frac{1}{m}\stackrel{}{}V`$ which might be expected on the basis of kinetic theory intuitions and moment identities directly derivable from the Karmers equation, but to the term $`+\stackrel{}{}\mathrm{\Omega }`$.

Moreover, instead of the standard pressure term, there appears a contribution from a probability density $`\rho `$-dependent potential $`Q(\stackrel{}{x},t)`$. It is given in terms of the so-called osmotic velocity field $`\stackrel{}{u}(\stackrel{}{x},t)=D\stackrel{}{}ln\rho (\stackrel{}{x},t)`$, (cf. ): $`Q(\stackrel{}{x},t)=\frac{1}{2}\stackrel{}{u}^2+D\stackrel{}{}\stackrel{}{u}`$ and is generic to a local momentum conservation law respected by isothermal Markovian diffusion processes, cf. -.

The Smoluchowski drift does not refer to any flows in the noise carrier. To analyze perturbations of the medium and then the resulting intrinsic (mean) flows, a more general function $`\stackrel{}{b}(\stackrel{}{X}(t),t)`$, must replace the Smoluchowski drift in Eqs. (1), (2). Forward drifts modify additively the pure noise (Wiener process entry) term in the Itô equations. Under suitable restrictions, we can relate probability measures corresponding to different (in terms of forward drifts !) Fokker-Planck equations and processes by means of the Cameron-Martin-Girsanov theory of measure transformations. The Radon-Nikodym derivative of measures is here involved and for suitable forward drifts which are are gradient fields that yields, , the most general form of an auxiliary potential $`\mathrm{\Omega }(\stackrel{}{x},t)`$ in Eq. (3):

$$\mathrm{\Omega }(\stackrel{}{x},t)=2D[_t\varphi +\frac{1}{2}(\frac{\stackrel{}{b}^2}{2D}+\stackrel{}{}\stackrel{}{b})].$$
$`(4)`$   
We denote $`\stackrel{}{b}(\stackrel{}{x},t)=2D\stackrel{}{}\varphi (\stackrel{}{x},t)`$.

Eq. (4) is a trivial identity, if we take for granted that all drifts are known from the beginning, like in case of typical Smoluchowski diffusions where the external force $`\stackrel{}{F}`$ is a priori postulated. We can proceed otherwise and, on the contrary, one can depart from a suitably chosen space-time dependent function $`\mathrm{\Omega }(\stackrel{}{x},t)`$. From this point of view, while developing the formalism, one should decide what is a quantity of a primary physical interest: the field of drifts $`\stackrel{}{b}(\stackrel{}{x},t)`$ or the potential $`\mathrm{\Omega }(\stackrel{}{x},t)`$.

Mathematical features of the formalism appear to depend crucially on the properties (like continuity, local and global boundedness, Rellich class) of the potential $`\mathrm{\Omega }`$, see e.g. . Let us consider a bounded from below (local boundedness from above is useful as well), continuous function $`\mathrm{\Omega }(\stackrel{}{x},t)`$, cf. ). Then, by means of the gradient field ansatz for the diffusion current velocity ($`\stackrel{}{v}=\stackrel{}{}S_t\rho =\stackrel{}{}[(\stackrel{}{}S)\rho ]`$) we can transform the momentum conservation law (3) of a Markovian diffusion process to the universal Hamilton-Jacobi form:

$$\mathrm{\Omega }=_tS+\frac{1}{2}|\stackrel{}{}S|^2+Q$$
$`(5)`$   
where $`Q(\stackrel{}{x},t)`$ was defined before. By applying the gradient operation to Eq. (5) we recover (3). In the above, the contribution due to $`Q`$ is a direct consequence of an initial probability measure choice for the diffusion process, while $`\mathrm{\Omega }`$ via Eq. (4) does account for an appropriate forward drift of the process.

Thus, in the context of Markovian diffusion processes, we can consider a closed system of partial differential equations which comprises the continuity equation $`_t\rho =\stackrel{}{}(\stackrel{}{v}\rho )`$ and the Hamilton-Jacobi equation (5), plus suitable initial (and/or boundary) data. The underlying diffusion process is specified uniquely, cf. .

Since the pertinent nonlinearly coupled system equations looks discouraging, it is useful to mention that a linearisation of this problem is provided by a time-adjoint pair of generalised diffusion equations in the framework of the so-called Schrödinger boundary data problem, . The standard heat equation appears as a very special case in this formalism.

The local conservation law (3) acquires a direct physical meaning (the rate of change of momentum carried by a locally co-moving with the flow volume, ), only if averaged with respect to $`\rho (\stackrel{}{x},t)`$ over a simply connected spatial area. If $`V`$ stands for a volume enclosed by a two-dimensional outward oriented surface $`V`$, we define a co-moving volume on small time scales, by deforming the boundary surface in accordance with the local current velocity field values. Then , let us consider at time $`t`$ the displacement of the boundary surface $`V(t)`$ defined as follows: $`\stackrel{}{x}V\stackrel{}{x}+\stackrel{}{v}(\stackrel{}{x},t)\mathrm{}t`$ for all $`\stackrel{}{x}V`$. Up to the first order in $`\mathrm{}t`$ this guarantees the conservation of mass (probability measure) contained in $`V`$ at time $`t`$ i. e. $`_{V(t+\mathrm{}t)}\rho (\stackrel{}{x},t+\mathrm{}t)d^3x_{V(t)}\rho (\stackrel{}{x},t)d^3x0`$.

The corresponding (to the leading order in $`\mathrm{}t`$) quantitative momentum rate-of-change measure reads, cf. , $`_V\rho \stackrel{}{}(\mathrm{\Omega }Q)d^3x`$.

For a particular case of the free Brownian expansion of an initially given $`\rho _0(\stackrel{}{x})=\frac{1}{(\pi \alpha ^2)^{3/2}}exp(\frac{x^2}{\alpha ^2})`$, where $`\alpha ^2=4Dt_0`$, we would have $`_V\rho \stackrel{}{}Qd^3x=_VP𝑑\stackrel{}{\sigma }`$, where $`Q(\stackrel{}{x},t)=\frac{\stackrel{}{x}^2}{8(t+t_0)^2}\frac{3D}{2(t+t_0)}`$, while the ”osmotic pressure” contribution reads $`P(\stackrel{}{x},t)=\frac{D}{2(t+t_0)}\rho (\stackrel{}{x},t)`$ for all $`\stackrel{}{x}R^3`$ and $`t0`$.

The current velocity $`\stackrel{}{v}(\stackrel{}{x},t)=\stackrel{}{}S(\stackrel{}{x},t)=\frac{\stackrel{}{x}}{2(t+t_0)}`$ is linked to the Hamilton-Jacobi equation $`_tS+\frac{1}{2}|\stackrel{}{}S|^2+Q=0`$ whose solution reads: $`S(\stackrel{}{x},t)=\frac{\stackrel{}{x}^2}{4(t+t_0)}+\frac{3}{2}Dln[4\pi D(t+t_0)]`$.

Let us observe that the initial data $`\stackrel{}{v}_0=D\stackrel{}{}ln\rho _0=\stackrel{}{u}_0`$ for the current velocity field indicate that we have totally ignored a crucial preliminary stage of the dynamics on the $`\beta ^1`$ time scale, when the Brownian expansion of an initially static ensemble has been ignited and so particles have been ultimately set in motion.

Notice also that our ”osmotic expansion pressure” $`P(\stackrel{}{x},t)`$ is not positive definite, in contrast to the familiar kinetic theory (equation of state) expression for the pressure $`P(\stackrel{}{x})=\alpha \rho ^\beta (\stackrel{}{x}),\alpha >0`$ appropriate for gases. The admissibility of the negative sign of the ”pressure” function encodes the fact that the Brownian evolving concentration of particles generically decompresses (blows up), instead of being compressed by the surrounding medium.

The loss (in view of the ”osmotic” migration) of momentum stored in a control volume at a given time, may be here interpreted in terms of an acceleration $`_V\rho \stackrel{}{}Qd^3x`$ induced by a fictituous ”attractive force”.

By invoking an explicit Hamilton-Jacobi connection (5), we may attribute to a diffusing Brownian ensemble the mean kinetic energy per unit of mass $`_V\rho \frac{1}{2}\stackrel{}{v}^2d^3x`$. In view of $`<\stackrel{}{x}^2>=6D(t+t_0)`$, we have also $`_{R^3}\rho \frac{1}{2}\stackrel{}{v}^2d^3x=\frac{3D}{4(t+t_0)}`$. Notice that the mean energy $`_V\rho (\frac{1}{2}\stackrel{}{v}^2+Q)d^3x`$ needs not to be positive. Indeed, this expression identically vanishes after extending integrations from $`V`$ to $`R^3`$. On the other hand the kinetic contribution, initially equal $`_{R^3}\frac{1}{2}\rho v^2d^3x=3D/\alpha ^2`$ and evidently coming from nowhere, continually diminishes and is bound to disappear in the asymptotic $`t\mathrm{}`$ limit, when Brownian particles become uniformly distributed in space.

Normally, diffusion processes yielding a nontrivial matter transport (diffusion currents) are observed for a non-uniform concentration of colloidal particles which are regarded as independent (non-interacting). We can however devise a thought (numerical) experiment that gives rise to a corresponding transport in terms of an ensemble of sample (and thus independent) Brownian motion realisations on a fixed finite time interval, instead of considering a multitude of them (migrating swarm of Brownian particles) simultaneously.

Let us assume that ”an effort” (hence, an energy loss/gain which implies local deviations from thermal equilibrium conditions) of the random medium, on the $`\beta ^1`$ scale, to produce a local Brownian diffusion current $`(\rho \stackrel{}{v})(\stackrel{}{x},t_0)`$ out of the initially static ensemble and thus to decompress (lower the blow-up tendency) an initial non-uniform probability distribution, results in the effective osmotic reaction of the random medium. This is the Brownian recoil effect of Ref. .

In that case, the particle swarm propagation scenario becomes entirely different from the standard Brownian one. First of all, the nonvanishing forward drift $`\stackrel{}{b}=\stackrel{}{u}`$ is generated as a dynamical (effective, statistical here !) response of the bath to the enforced by the bath particle transport with the local velocity $`\stackrel{}{v}=\stackrel{}{u}`$. Second, we need to account for a parellel inversion of the pressure effects (compression $`+\stackrel{}{}Q`$ should replace the decompression $`\stackrel{}{}Q`$) in the respective local momentum conservation law.

Those features can be secured through an explicit realization of the action-reaction principle which we promote to the status of the third Newton law in the mean.

On the level of Eq. (3), once averaged over a finite volume, we interpret the momentum per unit of mass rate-of-change $`_V\rho \stackrel{}{}(\mathrm{\Omega }Q)d^3x`$ which occurs exclusively due to the Brownian expansion mechanism, to generate a counterbalancing rate-of-change tendency in the random medium. To account for the emerging forward drift and an obvious modification of the subsequent expansion of an ensemble of particles, we re-define Eq. (3) by setting $`_V\rho \stackrel{}{}(\mathrm{\Omega }Q)d^3x`$ in its right-hand-side instead of $`+_V\rho \stackrel{}{}(\mathrm{\Omega }Q)d^3x`$ . That amounts to an instantaneous implementation of the third Newton law in the mean (action-reaction principle) in Eq. (3).

Hence, the momentum conservation law for the process with a recoil (where the reaction term replaces the decompressive ”action” term) would read:

$$_t\stackrel{}{v}+(\stackrel{}{v}\stackrel{}{})\stackrel{}{v}=\stackrel{}{}(Q\mathrm{\Omega })$$
$`(6)`$   
so that

$$_tS+\frac{1}{2}|\stackrel{}{}S|^2Q=\mathrm{\Omega }$$
$`(7)`$   
stands for the corresponding Hamilton-Jacobi equation, cf. , instead of Eq. (5). A suitable adjustment (re-setting) of the initial data is here necessary.

In the coarse-grained picture of motion we thus deal with a sequence of repeatable scenarios realised on the Smoluchowski process time scale $`\mathrm{}t`$: the Brownian swarm expansion build-up is accompanied by the parallel counterflow build-up, which in turn modifies the subsequent stage of the Brownian swarm migration (being interpreted to modify the forward drift of the process) and the corresponding built-up anew counterflow.

The new closed system of partial differential equations refers to Markovian diffusion-type processes again, . The link is particularly obvious if we observe that the new Hamilton-Jacobi equation (7) can be formally rewritten in the previous form (5) by introducing:

$$\mathrm{\Omega }_r=_tS+\frac{1}{2}|\stackrel{}{}S|^2+Q$$
$`(8)`$   
where $`\mathrm{\Omega }_r=2Q\mathrm{\Omega }`$ and $`\mathrm{\Omega }`$ represents the previously defined potential function of any Smoluchowski (or more general) diffusion process.

It is $`\mathrm{\Omega }_r`$ which via Eq. (4) would determine forward drifts of the Markovian diffusion process with a recoil. They must obey the Cameron-Martin-Girsanov identity $`\mathrm{\Omega }_r=2Q\mathrm{\Omega }=2D[_t\varphi +\frac{1}{2}(\frac{\stackrel{}{b}^2}{2D}+\stackrel{}{}\stackrel{}{b})]`$.

Our new closed system of equations is badly nonlinear and coupled, but its linearisation can be immediately given in terms of an adjoint pair of Schrödinger equations with a potential $`\mathrm{\Omega }`$, . Indeed, $`i_t\psi =D\mathrm{}\psi +\frac{\mathrm{\Omega }}{2D}\psi `$ with a solution $`\psi =\rho ^{1/2}exp(iS)`$ and its complex adjoint makes the job, if we regard $`\rho `$ together with $`S`$ to remain in conformity with the previous notations. The choice of $`\psi (\stackrel{}{x},0)`$ gives rise to a solvable Cauchy problem. (Notice that by setting $`D=\frac{\mathrm{}}{2m}`$ we recover the standard quantum mechanical notation.)

This feature we shall exploit in below. Notice that, for time-indepedent $`\mathrm{\Omega }`$, the total energy $`_{R^3}(\frac{v^2}{2}Q+\mathrm{\Omega })\rho d^3x`$ of the diffusing ensemble is a conserved quantity.

The general existence criterions for Markovian diffusion processes of that kind, were formulated in Ref. , see also .

Let us consider a simple one-dimensional example. In the absence of external forces, we solve the equations (in space dimension one) $`_t\rho =(v\rho )`$ and $`_tv+(v)v=+Q`$, with an initial probability density $`\rho _0(x)`$ chosen in correspondence with the previous free Brownian motion example. We denote $`\alpha ^2=4Dt_0`$. Then, $`\rho (x,t)=\frac{\alpha }{[\pi (\alpha ^4+4D^2t^2)]^{1/2}}exp[\frac{x^2\alpha ^2}{\alpha ^4+D^2t^2}]`$ and $`b(x,t)=v(x,t)+u(x,t)=\frac{2D(\alpha ^22Dt)x}{\alpha ^4+4D^2t^2}`$ are the pertinent solutions. Notice that $`u(x,0)=\frac{2Dx}{\alpha ^2}=b(x,0)`$ amounts to $`v(x,0)=0`$, while in the previous free Brownian case the initial current velocity was equal to $`Dln\rho _0`$. This re-adjustment of the initial data can be interpreted in terms of the counterbalancing (recoil) phenomenon: the would-be initial Brownian ensemble current velocity $`v_0=u_0`$ is here completely saturated by the emerging forward drift $`b_0=u_0`$, see e.g. also . We deal also with a fictituous ”repulsive” force, which corresponds to the compression (pressure upon) of the Brownian ensemble due to the counter-reaction of the surrounding medium. We can write things more explicitly. Namely, now: $`Q(x,t)=\frac{2D^2\alpha ^2}{\alpha ^4+4D^2t^2}(\frac{\alpha ^2x^2}{\alpha ^4+4D^2t^2}1)`$ and the corresponding pressure term ($`Q=\frac{1}{\rho }P`$) reads $`P(x,t)=\frac{2D^2\alpha ^2}{\alpha ^4+4D^2t^2}\rho (x,t)`$ giving a positive contribution $`+Q`$ to the local conservation law. The related Hamilton-Jacobi equation $`_tS+\frac{1}{2}|S|^2=+Q`$ is solved by $`S(x,t)=\frac{2D^2x^2t}{\alpha ^4+4D^2t^2}Darctan(\frac{2Dt}{\alpha ^2})`$. With the above form of $`Q(x,t)`$ one can readily check that the Cameron-Martin-Girsanov constraint euqation for the forward drift of the Markovian diffusion process with a recoil is automatically valid for $`\varphi =\frac{1}{2}ln\rho +S`$: $`2Q=2D[_t\varphi +\frac{1}{2}(\frac{b^2}{2D}+b)]`$.

In anology with our free Brownian motion discussion, let us observe that presently $`<x^2>=\frac{\alpha ^2}{2}+\frac{2D^2t^2}{\alpha ^2}`$. It is easy to demonstrate that the quadratic dependence on time persists for arbitrarily shaped initial choices of the probability distribution $`\rho _0(x)>0`$. That signalizes an anomalous behaviour (enhanced diffusion) of the pertinent Markovian process when $`\mathrm{\Omega }=0`$ i. e. $`\mathrm{\Omega }_r=2Q`$.

We can evaluate the kinetic energy contribution $`_R\rho \frac{v^2}{2}𝑑x=\frac{4D^4t^2}{\alpha ^2(\alpha ^4+4D^2t^2)}`$ which in contrast to the Brownian case shows up a continual growth up to the terminal (asymptotic) value $`\frac{D^2}{\alpha ^2}`$. This value was in turn an initial kinetic contribution in the previous free Brownian expansion example. In contrast to that case, the total energy integral is now finite (finite energy diffusions of Ref. ) and reads $`_R(\frac{1}{2}v^2Q)\rho 𝑑x=\frac{D^2}{\alpha ^2}`$ (it is a conservation law). The asymptotic value of the current velocity $`v\frac{x}{t}`$ is twice larger than this appropriate for the Brownian motion, $`v\frac{x}{2t}`$.

It is easy to produce an explicit solution to (7), (8) in case of $`\mathrm{\Omega }(x)=\frac{1}{2}\gamma ^2x^2D\gamma `$, with exactly the same inital probability density $`\rho _0(x)`$ as before. The forward drift of the corresponding diffusion-type process does not show up any obvious contribution from the harmonic Smoluchowski force. It is completely eliminated by the Brownian recoil scenario. One may check that $`b(x,0)=\frac{2Dx}{\alpha ^2}=u(x,0)`$, while obviously $`b=\frac{F}{m\beta }=\gamma x`$ would hold true for all times, in case of the Smoluchowski diffusion process. Because of the harmonic attraction and suitable initial probability measure choice, we have here wiped out all previously discussed enhanced diffusion features. Now, the dispersion is attentuated and actually the non-dispersive diffusion-type process is realised: $`<x^2>`$ does not spread at all despite of the intrinsically stochastic nature of the dynamics (finite-energy diffusions of Ref. ).

It is clear that stationary processes are the same both in case of the standard Brownian motion and the Brownian motion with a recoil. The respective propagation scenarios substantially differ in the non-stationary case only.

One may possibly argue, that a sign inversion of the right-hand-side of Eq. (3) which implies (6), can be accomplished by means of an analytic continuation in time (cf. the Euclidean quantum mechanics discussion in Ref. ). An examination of free dynamics cases given in the above proves that the third Newton law ansatz is an independent procedure.

Acknowledgements: I would like to thank Professor Eric Carlen for discussion about the scaling limits of the Boltzmann equation and related conservation laws.