51
|
Lara-Peña M, Domínguez H. Gas sorption in solid surfaces: a computational study using rigid and Einstein-solid models. Phys Chem Chem Phys 2015; 17:27894-9. [PMID: 26434638 DOI: 10.1039/c5cp04676a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The reactive Monte Carlo (RxMC) method was proposed to describe the sorption of gases in solid materials due to the chemical reaction A + B ⇌ C. Two models were used to simulate the solid; the first model considered simulations with rigid particles in the solid whereas in the second model the particles were allowed to vibrate inside the solid with a given spring constant, i.e. an Einstein solid was used to simulate the substrate. In both models not only physisorption but also chemisorption of the fluid was observed. Sorption curves, at different spring constants, were simulated and it was noted that sorption was always enhanced with the Einstein solid model. Moreover, an inverse dependent function of the spring constant with the temperature was found. Finally, the second model might be used to explain the unusual sorption behavior observed in actual experimental reactions such as CO2 + Li2O ⇌ Li2CO3.
Collapse
Affiliation(s)
- Mayra Lara-Peña
- Posgrado en Ciencias Físicas, Universidad Nacional Autónoma de México, México, D.F. 04510, Mexico
| | | |
Collapse
|
52
|
Balaji SP, Gangarapu S, Ramdin M, Torres-Knoop A, Zuilhof H, Goetheer EL, Dubbeldam D, Vlugt TJ. Simulating the Reactions of CO2 in Aqueous Monoethanolamine Solution by Reaction Ensemble Monte Carlo Using the Continuous Fractional Component Method. J Chem Theory Comput 2015; 11:2661-9. [DOI: 10.1021/acs.jctc.5b00160] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sayee Prasaad Balaji
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime, and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Satesh Gangarapu
- Laboratory
of Organic Chemistry, Wageningen University, Dreijenplein 8, 6703HB Wageningen, 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
| | - Ariana Torres-Knoop
- Van’t
Hoff Institute for Molecular Sciences, University of Amsterdam, Science
Park 904, 1098XH Amsterdam, The Netherlands
| | - Han Zuilhof
- Laboratory
of Organic Chemistry, Wageningen University, Dreijenplein 8, 6703HB Wageningen, The Netherlands
- Department
of Chemical and Materials Engineering, King Abdulaziz University, Jeddah, Saudi Arabia
| | | | - 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
|
53
|
Schnell SK, Skorpa R, Bedeaux D, Kjelstrup S, Vlugt TJH, Simon JM. Partial molar enthalpies and reaction enthalpies from equilibrium molecular dynamics simulation. J Chem Phys 2015; 141:144501. [PMID: 25318729 DOI: 10.1063/1.4896939] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
We present a new molecular simulation technique for determining partial molar enthalpies in mixtures of gases and liquids from single simulations, without relying on particle insertions, deletions, or identity changes. The method can also be applied to systems with chemical reactions. We demonstrate our method for binary mixtures of Weeks-Chandler-Anderson particles by comparing with conventional simulation techniques, as well as for a simple model that mimics a chemical reaction. The method considers small subsystems inside a large reservoir (i.e., the simulation box), and uses the construction of Hill to compute properties in the thermodynamic limit from small-scale fluctuations. Results obtained with the new method are in excellent agreement with those from previous methods. Especially for modeling chemical reactions, our method can be a valuable tool for determining reaction enthalpies directly from a single MD simulation.
Collapse
Affiliation(s)
- Sondre K Schnell
- Process and Energy Laboratory, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Ragnhild Skorpa
- Department of Chemistry, Faculty of Natural Science and Technology, Norwegian University of Science and Technology, 4791 Trondheim, Norway
| | - Dick Bedeaux
- Department of Chemistry, Faculty of Natural Science and Technology, Norwegian University of Science and Technology, 4791 Trondheim, Norway
| | - Signe Kjelstrup
- Department of Chemistry, Faculty of Natural Science and Technology, Norwegian University of Science and Technology, 4791 Trondheim, Norway
| | - Thijs J H Vlugt
- Process and Energy Laboratory, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Jean-Marc Simon
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303, CNRS-Université de Bourgogne, 9, av. Savary, 21000 Dijon, France
| |
Collapse
|
54
|
Mogo C, Brandão J. The READY program: Building a global potential energy surface and reactive dynamic simulations for the hydrogen combustion. J Comput Chem 2014; 35:1330-7. [PMID: 24760762 DOI: 10.1002/jcc.23621] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2013] [Revised: 04/10/2014] [Accepted: 04/10/2014] [Indexed: 11/09/2022]
Abstract
READY (REActive DYnamics) is a program for studying reactive dynamic systems using a global potential energy surface (PES) built from previously existing PESs corresponding to each of the most important elementary reactions present in the system. We present an application to the combustion dynamics of a mixture of hydrogen and oxygen using accurate PESs for all the systems involving up to four oxygen and hydrogen atoms. Results at the temperature of 4000 K and pressure of 2 atm are presented and compared with model based on rate constants. Drawbacks and advantages of this approach are discussed and future directions of research are pointed out.
Collapse
Affiliation(s)
- César Mogo
- Centro de Investigação em Química do Algarve Universidade do Algarve, Portugal
| | | |
Collapse
|
55
|
Dubbeldam D, Torres-Knoop A, Walton KS. On the inner workings of Monte Carlo codes. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.819102] [Citation(s) in RCA: 202] [Impact Index Per Article: 18.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
56
|
Lísal M, Předota M, Brennan JK. Molecular-level simulations of chemical reaction equilibrium and diffusion in slit and cylindrical nanopores: model dimerisation reactions. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.797576] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
57
|
Bourasseau E, Maillet JB, Desbiens N, Stoltz G. Microscopic Calculations of Hugoniot Curves of Neat Triaminotrinitrobenzene (TATB) and of Its Detonation Products. J Phys Chem A 2011; 115:10729-37. [DOI: 10.1021/jp2047739] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Emeric Bourasseau
- Commissariat à l’Energie Atomique, DAM, DIF, F-91297 Arpajon, France
| | | | - Nicolas Desbiens
- Commissariat à l’Energie Atomique, DAM, DIF, F-91297 Arpajon, France
| | - Gabriel Stoltz
- Université Paris Est, CERMICS, MICMAC Project-team, INRIA-Ecole des Ponts ParisTech, 6 et 8 Av. Blaise Pascal, 77455 Marne-la-Vallée Cedex 2, France
| |
Collapse
|
58
|
Moučka F, Lísal M, Škvor J, Jirsák J, Nezbeda I, Smith WR. Molecular Simulation of Aqueous Electrolyte Solubility. 2. Osmotic Ensemble Monte Carlo Methodology for Free Energy and Solubility Calculations and Application to NaCl. J Phys Chem B 2011; 115:7849-61. [DOI: 10.1021/jp202054d] [Citation(s) in RCA: 90] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Filip Moučka
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, ON L1H7K4, Canada
| | - Martin Lísal
- E. Hála Laboratory of Thermodynamics, Institute of Chemical Process Fundamentals of the ASCR, v. v. i., 165 02 Prague 6, Czech Republic
| | - Jiří Škvor
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, ON L1H7K4, Canada
| | - Jan Jirsák
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, ON L1H7K4, Canada
- E. Hála Laboratory of Thermodynamics, Institute of Chemical Process Fundamentals of the ASCR, v. v. i., 165 02 Prague 6, Czech Republic
| | - Ivo Nezbeda
- E. Hála Laboratory of Thermodynamics, Institute of Chemical Process Fundamentals of the ASCR, v. v. i., 165 02 Prague 6, Czech Republic
| | - William R. Smith
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, ON L1H7K4, Canada
| |
Collapse
|
59
|
Rosch TW, Maginn EJ. Reaction Ensemble Monte Carlo Simulation of Complex Molecular Systems. J Chem Theory Comput 2011; 7:269-79. [PMID: 26596150 DOI: 10.1021/ct100615j] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Acceptance rules for reaction ensemble Monte Carlo (RxMC) simulations containing classically modeled atomistic degrees of freedom are derived for complex molecular systems where insertions and deletions are achieved gradually by utilizing the continuous fractional component (CFC) method. A self-consistent manner in which to utilize statistical mechanical data contained in ideal gas free energy parameters during RxMC moves is presented. The method is tested by applying it to two previously studied systems containing intramolecular degrees of freedom: the propene metathesis reaction and methyl-tert-butyl-ether (MTBE) synthesis. Quantitative agreement is found between the current results and those of Keil et al. (J. Chem. Phys. 2005, 122, 164705) for the propene metathesis reaction. Differences are observed between the equilibrium concentrations of the present study and those of Lísal et al. (AIChE J. 2000, 46, 866-875) for the MTBE reaction. It is shown that most of this difference can be attributed to an incorrect formulation of the Monte Carlo acceptance rule. Efficiency gains using CFC MC as opposed to single stage molecule insertions are presented.
Collapse
Affiliation(s)
- Thomas W Rosch
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, 182 Fitzpatrick Hall, Notre Dame, Indiana 46556-5637, United States
| | - Edward J Maginn
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, 182 Fitzpatrick Hall, Notre Dame, Indiana 46556-5637, United States
| |
Collapse
|
60
|
Bourasseau E, Maillet JB. Coupling microscopic and mesoscopic scales to simulate chemical equilibrium between a nanometric carbon cluster and detonation products fluid. Phys Chem Chem Phys 2011; 13:7060-70. [DOI: 10.1039/c0cp02622k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
61
|
Gubbins KE, Liu YC, Moore JD, Palmer JC. The role of molecular modeling in confined systems: impact and prospects. Phys Chem Chem Phys 2011; 13:58-85. [DOI: 10.1039/c0cp01475c] [Citation(s) in RCA: 136] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
62
|
Long Y, Palmer JC, Coasne B, Śliwinska-Bartkowiak M, Gubbins KE. Pressure enhancement in carbon nanopores: a major confinement effect. Phys Chem Chem Phys 2011; 13:17163-70. [DOI: 10.1039/c1cp21407a] [Citation(s) in RCA: 113] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Phenomena that only occur at high pressures in bulk phases are observed in nanopores because of a pressure enhancement effect.
Collapse
Affiliation(s)
- Yun Long
- Department of Chemical and Biomolecular Engineering
- North Carolina State University
- Raleigh
- USA
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering
- North Carolina State University
- Raleigh
- USA
| | - Benoit Coasne
- Institut Charles Gerhardt Monpellier
- CNRS (UMR 5253)
- Université de Montpellier 2 and CNRS (UMR 5253)
- Montpellier
- France
| | | | - Keith E. Gubbins
- Department of Chemical and Biomolecular Engineering
- North Carolina State University
- Raleigh
- USA
| |
Collapse
|
63
|
Farah K, Karimi-Varzaneh HA, Müller-Plathe F, Böhm MC. Reactive Molecular Dynamics with Material-Specific Coarse-Grained Potentials: Growth of Polystyrene Chains from Styrene Monomers. J Phys Chem B 2010; 114:13656-66. [DOI: 10.1021/jp107001e] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Karim Farah
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Petersenstrasse 20, D-64287 Darmstadt, Germany
| | - Hossein A. Karimi-Varzaneh
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Petersenstrasse 20, D-64287 Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Petersenstrasse 20, D-64287 Darmstadt, Germany
| | - Michael C. Böhm
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Petersenstrasse 20, D-64287 Darmstadt, Germany
| |
Collapse
|
64
|
Hansen N, Krishna R, van Baten J, Bell A, Keil F. Reactor simulation of benzene ethylation and ethane dehydrogenation catalyzed by ZSM-5: A multiscale approach. Chem Eng Sci 2010. [DOI: 10.1016/j.ces.2009.12.028] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
65
|
Affiliation(s)
- Doros N. Theodorou
- Department of Materials Science and Engineering, School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Zografou Campus, Athens 15780, Greece
| |
Collapse
|
66
|
Bourasseau E, Lachet V, Desbiens N, Maillet JB, Teuler JM, Ungerer P. Thermodynamic behavior of the CO2 + NO2/N2O4 mixture: a Monte Carlo simulation study. J Phys Chem B 2009; 112:15783-92. [PMID: 19367990 DOI: 10.1021/jp8068255] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The thermodynamic behavior of the carbon dioxide + nitrogen dioxide (CO2 + NO2) mixture was investigated using a Monte Carlo molecular simulation approach. This system is a particularly challenging one because nitrogen dioxide exists as a mixture of monomers (NO2) and dimers (N2O4) under certain pressure and temperature conditions. The chemical equilibrium between N2O4 and 2NO2 and the vapor-liquid equilibrium of CO2 + NO2/N2O4 mixtures were simulated using simultaneously the reaction ensemble and the Gibbs ensemble Monte Carlo (RxMC and GEMC) methods. Rigid all atoms molecular potentials bearing point charges were proposed to model both NO2 and N2O4 species. Liquid-vapor coexistence properties of the reacting NO2/N2O4 system were first investigated. The calculated vapor pressures and coexisting densities were compared to experimental values, leading to an average deviation of 10% for vapor pressures and 6% for liquid densities. The critical region was also addressed successfully using the subcritical Monte Carlo simulation results and some appropriate scaling laws. Predictions of CO2 + NO2/N2O4 phase diagrams at 300, 313, and 330 K were then proposed. Derivative properties calculations were also performed in the reaction ensemble at constant pressure and temperature for both NO2/N2O4 and CO2 + NO2/N2O4 systems. The calculated heat capacities show a maximum in the temperature range where N2O4 dissociation occurs, in agreement with available experimental data.
Collapse
Affiliation(s)
- Emeric Bourasseau
- Commissariat à l'Energie Atomique, Centre DAM-Ile de France, Département de Physique Théorique et Appliquée, Bruyères-le-Châtel, 91297 Arpajon Cedex, France.
| | | | | | | | | | | |
Collapse
|
67
|
Coe JD, Sewell TD, Shaw MS. Optimal sampling efficiency in Monte Carlo simulation with an approximate potential. J Chem Phys 2009; 130:164104. [PMID: 19405558 DOI: 10.1063/1.3116788] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Building on the work of Iftimie et al. [J. Chem. Phys. 113, 4852 (2000)] and Gelb [J. Chem. Phys. 118, 7747 (2003)], Boltzmann sampling of an approximate potential (the "reference" system) is used to build a Markov chain in the isothermal-isobaric ensemble. At the end points of the chain, the energy is evaluated at a more accurate level (the "full" system) and a composite move encompassing all of the intervening steps is accepted on the basis of a modified Metropolis criterion. For reference system chains of sufficient length, consecutive full energies are statistically decorrelated and thus far fewer are required to build ensemble averages with a given variance. Without modifying the original algorithm, however, the maximum reference chain length is too short to decorrelate full configurations without dramatically lowering the acceptance probability of the composite move. This difficulty stems from the fact that the reference and full potentials sample different statistical distributions. By manipulating the thermodynamic variables characterizing the reference system (pressure and temperature, in this case), we maximize the average acceptance probability of composite moves, lengthening significantly the random walk between consecutive full energy evaluations. In this manner, the number of full energy evaluations needed to precisely characterize equilibrium properties is dramatically reduced. The method is applied to a model fluid, but implications for sampling high-dimensional systems with ab initio or density functional theory potentials are discussed.
Collapse
Affiliation(s)
- Joshua D Coe
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | | | |
Collapse
|
68
|
Malijevský A, Lísal M. Density functional study of chemical reaction equilibrium for dimerization reactions in slit and cylindrical nanopores. J Chem Phys 2009; 130:164713. [DOI: 10.1063/1.3125925] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|