1
|
Wasik D, Vicent-Luna JM, Rezaie S, Luna-Triguero A, Vlugt TJH, Calero S. The Impact of Metal Centers in the M-MOF-74 Series on Formic Acid Production. ACS APPLIED MATERIALS & INTERFACES 2024; 16:45006-45019. [PMID: 39141894 PMCID: PMC11367578 DOI: 10.1021/acsami.4c10678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/27/2024] [Revised: 08/01/2024] [Accepted: 08/08/2024] [Indexed: 08/16/2024]
Abstract
The confinement effect of porous materials on the thermodynamical equilibrium of the CO2 hydrogenation reaction presents a cost-effective alternative to transition metal catalysts. In metal-organic frameworks, the type of metal center has a greater impact on the enhancement of formic acid production than the scale of confinement resulting from the pore size. The M-MOF-74 series enables a comprehensive study of how different metal centers affect HCOOH production, minimizing the effect of pore size. In this work, molecular simulations were used to analyze the adsorption of HCOOH and the CO2 hydrogenation reaction in M-MOF-74, where M = Ni, Cu, Co, Fe, Mn, Zn. We combine classical simulations and density functional theory calculations to gain insights into the mechanisms that govern the low coverage adsorption of HCOOH in the surrounding of the metal centers of M-MOF-74. The impact of metal centers on the HCOOH yield was assessed by Monte Carlo simulations in the grand-canonical ensemble, using gas-phase compositions of CO2, H2, and HCOOH at chemical equilibrium at 298.15-800 K, 1-60 bar. The performance of M-MOF-74 in HCOOH production follows the same order as the uptake and the heat of HCOOH adsorption: Ni > Co > Fe > Mn > Zn > Cu. Ni-MOF-74 increases the mole fraction of HCOOH by ca. 105 times compared to the gas phase at 298.15 K, 60 bar. Ni-MOF-74 has the potential to be more economically attractive for CO2 conversion than transition metal catalysts, achieving HCOOH production at concentrations comparable to the highest formate levels reported for transition metal catalysts and offering a more valuable molecular form of the product.
Collapse
Affiliation(s)
- Dominika
O. Wasik
- Materials
Simulation and Modelling, Department of Applied Physics and Science
Education, Eindhoven University of Technology, 5600MB Eindhoven, The Netherlands
- Eindhoven
Institute for Renewable Energy Systems, Eindhoven University of Technology,
PO Box 513, 5600 MB Eindhoven, The Netherlands
| | - José Manuel Vicent-Luna
- Materials
Simulation and Modelling, Department of Applied Physics and Science
Education, Eindhoven University of Technology, 5600MB Eindhoven, The Netherlands
| | - Shima Rezaie
- Energy
Technology, Department of Mechanical Engineering, Eindhoven University of Technology, 5600MB Eindhoven, The Netherlands
| | - Azahara Luna-Triguero
- Eindhoven
Institute for Renewable Energy Systems, Eindhoven University of Technology,
PO Box 513, 5600 MB Eindhoven, The Netherlands
- Energy
Technology, Department of Mechanical Engineering, Eindhoven University of Technology, 5600MB Eindhoven, 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
| | - Sofía Calero
- Materials
Simulation and Modelling, Department of Applied Physics and Science
Education, Eindhoven University of Technology, 5600MB Eindhoven, The Netherlands
- Eindhoven
Institute for Renewable Energy Systems, Eindhoven University of Technology,
PO Box 513, 5600 MB Eindhoven, The Netherlands
| |
Collapse
|
2
|
Staňo R, van Lente J, Lindhoud S, Košovan P. Sequestration of Small Ions and Weak Acids and Bases by a Polyelectrolyte Complex Studied by Simulation and Experiment. Macromolecules 2024; 57:1383-1398. [PMID: 38370910 PMCID: PMC10867894 DOI: 10.1021/acs.macromol.3c01209] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 11/27/2023] [Accepted: 11/29/2023] [Indexed: 02/20/2024]
Abstract
Mixing of oppositely charged polyelectrolytes can result in phase separation into a polymer-poor supernatant and a polymer-rich polyelectrolyte complex (PEC). We present a new coarse-grained model for the Grand-reaction method that enables us to determine the composition of the coexisting phases in a broad range of pH and salt concentrations. We validate the model by comparing it to recent simulations and experimental studies, as well as our own experiments on poly(acrylic acid)/poly(allylamine hydrochloride) complexes. The simulations using our model predict that monovalent ions partition approximately equally between both phases, whereas divalent ones accumulate in the PEC phase. On a semiquantitative level, these results agree with our own experiments, as well as with other experiments and simulations in the literature. In the sequel, we use the model to study the partitioning of a weak diprotic acid at various pH values of the supernatant. Our results show that the ionization of the acid is enhanced in the PEC phase, resulting in its preferential accumulation in this phase, which monotonically increases with the pH. Currently, this effect is still waiting to be confirmed experimentally. We explore how the model parameters (particle size, charge density, permittivity, and solvent quality) affect the measured partition coefficients, showing that fine-tuning of these parameters can make the agreement with the experiments almost quantitative. Nevertheless, our results show that charge regulation in multivalent solutes can potentially be exploited in engineering the partitioning of charged molecules in PEC-based systems at various pH values.
Collapse
Affiliation(s)
- Roman Staňo
- Faculty
of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, 1090 Vienna, Austria
| | - Jéré
J. van Lente
- Department
of Molecules & Materials, University
of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands
| | - Saskia Lindhoud
- Department
of Molecules & Materials, University
of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands
| | - Peter Košovan
- Department
of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 40 Prague 2, Czech Republic
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Habibi P, Rahbari A, Blazquez S, Vega C, Dey P, Vlugt TJH, Moultos OA. A New Force Field for OH - for Computing Thermodynamic and Transport Properties of H 2 and O 2 in Aqueous NaOH and KOH Solutions. J Phys Chem B 2022; 126:9376-9387. [PMID: 36325986 PMCID: PMC9677430 DOI: 10.1021/acs.jpcb.2c06381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2022] [Revised: 10/20/2022] [Indexed: 11/05/2022]
Abstract
The thermophysical properties of aqueous electrolyte solutions are of interest for applications such as water electrolyzers and fuel cells. Molecular dynamics (MD) and continuous fractional component Monte Carlo (CFCMC) simulations are used to calculate densities, transport properties (i.e., self-diffusivities and dynamic viscosities), and solubilities of H2 and O2 in aqueous sodium and potassium hydroxide (NaOH and KOH) solutions for a wide electrolyte concentration range (0-8 mol/kg). Simulations are carried out for a temperature and pressure range of 298-353 K and 1-100 bar, respectively. The TIP4P/2005 water model is used in combination with a newly parametrized OH- force field for NaOH and KOH. The computed dynamic viscosities at 298 K for NaOH and KOH solutions are within 5% from the reported experimental data up to an electrolyte concentration of 6 mol/kg. For most of the thermodynamic conditions (especially at high concentrations, pressures, and temperatures) experimental data are largely lacking. We present an extensive collection of new data and engineering equations for H2 and O2 self-diffusivities and solubilities in NaOH and KOH solutions, which can be used for process design and optimization of efficient alkaline electrolyzers and fuel cells.
Collapse
Affiliation(s)
- Parsa Habibi
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628 CDDelft, The Netherlands
| | - Ahmadreza Rahbari
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| | - Samuel Blazquez
- Depto.
Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040Madrid, Spain
| | - Carlos Vega
- Depto.
Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040Madrid, Spain
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628 CDDelft, 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 CBDelft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| |
Collapse
|
6
|
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
|
7
|
Vo P, Forsman J, Woodward CE. Boundary-Monte Carlo Method for Neutral and Charged Confined Fluids. J Chem Theory Comput 2022; 18:3766-3780. [PMID: 35575645 DOI: 10.1021/acs.jctc.1c01146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we describe a new Monte Carlo (MC) simulation method to investigate highly coupled fluids in confined geometries at a constant chemical potential. This method is based on so-called multi-scale Hamiltonian methods, wherein the chemical potential is determined using a more amenable Hamiltonian for a fluid in an "outer" region, which facilitates standard methods, such as grand canonical MC simulations or Widom's particle insertion method. The (inner region) fluid of interest is placed in diffusive contact with the simpler outer fluid via a boundary zone wherein the Hamiltonian is transformed. The current method utilizes an ideal fluid for the outer regions, which allows for implicit rather than explicit simulations. Only the boundary and inner region need explicit consideration; hence, the nomenclature used is boundary-Monte Carlo. We illustrate the utility of the method for simple neutral and charged fluids in cylindrical and planar pores. In the latter case, we use a dense room-temperature ionic liquid model and illustrate how the boundary zone establishes a proper Donnan equilibrium between inner and outer fluids in the presence of charged planar electrodes. Thus, the method allows direct calculation of properties such as the differential capacitance, without the need for additional difficult calculations of the requisite Donnan potential.
Collapse
Affiliation(s)
- Phuong Vo
- School of Science, University of New South Wales, Canberra, Canberra ACT 2600, Australia
| | - Jan Forsman
- Department of Theoretical Chemistry, Chemical Centre, Lund University, Lund S-22100, Sweden
| | - Clifford E Woodward
- School of Science, University of New South Wales, Canberra, Canberra ACT 2600, Australia
| |
Collapse
|
8
|
Salehi HS, Polat HM, de Meyer F, Houriez C, Coquelet C, Vlugt TJH, Moultos OA. Vapor pressures and vapor phase compositions of choline chloride urea and choline chloride ethylene glycol deep eutectic solvents from molecular simulation. J Chem Phys 2021; 155:114504. [PMID: 34551525 DOI: 10.1063/5.0062408] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Despite the widespread acknowledgment that deep eutectic solvents (DESs) have negligible vapor pressures, very few studies in which the vapor pressures of these solvents are measured or computed are available. Similarly, the vapor phase composition is known for only a few DESs. In this study, for the first time, the vapor pressures and vapor phase compositions of choline chloride urea (ChClU) and choline chloride ethylene glycol (ChClEg) DESs are computed using Monte Carlo simulations. The partial pressures of the DES components were obtained from liquid and vapor phase excess Gibbs energies, computed using thermodynamic integration. The enthalpies of vaporization were computed from the obtained vapor pressures, and the results were in reasonable agreement with the few available experimental data in the literature. It was found that the vapor phases of both DESs were dominated by the most volatile component (hydrogen bond donor, HBD, i.e., urea or ethylene glycol), i.e., 100% HBD in ChClEg and 88%-93% HBD in ChClU. Higher vapor pressures were observed for ChClEg compared to ChClU due to the higher volatility of ethylene glycol compared to urea. The influence of the liquid composition of the DESs on the computed properties was studied by considering different mole fractions (i.e., 0.6, 0.67, and 0.75) of the HBD. Except for the partial pressure of ethylene glycol in ChClEg, all the computed partial pressures and enthalpies of vaporization showed insensitivity toward the liquid composition. The activity coefficient of ethylene glycol in ChClEg was computed at different liquid phase mole fractions, showing negative deviations from Raoult's law.
Collapse
Affiliation(s)
- 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
| | - H Mert Polat
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Frédérick de Meyer
- CCUS and Acid Gas Entity, Liquefied Natural Gas Department, Exploration Production, Total Energies S.E., 92078 Paris, 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
| | - 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
|
9
|
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
|
10
|
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
|
11
|
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
|
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
|
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
|
14
|
Vo P, Lu H, Ma K, Forsman J, Woodward CE. Local Grand Canonical Monte Carlo Simulation Method for Confined Fluids. J Chem Theory Comput 2019; 15:6944-6957. [PMID: 31665596 DOI: 10.1021/acs.jctc.9b00804] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
We describe a new local grand canonical Monte Carlo method to treat fluids in pores in chemical equilibrium with a reference bulk. The method is applied to Lennard-Jones particles in pores of different geometry and is shown to be much more accurate and efficient than other techniques such as traditional grand canonical simulations or Widom's particle insertion method. It utilizes a penalty potential to create a gas phase, which is in equilibrium with a more dense liquid component in the pore. Grand canonical Monte Carlo moves are employed in the gas phase, and the system then maintains chemical equilibrium by "diffusion" of particles. This creates an interface, which means that the confined fluid needs to occupy a large enough volume so that this is not an issue. We also applied the method to confined charged fluids and show how it can be used to determine local electrostatic potentials in the confined fluid, which are properly referenced to the bulk. This precludes the need to determine the Donnan potential (which controls electrochemical equilibrium) explicitly. Prior approaches have used explicit bulk simulations to measure this potential difference, which are significantly costly from a computational point of view. One outcome of our analysis is that pores of finite cross-section create a potential difference with the bulk via a small but nonzero linear charge density, which diminishes as ∼1/ln(L), where L is the pore length.
Collapse
Affiliation(s)
- Phuong Vo
- School of Science , University of New South Wales, Canberra , Canberra ACT 2600 , Australia
| | - Hongduo Lu
- Department of Theoretical Chemistry, Chemical Centre , Lund University P. O. Box 124, S-22100 Lund , Sweden
| | - Ke Ma
- School of Materials Science and Engineering , Tianjin University of Technology Tianjin 300384 , People's Republic of China
| | - Jan Forsman
- Department of Theoretical Chemistry, Chemical Centre , Lund University P. O. Box 124, S-22100 Lund , Sweden
| | - Clifford E Woodward
- School of Science , University of New South Wales, Canberra , Canberra ACT 2600 , Australia
| |
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
|
Noroozi J, Smith WR. An Efficient Molecular Simulation Methodology for Chemical Reaction Equilibria in Electrolyte Solutions: Application to CO2 Reactive Absorption. J Phys Chem A 2019; 123:4074-4086. [DOI: 10.1021/acs.jpca.9b00302] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Javad Noroozi
- Department of Chemical Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - William R. Smith
- Department of Chemical Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
- Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario N1G 2W1, Canada
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, Ontario L1H 7K4, Canada
| |
Collapse
|
18
|
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]
|
19
|
Affiliation(s)
- Braden Kelly
- Department of Mathematics and Statistics, University of Guelph, Guelph, Canada
| | - William R. Smith
- Department of Mathematics and Statistics, University of Guelph, Guelph, Canada
- Department of Chemistry, University of Guelph, Guelph, Canada
- Department of Chemical Engineering, University of Waterloo, Waterloo, Canada
- Faculty of Science, University of Ontario Institute of Technology, Oshawa Canada
| |
Collapse
|
20
|
Nikolaidis IK, Poursaeidesfahani A, Csaszar Z, Ramdin M, Vlugt TJH, Economou IG, Moultos OA. Modeling the phase equilibria of asymmetric hydrocarbon mixtures using molecular simulation and equations of state. AIChE J 2018. [DOI: 10.1002/aic.16453] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Ilias K. Nikolaidis
- National Center for Scientific Research “Demokritos”, Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory; Aghia Paraskevi Attikis Greece
- School of Chemical Engineering; National Technical University of Athens; Athens Greece
| | - Ali Poursaeidesfahani
- Engineering Thermodynamics, Process and Energy Dept., Faculty of Mechanical, Maritime and Materials Engineering; Delft University of Technology; Delft The Netherlands
| | - Zsolt Csaszar
- Engineering Thermodynamics, Process and Energy Dept., Faculty of Mechanical, Maritime and Materials Engineering; Delft University of Technology; Delft The Netherlands
| | - Mahinder Ramdin
- Engineering Thermodynamics, Process and Energy Dept., Faculty of Mechanical, Maritime and Materials Engineering; Delft University of Technology; Delft The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process and Energy Dept., Faculty of Mechanical, Maritime and Materials Engineering; Delft University of Technology; Delft The Netherlands
| | - Ioannis G. Economou
- National Center for Scientific Research “Demokritos”, Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory; Aghia Paraskevi Attikis Greece
- Chemical Engineering Program; Texas A&M University at Qatar; Doha Qatar
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process and Energy Dept., Faculty of Mechanical, Maritime and Materials Engineering; Delft University of Technology; Delft The Netherlands
| |
Collapse
|
21
|
Smith WR, Qi W. Molecular Simulation of Chemical Reaction Equilibrium by Computationally Efficient Free Energy Minimization. ACS CENTRAL SCIENCE 2018; 4:1185-1193. [PMID: 30276252 PMCID: PMC6161046 DOI: 10.1021/acscentsci.8b00361] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2018] [Indexed: 05/25/2023]
Abstract
The molecular simulation of chemical reaction equilibrium (CRE) is a challenging and important problem of broad applicability in chemistry and chemical engineering. The primary molecular-based approach for solving this problem has been the reaction ensemble Monte Carlo (REMC) algorithm [Turner et al. Molec. Simulation2008, 34, (2), 119-146], based on classical force-field methodology. In spite of the vast improvements in computer hardware and software since its original development almost 25 years ago, its more widespread application is impeded by its computational inefficiency. A fundamental problem is that its MC basis inhibits the implementation of significant parallelization, and its successful implementation often requires system-specific tailoring and the incorporation of special MC approaches such as replica exchange, expanded ensemble, umbrella sampling, configurational bias, and continuous fractional component methodologies. We describe herein a novel CRE algorithm (reaction ensemble molecular dynamics, ReMD) that exploits modern computer hardware and software capabilities, and which can be straightforwardly implemented for systems of arbitrary size and complexity by exploiting the parallel computing methodology incorporated within many MD software packages (herein, we use GROMACS for illustrative purposes). The ReMD algorithm utilizes these features in the context of a macroscopically inspired and generally applicable free energy minimization approach based on the iterative approximation of the system Gibbs free energy function by a mathematically simple convex ideal solution model using the composition at each iteration as a reference state. Finally, we additionally describe a simple and computationally efficient a posteriori method to estimate the equilibrium concentrations of species present in very small amounts relative to others in the primary calculation. To demonstrate the algorithm, we show its application to two classic example systems considered previously in the literature: the N2-O2-NO system and the ammonia synthesis system.
Collapse
Affiliation(s)
- William R. Smith
- Department
of Mathematics and Statistics, University
of Guelph, Guelph, Ontario N1G 2W1, Canada
- Department
of Chemistry, University of Guelph, Guelph, Ontario N1G 2W1, Canada
- Department
of Chemical Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
- Faculty
of Science, University of Ontario Institute
of Technology, Oshawa, Ontario L1H 7K4, Canada
| | - Weikai Qi
- Department
of Mathematics and Statistics, University
of Guelph, Guelph, Ontario N1G 2W1, Canada
| |
Collapse
|
22
|
Mullen RG, Corcelli SA, Maginn EJ. Reaction Ensemble Monte Carlo Simulations of CO 2 Absorption in the Reactive Ionic Liquid Triethyl(octyl)phosphonium 2-Cyanopyrrolide. J Phys Chem Lett 2018; 9:5213-5218. [PMID: 30136851 DOI: 10.1021/acs.jpclett.8b02304] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The absorption of CO2 into an aprotic heterocyclic anion ionic liquid (IL) is modeled using reaction ensemble Monte Carlo (RxMC) with the semigrand reaction move. RxMC has previously been unable to sample chemical equilibrium involving molecular ions in nanostructured liquids due to the high free-energy requirements to open and close cavities and restructure the surrounding environment. Our results are validated by experiments in the modeled IL, triethyl(octyl)phosphonium 2-cyanopyrrolide ([P2228][cnp]), and in a close analog with longer alkyl chains on the cation. Heats of absorption and reaction from both experiment and simulation are exothermic and of comparable magnitude. Replacing experimental Henry's constants with their simulated counterparts improves the accuracy of a Langmuir-type model at moderate pressures. Nonidealities that affect chemical equilibrium are identified and calculated with high precision.
Collapse
Affiliation(s)
- Ryan Gotchy Mullen
- Department of Chemical and Biomolecular Engineering , University of Notre Dame , Notre Dame , Indiana 46556 , United States
- Physical and Life Sciences Directorate , Lawrence Livermore National Laboratory , Livermore , California 94550 , United States
| | - Steven A Corcelli
- Department of Chemistry and Biochemistry , University of Notre Dame , Notre Dame , Indiana 46556 , United States
| | - Edward J Maginn
- Department of Chemical and Biomolecular Engineering , University of Notre Dame , Notre Dame , Indiana 46556 , United States
| |
Collapse
|
23
|
Rahbari A, Ramdin M, van den Broeke LJP, Vlugt TJH. Combined Steam Reforming of Methane and Formic Acid To Produce Syngas with an Adjustable H 2:CO Ratio. Ind Eng Chem Res 2018; 57:10663-10674. [PMID: 30270977 PMCID: PMC6156100 DOI: 10.1021/acs.iecr.8b02443] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2018] [Revised: 07/16/2018] [Accepted: 07/17/2018] [Indexed: 12/19/2022]
Abstract
![]()
Syngas
is an important intermediate in the chemical process industry. It
is used for the production of hydrocarbons, acetic acid, oxo-alcohols,
and other chemicals. Depending on the target product and stoichiometry
of the reaction, an optimum (molar) ratio between hydrogen and carbon
monoxide (H2:CO) in the syngas is required. Different technologies
are available to control the H2:CO molar ratio in the syngas.
The combination of steam reforming of methane (SRM) and the water-gas
shift (WGS) reaction is the most established approach for syngas production.
In this work, to adjust the H2:CO ratio, we have considered
formic acid (FA) as a source for both hydrogen and carbon monoxide.
Using thermochemical equilibrium calculations, we show that the syngas
composition can be controlled by cofeeding formic acid into the SRM
process. The H2:CO molar ratio can be adjusted to a value
between one and three by adjusting the concentration of FA in the
reaction feed. At steam reforming conditions, typically above 900
K, FA can decompose to water and carbon monoxide and/or to hydrogen
and carbon dioxide. Our results show that cofeeding FA into the SRM
process can adjust the H2:CO molar ratio in a single step.
This can potentially be an alternative to the WGS process.
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
| | - 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
|
24
|
Rahbari A, Hens R, Nikolaidis IK, Poursaeidesfahani A, Ramdin M, Economou IG, Moultos OA, Dubbeldam D, Vlugt TJH. Computation of partial molar properties using continuous fractional component Monte Carlo. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1451663] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - I. K. Nikolaidis
- Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research Demokritos, Institute of Nanoscience and Nanotechnology , Attikis, Greece
| | - A. Poursaeidesfahani
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - I. G. Economou
- Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research Demokritos, Institute of Nanoscience and Nanotechnology , Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar , Doha, Qatar
| | - O. A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - D. Dubbeldam
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam , Amsterdam, The Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| |
Collapse
|
25
|
Matito-Martos I, Rahbari A, Martin-Calvo A, Dubbeldam D, Vlugt TJH, Calero S. Adsorption equilibrium of nitrogen dioxide in porous materials. Phys Chem Chem Phys 2018; 20:4189-4199. [PMID: 29362749 DOI: 10.1039/c7cp08017d] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effect of confinement on the equilibrium reactive system containing nitrogen dioxide and dinitrogen tetroxide is studied by molecular simulation and the reactive Monte Carlo (RxMC) approach. The bulk-phase reaction was successfully reproduced and five all-silica zeolites (i.e. FAU, FER, MFI, MOR, and TON) with different topologies were selected to study their adoption behavior. Dinitrogen tetroxide showed a stronger affinity than nitrogen dioxide in all the zeolites due to size effects, but exclusive adsorption sites in MOR allowed the adsorption of nitrogen dioxide with no competition at these sites. From the study of the adsorption isotherms and isobars of the reacting mixture, confinement enhanced the formation of dimers over the full range of pressure and temperature, finding the largest deviations from bulk fractions at low temperature and high pressure. The channel size and shape of the zeolite have a noticeable influence on the dinitrogen tetroxide formation, being more important in MFI, closely followed by TON and MOR, and finally FER and FAU. Preferential adsorption sites in MOR lead to an unusually strong selective adsorption towards nitrogen dioxide, demonstrating that the topological structure has a crucial influence on the composition of the mixture and must be carefully considered in systems containing nitrogen dioxide.
Collapse
Affiliation(s)
- I Matito-Martos
- Department of Physical, Chemical and Natural Systems, University Pablo de Olavide, Sevilla 41013, Spain.
| | | | | | | | | | | |
Collapse
|
26
|
Mullen RG, Maginn EJ. Reaction Ensemble Monte Carlo Simulation of Xylene Isomerization in Bulk Phases and under Confinement. J Chem Theory Comput 2017; 13:4054-4062. [DOI: 10.1021/acs.jctc.7b00498] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ryan Gotchy Mullen
- Department of Chemical and
Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana, United States
| | - Edward J. Maginn
- Department of Chemical and
Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana, United States
| |
Collapse
|