Reactive Monte Carlo and grand-canonical Monte Carlo simulations of the propene metathesis reaction system

The influence of silicalite-1 pores on the reaction equilibria and the selectivity of the propene metathesis reaction system in the temperature range between 300 and 600 K and the pressure range from 0.5 to 7 bars has been investigated with molecular simulations. The reactive Monte Carlo (RxMC) technique was applied for bulk-phase simulations in the isobaric-isothermal ensemble and for two phase systems in the Gibbs ensemble. Additionally, Monte Carlo simulations in the grand-canonical ensemble (GCMC) have been carried out with and without using the RxMC technique. The various simulation procedures were combined with the configurational-bias Monte Carlo approach. It was found that the GCMC simulations are superior to the Gibbs ensemble simulations for reactions where the bulk-phase equilibrium can be calculated in advance and does not have to be simulated simultaneously with the molecules inside the pore. The confined environment can increase the conversion significantly. A large change in selectivity between the bulk phase and the pore phase is observed. Pressure and temperature have strong influences on both conversion and selectivity. At low pressure and temperature both conversion and selectivity have the highest values. The effect of confinement decreases as the temperature increases.     

Flexible polyelectrolyte simulations at the Poisson-Boltzmann level: a comparison of the kink-jump and multigrid configurational-bias Monte Carlo methods

We present a new approach for simulating the motions of flexible polyelectrolyte chains based on the continuous kink-jump Monte Carlo technique coupled to a lattice field theory based calculation of the Poisson-Boltzmann (PB) electrostatic free energy "on the fly." This approach is compared to the configurational-bias Monte Carlo technique, in which the chains are grown on a lattice and the PB equation is solved for each configuration with a linear scaling multigrid method to obtain the many-body free energy. The two approaches are used to calculate end-to-end distances of charged polymer chains in solutions with varying ionic strengths and give similar numerical results. The configurational-bias Monte Carlo/multigrid PB method is found to be more efficient, while the kink-jump Monte Carlo method shows potential utility for simulating nonequilibrium polyelectrolyte dynamics.     

Application of a molybdenum oxide catalyst supported on Al-pillared montmorillonite in the propene metathesis reaction

The catalytic behavior of a molybdenum oxide catalyst supported on Al-pillared montmorillonite in the propene metathesis reaction has been investigated. Comparison with two reference catalysts (Mo/SiO₂ and Mo/Al₂O₃) was also studied. Pillared clay catalysts show an important activity loss and activity recovery by air regeneration was obtained.     

Motion of a branched polymer chain in confinement: a Monte Carlo study

The aim of the study was a theoretical investigation of the polymer molecules located between two parallel and impenetrable surfaces which were also attractive for polymer segments. The chains were constructed of identical segments and were restricted to knots of a simple cubic lattice. Since the chains were at good solvent conditions the only interactions between the segments of the chain were the excluded volume. The properties of the model chains were determined by means of Monte Carlo simulations with a sampling algorithm based on the chain's local changes of conformation. The differences and similarities in the structure for different adsorption regimes and the size of the slit were shown and discussed. It was observed that at certain conditions the polymer chain was adsorbed at one of the confining surfaces, and then after a certain period of time it detached from this surface and approached the opposite wall; this switch was repeated many times. The influence of the strength of the adsorption, the size of the slit, and the chain's length on the frequency of these jumps were determined. The mechanism of the chain's motion during the switch was also shown.     

Adsorption of benzene and propene in zeolite MCM-22: a grand canonical Monte Carlo study

The GCMC (grand canonical Monte Carlo) simulation technique was used to predict the competition adsorption characteristics of benzene and propene in different pore systems of MCM-22. The nine-site model of benzene was used, which proved to be effective and efficient. The zeolite was divided into three adsorption sites following a simulated annealing method. It is found that benzene and propene have the same preferential adsorption site and a similar adsorption order in different sites. Moreover, the pure and mixture isotherms of the three sites are drawn. From the isotherms, we obtained a selectivity reversal of the mixture isotherms of benzene and propene in different sites. It is also noted that the competition adsorption in the three adsorption sites for the two adsorbates can fall into three successive steps and the adsorption order of propene in mixture in these three sites is S3→S1→S2. A new model is presented to predict the benzene and propene adsorption equilibrium in MCM-22. This approach yields better multicomponent equilibrium predictions than ideal adsorbed solution theory (IAST). Isotherms at different mole fraction of benzene in gas phase indicate an advantage to increase the feed radio of benzene and propene. Thus, this work is helpful for a better understanding of the adsorption mechanism of benzene and propene in MCM-22 and hence the relation of the catalytic properties of the zeolite to its structure.     

