A comparative study of the adsorption of water and methanol in zeolite BEA: a molecular simulation study

Grand canonical Monte Carlo simulations were carried out to study the equilibrium adsorption concentration of methanol and water in all-silica BEA zeolite and HBEA zeolites with different Si/Al ratios over a wide range of temperatures and loadings. These zeolites have oval-shaped channels with one side longer than the other. Water sorption into the hydrophobic BEA zeolite had a sharp transition with its sorption going from zero to near full capacity over a very small pressure range. Methanol sorption was much more gradual with respect to pressure. With the addition of hydrophilic sites for the HBEA zeolites by decreasing the Si/Al ratio, adsorption at lower pressures increased significantly for water and methanol. At higher loadings, water and methanol adsorption were found to behave in fundamentally different ways. Water structures in the zeolite channels formed hydrogen-bonded chains while maximising contact with the surfaces on the longer edges of the zeolite channels. Methanol molecules, in contrast, formed very few hydrogen bonds between themselves, with their hydroxyl groups primarily binding with surface of the shorter edge of the zeolite channels and their methyl groups located near the middle of the zeolite channels. The addition of hydrophilic groups in the HBEA zeolites strongly influenced positions of the methanol hydroxyl groups at high loadings, but did not have a significant effect on water structure.     

Methane adsorption in several series of newly synthesised metal-organic frameworks: a molecular simulation study

In this work, the adsorption of methane in several newly synthesised metal-organic frameworks (MOFs) is studied using grand canonical Monte Carlo simulations. The factors influencing the adsorption behaviour of methane in different series of materials were investigated at low, moderate and high pressures, respectively. The simulation results show that the isosteric heat of adsorption at infinite dilute, the accessible surface area and the free volume of materials play important roles in methane adsorption, while the main factors are different in different ranges of pressure. These are consistent with the rules obtained from well-studied MOFs by the Snurr group and by Wang, proving that they are applicable to the newly synthesised MOFs. Based on these, the influencing factors were further discussed by classifying the MOFs through the pore topology, and it was found that the pore topology should not be ignored when performing efficient modifications of MOFs at moderate and high pressures. This work may provide useful information to design novel MOFs with high methane uptake in the future.     

RASPA: molecular simulation software for adsorption and diffusion in flexible nanoporous materials

A new software package, RASPA, for simulating adsorption and diffusion of molecules in flexible nanoporous materials is presented. The code implements the latest state-of-the-art algorithms for molecular dynamics and Monte Carlo (MC) in various ensembles including symplectic/measure-preserving integrators, Ewald summation, configurational-bias MC, continuous fractional component MC, reactive MC and Baker's minimisation. We show example applications of RASPA in computing coexistence properties, adsorption isotherms for single and multiple components, self- and collective diffusivities, reaction systems and visualisation. The software is released under the GNU General Public License.     

Molecular simulation study of the quantum effects of hydrogen adsorption in metal-organic frameworks: influences of pore size and temperature

In this work a systematic molecular simulation study was performed to investigate the influence of pore size and temperature on the quantum effects of hydrogen adsorption in metal-organic frameworks (MOFs) with temperature varied from 40 to 120 K. To do this, three isoreticular MOFs (IRMOFs) with different pore sizes were adopted, and quadratic Feynman–Hibbs (FH) effective potential was introduced to consider the quantum effect. The results show that quantum effects diminish with increasing pore size of IRMOFs at lower pressure (loading), while the opposite trend appears at higher pressure (loading). Through the simulations it is also found that the quantum effects may be dominantly determined by the adsorbate–adsorbate or adsorbate–MOFs interactions with the varying of pressure (loading). In addition, the results also indicate how the temperature influences the quantum effects of H₂ adsorption in MOFs within the pressure range considered.     

Enantioseparation of phenylglycinol in chiral-modified zeolite HY: a molecular simulation study

