Influence of Nanoscale Intimacy and Zeolite Micropore Size on the Performance of Bifunctional Catalysts for n-Heptane Hydroisomerization.

In this study, Pt nanoparticles on zeolite/γ-Al2O3 composites (50/50 wt) were located either in the zeolite or on the γ-Al2O3 binder, hereby varying the average distance (intimacy) between zeolite acid sites and metal sites from "closest" to "nanoscale". The catalytic performance of these catalysts was compared to physical mixtures of zeolite and Pt/γ-Al2O3 powders, which provide a "microscale" distance between sites. Several beneficial effects on catalytic activity and selectivity for n-heptane hydroisomerization were observed when Pt nanoparticles are located on the γ-Al2O3 binder in nanoscale proximity with zeolite acid sites, as opposed to Pt nanoparticles located inside zeolite crystals. On ZSM-5-based catalysts, mostly monobranched isomers were produced, and the isomer selectivity of these catalysts was almost unaffected with an intimacy ranging from closest to microscale, which can be attributed to the high diffusional barriers of branched isomers within ZSM-5 micropores. For composite catalysts based on large-pore zeolites (zeolite Beta and zeolite Y), the activity and selectivity benefitted from the nanoscale intimacy with Pt, compared to both the closest and microscale intimacies. Intracrystalline gradients of heptenes as reaction intermediates are likely contributors to differences in activity and selectivity. This paper aims to provide insights into the influence of the metal-acid intimacy in bifunctional catalysts based on zeolites with different framework topologies.

Brick-CFCMC: Open Source Software for Monte Carlo Simulations of Phase and Reaction Equilibria Using the Continuous Fractional Component Method.

We present a new molecular simulation code, Brick-CFCMC, for performing Monte Carlo simulations using state-of-the-art simulation techniques. The Continuous Fractional Component (CFC) method is implemented for simulations in the NVT/NPT ensembles, the Gibbs Ensemble, the Grand-Canonical Ensemble, and the Reaction Ensemble. Molecule transfers are facilitated by the use of fractional molecules which significantly improve the efficiency of the simulations. With the CFC method, one can obtain phase equilibria and properties such as chemical potentials and partial molar enthalpies/volumes directly from a single simulation. It is possible to combine trial moves from different ensembles. This enables simulations of phase equilibria in a system where also a chemical reaction takes place. We demonstrate the applicability of our software by investigating the esterification of methanol with acetic acid in a two-phase system.

Hydride Transfer versus Deprotonation Kinetics in the Isobutane-Propene Alkylation Reaction: A Computational Study.

The alkylation of isobutane with light alkenes plays an essential role in modern petrochemical processes for the production of high-octane gasoline. In this study we have employed periodic DFT calculations combined with microkinetic simulations to investigate the complex reaction mechanism of isobutane-propene alkylation catalyzed by zeolitic solid acids. Particular emphasis was given to addressing the selectivity of the alkylate formation versus alkene formation, which requires a high rate of hydride transfer in comparison to the competitive oligomerization and deprotonation reactions resulting in catalyst deactivation. Our calculations reveal that hydride transfer from isobutane to a carbenium ion occurs via a concerted C-C bond formation between a tert-butyl fragment and an additional olefin, or via deprotonation of the tert-butyl fragment to generate isobutene. A combination of high isobutane concentration and low propene concentration at the reaction center favor the selective alkylation. The key reaction step that has to be suppressed to increase the catalyst lifetime is the deprotonation of carbenium intermediates that are part of the hydride transfer reaction cycle.

Efficient Application of Continuous Fractional Component Monte Carlo in the Reaction Ensemble.

