Effects of confinement on the molar enthalpy of argon adsorption in graphitic cylindrical pores: a grand canonical Monte Carlo (GCMC) simulation study

Using a grand canonical Monte Carlo simulation, we study argon adsorption in graphitic cylindrical pores to investigate the differences between the isosteric heat and the integral molar enthalpy under subcritical and supercritical conditions and compare these results against those for a flat graphite surface to investigate the role of confinement on the enthalpy change of adsorption. The isosteric heat curve is finite under subcritical conditions, but for supercritical adsorption, it becomes infinite at the pressure where the excess concentration versus pressure is maximum. This can be circumvented using the integral molar enthalpy, which is a better variable to describe the energy change for supercritical adsorption. Finally, the effects of pore geometry (radius and length) on argon adsorption under subcritical and supercritical conditions are discussed.     

Adsorption of argon from sub- to supercritical conditions on graphitized thermal carbon black and in graphitic slit pores: a grand canonical Monte Carlo simulation study

In this paper we consider the adsorption of argon on the surface of graphitized thermal carbon black and in slit pores at temperatures ranging from subcritical to supercritical conditions by the method of grand canonical Monte Carlo simulation. Attention is paid to the variation of the adsorbed density when the temperature crosses the critical point. The behavior of the adsorbed density versus pressure (bulk density) shows interesting behavior at temperatures in the vicinity of and those above the critical point and also at extremely high pressures. Isotherms at temperatures greater than the critical temperature exhibit a clear maximum, and near the critical temperature this maximum is a very sharp spike. Under the supercritical conditions and very high pressure the excess of adsorbed density decreases towards zero value for a graphite surface, while for slit pores negative excess density is possible at extremely high pressures. For imperfect pores (defined as pores that cannot accommodate an integral number of parallel layers under moderate conditions) the pressure at which the excess pore density becomes negative is less than that for perfect pores, and this is due to the packing effect in those imperfect pores. However, at extremely high pressure molecules can be packed in parallel layers once chemical potential is great enough to overcome the repulsions among adsorbed molecules.     

High-pressure adsorption capacity and structure of CO2 in carbon slit pores: theory and simulation

We present new simulation results for the packing of single-center and three-center models of carbon dioxide at high pressure in carbon slit pores. The former shows a series of packing transitions that are well described by our density functional theory model developed earlier. In contrast, these transitions are absent for the three-center model. Analysis of the simulation results shows that alternations of flat-lying molecules and rotated molecules can occur as the pore width is increased. The presence or absence of quadrupoles has negligible effect on these high-density structures.     

Monte Carlo modeling of chiral adsorption on nanostructured chiral surfaces and slit pores

Adsorptive separation of chiral molecules is a powerful technique that has long been used in the chemical and pharmaceutical industries. An important challenge in this field is to design and optimize new adsorbents to provide selective discrimination of enantiomers. In this article, we introduce an off-lattice model of chiral adsorption on nanostructured surfaces and slit pores with the aim of predicting their enantioslective properties. The concept presented here involves finding the optimal chiral pattern of active sites on the pore walls that maximizes the difference between the binding energies of the enantiomers. Our initial effort focuses on chiral molecules that do not have specific interactions with the pore surface. One candidate meeting this requirement is 1,2-dimethylcyclopropane (DMCP), a chiral hydrocarbon whose interaction with a model pore surface was described using the Lennard-Jones potential. To model the adsorption of DMCP, we used the Monte Carlo simulation method. It was demonstrated that the separation of the enantiomers of DMCP is hardly obtainable because of the smoothness of the potential energy surface for molecules physisorbed in the pore. However, the simulated results allowed the identification of key factors that influence the binding of the enantiomers of DMCP to the pore walls with a special distribution of active sites. This information will be useful in future considerations of the adsorption of more complex chiral molecules.     

Polymer nanodroplets forming liquid bridges in chemically structured slit pores: a computer simulation

