1. Introduction
To understand transport through, into, and along interfaces is of practical importance in physics, chemistry biology, and engineering [
1,
2,
3,
4,
5,
6,
7,
8]. Large efforts have been made to develop a rigorous and thermodynamically consistent description of surfaces at global equilibrium and away from this condition. Gibbs laid the foundation in his description of surfaces at equilibrium, using excess variables for energy and mass [
9]. The method of dynamic boundary conditions [
10,
11] was proposed in 1976 to deal with surfaces away from equilibrium. The method provides a description of interface transport phenomena that are compatible with the second law of thermodynamics. Current methods for interface transport often do not possess the symmetry required by this method (Onsager symmetry) [
3]. In spite of the agreement with symmetry requirements and the second law, the dynamic boundary method has not yet gained common use. One reason may be that there is still only scattered knowledge on the validity of its fundamental assumptions. This article aims to mend this situation by presenting rigorous evidence for the validity and applicability of the basic assumptions of the theory of non-equilibrium thermodynamics for interfaces.
As in the original method [
10,
11] and later work [
3], we shall use Gibbs definition of excess thermodynamic properties [
9] to describe the interface. It is well known that these properties depend on a mathematical construct, namely, the dividing surface. The dividing surface has a location that needs to be defined. But can a thermodynamic description depend on a mathematical construct, one may therefore ask? The intensity of reflected light, for instance, certainly does not depend on this location [
12], and neither does the surface tension. We shall see that this is also the case for other properly chosen surface excess variables.
In the dynamic boundary method [
10,
11], the crucial step was to proclaim the surface between two bulk phases to be an autonomous system, also away from global equilibrium. Gibbs was working with interfaces between phases in global equilibrium. In global equilibrium, the temperature and the chemical potentials are everywhere the same and well defined. But can we define such variables in a unique way, when the system is away from global equilibrium and the variables are everywhere different? By proclaiming the surface with its excess variables to be an autonomous system, we imply that the system is in local equilibrium; or that thermodynamic equations apply. The assumption is central in non-equilibrium thermodynamics. The theory builds on the assumption of local equilibrium and of independence of mathematical constructs. Evidence for both properties shall be presented and examined for an interface of multi-component fluids in this article.
Away from global equilibrium, the temperature and chemical potentials of the surface will in general differ from the values of the adjacent bulk phases. A measured surface tension may therefore disagree with the value that can be predicted from the variables of the adjacent homogeneous phase. This problem has lead to the introduction of the concept of a dynamic surface tension [
13]. Such a variable may turn out not to be necessary, with more information of the variables of the autonomous surface [
14].
The present paper addresses these and related questions. We shall see that a unique definition of surface variables exists, independent of the chosen location of the dividing surface. Furthermore, we shall propose that the property of local equilibrium can be actively used to find surface properties. From knowledge of proper thermodynamic relations for one choice of the location of the dividing surface, one can find the relations for another choice, through transformation of excess densities. This can also help clarify misunderstandings in the literature.
We aim to clarify definitions used in the thermodynamic description of a fluid-fluid interface of mixtures. We restrict ourselves to flat interfaces, but the analysis can be extended to curved ones [
10,
11]. We consider also only cases where gradients and fluxes are in the direction normal to the surface.
Savin et al. [
7] and Öttinger et al. [
9] addressed the same problems using the concept of gauge invariance. A gauge transformation is a displacement of the dividing surface over an atomistic distance. The temperature and chemical potentials of the surface were found to be gauge-invariant, whereas the excess densities were gauge-variant. Schweizer et al. [
8] verified the results for a one-component liquid-vapor interface, using molecular dynamics simulations. The conclusions of these papers support those discussed here. Nevertheless, we consider it useful to further strengthen the conclusions, following a somewhat different path, as outlined below.
The dividing surface can essentially be displaced in two ways. One can invoke a notional change only. This means that a different dividing surface is chosen without physically moving the original surface. The other way is to physically move the surface to a new location. In the present paper, we are only concerned with notional changes. A notional change is the same as a gauge transformation. This means that the results obtained earlier will be obtained again, but in a different way.
We start by defining in
Section 2 what we mean by the dividing surface and an excess density, and show how the concept of excess densities depend on the surface location. Molecular dynamic simulations or density functional theories can be used to obtain continuous profiles of energy- and mass-densities across interfaces. In the interfacial region these profiles change considerably over the thickness of the surface. Gibbs used the profiles to derive his concept of a surface variable.
Section 3 presents the Gibbs excess thermodynamic properties of a surface. We shall use the condition of global equilibrium to find an equivalent variable set, a set that is gauge-invariant or does not depend on the location of the surface.
From the set of gauge-invariant excess densities in
Section 4, we proceed in
Section 5 to show that a corresponding set of surface intensive variables is also gauge-invariant. This knowledge in combination with the assumption of local equilibrium can be actively used to find unknown properties of a surface. We discuss how we can measure thermodynamic quantities of the surface away from global equilibrium. A measurement of, for example, the surface tension can then be regarded as equivalent to a temperature measurement.
The property of local equilibrium has a remarkable consequence; away from global equilibrium, the equation of state for the surface, once the other variables are known, can be used to find dependent variables. Experiments to verify such relations will be proposed.
2. Excess Densities Defined
We recapitulate Gibbs’ definition of excess variables and identify the problem with location-dependent thermodynamic properties. The excess density of component
j in moles per unit of surface area is
where the integral is over the length of the system,
is the Heaviside function (1 for a positive and 0 for a negative argument),
x is the coordinate measured normal to the surface, and
is the chosen location of the dividing surface. On a molecular scale, the density
changes continuously throughout the interfacial region. The values of
and
are the extrapolated densities of the component
j in the two fluid phases. There are
n components. We use the symbol
to denote the excess number of moles per surface area [
15,
16,
17]. The symbol
has been used for the same property.
The variation in the densities of the bulk phases over the thickness of the interfacial region is negligible. All densities depend on the time
t, but this will not be explicitly indicated. The excess energy and entropy densities are also defined per unit of surface area. The value of
depends on the choice of the location,
, of the dividing surface. A common and convenient choice of this location is the equimolar surface
of component
j. This is defined by
The location of an equimolar surface refers to one component. When the equimolar surface of one component is used as a dividing surface, the excess densities of all the other components are not equal to zero. Excess molar densities can be both positive and negative.
The surface tension (in Pa·m = J·m
) is defined by minus the excess parallel pressure and is equal to [
5,
9]
The excess internal energy density (in J·m
) is
The excess Helmholtz energy,
, the excess entropy,
, and the excess enthalpy,
, of the surface are given by similar definitions. The excess Gibbs energy is
where
are the chemical potentials.
The extrapolated parallel pressures
and
of a flat surface away from the surface are equal to the extrapolated normal pressures
and
. The normal pressures are the same and constant in the interfacial region,
. This is the reason why the surface tension of a flat surface is independent of the location of the dividing surface [
5,
9]. All excess densities, except the surface tension, depend on the location
of the dividing surface.
The excess entropy production is standardly found from the Gibbs equation and the balance equations for the excess mass, energy, and entropy densities [
3]. The entropy production of the surface is an absolute quantity. It cannot depend on any reference value. The fact that the excess properties depend on the location of the dividing surface is therefore not convenient.
We need to deal with this problem in order to understand the role of the dividing surface, and how it disappears from the entropy production. This will be done in two steps. We introduce a special reference for the chemical potentials in the first step (
Section 3). We suppose in the second step that the validity of local equilibrium implies that the wanted invariance is also true. The needed relations can be tested experimentally (
Section 4).
3. Thermodynamic Properties of Gibbs Surface
3.1. The Interface at Global Equilibrium with the Surroundings
The purpose of this subsection is to show that a set of geometry-dependent, excess densities in global equilibrium can be replaced by a set of absolute intensive variables, provided that we introduce a convenient reference chemical potential for the components.
Take the excess functions of the liquid-vapor interface as an example. In order to denote
global equilibrium, we introduce the subscript
. For an arbitrary choice of the dividing surface, we have from Gibbs [
9]
Superscript s is used to indicate excess densities of the surface that depend on the location
of the dividing surface. The excess quantities depend on the temperature and all but one of the chemical potentials. In equilibrium, the temperature and the chemical potentials are constant throughout the system and are therefore accessible experimentally via the bulk phase values. The definitions of Equations (
3)–(
5) apply.
The thermodynamic excess densities
,
,
, and
as well as
need a reference value to be uniquely defined. In equilibrium, one chemical potential depends on the temperature and the other chemical potentials. We then renumber the components so that this component has number
n. We now use
to define absolute excess densities:
where
becomes absolute because identical standard chemical potentials can be chosen, and the difference can be found by experiments. Equation (
6) can now be alternatively written as
These redefined excess energy densities of the surface depend only on . They are all absolute excess densities.
In global equilibrium, the set is the same throughout the bulk phases and the interface. The values of the homogeneous phases can be found experimentally and applied to the interface. We can then tabulate the set of densities for the homogeneous phases, the set of excess densities for the surface, and as a function of . The tables of the excess densities can be used as “equations of state” for the interface. These lists will be further discussed below.
3.2. The Interface Away from Global Equilibrium
We explain here the meaning of local equilibrium for the interface, and find under these circumstances that there is a relation between the densities in and at the interface.
A fluid-fluid interface can be studied also away from global equilibrium. Experimentally we can, for instance, measure the surface tension. From molecular dynamics simulations or with density functional methods, we can obtain density- and other profiles through the surface when it is exposed to gradients. Such profiles have been written as [
15,
16,
17]:
where
. (We have not explicitly indicated a possible time dependence.) The densities
as well as the chemical potentials
need a reference value. We use the same reference as above and define the absolute densities
where
and
, which are also absolute. Note that
. Following the procedure outlined above, we take the excesses of these expressions and obtain
where
and
. Furthermore,
and
are defined similarly to the other excess densities. All excess densities with a hat are now absolute.
Local equilibrium means that normal thermodynamic equations apply. For a surface in local equilibrium, all surface excess densities have the equilibrium value given by temperature and relative chemical potentials . For a surface in local equilibrium, we can also convert the above expressions into Gibbs’ relations for the surface, also in the presence of external fields.
In equilibrium, we can identify the temperature and the chemical potentials of the surface:
By furthermore substituting Equations (
12) and (
13) into Equation (
11), we obtain
For the absolute excess densities, the equivalent equations are
These are the relations given originally for autonomous interfaces by Gibbs, cf. Equation (
8). While Gibbs gave these relations only for a fluid-fluid interface in global equilibrium, we have now shown, through the use of the identifications of Equation (
12), that they are also valid away from equilibrium.
The importance of this cannot be understated. It means that the theory of non-equilibrium thermodynamics can be used to describe transport along, into, and through the surface (see [
3]). From the Gibbs relation and the balance equations for energy, entropy, and masses, we can obtain an expression for the entropy production away from global equilibrium. Relations can also be formulated for surface excess density fluctuations [
11]. The entropy production can be used to identify relevant thermodynamic forces and fluxes, and their linear relations leads to a complete description of the time dependence of the properties of the surface.
Two important questions remain, however. How can we determine the temperature and the relative chemical potentials of a surface in a non-equilibrium system? Moreover, why are , independent of the choice of the location of the dividing surface (gauge-invariant)?
4. Excess Densities and the Choice of Dividing Surface
To answer the questions formulated in the last section, consider for the sake of argument the two most common choices for the dividing surface: the equimolar surfaces and the surface of tension. The position,
, of the equimolar surface for component
j is defined by
where the integral is over the whole box, cf. Equations (
1) and (
2). The position,
, of the surface of tension, on the other hand, is defined by setting the first moment of the excess parallel pressure equal to zero
Other locations are possible. One can, for instance, set the excess internal energy density equal to zero and use the equation to determine the position of the corresponding dividing surface.
For the surface of tension, all excess molar densities are not equal to zero. For the equimolar surfaces, the first moment of the excess parallel pressure is not equal to zero. With knowledge of excess densities for a particular choice of the dividing surface, one can find the set of corresponding excess densities for other choices of the dividing surface. Consider the two choices of the location of the dividing surface
and
. We take
. Using the definition in Equation (
3), we find
The variation of the extrapolated values of the homogeneous phases,
and
over the thickness of the surface, is negligible. This implies that
. For this reason,
can be taken to be independent of the location of the dividing surface, and Equation (
18) gives
Similar expressions are valid for the other excess densities. This then results in
All relations in Equation (
20a) were found to apply in very good approximation for a one-component fluid [
18] and for a two-component fluid mixture [
16].
Following Gibbs, we introduce the excess concentrations of component
j relative to component
n:
It follows from Equation (
20a) that
The relative excess concentrations of a flat surface therefore do not depend on the location of the dividing surface (they are gauge-invariant). This is precisely the reason why Gibbs [
9] introduced the relative excess concentrations in global equilibrium.
We can similarly define all other excess densities relative to component
n by
All the excess densities relative to component
n are gauge-invariant or independent of the location of the dividing surface:
We have shown that the excess densities relative to component
n do not depend on the location of the dividing surface. This is true in global as well as local equilibrium. Non-equilibrium molecular dynamic (NEMD) simulations [
18,
19,
20] and the square gradient analysis [
18] of the liquid-vapor interface for a one-component fluid have confirmed this property. For a two-component fluid mixture, this was also documented [
16] using the square gradient analysis.