Global phase diagrams for freezing in porous media
J. Chem. Phys. 116, 1147 (2002); doi:10.1063/1.1426412
Issue Date: 15 January 2002
You are not logged in to this journal. Log in
Using molecular simulations and free energy calculations based on Landau theory, we show that freezing/melting behavior of fluids of small molecules in pores of simple geometry can be understood in terms of two main parameters: the pore width H* (expressed as a multiple of the diameter of the fluid molecule) and a parameter
that measures the ratio of the fluid-wall to the fluidfluid attractive interaction. The value of the
parameter determines the qualitative nature of the freezing behavior, for example, the direction of change in the freezing temperature and the presence or absence of new phases. For slit-shaped pores, larger
values lead to an increase in the freezing temperature of the confined fluid, and to the presence of a hexatic phase. For pores that accommodate three or more layers of adsorbate molecules several kinds of contact layer phase (inhomogeneous phases in which the contact layer has a different structure than the inner layers) are observed. Smaller
values lead to a decrease in the freezing temperature. The parameter H* determines the magnitude of shift in the freezing temperature, and can also affect the presence of some of the new phases. Results are presented as plots of transition temperature vs
for a particular pore width. Experimental results are also presented for a variety of adsorbates in activated carbon fibers (ACF) covering a wide range of
values; the ACF have slit-shaped pores with average pore width 1.2 nm. The experimental and simulation results show qualitative agreement. ©2002 American Institute of Physics.
that measures the ratio of the fluid-wall to the fluidfluid attractive interaction. The value of the
parameter determines the qualitative nature of the freezing behavior, for example, the direction of change in the freezing temperature and the presence or absence of new phases. For slit-shaped pores, larger
values lead to an increase in the freezing temperature of the confined fluid, and to the presence of a hexatic phase. For pores that accommodate three or more layers of adsorbate molecules several kinds of contact layer phase (inhomogeneous phases in which the contact layer has a different structure than the inner layers) are observed. Smaller
values lead to a decrease in the freezing temperature. The parameter H* determines the magnitude of shift in the freezing temperature, and can also affect the presence of some of the new phases. Results are presented as plots of transition temperature vs
for a particular pore width. Experimental results are also presented for a variety of adsorbates in activated carbon fibers (ACF) covering a wide range of
values; the ACF have slit-shaped pores with average pore width 1.2 nm. The experimental and simulation results show qualitative agreement. ©2002 American Institute of Physics.
| History: | Received 23 August 2001; accepted 17 October 2001 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/116/1147/1 |
REFERENCES (45)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- L. D. Gelb, K. E. Gubbins, R. Radhakrishnan, and M. Sliwinska-Bartkowiak,
Rep. Prog. Phys. 62, 1573 (1999) . - M. Miyahara and K. E. Gubbins, J. Chem. Phys. 106, 2865 (1997).
- H. Dominguez, M. P. Allen, and R. Evans,
Mol. Phys. 96, 209 (1998) . - R. Radhakrishnan and K. E. Gubbins,
Mol. Phys. 96, 1249 (1999) . - M. A. Castro, S. M. Clarke, A. Inaba, and R. K. Thomas,
J. Phys. Chem. B 101, 8878 (1997) . - J. Klein and E. Kumacheva,
Science 269, 816 (1995) . - H. Hu, G. Carson, and S. Granick, Phys. Rev. Lett. 66, 2758 (1991).
- S. Granick,
Science 253, 1374 (1991) . - S. T. Cui, P. T. Cummings, and H. D. Cochran, J. Chem. Phys. 114, 7189 (2001).
- H. Christenson,
Colloids Surf., A 123, 355 (1997) . - K. Kaneko, A. Watanabe, T. Iiyama, R. Radhakrishnan, and K. E. Gubbins,
J. Phys. Chem. 103, 7061 (1999) . - A. Watanabe and K. Kaneko,
Chem. Phys. Lett. 305, 71 (1999) . - M. Sliwinska-Bartkowiak, G. Dudziak, R. Sikorski, R. Gras, K. E. Gubbins, and R. Radhakrishnan,
Phys. Chem. Chem. Phys. 3, 1179 (2001) . - R. Radhakrishnan, K. Gubbins, and M. Sliwinska-Bartkowiak, J. Chem. Phys. 112, 11048 (2000).
- M. Sliwinska-Bartkowiak, R. Radhakrishnan, and K. E. Gubbins,
Mol. Simul. 27, 323 (2001) . - R. Radhakrishnan, K. E. Gubbins, and M. Sliwinska-Bartkowiak, preprint (2001).
- E. Molz, A. P. Y. Wong, M. H. W. Chan, and J. R. Beamish, Phys. Rev. B 48, 5741 (1993).
- K. M. Unruh, T. E. Huber, and C. A. Huber, Phys. Rev. B 48, 9021 (1993).
- K. Overloop and L. V. Gerven,
J. Magn. Reson. 101, 179 (1993) . - K. Morishige and K. Kawano, J. Chem. Phys. 110, 4867 (1998).
- K. Morishige and K. Nabuoka, J. Chem. Phys. 107, 6965 (1997).
- H. F. Booth and J. H. Strange,
Mol. Phys. 93, 263 (1998) . - M. Sliwinska-Bartkowiak, J. Gras, R. Sikorski, R. Radhakrishnan, L. D. Gelb, and K. E. Gubbins,
Langmuir 15, 6060 (1999) . - M. Sliwinska-Bartkowiak, G. Dudziak, R. Sikorski, R. Gras, R. Radhakrishnan, and K. E. Gubbins, J. Chem. Phys. 114, 950 (2001).
- D. Morineau, G. Dosseh, C. Alba-Simionesco, and P. Llewellyn,
Philos. Mag. B 79, 1847 (1999) . - G. Dosseh, D. Morineau, and C. Alba-Simionesco, J. Phys. IV 10, Pr7-99 (2000).
- D. Morineau, F. Casas, C. Alba-Simionesco, A. Grosman, M.-C. Bellisent-Funel, and N. Ratovélomanana, J. Phys. IV 10, Pr7-95 (2000).
- W. A. Steele,
Surf. Sci. 36, 317 (1973) . - R. Radhakrishnan, K. Gubbins, A. Watanabe, and K. Kaneko, J. Chem. Phys. 111, 9058 (1999).
- R. M. Lynden-Bell, J. S. van Duijneveldt, and D. Frenkel,
Mol. Phys. 80, 801 (1993) . - P. M. Chaikin and T. C. Lubinski, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 1995).
- N. D. Mermin and H. Wagner,
Phys. Rev. Lett. 17, 1133 (1966) . - D. R. Nelson and B. I. Halperin, Phys. Rev. B 19, 2457 (1979).
- J. M. Kosterlitz and D. J. Thouless,
J. Phys. C 6, 1181 (1973) . - B. K. Peterson and K. E. Gubbins,
Mol. Phys. 62, 215 (1987) . - D. A. Kofke, J. Chem. Phys. 98, 4149 (1993).
- R. Agrawal and D. A. Kofke,
Mol. Phys. 85, 43 (1995) . - For strongly attractive pores, there exists a possibility of hexatic transitions in the layers other than the contact layers, although such transitions were not observed in our simulations. Since these are subtle phase transitions, they have a strong dependence on system size. Therefore their occurrence in the infinite system cannot be ruled out based on our calculations.
- For modeling the potential due to mica surfaces in the surface force apparatus (SFA), Cui et al. (Ref. 9) fitted the parameters of the 10-4-3 Steele potential [Eq. (2)], to reproduce the surface energy of 220 mJ/m2 reported in experiments (Ref. 8; see Table III). By using the 10-4-3 potential, Cui et al. (Ref. 9) assumed that the wall-potential due to mica surfaces is purely dispersive. However, as reported by Tsujimichi et al. (Ref. 40) in a atomic force microscopy (AFM) study of cleaved mica surfaces, there is a distribution of K+ ions in a background of SiO
and AlO
ions, so that induction forces may make a significant contribution.
- K. Tsujimichi, H. Tamura, A. Hirotani, M. Kubo, M. Komiyama, and A. Miyamoto,
J. Phys. Chem. B 101, 4260 (1997) . - J. O. Hirshfelder, C. F. Curtiss, and B. Bird, Molecular Theory of Gases and Liquids (Wiley, New York, 1964).
- W. L. Jorgensen, J. Chandrasekhar, and D. Madura, J. Chem. Phys. 79, 926 (1983).
- M. Maddox and K. E. Gubbins, J. Chem. Phys. 107, 9659 (1997).
- L. D. Gelb and K. E. Gubbins,
Langmuir 14, 2097 (1998) . - K. Thomson and K. Gubbins,
Langmuir 16, 5761 (2000) .








