Nanothermodynamic Description and Molecular Simulation of a Single-Phase Fluid in a Slit Pore

We have described for the first time the thermodynamic state of a highly confined single-phase and single-component fluid in a slit pore using Hill’s thermodynamics of small systems. Hill’s theory has been named nanothermodynamics. We started by constructing an ensemble of slit pores for controlled temperature, volume, surface area, and chemical potential. We have presented the integral and differential properties according to Hill, and used them to define the disjoining pressure on the new basis. We identified all thermodynamic pressures by their mechanical counterparts in a consistent manner, and have given evidence that the identification holds true using molecular simulations. We computed the entropy and energy densities, and found in agreement with the literature, that the structures at the wall are of an energetic, not entropic nature. We have shown that the subdivision potential is unequal to zero for small wall surface areas. We have showed how Hill’s method can be used to find new Maxwell relations of a confined fluid, in addition to a scaling relation, which applies when the walls are far enough apart. By this expansion of nanothermodynamics, we have set the stage for further developments of the thermodynamics of confined fluids, a field that is central in nanotechnology.


Introduction
The thermodynamic state of a fluid in confinement is important for the understanding of adsorption to walls, chemical reactions, film formation and transport in porous media [1][2][3][4][5][6]. The molecular structuring at the walls and the forces between particles and walls are central. The change in thermodynamic properties upon confinement is substantial. This has been known for a long time [1,7]. Derjaguin considered the measurable force that attracts or repels two walls that are close together, and defined from this the disjoining pressure [7]. The disjoining pressure has also been called the solvation pressure [1]. When the walls are far apart, the disjoining pressure vanishes. It is not well known, however, how size and shape, as variables, affect the disjoining pressure or other properties of the confined fluid.
Confinement is considered to be important, for instance, in the context of CO 2 separation and sequestration by metal-organic frameworks [2] or for adsorption in zeolites [3]. The disjoining pressure is of interest when studying aggregation of colloidal particles, suspended or adsorbed [4][5][6]. It is likely to be important also for film flow on the macroscale [8].
More knowledge of confined fluids on the nanoscale is therefore needed. It may, for instance, help us solve the well-known up-scaling problem in porous media science [9]. The central problem is to understand how to integrate properties on the pore scale to the macroscale where Darcy's law applies. In order to account for shape and size effects, it was recently proposed to use the four Minkowski functionals [10][11][12]. This simplifies the description of a representative elementary volume (REV). Another procedure using the grand potential for average variables in the REV was also proposed [13,14]. surface area 2Ω t = 2 ∑ N j=1 Ω j , and total number of particles N. The factor of two in the total surface area arises because there are two fluid-solid surfaces of equal area per slit pore. By construction, the ensemble variables U t , S t , V t , Ω t , N, and N are Euler homogeneous functions of the first order in the number N of slit pores. The total differential of the total internal energy is This type of equation for the total internal energy we call the Hill-Gibbs equation [25]. The last term was added by Hill. The partial derivatives of the total internal energy define the temperature T, the normal pressure p ⊥ , the surface tension γ and the chemical potential µ The control variables in subscripts are kept constant while taking the derivatives. The volume derivative is taken while keeping the total surface area constant, which implies that the volume is changed by changing the distance between the surfaces h j ≡ V j /Ω j . The surface derivative is taken while keeping the volume constant, which implies that both the pore heights and surface areas are changed in such a way that the change in the total volume is zero.
The new thermodynamic variable ε is the subdivision potential. It is The subdivision potential is defined here as the increase in the total internal energy as the number of slit pores N increases while keeping S t , V t , Ω t , and N t constant. This definition is different from the definition in a previous article by us [23], where only the entropy, volume and number of particles were kept constant and not the surface area. This led to a different expression for the subdivision potential. The subdivision potential is the work done on the system when adding a new slit pore while keeping the other control variables constant. The subdivision potential may be positive or negative, depending on whether work is needed or gained by adding new slit pore replicas.
We are aiming to describe the disjoining pressure, see Section 2.3, of an open slit pore when the volume, surface area, temperature and chemical potential are control variables. We will change to this ensemble. This ensemble is useful for describing experiments and simulations. We use the average volume per slit pore V = V t /N and the average surface area per slit pore 2Ω = 2Ω t /N , rather than the total volume V t and the total surface area 2Ω t . In order to obtain an appropriate Hill-Gibbs equation for this case, we substitute the total volume with the average volume V t = VN and total surface area with the average surface area 2Ω t = 2ΩN . The total differentials of the volume and surface area are By introducing this into the Hill-Gibbs equation, see Equation (1), we obtain The parenthesis define the replica energy The subdivision potential ε, normal pressure p ⊥ , and surface tension γ depend on the control variable set T, V, Ω, and µ. The replica energy density will be used to define the disjoining pressure. The grand partition function covers all microstates available to the slit pore. In order to calculate this partition function one chooses a volume V and a surface area Ω. Both can be varied independently. For large h and Ω only a change in the volume V = hΩ, and not in h and Ω separately, matters. For a small volume, it is necessary to use the volume V and the surface area Ω as independent variables. This is because the same volume change due to a change of the pore height h = V/Ω or due to a change of the surface area Ω produces different changes in the partition function, and therefore in the thermodynamic variables.
The replica energy was identified by Hill as the grand potential, here equal to minus the integral pressure times the volume [26], In the last equality, we have chosen to identify the integral pressure as the integral normal pressure times the volume minus the integral surface tension times the surface area. We will refer top ⊥ ,γ andp as the integral normal pressure, the integral surface tension and the integral pressure, respectively, while p ⊥ and γ are the differential normal pressure and the differential surface tension, respectively. The names integral and differential pressure were coined by Hill to reflect that the differential pressure involves the differential of the integral pressure. We have chosen a control variable set with volume and surface area, such that we do not have a differential pressure but a differential normal pressure p ⊥ and differential surface tension γ in its place. From Equations (6) and (7) it follows that the subdivision potential is The subdivision potential ε indicates that the integral normal pressure and the integral surface tension may be different from the corresponding differential variables.

