Surface tensions of linear and branched alkanes from Monte Carlo simulations using the anisotropic united atom model

The anisotropic united atoms (AUA4) model has been used for linear and branched alkanes to predict the surface tension as a function of temperature by Monte Carlo simulations. Simulations are carried out for n-alkanes ( n-C5, n-C6, n-C7, and n-C10) and for two branched C7 isomers (2,3-dimethylpentane and 2,4-dimethylpentane). Different operational expressions of the surface tension using both the thermodynamic and the mechanical definitions have been applied. The simulated surface tensions with the AUA4 model are found to be consistent within both definitions and in good agreement with experiments.     

Monte Carlo simulation of polymer adsorption

We developed and employed the incremental gauge cell method to calculate the chemical potential (and thus free energies) of long, flexible homopolymer chains of Lennard-Jones beads with harmonic bonds. The free energy of these chains was calculated with respect to three external conditions: in the zero-density bulk limit, confined in a spherical pore with hard walls, and confined in a spherical pore with attractive pores, the latter case being an analog of adsorption. Using the incremental gauge cell method, we calculated the incremental chemical potential of free polymer chains before and after the globual-random coil transitions. We also found that chains confined in attractive pores exhibit behaviors typical of low temperature physisorption isotherms, such as layering followed by capillary condensation.     

A Monte Carlo simulation study of branched polymers

Monte Carlo simulations are presented for the static properties of highly branched polymer molecules. The molecules consist of a semiflexible backbone of hard-sphere monomers with semiflexible side chains, also composed of hard-sphere monomers, attached to either every backbone bead or every other backbone bead. The conformational properties and structure factor of this model are investigated as a function of the stiffness of the backbone and side chains. The average conformations of the side chains are similar to self-avoiding random walks. The simulations show that there is a stiffening of the backbone as degree of crowding is increased, for example, if the branch spacing is decreased or side chain length is increased. The persistence length of the backbone is relatively insensitive to the stiffness of the side chains over the range investigated. The simulations reproduce most of the qualitative features of the structure factor observed in experiment, although the magnitude of the stiffening of the backbone is smaller than in experiment.     

Unusual chain length dependent adsorption of linear and branched alkanes on UiO-66

In this work, the zero coverage adsorption properties of C₅–C₁₀ n- and iso-alkanes on the UiO-66, UiO-66-Me and UiO-66-NO₂ metal–organic frameworks are studied by gas phase pulse chromatography. Analysis of enthalpy values, entropy values, Gibbs free energies and Henry constants reveals unusual chain length dependent adsorption behaviour of linear and branched alkanes, caused by the complex structure of the zirconium metal–organic framework UiO-66. The UiO-66 structure consists of a small, tetrahedral and large, octahedral cage. It is shown that at specific carbon chain lengths (e.g. C₆–C₇ for n-alkanes), distinctive jumps in adsorption enthalpy, entropy values and Henry constants occur. This chain length dependent effect is even more pronounced for 2- and 3-methyl alkanes and double branched alkanes. This distinctive shift in adsorption behaviour occurs at a molecular size that coincides with the cavity dimensions of the smallest, tetrahedral cage. The resulting selective adsorption arises from confinement effects and is function of both the molecular shape and size.     

Monte Carlo simulation of size exclusion chromatography for randomly branched and crosslinked polymers

The elution curves of size exclusion chromatography for nonlinear polymers formed through random branching and crosslinking of long polymer chains were simulated with a Monte Carlo method. We considered two types of measured molecular weight distributions (MWDs): (1) the MWD calibrated relative to standard linear polymers and (2) the MWD obtained with a light scattering (LS) photometer in which the weight‐average molecular weight of polymers within the elution volume is determined directly. The calibrated MWDs clearly underestimate the molecular weights for both randomly branched and crosslinked polymers, and this technique can be used to assess the degree of deviation from the true MWD. When the primary chains conform to the most probable distribution, the calibrated MWD can be estimated reasonably well with the Zimm–Stockmayer equation for the g factor with the help of the relationship between the average number of branch points per molecule and the degree of polymerization. However, the LS method gives good estimates of the true MWD for both randomly branched and crosslinked polymers, although the agreement is better for the branched ones. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 2009–2018, 2000     

Storage and separation of CO2 and CH4 in silicalite, C168 schwarzite, and IRMOF-1: a comparative study from Monte Carlo simulation

