Three-phase equilibria of CO$_2$ hydrate from computer simulation in presence of NaCl

Abstract: In this work, the cryoscopic decrease effect, as a function of the NaCl concentration, on the carbon dioxide (CO$_2$) hydrate dissociation line conditions has been determined through molecular dynamic simulations. In particular, we have determined the three-phase (solid hydrate-aqueous phase-liquid CO$_2$) coexistence temperature at 100, 400, and 1000 at several initial NaCl concentrations in the aqueous phase, from 0.0 to $3.0\,\text{m}$, using the direct-coexistence technique. We have used the well-known TIP4P/2005 and TraPPE force fields for water and CO$_2$ molecules, respectively. Also, the water-salt interactions are described using the Madrid-2019 force field, which has been specifically developed for various salts in combination with the TIP4P/2005 water model. According to the results obtained in this work, the dissociation temperature of the CO$_2$ hydrate decreases when the NaCl concentration in the initial aqueous phase increases. The results obtained are in excellent agreement with experimental data reported in the literature. We have also observed how the dynamic of melting and growth of the CO$_2$ hydrate becomes slower when the NaCl concentration is increased. As a consequence, larger simulation times (in the order of dozens of microseconds) are necessary when the NaCl concentration increases. Finally, we have also analyzed finite-size effects on the three-phase coexistence temperature of these systems by performing simulations, at 400 bar, with two different system sizes at two different NaCl concentrations (0.0 and $3.0\,\text{m}$). Non-negligible deviations have been found between the results obtained from both system sizes.