Multilayer approximation for a confined fluid in a slit pore

 

G. J. Zarragoicoechea,^{1, 2} A. G. Meyra,¹ V. A. Kuz¹

*Email: vasco@iflysib.unlp.edu.ar

IFLYSIBInstituto de Física de Líquidos y Sistemas Biológicos (CONICET, UNLP, CICPBA), 59 No. 789, 1900 La Plata, Argentina.

CICPBAComisión de Investigaciones Científicas de la Prov. de Buenos Aires.

A simple LennardJones fluid confined in a slit nanopore with hard walls is studied on the basis of a multilayer structured model. Each layer is homogeneous and parallel to the walls of the pore. The Helmholtz energy of this system is constructed following van der Waalslike approximations, with the advantage that the model geometry permits to obtain analytical expressions for the integrals involved. Being the multilayer system in thermodynamic equilibrium, a system of nonlinear equations is obtained for the densities and widths of the layers. A numerical solution of the equations gives the density profile and the longitudinal pressures. The results are compared with Monte Carlo simulations and with experimental data for Nitrogen, showing very good agreement.

 

I. Introduction

The effects on phase transition of confined fluids in a slitlike pore have been studied by simulation and different theories [111]. In a previous work, we constructed a generalized van der Waals equation for a fluid confined in a nanopore [12, 13]. The shift of the critical parameters was in good agreement with lattice model and numerical simulation results, and the predicted critical temperature remarkably reproduced the experiment. In that work, we concluded that the confined van der Waals fluid theory seemed to work better than the bulk one, maybe due to the fact that the higher virial contributions not considered in both theories were less important in the confined fluid than in the bulk. A similar treatment was used previously by Schoen and Diestler [14]. Following that line of reasoning, here we study a simple fluid confined between two infinite parallel hard walls (slit pore). The walls are at a distance L apart. To study the confined fluid, we propose a multilayer model [15]: the fluid is distributed in n thin layers, one beside the other. Each layer has a uniform density, and can be observed as a nonautonomous phase. A particle in a given layer interacts with its neighbors inside the layer, and with every particle in the other layers. Defay and Prigogine and Murakami et al. have shown that, in a liquid gas interface, the deviation from the Gibbs’ adsorption equation becomes practically negligible in the case of a two layer model [16], and that as the number of transition layers grows, the multilayer model becomes perfectly consistent with the Gibbs’ equation [17]. The van der Waalslike approximations made in developing this multilayer model theory limit its validity to the low density regime.

II. Theory

 

The configuration integral Z_N for a pair potential v_ij may be approximated as

and further expanded in function of two partióle integráis

The first term in Eq. (3) stands for particles in the layer i. The second term comes from the interaction of one particle in layer i with one particle in layer j. In a compact form, and assuming that a layer sees three nearest neighbor layers,

The integrals I_i and I_ij , for the slit pore geometry and after low density approximations, can be analytically solved to give

In the above expressions v_ij was taken to be the LennardJones pair interaction, being e and a the potential parameters. The integráis 1^ ¿₊2 and Ií,í+3 do not contain the excluded volume term because we suppose that the layer widths are a. The expressions for A and B in the preceding equations, functions of L_xiy are given in the Appendix A.

The Helmholtz energy Eq. (1) together with Eq. (4), has the final expression

The pressure tensor [12, 13] and chemical potentials are obtained from the following equations

If the system is in mechanical and chemical equilibrium, the xx components of the pressure tensor and the chemical potentials for each layer must be equal. From these equations, giving as input the wall separation L and the mean density p* = per³, it is constructed a system of (n 1) nonlinear equations with (n 1) unknowns (layer densities and widths) to be numerically solved. The low computational cost is taken for granted given that the code is easily written and the calculations are carried out on a Pentium 4 processor running at 2.66 GHz. At a temperature T*=kT/£=1, we have explored the cases with L=10a and L=15cr, at different mean densities. We have also compared the theoretical results with experimental data coming from studies of nitrogen adsorption in graphite slit pores at room temperature [19].

III. Monte Cario simulation

For numerical simulations, N LennardJones particles are confined between hard walls separated at a distance L. The unit cell is build up taken the walls to be of size L_y and L_z in the y and z directions respectively, directions on which the periodical boundary conditions are applied. The density profiles and pressures were obtained taking average valúes in fluid slabs parallel to the walls. The pressure tensor was used in the simple virial form, as indicated in references [20, 21].

