1
|
Puibasset J, Judeinstein P, Zanotti JM. Bulk supercooled water versus adsorbed films on silica surfaces: specific heat by Monte Carlo simulation. Phys Chem Chem Phys 2021; 23:2275-2285. [PMID: 33443254 DOI: 10.1039/d0cp05387b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Between 150 and 230.6 K, bulk supercooled water freezes upon cooling, and amorphous ice crystallizes upon heating: bulk water thus exists only in its stable ice form. To circumvent this problem, experiments are generally performed on water adsorbed in SiO2 based porous systems. In this work, we take advantage of Monte Carlo simulations to explore this metastable supercooled region inaccessible to experiments. Using three rigid, non-polarizable water models, namely SPC, TIP4P and TIP4P/2005, we investigate the isobaric specific heat capacity (Cp), between 100 and 300 K, of bulk water and water films of few monolayers adsorbed on different SiO2 surfaces: a smooth surface, a non-hydroxylated (0001) surface of quartz, and a fully hydroxylated (001) surface of cristobalite. As Cp is directly related to the entropy fluctuations and we focus on low temperatures, the convergence of the Monte Carlo simulations is a critical point of this work. Also, due to the small mass of the hydrogen atoms, quantum corrections are taken into account, and lead to an excellent agreement of the simulated and experimental Cp values at low temperature (100 K region). Altogether, we conclude that, in bulk, Cp is shown to exhibit a broad peak around 225 K for the SPC and TIP4P models, and around 250 K for the TIP4P/2005 model, in qualitative agreement with the experimentally observed features in Cp measurements. For interfacial water, in all cases, the broad Cp peak disappears. This result, at odds with experimental observations, suggests that disorder and hydrogen bonding at the interface (not yet taken into account) have a fundamental role in confined water transitions.
Collapse
Affiliation(s)
- J Puibasset
- ICMN, UMR7374, CNRS, Université d'Orléans, 1b, Rue de la Férollerie, 45071, Orléans Cedex 2, France.
| | | | | |
Collapse
|
2
|
Puibasset J, Judeinstein P, Zanotti JM. Molecular simulation study of the heat capacity of metastable water between 100 and 300 K. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1535179] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- J. Puibasset
- ICMN, UMR7374, CNRS, Université d’Orléans, Orléans, France
| | - P. Judeinstein
- Laboratoire Léon Brillouin, CEA, CNRS, Université Paris-Saclay, Gif-sur-Yvette Cedex, France
| | - J.-M. Zanotti
- Laboratoire Léon Brillouin, CEA, CNRS, Université Paris-Saclay, Gif-sur-Yvette Cedex, France
| |
Collapse
|
3
|
Terrón-Mejía KA, López-Rendón R, Gama Goicochea A. Mesoscopic modeling of structural and thermodynamic properties of fluids confined by rough surfaces. Phys Chem Chem Phys 2015; 17:26403-16. [DOI: 10.1039/c5cp03823e] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Simulations show that the ordering of particles confined by rough surfaces induces a structural phase transition while the interfacial tension is insensitive to it.
Collapse
Affiliation(s)
- Ketzasmin A. Terrón-Mejía
- Laboratorio de Bioingeniería Molecular a Multiescala
- Facultad de Ciencias
- Universidad Autónoma del Estado de México
- Toluca, Mexico
| | - Roberto López-Rendón
- Laboratorio de Bioingeniería Molecular a Multiescala
- Facultad de Ciencias
- Universidad Autónoma del Estado de México
- Toluca, Mexico
| | | |
Collapse
|
4
|
Puibasset J, Kierlik E, Tarjus G. Influence of system size on the properties of a fluid adsorbed in a nanopore: Physical manifestations and methodological consequences. J Chem Phys 2014; 141:044716. [PMID: 25084946 DOI: 10.1063/1.4891359] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Hysteresis and discontinuities in the isotherms of a fluid adsorbed in a nanopore in general hamper the determination of equilibrium thermodynamic properties, even in computer simulations. A way around this has been to consider both a reservoir of small size and a pore of small extent in order to restrict the fluctuations of density and approach a classical van der Waals loop. We assess this suggestion by thoroughly studying through Monte Carlo simulations and density functional theory the influence of system size on the equilibrium configurations of the adsorbed fluid and on the resulting isotherms. We stress the importance of pore-symmetry-breaking states that even for modest pore sizes lead to discontinuous isotherms and we discuss the physical relevance of these states and the methodological consequences for computing thermodynamic quantities.
Collapse
Affiliation(s)
- Joël Puibasset
- CRMD, CNRS FRE 3520, Université d'Orléans, 1b rue de la Férollerie, 45071 Orléans Cedex 02, France
| | - Edouard Kierlik
- LPTMC, CNRS UMR 7600, Université Pierre et Marie Curie, boîte 121, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Gilles Tarjus
- LPTMC, CNRS UMR 7600, Université Pierre et Marie Curie, boîte 121, 4 place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
5
|
Puibasset J. Fluid adsorption in linear pores: a molecular simulation study of the influence of heterogeneities on the hysteresis loop and the distribution of metastable states. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2013.829221] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
6
|
Singh SK, Khan S, Jana S, Singh JK. Vapour–liquid phase equilibria of simple fluids confined in patterned slit pores. MOLECULAR SIMULATION 2010. [DOI: 10.1080/08927022.2010.514778] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Sudhir K. Singh
- a Department of Chemical Engineering , Indian Institute of Technology , Kanpur, 208016, India
| | - Sandip Khan
- a Department of Chemical Engineering , Indian Institute of Technology , Kanpur, 208016, India
| | - Subimal Jana
- a Department of Chemical Engineering , Indian Institute of Technology , Kanpur, 208016, India
| | - Jayant K. Singh
- a Department of Chemical Engineering , Indian Institute of Technology , Kanpur, 208016, India
| |
Collapse
|
7
|
Puibasset J. Counting metastable states within the adsorption/desorption hysteresis loop: A molecular simulation study of confinement in heterogeneous pores. J Chem Phys 2010; 133:104701. [DOI: 10.1063/1.3483790] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Joël Puibasset
- Centre de Recherche sur la Matière Divisée, CNRS-Université d'Orléans, 1b rue de la Férollerie, Orléans Cedex 02 45071, France.
| |
Collapse
|
8
|
Puibasset J, Kierlik E, Tarjus G. Influence of reservoir size on the adsorption path in an ideal pore. J Chem Phys 2009; 131:124123. [DOI: 10.1063/1.3236510] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Bhattacharya S, Coasne B, Hung FR, Gubbins KE. Modeling micelle-templated mesoporous material SBA-15: atomistic model and gas adsorption studies. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2009; 25:5802-5813. [PMID: 19099416 DOI: 10.1021/la801560e] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
We report the development of a realistic molecular model for mesoporous silica SBA-15, which includes both the large cylindrical mesopores and the smaller micropores in the pore walls. The methodology for modeling the SBA-15 structure involves molecular and mesoscale simulations combined with geometrical interpolation techniques. First, a mesoscale model is prepared by mimicking the synthesis process using lattice Monte Carlo simulations. The main physical features of this mesoscale pore model are then carved out of an atomistic silica block; both the mesopores and the micropores are incorporated from the mimetic simulations. The calculated pore size distribution, surface area, and simulated TEM images of the model structure are in good agreement with those obtained from experimental samples of SBA-15. We then investigate the adsorption of argon in this structure using Grand Canonical Monte Carlo (GCMC) simulations. The adsorption results for our SBA-15 model are compared with those for a similar model that does not include the micropores; we also compare with results obtained in a regular cylindrical pore. The simulated adsorption isotherm for the SBA-15 model shows semiquantitative agreement with the experimental isotherm for a SBA-15 sample having a similar pore size. We observe that the presence of the micropores leads to increased adsorption at low pressure compared to the case of a model without micropores in the pore walls. At higher pressures, for all models, the filling proceeds via the monolayer-multilayer adsorption on the mesopore surface followed by capillary condensation, which is mainly controlled by the mesopore diameter and is not influenced by the presence of the micropores.
Collapse
Affiliation(s)
- Supriyo Bhattacharya
- Center for High Performance Computing and Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695-7905, USA
| | | | | | | |
Collapse
|
10
|
Kierlik E, Puibasset J, Tarjus G. Effect of the reservoir size on gas adsorption in inhomogeneous porous media. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:155102. [PMID: 21825355 DOI: 10.1088/0953-8984/21/15/155102] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
We study the influence of the relative size of the reservoir on the adsorption isotherms of a fluid in disordered or inhomogeneous mesoporous solids. We consider both an atomistic model of a fluid in a simple, yet structured pore, whose adsorption isotherms are computed by molecular simulation, and a coarse-grained model for adsorption in a disordered mesoporous material, studied by a density functional approach in a local mean-field approximation. In both cases, the fluid inside the porous solid exchanges matter with a reservoir of gas that is at the same temperature and chemical potential and whose relative size can be varied, and the control parameter is the total number of molecules present in the porous sample and in the reservoir. Varying the relative sizes of the reservoir and the sample within experimental range may change the shape of the hysteretic isotherms, leading to a 're-entrant' behavior compared to the grand-canonical isotherm when the latter displays a jump in density. We relate these phenomena to the organization of the metastable states that are accessible for the adsorbed fluid at a given chemical potential or density.
Collapse
Affiliation(s)
- E Kierlik
- Laboratoire de Physique Théorique de la Matière Condensée, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France
| | | | | |
Collapse
|
11
|
Puibasset J. Monte-Carlo multiscale simulation study of argon adsorption/desorption hysteresis in mesoporous heterogeneous tubular pores like MCM-41 or oxidized porous silicon. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2009; 25:903-911. [PMID: 19063620 DOI: 10.1021/la802474c] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
In a recent paper [J. Chem. Phys. 2007, 127, 154701] a multiscale approach was introduced which allowed calculation of adsorption/desorption hysteresis for fluid confined in a single mesoporous, heterogeneous tubular pore. The main interest in using such an approach is that it allows one to reconcile a molecular simulation approach generally limited to the nanometer scale (atomistic description of the confined fluid and pore roughness) with the much larger scale (micrometer) relevant to understand the complexity of adsorption/desorption hysteresis (the numerous metastable states in the hysteresis loop are a consequence of the large-scale disorder in the porous material). In this paper, this multiscale approach is used to study adsorption phenomena in mesoporous models made of a collection of disordered, noninterconnected tubular pores, as MCM-41 or porous silicon. A double distribution is introduced: one to characterize the disorder in a given pore, and the other to characterize the disorder between the pores. We consider two distribution shapes: Gaussian and uniform truncated and two cases of pores open at one or both ends. These models are expected to cover a wide variety of real materials made of independent pores, as MCM-41 and oxidized porous silicon. A large variety of hysteresis shapes is obtained, ranging from almost parallel adsorption/desorption branches typical of MCM-41 adsorption to triangular hysteresis typical of porous silicon. The structure of the metastable states inside the hysteresis (scanning adsorption/desorption curves) is also examined. The results are expected to be useful to experimentalists who want to infer pore structure and level of disorder from experimental adsorption/desorption experiments.
Collapse
Affiliation(s)
- Joël Puibasset
- Centre de Recherche sur la Matiere Divisee, CNRS-Universite d'Orleans, 1b rue de la Ferollerie, 45071 Orleans, Cedex 02, France.
| |
Collapse
|
12
|
Puibasset J. Pseudocritical or hysteresis temperature versus pore size for simple fluids confined in cylindrical nanopores. J Chem Phys 2008; 129:024705. [DOI: 10.1063/1.2948967] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Feng Z, Zhang X, Wang W. Adsorption of fluids in a pore with chemical heterogeneities: the cooperative effect. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:051603. [PMID: 18643075 DOI: 10.1103/physreve.77.051603] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2008] [Indexed: 05/26/2023]
Abstract
In this work, we study the cooperative adsorption of fluids in a heterogeneous pore, in which the pore walls are composed of homogeneous substrates with chemical groups (CGs) decorating them. The adsorption caused by the homogeneous substrates alone and that by CGs do not add up to the overall adsorption, indicating the existence of a cooperative effect. The cooperative effect is the source of cooperative adsorption, and is characterized in this work by the ratio of the overall adsorption to the sum of adsorption by the substrate only and that by CGs. It is found that the cooperative adsorption does not depend monotonically on the substrate or the CGs. Two different origins of the cooperative adsorption play different roles depending on which one dominates the overall adsorption. Our simulations reveal that, when the homogeneous substrate dominates the overall adsorption, weakening of the attractive fluid-substrate interaction or alternatively strengthening of the fluid-CGs interaction leads to a stronger cooperative effect and enhances the cooperative adsorption. However, when CGs dominate the overall adsorption, weakening of the attractive fluid-CG interaction or strengthening the fluid-substrate interaction results in strong cooperative adsorption. In order to investigate the effects of the distribution of CGs on cooperative adsorption, a design-test method is generalized and used in this work. Simulation results show that the overall adsorption can be significantly affected by the CG distribution.
Collapse
Affiliation(s)
- Zhikuan Feng
- Division of Molecular and Materials Simulation, Key Laboratory for Nanomaterials, Ministry of Education, Beijing University of Chemical Technology, Beijing 100029, China
| | | | | |
Collapse
|
14
|
Pellicane G, Vink RLC, Caccamo C, Löwen H. Colloid-polymer mixtures in the presence of quenched disorder: a theoretical and computer simulation study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2008; 20:115101. [PMID: 21694215 DOI: 10.1088/0953-8984/20/11/115101] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We use theory and computer simulation to study the structure and phase behavior of colloid-polymer mixtures in the presence of quenched disorder. The Asakura-Oosawa model (AO) (Asakura and Oosawa 1954 J. Chem. Phys. 22 1255) is used to describe the colloid-colloid, colloid-polymer, and polymer-polymer pair interactions. We then investigate the behavior of this model in the presence of frozen-in (quenched) obstacles. The obstacles will be placed according to two different scenarios, both of which are experimentally feasible. In the first scenario, polymers are distributed at positions drawn from an ideal gas configuration. In the second scenario, colloidal particles are distributed at positions drawn from an equilibrium hard sphere configuration. We investigate how the unmixing transition of the AO model is affected by the type of quenched disorder. The theoretical formalism is based on the replica method of Given and Stell (1994 Physica A 209 495). Our foremost aim is to test the accuracy of three common closures to the replica Ornstein-Zernike equations, namely the hypernetted chain, the Percus-Yevick, and the Martinov-Sarkisov equations. The accuracy is determined by comparison with grand canonical Monte Carlo simulations. We find that, for quenched polymer disorder, all three closures perform remarkably well. However, when quenched colloid disorder is considered, i.e. the second mentioned scenario, the predictions of all three closures worsen dramatically.
Collapse
Affiliation(s)
- Giuseppe Pellicane
- Dipartimento di Fisica Contrada Papardo, Università degli Studi di Messina, 98166 Messina, Italy
| | | | | | | |
Collapse
|
15
|
Puibasset J. Adsorption∕desorption hysteresis of simple fluids confined in realistic heterogeneous silica mesopores of micrometric length: A new analysis exploiting a multiscale Monte Carlo approach. J Chem Phys 2007; 127:154701. [DOI: 10.1063/1.2790423] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
16
|
Hung FR, Bhattacharya S, Coasne B, Thommes M, Gubbins KE. Argon and krypton adsorption on templated mesoporous silicas: molecular simulation and experiment. ADSORPTION 2007. [DOI: 10.1007/s10450-007-9034-4] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
17
|
Puibasset J. Thermodynamic pressure of simple fluids confined in cylindrical nanopores by isothermal-isobaric Monte Carlo: influence of fluid/substrate interactions. J Chem Phys 2007; 127:074702. [PMID: 17718622 DOI: 10.1063/1.2764477] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The thermodynamic pressure or grand potential density is calculated by isobaric-isothermal Monte Carlo algorithm for simple Lennard-Jones fluid confined in cylindrical pores presenting chemical heterogeneities along their axis. Heuristic arguments and simulation results show that the thermodynamic pressure of the confined fluid contains two contributions. The first term is the usual pressure of the bulk fluid for a density equal to the confined fluid density defined as the total number of confined particles divided by the accessible volume due to thermal agitation. A second term has to be added, which is empirically shown to be proportional to the fluid/wall interface area and almost constant along the adsorption and desorption branches. This interfacial contribution, calculated for various pore models, has small variations reminiscent of the fluid adsorption/desorption properties calculated in the various pores. In particular, it is shown that this interfacial quantity is maximum for a fluid/substrate interaction intensity of the same order as the fluid/fluid one, while the thermodynamic pressure at which rapid desorption occurs presents a minimum. Stronger or weaker fluid/wall affinity favors gas state nucleation on the desorption of confined fluids.
Collapse
Affiliation(s)
- Joël Puibasset
- Centre de Recherche sur la Matière Divisée, CNRS-Université d'Orléans, 1b rue de la Férollerie, 45071 Orléans, Cedex 02, France.
| |
Collapse
|
18
|
Puibasset J. Surface excess free energy of simple fluids confined in cylindrical pores by isothermal-isobaric Monte Carlo: influence of pore size. J Chem Phys 2007; 126:184701. [PMID: 17508818 DOI: 10.1063/1.2735582] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Confined fluid properties are mainly determined by interfacial phenomena characterized by surface quantities. Based on a simple model of Lennard-Jones particles confined in a cylindrical pore, this study introduces a grand potential surface quantity to quantify the difference in the thermodynamic pressure between the bulk and the confined fluids. The usual surface tension gamma defined as this grand potential difference for the same chemical potential in both confined and bulk states is generally strongly dependent on both the chemical potential and temperature. It is proposed here to introduce another surface quantity zeta which measures the thermodynamic pressure difference between confined and bulk states for identical densities. It is shown that this quantity is much less dependent on confined fluid density or chemical potential. It is actually constant along the gas-like and liquid-like adsorption/desorption branches for an irreversible isotherm (hysteresis), with a different value for each branch. For reversible supercritical isotherms, zeta is shown to remain constant in the low and high density parts of the isotherm. This independence on chemical potential (or equivalently fluid density) is believed to be of great interest for practical applications when one desires to calculate thermodynamic quantities such as the usual surface tension gamma or the thermodynamic pressure of a confined fluid for any given chemical potential and temperature. Such calculations are required to determine fundamental properties such as metastability or coexistence. The effects of temperature, fluid/substrate interaction strength, and pore size are studied.
Collapse
Affiliation(s)
- Joël Puibasset
- Centre de Recherche sur la Matière Divisée, CNRS-Université d'Orléans, 1b rue de la Férollerie, 45071 Orléans, Cedex 02, France.
| |
Collapse
|
19
|
Hung FR, Coasne B, Gubbins KE, Siperstein FR, Thommes M, Sliwinska-Bartkowiak M. A Monte Carlo study of capillary condensation of krypton within realistic models of templated mesoporous silica materials. ACTA ACUST UNITED AC 2007. [DOI: 10.1016/s0167-2991(07)80021-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
|
20
|
Puibasset J. Generalized isobaric–isothermal ensemble: application to capillary condensation and cavitation in heterogeneous nanopores. Mol Phys 2006. [DOI: 10.1080/00268970600938485] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
21
|
Puibasset J. Influence of surface chemical heterogeneities on adsorption/desorption hysteresis and coexistence diagram of metastable states within cylindrical pores. J Chem Phys 2006; 125:074707. [PMID: 16942364 DOI: 10.1063/1.2229193] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Grand canonical Monte Carlo simulations are performed to determine the adsorption/desorption isotherms at different temperatures of a Lennard-Jones fluid confined within a simple model of cylindrical pores presenting chemical heterogeneities. A complex hysteresis loop is observed, showing hysteresis subloops (scanning curves). This is shown to be consistent with the existence of several metastable states (local minima in the system free energy). A recent extension to the Gibbs ensemble technique is then used to calculate the complete coexistence diagram of these local minima.
Collapse
Affiliation(s)
- Joël Puibasset
- Centre de Recherche sur la Matière Divisée, CNRS, Université d'Orléans, 1b rue de la Férollerie, 45071 Orléans, Cedex 02, France.
| |
Collapse
|
22
|
Jiang J, Sandler SI. Capillary phase transitions of linear and branched alkanes in carbon nanotubes from molecular simulation. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2006; 22:7391-9. [PMID: 16893243 DOI: 10.1021/la0608720] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Capillary phase transitions of linear (from C(1) to C(12)) and branched (C(5) isomers) alkanes in single-walled carbon nanotubes have been investigated using the gauge-cell Monte Carlo simulation. The isotherm at a supercritical temperature increases monotonically with chemical potential and coincides with that from the traditional grand canonical Monte Carlo simulation, whereas the isotherm at a subcritical temperature exhibits a sigmoid van der Waals loop including stable, metastable, and unstable regions. Along this loop, the coexisting phases are determined using an Maxwell equal-area construction. A generic confinement effect is found that reduces the saturation chemical potential, lowers the critical temperature, increases the critical density, and shrinks the phase envelope. The effect is greater in a smaller diameter nanotube and is greater in a nanotube than in a nanoslit.
Collapse
Affiliation(s)
- Jianwen Jiang
- Department of Chemical and Biomolecular Engineering, National University of Singapore, 4 Engineering Drive 4, Singapore 117576.
| | | |
Collapse
|