A mechanism has been proposed for the separation of valinol enantiomers using a chiral-modified zeolite HY (i.e., zeolite HY containing (+)-(1R;2R)-hydrobenzoin) Molecular modeling of chiral-modified zeolite HY employed in enantioselective separation. Jirapongphan SS, Warzywoda J, Budil DE, Sacco A Jr. Chirality 2007; in press, which accurately predicted the experimentally measured enantioseparation. This methodology has been applied to predict the separation of an enantiomeric pair of phenylglycinol molecules (a precursor in the synthesis of HIV-1 protease inhibitors) using the modified zeolite HY as a CSP. Phenylglycinol and valinol molecules are similar in terms of the presence of polar (i.e., amine and hydroxyl) groups. These functional groups are important in the proposed chiral discrimination. Supercage-based docking simulations yielded an enantioselectivity of 1.3 with (+)-(S)-phenylglycinol molecule better retained in the zeolite. Also, the simulations predicted two binding modes that were the same as those in the valinol system. This suggests that specific structural features (i.e., number and type of polar groups), which generate the hypothesized binding modes, are required in an enantioseparation utilizing the chiral-modified zeolite HY.     

Cation redistribution upon dehydration of Na58Y faujasite zeolite: a joint neutron diffraction and molecular simulation study

Using neutron scattering and Monte Carlo simulation, we investigate the distribution of cations in Na₅₈Y faujasite upon (de)hydration. We introduce a new method for the assignment of cations to specific sites in molecular simulations from their local environment. This allows us to bypass the need of the coordinates of crystallographic sites, which vary as water adsorption induces changes in the zeolite framework structure. Although the agreement between experiments and simulation is excellent at high temperature, some differences are observed below 150°C. We show that these differences are due to the presence of water and that temperature itself as well as adsorption-induced deformation of the framework play a less important role. We demonstrate the migration of sodium to sites III upon water adsorption, not observed for other Si:Al ratios.     

On the existence of a hysteresis loop in open and closed end pores

We have studied the adsorption of argon at 87 K in slit pores of finite length with a smooth graphitic potential, open at both ends or closed at one end. Simulations were carried out using conventional GCMC (grand canonical Monte Carlo) or kMC (kinetic Monte Carlo) in the canonical ensemble with extremely long Markov chain, of at least 2 × 10⁸ configurations; selected simulations with much longer Markov chains do not show any change in the results. When the pore width is in the micropore range (0.65 nm), type I isotherms are obtained for both pore models and for both simulation methods. However, wider pores (1, 2 and 3 nm in width) all exhibit hysteresis loops in the GCMC simulations, while in the canonical ensemble simulations, the isotherms pass through a sigmoid van der Waals type loop in the transition region. This loop locates the true equilibrium transition. For the pores with one closed end, this transition is close to, or coincides with, the adsorption branch of the GCMC hysteresis loop, but for the open-ended pores, it is more closely associated with the desorption branch. In a separate study of adsorption hysteresis in an infinitely long slit pore, using both simulation techniques, the van der Waals loop follows the adsorption branch of the GCMC isotherm to the transition, then reverts to a long vertical section that falls midway between the two hysteresis branches and finally moves to the desorption transition close to the evaporation pressure. An examination of molecular distributions inside the pores reveals two coexisting phases in the canonical simulations, whereas in the grand canonical simulations, the molecules are uniformly distributed along the length of the pores.     

Effect of Si/Al ratio on tetralin adsorption on Y zeolite: a DFT study

Probing the adsorption of tetralin on zeolite is of prime scientific and industrial importance with the aim to upgrade the industrial process of tetralin cracking. In this work, the effect of Si/Al ratio ranging from 12 to 39 on tetralin adsorption property on Y zeolite is studied by DFT calculations. Tetralin adsorption on Y zeolite corresponds to a π-stacking adsorption mechanism between double bonds of aromatic ring and Brønsted acid sites. Therefore, the number of Brønsted acid sites influences the adsorption properties. Lower Si/Al ratio with more Brønsted acid sites interacting with the aromatic ring of tetralin leads to a higher adsorption energy. Furthermore, the charge and frontier molecular orbital analysis are also performed to understand the influence of Si/Al ratio on adsorption performance. Y zeolite with lower Si/Al ratio shows larger charge difference values and lower HOMO–LUMO gap, which directly manifests the stronger adsorption ability of tetralin and indicates bigger possibility of reacting.     

Monte Carlo study of sI hydrogen hydrates

