1
|
Drisko CR, Gezelter JD. A Reverse Nonequilibrium Molecular Dynamics Algorithm for Coupled Mass and Heat Transport in Mixtures. J Chem Theory Comput 2024; 20:4986-4997. [PMID: 38833377 DOI: 10.1021/acs.jctc.4c00182] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024]
Abstract
We present a new method for introducing stable nonequilibrium concentration gradients in molecular dynamics simulations of mixtures. This method extends earlier reverse nonequilibrium molecular dynamics (RNEMD) methods, which use kinetic energy scaling moves to create temperature or velocity gradients. In the new scaled particle flux (SPF-RNEMD) algorithm, energies and forces are computed simultaneously for a molecule existing in two nonadjacent regions of a simulation box, and the system evolves under a linear combination of these interactions. A continuously increasing particle scaling variable is responsible for the migration of the molecule between the regions as the simulation progresses, allowing for simulations under an applied particle flux. To test the method, we investigate diffusivity in mixtures of identical but distinguishable particles and in a simple mixture of multiple Lennard-Jones particles. The resulting concentration gradients provide Fick diffusion constants for mixtures. We also discuss using the new method to obtain coupled transport properties using simultaneous particle and thermal fluxes to compute the temperature dependence of the diffusion coefficient and activation energies for diffusion from a single simulation. Lastly, we demonstrate the use of this new method in interfacial systems by computing the diffusive permeability of a molecular fluid moving through a nanoporous graphene membrane.
Collapse
Affiliation(s)
- Cody R Drisko
- Department of Chemistry and Biochemistry, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| | - J Daniel Gezelter
- Department of Chemistry and Biochemistry, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| |
Collapse
|
2
|
Hao X, Li C, Meng Q, Sun J, Huang L, Bu Q, Li C. Molecular Dynamics Simulation of the Three-Phase Equilibrium Line of CO 2 Hydrate with OPC Water Model. ACS OMEGA 2023; 8:39847-39854. [PMID: 37901483 PMCID: PMC10601413 DOI: 10.1021/acsomega.3c05673] [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: 08/16/2023] [Accepted: 09/26/2023] [Indexed: 10/31/2023]
Abstract
The three-phase coexistence line of the CO2 hydrate was determined using molecular dynamics (MD) simulations. By using the classical and modified Lorentz-Berthelot (LB) parameters, the simulations were carried out at 10 different pressures from 3 to 500 MPa. For the OPC water model, simulations with the classic and the modified LB parameters both showed negative deviations from the experimental values. For the TIP4P/Ice water model, good agreement with experimental equilibrium data can be achieved when the LB parameter is adjusted based on the solubility of CO2 in water. Our results also show that the influence of the water model on the equilibrium prediction is much larger than the CO2 model. Current simulations indicated that the H2O-H2O and H2O-CO2 cross-interactions' parameters might contribute equally to the accurate prediction of T3. According to our simulations, the prediction of T3 values showed relatively higher accuracy while using the combination of TIP4P/Ice water and EPM2 CO2 with modified LB parameter. Furthermore, varied χ values are recommended for accurate T3 estimation over a wide pressure range. The knowledge obtained in this study will be helpful for further accurate MD simulation of the process of CO2/CH4 replacement.
Collapse
Affiliation(s)
- Xiluo Hao
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Chengfeng Li
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingguo Meng
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Jianye Sun
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Li Huang
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingtao Bu
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Congying Li
- Center
of Deep Sea Research, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| |
Collapse
|
3
|
Hao X, Li C, Liu C, Meng Q, Sun J. The performance of OPC water model in prediction of the phase equilibria of methane hydrate. J Chem Phys 2022; 157:014504. [DOI: 10.1063/5.0093659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular dynamics (MD) simulations were performed to determine the three-phase coexistence line of sI methane hydrates. The MD simulations were carried out at four different pressures (4, 10, 40 and 100 MPa) by using direct phase coexistence method. In current simulations, water was described by either TIP4P/Ice or OPC models and methane was described as a simple Lennard-Jones (LJ) interaction site. Lorentz-Berthelot combining rules were used to calculate the parameters of the cross interactions. For OPC model, positive deviations from the energetic Lorentz-Berthelot rule were also considered based on the solubility of methane in water. For TIP4P/Ice water model, the obtained three phase coexistence temperatures showed good agreement with experiment data at higher pressures, which is consistent with previous predictions. For OPC water model, simulations using the classic and the modified LB parameters both showed negative deviations to the experimental values. Our results also indicated that the deviation of the T3 prediction by OPC model not much correlated with the predicted melting point of ice. At 4 MPa, the modified OPC model showed outstanding prediction of hydrate equilibrium temperature, even better than the prediction by TIP4P/Ice. The relative higher accuracy in biomolecular MD of OPC model suggests that this model may have a better performance in hydrate MD simulations of biomolecule-based additives.
Collapse
Affiliation(s)
- Xiluo Hao
- Qingdao Institute of Marine Geology, China
| | | | | | | | - Jianye Sun
- Qingdao Institute of Marine Geology, China
| |
Collapse
|
4
|
Zhang Z, Qi N, Wu Y, Chen Z. Pressure-Induced Enhancement of Thermoelectric Performance in Rubrene. ACS APPLIED MATERIALS & INTERFACES 2021; 13:44409-44417. [PMID: 34515463 DOI: 10.1021/acsami.1c12832] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
In this work, the thermoelectric performance of a typical small-molecule organic semiconductor rubrene under different hydrostatic pressures was studied by first-principles calculation and molecular dynamics simulation. The ZT value of rubrene can reach 1.6 at 400 K due to an unprecedented increase in hole mobility under hydrostatic pressure. The underlying mechanism is ascribed to the suppression of low-frequency phonons (which weakens electron-phonon scattering) and the increase in the intermolecular electronic coupling. The effect of uniaxial stress has also been investigated to confirm this conclusion. Our results provide meaningful insights to understand the relationship between thermoelectric properties and hydrostatic pressure in organic semiconductors.
Collapse
Affiliation(s)
- Ziye Zhang
- Hubei Nuclear Solid Physics Key Laboratory, Department of Physics, Wuhan University, Wuhan 430072, China
| | - Ning Qi
- Hubei Nuclear Solid Physics Key Laboratory, Department of Physics, Wuhan University, Wuhan 430072, China
| | - Yichu Wu
- Hubei Nuclear Solid Physics Key Laboratory, Department of Physics, Wuhan University, Wuhan 430072, China
| | - Zhiquan Chen
- Hubei Nuclear Solid Physics Key Laboratory, Department of Physics, Wuhan University, Wuhan 430072, China
| |
Collapse
|
5
|
Gittus OR, Bresme F. Thermophysical properties of water using reactive force fields. J Chem Phys 2021; 155:114501. [PMID: 34551553 DOI: 10.1063/5.0057868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.
Collapse
Affiliation(s)
- Oliver R Gittus
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| |
Collapse
|
6
|
Celebi AT, Vlugt TJH, Moultos OA. Thermal conductivity of aqueous solutions of reline, ethaline, and glyceline deep eutectic solvents; a molecular dynamics simulation study. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1876263] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Alper T. Celebi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
7
|
Alosious S, Kannam SK, Sathian SP, Todd BD. Kapitza resistance at water-graphene interfaces. J Chem Phys 2020; 152:224703. [PMID: 32534537 DOI: 10.1063/5.0009001] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Heat transfer across fluid-solid interfaces in nanoconfinement has received significant attention due to its relevance in nanoscale systems. In this study, we investigate the Kapitza resistance at the water-graphene interface with the help of classical molecular dynamics simulation techniques in conjunction with our recently proposed equilibrium molecular dynamics (EMD) method [S. Alosious et al., J. Chem. Phys. 151, 194502 (2019)]. The size effect of the Kapitza resistance on different factors such as the number of graphene layers, the cross-sectional area, and the width of the water block was studied. The Kapitza resistance decreases slightly with an increase in the number of layers, while the influence of the cross-sectional area and the width of the water block is negligible. The variation in the Kapitza resistance as a function of the number of graphene layers is attributed to the large phonon mean free path along the graphene cross-plane. An optimum water-graphene system, which is independent of size effects, was selected, and the same was used to determine the Kapitza resistance using the predicted EMD method. The values obtained from both the EMD and the non-equilibrium molecular dynamics (NEMD) methods were compared for different potentials and water models, and the results are shown to be in good agreement. Our method allows us to compute the Kapitza resistance using EMD simulations, which obviates the need to create a large temperature gradient required for the NEMD method.
Collapse
Affiliation(s)
- Sobin Alosious
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sridhar Kumar Kannam
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| | - Sarith P Sathian
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - B D Todd
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| |
Collapse
|
8
|
Li J, Wang ZL. Fluctuation-dissipation analysis of nonequilibrium thermal transport at the hydrate dissociation interface. Phys Chem Chem Phys 2019; 21:23492-23500. [PMID: 31617505 DOI: 10.1039/c9cp04780h] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Herein, a nonequilibrium molecular-dynamics simulation in an NVE ensemble was performed to investigate the spontaneous dissociation of methane-hydrate when it came in contact with liquid water. The nonequilibrium in the interface region is linked to the dissociation process of the hydrate near the interface according to the Onsager's hypothesis. The simulated thickness of the interface was found to be close to the acoustic phonon mean path of methane hydrate and agreed with the reference value. The normalized heat flow autocorrelation function was introduced to study fluctuation-dissipation in terms of the thickness and moving velocity of the interface and the Stefan number. This helped to clearly identify three distinct hydrate-decomposition regimes dominated by sensible heat, latent heat and an intrinsically unstable lattice framework. It was found that the fluctuation-dissipation theory could express the nonequilibrium nature in the front two stages before the threshold was reached, and the dissociation rate increased in the latter stage; this was different from the case of thermal-driven dissociation. The Stefan number decreased rapidly with dissociation in the initial stage and then fluctuated in the intermediate stage; this was analogous to the fluctuation characteristics of the heat flow autocorrelation function. The Stefan number effect shows that thermal dissipation drives the hydrate dissociation and correlates fluctuation to the nonequilibrium nature. It was also found that a small Stefan number was enough to break up the residual hydrate soon after the threshold was achieved. The transient interfacial thermal resistance of the interfacial region was obtained as a typical value in the range of 10-7-10-9 m2 K W-1. This justifies that fluctuation-dissipation exists in the nonequilibrium process of hydrate dissociation either in terms of heat flux, as observed in this study, or the diffusion of guest molecules, as reported in other studies.
Collapse
Affiliation(s)
- Jia Li
- Department of Energy and Power Engineering, China University of Petroleum, Qingdao 266580, China.
| | | |
Collapse
|
9
|
Striolo A. Clathrate hydrates: recent advances on CH4 and CO2 hydrates, and possible new frontiers. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1646436] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Alberto Striolo
- Department of Chemical Engineering, University College London, London, UK
| |
Collapse
|
10
|
Energy and Environmental Analysis of Membrane-Based CH4-CO2 Replacement Processes in Natural Gas Hydrates. ENERGIES 2019. [DOI: 10.3390/en12050850] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Natural gas hydrates are the largest reservoir of natural gas worldwide. This paper proposes and analyzes the CH4-CO2 replacement in the hydrate phase and pure methane collection through the use of membrane-based separation. The investigation uses a 1 L lab reactor, in which the CH4 hydrates are formed in a quartz sand matrix partially saturated with water. CH4 is subsequently dissociated with a CO2 stream supplied within the sediment inside the reactor. An energy and environmental analysis was carried out to prove the sustainability of the process. Results show that the process energy consumption constitutes 4.75% of the energy stored in the recovered methane. The carbon footprint of the CH4-CO2 exchange process is calculated as a balance of the CO2 produced in the process and the CO2 stored in system. Results provide an estimated negative value, equal to 0.004 moles sequestrated, thus proving the environmental benefit of the exchange process.
Collapse
|
11
|
Ghafari H, Mohammadi-Manesh H. The thermal properties of binary structure sI clathrate hydrate from molecular dynamics simulation. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1572142] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Hakime Ghafari
- Department of Chemistry, Faculty of Science, Yazd University, Yazd, Iran
| | | |
Collapse
|
12
|
Xu J, Chen Z, Liu J, Sun Z, Wang X, Zhang J. A molecular dynamic study on the dissociation mechanism of SI methane hydrate in inorganic salt aqueous solutions. J Mol Graph Model 2017; 75:403-412. [PMID: 28666231 DOI: 10.1016/j.jmgm.2017.03.022] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2016] [Revised: 03/28/2017] [Accepted: 03/29/2017] [Indexed: 10/19/2022]
Abstract
Gas hydrate is not only a potential energy resource, but also almost the biggest challenge in oil/gas flow assurance. Inorganic salts such as NaCl, KCl and CaCl2 are widely used as the thermodynamic inhibitor to reduce the risk caused by hydrate formation. However, the inhibition mechanism is still unclear. Therefore, molecular dynamic (MD) simulation was performed to study the dissociation of structure I (SI) methane hydrate in existence of inorganic salt aqueous solution on a micro-scale. The simulation results showed that, the dissociation became stagnant due to the presence of liquid film formed by the decomposed water molecules, and more inorganic ions could shorten the stagnation-time. The diffusion coefficients of ions and water molecules were the largest in KCl system. The structures of ion/H2O and H2O/H2O were the most compact in hydrate/NaCl system. The ionic ability to decompose hydrate cells followed the sequence of: Ca2+>2K+>2Cl->2Na+.
Collapse
Affiliation(s)
- Jiafang Xu
- University of Petroleum (East China), Qingdao, China; National Engineering Laboratory for Subsea Equipment Testing and Detection Technology, China.
| | - Zhe Chen
- University of Petroleum (East China), Qingdao, China
| | | | - Zening Sun
- China United Coalbed Methane Corporation, Ltd., Taiyuan, China.
| | - Xiaopu Wang
- University of Petroleum (East China), Qingdao, China
| | - Jun Zhang
- University of Petroleum (East China), Qingdao, China.
| |
Collapse
|
13
|
Burnham CJ, English NJ. Study of clathrate hydrates via equilibrium molecular-dynamics simulation employing polarisable and non-polarisable, rigid and flexible water models. J Chem Phys 2016; 144:164503. [PMID: 27131553 DOI: 10.1063/1.4947039] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christian J. Burnham
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
14
|
Weak phonon scattering effect of twin boundaries on thermal transmission. Sci Rep 2016; 6:19575. [PMID: 26822675 PMCID: PMC4731792 DOI: 10.1038/srep19575] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2015] [Accepted: 12/15/2015] [Indexed: 11/20/2022] Open
Abstract
To study the effect of twin boundaries on thermal transmission, thermal conductivities of twinned diamond with different twin thicknesses have been studied by NEMD simulation. Results indicate that twin boundaries show a weak phonon scattering effect on thermal transmission, which is only caused by the additional twin boundaries’ thermal resistance. Moreover, according to phonon kinetic theory, this weak phonon scattering effect of twin boundaries is mainly caused by a slightly reduced average group velocity.
Collapse
|
15
|
Gillan MJ, Alfè D, Manby FR. Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations. J Chem Phys 2015; 143:102812. [PMID: 26374005 DOI: 10.1063/1.4926444] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.
Collapse
Affiliation(s)
- M J Gillan
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - D Alfè
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - F R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
16
|
Xu J, Gu T, Sun Z, Li X, Wang X. Molecular dynamics study on the dissociation of methane hydrate via inorganic salts. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1081708] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
17
|
Michalis VK, Costandy J, Tsimpanogiannis IN, Stubos AK, Economou IG. Prediction of the phase equilibria of methane hydrates using the direct phase coexistence methodology. J Chem Phys 2015; 142:044501. [DOI: 10.1063/1.4905572] [Citation(s) in RCA: 90] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Vasileios K. Michalis
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23847, Doha, Qatar
| | - Joseph Costandy
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23847, Doha, Qatar
| | - Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research NCSR “Demokritos,” Aghia Paraskevi, Attiki GR-15310, Greece
| | - Athanassios K. Stubos
- Environmental Research Laboratory, National Center for Scientific Research NCSR “Demokritos,” Aghia Paraskevi, Attiki GR-15310, Greece
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23847, Doha, Qatar
| |
Collapse
|
18
|
Ning FL, Glavatskiy K, Ji Z, Kjelstrup S, H. Vlugt TJ. Compressibility, thermal expansion coefficient and heat capacity of CH4 and CO2 hydrate mixtures using molecular dynamics simulations. Phys Chem Chem Phys 2015; 17:2869-83. [DOI: 10.1039/c4cp04212c] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Understanding the thermal and mechanical properties of CH4 and CO2 hydrates is essential for the replacement of CH4 with CO2 in natural hydrate deposits as well as for CO2 sequestration and storage.
Collapse
Affiliation(s)
- F. L. Ning
- Faculty of Engineering
- China University of Geosciences
- Wuhan
- China
| | - K. Glavatskiy
- School of Applied Sciences
- RMIT University
- Melbourne VIC 3001
- Australia
| | - Z. Ji
- Faculty of Materials Science and Chemistry
- China University of Geosciences
- Wuhan
- China
| | - S. Kjelstrup
- Department of Chemistry
- Norwegian University of Science and Technology
- 7491-Trondheim
- Norway
- Process & Energy Laboratory
| | - T. J. H. Vlugt
- Process & Energy Laboratory
- Delft University of Technology
- Delft
- The Netherlands
| |
Collapse
|
19
|
Perspectives on molecular simulation of clathrate hydrates: Progress, prospects and challenges. Chem Eng Sci 2015. [DOI: 10.1016/j.ces.2014.07.047] [Citation(s) in RCA: 144] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
20
|
Abstract
Equilibrium molecular dynamics simulations that use the Green–Kubo method for sI CO2-hydrate systems from medium to full occupancy were performed to estimate the corresponding thermal conductivities at temperatures that range from 233.15K to 278.15K and pressures that range from 3MPa to 100MPa. Specific potential models for water and CO2were adopted. The effects of guest occupancy ratios and outside thermobaric conditions on CO2hydrate thermal conductivity were studied. The thermal mechanism was also analyzed. The thermal conductivities of hydrates of CH4, C2H6, N2, and O2were estimated. The size ratio of guest diameter to cavity diameter provided an adequate basis for understanding the thermal conductivities of gas hydrates.
Collapse
|
21
|
Stocker KM, Gezelter JD. A Method for Creating Thermal and Angular Momentum Fluxes in Nonperiodic Simulations. J Chem Theory Comput 2014; 10:1878-86. [PMID: 26580518 DOI: 10.1021/ct500221u] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We present a new reverse nonequilibrium molecular dynamics method that can be used with nonperiodic simulation cells. This method applies thermal and/or angular momentum fluxes between two arbitrary regions of the simulation and is capable of creating stable temperature and angular velocity gradients while conserving total energy and angular momentum. One particularly useful application is the exchange of kinetic energy between two concentric spherical regions, which can be used to generate thermal transport between nanoparticles and the solvent that surrounds them. The rotational couple to the solvent (a measure of interfacial friction) is also available via this method. As tests of the new method, we have computed the thermal conductivities of gold nanoparticles and water clusters, the interfacial thermal conductivity (G) of a solvated gold nanoparticle, and the interfacial friction of a variety of solvated gold nanostructures.
Collapse
Affiliation(s)
- Kelsey M Stocker
- Department of Chemistry and Biochemistry, University of Notre Dame , 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| | - J Daniel Gezelter
- Department of Chemistry and Biochemistry, University of Notre Dame , 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| |
Collapse
|
22
|
Deible MJ, Tuguldur O, Jordan KD. Theoretical Study of the Binding Energy of a Methane Molecule in a (H2O)20 Dodecahedral Cage. J Phys Chem B 2014; 118:8257-63. [DOI: 10.1021/jp501592h] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Michael J. Deible
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Odbadrakh Tuguldur
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
23
|
Xu Y, Leitner DM. Communication maps of vibrational energy transport through Photoactive Yellow Protein. J Phys Chem A 2014; 118:7280-7. [PMID: 24552496 DOI: 10.1021/jp411281y] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We calculate communication maps for Photoactive Yellow Protein (PYP) from the purple phototropic eubacterium Halorhodospira halophile and use them to elucidate energy transfer pathways from the chromophore through the rest of the protein in the ground and excited state. The calculations reveal that in PYP excess energy from the chromophore flows mainly to regions of the surrounding residues that hydrogen bond to the chromophore. In addition, quantum mechanics/molecular mechanics and molecular dynamics (MD) simulations of the dielectric response of the protein and solvent environment due to charge rearrangement on the chromophore following photoexcitation are also presented, with both approaches yielding similar time constants for the response. Results of MD simulations indicate that the residues hydrogen bonding to the chromophore make the largest contribution to the response.
Collapse
Affiliation(s)
- Yao Xu
- Department of Chemistry and Chemical Physics Program, University of Nevada , Reno, Nevada 89557, United States
| | | |
Collapse
|
24
|
Review of CO2–CH4 clathrate hydrate replacement reaction laboratory studies – Properties and kinetics. J Taiwan Inst Chem Eng 2013. [DOI: 10.1016/j.jtice.2013.03.010] [Citation(s) in RCA: 90] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
25
|
Leitner DM. Thermal Boundary Conductance and Thermal Rectification in Molecules. J Phys Chem B 2013; 117:12820-8. [DOI: 10.1021/jp402012z] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David M. Leitner
- Department of Chemistry and Chemical
Physics Program, University of Nevada,
Reno, Nevada 89557, United States
| |
Collapse
|
26
|
Sirk TW, Moore S, Brown EF. Characteristics of thermal conductivity in classical water models. J Chem Phys 2013; 138:064505. [DOI: 10.1063/1.4789961] [Citation(s) in RCA: 98] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Chen J, Wang D, Shuai Z. First-Principles Predictions of Thermoelectric Figure of Merit for Organic Materials: Deformation Potential Approximation. J Chem Theory Comput 2012; 8:3338-47. [DOI: 10.1021/ct3004436] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jianming Chen
- Key Laboratory of Organic Solids,
Beijing National Laboratory for Molecular Science (BNLMS), Institute
of Chemistry, Chinese Academy of Sciences, 100190 Beijing, P. R. China
| | - Dong Wang
- MOE Key Laboratory of Organic
OptoElectronics and Molecular Engineering, Department of Chemistry,
Tsinghua University, 100084 Beijing, P. R. China
| | - Zhigang Shuai
- Key Laboratory of Organic Solids,
Beijing National Laboratory for Molecular Science (BNLMS), Institute
of Chemistry, Chinese Academy of Sciences, 100190 Beijing, P. R. China
- MOE Key Laboratory of Organic
OptoElectronics and Molecular Engineering, Department of Chemistry,
Tsinghua University, 100084 Beijing, P. R. China
| |
Collapse
|
28
|
Thermal conductivity of silicon and carbon hybrid monolayers: a molecular dynamics study. J Mol Model 2012; 18:4811-8. [DOI: 10.1007/s00894-012-1482-4] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2012] [Accepted: 05/21/2012] [Indexed: 11/25/2022]
|
29
|
English NJ, Gorman PD, MacElroy JMD. Mechanisms for thermal conduction in hydrogen hydrate. J Chem Phys 2012; 136:044501. [DOI: 10.1063/1.3677189] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Wan L, Liang D, Wu N, Guan J. Molecular dynamics simulations of the mechanisms of thermal conduction in methane hydrates. Sci China Chem 2012. [DOI: 10.1007/s11426-011-4473-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
31
|
Wang D, Shi W, Chen J, Xi J, Shuai Z. Modeling thermoelectric transport in organic materials. Phys Chem Chem Phys 2012; 14:16505-20. [DOI: 10.1039/c2cp42710a] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
32
|
Xu L, Wang X, Liu L, Yang M. First-principles investigation on the structural stability of methane and ethane clathrate hydrates. COMPUT THEOR CHEM 2011. [DOI: 10.1016/j.comptc.2011.09.035] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
33
|
Lin YS, Hsiao PY, Chieng CC. Constructing a force interaction model for thermal conductivity computation using molecular dynamics simulation: ethylene glycol as an example. J Chem Phys 2011; 134:154509. [PMID: 21513397 DOI: 10.1063/1.3578184] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This study aims to construct a force interaction model for thermal conductivity computation and to analyze the liquid properties in atomic level for liquid ethylene glycol (EG) using molecular dynamic simulation. The microscopic details of the molecular system and the macroscopic properties of experimental interest are connected by Green-Kubo relations. In addition, the major contributions of heat transfer modes for thermal conductivity due to convection, interaction, and torque are obtained quantitatively. This study reveals that the intramolecular interaction force fields result in different conformations of the EG in the liquid and thus the molecular shapes. The trans∕gauche ratio for EG's O-Me-Me-O torsional angle and the number of intermolecular∕intramolecular H-bonds are found to be important parameters affecting the thermal conductivity.
Collapse
Affiliation(s)
- Yung-Sheng Lin
- Department of Engineering and System Science, National Tsing Hua University, Hsinchu, Taiwan
| | | | | |
Collapse
|
34
|
|
35
|
Kuang S, Gezelter JD. A gentler approach to RNEMD: Nonisotropic velocity scaling for computing thermal conductivity and shear viscosity. J Chem Phys 2010; 133:164101. [DOI: 10.1063/1.3499947] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Conde MM, Vega C. Determining the three-phase coexistence line in methane hydrates using computer simulations. J Chem Phys 2010; 133:064507. [DOI: 10.1063/1.3466751] [Citation(s) in RCA: 167] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Guo Z, Zhang D, Zhai Y, Gong XG. The intriguing thermal conductivity of ice nanotubes. NANOTECHNOLOGY 2010; 21:285706. [PMID: 20585161 DOI: 10.1088/0957-4484/21/28/285706] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
We have investigated the thermal conductivity of various ice nanotubes (Ice-NTs) using the nonequilibrium molecular dynamics method. The results indicate that Ice-NTs have an unusually high thermal conductivity compared to that of the bulk ices. The thermal conductivity is sensitive to temperature, tube length and diameter, while being insensitive to polarization. We have also studied the confinement effect from single-walled carbon nanotubes (SWCNs). A very remarkable increase in the thermal conductivity is further observed after the Ice-NTs are confined in SWCNs.
Collapse
Affiliation(s)
- Zhixin Guo
- Key Laboratory for Computational Physics (MOE), and Surface Physics Laboratory (National Key), Fudan University, Shanghai 200433, People's Republic of China
| | | | | | | |
Collapse
|
38
|
|
39
|
Jacobson LC, Molinero V. A Methane−Water Model for Coarse-Grained Simulations of Solutions and Clathrate Hydrates. J Phys Chem B 2010; 114:7302-11. [DOI: 10.1021/jp1013576] [Citation(s) in RCA: 126] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Liam C. Jacobson
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850
| | - Valeria Molinero
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850
| |
Collapse
|
40
|
Conde MM, Vega C, McBride C, Noya EG, Ramírez R, Sesé LM. Can gas hydrate structures be described using classical simulations? J Chem Phys 2010; 132:114503. [DOI: 10.1063/1.3353953] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
41
|
Matanović I, Xu M, Moskowitz JW, Eckert J, Bačić Z. Methane molecule confined in the small and large cages of structure I clathrate hydrate: Quantum six-dimensional calculations of the coupled translation-rotation eigenstates. J Chem Phys 2009; 131:224308. [DOI: 10.1063/1.3268623] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
42
|
Jacobson LC, Hujo W, Molinero V. Thermodynamic Stability and Growth of Guest-Free Clathrate Hydrates: A Low-Density Crystal Phase of Water. J Phys Chem B 2009; 113:10298-307. [DOI: 10.1021/jp903439a] [Citation(s) in RCA: 230] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Liam C. Jacobson
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850
| | - Waldemar Hujo
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850
| | - Valeria Molinero
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850
| |
Collapse
|
43
|
English NJ, Tse JS. Mechanisms for thermal conduction in methane hydrate. PHYSICAL REVIEW LETTERS 2009; 103:015901. [PMID: 19659158 DOI: 10.1103/physrevlett.103.015901] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/25/2009] [Indexed: 05/28/2023]
Abstract
Crystalline clathrate hydrates exhibit an unusual thermal transport with glasslike thermal conductivity close to the Debye temperature but a crystal-like temperature dependence at low temperature. Molecular dynamics calculations on structure I methane clathrate hydrate reproduced the qualitative trend in the thermal conductivity. Analysis of the heat flux and local energy correlation functions shows that both the crystal structure of the clathrate framework and guest-host interactions contribute to thermal transport processes. The lower thermal conductivity relative to ice Ih is due to differences in crystal structures. The glasslike temperature dependence is governed by the guests and the guest-host interactions.
Collapse
Affiliation(s)
- Niall J English
- The SEC Strategic Research Cluster and the Centre for Synthesis and Chemical Biology, School of Chemical and Bioprocess Engineering, Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, Dublin 4, Ireland.
| | | |
Collapse
|
44
|
Schofield DP, Jordan KD. Molecular Dynamics Simulations of Bromine Clathrate Hydrates. J Phys Chem A 2009; 113:7431-8. [DOI: 10.1021/jp900237j] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Daniel P. Schofield
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| | - Kenneth D. Jordan
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| |
Collapse
|
45
|
Myshakin EM, Jiang H, Warzinski RP, Jordan KD. Molecular Dynamics Simulations of Methane Hydrate Decomposition. J Phys Chem A 2009; 113:1913-21. [DOI: 10.1021/jp807208z] [Citation(s) in RCA: 115] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Evgeniy M. Myshakin
- National Energy Technology Laboratory, U.S. Department of Energy, P.O. Box 10940, Pittsburgh, Pennsylvania 15236, Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, and Parsons, South Park, Pennsylvania 15129
| | - Hao Jiang
- National Energy Technology Laboratory, U.S. Department of Energy, P.O. Box 10940, Pittsburgh, Pennsylvania 15236, Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, and Parsons, South Park, Pennsylvania 15129
| | - Robert P. Warzinski
- National Energy Technology Laboratory, U.S. Department of Energy, P.O. Box 10940, Pittsburgh, Pennsylvania 15236, Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, and Parsons, South Park, Pennsylvania 15129
| | - Kenneth D. Jordan
- National Energy Technology Laboratory, U.S. Department of Energy, P.O. Box 10940, Pittsburgh, Pennsylvania 15236, Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, and Parsons, South Park, Pennsylvania 15129
| |
Collapse
|