Calculation of chemical potentials of chain molecules by the incremental gauge cell method

The gauge cell Monte Carlo method is extended to calculations of the incremental chemical potentials and free energies of linear chain molecules. The method was applied to chains of Lennard-Jones beads with stiff harmonic bonds up to 500 monomers in length. We show that the suggested method quantitatively reproduces the modified Widom particle insertion method of Kumar et al. [S. K. Kumar, I. Szleifer, and A. Z. Panagiotopoulos, Phys. Rev. Lett. 66(22), 2935 (1991)], and is by an order of magnitude more efficient for long chains in terms of the computational time required for the same accuracy of chemical potential calculations. The chain increment ansatz, which suggests that the incremental chemical potential is independent of the chain length, was tested at different temperatures. We confirmed that the ansatz holds only for coils above the θ temperature. Special attention is paid to the effects of the magnitude of adsorption potential and temperature on the behavior of single chains in confinements that are comparable in size with the free chain radius of gyration. At sufficiently low temperatures, the dependence of the incremental chemical potential on the chain length in wetting pores is superficially similar to a capillary condensation isotherm, reflecting monolayer formation following by pore volume filling, as the chain length increases. We find that the incremental gauge cell method is an accurate and efficient technique for calculations of the free energies of chain molecules in bulk systems and nanoconfinements alike. The suggested method may find practical applications, such as modeling polymer partitioning on porous substrates and dynamics of chain translocation into nanopores.     

A simulation method for the calculation of chemical potentials in small, inhomogeneous, and dense systems

We present a modification of the gauge cell Monte Carlo simulation method [A. V. Neimark and A. Vishnyakov, Phys. Rev. E 62, 4611 (2000)] designed for chemical potential calculations in small confined inhomogeneous systems. To measure the chemical potential, the system under study is set in chemical equilibrium with the gauge cell, which represents a finite volume reservoir of ideal particles. The system and the gauge cell are immersed into the thermal bath of a given temperature. The size of the gauge cell controls the level of density fluctuations in the system. The chemical potential is rigorously calculated from the equilibrium distribution of particles between the system cell and the gauge cell and does not depend on the gauge cell size. This scheme, which we call a mesoscopic canonical ensemble, bridges the gap between the canonical and the grand canonical ensembles, which are known to be inconsistent for small systems. The ideal gas gauge cell method is illustrated with Monte Carlo simulations of Lennard-Jones fluid confined to spherical pores of different sizes. Special attention is paid to the case of extreme confinement of several molecular diameters in cross section where the inconsistency between the canonical ensemble and the grand canonical ensemble is most pronounced. For sufficiently large systems, the chemical potential can be reliably determined from the mean density in the gauge cell as it was implied in the original gauge cell method. The method is applied to study the transition from supercritical adsorption to subcritical capillary condensation, which is observed in nanoporous materials as the pore size increases.     

Capillary condensation of short-chain molecules

A density-functional study of capillary condensation of fluids of short-chain molecules confined to slitlike pores is presented. The molecules are modeled as freely jointed tangent spherical segments with a hard core and with short-range attractive interaction between all the segments. We investigate how the critical parameters of capillary condensation of the fluid change when the pore width decreases and eventually becomes smaller than the nominal linear dimension of the single-chain molecule. We find that the dependence of critical parameters for a fluid of dimers and of tetramers on pore width is similar to that of the monomer fluid. On the other hand, for a fluid of chains consisting of a larger number of segments we observe an inversion effect. Namely, the critical temperature of capillary condensation decreases with increasing pore width for a certain interval of values of the pore width. This anomalous behavior is also influenced by the interaction between molecules and pore walls. We attribute this behavior to the effect of conformational changes of molecules upon confinement.     

The birth of a bubble: a molecular simulation study

We study the nucleation of a bubble in a metastable Lennard-Jones (LJ) fluid, confined to a spherical pore with wetting walls, by a combination of grand canonical, canonical ensemble, and gauge cell Monte Carlo simulation methods complemented by the Voronoi-Delaunay tessellation analysis of statistical geometry of intermolecular cavities. We construct the isotherm of confined fluid in the form of a continuous van der Waals' loop, in which the unstable backward trajectory between the spinodals corresponds to bubble states. We show that as the degree of metastability increases and the fluid becomes progressively stretched, the decrease of fluid density is associated with the evolution of a population of interstitial intermolecular cavities. At the spinodal, the fluid becomes mechanically unstable: Interstitial cavities partly coalesce into a larger cavity located due to the system symmetry around the pore center. This cavity represents a bubble embryo, which grows at the expense of interstitial cavities. The nucleation barrier is calculated by direct thermodynamic integration along the isotherm. We compare our simulation results to the predictions of the classical nucleation theory and experiments on capillary condensation-evaporation of nitrogen in pores of hybrid organic-inorganic mesoporous molecular sieve HMM-3.     