Maxwell Relations for a Slit Pore
Using that the total internal energy is Euler homogeneous of the first order in the number of slit pores, see Equation (5), we integrate it at constant T, V, µ and Ω We introduce the average internal energy, entropy and number of particles per slit pore By introducing the average properties into the total internal energy in Equation (9) we obtain the internal energy per slit pore Substituting the average properties per slit pore in Equation (10) into the corresponding Hill-Gibbs equation, see Equation (5), and using the internal energy in Equation (11), we obtain the total differential of the internal energy By differentiating the internal energy in Equation (11) and using the total differential of the internal energy in Equation (12), we obtain the total differential of the replica energy This equation was termed the Hill-Gibbs-Duhem Equation [25], because it reduces to the Gibbs-Duhem equation for a large system. It follows that the partial derivatives of the grand potential is Rather than the names replica energy or grand potential, we name from now X = −pV by minus the integral pressure times the volume. In Section 2.4 we will define the integral pressure in terms of the average tangential mechanical pressure which can be calculated from molecular simulations. In the simulations, we consider surface areas Ω much larger than the diameter of the fluid particles and heights h = V/Ω comparable to the diameter of the fluid particles. This implies thatp,p ⊥ ,γ, u ≡ U/V, s ≡ S/V, p ⊥ , γ, and ρ ≡ N/V do not depend on Ω, but the variables will depend on the height h and therefore on the volume V. The volume and surface derivative can be rewritten in terms of derivatives of the slit pore height and surface area. These derivatives are needed to calculate the differential normal pressure and differential surface tension. The differential normal pressure is The differential surface tension is Combining the equations for the differential normal pressure and surface tension it follows that the integral pressure isp The last term is the equal to minus the subdivision potential divided by volume, see Equations (7) and (8). The subdivision potential is This implies that in the grand canonical ensemble with T, V, Ω, µ as control variables, the thermodynamic description of the slit pore with a small height is the same as for the slit pore with a large height. However, it changes when the integral pressure depends on the surface area. This may not be the case for any other set of control variables. In general, the properties of a small system (confined fluid) depend on the set of control variables.
In this work, we deal with large surface areas, such that the integral pressure does not depend on it. As a consequence the subdivision potential is zero. We will find that the integral and differential normal pressure are equal. Using that ε = 0, it follows that the integral and differential surface tensions are also equal for a large surface area, see Equation (8), We will investigate the first identity for large surface areas using molecular simulations. From Equation (13) we obtain the following Maxwell relations of the differential surface tension, In the second identity, we have used that the entropy density, fluid number density and differential normal pressure are independent of the surface area for large surface areas.
Other Maxwell relations are possible, see the Hill-Gibbs-Duhem Equation (13). The last equality can be written as The important implication of this expression is that the constant nature of one variable implies the constant nature of the other for large pore heights. With the mechanical description of the integral properties given in Section 2.4 the consistency of all the above thermodynamic relations can be tested.

