Prediction of the univariant two-phase coexistence line of the tetrahydrofuran hydrate from computer simulation

Abstract

In this work, the univariant two-phase coexistence line of the hydrate of tetrahydrofuran (THF) is determined from 100 to 1000 bar by molecular dynamics simulations. The study is carried out by putting in contact a THF hydrate phase with a stoichiometric aqueous solution phase. Following the direct coexistence technique, the pressure has been fixed, and the coexistence line has been determined by analyzing if the hydrate phase grows or melts at different values of temperature. The model of water used is the well-known TIP4P/Ice model. We have used two different models of THF based on the transferable parameters for phase equilibria-united atom approach (TraPPE-UA), the original (flexible) TraPPe-UA model as well as a rigid and planar version of it. Overall, at high pressures, small differences have been observed in the results obtained by both models. Also, large differences have been observed in the computational efforts required by the simulations performed using both models, being the rigid and planar version much faster than the original one. The effect of the unlike dispersive interactions between the water and THF molecules has been also analyzed at 250 bar using the rigid and planar THF model. In particular, we have modified the Berthelot combining rule by adding a factor (ξO-THF) that modifies the unlike water-THF dispersive interactions and we have analyzed the effect on the dissociation temperature when ξO-THF is modified from 1.0 (original Berthelot combining rule) to 1.4 (modified Berthelot combining rule). We have extended the study using ξO-THF = 1.4 and the rigid THF model to the rest of the pressures considered in this work, finding an excellent agreement with the scarce experimental data taken from the literature.

10 pages, 5 figures

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Solubility of carbon dioxide in water: some useful results for hydrate nucleation
• Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate
• Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate
• Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates
• Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation

Cited by in corpus (3)

• Simulation of the THF hydrate-water interfacial free energy from computer simulation
• Rotationally invariant local bond order parameters for accurate determination of hydrate structures
• Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation. II. Effect of multiple occupancy