3. The Thermodiffusion and the Soret Equilibrium in the DML
In the first part of this section, we will shortly gather the main and relevant aspects of the DML avoiding repeating the concepts introduced and extensively discussed elsewhere. Readers interested in the details of the model are invited to refer to the literature [
17,
18,
19] where the theoretical developments, numerical simulations, experimental results, and applications are collected and discussed in depth. The second part is dedicated to showing that the DML may provide a completely new approach to the physical interpretation of the thermodiffusion phenomenon and of the Soret equilibrium occurring in liquid mixtures to which a temperature gradient is applied (when the phenomenon involves macromolecular solutes, it is generally named thermophoresis). A similar approach can be applied of course also to the Dufour effect, as well as to other phenomena typical of NET. In particular, we will show how the DML is able to model the thermodiffusion, taking into account even the sign variation that sometimes occurs for the same couple of solute-solvent of a mixture at varying parameters, such as the average temperature, the intensity of the temperature gradient or also the initial concentration. The capability of the DML to justify that the third law of dynamics is not falsified in such type of experiments, as some authors have supposed ([
12,
13,
14] and references therein), is dealt with in the
Section 5.
Two facts make challenging the experimental determination of the Soret coefficient, namely its small magnitude, usually of the order of
, and the presence of undesirable convective disturbances due to the coupling of the temperature gradient with gravity, sometimes even present in microgravity environments [
66,
68,
69,
70,
71,
72,
73,
74,
75,
76,
77]. It is then not surprising that the several theoretical models proposed through the years [
12,
13] often fail when dealing with peculiar characteristics of
, such as its sign inversion or the presence of minima. Experimental [
38,
39,
40,
41,
42,
43,
44,
45,
46], theoretical [
36,
45,
46,
47,
48,
49,
50,
51,
52,
53,
54,
55,
56,
57,
58,
59,
60,
61,
62,
63,
64] and computational [
12,
13,
14] studies on thermodiffusion and Soret equilibrium have significantly advanced in recent years, but several intriguing and fundamental questions still remain unanswered. One such question, maybe the most relevant, is the origin of the forces driving the observed phenomenology at the micro-mesoscopic scale. We will try to fill this gap by means of DML.
The DML rests on two hypotheses, both with an experimental background, namely: (1) experiments performed by means of Inelastic X-ray Scattering (IXS) and Inelastic Neutron Scattering (INS) techniques [
78,
79,
80,
81,
82,
83,
84,
85,
86,
87,
88,
89,
90,
91,
92,
93,
94,
95,
96,
97] allowed to discover that the mesoscopic structure of liquids is characterized by the ubiquitous presence of solid-like structures, whose size is of a few molecular diameters, through which elastic waves propagate as in the corresponding solid phase; the number and the size of these structures depend upon the temperature and pressure of the medium. (2) Elastic energy and momentum in liquids propagate by means of collective oscillations, or wave-packets, similarly to phonons in crystalline solids (this hypothesis is a consequence of the first). Although both arguments have been largely and deeply dealt with in [
17] where DML is introduced, it is the case to recall the attention of the reader on the relevance of the INS and IXS experiments findings, whose results are at the basis of a dual vision of the liquid modeling at the mesoscopic scale. In fact, even if the frequency ranges of these experiments are high, typically from less than 1 to 2–3 THz, however, the fact that the propagation speed of elastic waves in the pseudo-crystalline structures of liquids are higher than those typical of liquids, and closer to those of solids, ensures that these speeds are maintained over relatively long distances, covering several molecular diameters, i.e., the dimensions of the clusters (see also the
Section 6 for further recent findings on this topic). In fact, the experimental results show that, once these distances are exceeded, the propagation speeds of the elastic waves return to the range of values typical of liquids, thus providing a signature of the borders of the pseudo-crystalline structures. Among the many experimental works performed on such topic, other than that of Ruocco et al. [
81] who measured the speed of sound in water amounting to
on length scale of few molecular diameters, it is the case of citing that of Giordano and Monaco [
97]. They performed the first measurements across the glass-liquid and solid-liquid transitions and reported on a comparison of the collective excitations in liquid and polycrystalline sodium. As it concerns liquid sodium, it exhibits acoustic excitations of both longitudinal and transverse polarization at frequencies strictly related to those of the corresponding crystal. The relevant difference between the liquid and the crystal is the broadening of the excitations in the case of liquid because of an additional disorder-induced contribution coming into play. The authors deduce a direct connection between the structural and dynamic properties of liquids, with short-range order and overall structural disorder characterized by specific fingerprints. To confirm the validity of the DML, we evaluated the OoM of the phonon mean free path
and lifetime
within a
liquid particle [
17]. With reference to the experimental values,
will be a multiple of the phonon wavelength
,
, and
of
,
, with
. Using the data for the water of [
79,
94], typical values for the parameters characterizing a phonon (variation range is function of temperature, pressure, and
orientation) are: (central) frequency
, wave-length
and velocity
. Interestingly, this value fits very well with the experimental data obtained for the propagation velocity of thermal waves in water [
81]. The value of
represents also the typical OoM of the size of a
liquid particle at the exchanged momentum
, allowing us to confidentially assume a value of few units for
. A theoretical explanation of the occurrence of a solid-like value of sound propagation in liquids is based on the points of view of Brillouin and Frenkel. We may distinguish two limit regimes for the propagation of perturbations in liquids. One is the “normal” hydrodynamic regime, in which one considers the medium as a continuum and the excitations on a so long timescale that the system may be assumed to be in thermodynamic equilibrium. This regime is also usually referred to as viscous. The second is the “solid-like” or elastic regime, in which the dynamics becomes that of a free particle between successive elastic collisions. The investigation of collective dynamics becomes of particular interest when intermediate time- and length-scale are considered, i.e., distances comparable to those characterizing the structural correlations among particles, and times comparable to the lifetimes of such correlations. Although these two extreme behaviors, viscous and elastic, are well known in physics, the intermediate situation is still far from being fully characterized. This intermediate range is referred to as viscoelastic regime. It is just in this regime that the DML is applied, as is clear also from the characteristics of the elementary interactions described in the following.
According to this view, thermal energy is considered a form of elastic energy, as supposed by Debye [
98,
99], Brillouin [
100,
101], and Frenkel [
102]. At the mesoscopic level the liquid phase is modeled in the DML as constituted by two sub-systems, mutually interacting among them, the
liquid particles, i.e., a sort of solid-like aggregates of liquid molecules, and the wave-packets, or
lattice particles. Therefore we are led to define a “liquid” in the DML not as a liquid in the classical assumption, but as a mixture of solid icebergs (the solid-like part) and of an amorphous phase (the classical liquid part). It will not have escaped a very attentive and informed reader that this approach is not completely new in the world of the physics of liquids. First, it recalls that of Brillouin [
17,
100,
101] and Frenkel [
102]; Landau also described the HeII as made of “normal” and “superfluid” parts [
103]. Eyring in his Reaction Rate Theory [
104,
105,
106], also known as the “hole” theory, had supposed that a liquid could be assumed to have a quasi-crystalline structure. Similarly to a gas, the liquid should consist of molecules moving around empty spaces, or holes. A molecule is pictured as vibrating around an equilibrium position until it acquires sufficient energy to overtake the attracting forces “and” a hole is available in the nearest neighbor where the molecule can jump (see also Figure 2 in [
17]). More recently, Andreev [
107] in a very interesting paper assumed that a liquid consists of two weakly coupled systems, the phonons and the remainder of the liquid. He calculates the effects connected with the presence of weakly damped phonons in a normal liquid. In particular, he calculated the expressions for mechano-thermal and thermo-mechanical effects, as well as the propagation of shear oscillations (see also the
Section 6 for recent experimental findings on this topic).
Figure 2 shows in a pictorial way the model we have in mind. Solid-like icebergs, the
liquid particles, are dispersed in the amorphous phase of liquid, while thermo-elastic waves, the
lattice particles, travel through them. As far as the perturbations travel inside a solid-like structure, they show a solid-like character, i.e., they are (quasi) harmonic waves traveling at speeds close to that of the corresponding solid phase (3200 m/s in the case of water, [
17,
18,
81]). When however they leave the
liquid particle and travel through the amorphous environment, lose the harmonic outfit and become wave-packets. The interaction itself between these waves and the solid-like icebergs has an anharmonic character, allowing momentum and energy to be exchanged between the two reservoirs, the
liquid particles and the
lattice particles. The momentum transferred
generates the force
responsible for the energy and mass diffusion in liquids [
17,
18,
19];
represents the (finite) duration of the
lattice particle liquid particle interaction and
is the anharmonic interaction potential between the
lattice particle and the
liquid particle, given by [
17,
19] (for a discussion about the functional dependence of
upon distance and about its quantum/classical form, the reader may refer to the Discussion in [
17]):
Here
is the phase velocity associated with the wave-packet and
is the cross-section of the “obstacle”, the solid-like cluster, on the surface of which
is applied;
the energy flux carried by wave-packets [
17]. Some considerations on
, or
, are in order. First,
is supposed to be anharmonic because of the inelastic character of the interaction, which is not instantaneous but lasts
during which the particle is displaced by
(here and in the rest of the paper, the two brackets
indicate the average over a statistical ensemble of the quantity inside them). Second,
can either be positive or negative depending on whether the quantity
increases or decreases as a consequence of the interaction. This point is connected with the last equality in Equation (10) which deserves a more accurate analysis because it is pivotal for the purposes of the present manuscript. It is indeed well known that a wave impinging on an obstacle exerts on it a pressure. The radiation pressure was introduced for the first time by Balfour Stewart in 1871 [
20] for the case of thermal radiation. Maxwell took up this theme in Arts. 792–793 of his Treatise in 1873 [
21], arguing on the basis of his stress tensor. The first convincing experimental evidence for the radiation pressure of light was given by Lebedev in 1901 [
22] (the Crookes radiometer does not demonstrate electromagnetic radiation pressure). We have shown [
17] by means of the radiant vector
that this concept may be generalized and used for any type of wave (
is the oscillation velocity of the particles,
the unit vector along the direction of propagation and
the flux of elastic energy generating the pressure
). Here we exploit again this concept and apply it to the momentum transferred by elastic, or thermal, waves that impinge on an obstacle, such as for instance a solute molecule in a liquid mixture. Accordingly, we replace the gradient of elastic pressure with that of thermal energy, yielding:
Equation (12) contains two supplementary information; first, because thermal energy is carried by elastic wave-packets, the velocity of energy propagation remains the same. Second, the gradient of pressure gives rise to an inertial term,
ψth, dimensionally a force per unit of volume, responsible for momentum transport associated with the wave-packets propagation. Equation (12) is valid when the characteristics of the medium change continuously along the wave-packets propagation direction. The same argument may be used when wave-packets impinge on an “obstacle”, as for instance a
liquid particle, either of solute or of solvent: the radiant vector changes and a thermal radiation pressure Π
th is produced on the boundary. Observing that the inelastic character of the collision allows the exchange of momentum, this leads to the appearance of the force
fth, acting on molecular clusters. By rewriting Equation (12) in terms of discrete quantities, one yields the following expression for the net Π
th:
We conclude that the last equality in Equation (10) provides the algebraic sign for
according to the sign of the quantity in square brackets of Equation (14), i.e., on whether the energy associated with the propagation of elastic waves in the medium “1” is larger or smaller than that in the medium “2”. Expressions analogous to Equation (13) or Equation (14) have been previously derived by many authors following different approaches [
17,
20,
21,
22,
23,
24,
25,
26]. It derives from the Boltzmann–Ehrenfest’ Adiabatic Theorem [
27], later generalized by Smith [
28] and by Gaeta et al. [
29]; another approach gave analogous results by calculating the radiation pressure produced by acoustic waves on an idealized liquid-liquid interface [
30].
Let’s now go a little bit more into detail on the characteristics of
in the DML, i.e., on its intensity, orientation, functional dependence, etc. Everything begins by remembering that
is the internal energy per unit of volume of a liquid at temperature
,
the medium density,
the specific heat at constant volume per unit mass,
the Debye temperature of the liquid at temperature
. Since this dynamics occurs mainly at high frequencies and involves only the DoF of the lattice, we introduce the parameter
to account for the average number of DoF participating in the collective dynamics. Defining
as the ratio between the number of collective DoF surviving at temperature
and the total number of available collective DoF, the fraction
of
:
accounts for that part of the thermal energy transported by the wave-packets (the definition of
is easily deduced). Of course, it holds
[
17].
is the number of wave-packets per unit of volume and
their average energy.
is propagated through the liquid by means of inelastic interactions between the
liquid particles and the
lattice particles, i.e., the wave-packets that transport the thermal or, generally, the elastic energy.
Figure 3 exemplifies the elementary physical mechanism behind the DML. The
lattice particle ↔
liquid particle interaction allows the exchange of both energy and momentum between the phonon and the
liquid particle. To make more intuitive the model, we have ideally divided the elementary
liquid particle ↔
lattice particle interaction into two parts: one in which the
lattice particle collides with the
liquid particle and transfers to it momentum and energy (both kinetic and potential) and the other in which the
liquid particle relaxes returning the energy to the thermal reservoir through a
lattice particle, alike in a tunnel effect. Considering an event of type (a) of
Figure 3, an energetic wave-packet collides with a
liquid particle that acquires the momentum
(Equation (16)) and the energy
:
In Equation (17)
and
represent the wave-packet (central) frequency [
17] before and after the collision, respectively. Part of
becomes the kinetic energy
acquired by the
liquid particle, and the remaining part becomes
, i.e., the potential energy of internal DoF of the solid-like cluster. As a consequence of the collision (a), the wave-packet loses energy and momentum, which are acquired by the cluster and converted into kinetic and potential energy. The kinetic energy is dissipated as friction against the liquid. The potential energy
is relaxed by the
liquid particle once the interaction is completed, its dissipation lasts
during which the
liquid particle travels by
, at the end of which the residual energy stored into internal DoF returns to the liquid reservoir. The total displacement of the particle,
, and the total duration of the process,
, are:
There are several interesting observations at this point that we want to propose to the reader’s attention. First, as pointed out before, the elementary interaction looks like a tunnel effect, inasmuch as the energy subtracted from the phonon’s reservoir returns to it a time interval
later and a step
forward. We want to clarify however that the purely classical “tunnel effect” occurring in the
liquid particle ↔
lattice particle interaction should not be confused with tunneling in physics, which is ordinarily associated with quantum phenomenology. We have adopted such a terminology only to give the reader the idea of what happens in a liquid upon the interaction, but there are no hidden quantum phenomena beyond such a word. What matters here to the scope of DML is the dynamics linked to the relaxation time, i.e., to recognize that
is the time interval during which the energy disappears from the liquid thermal reservoir because it is trapped in the internal DoF; once
has elapsed, it reappears in a different place. In this way, relaxation times, introduced ad hoc by Frenkel [
102], find a clear physical interpretation (interestingly,
and
in Equations (18) and (19) have the same meaning as
and
defined by Frenkel [
102], see also Figure 2 in [
17]). The second point recalls us to the mind that in this manuscript we are dealing with the Soret equilibrium, i.e., a phenomenon characterized by a dynamical equilibrium between the diffusion of particles due to the heat crossing the mixture, that is balanced by the Fick back diffusion due to the concentration gradient that is being established by the thermodiffusion. In the papers published so far [
17,
18,
19], we have mainly been concerned with that part of the energy lost by the wave-packets and converted into internal (potential) energy of
liquid particles,
, and with the consequences of its propagation. At equilibrium, the variation of potential energy
is zero in the average because the thermal content of the system does not vary anymore, as well as the number of DoF excited. Accordingly, the only term surviving in the third member of Equation (17) is that accounting for the average kinetic energy acquired by the
liquid particle,
. Apart from this, in the frame of DML, we are instead concerned with a flux of momentum, in fact, we are analyzing the mechanism responsible for the solute, or solvent, molecules flux pushed by a flux of thermal (elastic) energy crossing the liquid and due to an external source, the temperature gradient applied to the system (a liquid mixture). We are then concerned here with that part of energy lost by wave-packets and converted into kinetic energy of liquid particles.
We want now to draw the attention of the reader to the event of type (b) of
Figure 3 which is the time-reversal of (a): an energetic cluster interacts with a wave-packet transferring to it energy and momentum. The outgoing wave-packet has momentum and energy larger than the incoming one so that the net effect of the interaction is to increase the liquid thermal energy carried by phonons at the expense of the kinetic and internal energy of the liquid particle. Another crucial point relevant to our purposes is that at (thermal or concentration) equilibrium conditions, events (a) will alternate over time with events (b), to keep statistically equivalent the balance of the two energy reservoirs. The effect of the interactions is null if integrated over the entire surface of the particle since its thermal regime does not change, nor does it move; consequently the mesoscopic equilibrium induces the macroscopic equilibrium, following the Onsager postulate [
1,
2,
3,
17]: events like those of
Figure 3 will be equally probable along any direction, with a null average over time and space. On the contrary, when the symmetry is broken by an external force field, for instance by a temperature (Soret) or a concentration (Dufour) gradient, type (a) events will prevail over type (b), or vice-versa, along the direction of the gradient; the gradient will generate an imbalance of the collisions and possibly also a variation of the energy of the phonons (their number does not change as a first approximation) and an imbalance of the number of collisions. We will return to this point when discussing the non-equilibrium phenomena which are the topic of the present manuscript.
Figure 4 represents a close-up of
Figure 3a and shows how the energy
is propagated inside a
liquid particle until it is given back to the thermal reservoir. The interaction begins with a wave-packet impinging on a
liquid particle, transferring to it energy and momentum.
is the extension of the wave-packet and
that of the liquid particle. Because the
liquid particle in the DML is considered a solid-like structure, following the experimental evidences obtained by means of IXS and INS experiments [
17,
78,
79,
80,
81,
82,
83,
84,
85,
86,
87,
88,
89,
90,
91,
92,
93,
94,
95,
96,
97] and the pioneering intuitions of Brillouin [
101] and Frenkel [
102], the phonon’s flight through the
liquid particle behaves like a (quasi-) harmonic wave, crossing it at a speed close to that of the corresponding solid phase. Once
has elapsed, the
liquid particle has travelled by
, it relaxes the energy stored in the internal DoF, then it travels again by
during
(not shown in the figure). The presence of a relaxation time is a signature of liquids allowing to justify, among others, some of their peculiar characteristics [
17,
19].
A question that the reader is certainly wondering at this point is “what pushes and drives the wave-packets in liquids between two successive collisions?”. To answer such question we should take in mind that in DML the wave-packet current is always present in liquid, even in isothermal conditions. During their free-flight between two successive interactions (
Figure 2) with the
liquid particle, the wave-packets are considered in the DML as local microscopic heat currents. Each
lattice particle propagates along an average distance
, during a time interval
, which is its average life-time. Accordingly, the ratio
defines the wave-packet average propagation velocity between two successive interactions,
, where
and
are the (central) frequency and wave-length of the wave-packet (although an elastic perturbation of the liquid lattice should be represented by a wave-packet, whose localization depends on the microscopic structure to which it is associated, here it will be represented, for the sake of mathematical simplification, before and after an interaction, as a monochromatic wave, longitudinally polarized, and hence identified by its average wave-length
and frequency
, whose product is equal to
). Considering a surface S arbitrarily oriented within the liquid, at thermal equilibrium an equal number of wave-packets will flow through it in opposite directions. Let’s assume now that locally every heat current, although isotropically distributed in the liquid, is driven by a
virtual temperature gradient applied over
. From Equation (16), the variation of the energy density within the liquid along the “
” direction is:
where the local virtual temperature gradient
in the last member represents the thermodynamic “
generalized” force (with the meaning of NET) driving the diffusion of thermal excitations along
z. The factor
accounts for the isotropy of the phonon current through a liquid at equilibrium. By applying this driving force over the distance
, a wave-packet diffusion, i.e., the heat current
, is generated along the “
” direction:
where
is the wave packets diffusion coefficient along “
”. Therefore, the propagation of an elementary heat current lasts the time interval
and is driven by a virtual temperature gradient
over the distance
. Of course, in an isothermal medium the number of scatterings along any direction is balanced by that in the opposite one, to get a null diffusion over time. On the contrary, when an external thermodynamic (generalized) force is applied to the system, as for instance a temperature gradient, a heat flux (in the opposite direction) is generated, and therefore, an increase of the number of “
wave-packet ↔
liquid particle” interactions in the same direction. In this case (Soret effect), events of type (a) of
Figure 3 have a larger probability to occur than events of type (b). Analogously, if the symmetry is broken by a particle concentration gradient (Dufour effect), the events having a larger probability of occurrence will be the (b). To quantify such reasoning, let
be the temperature gradient externally applied to a system. As before, we assume here as a first approximation that the application of an external temperature gradient does not change the average total number
of the wave-packets, (small gradient assumption). Let then
be the average number per second of “
wave-packet ↔
liquid particle” collisions at
in the isothermal liquid; due to the presence of the virtual temperature gradient, which is equally present in all three directions, for each direction this number is
. If an external (small) temperature gradient
is applied to the liquid, the total number of collisions per second
is assumed to stay constant because we have supposed that
does not change, but for collisions occurring in the direction of heat propagation, there will be
excess of collisions per second and an equal defect
in the opposite direction. If the intensity of the temperature gradient is not too high, we may assume that there will be no non-linear effects and that the initial state of the system at the microscopic level locally continues to be similar to the one at uniform temperature, except for thermal excitations along
. Therefore a
liquid particle at a given place in the liquid experiences the same number of collisions per second as if the temperature was uniform; the only difference with the isothermal case is that there is an imbalance in the number of collisions with wave-packets that originated in the two half-spaces along z, one hotter and the other cooler. In other words, while phonons in liquids are pushed by the local virtual temperature gradient
, the external temperature gradient
has the duty of orienting them along its own direction. In the linear range, we may assume that
is proportional to the temperature gradient
externally applied to the liquid. Consequently, there is a heat flux
superimposed to two equal and opposite heat fluxes
and
(
represent the thermal conductivity of the macroscopic system,
is that related to the collective DoF, [
17,
19]) giving:
Accordingly, an external temperature gradient
increases by
the phonon flux with respect to that corresponding to the local microscopic heat currents
due to spontaneous phonon diffusivity. Equation (22) gives the ratio of the number of “
wave-packet ↔
liquid particle” collisions per second due to the applied gradient, to that due to the random motions of collective thermal excitations. The clusters therefore will experience
collisions in excess along the direction of the heat flux and will execute as many jumps per second of average length
in excess in the same direction. Consequently, every particle travels the distance
per second along the direction of heat flux due to the external temperature gradient; this represents the drift velocity
of the
liquid particle along
due to the external temperature gradient:
Equation (23) allows defining the thermal diffusion coefficient
that we assume in DML to be the equivalent of
defined in NET (see Equation (6)) that by definition is the drift velocity of molecules in a unitary temperature gradient:
Because
is a signature of the “solute + solvent” mixture, it is correctly given in Equation (24) in terms of the
liquid particle drift
caused by the collisions
with the
lattice particle. Even more relevant to the purpose of the manuscript is its (reciprocal) dependence upon
. In fact
is the only term in Equation (24) that intrinsically contains an algebraic sign, determining in such a way the possibility for
to change sign when some particular conditions occur in a liquid mixture. The expression for the Soret coefficient
is easily obtained from its definition given by the ratio between the thermodiffusion coefficient and that of ordinary diffusion, for which we exploit the well-known Einstein relation (as specified before, we suppose the mixture sufficiently dilute so that the motion of molecules, either of solute or of solvent, does not interfere each other):
As for
,
in Equation (25) contains its algebraic sign hidden in the term
. We will discuss the consequences and implications of Equations (24) and (25) in the
Section 5.
A direct evaluation of
ST through Equation (25) is not a simple task, of course, for this reason, dedicated numerical evaluations are in due course and their results will be reported in a separate paper [Peluso, F. et al., in preparation]. However, Equation (25) contains very useful and important information about the meaning of
ST in the DML framework. Recalling that
is the
liquid particle drift caused by the
collisions per second with the
lattice particle, Equation (25) becomes:
Here is the average temperature difference experienced by a solute liquid particle over the distance along covered after a collision with a phonon. This is a direct consequence of the dynamics characterizing the DML and fully in line with the macroscopic meaning of the Soret coefficient in NET.
We now have all the ingredients to calculate the flux of solute (or solvent) molecules
due to the collisions with phonons in the presence of an external temperature gradient. If all solute particles have mass
and their density is
, we have:
To recover from the present lack of getting a numerical evaluation of Equation (25), we propose here an analysis of thermodiffusion from another perspective, still in the DML framework, namely that of the variation of Π
wp, or Π
th, along the gradient. Besides being an energy density,
represents in fact also the pressure Π
wp exerted by elastic wave-packets on the
liquid particles. Compiling Equation (13) with Equations (20) and (21) we get:
Of course also for
the time-averaged pressure generated by the ubiquitous wave-packets streams in a system without external gradients is zero. This is no longer true in the presence of a temperature (or a concentration) gradient working on the system. Assuming as always that the external gradient is small enough that linear deviations suffice to describe the induced variations due to its application, Equation (22) still holds for
, i.e., the increase in the number of collisions due to the external gradient. Such an increase in the number of collisions causes in turn an increase in the energy and momentum exchanged upon the collisions, namely:
The reader has certainly guessed that
and
are respectively the total power dissipated by collisions and the force exerted on the collision target, i.e., the
liquid particle or the
lattice particle, depending on which situation prevails in the system, (a) or (b) of
Figure 3. At the Soret equilibrium,
is the power needed to sustain the thermal and concentration gradients of the mixture (the Soret equilibrium is characterized by a positive rate of entropy production, the ratio
at the cold side being higher than that at the hot side,
, while
is the net force responsible for the separation of the two chemical species, for instance, cubes and spheres of
Figure 1. Recalling Equations (10), (13), (14), (22), and (30), and adapting Equation (28) for the continuum to the “discrete” situation involving the two distinct
liquid particles “1” and “2”, we get that:
Equation (31) provides a comparison between two expressions for the total force, or the pressure, exerted on the targets following the application of an external temperature gradient. Comparing in fact its third and last members, we understand how the imbalance of the number of collisions influences the dynamical behavior of the mixture. In other words, Equation (31) tells us which of the two species, “1” or “2”, is pushed by the collisions to the cold side and why: the chemical species which is pushed towards the cold side by the external thermal gradient is the one for which the difference is positive, and therefore generates a virtual thermal gradient parallel to the external one, determining in turn the algebraic sign of in Equation (25) or Equation (26). Because the concentration of the solute, i.e., of the less concentrated species, is usually observed in solutions, one refers to positive or negative Soret depending on its concentration variation. Actually, when there is a negative Soret coefficient, according to the DML it happens that , or , is positive for the solvent molecules, which will therefore concentrate on the cold side, while the solute will concentrate on the hot side due to Newton’s third law. In general, saying that , or , is positive (negative) means in the DML that the radiation pressure due to the lattice particle stream on species “1” is higher (lesser) than on species “2”. Because the system is closed, it will obviously happen that only one of the two species will be pushed towards the cold side by the stream of lattice particles, while the other chemical species will move to the hot side due to Newton’s third law (a bit like cork floats on water while also being attracted by gravity, or a foam floating on the liquid surface will concentrate towards the axis of rotation in a closed rotating system).
In order to get an expression for
that can be easily evaluated on an experimental or numerical basis, we write Equation (31) at equilibrium replacing the expression for the thermal flux
using the Fourier law; we get then:
where
is the thermal conductivity. Proceeding now in the reasoning, what happens at the equilibrium, i.e., when also the concentration gradient is established in the mixture, is that the power generated by the thermal engine is dissipated against the viscous forces:
The velocity of the liquid particle is
defined by Equation (23); accordingly:
In Equation (35) we have used Equation (24) and have indicated with the suffix
ext the external temperature gradient. It is trivial to get:
Introducing the well-known Stokes–Einstein expression for the diffusion coefficient
eventually we get for
:
Equations (25) and (37) provide for the first time a simple interpretation of the Soret effect at the mesoscopic level, they are without matching parameters and can be verified by means of experiments. Equation (25) is based on elementary principles applied at the mesoscopic level, while Equation (37) is based on macroscopic quantities. Both equations and their implications, as well as their capabilities in the comparison with experimental data on the Soret coefficient and on its sign, will be extensively analyzed in the
Section 5.
4. An “Unexpected” Mechano-Thermal Effect in Pure Liquids Explained by Means of DML
The other cross-effect we will frame in the DML is that discovered by the group of Noirez a few years ago [
6,
7,
8,
9,
10,
11]. The effect results from the coupling between an acceleration field and a thermal field in a confined liquid layer submitted to shear strain. The system consists of a pure liquid layer confined in between two coaxial disc plates, one is fixed while the other one is mobile following a Heaviside function or an oscillatory motion, see
Figure 5 for a schematic. The system is observed by far by means of a thermocamera, which allows for measurement of the thermal field that is being established inside the liquid, as well its time evolution. When the Heaviside motion law is applied to the mobile plate, typically the lower one, it undergoes a sudden acceleration in one direction, performs the expected movement, and then stops for a few seconds. Then it goes back to the starting point with an acceleration lower than that of the outward path. In the second type of experiment, the rotating plate oscillates according to a periodic law. The thermal behavior of the liquid in the two cases is not exactly the same, although the dependence of the intensity of the mechano-thermal phenomenon, in particular of the temperature difference that originates with respect to the acceleration field, remains almost the same. From the classical point of view, even a hardly detectable (faint) heating of the liquid layer facing the moving plate is not expected. Indeed, at these low frequencies and strain amplitudes, the energy stored cannot generate viscous heating following the classical empirical evaluation of the Nahme number (also known as the Brinkman number). Much faster velocities are needed for viscous liquids, closer to sound velocities. Experiments performed by Noirez and co-workers revealed instead that unexpectedly liquids confined to sub-millimeter scale show a thermal response consisting of the cooling of the layer facing the moving plate and a corresponding heating of that facing the fixed plate. The aim of this section is to provide a theoretical interpretation of such counter-intuitive response in the framework of DML.
Generally, the non-uniform temperature shear bands observed along the strain direction indicate the formation of a faint temperature gradient in phase with the deformation. The synchronization of the thermal field with the mechanical strain is a signature of a mechano-thermal effect of elastic origin indicating that viscous liquids might exhibit thermoelastic behavior, as pointed out also by the authors of the experiments [
108,
109].
Noirez and co-workers pointed out that this response to low-frequency mechanical stress is typical of amorphous solids, exhibiting a storage shear modulus
much higher than the loss viscous modulus
. This is a key point in both of the experiments and the theoretical interpretation of the DML framework that we are going to show in what follows. In fact, being the geometry of these experiments of shear strain type, it is implicitly demonstrated by the experimental results that liquids possess shear elasticity as solids, contrary to what is generally supposed due to the absence of rigidity, but in line with the duality of liquids as supposed in the DML. By increasing the frequency of the mechanical stress and/or the scale on which such stress is applied, the effect disappears, and the liquid returns to its classical full-viscous behavior (
). Nonaffine framework has provided the physical interpretation of such behavior, consisting of the cut-off of low-frequency transverse modes, that in turn allows the high-frequency modes to maximize their effect at the mesoscopic level [
16]. In this framework, the authors have assumed that transverse modes are carried by phonons. Accordingly the DML may provide the missing ingredient to nonaffine treatment by means of the phonon
liquid particle interactions, where phonons assume the role of transforming the momentum they receive from the liquid particles in a thermal field. Therefore,
is the best candidate for the nonaffine force introduced in the equation of motion of the microscopic building block [
16]. In fact, in the nonaffine framework, there is a competition between the rigid-like behavior of coordinated clusters and the relaxations via nonaffine motions due to unbalanced forces in the disordered molecular environment.
Let’s then proceed with the DML interpretation. The key is provided indeed by the elementary
lattice particle liquid particle interaction, described in
Figure 3b. The acceleration of the moving plate generates a velocity gradient through the liquid, which in turn gives rise to a momentum flux, directed towards the fixed plate, following the relation
, with
the liquid viscosity. According to the DML, and also following Rayleigh’s reasoning, the momentum flux is dimensionally an energy density and generates a temperature gradient. The
liquid particles, pushed by the momentum flux, transfer their excess impulse to the phonons of the liquid, following the (b) event of
Figure 3; the phonons emerge from each single impact with an increased energy, being frequency-shifted. They are then pushed and collected towards the fixed plate, giving rise to a temperature gradient directed towards the fixed plate. The rationale is the same adopted to justify Equation (22), where however in this case the imbalance of collisions
affects the wave packets instead of the
liquid particles, and is due to the ratio between velocity gradients instead of temperature:
The supposed linear dependence this time is between the external velocity gradient
and that virtual, normally present in mechanical equilibrium for a
liquid particle between two successive collisions with the
lattice particles,
. Each wave packet executes
jumps per second in the direction of the momentum flux, each long
; they accumulate towards the fixed plate and decrease near the moving one, generating the temperature field observed. The product
defines the drift velocity of the thermal front due to the velocity gradient. The phenomenon is therefore very rapid because it is determined by a succession of (almost) elastic collisions, which propagate very quickly. Obviously, the higher thermal content of the liquid in proximity to the fixed plate compensates for the lower thermal content near the rotating plate, keeping the thermal balance unchanged (neglecting the losses toward the environment). It is important to highlight that the
lattice particle liquid particle collision also explains why a temperature gradient (i.e., a generalized flux-force crossed effect) is generated inside the liquid instead of a simple heating due to the energy stored as a consequence of the mechanical stress
, as one would have expected on the basis of a classical interpretation of the elastic modulus. At the same time, the phenomenon can be observed only on the small length scale as that adopted in the experiments, the faintness of the thermal gradient would be in fact destroyed by the thermal unrest of the molecules in larger liquid volumes.
When the plate stops, the system relaxes and reaches the temperature it had at the beginning. During the return of the rotating plate to the starting position, the liquid is in exactly the same condition it had at the beginning of the experiment, and the temperature gradient forms again. This interpretation is in line with all the basic hypotheses made by the authors, including the one on the number of Grüneisen. In this particular case, the expansion/compression effect of the material medium would be a consequence of the density of the phononic contribution.
We will now try and get an expression for the (rate of) temperature variation due to the occurrence of events of type (b) of
Figure 3 which are at the base of the observed phenomenon. In fact, what happens inside the liquid is that the mechanical energy injected into the liquid by the rotating plate is (partially) commuted in thermal energy carried by wave-packets. This energy, as we have seen, accumulates close to the fixed plate. To do this, we begin by considering that the momentum flux
through the liquid due to the plate rotation represents also the pressure
exerted by the liquid molecules on the wave-packets, and is given by [
110]:
where
and
are the density and shear speed of liquid, respectively. Neglecting dissipation effects due to the liquid viscosity, and remembering Equation (13), the momentum flux is converted into the pressure that pushes the phonon’s stream and gives rise to a heat flux along the opposite direction:
where
is the wave-packet average speed. As is well known, the momentum flux corresponds even to the energy density
injected into the system by the mechanical stress. The duality of the DML provides us a key to (partially) convert at a mesoscopic level such mechanical energy in thermal energy, i.e., that carried by phonons that are excited at higher frequencies upon interactions of type (b) of
Figure 3 with the
liquid particles. Neglecting variations with respect to the temperature, we may use the average macroscopic value given by the ratio
, where
is the height of the liquid sample (parallel to the rotation axis), and
the time necessary to establish the observed temperature difference across the liquid. Such value will be a multiple of the ratio
. Getting now the temperature difference
(or the temperature gradient
) is trivial, being enough to divide Equation (41) by the specific heat (per unit of volume) of the substance, i.e.,:
Using the data contained in [
11], we give now a numerical estimation of
from Equation (42). First, we use the value of the momentum flux averaged over the disc surface,
, where
is the plate angular velocity and
its radius. The liquid used in the experiment was glycerol, for which
and
. Using a value of
for
and
for
, one gets
of the correct order of magnitude as measured and reported in [
11].
It is relevant to point out that the above interpretation is made possible because the DML foresees the presence of wave-packets that, interacting with the liquid particles, behave in this case not only as carriers of the thermal energy but also as converters of mechanical energy into thermal energy.
5. Discussion
In this paper, we have used the DML to explain two cross effects in liquid mixtures and in a pure liquid out of equilibrium, namely the thermodiffusion and a mechano-thermal effect recently detected in isothermal liquids in shear strain geometry. Both have been interpreted thanks to the duality of liquids, supposed in the DML as constituted by metastable solid-like icebergs fluctuating in the liquid amorphous phase [
17]. Elastic perturbations, hence the heat, are carried by elastic wave-packets that interact with the
liquid particles exchanging with them energy and momentum. In particular, the exchange of momentum allows to provide an interpretation of the selective migration of solute and solvent molecules, which are pushed along the thermal current, i.e., the current of wave-packets, following the dynamics shown in
Figure 3a. More in detail, the separation of the components of a liquid mixture upon the application of a thermal gradient, is due to the capability of the forces acting at the mesoscopic level in a liquid to convert the heat flux into a momentum flux. An external temperature gradient has the effect of orienting the local gradients along the same direction. When instead
ST or
Dth are negative means that the local virtual gradients are oriented along the direction opposite to that of the external gradient (see Equations (25) and (37)). This circumstance produces a thermodiffusive drift along the direction opposite to that of the external gradient, resulting in a negative value of
. In other words, the thermodiffusive drift of the solvent prevails on that of the solute, consequently we will observe the solute migrating to the hot side and the solvent to the cold one because the force due to the radiation pressure exerted by wave-packets on the solvent particles is higher than that produced on the solute particles.
Through the equations describing the DML and the
liquid particle lattice particle interactions, we have provided for the first time an expression (actually two) for the Soret coefficient
without matching parameters, falsifiable (in the Popperian meaning) because verifiable through experiments and manifesting even the sign dependence typical of the Soret coefficient. The order of magnitude of
as well the sign inversion will be compared in this section with experimental data available in the literature. The time-reversal mechanism shown in
Figure 3b is instead responsible, in our vision, for the Dufour effect and of the mechano-thermal effect described in
Section 4.
In this section, we will discuss the equations and deepen the concepts introduced in the previous two paragraphs. Let’s then start with the expressions for
. Equation (26) allows us to speculate on a very intriguing aspect related to the meaning and definition of
in NET. As widely discussed in [
17,
19], in the phonon-
liquid particle collisions of
Figure 3a, part of the phonon energy is transformed into kinetic and potential energies of the
liquid particle. The interaction definitively has an activation threshold for potential energy and, depending on how much energy is absorbed by internal DoF, the rest of the transferred energy will increase the kinetic energy of the
liquid particle. When an external temperature gradient is applied to a system, the first effect that occurs is the establishment of the temperature gradient across the system, according to the Cattaneo–Fourier equation [
19]. Therefore, once the internal DoF have reached the statistical equilibrium in terms of the distribution of their degree of excitation depending on the local temperature, the phenomenon of matter diffusion takes place due to the temperature gradient and at the expense of the kinetic energy reservoir, Equation (21). Inelastic effects dominate the dynamics during the transient phase, leaving the control to the kinetic reservoir at the steady-state, where only the elastic effects of the elementary interactions matter. Accepting the above reasoning
may be seen as the activation energy for the thermodiffusion process:
where
is the solute concentration: if the energy
available at the temperature
exceeds the activation energy
, or alternatively, the average temperature of the system exceeds the temperature difference over
following the elementary interaction (see Equations (25) and (26)), the phenomenon takes place. It should be interesting to theoretically speculate on whether, and how, such activation energy is related to the presence in liquids of
gapped momentum states [
111,
112]. We have already discussed in [
17] how the “
lattice particle liquid particle” interaction can be assumed as responsible for such phenomenon in liquids exhibiting
k-gap.
Of the two expressions obtained for , Equation (25) is based on elementary principles applied at the mesoscopic level, while Equation (37) is based on macroscopic quantities. A question that the reader has certainly wondered is what one expects from their comparison. Let’s then consider the third member of Equation (37); the numerator is the force exerted by the thermal field to displace the solute or solvent, molecules to generate and sustain the concentration gradient that gives origin to the Soret effect. At equilibrium the (average) energy available in the system is . The ratio between such two quantities is just the mean displacement of the solute (solvent) particle, , due to the collisions with the phonons, . Consequently, Equations (25) and (37) formally represent two similar expressions for the Soret coefficient (apart from a numerical factor due to statistical averages).
Equation (37) allows us to open a very interesting discussion, based on several points. The first and more immediate is that the last member is obtained by simplifying the temperature gradient in the two parentheses in the numerator with that in the denominator. Equation (37) has a huge relevance because it allows it to evaluate
more easily than Equation (25). Nevertheless, the difficulty in this case arises from the fact that the physical quantities involved are not, in principle, those of the bulk materials, but rather those in the effective configuration in which they are in the mixture. What helps is that in the DML, solvent, and solute
particles are not seen as isolated molecules but rather as islands, those we have dubbed
icebergs or
liquid particles. Incidentally, the proportionality of Equation (37) to the cross-section of the solute particle copes with the law
∝
to the solute mass early evidenced by Emery and Drickamer [
46] and successively by other authors [
40,
52] for solutes of large molecular weight. For such large particles, in fact, the surrounding layer of solvent molecules is relatively less relevant than for small solute particles.
In
Table 1 some examples of Soret coefficients evaluated from Equation (37) are provided for both positive and negative
.
Table 2 collects data for the rations
elaborated from thermodiffusion in liquid mixtures performed with a Klusius–Dickel device. Original data report only which of the two species involved in each single experiment concentrate on the hot or cold side, however, this is enough to verify how the direction of the thermodiffusion drift is correctly predicted from Equation (37).
Figure 6 shows another way to conceptually validate the expression of
given by Equation (37).
Figure 6a collects the data obtained by Bierlein, Finch, and Bowers for several mixtures of benzene and n-heptane [
116]. The interesting point is that
for such liquid mixtures increases with the average temperature for all the relative concentration values of the two components. Even more interesting is the fact that the three plots obtained as a function of temperature cross the “zero” value for
at about the same temperature, around 60 °C. More precisely,
is negative for temperature below such threshold, goes through zero, and becomes positive for average temperature above the threshold. In
Figure 6b we have plotted the values of the ratio
for the two components; it is straightforward to note that such a ratio assumes the same value around the temperature of 60 °C. Because the sign of
in Equation (37) depends just upon the difference of the ratios
for the two components, this observation confirms that Equation (37) is a good candidate to foresee even the sign inversion of
. As specified above and also in
Section 3, Equation (25), that is in principle the correct expression for
in the DML, contains the sign dependence hidden into the term
, more precisely in the orientation of
with respect to the external temperature gradient. Supporting the above reasoning is the dependence of
in Equation (25) upon
instead of the external gradient. This is not surprising because the thermodiffusion is a signature of the mixture, therefore it shall be dependent only upon the mixture’s parameters. Moreover, the reader must not be fooled by the inverse proportionality with
: it tells us that as large the local virtual gradient
, as difficult will be to observe the thermodiffusion, as one would have expected. Of course,
as well
or
will change according to the relative concentration of the components, or to the temperature of the system.
A more extensive collection of experimental data on thermodiffusion is in due course and their comparison with the theoretical predictions of Equations (25) and (37) will be shown elsewhere. Because we do not have access to a lab, the only possibility to provide further quantitative evaluations of Equation (37) is through experimental data available in the literature for which enough data should even be available for K’s and uϕ’s. As it concerns the validation of Equation (25), dedicated numerical analysis will be the core of another work.
The comparison between Equations (25) and (37) opens to a more wide reasoning concerning the displacement of the energy and mass barycenters of a mixture as a consequence of the thermodiffusion. Making reference here to
Figure 7, the drawing on the left shows a mixture, “0”, constituted by two chemical species, “1” and “2”, in a container in thermal equilibrium at
with the environment; only the exchange of thermal energy is allowed with the environment. For such a system the mass,
,
,
, and energy,
E, barycenters are all coincident and positioned in the cell middle.
Figure 7b on the right shows the same system once the Soret equilibrium is reached. In this situation, we have a heat flux
crossing the system, with
, so that the four barycenters are no longer coincident. In the drawing we have supposed “1” to be the solute, with a molecular mass
and ratio
higher than those of the solvent “2”,
and
; this explains the positioning of the three mass barycenters. As it concerns the energy barycenter
, it is shifted towards the warmer side due to the higher thermal content of the system. Several considerations are now in order. (i) The force acting on the mixture’s components pushes the solute particles toward the cold side; consequently, the solvent particles are pushed toward the warm side. The fact that
is shifted toward the opposite side with respect to
, is a direct consequence of the III Principle of dynamics, being the system closed, so that a net volume flux is not allowed,
. This is a relevant point because some authors in the past have argued that the phenomenon of thermodiffusion was characterized by a violation of the III Principle [
47]. (ii) The energy introduced into the system at the expense of the external thermal field ensures the stability of the mass and energy distributions obtained as effects of the thermodiffusion and gives also rise to a positive rate of entropy production. However, to keep such distributions stable with time, it is necessary that part of this energy is converted into momentum, which pushes the solute and solvent molecules towards their respective barycenters (see Equations (29) and (30) and related comments). Let’s then calculate this momentum rate, for which we will see that the duality of liquids of the DML is a necessary prerequisite, being the phonons the carriers of such momentum flux. If Δ
t is the time interval to reach the Soret equilibrium, the net mass current
will be given by (
Figure 7b):
where
is the cell cross area and
the coefficient of thermal diffusion of species “1”. The mechanical momentum flux,
at equilibrium, is then given by:
Equation (45) gives the net rate of transfer of momentum from the phonons to the solute particles in the course of the processes of type (a) shown in
Figure 3.
It is instructive and interesting to perform the same calculation but in a situation in which the system is closed and thermally isolated, and the starting situation is the one in which the two components are already separated, for instance with a solute concentrated solution on the bottom of the cell and the pure solvent on the top. We all know that the final distribution reached by the system will be that of a homogeneous mixture; we could thus think of a time reversal experiment with respect to thermal diffusion, known in NET as the ordinary diffusion giving rise to the Dufour effect. Of course, the constraint
= 0 is valid also in this case. Equation (44) now becomes:
is the density of the solution,
the solute diffusion coefficient and
the time averaged value of the solute concentration gradient. The momentum transfer accordingly is:
This time Equation (47) represents the momentum transfer from the
liquid particles to the phonons, following the events represented in
Figure 3b. The reader may now wonder how it is possible that the momentum of an isolated system has changed. The duality of the DML makes the answer straightforward; in fact, the phonon population was initially uniformly distributed in the system, which was indeed in thermal equilibrium, but not the population of solute particles, for which a concentration gradient is present at the beginning of the experiment. The mass flow due to the concentration gradient produced in turn a phonon concentration gradient opposite to the solute concentration gradient, due to an excess of collisions with solute particles in the direction of the mass flow, i.e., following just the type (b) processes of
Figure 3. This quantity is used by the
liquid particles to shift their barycenter mass. Of course, and as is well known, the Dufour effect produces a (faint) thermal gradient, that we identify as that of the phonon concentration gradient that is established in the system as a consequence of the
liquid particle current. Once again, the two events shown in
Figure 3 represent the elementary interactions occurring in liquids, either in equilibrium or out of equilibrium, their time-reversibility being the best candidate to justify the Onsager reciprocity law at the mesoscopic level [
17,
19].
The mechano-thermal effect analyzed in
Section 4 even deserves some comments. Equation (42), with all the limitations and simplifications adopted, shows that the temperature difference is proportional to the square of the shear speed, making the Δ
T orientation independent upon the shear versus, as expected. The energy and momentum input depend only upon the momentum of the moving disc, and the Δ
T observed is independent of the amount of substance or of the gap thickness of the layer (
is the specific heat per unit of volume). In the comparison of our results with the experimental observation, we should on one side keep in mind either the simplifications adopted in obtaining Equation (42), and the experimental difficulties in getting a stable and well-thermally isolated experimental setup. Nevertheless, on the other side, we could speculate on those aspects that may appear contradictory. One of these is that in some cases an apparent inversion of the thermal gradient was observed upon inversion of the shear motion. This effect could have the same origin as that observed in those cases where more than a single gradient was observed in the liquid, a sort of alternation between hot and cold layers. A possible explanation, although on a purely speculative basis, could be that the stationary plate reacts to the inversion of the shear direction with a sort of “mirror” effect, pushing back the thermal excitations. The liquid near the stationary plate perceives the shear inversion as a sort of mirror effect of the stationary plate. Very interestingly, upon the revision process, one of the Reviewers called my attention to the non-Fourier processes recently identified in various thermal mechanisms. The review by G. Chen [
117] may help us to understand the origin of such effects. Generally, the momentum of normal phonons is conserved in scattering processes; in systems where normal phonon scattering dominates, phonons under temperature gradient can acquire a non-vanishing drift velocity similar to fluid flow driven by a pressure gradient: this is named the phonon hydrodynamic regime, which originates from the Landau treatment of superfluid Helium [
103], and has been already dealt with in detail in a previous paper on the DML [
17]. The Casimir–Knudsen is one among the several heat conduction phenomena associated with the phonon hydrodynamic regime. A necessary (but not sufficient) condition for the Casimir–Knudsen regime to take place is that the phonon mean free path
is comparable to, or larger than, the characteristic length of the system, say
. In the system under study, the phonons propagating in the pure liquid interact with the material constituting the plates. Hence the problem is restricted to considering such interface, of thickness
, and the related highly uncertain boundary conditions. These uncertainties arise because of the atomic-scale defects present at the interface affected by physical roughness, in turn responsible for how phonons are scattered because the phonon wavelength amounts to several nanometers, as estimated also in DML ([
17],
Table 1). To bypass the problem a “specularity” parameter is usually introduced, accounting for the fraction of “specularly” reflected phonons, which should depend on the phonon wavelength and surface roughness [
118]. The speculative explanation proposed above could then be a consequence of the shear strain and may represent coherence bands of phonons originated just at the interface between the plate and the liquid. The incoming phonons experience reflection and transmission, being responsible for the temperature drop across the interface. This phenomenon, quantified in
Section 4, could explain either the inversion of the temperature drop upon inversion of the plate rotation or the formation of multiple temperature layers.
One of the characteristics of the Casimir–Knudsen regime is that the thermal conductivity is linearly proportional to the specific heat, as is the case just for the thermal conductivity
in the DML (see Equation (12) in reference [
19]). During the thermal transient, the heat pulse generates heat wave propagation (second sound) similar to a mechanical pulse generating a sound wave in a solid. Thus, as extensively described in [
19], phonon heat propagation can be described using a Cattaneo-type (telegraph) equation [
117]. In addition to thermal transport, it has also been predicted that signatures of phonon hydrodynamic transport can be observed from neutron and Brillouin scattering experiments [
17,
119]. In conclusion, the DML allows a direct comparison with many other interesting aspects characterizing either the solid-like behavior of liquids or the presence of phonons as carriers of elastic (and thermal) excitations.
Returning to Equation (28) or Equation (31), it is interesting to note that
represents the radiation pressure exerted by the current of wave-packets on the
liquid particles following the interaction. This pressure may also be regarded as a sort of osmotic pressure determining the displacement of
liquid particles following the collisions with the wave-packets. Indeed, to a pressure difference,
the pure liquid responds with the self-diffusion of the
liquid particles [
17,
19]. The tunnel effect could be regarded as a semi-permeable membrane, that works allowing only the passage of
liquid particles and preventing that of
lattice particles. A similar argument was raised by Ward and Wilks [
120] in dealing with the mechano-thermal effect observed in HeII.
Finally, some considerations could be dedicated to a comparison between the capabilities of DML applied to cross effects in NET and the previous theoretical models for thermodiffusion. As anticipated in
Section 2, we do not want here to enter into the debates of the models or make a historical excursus, that is however the topic of other interesting papers [
12,
13], but only propose a critical comparison of the DML with some of the most relevant alternative models. The first comment is dedicated to the interesting paper by Najafi and Golestanian [
47], in which they suppose the thermodiffusion to be due to the Soret-Casimir effect. They in fact calculate such force in several configurations, however this method is intrinsically unable to foresee negative values for
. Besides, they have also argued that in some cases the Soret effect violates the III Principle of dynamics, a circumstance that has been definitively excluded in the present manuscript because of the liquid’s duality. Duhr and Braun [
52] have modeled the thermodiffusion due to the gradients of thermodynamic potentials inside mixtures. They have applied the model to macromolecular solutions, in particular DNA, finding good correspondence also with the mass and size of the solute, however only in those cases with positive
.
Two very interesting historical reviews are represented by the works of Eslamian and Saghir [
12] and Rahman and Saghir [
13], from which a large amount of models of thermodiffusion jumps to the eye, very different from one another and often from experimental data. These reviews are united by similar conclusions on the several models they deal with. In particular, it is affirmed that many theoretical models, initially accepted by the scientific community, were found wrong after more careful analysis [
49,
121]; many others correctly reproduce the thermodiffusion only in gases, while “in liquids, reliable and satisfactory predictive theories are still lacking”; furthermore “the largest part of models necessitate of external matching parameters to be defined, making the models dependent on the experiments” [
12]. A relevant point is the fact that quite all the models either do not predict negative
or fail in its prediction [
14].
All the above points of weakness are overcome by the DML, as discussed above.