Using a coarse-grained bead-spring model of flexible polymer chains, the structure of a polymeric nanodroplet adsorbed on a chemically decorated flat wall is investigated by means of molecular dynamics simulation. We consider sessile drops on a lyophilic (attractive for the monomers) region of circular shape with radius R(D) while the remaining part of the substrate is lyophobic. The variation of the droplet shape, including its contact angle, with R(D) is studied, and the density profiles across these droplets also are obtained. In addition, the interaction of droplets adsorbed on two walls forming a slit pore with two lyophilic circular regions just opposite of one another is investigated, paying attention to the formation of a liquid bridge between both walls. A central result of our study is the measurement of the force between the two substrate walls at varying wall separation as well as the kinetics of droplet merging. Our results are compared to various phenomenological theories developed for liquid droplets of mesoscopic rather than nanoscopic size.     

Modeling of adsorption on nongraphitized carbon surface: GCMC simulation studies and comparison with experimental data

We model nongraphitized carbon black surfaces and investigate adsorption of argon on these surfaces by using the grand canonical Monte Carlo simulation. In this model, the nongraphitized surface is modeled as a stack of graphene layers with some carbon atoms of the top graphene layer being randomly removed. The percentage of the surface carbon atoms being removed and the effective size of the defect (created by the removal) are the key parameters to characterize the nongraphitized surface. The patterns of adsorption isotherm and isosteric heat are particularly studied, as a function of these surface parameters as well as pressure and temperature. It is shown that the adsorption isotherm shows a steplike behavior on a perfect graphite surface and becomes smoother on nongraphitized surfaces. Regarding the isosteric heat versus loading, we observe for the case of graphitized thermal carbon black the increase of heat in the submonolayer coverage and then a sharp decline in the heat when the second layer is starting to form, beyond which it increases slightly. On the other hand, the isosteric heat versus loading for a highly nongraphitized surface shows a general decline with respect to loading, which is due to the energetic heterogeneity of the surface. It is only when the fluid-fluid interaction is greater than the surface energetic factor that we see a minimum-maximum in the isosteric heat versus loading. These simulation results of isosteric heat agree well with the experimental results of graphitization of Spheron 6 (Polley, M. H.; Schaeffer, W. D.; Smith, W. R. J. Phys. Chem. 1953, 57, 469; Beebe, R. A.; Young, D. M. J. Phys. Chem. 1954, 58, 93). Adsorption isotherms and isosteric heat in pores whose walls have defects are also studied from the simulation, and the pattern of isotherm and isosteric heat could be used to identify the fingerprint of the surface.     

Low temperature adsorption of CO2 in carbonaceous wedge pores: a Monte Carlo simulation study

Adsorption - A systematic study of carbon dioxide in wedge pores under subcritical conditions were conducted with Grand Canonical Monte Carlo simulation. The effects of various factors:...     

The effects of energy sites on adsorption of Lennard-Jones fluids and phase transition in carbon slit pore of finite length a computer simulation study

A Monte Carlo simulation method is used to study the effects of adsorption strength and topology of sites on adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length. Argon is used as a model adsorbate, while the adsorbent is modeled as a finite carbon slit pore whose two walls composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. Impurities having well depth of interaction greater than that of carbon atom are assumed to be grafted onto the surface. Different topologies of the impurities; corner, centre, shell and random topologies are studied. Adsorption isotherms of argon at 87.3 K are obtained for pore having widths of 1, 1.5 and 3 nm using a Grand Canonical Monte Carlo simulation (GCMC). These results are compared with isotherms obtained for infinite pores. It is shown that the surface heterogeneity affects significantly the overall adsorption isotherm, particularly the phase transition. Basically it shifts the onset of adsorption to lower pressure and the adsorption isotherms for these four impurity models are generally greater than that for finite pore. The positions of impurities on solid surface also affect the shape of the adsorption isotherm and the phase transition. We have found that the impurities allocated at the centre of pore walls provide the greatest isotherm at low pressures. However when the pressure increases the impurities allocated along the edges of the graphene layers show the most significant effect on the adsorption isotherm. We have investigated the effect of surface heterogeneity on adsorption hysteresis loops of three models of impurity topology, it shows that the adsorption branches of these isotherms are different, while the desorption branches are quite close to each other. This suggests that the desorption branch is either the thermodynamic equilibrium branch or closer to it than the adsorption branch.     

Dynamics of benzene guest inside a self-assembled cylindrical capsule: a combined solid-state 2H NMR and molecular dynamics simulation study

