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.     

Mean field kinetic theory for a lattice gas model of fluids confined in porous materials

We consider the mean field kinetic equations describing the relaxation dynamics of a lattice model of a fluid confined in a porous material. The dynamical theory embodied in these equations can be viewed as a mean field approximation to a Kawasaki dynamics Monte Carlo simulation of the system, as a theory of diffusion, or as a dynamical density functional theory. The solutions of the kinetic equations for long times coincide with the solutions of the static mean field equations for the inhomogeneous lattice gas. The approach is applied to a lattice gas model of a fluid confined in a finite length slit pore open at both ends and is in contact with the bulk fluid at a temperature where capillary condensation and hysteresis occur. The states emerging dynamically during irreversible changes in the chemical potential are compared with those obtained from the static mean field equations for states associated with a quasistatic progression up and down the adsorption/desorption isotherm. In the capillary transition region, the dynamics involves the appearance of undulates (adsorption) and liquid bridges (adsorption and desorption) which are unstable in the static mean field theory in the grand ensemble for the open pore but which are stable in the static mean field theory in the canonical ensemble for an infinite pore.     

Adsorption/desorption hysteresis in inkbottle pores: a density functional theory and Monte Carlo simulation study

The mechanisms of adsorption and desorption in inkbottle-shaped pores are considered for lattice models using grand canonical mean field density functional theory and Monte Carlo simulation. We find that they depend significantly on the particular pore geometry, the nature of the fluid-solid interaction, and the temperature. We find two mechanisms for desorption. One mechanism involves the emptying of the main cavity even as the density of fluid in the necks remains high, a mechanism observed recently in studies of an off-lattice model of an inkbottle. The other is a simultaneous desorption from the entire pore space, behavior that is more closely related to the traditional picture of pore blocking in the inkbottle system.     

Dynamic aspects of mercury porosimetry: a lattice model study

Grand canonical Monte Carlo simulations using both Glauber dynamics and Kawasaki dynamics have been carried out for a recently developed lattice model of a nonwetting fluid confined in a porous material. The calculations are aimed at investigating the molecular scale mechanisms leading to mercury retention encountered during mercury porosimetry experiments. We first describe a set of simulations on slit and ink-bottle pores. We have studied the influence of the pore width parameter on the intrusion/extrusion curve shapes and investigated the corresponding mechanisms. Entrapment appears during Kawasaki dynamics simulations of extrusion performed on ink-bottle pores when the system is studied for short relaxation times. We then consider the more realistic and complex case of a Vycor glass building on recent work on the dynamics of adsorption of wetting fluids (Woo, H. J.; Monson, P. A. Phys. Rev. E 2003, 67, 041207). Our results suggest that mercury entrapment is caused by a decrease in the rate of mass transfer associated with the fragmentation of the liquid during extrusion.     

Analysis of nucleation ability of cluster configurations with Monte Carlo simulations of argon

We determine the nucleation ability of argon clusters from Monte Carlo simulations. The nucleation rate appears to be defined by a sole characteristic of the clusters, namely, the stability. The stability is calculated as the ratio of grand canonical growth and decay rates and can be assigned to individual cluster configurations. We study the connection between the stability of the cluster configurations and their volume and total potential energy. Neither the potential energy nor the volume of a cluster configuration has a clear relation to its stability, and thus to the nucleation ability. On the other hand, we show that it is possible to use a specific volume for each cluster size to calculate the work of the cluster formation. These clusters with a unique volume have the same average stability as the full set of clusters. Our simulation method allows us to study the effect of possible deviations from equilibrium in the cluster configuration distributions. We argue that the nucleation process itself can produce a source for such a deviation. We show that even a small deviation from equilibrium in the cluster configuration distribution can lead to a dramatic deceleration of the nucleation rate. Although our simulations may overestimate the magnitude of the effect, they give qualitative estimates for its importance.     

Computer simulation of Lennard-Jones fluid adsorption in wedge-shaped pores with smooth walls

The adsorption isotherms and local density distributions in wedge-shaped pores with smooth solid walls are obtained for a molecular fluid by the grand canonical Monte Carlo and multiple histogram reweighting methods. The effect of wedge opening angle on the phase behavior of the fluid is considered, and the peculiarities of the mechanism of wedge-shaped pore filling are established.     

Monte Carlo simulations of hydrogen adsorption in alkali-doped single-walled carbon nanotubes