The Disjoining Pressure
We define the excess replica energy as the replica energy of a slit pore of height h minus the replica energy of the slit pore where the slit pore height approaches infinity, We denote the thermodynamic variables where the slit pore height approaches infinity with ∞ in superscript, The excess replica energy can be written in terms of the excess integral pressurep ex , excess integral normal pressurep ex ⊥ , and excess integral surface tensionsγ ex We define now the disjoining pressure as the excess normal pressure, Other possible names are minus the excess replica density, or the excess grand potential density. This definition of the disjoining pressure is different from the typical definition found in the literature [7,29,30]. Typically the disjoining pressure is defined to be equal to what we in this work call the excess differential normal pressure, where the excess is relative to a bulk fluid. For large surface areas we will show with molecular simulations that p ⊥ = p ⊥ . We will furthermore show that the integral normal pressure, as the slit pore height approaches infinite separation, equals the bulk fluid pressure p ∞ ⊥ = p b . This shows that our choice is equivalent to the usual definition.

A Mechanical Description of the Slit Pore
The thermodynamic description of integral and differential pressures and surface tensions has its mechanical equivalent description in terms of components of the mechanical pressure tensor. A recent discussion on this topic clarified the challenge of translating the mechanical pressure tensor into a thermodynamic scalar variable in a meaningful manner [31][32][33]. We will identify the integral normal pressure, integral surface tension, and integral pressure in such a way that thermodynamic framework is self-consistent. However, we do not claim that this is the only valid choice of thermodynamic pressures and tensions in terms of the mechanical pressure tensor.
The mechanical pressure tensor of a heterogeneous system is ambiguously defined. This has been known for a long time, at least since the work by Irving and Kirkwood in the 1940s. It was shown by Schofield and Henderson [34] that the ambiguity is due to the arbitrary choice of the integration contour C ij , which is needed to calculate the configurational contribution to the pressure tensor. The local mechanical pressure tensor is calculated in a subvolume V l as a sum of the kinetic and the configurational contributions, Upper case P is used to denote mechanical pressure tensor components in order to distinguish them from the thermodynamic pressures, which are denoted by the lower case p. The kinetic contribution is the ideal contribution to the mechanical pressure and is calculated as The sum with subscript i ∈ V l represents a sum over all particles in the subvolume V l . The particle mass is m i and v i,α is the velocity in the α-direction. The solid walls do not have a velocity and consequently do not directly contribute to the kinetic pressure. The brackets . . . represent ensemble average. The configurational contribution is the non-ideal contribution to the mechanical pressure and is calculated as The sum represents a sum over all particle pairs. The α-component of the force vector acting on particle i due to particle j is f ij,α . The fluid-fluid and fluid-solid interactions contribute to the configurational pressure, the solid-solid interaction is zero and does not contribute to the pressure. The line integral is the β-component of the part of the contour C ij contained in the subvolume V l .
The contour C ij is the source of the ambiguity of the mechanical pressure tensor, it can be any continuous line from the centers of particles i to j. The Harasima [35] and the Irving-Kirkwood [36] contours are two common choices for C ij . The Harasima contour is defined as two continuous line segments, a line from the center of particle i parallel to the surface and a line normal to the surface to the center of particle j. The Irving-Kirkwood is the straight line from particle i to j. For flat surfaces they are equal. However, for spherical surfaces, the Harasima contour does not obey momentum balance [37]. There are cases where the Harasima contour is more useful than the Irving-Kirkwood contour [38]. In this work, we will use the equations by Ikeshoji et al. [39] with the Irving-Kirkwoood contour to calculate the mechanical pressure tensor.
We will only consider surface areas Ω much larger in both directions than the fluid particle diameter. This implies the mechanical pressure tensor does not depend on Ω. It does, however, depend on the height h ≡ V/Ω and therefore on the volume V.
Using the translational symmetry of the slit pore in the yand z-direction, where the x-direction is normal to the solid surface, the equilibrium mechanical pressure tensor in the slit pore has the form where e x , e y and e z are the unit vectors in the x-, yand z-directions. The normal pressure tensor component is equal to the xx-component, and tangential pressure tensor component is the average of the yyand zz-components, Mechanical equilibrium requires that the tangential pressure is independent of the yand z-coordinates, but depends on the x-coordinate. The normal pressure is independent of all spatial coordinates. We identify the thermodynamic integral normal pressure in terms of the volume integral of the normal mechanical pressure divided by the volume. However, since the normal mechanical pressure and the area are constant everywhere, this simplifies tô The integral normal pressure in a large pore is equal to the pressure in a bulk phase in equilibrium with the pore. The difference in the integral normal pressure in a liquid and vapor phase in a slit pore is described the Young-Laplace equation, while the integral pressure is the same in both the liquid and vapor phase [23]. An alternative route to the integral normal pressure is via the local fluid number density profile [40] where f f s (x) is the fluid-solid force and ρ(x, h) is the local fluid number density. The fluid density and fluid-solid forces are uniquely defined, and do not have the inherent problem that the mechanical pressure tensor has. We can use this equation to validate our method of calculating the integral normal pressure. We identify further the integral surface tension as the integral of the normal minus tangential pressure tensor components, where the factor half is due to the fact that there are two fluid-solid surfaces. It follows from Equation (7) together with Equations (31) and (33) that the integral pressure is, As shown by Harasima [35] and Schofield and Henderson [34] a sufficiently large volume integral of the mechanical pressure tensor components does not depend on the choice of the integral contour C ij . We have identified the integral normal pressure, integral surface tension and integral pressure in terms of the mechanical pressure tensor. The local mechanical pressure tensor has an inherent problem, specifically that the contour C ij can be any continuous line from i to j. However, the thermodynamic variables do not have this inherent problem. This is because we integrate the local mechanical pressure tensor across the whole volume V. This volume integral includes all interactions. The internal energy can be calculated as the sum of the kinetic and potential energy, By dividing the internal energy in Equation (11) by the volume we obtain the internal energy density u = U/V. We also use the entropy density s = S/V, and fluid number density ρ = N/V. By rearranging the equation we obtain the entropy density as where the internal energy density, integral pressure and fluid number density are known. The volume, chemical potential and temperature are imposed on the system.