Monte Carlo simulation of cavitation in pores with nonwetting defects

We investigate the onset of cavitation in a metastable fluid confined to nanoscale pores with nonwetting defects present. Using grand canonical and gauge cell mesocanonical Monte Carlo simulations, we study the degree of metastability (relative vapor pressure), at which the critical bubble forms in a spherical pore with a circular nonwetting defect. It is shown that an increase of the defect size leads to a transition from homogeneous to heterogeneous nucleation of critical bubbles formed at the defect site. In this case, the desorption process may be initiated at larger relative vapor pressures than those predicted by the theories of homogeneous cavitation.     

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.     

Grand potential, helmholtz free energy, and entropy calculation in heterogeneous cylindrical pores by the grand canonical Monte Carlo simulation method

The adsorption of fluids in porous media is still an open area of research, since no model is able to explain all experimental features. The difficulties rise from the complexity of the real porous materials which present surface heterogeneities, large pore size distributions, and complex networks of interconnected pores. In parallel to experimental efforts trying to produce more ordered porous materials, theoreticians try to introduce more disorder in their models, with the help of molecular simulation for instance. This grand canonical Monte Carlo simulation study concentrates on the adsorption of a simple Lennard-Jones fluid in three porous substrates, to compare the effect of purely geometric heterogeneity (spatial deformation of the external potential) as opposed to purely chemical heterogeneity (amplitude variations of the external potential). This separation is unrealistic, since geometric fluctuations of a real pore diameter along its axis generally induce variations in the amplitude of the external potential created by the pore. However it enables one to compare both effects. In this paper, a thermodynamic integration scheme is applied to a complete set of adsorption/desorption isotherms. The grand potential, free energy, and entropy are calculated, which allows one to discuss the features of the phase diagrams. It is shown that a purely geometric deformation (undulation) of the external potential does not affect the thermodynamic characteristics of the confined fluid. On the other hand, amplitude modulation of the external potential (chemical heterogeneity) strongly distorts the phase diagram. This heterogeneity is actually able to stabilize a "bridgelike" phase which corresponds to an accumulation of molecules in the most attractive region of the pore.     

Capillary condensation in a geometrically and a chemically heterogeneous pore: a molecular simulation study

A computer simulation study has been carried out, using an extended Gibbs ensemble Monte Carlo technique, to examine the influence of so-called geometric and chemical disorder on the thermodynamic behavior of simple fluids confined in porous media. The technique allows the equilibrium coexistence of gas and liquid phases to be calculated in a single run. The phase diagram of Lennard-Jones fluid has been calculated in a perfectly cylindrical pore as a reference. Some disorder is then introduced in the porous material, first by spatially modifying the external potential of the initially cylindrical pore, to imitate the geometric disorder of a more realistic pore (undulation, constrictions, etc.) and second by modulating the amplitude of the same initially cylindrical potential to reproduce the energetic disorder of realistic pores due to chemical variations along it. It is shown that the chemical disorder has a much stronger effect on the phase diagram of the confined fluid. The complete adsorption/desorption isotherms are also calculated to help in understanding the large effects of chemical disorder.     

Structure and phase behavior of a confined nanodroplet composed of the flexible chain molecules

A polymer density functional theory has been employed for investigating the structure and phase behaviors of the chain polymer, which is modelled as the tangentially connected sphere chain with an attractive interaction, inside the nanosized pores. The excess free energy of the chain polymer has been approximated as the modified fundamental measure-theory for the hard spheres, the Wertheim's first-order perturbation for the chain connectivity, and the mean-field approximation for the van der Waals contribution. For the value of the chemical potential corresponding to a stable liquid phase in the bulk system and a metastable vapor phase, the flexible chain molecules undergo the liquid-vapor transition as the pore size is reduced; the vapor is the stable phase at small volume, whereas the liquid is the stable phase at large volume. The wide liquid-vapor coexistence curve, which explains the wide range of metastable liquid-vapor states, is observed at low temperature. The increase of temperature and decrease of pore size result in a narrowing of liquid-vapor coexistence curves. The increase of chain length leads to a shift of the liquid-vapor coexistence curve towards lower values of chemical potential. The coexistence curves for the confined phase diagram are contained within the corresponding bulk liquid-vapor coexistence curve. The equilibrium capillary phase transition occurs at a higher chemical potential than in the bulk phase.     

Density functional description of adsorption in slitlike pores modified with chain molecules: a simple model for pillaredlike materials

