1
|
Raju D, Ramdin M, Vlugt TJH. Thermophysical Properties and Phase Behavior of CO 2 with Impurities: Insight from Molecular Simulations. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2024; 69:2735-2755. [PMID: 39139986 PMCID: PMC11318637 DOI: 10.1021/acs.jced.4c00268] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/13/2024] [Revised: 06/19/2024] [Accepted: 07/02/2024] [Indexed: 08/15/2024]
Abstract
Experimentally determining thermophysical properties for various compositions commonly found in CO2 transportation systems is extremely challenging. To overcome this challenge, we performed Monte Carlo (MC) and Molecular Dynamics (MD) simulations of CO2 rich mixtures to compute thermophysical properties such as densities, thermal expansion coefficients, isothermal compressibilities, heat capacities, Joule-Thomson coefficients, speed of sound, and viscosities at temperatures of (235-313) K and pressures of (20-200) bar. We computed thermophysical properties of pure CO2 and CO2 rich mixtures with N2, Ar, H2, and CH4 as impurities of (1-10) mol % and showed good agreement with available Equations of State (EoS). We showed that impurities decrease the values of thermal expansion coefficients, isothermal compressibilities, heat capacities, and Joule-Thomson coefficients in the gas phase, while these values increase in the liquid and supercritical phases. In contrast, impurities increase the value of speed of sound in the gas phase and decrease it in the liquid and supercritical phases. We present an extensive data set of thermophysical properties for CO2 rich mixtures with various impurities, which will help to design the safe and efficient operation of CO2 transportation systems.
Collapse
Affiliation(s)
- D. Raju
- Engineering Thermodynamics, Process & Energy Department, Faculty of
Mechanical Engineering, Delft University of Technology,
Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of
Mechanical Engineering, Delft University of Technology,
Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of
Mechanical Engineering, Delft University of Technology,
Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| |
Collapse
|
2
|
Broad J, Wheatley RJ, Graham RS. Parallel Implementation of Nonadditive Gaussian Process Potentials for Monte Carlo Simulations. J Chem Theory Comput 2023. [PMID: 37368843 DOI: 10.1021/acs.jctc.3c00113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
A strategy is presented to implement Gaussian process potentials in molecular simulations through parallel programming. Attention is focused on the three-body nonadditive energy, though all algorithms extend straightforwardly to the additive energy. The method to distribute pairs and triplets between processes is general to all potentials. Results are presented for a simulation box of argon, including full box and atom displacement calculations, which are relevant to Monte Carlo simulation. Data on speed-up are presented for up to 120 processes across four nodes. A 4-fold speed-up is observed over five processes, extending to 20-fold over 40 processes and 30-fold over 120 processes.
Collapse
Affiliation(s)
- Jack Broad
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Richard J Wheatley
- School of Chemistry, University of Nottingham, Nottingham, NG7 2RD, England
| | - Richard S Graham
- School of Mathematical Sciences, University of Nottingham, Nottingham, NG7 2RD, England
| |
Collapse
|
3
|
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]
|
4
|
Graham RS, Wheatley RJ. Machine learning for non-additive intermolecular potentials: quantum chemistry to first-principles predictions. Chem Commun (Camb) 2022; 58:6898-6901. [PMID: 35642644 DOI: 10.1039/d2cc01820a] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Prediction of thermophysical properties from molecular principles requires accurate potential energy surfaces (PES). We present a widely-applicable method to produce first-principles PES from quantum chemistry calculations. Our approach accurately interpolates three-body non-additive interaction data, using the machine learning technique, Gaussian Processes (GP). The GP approach needs no bespoke modification when the number or type of molecules is changed. Our method produces highly accurate interpolation from significantly fewer training points than typical approaches, meaning ab initio calculations can be performed at higher accuracy. As an exemplar we compute the PES for all three-body cross interactions for CO2-Ar mixtures. From these we calculate the CO2-Ar virial coefficients up to 5th order. The resulting virial equation of state (EoS) is convergent for densities up to the critical density. Where convergent, the EoS makes accurate first-principles predictions for a range of thermophysical properties for CO2-Ar mixtures, including the compressibility factor, speed of sound and Joule-Thomson coefficient. Our method has great potential to make wide-ranging first-principles predictions for mixtures of comparably sized molecules. Such predictions can replace the need for expensive, laborious and repetitive experiments and inform the continuum models required for applications.
Collapse
Affiliation(s)
- Richard S Graham
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, UK.
| | | |
Collapse
|
5
|
|
6
|
Broad J, Preston S, Wheatley RJ, Graham RS. Gaussian process models of potential energy surfaces with boundary optimization. J Chem Phys 2021; 155:144106. [PMID: 34654292 DOI: 10.1063/5.0063534] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A strategy is outlined to reduce the number of training points required to model intermolecular potentials using Gaussian processes, without reducing accuracy. An asymptotic function is used at a long range, and the crossover distance between this model and the Gaussian process is learnt from the training data. The results are presented for different implementations of this procedure, known as boundary optimization, across the following dimer systems: CO-Ne, HF-Ne, HF-Na+, CO2-Ne, and (CO2)2. The technique reduces the number of training points, at fixed accuracy, by up to ∼49%, compared to our previous work based on a sequential learning technique. The approach is readily transferable to other statistical methods of prediction or modeling problems.
Collapse
Affiliation(s)
- Jack Broad
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Simon Preston
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Richard J Wheatley
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Richard S Graham
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|
7
|
Liu Q, Liu L, An F, Huang J, Zhou Y, Xie D. A full-dimensional ab initio intermolecular potential energy surface and rovibrational spectra for OC-HF and OC-DF. J Chem Phys 2021; 155:084302. [PMID: 34470366 DOI: 10.1063/5.0061291] [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
We present a full-dimensional ab initio intermolecular potential energy surface (IPES) for the OC-HF van der Waals complex. 3167 ab initio points were computed at the frozen-core (FC) explicitly correlated coupled cluster [FC-CCSD(T)-F12b] level, with the augmented correlation-consistent polarized valence quadruple-zeta basis set plus bond functions. Basis set superposition error correction was also considered by the full counterpoise procedure. Gaussian process regression (GPR) was used to map out the potential energy surface, while a multipole expansion method was employed to smooth the ab initio noise of intermolecular potential in the long range. The global minimum of -1248.364 cm-1 was located at the linear configuration with the C atom pointing toward the H atom of the HF molecule. In addition, a local minimum of -602.026 cm-1 was found at another linear configuration with the O atom pointing toward the H atom of the HF molecule. The eigenstates were calculated on the vibrational averaged four-dimensional IPESs with the mixed radial discrete variable representation/angular finite basis representation method and Lanczos propagation algorithm. The dissociation energy D0 was calculated to be 701.827 cm-1, well reproducing the experimental value of 732 ± 2 cm-1. The dipole moment surfaces were also fitted by GPR from 3132 ab initio points calculated using the coupled cluster method [CCSD(T)] with AVTZ basis set plus bond functions. The frequencies and relative line intensities of rovibrational transitions in the HF (DF) and CO stretching bands were further calculated and compared well with the experimental results. These results indicate the high fidelity of the new IPES.
Collapse
Affiliation(s)
- Qiong Liu
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Lu Liu
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Feng An
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Jing Huang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yanzi Zhou
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
8
|
Deringer VL, Bartók AP, Bernstein N, Wilkins DM, Ceriotti M, Csányi G. Gaussian Process Regression for Materials and Molecules. Chem Rev 2021; 121:10073-10141. [PMID: 34398616 PMCID: PMC8391963 DOI: 10.1021/acs.chemrev.1c00022] [Citation(s) in RCA: 245] [Impact Index Per Article: 81.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Indexed: 12/18/2022]
Abstract
We provide an introduction to Gaussian process regression (GPR) machine-learning methods in computational materials science and chemistry. The focus of the present review is on the regression of atomistic properties: in particular, on the construction of interatomic potentials, or force fields, in the Gaussian Approximation Potential (GAP) framework; beyond this, we also discuss the fitting of arbitrary scalar, vectorial, and tensorial quantities. Methodological aspects of reference data generation, representation, and regression, as well as the question of how a data-driven model may be validated, are reviewed and critically discussed. A survey of applications to a variety of research questions in chemistry and materials science illustrates the rapid growth in the field. A vision is outlined for the development of the methodology in the years to come.
Collapse
Affiliation(s)
- Volker L. Deringer
- Department
of Chemistry, Inorganic Chemistry Laboratory, University of Oxford, Oxford OX1 3QR, United Kingdom
| | - Albert P. Bartók
- Department
of Physics and Warwick Centre for Predictive Modelling, School of
Engineering, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Noam Bernstein
- Center
for Computational Materials Science, U.S.
Naval Research Laboratory, Washington D.C. 20375, United States
| | - David M. Wilkins
- Atomistic
Simulation Centre, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - Michele Ceriotti
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, Lausanne, Switzerland
| | - Gábor Csányi
- Engineering
Laboratory, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
9
|
Uteva E, Graham RS, Wilkinson RD, Wheatley RJ. Active learning in Gaussian process interpolation of potential energy surfaces. J Chem Phys 2018; 149:174114. [DOI: 10.1063/1.5051772] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Elena Uteva
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Richard S. Graham
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Richard D. Wilkinson
- School of Mathematics and Statistics, University of Sheffield, Sheffield S10 2TN, United Kingdom
| | - Richard J. Wheatley
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|