Simulation Details
The thermodynamic state of slit pores of varying heights h was investigated by using grand-canonical Monte Carlo (GCMC) [41] in combination with molecular dynamic (MD) simulations with the Nosé-Hoover thermostat [42]. This produced the grand canonical ensemble, i.e., constant chemical potential, temperature, volume and surface area. The GCMC method inserted and deleted fluid particles to and from the simulation box from an imaginary fluid particle reservoir at the same temperature and chemical potential. This controlled the chemical potential of the fluid in the slit pore. The MD procedure updated the positions and velocities of the fluid particles and controlled the temperature with the Nosé-Hoover thermostat.
The simulations were carried out using LAMMPS [43]. The local mechanical pressure tensor was calculated by post-processing the particle trajectories with in-house software (available at D.O.I. 10.5281/zenodo.4405267). The chemical potential and temperature were kept constant at µ * = 1 and T * = 2. The critical temperature of the Lennard-Jones/spline fluid is T * c = 0.885 [44]. All units in this work are in reduced Lennard-Jones units, see Table 1 for a definition. Table 1. The reduced units are denoted with an asterisk in superscript, for example T * . The variables are reduced using the molecular diameter σ, potential well depth , fluid particle mass m and Boltzmann constant k B .

Description
Definition The simulation box was a rectangular cuboid of side lengths L x , L y = L z . The side lengths L y = L z were chosen such that the surface area was large. Large in this context indicates large enough forp,p ⊥ ,γ, p ⊥ , and γ to be independent of the surface area Ω = L y L z . The simulation box size was decided such that the average number of fluid particles was approximately 2 × 10 4 . The simulation box was periodic in the yand zdirections, and non-periodic in the x-direction. This implies that the particles did not interact across the simulation box boundary in the x-direction.
The fluid-fluid and fluid-solid interaction was modeled with the Lennard-Jones/spline potential [44]. The fluid-fluid and fluid-solid interactions were equal. The potential energy of a fluid-fluid or fluid-solid pair separated by a distance r was where r s , a, b and r c were chosen such that the potential energy and the force were continuous at the inflection point r = r s and the cut-off r = r c . The solid walls were placed at the simulation box boundaries x = −L x /2 and x = L x /2. The distance between the fluid and solids were r = |x f − L x /2| and r = |x f + L x /2|, where x f is the x-position of the fluid particle. The dividing surfaces of the fluid-solid surfaces were chosen to be at x = −L x /2 and x = L x /2. The slit pore height was consequently determined to be h = L x . Other choices of the dividing surface are possible, for example, the Gibbs dividing surface or the surface of tension. When L x < 2r c the fluid particle can interact with both solid walls. See Figure 1 for an visualization of the simulation box for the case L x = 4σ. The mechanical pressure tensor was calculated in thin rectangular cuboids, called layers l, of side lengths ∆x, L y , L z . The thickness of the layers was ∆x = 0.005σ and the number of layers was n l = h/∆x. The diagonal components of the mechanical pressure tensor was calculated using Equations (27)- (29).
The kinetic energy was calculated as the sum of the kinetic energy for each fluid particle and the potential energy was calculated as the sum of the potential energy of each fluid-fluid and fluid-solid pairs, The sums of the kinetic and potential energies were used to calculate the internal energy and entropy densities.