With T*=l, and for both slit pore widths L=10a and L=15<7, the size of the unit cell was set to L_y=L_z =30a, taking the number of particles N to correspond with the mean density. The range of the LennardJones interactions was considered with a cutoff radius of 5<7.

Figure 1: Density profiles for a n=9 layer model of a confined fluid in a slit pore (solid symbols). The temperature is T*=1.0 and the wall separation is L=10<7, with mean densities p*=l/20 (circles), 1/15 (squares), 1/10 (triangles), 1/5 (diamonds), and ^ (stars). Open symbols represent the Monte Cario simulations.

Figure 2: zz pressure tensor component. Captions as in Fig. 1.

IV. Results

In Figs. 1 and 2, the density profiles and zz components of the pressure tensor are shown for T*=l and L=10a. The mean densities studied are p*=l/20, 1/15, 1/10, 1/5, and 1/4. The agreement of the theoretical density profiles with the Monte Cario simulations is very good. For the pressure there is a rather good correspondence for low densities, up to p*=l/10. For the higher densities, differences appear, though the tendencies are similar. The discrepancies come first from the low density approximations done to get the Helmholtz energy. But, while in the simulation slab partióles fluctuate and at higher densities some clusterization occurs, in the theory each layer is supposed to have a homogeneous density which makes it hard for the theoretical pressures to follow those obtained by simulation. For the density profiles, averaging the number of partióles in each slab evidently compensates the clusterization, and the theory gives good results, at least for the rather low densities studied. The same picture applies to the behavior of the system for T*=l and L=15cr, at mean densities p* = 1/10, and 1/5, represented in Figs. 3 and 4.

 

Figure 3: Density profiles for a n=13 layer model of a confined fluid in a slit pore (solid symbols). Trie temperature is T*=1.0 and trie wall separation is L=15<7, with mean densities p*=l/10 (triangles), and 1/5 (circles). Open symbols represent the Monte Cario simulations.

The results, as expected for hard repulsive walls, show a low density región next to the walls and an increasing density profile, with a máximum at the center of the slit pore. This behavior is also shown with density functional theory [1] and in other Monte Cario simulations [2].

Finally, the good agreement of the theory with the experiment can be seen in the results shown in Fig. 5. In this figure, the excess number of molecules per unit área of pore surface T is plotted in function of the external pressure, at T*=3.18 and L=4cr. These parameters approximate the experimental valúes [19] T=303 K and L=1.45 nm, e/k = 95.2 K and a = 3.75 Á are used to characterize the nitrogen. In this case, due to the size of the sample, n=3 layers have been used for calculation. r is defined as

Figure 4: zz pressure tensor component. Captions as in Fig. 3.

 

where N_g/p* is the number/density of particles which would occupy the slit pore in the absence of the adsorption forces. p* and the external pressure are determined equating the chemical potential inside the slit pore (Eq. 10) to the chemical potential coming from the bulk van der Waals equation at the same temperature. The theoretical results presented here are similar to the numerical simulation results obtained by the same authors who have done the experiment [19]. They assume that the differences at higher pressures could be a consequence of the uncertainty in the determination of the pore geometry.

V. Conclusions

The application of a simple theory, with van der Waalslike approximations to the Helmholtz energy, to a particular model of spatial distribution makes it possible to obtain analytical expressions for the thermodynamic quantities. The study of a confined fluid in a slit pore geometry with a multilayer approximation produces good results when compared with Monte Cario simulations at low densities. The agreement with a particular experiment on nitrogen confined in a graphite slit pore is remarkable, even though an excess quantity is in study. It may be concluded that the confinement reduces the importance that higher virial contributions have on the equation of the state of the confined fluid. Classical density functional theory [22] can also be applied to study the slit pore geometry, with very good agreement with experiments and simulations. Though the theoretical work developed in these pages is not a competitor of density functional theory, it has the advantages of having analytical expressions, and the possibility of easily introducing two immiscible components: for instance one or two layer lubricants wetting the walls and a gas or a liquid filling the rest of layers forming the capillary volume.

 

Figure 5: Excess number of molecules per unit area of pore surface G as function of the external pressure. The full line represents the experiment (digitalized from Ref. [19]), and the dots are our theoretical results.

A correction has been made to get good critical parameters for the bulk (L - oo). For Argón a\ = 5.7538 and 6=1.3538, and for nitrogen a\= 1.5955 and 6=1.0349.

 

