1
|
Morishita T, Shiga M. Ab Initio Characterization of the CO 2-Water Interface Using Unsupervised Machine Learning for Dimensionality Reduction. J Phys Chem B 2024; 128:5781-5791. [PMID: 38829554 DOI: 10.1021/acs.jpcb.4c01526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
Precise characterization of the supercritical CO2-water interface under high pressure and temperature conditions is crucial for the geological storage of carbon dioxide (CO2) in deep saline aquifers. Molecular dynamics (MD) simulations offer a valuable approach to gaining insight into the CO2-water interface at a microscopic level. However, no attempt has been made to characterize the CO2-water interface with the accuracy afforded by ab initio calculations. In this study, we performed ab initio MD (AIMD) simulations to investigate the structural and dynamical properties of the CO2-water interface, comparing the results with those obtained from classical force-field MD (FF-MD) simulations. Molecular orientation at the interface was well reproduced in both AIMD and FF-MD simulations. Characteristic structural fluctuations of water at the interface were unveiled by applying multidimensional scaling and time-dependent principal component analysis to the AIMD trajectories; however, they were not prominent in the FF-MD simulations. Furthermore, our study demonstrated a marked difference in the residence time of molecules in the interface region between AIMD and FF-MD simulations, indicating that time-dependent properties of the CO2-water interface strongly depend on the description of the intermolecular forces.
Collapse
Affiliation(s)
- Tetsuya Morishita
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Masashige Shiga
- Geological Survey of Japan, National Institute of Advanced Industrial Science and Technology (AIST), Central 7, 1-1-1 Higashi, Tsukuba, Ibaraki 305-8567, Japan
| |
Collapse
|
2
|
Hulikal Chakrapani T, Hajibeygi H, Moultos OA, Vlugt TJH. Mutual Diffusivities of Mixtures of Carbon Dioxide and Hydrogen and Their Solubilities in Brine: Insight from Molecular Simulations. Ind Eng Chem Res 2024; 63:10456-10481. [PMID: 38882502 PMCID: PMC11177264 DOI: 10.1021/acs.iecr.4c01078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 05/14/2024] [Accepted: 05/20/2024] [Indexed: 06/18/2024]
Abstract
H2-CO2 mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H2-CO2 mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H2 due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO2, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO2 content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO2-H2-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO2 and H2 in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO2 in NaCl brine decreased in the ternary system compared to the binary CO2-NaCl brine system, the solubility of H2 in NaCl brine increased less in the ternary system compared to the binary H2-NaCl brine system. The cooperative effect of H2-CO2 enhances the H2 solubility while suppressing the CO2 solubility. The water content in the gas phase was found to be intermediate between H2-NaCl brine and CO2-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO2-H2 mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H2-CO2 mixtures.
Collapse
Affiliation(s)
- Thejas Hulikal Chakrapani
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Hadi Hajibeygi
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| |
Collapse
|
3
|
Shiga M, Morishita T, Nishiyama N, Sorai M, Aichi M, Abe A. Atomic-Scale Insights into the Phase Behavior of Carbon Dioxide and Water from 313 to 573 K and 8 to 30 MPa. ACS OMEGA 2024; 9:20976-20987. [PMID: 38764624 PMCID: PMC11097351 DOI: 10.1021/acsomega.4c00133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 03/17/2024] [Accepted: 04/10/2024] [Indexed: 05/21/2024]
Abstract
We performed molecular dynamics (MD) simulations of CO2 + H2O systems by employing widely used force fields (EPM2, TraPPE, and PPL models for CO2; SPC/E and TIP4P/2005 models for H2O). The phase behavior observed in our MD simulations is consistent with the coexistence lines obtained from previous experiments and SAFT-based theoretical models for the equations of state. Our structural analysis reveals a pronounced correlation between phase transitions and the structural orderliness. Specifically, the coordination number of Ow (oxygen in H2O) around other Ow significantly correlates with phase changes. In contrast, coordination numbers pertaining to the CO2 molecules show less sensitivity to the thermodynamic state of the system. Furthermore, our data indicate that a predominant number of H2O molecules exist as monomers without forming hydrogen bonds, particularly in a CO2-rich mixture, signaling a breakdown in the hydrogen bond network's orderliness, as evidenced by a marked decrease in tetrahedrality. These insights are crucial for a deeper atomic-level understanding of phase behaviors, contributing to the well-grounded design of CO2 injection under high-pressure and high-temperature conditions, where an atomic-scale perspective of the phase behavior is still lacking.
Collapse
Affiliation(s)
- Masashige Shiga
- Geological
Survey of Japan, National Institute of Advanced
Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8567, Japan
| | - Tetsuya Morishita
- Research
Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science
and Technology (AIST), Tsukuba, Ibaraki 305-8568, Japan
| | - Naoki Nishiyama
- Geological
Survey of Japan, National Institute of Advanced
Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8567, Japan
| | - Masao Sorai
- Geological
Survey of Japan, National Institute of Advanced
Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8567, Japan
| | - Masaatsu Aichi
- Department
of Environment Systems, Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba 277-8563, Japan
| | - Ayaka Abe
- Japan
Organization for Metals and Energy Security (JOGMEC), Minato-ku, Tokyo 105-0001, Japan
| |
Collapse
|
4
|
Coelho FM, Franco LFM, Firoozabadi A. Thermodiffusion of CO 2 in Water by Nonequilibrium Molecular Dynamics Simulations. J Phys Chem B 2023; 127:2749-2760. [PMID: 36930893 DOI: 10.1021/acs.jpcb.2c08260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023]
Abstract
The components of a fluid mixture may segregate due to the Soret effect, a coupling phenomenon in which mass flux can be induced by a thermal gradient. In this work, we evaluate systematically the thermodiffusion of the CO2-H2O mixture, and the influence of the geothermal gradient on CO2 segregation in deep saline aquifers in CO2 storage. The eHeX method, a nonequilibrium molecular dynamics simulation approach, is judiciously selected to simulate the phenomenon. At 350 K, 400 bar, and CO2 mole fraction of 0.02 (aquifer conditions), CO2 accumulates on the cold side, and the thermal diffusion factor is close to 1 in a number of force fields. The lower the temperature, the higher is the separation and the thermal diffusion factor. In colder regions, water self-association is stronger, whereas the CO2-H2O cross-association and the CO2-CO2 interactions enhance at higher temperatures. Thermodiffusion and gravitational segregation have opposite effects on CO2 segregation. At typical subsurface conditions, the Soret effect is more pronounced than gravity segregation, and CO2 concentrates in the top (colder region). Our work sets the stage to model the effect of electrolytes on CO2 segregation in subsurface aquifers and other areas of interest.
Collapse
Affiliation(s)
- Felipe M Coelho
- School of Chemical Engineering, University of Campinas, Campinas, SP 13083-852, Brazil.,Reservoir Engineering Research Institute (RERI), Palo Alto, California, 94306, United States
| | - Luís F M Franco
- School of Chemical Engineering, University of Campinas, Campinas, SP 13083-852, Brazil
| | - Abbas Firoozabadi
- Reservoir Engineering Research Institute (RERI), Palo Alto, California, 94306, United States
| |
Collapse
|
5
|
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]
|
6
|
Thermodynamics and Transport Properties of CBD and Δ9-THC: A first attempt using Molecular Dynamics. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.121048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
7
|
Algaba J, Acuña E, Míguez JM, Mendiboure B, Zerón IM, Blas FJ. Simulation of the carbon dioxide hydrate-water interfacial energy. J Colloid Interface Sci 2022; 623:354-367. [PMID: 35594594 DOI: 10.1016/j.jcis.2022.05.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 05/04/2022] [Accepted: 05/05/2022] [Indexed: 10/18/2022]
Abstract
HYPOTHESIS Carbon dioxide hydrates are ice-like nonstoichiometric inclusion solid compounds with importance to global climate change, and gas transportation and storage. The thermodynamic and kinetic mechanisms that control carbon dioxide nucleation critically depend on hydrate-water interfacial free energy. Only two independent indirect experiments are available in the literature. Interfacial energies show large uncertainties due to the conditions at which experiments are performed. Under these circumstances, we hypothesize that accurate molecular models for water and carbon dioxide combined with computer simulation tools can offer an alternative but complementary way to estimate interfacial energies at coexistence conditions from a molecular perspective. CALCULATIONS We have evaluated the interfacial free energy of carbon dioxide hydrates at coexistence conditions (three-phase equilibrium or dissociation line) implementing advanced computational methodologies, including the novel Mold Integration methodology. Our calculations are based on the definition of the interfacial free energy, standard statistical thermodynamic techniques, and the use of the most reliable and used molecular models for water (TIP4P/Ice) and carbon dioxide (TraPPE) available in the literature. FINDINGS We find that simulations provide an interfacial energy value, at coexistence conditions, consistent with the experiments from its thermodynamic definition. Our calculations are reliable since are based on the use of two molecular models that accurately predict: (1) The ice-water interfacial free energy; and (2) the dissociation line of carbon dioxide hydrates. Computer simulation predictions provide alternative but reliable estimates of the carbon dioxide interfacial energy. Our pioneering work demonstrates that is possible to predict interfacial energies of hydrates from a truly computational molecular perspective and opens a new door to the determination of free energies of hydrates.
Collapse
Affiliation(s)
- Jesús Algaba
- Department of Chemical Engineering, South Kensington Campus, Imperial College London, SW7 2AZ London, United Kingdom
| | - Esteban Acuña
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - José Manuel Míguez
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - Bruno Mendiboure
- Laboratoire des Fluides Complexes et Leurs Reservoirs, UMR5150, Universite de Pau et des Pays de l'Adour, B. P. 1155, Pau Cedex 64014, France
| | - Iván M Zerón
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - Felipe J Blas
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain.
| |
Collapse
|
8
|
Narayanan Nair AK, Anwari Che Ruslan MF, Ramirez Hincapie ML, Sun S. Bulk and Interfacial Properties of Brine or Alkane in the Presence of Carbon Dioxide, Methane, and Their Mixture. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c00249] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Arun Kumar Narayanan Nair
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| | - Mohd Fuad Anwari Che Ruslan
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| | - Marcia Luna Ramirez Hincapie
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| | - Shuyu Sun
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| |
Collapse
|
9
|
Yue S, Riera M, Ghosh R, Panagiotopoulos AZ, Paesani F. Transferability of data-driven, many-body models for CO2 simulations in the vapor and liquid phases. J Chem Phys 2022; 156:104503. [DOI: 10.1063/5.0080061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shuwen Yue
- Princeton University, United States of America
| | - Marc Riera
- Chemistry and Biochemistry, University of California San Diego Department of Chemistry and Biochemistry, United States of America
| | - Raja Ghosh
- University of California San Diego, United States of America
| | | | | |
Collapse
|
10
|
Mitev P, Briels WJ, Hermansson K. Space-Resolved OH Vibrational Spectra of the Hydration Shell around CO 2. J Phys Chem B 2021; 125:13886-13895. [PMID: 34927438 PMCID: PMC8724796 DOI: 10.1021/acs.jpcb.1c06123] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2021] [Revised: 11/01/2021] [Indexed: 11/30/2022]
Abstract
The CO2 molecule is weakly bound in water. Here we analyze the influence of a dissolved CO2 molecule on the structure and OH vibrational spectra of the surrounding water. From the analysis of ab initio molecular dynamics simulations (BLYP-D3) we present static (structure, coordination, H-bonding, tetrahedrality) and dynamical (OH vibrational spectra) properties of the water molecules as a function of distance from the solute. We find a weakly oscillatory variation ("ABBA") in the 'solution minus bulk water' spectrum. The origin of these features can largely be traced back to solvent-solute hard-core interactions which lead to variations in density and tetrahedrality when moving from the solute's vicinity out to the bulk region. The high-frequency peak in the solute-affected spectra is specifically analyzed and found to originate from both water OH groups that fulfill the geometric H-bond criteria, and from those that do not (dangling ones). Effectively, neither is hydrogen-bonded.
Collapse
Affiliation(s)
- Pavlin
D. Mitev
- Department
of Chemistry-Ångström Laboratory, Uppsala University, Box 538, S-751 21, Uppsala, Sweden
- Uppsala
Multidisciplinary Center for Advanced Computational Science, Uppsala University, Uppsala, SE-751 05, Sweden
| | - W. J. Briels
- MESA+
Institute for Nanotechnology, University
of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
- IBI-4, Forschungszentrum Jülich, D-52425 Jülich, Germany
| | - Kersti Hermansson
- Department
of Chemistry-Ångström Laboratory, Uppsala University, Box 538, S-751 21, Uppsala, Sweden
| |
Collapse
|
11
|
Heijmans K, Tranca IC, Smeulders DMJ, Vlugt TJH, Gaastra-Nedea SV. Gibbs Ensemble Monte Carlo for Reactive Force Fields to Determine the Vapor-Liquid Equilibrium of CO 2 and H 2O. J Chem Theory Comput 2021; 17:322-329. [PMID: 33350819 PMCID: PMC7808213 DOI: 10.1021/acs.jctc.0c00876] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
![]()
Absorption and reactive properties
of fluids in porous media are
key to the design and improvement of numerous energy related applications.
Molecular simulations of these systems require accurate force fields
that capture the involved chemical reactions and have the ability
to describe the vapor–liquid equilibrium (VLE). Two new reactive
force fields (ReaxFF) for CO2 and H2O are developed,
which are capable of not only modeling bond breaking and formation
in reactive environments but also predicting their VLEs at saturation
conditions. These new force fields include extra terms (ReaxFF-lg)
to improve the long-range interactions between the molecules. For
validation, we have developed a new Gibbs ensemble Monte Carlo (GEMC–ReaxFF)
approach to predict the VLE. Computed VLE data show good agreement
with National Institute of Standards and Technology reference data
as well as existing nonreactive force fields. This validation proves
the applicability of the GEMC–ReaxFF method to test new reactive
force fields, and simultaneously it proves the applicability to extend
newly developed ReaxFF force fields to other more complex reactive
systems.
Collapse
Affiliation(s)
- Koen Heijmans
- Department of Mechanical Engineering, Eindhoven University of Technology, Groene Loper 15, 5600MB Eindhoven, The Netherlands
| | - Ionut C Tranca
- Department of Mechanical Engineering, Eindhoven University of Technology, Groene Loper 15, 5600MB Eindhoven, The Netherlands
| | - David M J Smeulders
- Department of Mechanical Engineering, Eindhoven University of Technology, Groene Loper 15, 5600MB Eindhoven, The Netherlands
| | - Thijs J H Vlugt
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Silvia V Gaastra-Nedea
- Department of Mechanical Engineering, Eindhoven University of Technology, Groene Loper 15, 5600MB Eindhoven, The Netherlands
| |
Collapse
|
12
|
Zhang M, Chen L, Ma J. Multiscale Simulation of Vinyl Acetate Systems Applied in the Industrial Gas Separation Column. Ind Eng Chem Res 2020. [DOI: 10.1021/acs.iecr.0c03969] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Minhua Zhang
- Key Laboratory for Green Chemical Technology of Ministry of Education, R&D Center for Petrochemical Technology, Tianjin University, Tianjin 300072, China
- Collaborative Innovation Center of Chemical Science and Engineering, Tianjin University, Tianjin 300072, China
| | - Lihang Chen
- Key Laboratory for Green Chemical Technology of Ministry of Education, R&D Center for Petrochemical Technology, Tianjin University, Tianjin 300072, China
- Collaborative Innovation Center of Chemical Science and Engineering, Tianjin University, Tianjin 300072, China
| | - Jing Ma
- Key Laboratory for Green Chemical Technology of Ministry of Education, R&D Center for Petrochemical Technology, Tianjin University, Tianjin 300072, China
- Collaborative Innovation Center of Chemical Science and Engineering, Tianjin University, Tianjin 300072, China
| |
Collapse
|
13
|
Yang Y, Narayanan Nair AK, Anwari Che Ruslan MF, Sun S. Bulk and Interfacial Properties of the Decane + Water System in the Presence of Methane, Carbon Dioxide, and Their Mixture. J Phys Chem B 2020; 124:9556-9569. [PMID: 33059452 DOI: 10.1021/acs.jpcb.0c05759] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular dynamics simulations are carried out to study the two-phase behavior of the n-decane + water system in the presence of methane, carbon dioxide, and their mixture under reservoir conditions. The simulation studies were complemented by theoretical modeling using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state (EoS) and density gradient theory. Our results show that the presence of methane and carbon dioxide decreases the interfacial tension (IFT) of the decane + water system. In general, the IFT increases with increasing pressure and decreasing temperature for the methane + decane + water and carbon dioxide + decane + water systems, similar to what has been found for the corresponding decane + water system. The most important finding of this study is that the presence of carbon dioxide decreases the IFT of the methane + decane + water system. The atomic density profiles provide evidence of the local accumulation of methane and carbon dioxide at the interface, in most of the studied systems. The results of this study show the preferential dissolution in the water-rich phase and enrichment at the interface for carbon dioxide in the methane + carbon dioxide + decane + water system. This indicates the preferential interaction of water with carbon dioxide relative to methane and decane. Notably, there is an enrichment of the interface by decane at high mole fractions of methane in the methane/decane-rich or methane/carbon dioxide/decane-rich phase. Overall, the solubility of methane and carbon dioxide in the water-rich phase increases with increasing pressure and temperature. Additionally, we find that the overall performance of the PC-SAFT EoS and the cubic-plus-association EoS is similar with respect to the calculation of bulk and interfacial properties of these systems.
Collapse
Affiliation(s)
- Yafan Yang
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Arun Kumar Narayanan Nair
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Mohd Fuad Anwari Che Ruslan
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Shuyu Sun
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
14
|
Riera M, Yeh EP, Paesani F. Data-Driven Many-Body Models for Molecular Fluids: CO2/H2O Mixtures as a Case Study. J Chem Theory Comput 2020; 16:2246-2257. [DOI: 10.1021/acs.jctc.9b01175] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Eric P. Yeh
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
15
|
Fuentes-Azcatl R, Domínguez H. Prediction of experimental properties of CO2: improving actual force fields. J Mol Model 2019; 25:146. [DOI: 10.1007/s00894-019-4034-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 04/15/2019] [Indexed: 12/31/2022]
|
16
|
A review of molecular simulation applied in vapor-liquid equilibria (VLE) estimation of thermodynamic cycles. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.05.101] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
17
|
Yang Y, Narayanan Nair AK, Sun S. Molecular Dynamics Simulation Study of Carbon Dioxide, Methane, and Their Mixture in the Presence of Brine. J Phys Chem B 2017; 121:9688-9698. [DOI: 10.1021/acs.jpcb.7b08118] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Yafan Yang
- Physical Science and Engineering Division
(PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia
| | - Arun Kumar Narayanan Nair
- Physical Science and Engineering Division
(PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia
| | - Shuyu Sun
- Physical Science and Engineering Division
(PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia
| |
Collapse
|
18
|
13 The Role of Molecular Thermodynamics in Developing Industrial Processes and Novel Products That Meet the Needs for a Sustainable Future. ACTA ACUST UNITED AC 2017. [DOI: 10.1201/9781315153209-14] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
19
|
Waage MH, Vlugt TJH, Kjelstrup S. Phase Diagram of Methane and Carbon Dioxide Hydrates Computed by Monte Carlo Simulations. J Phys Chem B 2017; 121:7336-7350. [DOI: 10.1021/acs.jpcb.7b03071] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | - 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
|
20
|
Jiang H, Economou IG, Panagiotopoulos AZ. Molecular Modeling of Thermodynamic and Transport Properties for CO 2 and Aqueous Brines. Acc Chem Res 2017; 50:751-758. [PMID: 28234455 DOI: 10.1021/acs.accounts.6b00632] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular simulation techniques using classical force-fields occupy the space between ab initio quantum mechanical methods and phenomenological correlations. In particular, Monte Carlo and molecular dynamics algorithms can be used to provide quantitative predictions of thermodynamic and transport properties of fluids relevant for geologic carbon sequestration at conditions for which experimental data are uncertain or not available. These methods can cover time and length scales far exceeding those of quantum chemical methods, while maintaining transferability and predictive power lacking from phenomenological correlations. The accuracy of predictions depends sensitively on the quality of the molecular models used. Many existing fixed-point-charge models for water and aqueous mixtures fail to represent accurately these fluid properties, especially when descriptions covering broad ranges of thermodynamic conditions are needed. Recent work on development of accurate models for water, CO2, and dissolved salts, as well as their mixtures, is summarized in this Account. Polarizable models that can respond to the different dielectric environments in aqueous versus nonaqueous phases are necessary for predictions of properties over extended ranges of temperatures and pressures. Phase compositions and densities, activity coefficients of the dissolved salts, interfacial tensions, viscosities and diffusivities can be obtained in near-quantitative agreement to available experimental data, using relatively modest computational resources. In some cases, for example, for the composition of the CO2-rich phase in coexistence with an aqueous phase, recent results from molecular simulations have helped discriminate among conflicting experimental data sets. The sensitivity of properties on the quality of the intermolecular interaction model varies significantly. Properties such as the phase compositions or electrolyte activity coefficients are much more sensitive than phase densities, viscosities, or component diffusivities. Strong confinement effects on physical properties in nanoscale media can also be directly obtained from molecular simulations. Future work on molecular modeling for CO2 and aqueous brines is likely to be focused on more systematic generation of interaction models by utilizing quantum chemical as well as direct experimental measurements. New ion models need to be developed for use with the current generation of polarizable water models, including ion-ion interactions that will allow for accurate description of dense, mixed brines. Methods will need to be devised that go beyond the use of effective potentials for incorporation of quantum effects known to be important for water, and reactive force fields developed that can handle bond creation and breaking in systems with carbonate and silicate minerals. Another area of potential future work is the integration of molecular simulation methods in multiscale models for the chemical reactions leading to mineral dissolution and flow within the porous media in underground formations.
Collapse
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
21
|
Jiang H, Economou IG, Panagiotopoulos AZ. Phase Equilibria of Water/CO2 and Water/n-Alkane Mixtures from Polarizable Models. J Phys Chem B 2017; 121:1386-1395. [DOI: 10.1021/acs.jpcb.6b12791] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
22
|
Orozco GA, Lachet V, Mackie AD. Physical Absorption of Green House Gases in Amines: The Influence of Functionality, Structure, and Cross-Interactions. J Phys Chem B 2016; 120:13136-13143. [PMID: 27966955 DOI: 10.1021/acs.jpcb.6b09819] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Monte Carlo simulations were performed in the isothermal-isobaric ensemble (NPT) to calculate the Henry constants of methane (CH4), nitrous oxide (N2O), and carbon dioxide (CO2) in pure H2O, amines, and alkanolamines using the classical Lorentz-Berthelot combining rules (L-B). The Henry constants of N2O and CO2 in water are highly overestimated and motivated us to propose a new set of unlike interactions. Contrarily, the Henry constant of N2O in MEA is underestimated by around 40%, and again, a new reoptimized cross unlike parameter is able to reproduce the constant to within 10%. An analysis is given of the relationship between the physical absorption of these gases and the chemical structure or functionality of 12 molecules including amines and alkanolamines using the anisotropic united atom intermolecular potential (AUA4). Finally, the solubility of N2O in an aqueous solution of monoethanolamine (MEA) at 30% (wt) was also studied. A Henry constant within 7% of the experimental value was found by using the reoptimized parameters along with L-B to account for the MEA + H2O unlike interactions. This very good agreement without additional adjustments for the MEA + H2O system may be attributed to the good excess properties predictions found in previous works for the binary mixture (MEA + H2O). However, further work, including additional alkanolamines in aqueous solutions at several concentrations, is required to verify this particular point.
Collapse
Affiliation(s)
- Gustavo A Orozco
- Facultad de Ciencias, Departamento de Física, Universidad Antonio Nariño , Carrera 3 Este No. 47 A-15, 110211 Bogotá, Colombia
| | - Véronique Lachet
- IFP Energies nouvelles , 1-4 Avenue de Bois-Préau, 92852 Rueil-Malmaison, France
| | - Allan D Mackie
- Departament d'Enginyeria Química, ETSEQ, Universitat Rovira i Virgili , Av. dels Països Catalans 26, 43007 Tarragona, Spain
| |
Collapse
|
23
|
Jiang H, Moultos OA, Economou IG, Panagiotopoulos AZ. Gaussian-Charge Polarizable and Nonpolarizable Models for CO2. J Phys Chem B 2016; 120:984-94. [DOI: 10.1021/acs.jpcb.5b11701] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
24
|
Floudas CA, Niziolek AM, Onel O, Matthews LR. Multi‐scale systems engineering for energy and the environment: Challenges and opportunities. AIChE J 2016. [DOI: 10.1002/aic.15151] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Christodoulos A. Floudas
- Artie McFerrin Dept. of Chemical EngineeringTexas A&M UniversityCollege Station TX77843 USA
- Texas A&M Energy Institute, 302D Williams Administration Building, 3372 Texas A&M UniversityCollege Station TX77843USA
| | - Alexander M. Niziolek
- Dept. of Chemical and Biological EngineeringPrinceton UniversityPrinceton NJ08544 USA
- Artie McFerrin Dept. of Chemical EngineeringTexas A&M UniversityCollege Station TX77843 USA
- Texas A&M Energy Institute, 302D Williams Administration Building, 3372 Texas A&M UniversityCollege Station TX77843USA
| | - Onur Onel
- Dept. of Chemical and Biological EngineeringPrinceton UniversityPrinceton NJ08544 USA
- Artie McFerrin Dept. of Chemical EngineeringTexas A&M UniversityCollege Station TX77843 USA
- Texas A&M Energy Institute, 302D Williams Administration Building, 3372 Texas A&M UniversityCollege Station TX77843USA
| | - Logan R. Matthews
- Dept. of Chemical and Biological EngineeringPrinceton UniversityPrinceton NJ08544 USA
- Artie McFerrin Dept. of Chemical EngineeringTexas A&M UniversityCollege Station TX77843 USA
- Texas A&M Energy Institute, 302D Williams Administration Building, 3372 Texas A&M UniversityCollege Station TX77843USA
| |
Collapse
|
25
|
Costandy J, Michalis VK, Tsimpanogiannis IN, Stubos AK, Economou IG. The role of intermolecular interactions in the prediction of the phase equilibria of carbon dioxide hydrates. J Chem Phys 2015; 143:094506. [DOI: 10.1063/1.4929805] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joseph Costandy
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Vasileios K. Michalis
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos,” GR-15310 Aghia Paraskevi Attikis, Greece
| | - Athanassios K. Stubos
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos,” GR-15310 Aghia Paraskevi Attikis, Greece
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| |
Collapse
|
26
|
Moultos OA, Orozco GA, Tsimpanogiannis IN, Panagiotopoulos AZ, Economou IG. Atomistic molecular dynamics simulations of H2O diffusivity in liquid and supercritical CO2. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1023224] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Othonas A. Moultos
- Chemical Engineering Program , Texas A&M University at Qatar , Doha, Qatar
| | - Gustavo A. Orozco
- Department of Chemical and Biological Engineering, Princeton University , Princeton, United States
| | - Ioannis N. Tsimpanogiannis
- Chemical Engineering Program , Texas A&M University at Qatar , Doha, Qatar
- Environmental Research Laboratory, Institute of Nuclear and Radiology Sciences and Technology, Energy and Safety, National Center for Scientific Research “Demokritos” , Aghia Paraskevi, Greece
| | | | | |
Collapse
|
27
|
Orozco GA, Moultos OA, Jiang H, Economou IG, Panagiotopoulos AZ. Molecular simulation of thermodynamic and transport properties for the H2O+NaCl system. J Chem Phys 2014; 141:234507. [DOI: 10.1063/1.4903928] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Gustavo A. Orozco
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Hao Jiang
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|