Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation

Abstract

In this work, we determine the dissociation line of the nitrogen (N2) hydrate by computer simulation using the TIP4P/Ice model for water and the TraPPE force field for N2. We use the solubility method proposed recently by some of us to evaluate the dissociation temperature of the hydrate at different pressures, from 500 to 1500 bar. Particularly, we calculate the solubility of N2 in the aqueous solution when it is in contact with a N2-rich liquid phase and when in contact with the hydrate phase via planar interfaces as functions of temperature. Since the solubility of N2 decreases with temperature in the first case and increases with temperature in the second case, both curves intersect at a certain temperature that determines the dissociation temperature at a given pressure. We find a good agreement between the predictions obtained in this work and the experimental data taken from the literature in the range of pressures considered in this work. From our knowledge of the solubility curves of N2 in the aqueous solution, we also determine the driving force for nucleation of the hydrate, as a function of temperature, at different pressures. In particular, we use two different thermodynamic routes to evaluate the change in chemical potential for hydrate formation. Although the driving force for nucleation slightly decreases (in absolute value) when the pressure is increased, our results indicate that the effect of pressure can be considered negligible in the range of pressures studied in this work. To the best of our knowledge, this is the first time the driving force for nucleation of a hydrate that exhibits crystallographic structure sII, along its dissociation line, is studied from computer simulation.