We propose a density functional theory to describe adsorption of Lennard-Jones fluid in slitlike pores modified by chain molecules. Specifically, the chains are bonded by their ends to the opposite pore walls, so they can form pillaredlike structure. Two models are studied. In the first model, the nonterminating segments of chains can change their configuration inside the pore upon adsorption of spherical species. In the second model, the chains configuration remains fixed, so that the system is similar to a nonuniform quenched-annealed mixture. We study capillary condensation of fluid species inside such modified pores and compare the results obtained for two models.     

Capillary condensation hysteresis in overlapping spherical pores: a monte carlo simulation study

The mechanisms of hysteretic phase transformations in fluids confined to porous bodies depend on the size and shape of pores, as well as their connectivity. We present a Monte Carlo simulation study of capillary condensation and evaporation cycles in the course of Lennard-Jones fluid adsorption in the system of overlapping spherical pores. This model system mimics pore shape and connectivity in some mesoporous materials obtained by templating cubic surfactant mesophases or colloidal crystals. We show different mechanisms of capillary hysteresis depending on the size of the window between the pores. For the system with a small window, the hysteresis cycle is similar to that in a single spherical pore: capillary condensation takes place upon achieving the limit of stability of adsorption film and evaporation is triggered by cavitation. When the window is large enough, the capillary condensation shifts to a pressure higher than that of the isolated pore, and the possibility for the equilibrium mechanism of desorption is revealed. These finding may have important implications for practical problems of assessment of the pore size distributions in mesoporous materials with cagelike pore networks.     

Anisotropic dynamics of dipolar liquids in narrow slit pores

We report molecular dynamics simulation results for Stockmayer fluids confined to narrow slitlike pores with structureless, nonconducting walls. The translational and rotational dynamics of the dipolar particles have been investigated by calculating autocorrelation functions, diffusion coefficients, and relaxation times for various pore widths (five or less particle diameters) and directions parallel and perpendicular to the walls. The dynamic properties of the confined systems are compared to bulk properties, where corresponding bulk and pore states at the same temperature and chemical potential are determined in parallel grand canonical Monte Carlo simulations. We find that the dynamic behavior inside the pore depends on the distance from the walls and can be strongly anisotropic even in globally isotropic systems. This concerns especially the particles in the surface layers close to the walls, where the single particle and collective dipolar relaxation resemble that of true two-dimensional dipolar fluids with different in-plane and out-of-plane relaxations. On the other hand, bulklike relaxation is observed in the pore center of sufficiently wide pores.     

Thermodynamic pressure of simple fluids confined in cylindrical nanopores by isothermal-isobaric Monte Carlo: influence of fluid/substrate interactions

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.     

Monte Carlo Calculations of Polymer Adsorption at the Entrance of Cylindrical Pores in Flat Adsorbing Surfaces

The influence of surface interactions on the conformation of flexible polymers partially confined inside narrow cylindrical pores in a flat surface is studied above the critical adsorption energy in a good solvent. We use a static configurational bias computational sampling method to calculate the adsorption free energy and the radius of gyration components parallel and perpendicular to the pore axis as a function of the polymer center of mass position at different degrees of confinement. We find strong free‐energy minima just in front of the pore entry for all degrees of confinement studied. At the location of the free‐energy minimum, polymers are partially adsorbed inside the pore and on the outer solid surface and adopt “drawing pin”‐like conformations. A distinct maximum in the average loop length at the pore entry indicates that the polymer bridges the pore entry of small pores.     

Behavior of acetonitrile confined to mesoporous silica gels as studied by 129Xe NMR: a novel method for determining the pore sizes

129Xe NMR spectra of xenon dissolved in acetonitrile confined into three mesoporous silica gels with nominal pore diameters of 40, 60, and 100 A have been measured over the temperature range 170-245 K. The spectra consist of a number of lines, which contain detailed information on the system. The most interesting result is that the chemical shift of a particular signal observed below the melting point of confined acetonitrile is highly sensitive to the pore size, and hence its shape is sensitive to the pore size distribution function. This signal originates from the xenon atoms sited in very small cavities built up inside the pores during the freezing transition. It can be used to determine the size or even the size distribution function of the pores. In addition, the emergence of this signal reveals the phase transition temperature of acetonitrile inside the pores, which can also be used to determine the size of the pores. The difference in the chemical shifts of two other signals, which arise from xenon dissolved in bulk and confined acetonitrile, provides still another novel method for determining the size of the pores.     

Capillary condensation in pores with rough walls: a density functional approach