Storage of pure CO2 and CH4 and separation of their binary mixture in three different classes of nanostructured adsorbents--silicalite, C168 schwarzite, and IRMOF-1--have been compared at room temperature using atomistic simulation. CH4 is represented as a spherical Lennard-Jones molecule, and CO2 is represented as a rigid linear molecule with a quadrupole moment. For pure component adsorption, CO2 is preferentially adsorbed than CH4 in all the three adsorbents over the pressure range under this study, except in C168 schwarzite at high pressures. The simulated adsorption isotherms and isosteric heats match closely with available experimental data. A dual-site Langmuir-Freundlich equation is used to fit the isotherms satisfactorily. Compared to silicalite and C168 schwarzite, the gravimetric adsorption capacity of pure CH4 and CO2 separately in IRMOF-1 is substantially larger. This implies that IRMOF-1 might be a potential storage medium for CH4 and CO2. For adsorption from an equimolar binary mixture, CO2 is preferentially adsorbed in all three adsorbents. Predictions of mixture adsorption with the ideal-adsorbed solution theory on the basis of only pure component adsorption agree well with simulation results. Though IRMOF-1 has a significantly higher adsorption capacity than silicalite and C168 schwarzite, the adsorption selectivity of CO2 over CH4 is found to be similar in all three adsorbents.     

Capillary phase transitions of linear and branched alkanes in carbon nanotubes from molecular simulation

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.     

Monte Carlo simulation studies of the size and shape of linear and star‐branched polymers embedded in the tetrahedral lattice

In recent years the author investigated the size and shape of linear and star‐branched model chains as a function of functionality (number of arms) F and chain length (total number of beads) N for several lattice and off‐lattice random walk models as well as for self‐avoiding random walks (embedded in the tetrahedral lattice) subjected to various thermodynamic conditions. Not only mean values – characteristic quantities averaged over a large number of configurations of given length, functionality and thermodynamic conditions – have been computed, but also distributions of shape parameters, and the correlation operative between size and shape has been explored. The present feature article in principle summarizes these results, however, the data given in part are recomputed for still larger sample sizes and chain‐lengths as in the underlying papers. In addition to stars with F = 4, 8 and 12 arms, so far unpublished investigations on stars with F = 3, F = 6 and F = 10 arms are presented and discussed.     

Molecular simulation of adsorption and separation of mixtures of short linear alkanes in pillared layered materials at ambient temperature

Grand canonical Monte Carlo and configurational-bias Monte Carlo techniques are carried out to simulate the adsorption of ternary and quaternary mixtures of short linear alkanes, involving methane, ethane, propane, and n-butane, in pillared layered materials at ambient temperature, T=300 K. In the simulation, a pillared layered pore is modeled by a uniform distribution of pillars between two layered walls built by making two separate talc lamellas parallel each other with a given size of interlayer distance. The interaction between fluid molecules and two layered walls is measured by storing potentials calculated in advance at a series of grid points. The interaction between fluid molecules and pillars is also calculated by a site-to-site method. The potential model proposed in this work is proved to be effective because of the simulation result being good agreement with the experimental data for the adsorption of nitrogen at 77 K. Then, the adsorption isotherms of mixtures of short linear alkanes in pillared layered pores with three different porosities psi=0.98, 0.93 and 0.85, and three pore widths H=1.02, 1.70 and 2.38 nm at 300 K are obtained by taking advantage of the model. The simulation results tell us that the longer chain component is preferentially adsorbed at low pressures, and its adsorption increases and then decreases as the pressure increases while the shorter chain component is still adsorbed at high pressures. Moreover, the sorption selectivity of pillared layered materials for the longest chain component in alkane mixtures increases as the mole fraction of methane in the gas phase increases. The selectivity of pillared layered materials for the longest chain component in alkane mixtures also increases as the pore width decreases and the porosity increases.     

Monte Carlo simulation of n-alkane adsorption isotherms in carbon slit pores

The single component adsorption of alkanes in carbon slit pores was studied using configurational-biased grand canonical Monte Carlo simulations. Wide ranges of temperature, pressure, alkane chain length, and slit height were studied to evaluate their effects on adsorption. Adsorption isotherms and density and orientation profiles were calculated. The behavior of long alkanes at high temperatures was found to be similar to short alkanes at lower temperatures. This suggests that the isotherms may be related through the Polanyi potential theory.     

Density functional steric analysis of linear and branched alkanes