The reorientational dynamics of benzene-d(6) molecules hosted into the cavity of a cavitand-based, self-assembled capsule was investigated by Molecular Dynamics (MD) simulations and temperature-dependent solid-state (2)H NMR spectroscopy. MD simulations were preliminarily performed to assess the motional models of the guest molecules inside the capsules. An in-plane fast reorientation of the benzene guest around the C(6) symmetry axis (B1 motion), characterized by correlation times of the order of picoseconds, was predicted with an activation barrier ( approximately 8 kJ/mol) very similar to that found for neat benzene in the liquid state. An out-of-plane reorientation corresponding to a nutation of the C(6) symmetry axis in a cone angle of 39 degrees (B2 motion, 373 K) with an activation barrier ( approximately 39 kJ/mol) definitely larger than that of liquid benzene was also anticipated. In the temperature range 293-373 K correlation times of the order of a nanosecond have been calculated and a transition from fast to slow regime in the (2)H NMR scale has been predicted between 293 and 173 K. (2)H NMR spectroscopic analysis, carried out in the temperature range 173-373 K on the solid capsules containing the perdeuterated guest (two benzene molecules/capsule), confirmed the occurrence of the B1 and B2 motions found in slow exchange in the (2)H NMR time scale. Line shape simulation of the (2)H NMR spectral lines permitted defining a cone angle value of 39 degrees at 373 K and 35 degrees at 173 K for the nutation axis. The T(1) values measured for the (2)H nuclei of the encapsulated aromatic guest gave correlation times and energetic barrier for the in-plane motion B1 in fine agreement with theoretical calculation. The experimental correlation time for B2 as well as the corresponding energetic barrier are in the same range found for B1. A molecular mechanism for the encapsulated guest accounting for the B1 and B2 motions was also provided.     

Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: for very long pores, the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place, fluctuations occur between vapor- and liquidlike states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality), the system at phase coexistence is in an axially inhomogeneous multidomain state, where long cylindrical liquid- and vaporlike domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures, the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in the axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multidomain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical temperature should increase logarithmically with the length of the pore.     

Packing density and structure effects on energy-transfer dynamics in argon collisions with organic monolayers

A combined experimental and molecular-dynamics simulation study has been used to investigate energy-transfer dynamics of argon atoms when they collide with n-alkanethiols adsorbed to gold and silver substrates. These surfaces provide the opportunity to explore how surface structure and packing density of alkane chains affect energy transfer in gas-surface collisions while maintaining the chemical nature of the surface. The chains pack standing up with 12 degrees and 30 degrees tilt angles relative to the surface normal and number densities of 18.9 and 21.5 A(2)molecule on the silver and gold substrates, respectively. For 7-kJmol argon scattering, the two surfaces behave equivalently, fully thermalizing all impinging argon atoms. In contrast, these self-assembled monolayers (SAMs) are not equally efficient at absorbing the excess translational energy from high-energy, 35 and 80 kJmol, argon collisions. When high-energy argon atoms are scattered from a SAM on silver, the fraction of atoms that reach thermal equilibrium with the surface and the average energy transferred to the surface are lower than for analogous SAMs on gold. In the case of argon atoms with 80 kJmol of translational energy scattering from long-chain SAMs, 60% and 45% of the atoms detected have reached thermal equilibrium with the monolayers on gold and silver surfaces, respectively. The differences in the scattering characteristics are attributed to excitation efficiencies of different types of surface modes. The high packing density of alkyl chains on silver restricts certain low-energy degrees of freedom from absorbing energy as efficiently as the lower-density monolayers. In addition, molecular-dynamics simulations reveal that the extent to which argon penetrates into the monolayer is related to packing density. For argon atoms with 80-kJmol incident energy, we find 16% and 7% of the atoms penetrate below the terminal methyl groups of C(10) SAMs on gold and silver, respectively.     

Model liquid crystalline dimers in a slit pore: Monte Carlo simulation study

ABSTRACT