Permeation and gating in proteins: kinetic Monte Carlo reaction path following

We present a new Monte Carlo technique, kinetic Monte Carlo reaction path following (kMCRPF), for the computer simulation of permeation and large-scale gating transitions in protein channels. It combines ideas from Metropolis Monte Carlo (MMC) and kinetic Monte Carlo (kMC) algorithms, and is particularly suitable when a reaction coordinate is well defined. Evolution of transition proceeds on the reaction coordinate by small jumps (kMC technique) toward the nearest lowest-energy uphill or downhill states, with the jumps thermally activated (constrained MMC). This approach permits navigation among potential minima on an energy surface, finding the minimum-energy paths and determining their associated free-energy profiles. The methodological and algorithmic strategies underlying the kMCRPF method are described. We have tested it using an analytical model and applied it to study permeation through the curvilinear ClC chloride and aquaporin pores and to gating in the gramicidin A channel. These studies of permeation and gating in real proteins provide extensive procedural tests of the method.     

Gibbs ensemble Monte Carlo simulation of supercritical CO2 adsorption on NaA and NaX zeolites

Adsorption of supercritical carbon dioxide on two kinds of zeolites with identical chemical composition but different pore structure (NaA and NaX) was studied using the Gibbs ensemble Monte Carlo simulation. The model frameworks for the two zeolites with SiAl ratio being unity have been chosen as the solid structures in the simulation. The adsorption behaviors of supercritical CO2 on the NaA and NaX zeolites, based on the adsorption isotherms and isosteric heats of adsorption, were discussed in detail and were compared with the available experimental results. A good agreement between the simulated and experimental results is obtained for both the adsorbed amount and the bulk phase density. The intermediate configurational snapshots and the radial distribution functions between zeolite and adsorbed CO2 molecules were collected in order to investigate the preferable adsorption locations and the confined structure behavior of CO2. The structure behaviors of the adsorbed CO2 molecules show various performances, as compared with the bulk phase, due to the confined effect in the zeolite pores.     

Adsorption and separation of carbon dioxide and methane in new zeolites using the Grand Canonical Monte Carlo method

The adsorption and separation behaviors of CO₂ and CH₄ in new siliceous zeolites (IFO, JSR, OKO, SEW) were simulated using the Grand Canonical Monte Carlo method in the paper. The adsorption isotherms for pure components and binary mixtures of CO₂ and CH₄ in four siliceous zeolites were obtained. The adsorption thermodynamic properties including Gibb’s free energy change, enthalpy change and entropy change were investigated. The results demonstrate that the adsorbed amount of pure components increases with an increase in pressure, and larger pore volume and surface area are beneficial to improve the adsorption capacity. The adsorption amount of CO₂ and CH₄ in the JSR zeolite is 7.08 and 2.27 mmol g^?1 at 1000 kPa, respectively. In view of the thermodynamic results, the new siliceous zeolites show a higher affinity for CO₂. The adsorption capacities of CO₂ in all zeolites were five times more than those of CH₄ in binary mixtures based on the ratios of equilibrium adsorption capacity. Considering the adsorption uptake and selectivity for CO₂/CH₄, the JSR zeolite is a good candidate for the separation of CO₂/CH₄ at low pressure.     

Monte Carlo simulation of the influence of solvent on nucleic acid base associations

The results of Monte Carlo calculations of the association between nucleic acid bases in a nonpolar solvent (CCl4) are described. The influence of the solvent on planar and stacked associations of bases was examined by analyzing the total energy of the system, including solute-solute, solute-solvent, and solvent-solvent contributions. Good quantitative agreement with the available experimental data was obtained. Solute-solvent interactions are primarily determined by dispersion forces; consequently, solute-solvent interactions vertical to the solute plane that maximize dispersion interactions are most favored, and a rough proportionally between solute-solvent energy and the surface of the solute was observed. Analysis of solvent-solvent energy is not necessarily reduced when surface area decreases, contrary to the simple cavity concept. "Single molecule probe" calculations were performed to explain the differences in base associations in H2O and CCl4. In CCl4 dispersion forces dominate and planar complexes are stabilized by maximum exposure of molecular planes to the solvent. In H2O electrostatic forces dominate so that the most stable structures are stacked association that allow the maximum number of hydrophilic centers to be exposed to the solvent.     

