J. Algaba, S. Blazquez, E. Feria, J. M. Míguez, M. M. Conde, F. J. Blas
In this work, the effects of finite size on the determination of the three-phase coexistence temperature ($T_3$) of carbon dioxide (CO$_2$) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved are placed together in the same simulation box. By varying the number of molecules of each phase it is possible to analyze the effect of simulation size and stoichiometry on the $T_3$ determination. In this work, we have determined the $T_3$ value at 8 different pressures and using 6 different simulation boxes with different numbers of molecules and sizes. In 2 of these configurations, the ratio of the number of water and CO$_2$ molecules in the aqueous solution and the liquid CO$_2$ phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO$_2$ in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO$_2$ available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above $T_3$, overestimating the $T_3$ value obtained from direct coexistence simulations. The simulation results obtained for the CO$_{2}$ hydrate confirm the sensitivity of $T_{3}$ depending on the size and composition of the system, explaining the discrepancies observed in the original work by M'iguez emph{et al.} Non-stoichiometric configurations with larger unit cells show convergence of $T_{3}$ values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.
{"title":"Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate","authors":"J. Algaba, S. Blazquez, E. Feria, J. M. Míguez, M. M. Conde, F. J. Blas","doi":"arxiv-2408.02069","DOIUrl":"https://doi.org/arxiv-2408.02069","url":null,"abstract":"In this work, the effects of finite size on the determination of the\u0000three-phase coexistence temperature ($T_3$) of carbon dioxide (CO$_2$) hydrate\u0000have been studied by molecular dynamic simulations and using the direct\u0000coexistence technique. According to this technique, the three phases involved\u0000are placed together in the same simulation box. By varying the number of\u0000molecules of each phase it is possible to analyze the effect of simulation size\u0000and stoichiometry on the $T_3$ determination. In this work, we have determined\u0000the $T_3$ value at 8 different pressures and using 6 different simulation boxes\u0000with different numbers of molecules and sizes. In 2 of these configurations,\u0000the ratio of the number of water and CO$_2$ molecules in the aqueous solution\u0000and the liquid CO$_2$ phase is the same as in the hydrate (stoichiometric\u0000configuration). In both stoichiometric configurations, the formation of a\u0000liquid drop of CO$_2$ in the aqueous phase is observed. This drop, which has a\u0000cylindrical geometry, increases the amount of CO$_2$ available in the aqueous\u0000solution and can in some cases lead to the crystallization of the hydrate at\u0000temperatures above $T_3$, overestimating the $T_3$ value obtained from direct\u0000coexistence simulations. The simulation results obtained for the CO$_{2}$\u0000hydrate confirm the sensitivity of $T_{3}$ depending on the size and\u0000composition of the system, explaining the discrepancies observed in the\u0000original work by M'iguez emph{et al.} Non-stoichiometric configurations with\u0000larger unit cells show convergence of $T_{3}$ values, suggesting that\u0000finite-size effects for these system sizes, regardless of drop formation, can\u0000be safely neglected. The results obtained in this work highlight that the\u0000choice of a correct initial configuration is essential to accurately estimate\u0000the three-phase coexistence temperature of hydrates by direct coexistence\u0000simulations.","PeriodicalId":501146,"journal":{"name":"arXiv - PHYS - Soft Condensed Matter","volume":"1 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141969019","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
J. Algaba, S. Blazquez, J. M. Míguez, M. M. Conde, F. J. Blas
In this work, the effect of the range of the dispersive interactions in the determination of the three-phase coexistence line of the CO$_2$ and CH$_4$ hydrates has been studied. In particular, the temperature ($T_3$) at which solid hydrate, water, and liquid CO$_2$/gas CH$_4$ coexist has been determined through molecular dynamics simulations using different cut-off values (from 0.9 to 1.6 nm) for the dispersive interactions. The $T_3$ of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures $T$. If the hydrate melts, then $T>T_3$. Contrary, if the hydrate grows, then $T