In this work, we present results from (isobaric–isothermal) Monte Carlo Simulation studies of liquid crystalline dimer systems confined in a slit pore. Liquid crystalline dimer systems of various spacer numbers have been considered. Surface-induced conformational and alignment properties of these systems at different pressures under homeotropic anchoring condition have been investigated. We have used easily manageable coarse grained force fields to model both monomer–monomer and monomer–substrate interaction potentials. According to the simulated result, the anchoring of dimers to the surface and orientation of mesogenic units with respect to the surface normal seem to depend on the spacer number for messogen attractive confinement. Dimers with lower spacer number are able be adsorbed to the surface and most of their mesogens are oriented along the surface normal even at lower pressure. Those with larger spacer number are distributed throughout the volume at lower pressure. In the case of mesogen repulsive confinement, most of the dimers are adsorbed to the surface and most mesogens are randomly oriented at low pressure. As the pressure gets higher, the adsorption and orientability increase depending on the type of confinement and spacer number. As a result, clear submolecular partitioning and smectic A like structure have been identified.     

Sr掺杂对CaO(100)表面吸附甲醇影响的分子模拟

借助分子模拟手段,研究了锶掺杂对氧化钙表面甲醇吸附行为的影响。构建了甲醇在CaO(100)和CaO(100)-Sr表面吸附的模型,计算了甲醇在氧化钙表面的吸附能和解离活化能,分析了甲醇在氧化钙表面成键的态密度以及锶掺杂前后甲醇在氧化钙表面电荷布局和差分电荷密度,评估了锶掺杂量对氧化钙表面甲醇吸附性能的影响。结果表明,锶掺杂能够显著强化氧化钙对甲醇的吸附性能,降低甲醇的解离活化能,且吸附性能随锶掺杂量的增加而增强;甲醇在氧化钙表面吸附时活化,锶掺杂后活化程度增加。     

Simulation study of two-dimensional phase transitions of argon on graphite surface and in slit micropores

Molecular simulation has been increasingly used in the analysis and modeling of gas adsorption on open surfaces and in porous materials because greater insight could be gained from such a study. In case of homogeneous surfaces or pore walls the adsorption behavior is often complicated by the order–disorder transition. It is shown in our previous publications (Ustinov and Do, Langmuir 28:9543–9553, 2012a; Ustinov and Do, Adsorption 19:291–304, 2013) that once an ordered molecular layer has been formed on the surface, the lattice constant depends on the simulation box size, which requires adjusting the box dimensions parallel to the surface for each value of loading. It was shown that this can be accomplished with the Gibbs–Duhem equation, which results in decreasing lattice constant with an increase of the amount adsorbed. The same feature is expected to be valid for gas adsorption in narrow pores, but this has not been analyzed in the literature. This study aims at an extension of our approach to adsorption in slit graphitic pores using kinetic Monte Carlo method (Ustinov and Do, J Colloid Interface Sci 366:216–223, 2012b). The emphasis rests on the thermodynamic analysis of the two-dimensional (2D) ordering transition and state of the ordered phase; if the ordered phase exists in narrow slit pores, simulation with constant volume box always leads to erroneous results, for example, seemingly incompressible adsorbed phase. We proposed a new approach that allows for modeling thermodynamically consistent adsorption isotherms, which can be used as a basis for further refinement of the pore size distribution analysis of nanoporous materials.     

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.     

Molecular simulation study on adsorption of methanol/water mixtures in mesoporous silicas modified pore surface silylation

Two types of molecular simulation techniques have been utilized to investigate adsorption of methanol/water mixtures in a mesoporous silica with a hydrophobic pore surface: the NVT-ensemble Molecular Dynamics method with the melt-quench algorithm for modeling a fully-silylated mesoporous silica and the μVT-ensemble Orientaional-Biased Monte Carlo method for calculating adsorption isotherms. Adsorption isotherms of methanol and water at 333 K are calculated for an equi-relative-pressure mixture (each component has the same relative pressure which is defined as the ratio of the partial pressure to the saturation pressure of the pure gas) together with pure gases. In the case of the pure gas, water hardly adsorb even at elevated pressures, while the adsorption isotherm for methanol shows the condensable adsorption. On the other hand, in the case of the mixture, water molecules are substantially adsorbed along with methanol molecules, showing an isotherm representing the condensation mechanism. In addition, it is found that the separation factor of methanol to water is the highest in the case of monolayer adsorption from a liquid mixture.