Branched alkane hydrocarbons are thermodynamically more stable than straight-chain linear alkanes. This thermodynamic stability is also manifest in alkane bond separation energies. To understand the physical differences between branched and linear alkanes, we have utilized a novel density functional theory (DFT) definition of steric energy based on the Weiz?cker kinetic energy. Using the M06-2X functional, the total DFT energy was partitioned into a steric energy term (E(s)[ρ]), an electrostatic energy term (E(e)[ρ]), and a fermionic quantum energy term (E(q)[ρ]). This analysis revealed that branched alkanes have less (destabilizing) DFT steric energy than linear alkanes. The lower steric energy of branched alkanes is mitigated by an equal and opposite quantum energy term that contains the Pauli component of the kinetic energy and exchange-correlation energy. Because the steric and quantum energy terms cancel, this leaves the electrostatic energy term that favors alkane branching. Electrostatic effects, combined with correlation energy, explains why branched alkanes are more stable than linear alkanes.     

Modified Monte Carlo simulation of adsorption in a spatially inhomogeneous porous system

A modified Monte Carlo method in conjunction with the canonical and grand canonical ensembles is proposed for simulating adsorption in spatially inhomogeneous porous systems. Unlike the traditional Monte Carlo simulation in terms of the grand canonical ensemble, the simulation for the regions of pore space having no direct communication with the bulk phase is performed in local conditions of the canonical ensemble.     

Efficient simulation of binary adsorption isotherms using transition matrix Monte Carlo

Molecular simulations of binary adsorption in porous materials are a useful complement to experimental studies of mixture adsorption. Most molecular simulations of binary adsorption are performed using grand canonical Monte Carlo (GCMC) to independently examine a range of state points of interest. A disadvantage of this approach is that it only yields information at a discrete set of state points; therefore, if a complete isotherm is required for arbitrary conditions, some type of data fitting or interpolation must be used in combination with the GCMC data. We show that the transition matrix Monte Carlo (TMMC) method of Shen and Errington (Shen, V. K.; Errington, J. R. J. Chem.Phys. 2005, 122, 064508) is well-suited to simulation of binary adsorption in porous materials. At the completion of a TMMC simulation, the adsorption isotherm for all possible bulk phase compositions and pressures is available without data fitting or interpolation. It is also straightforward to use results from TMMC to compute derivatives of the isotherm such as the mixture thermodynamic correction factors, partial differential ln f(i)/partial differential ln c(j), again without data fitting or interpolation. This approach should be useful in contexts where information on the full adsorption isotherm is needed, such as the design of adsorption- or membrane-based separations.     

Monte Carlo simulation of phase separation kinetics in a concentrated polymer solution

Ｔｈｅｐｈａｓｅｂｅｈａｖｉｏｒｉｎｍｕｌｔｉｐｌｅｃｏｍｐｏｎｅｎｔｐｏｌｙｍｅｒｓｃｏｎｓｔｉｔｕｔｅｓａｌｏｎｇｓｔａｎｄｉｎｇａｃｔｉｖｅａｃａｄｅｍｉｃｓｕｂｊｅｃｔｂｏｔｈｉｎｐｏｌｙｍｅｒｓｃｉｅｎｃｅａｎｄｃｏｎｄｅｎｓｅｄｓｔａｔｅｐｈｙｓｉｃｓ.Ｉｔｉｓｅｓｐｅｃｉａｌｌｙｓｉｇｎｉｆｉｃａｎｔｉｎｇｕｉｄｉｎｇｔｈｅｆａｂｒｉｃａｔｉｏｎｏｆｐｏｌｙｍｅｒａｌｌｏｙｓ[1].Ｄｕｒｉｎｇｔｈｅｌａｓｔｄｅｃａｄｅｓｍｕｃｈａｔｔｅｎｔｉｏｎｈａｓｂｅｅｎｐａｉｄｔｏｔｈｅｃｏｍ…     

Lattice Monte Carlo simulation of polymer adsorption at an interface, 1. Monodisperse polymer

Adsorption of a monodisperse polymer at a solid-liquid interface is comprehensively studied by Monte Carlo simulation. The distributions of total segment density and different adsorption configurations including trains, loops and tails are obtained. Effects of reduced exchange interaction energies $ \tilde \varepsilon $, bulk concentrations ϕ^*, reduced adsorption energies $ \tilde \varepsilon_a $ and chain lengths r on those distributions are studied. Comparisons with predictions of the Scheutjens-Fleer (SF) theory are also provided. Generally, the chain molecules are more easily adsorbed at an interface in non-solvents than in good solvents. Longer chains are more likely to be adsorbed than shorter ones. The reduced adsorption energy and the bulk concentration have shown strong effects on the segment-density distributions. In addition, the thickness of the adsorption layer is mainly determined by the extension of tails into the bulk solution, which are in turn determined by the chain length. The trains, loops and tails are overwhelmingly short. On the other hand, the amounts of trains and loops are usually much greater than that of tails. Though not perfect, satisfactory agreement is found in comparison with the theoretical predictions of the SF theory.     