A new formulation of the Reaction Ensemble Monte Carlo technique (RxMC) combined with the Continuous Fractional Component Monte Carlo method is presented. This method is denoted by serial Rx/CFC. The key ingredient is that fractional molecules of either reactants or reaction products are present and that chemical reactions always involve fractional molecules. Serial Rx/CFC has the following advantages compared to other approaches: (1) One directly obtains chemical potentials of all reactants and reaction products. Obtained chemical potentials can be used directly as an independent check to ensure that chemical equilibrium is achieved. (2) Independent biasing is applied to the fractional molecules of reactants and reaction products. Therefore, the efficiency of the algorithm is significantly increased, compared to the other approaches. (3) Changes in the maximum scaling parameter of intermolecular interactions can be chosen differently for reactants and reaction products. (4) The number of fractional molecules is reduced. As a proof of principle, our method is tested for Lennard-Jones systems at various pressures and for various chemical reactions. Excellent agreement was found both for average densities and equilibrium mixture compositions computed using serial Rx/CFC, RxMC/CFCMC previously introduced by Rosch and Maginn (Journal of Chemical Theory and Computation, 2011, 7, 269-279), and the conventional RxMC approach. The serial Rx/CFC approach is also tested for the reaction of ammonia synthesis at various temperatures and pressures. Excellent agreement was found between results obtained from serial Rx/CFC, experimental results from literature, and thermodynamic modeling using the Peng-Robinson equation of state. The efficiency of reaction trial moves is improved by a factor of 2 to 3 (depending on the system) compared to the RxMC/CFCMC formulation by Rosch and Maginn.

Behavior of the Enthalpy of Adsorption in Nanoporous Materials Close to Saturation Conditions.

Many important industrial separation processes based on adsorption operate close to saturation. In this regime, the underlying adsorption processes are mostly driven by entropic forces. At equilibrium, the entropy of adsorption is closely related to the enthalpy of adsorption. Thus, studying the behavior of the enthalpy of adsorption as a function of loading is fundamental to understanding separation processes. Unfortunately, close to saturation, the enthalpy of adsorption is hard to measure experimentally and hard to compute in simulations. In simulations, the enthalpy of adsorption is usually obtained from energy/particle fluctuations in the grand-canonical ensemble, but this methodology is hampered by vanishing insertions/deletions at high loading. To investigate the fundamental behavior of the enthalpy and entropy of adsorption at high loading, we develop a simplistic model of adsorption in a channel and show that at saturation the enthalpy of adsorption diverges to large positive values due to repulsive intermolecular interactions. However, there are many systems that can avoid repulsive intermolecular interactions and hence do not show this drastic increase in enthalpy of adsorption close to saturation. We find that the conventional grand-canonical Monte Carlo method is incapable of determining the enthalpy of adsorption from energy/particle fluctuations at high loading. Here, we show that by using the continuous fractional component Monte Carlo, the enthalpy of adsorption close to saturation conditions can be reliably obtained from the energy/particle fluctuations in the grand-canonical ensemble. The best method to study properties at saturation is the NVT energy (local-) slope methodology.

Direct Free Energy Calculation in the Continuous Fractional Component Gibbs Ensemble.

A new formulation of the Gibbs ensemble (GE) combined with the continuous fractional component Monte Carlo method is presented. In the proposed formulation, only a single fractional molecule per component is used instead of two in the original formulation by Shi and Maginn ( J. Comput. Chem. 2008 , 29 , 2520 - 2530 ). This has the following advantages: (1) one directly obtains chemical potentials, without using test particles. We show analytically that the expressions for the chemical potential are identical to those in the conventional Gibbs ensemble; (2) biasing is applied to each simulation box independently; (3) maximum allowed changes in the scaling parameter of intermolecular interactions can be chosen differently in each simulation box. Obtaining chemical potentials directly facilitates thermodynamic modeling using equations of state, and it can be used as an independent check to ensure that chemical equilibrium is achieved. As a proof of principle, our method is tested for Lennard-Jones (LJ) particles and the TIP3P-Ew water model. Results are compared with the conventional GE. Excellent agreement was found both for average densities and chemical potentials. In our new approach, the acceptance probability for molecule exchanges between the boxes is much higher (typically larger than 40% for LJ particles) than for the conventional GE (typically lower than 2% for LJ particles). It is also shown that the contribution of the fractional molecule should be disregarded when computing ensemble averages such as the average energy per molecule and the average densities. The algorithm can be easily extended to mixtures and molecules with intramolecular interactions.