Acknowledgements This work was partially supported by Universidad Nacional de La Plata and CICPBA. G. J. Z. is member of “Carrera del Investigador Científico” CICPBA.

Appendix A

Expressions of quantities used in Eqs. 58:

 

 

[1] S A Sartarelli, L Szybisz, Correlation between asymmetric profiles in slits and standard prewetting lines, Pap. Phys. 1, 010001 (2009); L Szybisz, S A Sartarelli, Density profiles of Ar adsorbed in slits of CO₂: Spontaneous symmetry breaking revisited, J. Chem. Phys. 128, 124702 (2008).         [ Links ]

[2] M Schoen, Computer Simulation of Condensed Phases in Complex Geometries (Lecture Notes in Physics), Springer, Berlin (1993).         [ Links ]

[3] M Schoen, Structure and phase behavior of confined soft condensed matter, In: Computational Methods in Surface and Colloid Science, Ed. M Borowko, Pag. 1, Marcel Dekker, New York (2000).         [ Links ]

[4] S Dietrich, Fluids in contact with structured substrates, In: New Approaches to Problems in Liquid State Theory, Eds. C Caccamo, J P Hansen, G Stell, Pag. 197, Kluwer, Dordrecht (1999).         [ Links ]

[5] J P R B Walton, N. Quirke, Capillary Condensation: A Molecular Simulation Study, Mol. Simul. 2, 361 (1989).         [ Links ]

[6] L D Gelb, K E Gubbins, R Radhakrishnan, M SliwinskaBartkowiak, Phase separation in confined systems, Rep. Prog. Phys. 62, 1573 (1999).         [ Links ]

[7] A Maciolek, A Ciach, R Evans, Critical depletion of fluids in pores: Competing bulk and surface fields, J. Chem. Phys. 108, 9765 (1998).         [ Links ]

[8] A Maciolek, R Evans, N B Wilding, Effects of weak surface fields on the density profiles and adsorption of a confined fluid near bulk criticality, J. Chem. Phys. 119, 8663 (2003).         [ Links ]

[9] P B Balbuena, K E Gubbins, Classification of adsorption behavior: simple fluids in pores of slitshaped geometry, Fluid Phase Equilib. 76, 21 (1992).         [ Links ]

[10] P B Balbuena, K E Gubbins, Theoretical interpretation of adsorption behavior of simple fluids in slit pores, Langmuir 9, 1801 (1993).         [ Links ]

[11] O Pizio, A Patrykiejew, S Sokolowski, Phase Behavior of LennardJones Fluids in Slitlike Pores with Walls Modified by Preadsorbed Molecules: A Density Functional Approach, J. Phys. Chem. C 111, 15743 (2007).         [ Links ]

[12] G J Zarragoicoechea, V A Kuz, van der Waals equation of state for a fluid in a nanopore, Physical Review E 65, 021110 (2002).         [ Links ]

[13] G J Zarragoicoechea,V A Kuz, Critical shift of a confined fluid in a nanopore, Fluid Phase Equilib. 220, 7 (2004).         [ Links ]

[14] M Schoen, D J Diestler, Liquidvapor coexistence in a chemically heterogeneous slitnanopore, Chem. Phys. Letters 270, 339 (1997).         [ Links ]

[15] R Defay, I Prigogine, Surface Tension and Adsorption, Longmans, London (1966).

[16] R Defay, I Prigogine, Surface tension of regular solutions, Trans. Faraday Soc. 46, 199 (1950)).

[17] T Murakami, S Ono, M Tamura, M Kurata, On the theory of surface tension of regular solution, J. Phys. Soc. Japan 6, 309 (1951).

[18] L D Landau, E M Lifshitz, Física Estadística, Revert´e, Barcelona (1969), pp. 269273.

[19] K Kaneko, R F Cracknell and D Nicholson, Nitrogen Adsorption in Slit Pores at Ambient Temperatures: Comparison of Simulation and Experiment, Langmuir 10, 4606 (1994).

[20] M Schoen, D J Diestler, Analytical treatment of a simple fluid adsorbed in a slitpore, J. Chem. Phys. 109, 5596 (1998).

[21] M P Allen, D J Tildesley, Computer Simulation of Liquids, Oxford University Press, London (1987), Pag. 4647.

[22] J Wu, Density Functional Theory for chemical engineering: from capillarity to soft materials, AIChe J. 52, 1169 (2006).