In this study, we perform grand canonical Monte Carlo simulations to evaluate the hydrogen storage capacity of structure I (sI) hydrogen hydrates at pressures up to 500 MPa. Initially, we calculate the upper limit of H₂ content of sI hydrates by studying the hypothetical sI hydrate, where H₂ is the single guest component. It is found that the storage capacity of the hypothetical pure H₂ sI hydrate could reach 3.5 wt% at 500 MPa and 274 K. Depending on pressure, the large cavities of the pure H₂ hydrate can accommodate up to three H₂ molecules while the small ones are singly occupied at most, even at pressures as high as 500 MPa, without any double occupancy being observed. Subsequently, the binary H₂–ethylene oxide (EO) hydrate is examined. In this case, the large cavities are occupied by a single EO molecule while the small cavities can accommodate at most a single H₂ molecule. Such configuration results in a maximum H₂ content of only 0.37 wt%. The hydrogen storage capacity does not improve significantly even in case when EO is replaced by a component with smaller molecular weight.     

Influence of substrate pattern on the adsorption of HP lattice proteins

With the highly simplified hydrophobic-polar model representation of a protein, we can study essential qualitative physics without an unnecessarily large computational overhead. Using Wang-Landau sampling in conjunction with a set of efficient Monte Carlo trial moves, we studied the adsorption of short HP lattice proteins on various simple patterned substrates and in particular for checkered patterned surfaces. A set of single-site mutated HP proteins is used to investigate the role of hydrophobicity of a protein chain and surface pattern for substrates with various pattern cell sizes relative to the protein’s native configuration. For most cases, we found that the adsorption transition occurs at a lower temperature, while the hydrophobic core formation is less affected. The flattening procedure after the HP protein is adsorbed is more sensitive to the change in surface patterns and single-site mutations. These observations stay valid for both strongly and weakly attractive surfaces.     

The self-diffusion and hydrogen bond interaction in neat liquid alkanols: a molecular dynamic simulation study

Self-diffusion of methanol, ethanol, 1-propanol and 2-propanol has been studied by molecular dynamics simulation in the temperature range between the melting pressure curve and 478 K at pressures up to 300 MPa. The simulation results on self-diffusion of methanol, ethanol and 2-propanol (for 2-propanol, at high temperatures) agree well with experiment, which suggests that the simulation method is a powerful tool to obtain self-diffusion coefficients over wide range of temperature and pressure, under which it is rather difficult for experiments. The local structures of methanol, ethanol and 2-propanol are investigated by calculating the radial distribution functions, H-bond numbers, coordination numbers and the ratios of H-bond number divided by coordination number. The correlation between self-diffusion and structural properties, and the influence of temperature and pressure on them are discussed. The degree of forming H-bond space network in methanol, ethanol and water is higher than that in 2-propanol, and they are all higher than those in ammonia and methylamine. The simulation results demonstrate that the effect of hydrogen bonding on the translational dynamics in methanol and ethanol is more pronounced than that in 2-propanol.     

Reconciliation of different simulation methods in the determination of the equilibrium branch for adsorption in pores

The recently proposed mid-density scheme [Liu Z, Herrera L, Nguyen VT, Do DD, Nicholson D. A Monte Carlo scheme based on mid-density in a hysteresis loop to determine equilibrium phase transition. Mol Simul. 2011; 37(11):932–939, Liu Z, Do DD, Nicholson D. A thermodynamic study of the mid-density scheme to determine the equilibrium phase transition in cylindrical pores. Mol Simul. 2012; 38(3):189–199] is tested against a method 2V-NVT (similar to the well-established gauge cell method) and the canonical ensemble (CE) method, using argon adsorption at 87 K in graphitic slit pores of infinite and finite length. In infinitely long pores, the equilibrium transition is vertical that is expected for an infinite system to have a first-order transition and this vertical transition was found to lie at the middle of the hysteresis loop and satisfies the well-known Maxwell rule of equal area. For pores of finite length, the equilibrium transitions are steep and are close to, but not exactly identical to, the desorption branch. This lends support to the conventional view that the desorption branch is nearest to equilibrium, although both adsorption and desorption branches are strictly speaking metastable; a view proposed originally by Everett [Everett DH. Capillary condensation and adsorption hysteresis. Berichte Der Bunsen-Gesellschaft [Phys Chem Chem Phys]. 1975; 79(9):732–734]. As a consequence, the Maxwell rule of equal area does not apply to finite systems. As the widely accepted CE and gauge cell methods do not falsify the mid-density scheme, this study lends strong support to the validity of this technique for the study of equilibria.     

基于生命表参数的棉铃虫种群动态模拟研究

