1
|
Cortés Morales AD, Diamantonis N, Economou IG, Peters CJ, Siepmann JI. Molecular Modeling of Double Retrograde Vaporization Using Monte Carlo Simulations and Equations of State. J Phys Chem B 2023; 127:3672-3681. [PMID: 37067787 DOI: 10.1021/acs.jpcb.3c00706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2023]
Abstract
Vapor-liquid equilibria of binary systems consisting of a low-boiling (i.e., more volatile) and a high-boiling compound may exhibit unexpected behavior near the critical point of the low-boiling compound. Near the critical temperature of the low-boiling compound and for compositions rich in the low-boiling compound, increasing the pressure may result in multiple crossings of the dew- and bubble-point curves. This phenomenon is often called double retrograde vaporization (or condensation) and may play a role in oil field operations and gas transport through pipelines, but the microscopic driving forces for the unusual shape of the dew-point curve are not well understood. Monte Carlo simulations in the constant-pressure, constant-temperature Gibbs ensemble using the united-atom version of the TraPPE force field were carried out for the methane/n-butane mixture at temperatures ranging from 0.95 to 1.05 of the reduced (T/Tc) temperature of methane. The simulations predict a wealth of additional thermodynamic data (densities and free energies of transfer) and structural data that are used to provide much needed molecular-level insights into the fluid properties associated with double retrograde vaporization. Simulated thermodynamic data are also compared with calculations using the Peng-Robinson and PC-SAFT equations of state.
Collapse
Affiliation(s)
- Angel D Cortés Morales
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Nikolaos Diamantonis
- Department of Chemical Engineering, The Petroleum Institute, Khalifa University of Science and Technology, P.O. Box 2533, Abu Dhabi, UAE
- Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research "Demokritos", GR-15310 Aghia Paraskevi Attikis, Greece
| | - Ioannis G Economou
- Department of Chemical Engineering, The Petroleum Institute, Khalifa University of Science and Technology, P.O. Box 2533, Abu Dhabi, UAE
- Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research "Demokritos", GR-15310 Aghia Paraskevi Attikis, Greece
| | - Cornelis J Peters
- Department of Chemical Engineering, The Petroleum Institute, Khalifa University of Science and Technology, P.O. Box 2533, Abu Dhabi, UAE
| | - J Ilja Siepmann
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
- Department Chemical Engineering and Materials Science, University of Minnesota, 421 Washington Avenue SE, Minneapolis, Minnesota 55455-0132, United States
| |
Collapse
|
2
|
Comprehensive review on physical properties of supercritical carbon dioxide calculated by molecular simulation. KOREAN J CHEM ENG 2023. [DOI: 10.1007/s11814-022-1316-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
3
|
Scheiner B, Yoon TJ. Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory. J Chem Phys 2021; 154:134101. [PMID: 33832259 DOI: 10.1063/5.0045211] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. The first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO2 is the C-C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.
Collapse
Affiliation(s)
- Brett Scheiner
- Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Tae Jun Yoon
- Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| |
Collapse
|
4
|
Mi W, Ramos P, Maranhao J, Pavanello M. Ab Initio Structure and Dynamics of CO 2 at Supercritical Conditions. J Phys Chem Lett 2019; 10:7554-7559. [PMID: 31738567 DOI: 10.1021/acs.jpclett.9b03054] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Green technologies rely on green solvents and fluids. Among them, supercritical CO2 already finds many important applications. The molecular-level understanding of the dynamics and structure of this supercritical fluid is a prerequisite for rational design of future green technologies. Unfortunately, the commonly employed Kohn-Sham density functional theory (DFT) is too computationally demanding to produce meaningfully converged dynamics within a reasonable time and with a reasonable computational effort. Thanks to subsystem DFT, we analyze finite-size effects by considering simulation cells of varying sizes (up to 256 independent molecules in the cell) and finite-time effects by running 100 ps trajectories. We find that the simulations are in reasonable and semiquantitative agreement with the available neutron diffraction experiments and that, as opposed to the gas phase, the CO2 molecules in the fluid are bent with an average OCO angle of 175.8°. Our simulations also confirm that the dimer T-shape is the most prevalent configuration. Our results further strengthen the experiment-simulation agreement for this fluid when comparing radial distribution functions and diffusion coefficient, confirming subsystem DFT as a viable tool for modeling structure and dynamics of condensed-phase systems.
Collapse
Affiliation(s)
- Wenhui Mi
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Pablo Ramos
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Jack Maranhao
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Michele Pavanello
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
- Department of Physics , Rutgers University , Newark , New Jersey 07102 , United States
| |
Collapse
|
5
|
Goel H, Windom ZW, Jackson AA, Rai N. Performance of density functionals for modeling vapor liquid equilibria of CO 2 and SO 2. J Comput Chem 2018; 39:397-406. [PMID: 29164642 DOI: 10.1002/jcc.25123] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2017] [Revised: 10/12/2017] [Accepted: 11/07/2017] [Indexed: 01/16/2023]
Abstract
Vapor liquid equilibria (VLE) and condensed phase properties of carbon dioxide and sulfur dioxide are calculated using first principles Monte Carlo (FPMC) simulations to assess the performance of several density functionals, notably PBE-D3, BLYP-D3, PBE0-D3, M062X-D3, and rVV10. GGA functionals were used to compute complete vapor liquid coexistence curves (VLCCs) to estimate critical properties, while the hybrid and nonlocal van der Waals functionals were used only for computing density at a single state point due to the high computational cost. Our results show that the BLYP-D3 functional performs well in predicting VLE properties for both molecules when compared with other functionals. In the liquid phase, pair correlation functions reveal that there is not a significant difference in the location of the peak for the first solvation shell while the peak heights are different for different functionals. Overall, the BLYP-D3 functional is a good choice for modeling VLE of acidic gases with significant environmental implications such as CO2 and SO2 . © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| | - Zachary W Windom
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| | - Amber A Jackson
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| |
Collapse
|
6
|
Fetisov EO, Shah MS, Long JR, Tsapatsis M, Siepmann JI. First principles Monte Carlo simulations of unary and binary adsorption: CO2, N2, and H2O in Mg-MOF-74. Chem Commun (Camb) 2018; 54:10816-10819. [DOI: 10.1039/c8cc06178e] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Dative bonding of adsorbate molecules onto coordinatively-unsaturated metal sites in metal–organic frameworks (MOFs) can lead to unique adsorption selectivities.
Collapse
Affiliation(s)
- Evgenii O. Fetisov
- Department of Chemistry and Chemical Theory Center
- University of Minnesota
- Minneapolis
- USA
| | - Mansi S. Shah
- Department of Chemistry and Chemical Theory Center
- University of Minnesota
- Minneapolis
- USA
- Department of Chemical Engineering and Materials Science
| | - Jeffrey R. Long
- Department of Chemistry and Chemical & Biomolecular Engineering
- University of California
- Berkeley
- USA
- Materials Sciences Division
| | - Michael Tsapatsis
- Department of Chemical Engineering and Materials Science
- University of Minnesota
- Minneapolis
- USA
| | - J. Ilja Siepmann
- Department of Chemistry and Chemical Theory Center
- University of Minnesota
- Minneapolis
- USA
- Department of Chemical Engineering and Materials Science
| |
Collapse
|
7
|
|
8
|
Mielke SL, Truhlar DG. Improved methods for Feynman path integral calculations and their application to calculate converged vibrational–rotational partition functions, free energies, enthalpies, entropies, and heat capacities for methane. J Chem Phys 2015; 142:044105. [PMID: 25637967 DOI: 10.1063/1.4905526] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Steven L. Mielke
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. S.E., Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
9
|
Cabral BJC, Rivelino R, Coutinho K, Canuto S. A first principles approach to the electronic properties of liquid and supercritical CO2. J Chem Phys 2015; 142:024504. [DOI: 10.1063/1.4905256] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Orozco GA, Economou IG, Panagiotopoulos AZ. Optimization of Intermolecular Potential Parameters for the CO2/H2O Mixture. J Phys Chem B 2014; 118:11504-11. [DOI: 10.1021/jp5067023] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Gustavo A. Orozco
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | | | |
Collapse
|
11
|
Gurina DL, Antipova ML, Petrenko VE. Structural features of binary mixtures of supercritical CO2 with polar entrainers by molecular dynamics simulation. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2013. [DOI: 10.1134/s0036024413100099] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
12
|
Pérez-Sánchez G, González-Salgado D, Piñeiro MM, Vega C. Fluid-solid equilibrium of carbon dioxide as obtained from computer simulations of several popular potential models: the role of the quadrupole. J Chem Phys 2013; 138:084506. [PMID: 23464159 DOI: 10.1063/1.4792443] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
In this work the solid-fluid equilibrium for carbon dioxide (CO2) has been evaluated using Monte Carlo simulations. In particular the melting curve of the solid phase denoted as I, or dry ice, was computed for pressures up to 1000 MPa. Four different models, widely used in computer simulations of CO2 were considered in the calculations. All of them are rigid non-polarizable models consisting of three Lennard-Jones interaction sites located on the positions of the atoms of the molecule, plus three partial charges. It will be shown that although these models predict similar vapor-liquid equilibria their predictions for the fluid-solid equilibria are quite different. Thus the prediction of the entire phase diagram is a severe test for any potential model. It has been found that the Transferable Potentials for Phase Equilibria (TraPPE) model yields the best description of the triple point properties and melting curve of carbon dioxide. It is shown that the ability of a certain model to predict the melting curve of carbon dioxide is related to the value of the quadrupole moment of the model. Models with low quadrupole moment tend to yield melting temperatures too low, whereas the model with the highest quadrupole moment yields the best predictions. That reinforces the idea that not only is the quadrupole needed to provide a reasonable description of the properties in the fluid phase, but also it is absolutely necessary to describe the properties of the solid phase.
Collapse
Affiliation(s)
- G Pérez-Sánchez
- Departamento de Física Aplicada, Universidade de Vigo, 36310 Vigo, Spain
| | | | | | | |
Collapse
|
13
|
Mielke SL, Dinpajooh M, Siepmann JI, Truhlar DG. Efficient methods for including quantum effects in Monte Carlo calculations of large systems: Extension of the displaced points path integral method and other effective potential methods to calculate properties and distributions. J Chem Phys 2013; 138:014110. [PMID: 23298031 DOI: 10.1063/1.4772667] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Steven L Mielke
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA.
| | | | | | | |
Collapse
|
14
|
El-Hendawy MM, Garate JA, English NJ, O'Reilly S, Mooney DA. Diffusion and interactions of carbon dioxide and oxygen in the vicinity of the active site of Rubisco: Molecular dynamics and quantum chemical studies. J Chem Phys 2012; 137:145103. [DOI: 10.1063/1.4757021] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Morad M El-Hendawy
- SFI Strategic Research Cluster in Solar Energy Conversion, University College Dublin, Belfield, Dublin 4, Ireland
| | | | | | | | | |
Collapse
|
15
|
Wang Z, Liu J, Zhang M, Cukier RI, Bu Y. Solvation and evolution dynamics of an excess electron in supercritical CO2. PHYSICAL REVIEW LETTERS 2012; 108:207601. [PMID: 23003186 DOI: 10.1103/physrevlett.108.207601] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2011] [Revised: 02/08/2012] [Indexed: 06/01/2023]
Abstract
We present an ab initio molecular dynamics simulation of the dynamics of an excess electron solvated in supercritical CO2. The excess electron can exist in three types of states: CO2-core localized, dual-core localized, and diffuse states. All these states undergo continuous state conversions via a combination of long lasting breathing oscillations and core switching, as also characterized by highly cooperative oscillations of the excess electron volume and vertical detachment energy. All of these oscillations exhibit a strong correlation with the electron-impacted bending vibration of the core CO2, and the core-switching is controlled by thermal fluctuations.
Collapse
Affiliation(s)
- Zhiping Wang
- The Center of Molecular Modeling & Simulation, Institute of Theoretical Chemistry, Shandong University, Jinan, 250100, People's Republic of China
| | | | | | | | | |
Collapse
|
16
|
Huang SN, Pascal TA, Goddard WA, Maiti PK, Lin ST. Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model. J Chem Theory Comput 2011; 7:1893-901. [DOI: 10.1021/ct200211b] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Shao-Nung Huang
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
| | - Tod A. Pascal
- Materials and Process Simulation Center, Beckman Institute, California Institute of Technology, Pasadena, California 91125, United States
- Graduate School of EEWS (WCU), Korea Advanced Institute of Science and Technology, Daejeon 305-701, Republic of Korea
| | - William A. Goddard
- Materials and Process Simulation Center, Beckman Institute, California Institute of Technology, Pasadena, California 91125, United States
- Graduate School of EEWS (WCU), Korea Advanced Institute of Science and Technology, Daejeon 305-701, Republic of Korea
| | - Prabal K. Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Shiang-Tai Lin
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|
17
|
Sarupria S, Debenedetti PG. Molecular Dynamics Study of Carbon Dioxide Hydrate Dissociation. J Phys Chem A 2011; 115:6102-11. [DOI: 10.1021/jp110868t] [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)
- Sapna Sarupria
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
18
|
Xu W, Yang J. A Computer Simulation Study on Self- and Cross-Aggregation of Multiple Polar Species in Supercritical Carbon Dioxide. J Phys Chem A 2010; 114:5414-28. [DOI: 10.1021/jp101092v] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Wenhao Xu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, PR China
| | - Jichu Yang
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, PR China
| |
Collapse
|
19
|
Herrebout WA, Nagels N, van der Veken BJ. On the ν1CO2/2âν2CO2Resonance in the Complex of Carbon Dioxide with Dimethyl Ether. Chemphyschem 2009; 10:3054-60. [DOI: 10.1002/cphc.200900414] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
20
|
Balasubramanian S, Kohlmeyer A, Klein ML. Ab initio molecular dynamics study of supercritical carbon dioxide including dispersion corrections. J Chem Phys 2009; 131:144506. [DOI: 10.1063/1.3245962] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|