1
|
van Rooijen WA, Habibi P, Xu K, Dey P, Vlugt TJH, Hajibeygi H, Moultos OA. Interfacial Tensions, Solubilities, and Transport Properties of the H 2/H 2O/NaCl System: A Molecular Simulation Study. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2024; 69:307-319. [PMID: 38352074 PMCID: PMC10859954 DOI: 10.1021/acs.jced.2c00707] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Accepted: 12/23/2022] [Indexed: 02/16/2024]
Abstract
Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H2), such as underground H2 storage and H2O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H2O, (b) self-diffusivities of infinitely diluted H2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O, and (c) solubilities of H2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O. The force fields used are the TIP4P/2005 for H2O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H2. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.
Collapse
Affiliation(s)
- W. A. van Rooijen
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, Stevinweg 1, 2628CN, Delft, The Netherlands
| | - P. Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - K. Xu
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - P. Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - T. J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
| | - H. Hajibeygi
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, Stevinweg 1, 2628CN, Delft, The Netherlands
| | - O. A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
| |
Collapse
|
2
|
Hulikal Chakrapani T, Hajibeygi H, Moultos OA, Vlugt TJH. Calculating Thermodynamic Factors for Diffusion Using the Continuous Fractional Component Monte Carlo Method. J Chem Theory Comput 2024; 20:333-347. [PMID: 38113860 PMCID: PMC10782482 DOI: 10.1021/acs.jctc.3c01144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 11/28/2023] [Accepted: 11/28/2023] [Indexed: 12/21/2023]
Abstract
Thermodynamic factors for diffusion connect the Fick and Maxwell-Stefan diffusion coefficients used to quantify mass transfer. Activity coefficient models or equations of state can be fitted to experimental or simulation data, from which thermodynamic factors can be obtained by differentiation. The accuracy of thermodynamic factors determined using indirect routes is dictated by the specific choice of an activity coefficient model or an equation of state. The Permuted Widom's Test Particle Insertion (PWTPI) method developed by Balaji et al. enables direct determination of thermodynamic factors in binary and multicomponent systems. For highly dense systems, for example, typical liquids, it is well known that molecular test insertion methods fail. In this article, we use the Continuous Fractional Component Monte Carlo (CFCMC) method to directly calculate thermodynamic factors by adopting the PWTPI method. The CFCMC method uses fractional molecules whose interactions with their surrounding molecules are modulated by a coupling parameter. Even in highly dense systems, the CFCMC method efficiently handles molecule insertions and removals, overcoming the limitations of the PWTPI method. We show excellent agreement between the results of the PWTPI and CFCMC methods for the calculation of thermodynamic factors in binary systems of Lennard-Jones molecules and ternary systems of Weeks-Chandler-Andersen molecules. The CFCMC method applied to calculate the thermodynamic factors of realistic molecular systems consisting of binary mixtures of carbon dioxide and hydrogen agrees well with the NIST REFPROP database. Our study highlights the effectiveness of the CFCMC method in determining thermodynamic factors for diffusion, even in densely packed systems, using relatively small numbers of molecules.
Collapse
Affiliation(s)
- Thejas Hulikal Chakrapani
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, 2628 CN Delft, The
Netherlands
| | - Hadi Hajibeygi
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, 2628 CN Delft, The
Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, 2628 CB Delft, The
Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, 2628 CB Delft, The
Netherlands
| |
Collapse
|
3
|
Kobayashi K, Firoozabadi A. Effect of Branching on Mutual Solubility of Alkane-CO 2 Systems by Molecular Simulations. J Phys Chem B 2022; 126:8300-8308. [PMID: 36197719 DOI: 10.1021/acs.jpcb.2c05774] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Mutual solubilities of hydrocarbon-CO2 systems are important in a broad range of applications. Experimental data and theoretical understanding of phase behavior of large hydrocarbon molecules and CO2 are limited. This is especially true in relation to the molecular structure of hydrocarbons when the carbon number exceeds 12. In this work, the continuous fractional component Gibbs ensemble Monte Carlo simulations are used to investigate mutual solubility of different alkane and CO2 systems and the molecular structure. We investigate the mutual solubility of n-decane, n-hexadecane, n-eicosane, and the corresponding structural isomers in the CO2-rich and hydrocarbon-rich phase. The focus will be solubility of the heavy normal alkanes and their structural isomers in CO2. The simulation results are verified by comparing the experimental data when measurements are available. The simulation of phase behavior of the n-decane-CO2 system agrees with the experiments. We also present simulation results of n-hexadecane-CO2 and n-eicosane-CO2 systems away from the critical region partly due to the finite size effect. We establish that solubility of the hydrocarbons in CO2 is improved by change of the molecular structure in heavier alkanes. The enhanced solubility is limited in decane isomers, but the isomers of hexadecane and eicosane show 2- to 3-time solubility enhancement. The molecular dynamics simulations suggest that the improvement is from a higher coordination number of CO2 for methyl (CH3) rather than for methylene (CH2) groups. This study sets the stage for molecular engineering and synthesis of hydrocarbons that are soluble in CO2 not only by considering functionality but also by changing the molecular structure. The solubility enhancement is the first step in viscosification of CO2 which broadens the use of CO2.
Collapse
Affiliation(s)
- Kazuya Kobayashi
- INPEX Corporation, Akasaka Biz Tower 5-3-1 Akasaka, Minato-ku, Tokyo107-6332, Japan.,Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas77005, United States
| | - Abbas Firoozabadi
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas77005, United States
| |
Collapse
|
4
|
Vlugt TJH. Partial molar properties from single molecular dynamics simulations. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2126509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
Affiliation(s)
- Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft, The Netherlands
| |
Collapse
|
5
|
Dawass N, Langeveld J, Ramdin M, Pérez-Gallent E, Villanueva AA, Giling EJM, Langerak J, van den Broeke LJP, Vlugt TJH, Moultos OA. Solubilities and Transport Properties of CO 2, Oxalic Acid, and Formic Acid in Mixed Solvents Composed of Deep Eutectic Solvents, Methanol, and Propylene Carbonate. J Phys Chem B 2022; 126:3572-3584. [PMID: 35507866 PMCID: PMC9125562 DOI: 10.1021/acs.jpcb.2c01425] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
![]()
Recently, deep eutectic
solvents (DES) have been considered as
possible electrolytes for the electrochemical reduction of CO2 to value-added products such as formic and oxalic acids.
The applicability of pure DES as electrolytes is hindered by high
viscosities. Mixtures of DES with organic solvents can be a promising
way of designing superior electrolytes by exploiting the advantages
of each solvent type. In this study, densities, viscosities, diffusivities,
and ionic conductivities of mixed solvents comprising DES (i.e., reline
and ethaline), methanol, and propylene carbonate were computed using
molecular simulations. To provide a quantitative assessment of the
affinity and mass transport of CO2 and oxalic and formic
acids in the mixed solvents, the solubilities and self-diffusivities
of these solutes were also computed. Our results show that the addition
of DES to the organic solvents enhances the solubilities of oxalic
and formic acids, while the solubility of CO2 in the ethaline-containing
mixtures are in the same order of magnitude with the respective pure
organic components. A monotonic increase in the densities and viscosities
of the mixed solvents is observed as the mole fraction of DES in the
mixture increases, with the exception of the density of ethaline-propylene
carbonate which shows the opposite behavior due to the high viscosity
of the pure organic component. The self-diffusivities of all species
in the mixtures significantly decrease as the mole fraction of DES
approaches unity. Similarly, the self-diffusivities of the dissolved
CO2 and the oxalic and formic acids also decrease by at
least 1 order of magnitude as the composition of the mixture shifts
from the pure organic component to pure DES. The computed ionic conductivities
of all mixed solvents show a maximum value for mole fractions of DES
in the range from 0.2 to 0.6 and decrease as more DES is added to
the mixtures. Since for most mixtures studied here no prior experimental
measurements exist, our findings can serve as a first data set based
on which further investigation of DES-containing electrolyte solutions
can be performed for the electrochemical reduction of CO2 to useful chemicals.
Collapse
Affiliation(s)
- Noura Dawass
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Jilles Langeveld
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Mahinder Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Elena Pérez-Gallent
- Department of Sustainable Process and Energy Systems, TNO, Delft, Zuid-Holland 2628CA, The Netherlands
| | - Angel A Villanueva
- Department of Sustainable Process and Energy Systems, TNO, Delft, Zuid-Holland 2628CA, The Netherlands
| | - Erwin J M Giling
- Department of Sustainable Process and Energy Systems, TNO, Delft, Zuid-Holland 2628CA, The Netherlands
| | - Jort Langerak
- Research and Development Department, DMT Environmental Technology, Yndustrywei 3, 8501SN Joure, The Netherlands
| | - Leo J P van den Broeke
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| |
Collapse
|
6
|
Polat HM, Salehi HS, Hens R, Wasik DO, Rahbari A, de Meyer F, Houriez C, Coquelet C, Calero S, Dubbeldam D, Moultos OA, Vlugt TJH. New Features of the Open Source Monte Carlo Software Brick-CFCMC: Thermodynamic Integration and Hybrid Trial Moves. J Chem Inf Model 2021; 61:3752-3757. [PMID: 34383501 PMCID: PMC8385706 DOI: 10.1021/acs.jcim.1c00652] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present several
new major features added to the Monte Carlo
(MC) simulation code Brick-CFCMC for phase- and reaction equilibria
calculations (https://gitlab.com/ETh_TU_Delft/Brick-CFCMC). The first one
is thermodynamic integration for the computation of excess chemical
potentials (μex). For this purpose, we implemented
the computation of the ensemble average of the derivative of the potential
energy with respect to the scaling factor for intermolecular interactions
(). Efficient bookkeeping is implemented
so that the quantity is updated after every MC trial
move with
negligible computational cost. We demonstrate the accuracy and reliability
of the calculation of μex for sodium chloride in
water. Second, we implemented hybrid MC/MD translation and rotation
trial moves to increase the efficiency of sampling of the configuration
space. In these trial moves, short Molecular Dynamics (MD) trajectories
are performed to collectively displace or rotate all molecules in
the system. These trajectories are accepted or rejected based on the
total energy drift. The efficiency of these trial moves can be tuned
by changing the time step and the trajectory length. The new trial
moves are demonstrated using MC simulations of a viscous fluid (deep
eutectic solvent).
Collapse
Affiliation(s)
- H Mert Polat
- CCUS and Acid Gas Entity, Liquefied Natural Gas Department, Exploration Production, TotalEnergies S.E., 92078 Paris, France.,Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands.,CTP - Centre of Thermodynamics of Processes, Mines ParisTech, PSL University, 35 rue Saint Honoré, 77305 Fontainebleau, France
| | - Hirad S Salehi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Remco Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Dominika O Wasik
- Materials Simulation and Modelling, Department of Applied Physics, Eindhoven University of Technology, 5600MB Eindhoven, The Netherlands
| | - Ahmadreza Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Frédérick de Meyer
- CCUS and Acid Gas Entity, Liquefied Natural Gas Department, Exploration Production, TotalEnergies S.E., 92078 Paris, France.,CTP - Centre of Thermodynamics of Processes, Mines ParisTech, PSL University, 35 rue Saint Honoré, 77305 Fontainebleau, France
| | - Céline Houriez
- CTP - Centre of Thermodynamics of Processes, Mines ParisTech, PSL University, 35 rue Saint Honoré, 77305 Fontainebleau, France
| | - Christophe Coquelet
- CTP - Centre of Thermodynamics of Processes, Mines ParisTech, PSL University, 35 rue Saint Honoré, 77305 Fontainebleau, France
| | - Sofia Calero
- Materials Simulation and Modelling, Department of Applied Physics, Eindhoven University of Technology, 5600MB Eindhoven, The Netherlands.,Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, Ctra. Utrera Km. 1, Seville ES-41013, Spain
| | - David Dubbeldam
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, 1098XH Amsterdam, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| |
Collapse
|
7
|
Rahbari A, Garcia-Navarro JC, Ramdin M, van den Broeke LJP, Moultos OA, Dubbeldam D, Vlugt TJH. Effect of Water Content on Thermodynamic Properties of Compressed Hydrogen. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2021; 66:2071-2087. [PMID: 34054140 PMCID: PMC8154567 DOI: 10.1021/acs.jced.1c00020] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Accepted: 03/16/2021] [Indexed: 06/12/2023]
Abstract
Force field-based molecular simulations were used to calculate thermal expansivities, heat capacities, and Joule-Thomson coefficients of binary (standard) hydrogen-water mixtures for temperatures between 366.15 and 423.15 K and pressures between 50 and 1000 bar. The mole fraction of water in saturated hydrogen-water mixtures in the gas phase ranges from 0.004 to 0.138. The same properties were calculated for pure hydrogen at 323.15 K and pressures between 100 and 1000 bar. Simulations were performed using the TIP3P and a modified TIP4P force field for water and the Marx, Vrabec, Cracknell, Buch, and Hirschfelder force fields for hydrogen. The vapor-liquid equilibria of hydrogen-water mixtures were calculated along the melting line of ice Ih, corresponding to temperatures between 264.21 and 272.4 K, using the TIP3P force field for water and the Marx force field for hydrogen. In this temperature range, the solubilities and the chemical potentials of hydrogen and water were obtained. Based on the computed solubility data of hydrogen in water, the freezing-point depression of water was computed ranging from 264.21 to 272.4 K. The modified TIP4P and Marx force fields were used to improve the solubility calculations of hydrogen-water mixtures reported in our previous study [Rahbari A.;J. Chem. Eng. Data2019, 64, 4103-4115] for temperatures between 323 and 423 K and pressures ranging from 100 to 1000 bar. The chemical potentials of ice Ih were calculated as a function of pressure between 100 and 1000 bar, along the melting line for temperatures between 264.21 and 272.4 K, using the IAPWS equation of state for ice Ih. We show that at low pressures, the presence of water has a large effect on the thermodynamic properties of compressed hydrogen. Our conclusions may have consequences for the energetics of a hydrogen refueling station using electrochemical hydrogen compressors.
Collapse
Affiliation(s)
- Ahmadreza Rahbari
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | | | - Mahinder Ramdin
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Leo J. P. van den Broeke
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - David Dubbeldam
- Van’t
Hoff Institute for Molecular Sciences, University
of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| |
Collapse
|
8
|
Caro-Ortiz S, Zuidema E, Rigutto M, Dubbeldam D, Vlugt TJH. Competitive Adsorption of Xylenes at Chemical Equilibrium in Zeolites. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2021; 125:4155-4174. [PMID: 33841605 PMCID: PMC8025683 DOI: 10.1021/acs.jpcc.0c09411] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/17/2020] [Revised: 01/28/2021] [Indexed: 06/12/2023]
Abstract
The separation of xylenes is one of the most important processes in the petrochemical industry. In this article, the competitive adsorption from a fluid-phase mixture of xylenes in zeolites is studied. Adsorption from both vapor and liquid phases is considered. Computations of adsorption of pure xylenes and a mixture of xylenes at chemical equilibrium in several zeolite types at 250 °C are performed by Monte Carlo simulations. It is observed that shape and size selectivity entropic effects are predominant for small one-dimensional systems. Entropic effects due to the efficient arrangement of xylenes become relevant for large one-dimensional systems. For zeolites with two intersecting channels, the selectivity is determined by a competition between enthalpic and entropic effects. Such effects are related to the orientation of the methyl groups of the xylenes. m-Xylene is preferentially adsorbed if xylenes fit tightly in the intersection of the channels. If the intersection is much larger than the adsorbed molecules, p-xylene is preferentially adsorbed. This study provides insight into how the zeolite topology can influence the competitive adsorption and selectivity of xylenes at reaction conditions. Different selectivities are observed when a vapor phase is adsorbed compared to the adsorption from a liquid phase. These insight have a direct impact on the design criteria for future applications of zeolites in the industry. MRE-type and AFI-type zeolites exclusively adsorb p-xylene and o-xylene from the mixture of xylenes in the liquid phase, respectively. These zeolite types show potential to be used as high-performing molecular sieves for xylene separation and catalysis.
Collapse
Affiliation(s)
- Sebastián Caro-Ortiz
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Erik Zuidema
- Shell
Global Solutions International B.V., PO Box 38000, 1030 BN Amsterdam, The Netherlands
| | - Marcello Rigutto
- Shell
Global Solutions International B.V., PO Box 38000, 1030 BN Amsterdam, The Netherlands
| | - David Dubbeldam
- Van’t
Hoff Institute of Molecular Sciences, University
of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
9
|
Dawass N, Wanderley RR, Ramdin M, Moultos OA, Knuutila HK, Vlugt TJH. Solubility of Carbon Dioxide, Hydrogen Sulfide, Methane, and Nitrogen in Monoethylene Glycol; Experiments and Molecular Simulation. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2021; 66:524-534. [PMID: 33487733 PMCID: PMC7818648 DOI: 10.1021/acs.jced.0c00771] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/21/2020] [Accepted: 11/24/2020] [Indexed: 06/12/2023]
Abstract
Knowledge on the solubility of gases, especially carbon dioxide (CO2), in monoethylene glycol (MEG) is relevant for a number of industrial applications such as separation processes and gas hydrate prevention. In this study, the solubility of CO2 in MEG was measured experimentally at temperatures of 333.15, 353.15, and 373.15 K. Experimental data were used to validate Monte Carlo (MC) simulations. Continuous fractional component MC simulations in the osmotic ensemble were performed to compute the solubility of CO2 in MEG at the same temperatures and at pressures up to 10 bar. MC simulations were also used to study the solubility of methane (CH4), hydrogen sulfide (H2S), and nitrogen (N2) in MEG at 373.15 K. Solubilities from experiments and simulations are in good agreement at low pressures, but deviations were observed at high pressures. Henry coefficients were also computed using MC simulations and compared to experimental values. The order of solubilities of the gases in MEG at 373.15 K was computed as H2S > CO2 > CH4 > N2. Force field modifications may be required to improve the prediction of solubilities of gases in MEG at high pressures and low temperatures.
Collapse
Affiliation(s)
- Noura Dawass
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Ricardo R. Wanderley
- Department
of Chemical Engineering, Norwegian University
of Science and Technology, 7034 Trondheim, Norway
| | - Mahinder Ramdin
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Hanna K. Knuutila
- Department
of Chemical Engineering, Norwegian University
of Science and Technology, 7034 Trondheim, Norway
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
10
|
Rahbari A, Hens R, Ramdin M, Moultos OA, Dubbeldam D, Vlugt TJH. Recent advances in the continuous fractional component Monte Carlo methodology. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1828585] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - O. A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - D. Dubbeldam
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| |
Collapse
|
11
|
Werdehausen D, Santiago XG, Burger S, Staude I, Pertsch T, Rockstuhl C, Decker M. Modeling Optical Materials at the Single Scatterer Level: The Transition from Homogeneous to Heterogeneous Materials. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000192] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Daniel Werdehausen
- Corporate Research & Technology Carl Zeiss AG Carl Zeiss Promenade 10 07745 Jena Germany
- Institute of Applied Physics Abbe Center of Photonics Friedrich Schiller University Jena Albert‐Einstein‐Str. 15 07745 Jena Germany
| | - Xavier Garcia Santiago
- JCMwave GmbH Bolivarallee 22 14050 Berlin Germany
- Zuse Institute Berlin Takustr. 7 14195 Berlin Germany
- Institut für Nanotechnology Karlsruher Institut für Technologie PO‐Box 3640 76021 Karlsruhe Germany
| | - Sven Burger
- JCMwave GmbH Bolivarallee 22 14050 Berlin Germany
- Zuse Institute Berlin Takustr. 7 14195 Berlin Germany
| | - Isabelle Staude
- Institute of Applied Physics Abbe Center of Photonics Friedrich Schiller University Jena Albert‐Einstein‐Str. 15 07745 Jena Germany
- Institute for Solid State Physics Friedrich Schiller University Jena Max‐Wien‐Platz 1 07743 Jena Germany
| | - Thomas Pertsch
- Institute of Applied Physics Abbe Center of Photonics Friedrich Schiller University Jena Albert‐Einstein‐Str. 15 07745 Jena Germany
- Fraunhofer Institute for Applied Optics and Precision Engineering Albert‐Einstein‐Str. 7 07745 Jena Germany
- Max Planck School of Photonics Germany
| | - Carsten Rockstuhl
- Institut für Nanotechnology Karlsruher Institut für Technologie PO‐Box 3640 76021 Karlsruhe Germany
- Institut für Theoretische Festkörperphysik Karlsruher Institut für Technologie Wolfgang‐Gaede‐Str. 1 76131 Karlsruhe Germany
- Max Planck School of Photonics Germany
| | - Manuel Decker
- Corporate Research & Technology Carl Zeiss AG Carl Zeiss Promenade 10 07745 Jena Germany
| |
Collapse
|
12
|
Hens R, Rahbari A, Caro-Ortiz S, Dawass N, Erdős M, Poursaeidesfahani A, Salehi HS, Celebi AT, Ramdin M, Moultos OA, Dubbeldam D, Vlugt TJH. Brick-CFCMC: Open Source Software for Monte Carlo Simulations of Phase and Reaction Equilibria Using the Continuous Fractional Component Method. J Chem Inf Model 2020; 60:2678-2682. [PMID: 32275829 PMCID: PMC7312392 DOI: 10.1021/acs.jcim.0c00334] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2020] [Indexed: 12/02/2022]
Abstract
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.
Collapse
Affiliation(s)
- Remco Hens
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Ahmadreza Rahbari
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Sebastián Caro-Ortiz
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Noura Dawass
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Máté Erdős
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Ali Poursaeidesfahani
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Hirad S. Salehi
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Alper T. Celebi
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Mahinder Ramdin
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - David Dubbeldam
- Van’t
Hoff Institute for Molecular Sciences, University
of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| |
Collapse
|
13
|
Li D, Gao Z, Vasudevan NK, Li H, Gao X, Li X, Xi L. Molecular Mechanism for Azeotrope Formation in Ethanol/Benzene Binary Mixtures through Gibbs Ensemble Monte Carlo Simulation. J Phys Chem B 2020; 124:3371-3386. [PMID: 32250637 DOI: 10.1021/acs.jpcb.9b12013] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Azeotropes have been studied for decades due to the challenges they impose on separation processes but fundamental understanding at the molecular level remains limited. Although molecular simulation has demonstrated its capability of predicting mixture vapor-liquid equilibrium (VLE) behaviors, including azeotropes, its potential for mechanistic investigation has not been fully exploited. In this study, we use the united atom transferable potentials for phase equilibria (TraPPE-UA) force field to model the ethanol/benzene mixture, which displays a positive azeotrope. Gibbs ensemble Monte Carlo (GEMC) simulation is performed to predict the VLE phase diagram, including an azeotrope point. The results accurately agree with experimental measurements. We argue that the molecular mechanism of azeotrope formation cannot be fully understood by studying the mixture liquid-state stability at the azeotrope point alone. Rather, azeotrope occurrence is only a reflection of the changing relative volatility between the two components over a much wider composition range. A thermodynamic criterion is thus proposed on the basis of the comparison of partial excess Gibbs energy between the components. In the ethanol/benzene system, molecular energetics shows that with increasing ethanol mole fraction, its volatility initially decreases but later plateaus, while benzene volatility is initially nearly constant and only starts to decrease when its mole fraction is low. Analysis of the mixture liquid structure, including a detailed investigation of ethanol hydrogen-bonding configurations at different composition levels, reveals the underlying molecular mechanism for the changing volatilities responsible for the azeotrope.
Collapse
Affiliation(s)
- Dongyang Li
- School of Chemical Engineering and Technology, National Engineering Research Center of Distillation Technology, and Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin University, Tianjin 300072, China.,Department of Chemical Engineering, McMaster Universtiy, Hamilton, Ontario L8S 4L7, Canada
| | - Ziqi Gao
- Department of Chemical Engineering, McMaster Universtiy, Hamilton, Ontario L8S 4L7, Canada
| | - Naveen Kumar Vasudevan
- Department of Chemical Engineering, McMaster Universtiy, Hamilton, Ontario L8S 4L7, Canada
| | - Hong Li
- School of Chemical Engineering and Technology, National Engineering Research Center of Distillation Technology, and Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin University, Tianjin 300072, China
| | - Xin Gao
- School of Chemical Engineering and Technology, National Engineering Research Center of Distillation Technology, and Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin University, Tianjin 300072, China
| | - Xingang Li
- School of Chemical Engineering and Technology, National Engineering Research Center of Distillation Technology, and Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin University, Tianjin 300072, China
| | - Li Xi
- Department of Chemical Engineering, McMaster Universtiy, Hamilton, Ontario L8S 4L7, Canada
| |
Collapse
|
14
|
Rahbari A, Hens R, Moultos OA, Dubbeldam D, Vlugt TJH. Multiple Free Energy Calculations from Single State Point Continuous Fractional Component Monte Carlo Simulation Using Umbrella Sampling. J Chem Theory Comput 2020; 16:1757-1767. [PMID: 31999461 PMCID: PMC7066647 DOI: 10.1021/acs.jctc.9b01097] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
We introduce an alternative method
to perform free energy calculations
for mixtures at multiple temperatures and pressures from a single
simulation, by combining umbrella sampling and the continuous fractional
component Monte Carlo method. One can perform a simulation of a mixture
at a certain pressure and temperature and accurately compute the chemical
potential at other pressures and temperatures close to the simulation
conditions. This method has the following advantages: (1) Accurate
estimates of the chemical potential as a function of pressure and
temperature are obtained from a single state simulation without additional
postprocessing. This can potentially reduce the number of simulations
of a system for free energy calculations for a specific temperature
and/or pressure range. (2) Partial molar volumes and enthalpies are
obtained directly from the estimated chemical potentials. We tested
our method for a Lennard-Jones system, aqueous mixtures of methanol
at T = 298 K and P = 1 bar, and
a mixture of ammonia, nitrogen, and hydrogen at T = 573 K and P = 800 bar. For pure methanol (N = 410 molecules), we observed that the estimated chemical
potentials from umbrella sampling are in excellent agreement with
the reference values obtained from independent simulations, for ΔT = ±15 K and ΔP = 100 bar (with
respect to the simulated system). For larger systems, this range becomes
smaller because the relative fluctuations of energy and volume become
smaller. Without sufficient overlap, the performance of the method
will become poor especially for nonlinear variations of the chemical
potential.
Collapse
Affiliation(s)
- Ahmadreza Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Remco Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - David Dubbeldam
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
15
|
Josephson TR, Singh R, Minkara MS, Fetisov EO, Siepmann JI. Partial molar properties from molecular simulation using multiple linear regression. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1648898] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Tyler R. Josephson
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, USA
| | - Ramanish Singh
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, USA
| | - Mona S. Minkara
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
| | - Evgenii O. Fetisov
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Chemical Physics and Analysis, Pacific Northwest National Laboratory, Richland, Washington, USA
| | - J. Ilja Siepmann
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, USA
| |
Collapse
|
16
|
Rahbari A, Hens R, Dubbeldam D, Vlugt TJH. Improving the accuracy of computing chemical potentials in CFCMC simulations. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1631497] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft, Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft, Netherlands
| | - D. Dubbeldam
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft, Netherlands
| |
Collapse
|
17
|
Effect of truncating electrostatic interactions on predicting thermodynamic properties of water–methanol systems. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1547824] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|