Results and Discussion
The results are presented in Figures 2-10 and discussed in that order before general remarks are offered.
The normal mechanical pressure P ⊥ is presented in Figure 2a. It does not depend on the position x, but it depends strongly on the slit pore height h. The figure shows a straight line of various lengths for each of the three heights, which reflect the slit pore height h. The integral normal pressure was identified as this component,p ⊥ = P ⊥ . We see that it is always constant, as demanded by Equation (31). For slit pore heights h > 7σ we find that the normal mechanical pressure is equal to the bulk pressure. At h = 2.04σ the normal mechanical pressure is at a global maximum and at h = 2.59σ it has a local minimum for the given temperature and chemical potential. The normal mechanical pressures divided by the bulk pressure for the two cases are P ⊥ /p b = 2.5135 ± 0.0008 and P ⊥ /p b = 0.5539 ± 0.0002, respectively. The global minimum, which is zero, is at h < 1.8σ when no fluid particles fit in the slit pore. This is because solid-solid interactions and quantum effects are not considered in this work. The tangential mechanical pressure P , illustrated in Figure 2b, depends in contrast on the position x as well as on the slit pore height h. The integral pressure is the average of P , see Equation (34). The tangential mechanical pressure follows the trend of the fluid number density, compare Figures 2b,c. For pore sizes h > 7σ the tangential mechanical pressure is constant and equal to the bulk pressure p b in the center of the pore. This indicates that the pore is large enough to accommodate bulk liquid in the center. The fluid is highly structured close to the fluid-solid surface [1]. As the slit pore height is decreased, fluid structures on the two sides overlap. When regions of structured fluids overlap, repulsive and attractive forces between the surfaces appear, and the disjoining pressure becomes non-zero.
The fluid number density ρ = N/V is presented in Figure 3 as a function of the slit pore height h. The bulk fluid number density ρ b in equilibrium with the slit pore is shown as a dashed line. The fluid number density converges to the bulk value as the slit pore height approaches infinity. The volume V depends on the choice of the fluid-solid dividing surfaces. We have chosen the dividing surfaces to be at x = −L x /2 and x = L x /2. The choice of the dividing surface determines how rapidly the slit pore values converge to the bulk values. A dividing surface closer to the fluid phase will reduce the volume and consequently the slit pore values will converge faster to the bulk values. Other choices of the dividing surface are possible. The entropy density is presented in Figure 4 as a function of the slit pore height h. The entropy density is a monotonically increasing function of the height h. This confirms the observation by Israelachvili [1] that the origin of the oscillations of the disjoining pressure as a function of the height is not entropic. As a further confirmation of this point, we find that the internal energy density oscillates with a period equal to the particle diameter, see Figure 5. The oscillating forces or pressures are thus of energetic origin. The bulk entropy and internal energy densities are shown as dashed lines. The entropy and internal energy densities of the slit pore converge to the bulk values as the slit pore height is increased.  The integral pressure is equal to the volume average tangential mechanical pressurê p = h −1 h 0 P dx. It is of special interest because it is equal to minus the grand potential divided by volumep = −X/V, or the replica energy density. The grand potential is the starting point for the definition of the average thermodynamic properties of the REV [13]. The integral pressure is presented in Figure 6 as a function of the slit pore height h. The integral pressure converges to the bulk pressure p b as the slit pore height approaches infinity, as expected. The bulk pressure is shown as a dashed line.
In previous works [13,14,21] we argued that the gradient of the integral pressure is the driving force for mass flux. In another work [23] we found the integral pressure of a two-phase system in a slit pore to be equal in the liquid and vapor in equilibrium. The identification of the integral pressure in this work is consistent with this interpretation. The gradient of the integral pressure is the driving force of the mass flux. For fluid flows tangential to the slit pore surfaces it is the gradient in the tangential mechanical pressure tensor component that gives the driving force when the system is out of equilibrium. In this work we identify the integral pressure as the average of the tangential mechanical pressure tensor components.
As stated in Section 2.2, the integral and differential normal pressures and integral and differential surface tensions are expected to be equal when the surface area is large. If this is correct, the integral pressure can be computed aŝ The lower integration limit is h 0 = 1.8σ, at which point the integral pressure is in good approximation zero. The integral pressure computed from Equation (39) is shown in Figure  6. The curves are identical. This implies that the integral and differential normal pressures are equal. As we have already shown that the subdivision potential is zero ε = 0 for large surface areas with this set of control variables, it follows that the integral and differential surface tensions are also equal. We will from now on refer to the integral normal pressure and integral surface tensions as the normal pressure and surface tension. The normal pressure was identified as the normal mechanical pressure p ⊥ = P ⊥ in Equation (31). It is presented in Figure 7 as a function of the slit pore height h. The normal pressure was also calculated from the local fluid density and the fluid-solid force using Equation (32). The two methods of calculating the normal pressure agree, which indicates that we have calculated the mechanical pressure tensor correctly.
The normal pressure oscillates with a period equal to the fluid particle diameter at small heights h. The oscillations decay as the height h increases. Such oscillations have been observed in experiments and are well known, see for example Israelachvili [1]. The oscillations are caused by the structuring of the fluid particles between the surfaces, and by the fact that the fluid particles all have the same diameter. As the height is increased above h > 7σ the oscillations vanish and the normal pressure is constant and equal to the bulk pressure. The bulk pressure is shown as a dashed line in the figure. At heights h > 7σ the fluid structuring near the walls do not overlap. At lower densities, smaller oscillations are expected with a faster decay. The normal pressure shows a similar trend to previous works [15,16,40]. At very small heights, the solid-solid interaction will dominate and completely overshadow the fluid-fluid and fluid-solid interactions presented here. We have not included any solid-solid interaction in this work, and as a consequence the normal pressure approaches zero because there is no room for any fluid particles to enter the slit pore. When p ⊥ =p ⊥ and ε = 0, it follows from Equation (8) that the integral and differential surface tensions are equal γ =γ. The surface tension as a function of the slit pore height is presented in Figure 8, see Equation (33). The surface tension at infinite separation γ ∞ is computed as the average surface tension of the slit pore, with height h > 10σ, at which point the surface tension is independent of h. The disjoining pressure Π was computed from Equation (25), and is shown as a function of the slit pore height h in Figure 9. The disjoining pressure was here defined to be equal to the integral normal pressure minus the normal pressure at infinite separation. Because the integral and the differential normal pressures are the same and because the normal pressure at infinite separation is equal to the bulk pressure, in this case, the definition contain the commonly used definition of the disjoining pressure [7,29]. The normal pressure is constant for slit pore heights h > 7σ. Consequently, the normal pressure at infinite separation can be calculated as the normal pressure when h > 7σ. The normal pres-sure at infinite separation is equal to the pressure in a bulk fluid with the same temperature and chemical potential.
For the present case, we can claim that our definition of the disjoining pressure is equivalent to the common definition. Our definition is general, as it also covers the cases where the integral and the differential normal pressures are unequal. Examples where this is the case are given below. Figure 10 shows how the normal pressure minus the integral pressure scales with inverse slit pore height. It is equal to the scaling law presented in a previous work [23]. The slope of this ideal curve is equal to two times the surface tension. When the inverse slit pore height approaches zero, i.e., when the walls are far apart, the normal pressure and integral pressures are equal as predicted. For pores that are so small that no bulk fluid can form in the center, the structuring at the walls starts to overlap, and a fluctuating difference is seen in the difference of the normal and integral pressures. In the region of the straight line, we have a scaling law, that relates states of different heights. At heights smaller than approximately h < 5σ the scaling law breaks down, the difference of the two pressures starts to oscillate. There are positive and negative deviations from the law.  In our earlier work [23], we studied liquid-vapor coexistence in a slit pore. In that work the surface area was not a control variable, and as a consequence the subdivision potential was found to be equal to two times the surface tension divided by the slit pore height, see Figure 10. Another reasonable set of control variables is the height h instead of the volume V. For a control variable set consisting of temperature, height, surface area and chemical potential, the differential normal pressure is equal to The differential surface tensions is For this set of control variables, the height h is kept constant instead of the volume V.
The subdivision potential is accordingly These relations also help us characterize the smallness of the slit pore with the large walls. The subdivision potential, introduced by Hill as a measure of smallness, deviates from zero in the last relation, also when the integral pressure does not depend on the size of the area, with height h and surface area Ω as control variables. A dependency on the area is relevant when adsorption takes place on small spheres [24]. A slit pore with large walls may be expected to be small for small heights h, since the confined fluid is not bulk-like. However, we have seen that when we use volume V and surface area Ω as control variables, the subdivision potential is zero for large surface areas. A zero subdivision potential means that the system also can be described perfectly using regular thermodynamics [18]. It is nevertheless meaningful to define a non-zero integral pressure, because the integral pressure enters the grand potential. It, therefore, determines the thermodynamic properties of a REV. Away from equilibrium, it will create a driving force. The grand potential or minus the integral pressure times the volume are equal to the replica energy. The replica energy is not zero in the present case. Clearly, we have here an example where smallness is brought out in Hill's terms through the replica potential.