The influence of side chains on the properties of simple model biopolymers: Monte Carlo simulations

A model for the simulation of proteins is introduced which is based on a new set of bond vectors and a new method for modeling the side chains of proteins. The drawbacks of united atoms models are summarized and the motivation for this new model is given. Some preliminary results are shown which shall demonstrate the suitability of the model proposed.     

Monte Carlo simulation of polymer analogous reaction in confined conditions: effects of ordering

The influence of energetic parameters of the interchain homo- and heterocontacts on a local ordering of Bernoullian copolymers has been studied using Monte Carlo simulations and probabilistic analysis. The results of both methods are in a good agreement. Then simple Monte Carlo procedure was employed to study the ordering in products of a polymeranalogous reaction with accelerating effect of neighboring groups. When the reaction with intra- and interchain acceleration and local ordering proceed simultaneously in confined conditions, the ordering might affect the process so that the formation of certain nano-structures (in particular, not trivial strip-like ones) is possible.     

Monte Carlo simulation of adsorption of binary and quaternary alkane isomers mixtures in zeolites: Effect of pore size and structure

Grand canonical Monte Carlo and configurational bias Monte Carlo techniques were employed to simulate the adsorption of binary mixtures of butane isomers and quaternary mixtures in nine zeolites at 300 K. For binary mixtures the results show there is a critical pore size, which is 10-membered-ring about 5.6 Å. The channel sizes of BEA, ISV, MOR and CFI are larger than this critical pore size, they prefer i-butane than n-butane, whereas TON with smaller channel size than critical pore size prefers n-butane than i-butane, but its selectivity decreases with pressure increasing. MFI, MEL and TER prefer i-butane than n-butane at low pressure, but with pressure increasing, the selectivity is reversed. BOG prefers i-butane than n-butane but the selectivity decreased with pressure increasing. It demonstrates that the adsorption and selectivity are controlled by both pore size and pore structure. The n-butane–i-butane–n-pentane–2-methylbutane quaternary mixtures adsorbed in these nine zeolites were studied, and the results show alkane chain length dependence at low pressure, but the adsorption is controlled by pore size and structure with pressure increasing in all the zeolites except for TON and BOG.     

Effects of confinement on the adsorption behavior of methane in high-silica zeolites

Adsorption isotherms of methane on high-silica zeolites at ambient temperature and up to high pressure were experimentally measured. The isotherms were analyzed on the basis of thermodynamics and statistical mechanics, and the relationship between microscopic properties and the macroscopic adsorption behavior was investigated. A comparison between the confined and unconfined phases revealed that molecular motion is restricted in the pores. As a result, the adsorbed phase is entropically destabilized, which cannot be neglected in comparison with the energetical stabilization that occurs as a result of the solid-molecule interaction. Our findings also indicate that the smaller slope (drho/d ln p) of the adsorption isotherms compared to that of the isotherm of the bulk at the same density is due to the smaller intermolecular interaction in the pores.The pore-size dependence is indicated not only in solid-molecule interactions but also in intermolecular interactions and molecular motion. Of these, the solid-molecule interactions strongly influence the adsorption behaviors in pores of different sizes. The origin of the restriction of molecular motion in the pores is well-explained by the one-dimensional transition and two-dimensional vibration (1D-trans, 2D-vib) model.     

The influence of correlation on the interpretation of Hund's multiplicity rule: A quantum Monte Carlo study

A systematic quantum Monte Carlo study of 2p atoms (C, N, O) and 3p atoms (Si, P, S) is performed to investigate the influence of correlation on the interpretation of Hund's multiplicity rule, which is an extension of our previous study of the carbon atom [J. Chem. Phys. 121, 7144 (2004)] to heavier atoms. The accuracy in the present study is significantly improved as compared with the previous study. A detailed analysis of the correlation contribution to individual energy components of the total energy is given beyond the self-consistent Hartree-Fock calculation. The stability of the highest spin-multiplicity state of all the atoms is ascribed to the greater electron-nucleus attraction energy that is gained at the cost of increasing the electron-electron repulsion energy as well as the kinetic energy. The present study demonstrates that correlation does not change the above conclusion due to the Hartree-Fock theory to support Boyd's less screening mechanism.     

Isomer-specific influences on the composition of reaction intermediates in dimethyl ether/propene and ethanol/propene flame