Monte Carlo simulation for the adsorption of symmetric triblock copolymers II. Adsorption layer information

By using Monte Carlo simulation, adsorption of both end-adsorbed and middle-adsorbed symmetric triblock copolymers from a non-selective solvent on an impenetrable surface has been studied. Influences of the adsorption energy, the bulk concentration, the chain composition and the chain length on the adsorption behavior including the surface coverage, the adsorption amount and the layer thickness are presented. It is shown that the total surface coverage for both end-adsorbed and middle-adsorbed copolymers increases monotonically as the bulk concentration increases. The higher the adsorption energy and the more the attractive segments, the higher the total surface coverage is exhibited. Surface coverage θ decreases with increasing the length of the non-attractive segments, but the product of θ and the proportion of the non-attractive segments in a triblock copolymer chain is nearly independent of the chain length. The adsorption amount increases almost monotonically with the bulk concentration. The logarithm of the adsorption amount is a linear function of the reciprocal of the reduced temperature. When the adsorption energy is large, the adsorption amount exhibits a maximum as the composition of the attractive segment increases. The adsorption isotherms of copolymers with different length of the non-attractive segments can be mapped onto a single curve under certain energy indicating that copolymers with different chain length have the same adsorption amount. The adsorption layer thickness for the end-adsorbed copolymers decreases as the energy and the number of adsorbing segments increases. The longer non-attractive segments, the larger adsorbed layer thickness is found. The tails mainly governs the adsorption layer thickness.     

Gibbs ensemble Monte Carlo simulation of adsorption for model surfactant solution in confined slit pores

An isobaric-isothermal Gibbs ensemble Monte Carlo simulation has been carried out to study the adsorption of a model surfactant/solvent mixture in slit nanopores. The adsorption isotherms, the density distributions, and the configuration snapshots were simulated to illustrate the adsorption and self-assembly behaviors of the surfactant in the confined pores. The adsorption isotherms are stepwise: a two-step curve for the smaller (30 A) pore and a three-step one for the larger (50 A) pore. The adsorption isotherms and the interfacial aggregate structure of the surfactants in the pores with various sizes show a qualitatively consistent performance with the previous experimental observation. The micelle size distributions of the adsorbed surfactant aggregates have been analyzed in order to understand the adsorption mechanism, which suggests that the step rise in the surfactant adsorption is associated with the considerable formation of the micelle aggregates in the confined pores. The effect of the interaction between the pore surface and the surfactant on the adsorption behavior has also been investigated. The simulation results indicate that a change in the interaction can modify the shape of adsorption isotherms. A nonlinear mathematical model was used to represent the multistep adsorption isotherms. A good agreement between the model fitting and the simulation data was obtained for both the amount of adsorption and the jump point concentration.     

Monte Carlo simulation of water adsorption in hydrophobic MFI zeolites with hydrophilic sites

The effect of strong and weak hydrophilic sites, Al atoms with associated extraframework Na cations and silanol nests, respectively, in high-silica MFI zeolites on water adsorption was investigated using Monte Carlo simulations. For this purpose, a new empirical model to represent potential energy interactions between water molecules and the MFI framework was developed, which reproduced the hydrophobic characteristics of a siliceous MFI-type zeolite, silicalite-1, with both the vapor-phase adsorption isotherm and heats of adsorption at 298 K being in good agreement with experimental data. The proposed model is also compatible with previous hydrocarbon potential models and can be used in the adsorption simulations of VOC-water mixtures. Adsorption simulations revealed that strongly hydrophilic Al sites in Na-ZSM-5 zeolites coordinate two water molecules per site at low coverage, which promotes water clustering in the vicinity of these sites. However, weakly hydrophilic silanol nests in silicalite-1 are in coordination with a single water molecule per site, which does not affect the adsorption capacity significantly as expected. However, even in the presence of 0.125 silanol nest per unit cell, the increase in the heat of adsorption at low coverage is drastic.     

On the efficiency of azimuthal and rotational splitting for Monte Carlo simulation of clinical linear accelerators

This article presents rotational splitting as a modification to the sampling process of the azimuthal angle used in the variance-reduction technique of azimuthal particle redistribution, with the goal to improve the efficiency of this variance-reduction technique for the Monte Carlo simulation of radiation transport in clinical linear accelerators. Using a constant azimuthal angle, instead of a random one, in the azimuthal particle redistribution technique, increases the efficiency of the simulation of a clinical linear accelerator by about 30% and reduces the latent variance of a 10×10 cm² phase space by about 40%.