Conclusions
We have developed a nanothermodynamic description based on the ideas of Hill to describe single-phase and single-component fluids in slit pores in a new manner. As environmental control variables we chose the chemical potential, temperature, volume and surface area. We have seen that the outcome varies with the set chosen, but the procedure can be used for complex geometries and different sets.
Following Hill, we introduced the subdivision potential. It is non-zero only when the integral pressure depends on the surface area Ω. For large surface areas we have shown that the subdivision potential is zero, and that p ⊥ =p ⊥ and γ =γ. In this sense, we have found that nanothermodynamics is equivalent to the usual thermodynamic description for all heights. The replica energy, and therefore the integral pressure, were shown to be non-zero. This allowed us to identify a scaling law, which confirms earlier results [22,23]. By choosing height h rather than volume V among the control variables, a non-zero subdivision potential appears.
We have identified the thermodynamic properties by their mechanical counterparts in a consistent manner. The integral pressure, which is equal to minus the grand potential divided by volume, can be understood as the average tangential mechanical pressure. The normal pressure is the normal mechanical pressure, and the surface tension is the integral of the normal minus the tangential mechanical pressure. The entropy and internal energy densities vary with the slit pore height, confirming the observations of Israelachvili [1]. The entropy density increases monotonically with increasing height, while the energy density oscillates. This confirms that the disjoining pressure is not of entropic origin, it is of energetic origin [1].
By these investigations of the nanothermodynamic theory and the subsequent simulations, we hope to have expanded on the knowledge on Hill's method, making it more available for further studies, for instance of flow and reactions in porous media.
Author Contributions: O.G. contributed to investigation, methodology, software and visualization. D.B. and S.K. contributed to supervision. All authors contributed to conceptualization, formal analysis, writing original draft, reviewing and editing. All authors have read and agreed to the published version of the manuscript.