Ali Poursaeidesfahani

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.

Computation of thermodynamic properties in the continuous fractional component Monte Carlo Gibbs ensemble

It is shown that ensemble averages computed in the Gibbs Ensemble with Continuous Fractional Component Monte Carlo (CFCMC GE) are different from those computed in the conventional Gibbs Ensemble (GE). However, it is possible to compute averages corresponding to the conventional GE while performing simulations in the CFCMC GE. In this way, one can benefit from the nice features of CFCMC GE (e.g. more efficient particle exchange) and at the same time compute the ensemble averages that correspond to the conventional GE. As a case study, the equilibrium pressure and densities of the systems of 256 and 512 LJ particles at different reduced temperatures () are computed in the conventional GE and CFCMC GE. The validity of the expressions derived for computation of the thermodynamic pressure and densities corresponding to the conventional GE and computed in the CFCMC GE is examined numerically. The thermodynamic pressure in the conventional GE and CFCMC GE typically differs by at most 4%. It is shown that a very good estimate of the average pressure and densities corresponding to the conventional GE can be obtained by performing simulation in CFCMC GE and ignoring the contributions of the fractional molecule. It is also shown that the fractional molecule does not have an influence on the structure of the liquid, even for very small system sizes (e.g. 40 particles). The approach used here to compute the equilibrium pressure and densities of the conventional GE using the CFCMC GE can be easily extended to other thermodynamic properties and other ensembles.

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.

Chemical potentials of water, methanol, carbon dioxide and hydrogen sulphide at low temperatures using continuous fractional component Gibbs ensemble Monte Carlo

Chemical potentials of coexisting gas and liquid phases for water, methanol, hydrogen sulphide and carbon dioxide for the temperature range K to K are computed using two different methodologies: (1) Widom’s test particle insertion (WTPI) method in the conventional Gibbs Ensemble (GE), and (2) the Continuous Fractional Component Gibbs Ensemble Monte Carlo (CFCGE MC) method. It is shown that the WTPI method fails to accurately compute the chemical potentials of water and methanol in the liquid phase at low temperatures, while accurate chemical potentials in the liquid phase are computed using CFCGE MC method. For the CFCGE MC method, the statistical uncertainty for computed chemical potentials of water and methanol in the liquid phase are considerably smaller compared to the WTPI method. For the water models considered in this study (SPC, TIP3P-EW, TIP4P-EW, TIP5P-EW), computed excess chemical potentials based on three-site models are in better agreement with the chemical potentials computed from an empirical equation of state from the NIST database. For water, orientational biasing is applied during test particle insertion to check whether certain orientations of test particle are energetically unfavourable. A two-dimensional Overlapping Distribution Method (ODM) in the NVT ensemble is derived for this purpose. It is shown that failure of the WTPI method for systems with a strong hydrogen bonding network does not depend on orientation of the test molecule in that system. For all systems in this study, the WTPI method breaks down when the void fraction of the system drops below approximately 0.50.

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.

Improving Olefin Purification Using Metal Organic Frameworks with Open Metal Sites

The separation and purification of light hydrocarbons is challenging in the industry. Recently, a ZJNU-30 metal-organic framework (MOF) has been found to have the potential for adsorption-based separation of olefins and diolefins with four carbon atoms [H. M. Liu et al. Chem.-Eur. J. 2016, 22, 14988-14997]. Our study corroborates this finding but reveals Fe-MOF-74 as a more efficient candidate for the separation because of the open metal sites. We performed adsorption-based separation, transient breakthrough curves, and density functional theory calculations. This combination of techniques provides an extensive understanding of the studied system. Using this MOF, we propose a separation scheme to obtain a high-purity product.