The effect of surface roughness of slit-like pore walls on the capillary condensation of a spherical particles and short chains is studied. The gas molecules interact with the substrate by a Lennard-Jones (9,3) potential. The rough layer at each pore wall has a variable thickness and density and consists of a disordered quenched matrix of spherical particles. The system is described in the framework of a density functional approach and using computer simulations. The contribution due to attractive van der Waals interactions between adsorbate molecules is described by using first-order mean spherical approximation and mean-field approximation.     

Sol-Gel Preparation of Fast Proton-Conducting P2O5-SiO2 Glasses

We have investigated the proton conductivities of the sol-gel-derived P₂O₅-SiO₂ glass at –50 to 120°C. The obtained glass is porous, where the surface area, pore volume and pore diameter are 740 m²/g, 0.5 cm³/g and <5 nm, respectively. The freezing temperature of water molecules adsorbed in the pores was –20°C, which is much lower than that of free liquid water due to the quantum size effect of the water confined in the pores. The electrical conductivities followed the Arrhenius equation in the temperatures between –20 and 120°C. Below –20°C, the adsorbed-water molecules were frozen, resulting in a rapid decrease of the proton conductivity. Considering the high conductivity, chemical and thermal stability, this oxide glass membranes have potential for the fuel cell membrane.     

An alternative NMR method to determine nuclear shielding anisotropies for molecules in liquid-crystalline solutions with (13)C shielding anisotropy of methyl iodide as an example

An alternative NMR method for determining nuclear shielding anisotropies in molecules is proposed. The method is quite simple, linear and particularly applicable for heteronuclear spin systems. In the technique, molecules of interest are dissolved in a thermotropic liquid crystal (LC) which is confined in a mesoporous material, such as controlled pore glass (CPG) used in this study. CPG materials consist of roughly spherical particles with a randomly oriented and connected pore network inside. LC Merck Phase 4 was confined in the pores of average diameter from 81 to 375 A and LC Merck ZLI 1115 in the pores of average diameter 81 A. In order to demonstrate the functionality of the method, the (13)C shielding anisotropy of (13)C-enriched methyl iodide, (13)CH(3)I, was determined as a function of temperature using one dimensional (13)C NMR spectroscopy. Methane gas, (13)CH(4), was used as an internal chemical shift reference. It appeared that methyl iodide molecules experience on average an isotropic environment in LCs inside the smallest pores within the whole temperature range studied, ranging from bulk solid to isotropic phase. In contrast, in the spaces in between the particles, whose diameter is approximately 150 microm, LCs behave as in the bulk. Consequently, isotropic values of the shielding tensor can be determined from spectra arising from molecules inside the pores at exactly the same temperature as the anisotropic ones from molecules outside the pores. Thus, for the first time in the solution state, shielding anisotropies can easily be determined as a function of temperature. The effects of pore size as well as of different LC media on the shielding anisotropy are examined and discussed.     

Determination of pore sizes and volumes of porous materials by 129Xe NMR of xenon gas dissolved in a medium

In our previous paper (J. Phys. Chem. B 2005, 109, 757) it was illustrated that the 129Xe NMR spectra of xenon dissolved in acetonitrile confined into mesoporous materials give detailed information on the system, especially about the pore sizes. A resonance signal originating from xenon atoms sited in very small cavities built up inside the pores during the freezing transition (referred to as signal D) turned out to be highly sensitive to the pore size. The emergence of this signal reveals the phase transition temperature of acetonitrile inside the pores, which can also be used to determine the size of the pores. In addition, the difference in the chemical shifts of two other signals arising from xenon dissolved in bulk and confined acetonitrile (B and C) provides another method for determining the pore sizes. In the present work, the observed correlations have been investigated using an extensive set of measurements with a variety of porous materials (silica gels and controlled pore glasses) with the mean pore diameters ranging from 43 to 2917 A. The usefulness of the correlations has been demonstrated by calculating the pore size distributions from the spectral data. The distributions are in agreement with those reported by the manufacturers, when the mean pore diameter is smaller than approximately 500 A. In addition, it has been shown that the porosity of the materials can be determined by comparing the intensities of the signals arising from the bulk and confined liquid. When acetonitrile is replaced by cyclohexane in the sample, the dependence of the chemical shift difference between the B and C signals on the pore size becomes more sensitive, but no D signal appears below the freezing point. In addition, the influence of xenon gas on the melting points of bulk and confined acetonitrile has been studied by 1H NMR cryoporometry. The measurements show that the temperature of the latter transition lowers slightly more, and consequently affects the pore sizes calculated by means of the difference in the phase transition temperatures. Hysteresis in the phase transitions in a cooling-warming cycle has also been studied as a function of the temperature stabilization time by 129Xe NMR of xenon dissolved in acetonitrile.