依据对棉铃虫自然种群生命表研究所获得的种群统计参数,采用蒙特卡罗方法,模拟了随机环境条件下的棉铃虫种群发生动态。模拟结果显示,同一输入参数(模拟的1代残虫量),经过确定世代或生活史阶段后,棉铃虫的种群密度既可能为轻发生,也可能种群暴发。环境条件有利时,输入变量(模拟的1代残虫量)影响棉铃虫最大可能的发生程度。当棉铃虫为中等或轻发生时,1代残虫量(模型输入参数)和最终的棉铃虫发生程度(模型输出)没有明显的依赖关系。相同的最终棉铃虫发生程度所对应的模型的输入变量可以有多种水平。模拟结果显示,在环境条件比较有利情况下,棉铃虫只需要1个世代或某一个生活史阶段(如蛹期),种群密度即可达暴发水平。本文的研究结果说明,对棉铃虫中长期发生态势的预测,存在一定程度的不确定性。     

Effect of chirality,length and diameter of carbon nanotubes on the adsorption of 20 amino acids: a molecular dynamics simulation study

We carried out molecular dynamics simulations to study the adsorption of all the 20 amino acids (AAs; aromatic, polar and non-polar) on the surface of chiral, zigzag and armchair single-walled carbon nanotubes. The adsorption was occurring in all systems. In the aromatic AAs, the π–π stacking and the semi-hydrogen bond formation cause a strong interaction with the carbon nanotubes (CNTs). We also investigated the chirality, length and diameter dependencies on adsorption energies. We found that all AAs have more tendency to adsorption on the chiral and zigzag CNTs over the armchair. The results show that increasing both the diameter and the length causes the enhancement of the adsorption energy. But, the effect of the length is more than of the diameter. For example, the adsorption energy of Trp on the surface of CNT (4,1), with 2 nm length, is 20.4 kcal/mol. When the length of CNT becomes twice, the adsorption energy increases by 24 ± 0.3%. But by doubling the diameter, the adsorption energy increased only by 9.8 ± 0.25%.     

Monte Carlo simulation of damascene electroplating: effects of additives

The influence of additives on the filling process of via holes in damascene electroplating is investigated with use of a kinetic Monte Carlo method. The basic system is the solid-by-solid model for crystal growth which includes the vacancy formation during the growth of thin film. Three kinds of additives are included in the model to control the local surface growth rate. Inhibitors and levelers have the effect of preventing the deposition, while accelerators increase the local growth rate. Levelers are modeled to stick to the tips of the surface. We performed a series of simulations by changing the parameters which characterize the additives to see their influence on the filling mechanism. It is shown that void-free filling is possible by the combination of the effects of the additives.     

Effect of hydrogen coverage on the buckling of penta-graphene by molecular dynamics simulation

By employing a series of MD simulations, buckling behaviour of penta-graphene and functionalised penta-graphene having different hydrogen (H) coverage is presented in this study. To this end, the buckling onset strain is determined for different systems. The results reveal that the new allotrope is slightly stiffer than graphene. Moreover, the effect of H adatoms in the range 0–56% on buckling behaviour is investigated. Finally it is shown that the H coverage has an influence on stability, and ripple-type destortion of the penta-graphene under compression.     

Evaluation of three scatter correction methods based on estimation of photopeak scatter spectrum in SPECT imaging: A simulation study

Three practical methods for scatter correction of Tc-99m SPECT images are evaluated. Among these, two methods, three-energy window (TEW) methods using the trapezoidal and triangular approximations, have been described previously by investigators, and a new approximation is offered in this work. The SIMIND (SIMulation of Imaging Nuclear Detectors) Monte Carlo program is used to simulate a line source placed at on-axis and 5 cm off-axis locations, a cold-sphere/hot-background phantom, a hot-sphere/cold-background phantom, and a more clinically realistic NCAT (Nonuniform Rational B-spline-based CArdiac-Torso) phantom. For evaluation of these methods, the scatter line-spread functions and scatter fractions for the on- and off-axis line source, image contrast, signal-to-noise ratio and relative noise for the cold spheres, and recovery coefficient for the hot spheres of different diameters are compared. For the NCAT phantom, a line profile through a slice of the reconstructed image is considered before and after scatter correction, and also image contrast defined by this profile is used to compare the correction methods. The results of this study indicate that for the line source simulation the scatter fractions obtained from the proposed method are a better estimation of true scatter fractions. Also, for both the sphere simulation and NCAT simulation, the proposed method improves the image contrast as compared to the two other methods.