1. Introduction
In the thermodynamic theory of small systems, constructed by Terrel L. Hill in the 1960’ies [
1,
2,
3], the subdivision potential
is central. In spite of its central position in a Gibbs type of equation, its meaning has not been made clear over the years. It still is referred to as confusing, ill defined, unmeasurable and even unnecessary [
4]. It has been claimed to be zero in equilibrium under certain cases. While many aspects of the subdivision potential are obvious, it is also difficult to get used to. The equation of state, for instance, is different for different sets of variables. The shape of the system can be a control variable. This is contrary to classical thermodynamics. The subdivision potential depends in other words on the ensemble in use. For some cases one needs two pressures rather than one, or two chemical potentials rather than one [
4]. In addition to the common pressure (chemical potential) the so-called integral pressure (chemical potential) has been defined.
The aim of this paper is to define, clarify and defend this property as it is used to describe various ensembles in equilibrium. We shall clarify how the subdivision potential enters can be accessed by measurements. As most other thermodynamic properties it does not have an absolute value; only differences can be obtained. In this sense the subdivision potential is not directly measurable, but can be computed from measured values. We shall use the word measurable in this sense. We shall also present statistical mechanical models for the subdivision potential, accessible by molecular simulations. As example we consider an ideal gas (mixture) of a sphere with surface adsorption, cf.
Figure 1. The adsorption energy of each gas is taken to be independent of the control variables used of the ensemble of replicas of the small system. We shall give analytic expression for some ensembles explicitly.
We shall see how changes when the temperature, pressure and chemical potentials are changed in six relevant ensembles. These changes are all measurable.
Hill introduced the subdivision potential
, as the energy needed to subdivide a larger system over one more replica. It is a characteristic of the dependence of the energy on the size of the small system. As stated above
depends on the choice of environmental and other variables. To add variables like the surface area and its curvature [
5] to Gibbs equation may reduce the need for
.
The systematic approach to include size and shape as variables in a thermodynamic description, is to first choose the set of control variables. These variables for a similar, but large system are often the same. The important difference is that the variables for the small system are no longer Euler homogeneous functions of the first order in the size. The next step is to consider an ensemble of replicas of the small system. The total properties of the ensemble are now Euler homogeneous functions of the first order in the number of replicas. This in a sense restores the validity of the Euler homogeneity and makes the thermodynamic treatment of small systems possible, as we will discuss in detail below.
This work builds heavily on a recently published monograph on the subject [
6]. The general basis of the formulas for the subdivision potential are given in this book, but are here extended and strengthened. We provide explicit formulas from the thermodynamic and the statistical mechanical points of view. We start by recapitulating the general methodology of Hill’s work. In the subsequent sections we treat the specifics of the small system for the nanocanonical, microcanonical, canonical, isobaric and grand canonical ensembles. We end with a discussion of the results and a conclusion.
2. Systems That Are No longer Euler Homogeneous in V
When a system is small, properties like the internal energy,
U, entropy,
S, and mole numbers of the components,
, are no longer Euler homogeneous functions of the first order in the volume,
V, of the system. In short we may say that they are no longer extensive in
V. Contributions to thermodynamic variables due to non-negligible surface energy, curvature (shape) and three-phase contact lines, are the origin of this behavior. Hill [
1,
2] resolved the size- and shape dependence by considering an ensemble of replicas of the small system. The number of replicas,
, is then very large. For the total internal energy of all replicas,
, he wrote
Subscript
t indicates the sum of the corresponding single system properties over all the replicas of the ensemble. The ensemble of replica properties becomes Euler homogeneous functions of the replica number, by this choice.
Compared to the standard Gibbs equation there is a new term,
, which gives the internal energy needed to change the number of replicas, while
are kept constant
In other words
is the internal energy needed to subdivide the same
over one more replica. Because of the relation to this process, Hill suggested the name
subdivision potential, a suggestion we follow. Equation
1 was therefore called the Hill-Gibbs equation by us. The ensemble totals
are Euler homogeneous functions of the first order in the number of replicas,
. It follows that
The subdivision potential is, as we shall see, a measure of how the energy of the replicas depends on their size and shape. It several ways it is similar to the chemical potential, as it describes the deviation of properties from those of bulk materials.
The ensemble averages are given by
Together with Eq.
3 this implies that
From this relation we see that the subdivision potential
is a measure of how much
deviate from being Euler homogeneous functions of the first order in the volume,
V, of the small system. We substitute
times the averages into Eq.
1, use Eq.
5, and obtain
The averages satisfy the usual Gibbs equation. The difference to the description with the standard large system case, is that the small system is not Euler homogeneous in
V. This implies that for small systems
. The finite subdivision potential
is a consequence of this.
Related to
is the socalled hat-variables. A hat-variable describes the effect of a finite volume changes, introduced by adding or subtracting one replica to the ensemble. For the nanocanonical ensemble, there is such an effect on the internal energy called the integral internal energy
or
U-hat. It was found that
is equal to
[
6]. Similar hat variables will be introduced for more ensembles belowe. For small systems, the variables of the ensemble of replicas (variables with subscript
t) are Euler homogeneous functions of degree one in the replica number.
For small systems therefore, the thermodynamic properties depend on the variables that are controlled, and the shape and size of the volume. This is unlike the situation in standard thermodynamics. With small system thermodynamics, each ensemble has a different equation of state.
These differences will now be illustrated the particular subdivision potential for a multi-component ideal gas in a sphere with surface adsorption. We examine this system in six ensemble settings.
3. Nanocanonical Ensemble
In the nanocanonical ensemble the temperature, T, pressure, chemical potentials, , and the shape of the volume are controlled by putting the system in a reservoir. The system wall is open for exchange of particles with chemical potential . This particular choice is only appropriate for small systems. For large systems the pressure is a unique function of T and . We add a superscript NC to indicate the variable choice.
For small systems it is of special interest to introduce the so-called integral internal energy
which is also called
U-hat. It follows from Eq.
5 that
For large systems the subdivision potential is zero so that
.
By differentiating Eq.
5 and using Eq.
6 we obtain the Hill-Gibbs-Duhem equation
It follows that
These three equations can in principle all be used to determine the subdivision potential. The top equation involves the entropy which is not too easy to find. The center equation expresses that we can obtain
by integrating the volume
as a function of pressure at constant
. From the bottom equation, we can find
by integrating
as a function of
at constant
and
. The central equation is analogous to the standard way to obtain the chemical potential. It provides possibly the best experimental route to the subdivision potential out of the three. With reference to some standard value, we can obtain the subdivision potential for a new state by the integration:
Similar equations are valid for the dependence on the temperature and the chemical potentials. An important implication is that the subdivision potential is measurable using these integral relations.
A statistical mechanical expression for the variable is of interest. For the example of an ideal multi-component mixture in a cube with surface adsorption, we calculated earlier [
6] the subdivision potential and found:
The equation was derived, with the assumption that the adsorption energy
is independent of
, see [
6] for details
1. The expression gives the ideal gas law as a leading term. Deviations from this law are a function of the chemical potential of all particles,
, the particle’s adsorption energies,
, as well as of
, the Broeglie wavelength. The factor
is the gamma-function.
For the sphere in
Figure 1 we find a similar expression:
The results for the cube and the sphere are rather similar. For the cube, we obtain the factor 6
, while the corresponding prefactor for the sphere is
. The results reflects the fact that the surface of a sphere is smaller than of a cube with the same volume. The two formulae lend themselves available for tests by molecular dynamics simulations. Through such simulations one may be able to verify the prefactor of the size-dependent term, and obtain theoretical support for the theory [
6].
4. Microcanonical Ensemble
In the microcanonical ensemble one controls the internal energy,
U, the volume,
V, and the mole numbers,
, all per replica. Their total values in the ensemble of replicas are therefore,
,
and
. Substitution in Eq.
1 gives
where the superscript MC indicates that the corresponding variable depends on
and the shape of the system. The so-called
replica energy (the prefactor of
) is therefore given by
It follows from Eq.
14 that
Using that
is an Euler homogeneous function of the first order of
it follows from Eq.
14 that
By substitution of Eqs.
15 and
17 into Eq.
14 we obtain the Gibbs equation for small systems
It has exactly the same form as for large systems. The important distinction is that for small systems
are not Euler homogeneous functions of
V. They have this property for large systems and in that case
, which implies that
. We define the integral entropy for small systems by
Together with Eqs.
15 and
17 it follows that
This property is typical for small systems.
It follows from Eq.
18 that
Using Eqs.
18 and
19 one obtains
This implies that
In Eqs.
22 and
23, the proper variables to use are the Planck potentials
. In order to write changes of
as a function of these Planck potentials, similar to Eq.
11, one should introduce new symbols for these potentials to avoid confusion. We will not do this explicitly.
Equation
23 prescribes how the subdivision potential depends on the Planck potentials, which again depend on the temperature, pressure and chemical potentials. Like explained for the nanocanonical ensemble, it follows that the subdivision potential over the temperature can be measured,
i.e. by integrating the volume as a function over the pressure.
5. Canonical Ensemble:
In this case one controls the temperature, the volume, and the mole numbers of each replica. The totals of the volume and the mole numbers of the ensemble of replicas are therefore,
and
. Substitution in Eq.
1 gives
where the superscript C indicates that the corresponding variable depends on
and the shape of the system. The so-called
replica energy (the prefactor of
) is therefore given by
It follows from Eq.
24 that
This implies that
is the internal energy needed to subdivide the total entropy
over one more replica. Using that
and
are Euler homogeneous functions of the first order in
it follows from Eq.
24 that
The replica energy is therefore equal to the Helmholtz energy per replica. By substitution of Eqs:
25 and
27 into Eq.
24 we obtain the Gibbs equation for small systems
It has exactly the same form as for large systems. The important distinction is that for small systems
are not Euler homogeneous functions of
V. They have this property for large systems and in that case
, which implies that
. We define the integral Helmholtz energy for small systems by
This expression defines
as the sum of the compression and the chemical energies per replica. Together with Eqs.
25 and
27 it follows that
This property is typical for small systems as pointed out by Dong [
4]. Using Eqs.
28-
30 we obtain
It follows from Eq.
28 that
and from Eq.
31 that
The equations describe how the subdivision potential changes as a function of the temperature, pressure and chemical potentials. as discussed for the nanocanonical ensemble, a practical route may be to measure changes as a function of pressure. This can be found from the following integral
Similar relations are true for the dependence on the temperature and the chemical potentials. The important implication is again that the subdivision potential can be determined from measurements.
We have earlier [
6] calculated the subdivision potential of a multi-component ideal gas in a cube with surface adsorption. The rib of the cube had the length
L so that the volume equals
and the surface area is
. The result was
The corresponding value for
for a sphere is obtained by substituting 2
R for
L, where
R is the radius of the sphere. Again, this gives expressions which can be tested by molecular dynamics simulations.
In the thermodynamic limit the
in both expressions goes to zero [
6]. Hadwiger’s theorem [
7], which predicts that
is proportional to
in the case of adsorption on a flat surface, is clearly incorrect. We find that all powers of
contribute.
6. Isobaric Ensemble:
In this case we control the temperature, the pressure and the mole numbers of each replica. The totals of the mole numbers of the ensemble of replicas are therefore,
. Substitution in Eq.
1 gives
where the superscript P indicates that the corresponding variable depends on
and the shape of the system. The
replica energy (the prefactor of
) is therefore given by
It follows from Eq.
37 that
This implies that
is the internal energy needed to subdivide
and
over one more replica. Using that
and
are Euler homogeneous functions of the first order in
it follows from Eq.
37 that
The replica energy is therefore equal to the Gibbs energy per replica. By substitution of Eqs:
38 and
40 into Eq.
37 we obtain the Gibbs equation for small systems
It has exactly the same form as for large systems. The important distinction is that for small systems
are not Euler homogeneous functions of
. They have this property for large systems and in that case
, which implies that
.
We define the integral Gibbs energy for small systems like nuclei by
This defines
as the chemical energy per replica. Defining furthermore the integral chemical potentials by
, and using Eqs.
38 and
42 it follows that
This property is typical for small systems. Using Eqs.
41-
43 one obtains
It follows from Eq.
41 that
and from Eq.
44 that
Equation
46 expresses how the subdivision potential changes as a function of temperature, pressure and the chemical potentials. By integration we find, for instancė:
Similar expressions can be given for the dependence on the temperature and the chemical potentials. The subdivision potential is as a consequence clearly measurable.
7. Grand Canonical Ensemble:
We control the temperature, the volume, and the chemical potentials in the grand canonical ensemble. The total volume of the ensemble of replicas is,
. By substitution in Eq.
1 we obtain
where the superscript GC indicates that the corresponding variable depends on
and the shape of the volume. The so-called
replica energy (the prefactor of
) is therefore given by
It follows from Eq.
48 that
This implies that
is the internal energy needed to subdivide
and
over one more replica. Here
,
and
are Euler homogeneous functions of the first order in
. It follows from Eq.
48 that
In the last identity we introduced, following Hill, the integral pressure
, also called
p-hat. In the grand canonical ensemble there are two pressures: the integral pressure
which corresponds to the internal energy needed to increase the volume by adding a replica, and the common (differential) pressure which corresponds to the internal energy needed to increase the volume by increasing the volume of all the replicas keeping the total number of replicas constant.
By substitution of Eqs.
49 and
51 into Eq.
48 we obtain
which is Gibbs equation for small systems. It has exactly the same form as for large systems. The important distinction is that for small systems
are not Euler homogeneous functions of order one in
V. For large systems, however,
, which implies that
. It is convenient to define
By using Eq.
52, we obtain
It follows from Eqs.
49,
51 and
53 that
This property is typical for small systems. By using Eqs.
52-
54 we obtain
It follows from Eq.
54 that
and from Eq.
56 that
Equation
58 gives how the subdivision potential changes as a function of the temperature, pressure and chemical potentials. Equation
46 gives how the subdivision potential changes as a function of the temperature, pressure and the chemical potentials. By integration we find, for instancė:
Similar expressions can be given for the dependence on the temperature and the chemical potentials. It again follows that the subdivision potential can be obtained from measurements.
For the subdivision potential of a multi-component gas in a sphere with surface adsorption we obtain [
6]
In the thermodynamic limit
goes to zero. The subdivision potential divided by
V is proportional to
. The result for
is a scaling law. The result is in this particular case expected from Hadwiger’s theorem [
7].
We can find the results for adsorption on the walls of a sphere by substituting 2R for L, where we recall that R is the radius of the sphere. This leads to a different slope in the scaling law.
8. Small System Examples in the Grand Canonical Ensemble
We give below some small system examples in the grand canonical ensemble. We show how we can find the subdivision potential and its corresponding scaling law. We study three simple systems, which also can be and have been described without the subdivision potential. These are spherical and cylindrical systems, as well as a slit pore. Gibbs’ equation is available without the subdivision potential in these cases [
8], so we can compare the common to the alternative description. We continue with a system which can serve the purpose of a simple representative elementary volume of a porous medium. In this case, the subdivision potential provides a possible description of a control volume or representative elementary volume on the continuum scale.
Spherical systems
The radius of a sphere or cylinder decides on the energy of the fluid(s) inside the sphere [
6,
8,
9,
10]. When the surface is flat, the system is large, and the subdivision potential is zero [
9]. When the energy depends on the area and its curvatures, we can introduce the curvatures as additional variables. We can use the subdivision potential to describe geometrical effects [
6,
11]. The energy the arises from adsorption on a curved surface can also be described using the subdivision potential, as we shall see. This allows for a more general description of the surface energy, including heterogeneous materials.
The integral pressure of a spherical system in the grand canonical ensemble is given by
where we have assumed that the bulk pressure
is independent of
The integral surface tension
depends on the curvature
, and therefore on the volume. The surface area is
. The differential pressure is then
The resulting subdivision potential is now given by
Cylindrical systems
The integral pressure of a cylinder with a length
L and a radius
R in the grand canonical ensemble is likewise given by
where we assume that the bulk pressure
is independent of
The integral surface tension
depends on the curvature
, and therefore on the volume. The surface area is
. The differential pressure is
The resulting subdivision potential is given by
which is identical to the expression for the sphere.
Slit pore
A slit pore system with one fluid in the two-phase region was studied with two flat walls at a distance
d [
6,
9,
10]. For small
d this is a small system. When
d is large and their volume fractions sufficiently substantial the properties of the fluids and their equations of state are the same as for a bulk phase. There should not be complete wetting of one of the fluid phases. When one introduces the surface area as additional variable the description can be done without the use of a subdivision potential.
The integral pressure in a square section with two surface areas,
, of a slit pore with a distance
in the grand canonical ensemble is given by
where we assume that the bulk pressure
is independent of
The integral surface tension
depends on
h, and therefore on the volume. The total surface area is
. The differential pressure is
The resuling subdivision potential is given by
Again we find an expression identical to the expressions for the sphere. This is found when
R replaces
h.
A porous medium
Porous media have normally a pore size distribution. Pore diameters may be nanoscopic or microscopic, regular or not. The pores are by no means straight in reality. Numerous examples exist, where the fluid transport in such materials is of interest. On the pore scale, the porous medium often is a small thermodynamic system. We have proposed that a representative elementary volume (REV) for a continuum scale description can be constructed using a sufficiently large, but not larger volume element of the porous medium [
6].
Consider therefore a porous medium with a single phase fluid in a representative number of pores. The porosity is
. The pores are assumed to be in good approximation cylindrical with a radius
. In a representative elementary volume
V (REV) the pore volume is
. Here
is by definition the total length of the pores in the REV. The surface area of the pores is
. The integral pressure in the REV is
We assume that both
and
are independent of
V. Furthermore we assume that the porosity
is independent of
V. The changes in
V are isomorphic, i.e. if
V increases to
both
and
increase by a factor
. The differential pressure is
The resulting subdivision potential is given by
The expression is similar to those for the sphere, cylinder and the slit pore. The factor 2/3 accounts for the fact that the isomorphic expansion is in three dimensions, while the cross section of the pores is two-dimensional. The factor
is there because epsilon is proportional to the volume fraction of the pores.
We see that the subdivision potential contains typical terms that are major contributions to the capillary pressure, , in the presence of two fluid phases. This makes an interesting identification. The capillary pressure is accessible from two bulk liquid pressures and the surface tension.
The capillary pressure is defined in the case of mechanical equilibrium, by Young’s law. Here we obtain a thermodynamic quantity which is immediately can be useful in a non-equilibrium thermodynamic description of transport through the system. An overview of the characteristic length of each geometry is given in
Table 2.
9. Discussion
This paper explains Hill’s [
1,
2,
3,
6] thermodynamic method for small systems of five ensembles studied here. The focus has been on the subdivision potential, a quantity that plays a central role in the theory. We have shown that for all five ensembles, this is a measurable quantity, contrary to what has been claimed in the literature.
For three of the ensembles we have given statistical mechanical expressions for the subdivision potential, in the case of a multi-component ideal gas in a sphere or in a cylinder with adsorptions on the walls. As stated before [
6], these expressions show that Hadwiger’s theorem is not generally valid. They can apply under certain circumstances. Therefore, in general we cannot replace the subdivision potential by a Minkowski set of variables [
12]. The Minkowski set of variables contains the system’s volume, surface areas, and two integrals over surface curvatures [
12].
For large systems is small. This statement must not be confused with the statement that the subdivision potential is zero at equilibrium in small systems. Small systems in equilibrium have a finite subdivision potential, which depends on the choice of the control variables. The statistical mechanical expressions confirm that the subdivision potential in principle is not zero in equilibrium.
We use the name "nanothermodynamics" to name the extension of classical thermodynamics introduced by Hill. A summary of the ensemble dependency of the subdivision potential from this theory is given here in
Table 1. From the expressions in the table we conclude that subdivision potentials can be defined for each choices of control variables following Hill’s systematic procedure. They are particular for the ensemble or set of control variables, as shown, and their construction allow us to deal with shape and size as thermodynamic variables.
Table 1.
The subdivision potential for five ensembles of small systems [
6]. The shape of the small system is controlled in all ensembles
Table 1.
The subdivision potential for five ensembles of small systems [
6]. The shape of the small system is controlled in all ensembles
Ensemble |
Subdivision potential |
Expression |
Section |
NC |
|
|
2 |
MC |
|
|
3 |
C |
|
|
4 |
P |
|
|
5 |
GC |
|
|
6 |
Table 2.
The characteristic length for some small systems in the common scaling law of the subdivision potential.
Table 2.
The characteristic length for some small systems in the common scaling law of the subdivision potential.
Small system |
Characteristic length |
Symbol |
sphere |
Radius |
R |
cylinder |
Radius |
R |
cube |
Rib length |
L |
slit pore |
Half distance of walls |
h |
porous medium |
Average pore radius |
|
9.1. How Can We Measure and Compute the Subdivision Potential?
The subdivision potential is not a property that exists in theory only. It is equally practical as the chemical potential. Like the chemical potential, it can be used to describe how the small system deviate from ideal bulk behaviour. This has not been done systematically, however, for size and shape as variables. A purpose of this work has been to point at this possibility to find new information. For all ensembles, we have therefore presented expressions that allow us to find and use the subdivision potential from either measurements or computations. The subdivision potential is a function of temperature, pressure and chemical potential and must be determined as such.
We therefore claim that subdivision potential of a system is measurable for all ensembles. Well known procedures are relevant as discussed. The Maxwell relations [
6], that apply to the subdivision potential, are also useful for construction of experiments. We have not elaborated on these, because it carries too far. The subdivision potential divided by the volume will always become zero in the large system limit. This may help choose a practical standard state for the subdivision potential.
We have seen that the subdivision potential of a small system can be found from a statistical mechanical description of the small system, which is in essence the same as for larger systems. The partition function for a small systems includes all thermodynamic variables, and gives thus the equations of state for the small system [
13]. Such equations of state for small systems have been obtained already using molecular dynamics simulations.
9.2. Potential Applications
Nanothermodynamics has already gained some popularity. In solution theory, the common scaling laws for the thermodynamic variables and finite volume Kirkwood-Buff integrals [
6,
14,
15] have been recognized. This integral follows the scaling law
applicable to small systems, where
L is the characteristic length of the system [
16]. The origin of the scaling law is the dominant surface area contribution to the finite-volume Kirkwood-Buff integrals [
15]. The smaller the volume of the small system becomes, the higher powers of
must be included in the scaling law, due curvature and contact lines contributions.
We have seen above that it is possible to find the thermodynamic variables and corresponding equations of state for a large system, by extrapolating results for the small system as a function of
. This has been done for a few cases only [
16]. The equations presented for the subdivision potential in the examples above for the grand canonical ensemble are opening up an interesting road in solution theory.
10. Conclusions and Perspectives
We have seen in this article that the subdivision potential of Hill [
1,
2,
3] holds an interesting potential for future ways to describe energies on the small scale in a systematic manner. We have seen that we will be able to quantify small system’s properties, and discriminate between effect of sets of control variables by measuring or computing the subdivision potential. Thereby, we may understand how we can predict the system’s stability and energy conversion properties. Such predictions rest, however, on the ability to proceed with measurements and computations. It is the hope of the authors to have contributed to open up such a line of research, by clearing away some misconceptions in the literature.
Acknowledgments
This work grew out of an inspiring visit of Wei Dong to PoreLab, Trondheim, Norway. Discussions with Sondre Kvalvåg Schnell are appreciated. The authors are grateful to the Research Council of Norway for their Center of Excellence Funding Scheme, project no 262644 PoreLab.
References
- Hill, T. Thermodynamics of Small Systems. J. Chem. Phys. 1962, 36, 3182–3197. [Google Scholar] [CrossRef]
- Hill, T. Thermodynamics of Small Systems, Part 1 and 2; Benjamin, New York, 1963. [Google Scholar]
- Hill, T. Thermodynamics of Small Systems; Dover, New York, 1994. [Google Scholar]
- Dong, W. Nanoscale thermodynamics needs the concept of a disjoining chemical potential. Phys. Rev. Letters. 2023, 10, 1038. [Google Scholar] [CrossRef] [PubMed]
- Bedeaux, D.; Kjelstrup, S. Hill’s nano-thermodynamics is equivalent with Gibbs’ thermodynamics for surfaces of constant curvatures. Chem. Phys. Letters 2018, 707, 40–43. [Google Scholar] [CrossRef]
- Bedeaux, D.; Kjelstrup, S.; Schnell, S. Nanothermodynamics. Theory and Applications; World Scientific, 2024. [Google Scholar]
- Hadwiger, H. Vorlesungen über Inhalt, Oberflache und Isoperimentrie; Springer: Berlin, 1957. [Google Scholar]
- Bedeaux, D.; Kjelstrup, S. Hill’s nano-thermodynamics is equivalent with Gibbs’ thermodynamics for curved surfaces. Chem. Pys. Letters 2018, 707, 40–43. [Google Scholar] [CrossRef]
- Rauter, M.; Galteland, O.; Erdos, M.; Moultos, O.; Vlugt, T.; Schnell, S.; Bedeaux, D.; Kjelstrup, S. Two-Phase Equilibrium Conditions in Nanopores. Nanomaterials 2020, 10, 608. [Google Scholar] [CrossRef] [PubMed]
- Erdos, M.; Galteland, O.; Bedeaux, D.; Kjelstrup, S.; Moultos, O.; Vlugt, T. Two-Phase Equilibrium Conditions in Nanopores.
- Olav Galteland, D.B.; Kjelstrup, S. Nanothermodynamic description and molecular simulation of a single-phase fluid in a slit pore. 2021; 11, 165. [Google Scholar]
- Armstrong, R.; McClure, J.; Robins, V.; Liu, Z.; Arns, C.; Schlüter, S.; Berg, S. Porous media characterization using Minkowski functionals: Theories, applications and future directions. Transport in Porous Media 2019, 30, 305–355. [Google Scholar] [CrossRef]
- Bråten, V.; Bedeaux, D.; Wilhelmsen, ∅.; Schnell, S. Small size effects in open and closed systems: What can we learn from ideal gases about systems with interacting particles? J. Chem. Phys. 2021, 155, 244504–17. [Google Scholar] [CrossRef] [PubMed]
- Krüger, P.; Scnell, S.; Bedeaux, D.; Kjelstrup, S.; T.J.H, V.; J.M.Simon. Kirkwood-Buff integrals for finite volumes. J. Phys. Chem. Letters 2013, 4, 235–238. [Google Scholar] [CrossRef] [PubMed]
- Simon, J.; Krüger, P.; Schnell, S.; Vlugt, T.; Kjelstrup, S.; Bedeaux, D. Kirkwood-Buff Integrals: from fluctuations in finite volumes to the thermodynamic limit. J. Chem. Phys. 2022, 157. [Google Scholar] [CrossRef] [PubMed]
- Bjørn Strøm, Sondre K. Schnell, J.M.S.J.H.S.K. Size and shape effects on the thermodynamic properties of nanoscale volumes of water. Phys. Chem. Chem. Phys. 2017, 19, 9016–9027. [CrossRef] [PubMed]
1 |
Equation 13 in [ 6] is obtained from Eq.6.26. We take the opportunity to correct the wrong formula in Eq.6.27 in [ 6]. In the thermodynamic limit goes to zero because goes to infinity in that limit. |
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).