Monte Carlo simulations and Widom's test particle insertion method have been used to calculate the solubility coefficients (S) and the adsorption equilibrium constants (K) in single-walled (10,10) armchair carbon nanotubes including single nanotubes, and nanotube bundles with various configurations with and without alkali dopants. The hydrogen adsorption isotherms at room temperature were predicted by following the Langmuir adsorption model using the calculated constants S and K. The simulation results were in good agreement with experimental data as well as the grand canonical Monte Carlo simulation results reported in the literature. The simulations of nanotube bundle configurations suggest that the gravimetric hydrogen adsorption increases with internanotube gap size. It may be attributed to favorable hydrogen-nanotube interactions outside the nanotubes. The effect of alkali doping on hydrogen adsorption was studied by incorporating K+ or Li+ ions into nanotube arrays using a Monte Carlo simulation. The results on hydrogen adsorption isotherms indicate hydrogen adsorption of 3.95 wt% for K-doping, and 4.21 wt% for Li-doping, in reasonable agreement with the experimental results obtained at 100 atm and room temperature.     

Enabling grand‐canonical Monte Carlo: Extending the flexibility of GROMACS through the GromPy python interface module

We report on a python interface to the GROMACS molecular simulation package, GromPy (available at https://github.com/GromPy ). This application programming interface (API) uses the ctypes python module that allows function calls to shared libraries, for example, written in C. To the best of our knowledge, this is the first reported interface to the GROMACS library that uses direct library calls. GromPy can be used for extending the current GROMACS simulation and analysis modes. In this work, we demonstrate that the interface enables hybrid Monte‐Carlo/molecular dynamics (MD) simulations in the grand‐canonical ensemble, a simulation mode that is currently not implemented in GROMACS. For this application, the interplay between GromPy and GROMACS requires only minor modifications of the GROMACS source code, not affecting the operation, efficiency, and performance of the GROMACS applications. We validate the grand‐canonical application against MD in the canonical ensemble by comparison of equations of state. The results of the grand‐canonical simulations are in complete agreement with MD in the canonical ensemble. The python overhead of the grand‐canonical scheme is only minimal. © 2012 Wiley Periodicals, Inc.     

Local structure and thermodynamics of a core-softened potential fluid: theory and simulation.

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.     

Density functional theory of adsorption in pillared slit-like pores

We propose a density functional theory to describe adsorption of Lennard-Jones fluid in pillared slit like pores. Specifically, the pillars are built of chains that are bonded by their ends to the opposite pore walls. The approach we propose combines theory of quenched-annealed systems and theory of nonuniform fluids involving chain molecules. We compare the results of theoretical predictions with grand canonical ensemble Monte Carlo simulations and compute theoretical capillary condensation phase diagrams for several model systems.     

Phase equilibria and interfacial tension of fluids confined in narrow pores

Correlation between phase behaviors of a Lennard-Jones fluid in and outside a pore is examined over wide thermodynamic conditions by grand canonical Monte Carlo simulations. A pressure tensor component of the confined fluid, a variable controllable in simulation but usually uncontrollable in experiment, is related with the pressure of a bulk homogeneous system in equilibrium with the confined system. Effects of the pore dimensionality, size, and attractive potential on the correlations between thermodynamic properties of the confined and bulk systems are clarified. A fluid-wall interfacial tension defined as an excess grand potential is evaluated as a function of the pore size. It is found that the tension decreases linearly with the inverse of the pore diameter or width.     

Molecular simulations of aqueous electrolyte solubility: 1. The expanded-ensemble osmotic molecular dynamics method for the solution phase

We have developed a molecular-level simulation technique called the expanded-ensemble osmotic molecular dynamics (EEOMD) method, for studying electrolyte solution systems. The EEOMD method performs simulations at a fixed number of solvent molecules, pressure, temperature, and overall electrolyte chemical potential. The method combines elements of constant pressure-constant temperature molecular dynamics and expanded-ensemble grand canonical Monte Carlo. The simulated electrolyte solution systems contain, in addition to solvent molecules, full and fractional ions and undissociated electrolyte molecular units. The fractional particles are coupled to the system via a coupling parameter that varies between 0 (no interaction between the fractional particle and the other particles in the system) and 1 (full interaction between the fractional particle and the other particles in the system). The time evolution of the system is governed by the constant pressure-constant temperature equations of motion and accompanied by random changes in the coupling parameter. The coupling-parameter changes are accepted with a probability derived from the expanded-ensemble osmotic partition function corresponding to the prescribed electrolyte chemical potential. The coupling-parameter changes mimic insertion/deletion of particles as in a crude grand canonical Monte Carlo simulation; if the coupling parameter becomes 0, the fractional particles disappear from the system, and as the coupling parameter reaches unity, the fractional particles become full particles. The method is demonstrated for a model of NaCl in water at ambient conditions. To test our approach, we first determine the chemical potential of NaCl in water by the thermodynamic integration technique and by the expanded-ensemble method. Then, we carry out EEOMD simulations for different specified values of the overall NaCl chemical potential and measure the concentration of ions resulting from the simulations. Both computations give consistent results, validating the EEOMD methodology.     

Spatial updating grand canonical Monte Carlo algorithms for fluid simulation: generalization to continuous potentials and parallel implementation

Spatial updating grand canonical Monte Carlo algorithms are generalizations of random and sequential updating algorithms for lattice systems to continuum fluid models. The elementary steps, insertions or removals, are constructed by generating points in space either at random (random updating) or in a prescribed order (sequential updating). These algorithms have previously been developed only for systems of impenetrable spheres for which no particle overlap occurs. In this work, spatial updating grand canonical algorithms are generalized to continuous, soft-core potentials to account for overlapping configurations. Results on two- and three-dimensional Lennard-Jones fluids indicate that spatial updating grand canonical algorithms, both random and sequential, converge faster than standard grand canonical algorithms. Spatial algorithms based on sequential updating not only exhibit the fastest convergence but also are ideal for parallel implementation due to the absence of strict detailed balance and the nature of the updating that minimizes interprocessor communication. Parallel simulation results for three-dimensional Lennard-Jones fluids show a substantial reduction of simulation time for systems of moderate and large size. The efficiency improvement by parallel processing through domain decomposition is always in addition to the efficiency improvement by sequential updating.     

Accelerating molecular simulations by reversible mapping between local minima

A new framework is presented for performing Monte Carlo simulations of condensed matter based on a recently developed bijective mapping between local energy minima. The framework is used to implement a range of new multiparticle Monte Carlo moves, which are investigated by simulating atomic Lennard-Jones fluids in the canonical and grand canonical ensembles. Important aspects of the simulation protocol and their effect on performance are analyzed in detail. Using the mapping accelerates the simulations by many orders of magnitude when compared to the equivalent moves without the mapping, and leads to particularly efficient configurational sampling at low temperatures and high densities. The method appears to be suitable for adapting to quantitative simulations of more complex molecular systems over long effective time scales.     

Solvation force induced by short range, exact dissipative particle dynamics effective surfaces on a simple fluid and on polymer brushes

The thermodynamic properties of a simple fluid confined by effective wall forces are calculated using Monte Carlo simulations in the grand canonical ensemble. The solvation force produced by polymer brushes of two different lengths is obtained also. For the particular type of model interactions used, known as the dissipative particle dynamics method, we find that it is possible to obtain an exact, simple expression for the effective force induced by a planar wall composed of identical particles that interact with those in the fluid. We show that despite the short range of all forces in the model, the solvation force can be finite at relatively large distances and therefore does not depend only on the range of the interparticle or solvent-surface forces. As for the polymer brushes, we find that the shape of the solvation force profiles is in fair agreement with scaling and self-consistent field theories. The applications and possible extensions of this work are discussed.     

Chemical potential of a hard sphere fluid adsorbed in model disordered polydisperse matrices

We consider a model for adsorption of a simple fluid in disordered polydisperse adsorbents. The fluid consists of hard sphere particles. On the other hand, the adsorbents of this study are modeled as a collection of hard spheres with their diameter obeying a certain distribution function. Our focus is in the evaluation of the chemical potential of the fluid immersed in such a polydisperse material. It permits us to obtain porosity and pore size distribution for the adsorbent, as well as a set of adsorption isotherms. The latter have been calculated theoretically and by grand canonical Monte Carlo simulations. We observe that the width of assumed polydispersity distribution affects all the properties of the system. Nevertheless, the effect of matrix packing is dominant in determining adsorption for this class of models. We are convinced that the matrix structures generated via more sophisticated algorithms would exhibit stronger effects of polydispersity on the entire set of properties of adsorbed simple fluids.     

Boltzmann bias grand canonical Monte Carlo

We derive an efficient method for the insertion of structured particles in grand canonical Monte Carlo simulations of adsorption in very confining geometries. We extend this method to path integral simulations and use it to calculate the isotherm of adsorption of hydrogen isotopes in narrow carbon nanotubes (two-dimensional confinement) and slit pores (one-dimensional confinement) at the temperatures of 20 and 77 K, discussing its efficiency by comparison to the standard path integral grand canonical Monte Carlo algorithm. We use this algorithm to perform multicomponent simulations in order to calculate the hydrogen isotope selectivity for adsorption in narrow carbon nanotubes and slit pores at finite pressures. The algorithm described here can be applied to the study of adsorption of real oligomers and polymers in narrow pores and channels.     

Water condensation in hydrophobic silicalite-1 zeolite: a molecular simulation study

We report grand canonical Monte Carlo simulations of the gas and liquid phase adsorption of water in silicalite-1 zeolite. Simple but effective models and simulation methods, found useful for studying gas adsorption in nanoporous materials, have been extended to describe the intrusion/extrusion cycle of water in this hydrophobic solid. The picture of water confined to hydrophobic spaces of nanoscopic dimensions that emerges from this study is one of a strongly depleted and highly inhomogeneous fluid.