This work provides experimental evidence on how the molecular compositions of fuel-rich low-pressure premixed flames are influenced as the oxygenates dimethyl ether (DME) or ethanol are incrementally blended into the propene fuel. Ten different flames with a carbon-to-oxygen ratio of 0.5, ranging from 100% propene (phi = 1.5) to 100% oxygenated fuel (phi = 2.0), are analyzed with flame-sampling molecular-beam mass spectrometry employing electron- or photoionization. Absolute mole fraction profiles for flame species with masses ranging from m/z = 2 (H2) to m/z = 80 (C6H8) are analyzed with particular emphasis on the formation of harmful emissions. Fuel-specific destruction pathways, likely to be initiated by hydrogen abstraction, appear to lead to benzene from propene combustion and to formaldehyde and acetaldehyde through DME and ethanol combustion, respectively. While the concentration of acetaldehyde increases 10-fold as propene is substituted by ethanol, it decreases as propene is replaced with DME. In contrast, the formaldehyde concentration rises only slightly with ethanol replacement but increases markedly with addition of DME. Allyl and propargyl radicals, the dominant precursors for benzene formation, are likely to be produced directly from propene decomposition or via allene and propyne. Benzene formation through propargyl radicals formed via unsaturated C2 intermediates in the decomposition of DME and ethanol is negligibly small. As a consequence, DME and ethanol addition lead to similar reductions of the benzene concentration.     

Monte Carlo simulation of the interchain exchange reaction in a blend of incompatible polymers

The interchain exchange reaction in a blend composed of two contacting layers of incompatible A and B homopolymers was simulated by means of the dynamic off-lattice Monte Carlo method. The evolution of local molecular-mass and block mass distributions, depending on the effective temperature and the reaction rate, was studied for the first time. It was shown that the components interpenetrate as the copolymer forms in the interphase layer and the average block length decreases below a certain, temperature-dependent value. The state of dynamic equilibrium, whose characteristics are determined mainly by temperature, is established in the system. The time of establishment of equilibrium and the intensity of compatibilization at the early steps of the process are controlled by the rate of the reaction. The results of the study allow the contribution of the reaction to the interchange processes to be evaluated.     

Dynamic Monte Carlo simulation of the NO + CO reaction on Rh(111)

The kinetics of the catalytic reduction of NO by CO on Rh(111) surfaces was investigated by using dynamic Monte Carlo simulations. Our model takes into account recent experimental findings and introduces relevant modifications to the classical reaction scheme, including an alternative pathway to produce N2 through an (N-NO)* intermediate, the formation of atomic nitrogen islands in the adsorbed phase, and the influence of coadsorbed species on the dissociation of NO. All elementary steps are expressed as activated processes with temperature-dependent rates and realistic values dictated by experiments. Calculated steady-state phase diagrams are presented for the NO + CO reaction showing the windows for the conditions under which the reaction is viable. The model predicts variations in both production rates and adsorbate coverages with temperature consistent with experimental data. The effect of varying the individual kinetic parameters and the importance of each step in the reaction scheme were probed.     

Insight into the topology effect on the diffusion of ethene and propene in zeolites: A molecular dynamics simulation study

Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of products.Besides the topology effect on reaction kinetics,the topology influences the diffusion of reactants and products in catalysts as well.In this work,by using COMPASS force-field molecular dynamics method,we investigated the intracrystalline diffusion of ethene and propene in four different zeolites,CHA,MFI,BEA and FAU,at different temperatures.The self-diffusion coefficients and diffusion activation barriers were calculated.A strong restriction on the diffusion of propene in CHA was observed because the self-diffusion coefficient ratio of ethene to propene is larger than 18 and the diffusion activation barrier of propene is more than 20 kJ/mol in CHA.This ratio decreases with the increase of temperature in the four investigated zeolites.The shape selectivity on products from diffusion perspective can provide some implications on the understanding of the selectivity difference between HSAPO-34 and HZSM-5 catalysts for the MTO conversion.     

Monte Carlo simulation for inhomogeneous chemical kinetics: Application to the Belousov–Zhabotinskii reaction

A Monte Carlo algorithm, capable of simulating numerically the time and space dependence of chemical concentrations in a reacting system, is presented. This method is used to study the phenomenon of trigger waves in the Oregonator model of the Belousov–Zhabotinskii reaction, including the diffusion of species X and Y in one dimension. The results show that a small disturbance in a homogeneous mixture can grow into a chemical (trigger) wave propagating in space at constant velocity. The dependence of this velocity on several factors is studied, namely, initial concentrations, the diffusion of Y, and the stoichiometry of the autocatalytic step of the model. A comparison of the Monte Carlo results with a previous simulation also is discussed.