1. Introduction
In a recent article [
1], we have shown that in a perfect fluid the Bohm potential can be obtained without any quantum mechanical prior assumption. In fact, unlike previous works (see the comprehensive review of Dürr et al. [
2]) where the Madelung equation is derived starting from the Schrödinger equation, here the same result is obtained assuming a classical behavior, with a non-local logarithmic dependence of the free energy from the mass (or number) density. Let us start by reviewing this earlier work.
Consider the Hamilton–Jacobi equation for a particle of unit mass m = 1 located at position
and time
t in a potential energy field,
,
where
S is the action, defined as the time integral of the Lagrangian function. Here,
S depends on position
, time
t, and momentum
, i.e.,
, where
, showing that the momentum (and the velocity
as well) can be expressed as the gradient of a scalar function. Note that by taking the gradient of Equation (
1) we obtain Newton’s equation of motion:
denoting the material derivative.
Now, let us interpret the action
S as the velocity potential of an inviscid fluid (in fact, its momentum is
), assuming that this particle is one of the very many identical particles, not interacting with each other, that constitute an ideal fluid in isothermal conditions. Then, the governing equations, expressing the conservation of mass and momentum, read [
3]
where
is the mass density, while
is the force exerted on a particle, defined in Equation (
2). Note that
denotes both the mass density and the number density since, without loss of generality, we have assumed that particles have unit mass.
It should be stressed that
and
are mean values, defined within an elementary point volume, assuming local equilibrium. From a different point of view, we are neglecting fluctuations, i.e., the Ginzburg criterion is satisfied, so the mean field approximation can be applied (see Landau et al. [
4]).
Classically, in Equation (
4), the force
, denoted as
, is often expressed as the sum of a pressure term plus an external conservative force, i.e.,
where
P is the pressure which, for an ideal fluid, is readily written as
, so that
This expression, showing that the classical potential energy is the sum of a pressure-dependent inner part and an external part, assumes that the gas evolves under isothermal conditions. If we remove this assumption, an energy conservation equation should be added to Equations (
3) and (
4), as will be shown in future works.
Now, define a wave-like complex function,
,
where
,
, and
h is the Planck constant, imposing that both its real and imaginary components satisfy the Schrödinger equation,
with
denoting the classical potential energy (
6) acting upon the particle. The real part of Equation (
8) coincides with Equation (
3), showing that
can be interpreted as the probability to find a particle at location
and at time
t. In addition, the imaginary part of Equation (
8) yields the Hamilton–Jacobi Equation (
1), where
is the Bohm potential [
5,
6] that is a central concept of the de Broglie–Bohm formulation of quantum mechanics (see reviews in [
2,
7]). Therefore, we conclude that the classical Hamilton–Jacobi Equation (
1) is equivalent to the Schrödinger Equation (
6), provided that
, that is, the potential energy is the sum of the classical potential energy and its quantum counterpart, i.e., the Bohm potential. To observe this, first consider that when the Hamiltonian
H does not depend on time explicitly, then, applying separation of variables, we find that
, where
E is the energy, which is constant when the system is isolated or at steady state. Here,
W is called Hamilton’s characteristic function, and the quantum Hamilton–Jacobi equation becomes
At steady state, momentum is nul, that is
. Then,
W is uniform (we may assume it to be zero, as it is determined within an arbitrary constant), so that
and therefore Equation (
10) reduces to
. Then, considering Equation (
9) and denoting
, we find
which coincides with the steady solution of the Schrödinger equation.
It should be observed that the Bohmian formulation, involving a second-order differential equation of motion in the presence of the quantum potential, reproduces the classical predictions of quantum mechanics only if the initial position
is distributed according to the Born rule, i.e., according to the density
and if the particle momentum is restricted to the value
[
8].
Finally, summarizing, from Equations (
2) and (
4), the force balance on each particle yields
These results were obtained in 1927 by Madelung [
9], showing that the Schrödinger equation for one-electron problems can be transformed into hydrodynamical equations of an ideal (i.e., non-viscous) gas in isothermal conditions, subjected to the action of both classical and quantum potentials.
2. The Quantum Potential
Although many researchers have stated that the Bohm quantum potential is due to non-local effects, there is no clear explanation of how that happens. Recently [
1], recalling that, as
is an additive integral of motion it must be proportional to the energy (see Landau et al. [
4]), we propose the following non-local expression for the free energy per unit mass,
:
Here,
is the potential energy per unit mass resulting from the action of an external conservative force, while
is an interaction kernel between particles located at a distance
, with the normalization condition
. Dropping for convenience the time dependence, expanding
in Taylor series,
and truncating the series after the second term (there can be no
term, due to the isotropy of the fluid), we find
Here, the first term on the RHS is the usual classic (i.e., thermodynamic) free energy (per unit mass) of an ideal gas, coinciding with the classical potential energy (
6),
while
is the non-local part, with
denoting the square of a characteristic length,
a. Note the negative sign in Equation (
18), revealing that particles are assumed to attract each other. Higher-order terms in the expansion (
14) can be neglected under quite general conditions, as discussed in Sections 3 and 4 of Ref. [
1]. In fact, assuming an exponentially decreasing interaction kernel
in [
1], we show that the additional terms, correcting the Madelung equation, are very small.
Now, define the free energy per unit volume,
, so that the total free energy is given by the following functional:
where
V is the total volume, that we assume to be infinite. Integrating parts with the assumption that
exponentially as
, we find
where
f denotes the effective free energy per unit volume.
Imposing that the free energy is minimal, under the constraint of mass conservation, i.e.,
, we obtain the following generalized chemical potential:
where
and
Here,
is the thermodynamic chemical potential (an energy per unit mass) which, apart from an irrelevant additive constant (i.e.,
), coincides with the thermodynamic free energy (i.e., the classical potential energy), while
is the non-local contribution. After a straightforward calculation (see Ref. [
1]), it can be shown that the non-local potential (
23) can be expressed in the following equivalent form:
In particular, imposing that at leading order the non-local potential (
24) reduces to Bohm’s quantum potential (
9), i.e.,
then
a coincides with the thermal de Broglie wavelength,
The equations of motion can be determined by applying a variational principle [
10] to derive the Euler Equation (
2) for a compressible, inhomogeneous fluid where
is a body force, driven by density gradients in the fluid. The same result can be obtained more heuristically by applying Noether’s theorem [
11] or by the least-action principle [
12]. At the end, we find Equations (
3) and (
4), with the following force per unit mass:
Let us consider the two terms on the RHS separately. On one hand, we have
and therefore, considering that
,
On the other hand, from Equations (
23) and (
26), we define a non-local reversible body force, which is usually referred to as the Korteweg force [
10]:
Finally, summarizing, from Equation (
2), the momentum balance per unit mass yields the Madelung equation:
Equation (
31) can also be written as a quantum Bernoulli equation as follows [
13]:
where the RHS denotes the gradient of the total energy,
E, that is the sum of kinetic energy, potential energy, and quantum potential energy.
It should be stressed that the Madelung Equation (
31) (as well as the quantum Bernoulli Equation (
32)) is a leading-order approximation of the more complete equation of motion where the full gradient expansion (
14) is accounted for (see Equation (
38) of Ref. [
1]).
3. Retarded Potential and Covariant Formulation
The non-local constitutive Equation (
13) implies that density changes in
determine an instantaneous variation of the free energy at
, with an infinite propagation velocity. To correct this obviously unrealistic assumption, Equation (
13) is replaced with the following expression:
where
is a retarded time and
c denotes the propagation velocity.
Expanding
in the Taylor series,
and truncating the series after the quadratic term (the linear terms give no contribution due to simple spatial and temporal symmetry considerations), we find Equation (
15), i.e.,
, where
is again the classical potential energy (
16), while the non-local contribution reads
with
denoting the d’Alambertian wave operator. Equation (
37) is the Lorentz-invariant form of the non-local energy (
17). This result has been obtained previously by Holland [
14] or, more recently, by Nicolić [
15,
16] and Rahmani et al. [
17], using a variational approach.
Consider a flat Lorentz-invariant Minkowski space-time,
, assuming the invariance of the element
, where
is the Minkowski metric, i.e.,
. Consequently, defining the covariant and contravariant gradient operators [
18],
shows that the d’Alambertian corresponds indeed to a four-dimensional Laplacian, i.e.,
Here, the relative time t is related to the proper time (i.e., measured in the co-moving reference frame) as , where , with and is the 3D velocity. Note that should not be confused with the 4-velocity , with .
Finally, proceeding as in Equations (
19)–(
25), we find the chemical potential to be
, where
is unchanged from the non-relativistic case, while
which generalizes Equation (
25).
Now we can write the governing equations. First, the continuity Equation (
3) in covariant form reads
where
is the 4-flux (that is, the 4-momentum per unit volume). Additionally,
, where
is the velocity (that is, the 4-momentum per unit mass).
As for the momentum balance equation, the equation of motion for an ideal fluid, that is free of any dissipative effects, can be obtained assuming that Equations (
4) and (
31) must be reduced in the non-relativistic limit, obtaining, refs. [
18,
19,
20]
where
with
, and
is the energy-momentum 4-tensor, while
is the force density. Finally, the following governing equation is obtained:
where
h is the enthalpy per unit volume, i.e.,
Thus, in agreement with Tavernelli [
21], we see that particles move deterministically along geodesic lines in a curved space whose curvature depends nonlocally on Bohm’s quantum potential.