Enthalpies of Hydrate Formation from Hydrate Formers Dissolved in Water

: The international interest in the energy potential related to the huge amounts of methane trapped in the form of hydrates is rapidly increasing. Unlike conventional hydrocarbon sources these natural gas hydrate deposits are widely spread around the world. This includes countries which have limited or no conventional hydrocarbon sources, like for instance Japan. A variety of possible production methods have been proposed during the latest four decades. The pressure reduction method has been dominant in terms of research efforts and associated investments in large scale pilot test studies. Common to any feasible method for producing methane from hydrates is the need for transfer of heat. In the pressure reduction method necessary heat is normally expected to be supplied from the surrounding formation. It still remain, however, unveriﬁed whether the capacity, and heat transport capabilities of surrounding formation, will be sufﬁcient to supply enough heat for a commercial production based on reduction in pressure. Adding heat is very costly. Addition of limited heat in critical areas (regions of potential freezing down) might be economically feasible. This requires knowledge about enthalpies of hydrate dissociation under various conditions of temperature and pressure. When hydrate is present in the pores then it is the most stable phase for water. Hydrate can then grow in the concentration range in between liquid controlled solubility concentrations, and the minimum concentration of hydrate in water needed to keep the hydrate stable. Every concentration in that range off concentrations results unique free energy and enthalpy of the formed hydrate. Similarly for hydrate dissociation towards water containing less hydrate former than


Introduction
Industrial problems related to hydrate formation in pipelines and process equipment has historically motivated much of the hydrate research. In many of these systems the heat released during hydrate formation is rapidly distributed through liquid water and solid metal. Related kinetic models have therefore typically focused on the thermodynamic driving force and related mass transport. During the few latest decades the interest in natural gas hydrates as an energy source has increased substantially. Any method utilized for production of methane gas from methane gas hydrates will involve transfer of heat. Most efforts has been on reducing the pressure to outside hydrate stability, which addresses the thermodynamic driving force but the heat still needs to be supplied from the surroundings, or added in some other way. Some very short pilot tests have been conducted in Alaska [1,2]. In more recent years two tests were conducted off the shore of Japan. The first one lasted for six days and was stopped due to problems with production of sand and water, as well as a freezing down due to limited heat supply capacity from the surroundings [3]. The second test was designed and planned for six months of continuous production but closed down after 24 days due to freezing down problems [4] Thermal stimulation, through for instance steam or hot water, is also possible but considered far too expensive as the primary method. A more novel technology is to inject carbon dioxide. A solid state conversion mechanism has been proven by Ripmeester et al. [5] and Kuhs et al. [6] for temperatures far below zero Celsius. For temperatures above zero Celsius a substantially faster mechanism [7,8] involves the formation of a new CO 2 dominated hydrate inside the pores using free pore water. Even in permafrost areas it is hard to find any natural gas hydrate deposits with higher hydrate saturations than 85% of the available pore volume. The more typical is 75% or lower like also in the Ignik Sikimu pilot test [9,10]. Injection of CO 2 will lead to the formation of a new CO 2 hydrate from the free water in the pores and the injected CO 2 . Addition of small amounts of surfactant will keep the interface free of hydrate [11] and stimulate to hydrate nucleation and growth below the CO 2 /water interface region. The effects of surfactant on interface free energy and increased solubility of CO 2 below the surfactant stimulated interface are additional effects that will stimulate hydrate growth from solution. Practically it has been observed that small amounts of methanol will act as promotor for hydrate formation in a pipeline [12].
In view of the above there is a need for consistent and reliable data for enthalpies of hydrate formation or dissociation. As will be discussed in Section 2 available experimental data are very limited and normally incomplete in various ways like for instance missing hydrations number, and unprecise information on thermodynamic conditions. Pressure temperature hydrate equilibrium curves can be measured, or calculated using one of the many codes around the world. With additional data on volumetric changes it is possible to calculate enthalpies of phase transitions. This approach, the so-called Clapeyron method, becomes more complex for mixtures. And the method is not directly usable for hydrate formation from solution of hydrate formers in water. The use of residual thermodynamics, as proposed here, is simple and has no limitations as to which phases water and hydrate formers come from.
A brief review of available experimental data for methane hydrate and carbon dioxide hydrate is given in Section 2. The Clapeyron of using the pressure/temperature gradient, and the simplified Claussius-Clapeyron, have been typical methods for estimating enthalpies of hydrate formation. These methods are briefly discussed in Section 2. The main focus of Section 2, however, is the description of the proposed method for calculating enthalpies of hydrate phase transitions. The method is illustrated and compared to available experimental data. Additional results are presented and discussed in Section 3, followed by discussions in Section 4 and our conclusions in Section 5.

Methodologies
The use of residual thermodynamics also for water as liquid, ice or in empty clathrates of structures I and II makes it feasible to calculate thermodynamic properties related to hydrate phase transitions in a consistent way since the reference state for all components, in all phases is the same. The use of model molecules and Molecular Dynamics (MD) simulations will directly provide residual properties through samplings in configurational space. Average ideal gas properties are also available through samplings in momentum space. For rigid molecules like those utilized in our work ideal gas properties are trivially available from atomic masses and moments of inertia. In this work we utilize this approach and evaluate calculated results for various properties towards experimental data from open literature. Enthalpies of hydrate formation along pressure temperature equilibrium curves will of course have same values as dissociation enthalpies for the same conditions, but with opposite sign. Extensions to conditions outside equilibrium follow standard thermodynamic relationships and are the subjects of follow-up papers. The first section of methodologies describes the theory for situations in which temperature and pressure are the independent thermodynamic variables. This is the most commonly studied route to hydrate formation from water and a separate hydrate former phase.
For hydrate formation from dissolved hydrate formers the main modifications of the calculation scheme are for the change in hydrate former enthalpy, from a rather dilute solution in water to being captured inside a hydrate cavity. Unlike hydrate formers from a gas (or liquid) the enthalpy changes for hydrate formers in this case is smaller. This is expected since water have to optimize entropy in hydrate bonding structures around these solutes. This results in very cavity like water structures that we, and many other research groups, have reported during the latest three decades. In terms of enthalpies of hydrate formation for homogeneous hydrate formation from water solution we accordingly expect smaller absolute values than for the heterogeneous case.

Residual Thermodynamic Modeling of Hydrate Formation from Water and Separate Hydrate Former Phase
By definition a hydrate phase transition is reversible along the equilibrium curve. This is also utilized for the Clapeyron based methods for the formation of hydrate from a separate gas (or liquid) hydrate former phase and a free water phase. The free energy change for this phase transition can be written as: The superscript H 1 is used to distinguish the specific heterogeneous phase transition from other hydrate formation phase transitions. T is temperature, P is pressure. x is the mole-fraction in either liquid or hydrate (denoted with a superscript H) while y is the mole-fraction in gas (or liquid) hydrate former phase. i is an index for hydrate formers. The superscript water denotes a water phase that is converted into hydrate. Generally this is ice or liquid but in this work we only consider liquid water. µ is the chemical potential. In subsequent papers this can also be chemical potential for water adsorbed in various distances from a solid mineral surface.
Liquid water chemical potential is calculated from the symmetric excess conventions as: The approximation on the right-hand side is strictly not necessary. We could utilize a theoretical model for the activity coefficient or we could also use Gibbs-Duhem [13] to relate this to the activity of CO 2 in water as we have done in our Phase Field Theory modelling of CO 2 hydrate phase transition dynamics studies [14,15]. Since the focus here is to illustrate the complexity of multiple hydrate formation in systems of water and CO 2 we use a simpler kinetic model which is more visible in terms of the various contributions to the hydrate phase transition dynamics. As such the approximation on the right-hand side of (2) is accurate enough for the purpose. The solubility of CH 4 in water is small and right hand side will be close to pure water chemical potential.
Chemical potential for water in the hydrate structure is given by [16]: in which H denotes hydrate and 0 in the superscript on first term on right hand side means empty clathrate. These chemical potentials are readily available from model water (TIP4P) simulations [17]. The number of cavities per water ν k is 1/23 for small cavities of structure I and 3/23 for large cavities.
With CO 2 as only guest i is 1 in the sum over canonical partition functions for small and large cavities.
where β is the inverse of the universal gas constant times temperature. At equilibrium chemical potential of guest molecules i in hydrate cavity k is equal to the chemical potential of molecules i in the co-existing phase it comes from. For non-equilibrium the chemical potential is adjusted for distance from equilibrium through a Taylor expansion as discussed later. The free energies of inclusion (latter term in the exponent) are reported elsewhere [11,[18][19][20][21]. At thermodynamic equilibrium between a free hydrate former phase and the hydrate, then µ ki is the chemical potential of the guest molecule in the hydrate former phase (gas, liquid or fluid) at hydrate equilibrium temperature and pressure. The corresponding filling fracrions and mole-fractions of methane in the hydrate is given by: θ ki is the filling fraction of component i in cavity type k: where ν is the fraction of cavity per water for the actual cavity type as indicated by subscripts. Corresponding mole-fraction water is then given by: x H and the associated hydrate free energy is then: in which µ denote chemical potential. Subscripts H 2 O and i denote water and hydrate formers respectively. Superscripts H, water and gas denote hydrate, liquid water and gas phases respectively.
x is mole-fraction in liquid or hydrate (superscript H) and y is mole-fraction in hydrate former phase. T and P are temperature and pressure respectively and G is free energy. The ∆ symbol the change in free energy and superscript H 1 indicate this hydrate formation route. The chemical potential for guest molecule j (in the case of this work either CO 2 or CH 4 ) which enters Equations (4) and (8) at equilibrium is, according to residual thermodynamics: where y i is mole-fraction of component i in the gas mixture. φ i is the fugacity coefficient for i. Ideal gas chemical potential for pure i can be trivially calculated for any model molecule via statistical mechanics from mass and intramolecular structure (bond lengths and bond angles). For rigid model molecules this is particularly simple and analytic though molecular weight, centre of mass and the moments of inertia. Together with density and temperature the ideal gas chemical potential is available from the momentum space canonical partition function. We have utilized the SRK [22] equation of state for calculating the fugacity coefficient, and also the density needed for the ideal gas chemical potential calculations. The enthalpy change is trivially related to the corresponding free energy change by the thermodynamic relationship: The superscript "Total" is introduced to also include the penalty of pushing aside the old phases. Practically the total free energy change will be the extensive version of Equation (1) plus the work needed to push away old phases in order to give space for the new hydrate phase. This penalty during the nucleation process of hydrate is interface free energy times contact area between water and hydrate forming phase. See Kvamme et al. [11] for more details. Since critical nuclei sizes are small [23][24][25] the whole particle can be considered as covered with water due to capillary forces. As such it is expected that interface free energy between hydrate and liquid water can be used with confidence. Above critical core size the penalty diminishes rapidly relative to the free energy change benefits as given by Equation (1).
For the liquid water phase in (1), as well as for the empty hydrate chemical potential on right hand side of (3) results are trivially obtained from Kvamme and Tanaka [16]. The second term on right hand side is reorganized as: and the derivatives of the cavity partition functions can be written as: The partial derivatives in the last term on right hand side is numerically differentiated from the polyonomial fits of Kvamme and Tanaka [16]: H H For liquid water the enthalpy is even more trivially obtained by numerical differentiation of the polyonomial fit of chemical potential as function of T given by Kvamme and Tanaka [16].
In an equilibrium situation the chemical potential of same guest in the two cavity types must be the same and these have to be equal to the chemical potential of the same molecule in the phase that it came from. For the heterogeneous case this means chemical potential of the molecule in gas (or liquid) hydrate former phase, but outside of equilibrium the gradients in chemical potentials as function of T, P and mole-fractions have to reflect how the molecule behaves in the cavity.
Enthalpies for various guest molecules in the two types of cavities can be evaluated by Monte Carlo simulations along the lines described by Kvamme and Lund [26] and Kvamme and Førrisdahl [27] by sampling guest water interaction energies and efficient volumes from the movements of the guest molecules, i.e., where U is energy and superscript R denote residual (interaction) contribution. z ki is compressibility factor for the guest molecule i in cavity k. Consistent ideal gas values for the same interaction models that were applied in calculation of the residual values is trivial: In which k B is Boltzmann's constant and V ki is the excluded volume of a molecule of type i in cavity of type k. This latter volume is calculated from the sampled volume of center of mass movements plus the excluded volume due to water/guest occupation. Slightly more complex sampling and calculation for molecules which are not monoatomic (or approximated as monoatomic like methane) but still fairly standard [26,27] and explicit discussion on this is not needed here.
For a relevant temperature span in the order of 10 K (273 K-283 K) the differences in enthalpies as calculated from (16) using Monte Carlo sampled data does not vary substantially and can even be approximated as constant for the purpose of this work. This is as expected since the hydrate water lattice is fairly rigid and the average movements almost the same for the limited temperature range. Sampled cavity partition functions will of course vary significantly over the same temperature range due to the direct exponential (Boltzmann factor) dependency. The interaction models for CH 4 and CO 2 utilized are the same as those utilized by Kvamme and Tanaka [16]. Also note that while there is an average attraction also for CO 2 in small cavity, but the sampled Langmuir constant is very small and not significant. This is also confirmed by the Molecular Dynamics studies along the lines of Kvamme and Tanaka [16], for which the movements of CO 2 in the small cavity interferes with several water libration frequencies and the resulting free energy of inclusion is not beneficial for CO 2 in the small cavity. While small cavity occupation of CO 2 has been found at extreme conditions in the ice range of temperatures in some studies [28] it remains unclear if there would be any significant small cavity filling of CO 2 at all for temperatures above zero degrees Celsius. The sampled values are given in Table 1 below. The energies for CO 2 in small cavity was sampled after a structure I containing only CO 2 in large cavities was stabilized. Then small cavities were gradually filled and simulations were run until the average fluctuations in sampled interaction energies were symmetrical and on average less than 0.5 per cent of the average energy for the hydrate crystal. Practically this energy value does not have any implication on the enthalpy since the canonical partition function for CO 2 in small cavity is practically zero and as such also the filling fraction of CO 2 in small cavity is almost zero. At pressures below 95 bars along the equilibrium curve it is zero to third digit in mole-fraction while the maximum filling fraction in small cavity contributes with 0.006 to the mole-fraction at (290.00K, 403.0 bars). In contrast the calculated mole-fractions of methane in structure I varied between 0.134 at (276.16 K, 25.2 bars) along the equilibrium curve to 0.138 at (290.00 K, 164.7 bars).
The derivative of the chemical potential of a guest molecule i in cavity type k with respect to temperature as needed in Equation (26) is the negative of partial molar entropy for the same guest molecule and can be calculated according to: Equation (15) can then be rearranged into: The sampled temperature derivatives of guest inclusion free energies for methane in small cavity exhibits a very different behaviour due to the slightly higher density (roghly 85% higher density as compared to CH 4 in large cavity). CO 2 in small cavity is not relevant since the canonocal partition functions are practically zero within significance, as discussed above. The samples values for the gradient in the last term is plotted in Figure 1 below. The energies for CO2 in small cavity was sampled after a structure I containing only CO2 in large cavities was stabilized. Then small cavities were gradually filled and simulations were run until the average fluctuations in sampled interaction energies were symmetrical and on average less than 0.5 per cent of the average energy for the hydrate crystal. Practically this energy value does not have any implication on the enthalpy since the canonical partition function for CO2 in small cavity is practically zero and as such also the filling fraction of CO2 in small cavity is almost zero. At pressures below 95 bars along the equilibrium curve it is zero to third digit in mole-fraction while the maximum filling fraction in small cavity contributes with 0.006 to the mole-fraction at (290.00K, 403.0 bars). In contrast the calculated mole-fractions of methane in structure I varied between 0.134 at (276.16 K, 25.2 bars) along the equilibrium curve to 0.138 at (290.00 K, 164.7 bars).
The derivative of the chemical potential of a guest molecule i in cavity type k with respect to temperature as needed in equation (26) is the negative of partial molar entropy for the same guest molecule and can be calculated according to: Equation (15) can then be rearranged into: The sampled temperature derivatives of guest inclusion free energies for methane in small cavity exhibits a very different behaviour due to the slightly higher density (roghly 85% higher density as compared to CH4 in large cavity). CO2 in small cavity is not relevant since the canonocal partition functions are practically zero within significance, as discussed above. The samples values for the gradient in the last term is plotted in figure 1 below. Residual enthalpies for hydrate former in a separate hydrate former phase are trivially given by: In which the same equation of state (SRK) is utilized as the one used for calculating fugacity coefficients for the chemical potentials. As examples we plot the calculated equilibrium curves for CH 4 and CO 2 in Figure 2a and the corresponding calculated heats of hydrate formation along the equilibrium curves in Figure 2b. The calculated results for CH 4 in 1a) is compared to experimental data in Kvamme et al. [11] and comparison between calculated results and experimental data for CO 2 equilibrium curve can be found in Kvamme [10]. What is frequently missing in reported experimental data on enthalpies of hydrate dissociation is the hydrate number. Calculated mole-fraction hydrate former along the equilibrium curves for CH 4 and CO 2 is plotted in Figure 3a and the corresponding hydrate number (mole water per mole hydrate former) is plotted in Figure 3b. As expected the trend for CH 4 is smooth and development as expected in terms of temperature and pressure changes along the equilibrium curve. For CO 2 there are distinct and expected changes after the change to high density.

by:
In which the same equation of state (SRK) is utilized as the one used for calculating fugacity coefficients for the chemical potentials. As examples we plot the calculated equilibrium curves for CH4 and CO2 in Figure 2a and the corresponding calculated heats of hydrate formation along the equilibrium curves in Figure 2b. The calculated results for CH4 in 1a) is compared to experimental data in Kvamme et al. [11] and comparison between calculated results and experimental data for CO2 equilibrium curve can be found in Kvamme [10]. What is frequently missing in reported experimental data on enthalpies of hydrate dissociation is the hydrate number. Calculated mole-fraction hydrate former along the equilibrium curves for CH4 and CO2 is plotted in Figure 3a and the corresponding hydrate number (mole water per mole hydrate former) is plotted in Figure  3b. As expected the trend for CH4 is smooth and development as expected in terms of temperature and pressure changes along the equilibrium curve. For CO2 there are distinct and expected changes after the change to high density.
In which the same equation of state (SRK) is utilized as the one used for calculating fugacity coefficients for the chemical potentials. As examples we plot the calculated equilibrium curves for CH4 and CO2 in Figure 2a and the corresponding calculated heats of hydrate formation along the equilibrium curves in Figure 2b. The calculated results for CH4 in 1a) is compared to experimental data in Kvamme et al. [11] and comparison between calculated results and experimental data for CO2 equilibrium curve can be found in Kvamme [10]. What is frequently missing in reported experimental data on enthalpies of hydrate dissociation is the hydrate number. Calculated mole-fraction hydrate former along the equilibrium curves for CH4 and CO2 is plotted in Figure 3a and the corresponding hydrate number (mole water per mole hydrate former) is plotted in Figure  3b. As expected the trend for CH4 is smooth and development as expected in terms of temperature and pressure changes along the equilibrium curve. For CO2 there are distinct and expected changes after the change to high density. Outside of equilibrium the only phase that needs special attention is the hydrate phase since description of all fluids phase are continuous while the theory for hydrate is based on a Langmuir type of equilibrium adsorption theory: in which the equilibrium properties, as reference state, is available from the calculations leading to Figure 2a. The third term on the right hand side is a trivial volume term and the final temperature term on the right hand side can be avoided by expanding in a constant temperature from the equilibrium curve. The second term involves an iteration in which a new set of filling fractions are calculated from the properties at the new pressure. the correction is very limited and no further iteration is needed. The corresponding chemical potentials needed for Equation (8) in the non-equilibrium situation is then trivially given by: which then in turn can be used for calculation of phase transition enthalpies outside of hydrate equilibrium.

Residual Thermodynamic Modeling of Hydrate Formation from Water and Dissolved Hydrate Former in Water
where the superscript H 2 now refers to homogeneous hydrate formation from liquid water and dissolved hydrate former in water according to a similar notation in Kvamme et al. [11]. The chemical potential for dissolved methane in water can then [11] be formulated as: with: where T R is temperature divided by critical temperature for methane. The maximum temperature used in the fitting is 325 K. Ideal gas as function of temperature and density is trivial to consistently calculate using the TIP4P model moments in inertia for the rotational contribution [16]. The activity coefficient for methane in water, based on the asymmetric excess convention (activity coefficient for CH 4 in water unity as mole-fraction CH 4 goes to zero) has been fitted to the following function: where T R is reduced temperature and defined as actual T in Kelvin divided by critical temperature for methane (190.6 K). The lower summation 1,2 indicates starting from 1 and counting in steps of 2. Parameters are given by Kvamme et al. [11] ( Table 3 in that paper). For CO 2 a slightly different approach is utilized. The density of CO 2 as dissolved in water will correspond to the partial molar volume of CO 2 at infinite dilution. The infinite dilution ideal gas chemical potential is not very sensitive to pressure so the following approximation to only temperature dependency is considered as adequate [29]: where T 0,R is 273.15 K divided by the actual temperature. Equation (11) does not apply to temperatures above 303, due to the limited range of temperatures for which infinite partial molar volumes are used, and for temperatures above 273.15 K. The fugacity coefficient for CO 2 in water is fitted to the following function: where T R is reduced temperature and defined as actual T in Kelvin divided by critical temperature for CO 2 (304.35 K). The lower summation 1, 2 indicates starting from 1 and counting in steps of 2. Parameters are given in Table 2 below. The vector sign on mole-fraction x denote the vector of mole-fractions. The chemical potential for CO 2 which applies to Equations (10) and (4) for an equilibrium case is the given by: Since the chemical potential of CO 2 is not necessarily, the same for dissolved CO 2 in water, and CO 2 in gas (or liquid) in a non-equilibrium situation, then hydrate formed according to Equation (23) will be different from the first hydrate and accordingly denoted H 2 . The composition of this hydrate will be different as seen from the corresponding compositions, which follows from Equations (4) to (8).
For homogeneous hydrate formation from liquid water and dissolved hydrate former then the number of degrees of freedom (Gibbs phase rule) is 2 and each specific set of temperature and pressure will have a maximum dissolved content of hydrate former according a Henry's law type of calculation [11], and a minimum content of hydrate former needed to keep contacting hydrate at stability. Examples for CH 4 are given in Figure 4 below for CH 4 for a larger range of temperatures and pressures than those reported by Kvamme et al. [11] (Figure 16 in that paper). Similar results for CO 2 are given in Figure 5.
will be different from the first hydrate and accordingly denoted H2. The composition of this hydrate will be different as seen from the corresponding compositions, which follows from equations (4) to (8).
For homogeneous hydrate formation from liquid water and dissolved hydrate former then the number of degrees of freedom (Gibbs phase rule) is 2 and each specific set of temperature and pressure will have a maximum dissolved content of hydrate former according a Henry's law type of calculation [11], and a minimum content of hydrate former needed to keep contacting hydrate at stability. Examples for CH4 are given in figure 4 below for CH4 for a larger range of temperatures and pressures than those reported by Kvamme et.al [11] (figure 16 in that paper). Similar results for CO2 are given in Figure 5.  One should note the strange frames that are drawn between the corners of figures 3b) and 4b) which seem to be artificial drawing in Matlab, as some straight lines drawn between the four corners in the fairly symmetric data. This might be due to the arrangement of the data from the calculations,  For homogeneous hydrate formation from liquid water and dissolved hydrate former then the number of degrees of freedom (Gibbs phase rule) is 2 and each specific set of temperature and pressure will have a maximum dissolved content of hydrate former according a Henry's law type of calculation [11], and a minimum content of hydrate former needed to keep contacting hydrate at stability. Examples for CH4 are given in figure 4 below for CH4 for a larger range of temperatures and pressures than those reported by Kvamme et.al [11] (figure 16 in that paper). Similar results for CO2 are given in Figure 5.  One should note the strange frames that are drawn between the corners of figures 3b) and 4b) which seem to be artificial drawing in Matlab, as some straight lines drawn between the four corners in the fairly symmetric data. This might be due to the arrangement of the data from the calculations, One should note the strange frames that are drawn between the corners of Figures 3b and 4b which seem to be artificial drawing in Matlab, as some straight lines drawn between the four corners in the fairly symmetric data. This might be due to the arrangement of the data from the calculations, which are in a long row of 50 temperatures for each of 50 pressures in the range from 50 to 400 bars. This is of course not a physical part of the calculated data and accordingly must be neglected in the reading of the curves.
In Figure 5b there is a fairly steep change in liquid mole-fraction for lower limit of hydrate stability, corresponding to the rapid change in density as also reflected in Figure 2a and also the rest of results illustrated in Figures 2 and 3. Methane is supercritical in these ranges of temperatures and do have any sharp density increase in the same range. The contours in Figure 4b is slightly concave, and again-the outer strange "side-walls" are not physical.
The only change in Equation (19) for this case is in the cavity partition functions so that Equation (4), now with more specific notations, can be written at equilibrium as: where the chemical potentials are either from Equation (24) for CH 4 or (29) for CO 2 . For the case of equilibrium it is assumed that chemical potential for the guest molecule has same chemical potential in both cavities. Outside of equilibrium then Equation (21) applies, and the Taylor expansion will also imply that the chemical potentials of guest in the two types of cavities are not neccesarily the same.
A second set of changes, relative to the case of separate phases for water and hydrate former, lies in the enthalpy for the hydrate former that enters in the changes associated with the phase transition: and:

Results
It is outside the scope of this work to provide a comprehensive review and discussion of available experimental data. That could have been natural and easy if the experimental data available were following certain standards. The problem is that all the experimental data that we have found so far in open literature is lacking some critical information and the associated error bars are uncertain. Frequent reasons for error bars are lack of information on hydrate number, no indication of pressure and lack of information on degree of super heating. For these reasons there are substantial differences in reported values for heats of dissociation between different publications. These differences are reflected in differences between reported results from various research groups even though they utilize the same experimental method. Some experimental methods are indirect in the sense that heat is not measured directly but interpreted indirectly via measurements of the pressure temperature hydrate stability limits. One example of the use of Clausius-Clapeyron is plotted in Figure 6a below. The only experimental information they utilize is the gradient of their measured equilibrium curve. Their experimental data for the pressure temperature hydrate equilibrium curve is very close to recent data, and also in accordance with Figure 2a. They utilized the calculated compressibility factor, z, from the SRK in the following classical version of the Clausius simplification of Clapeyron's equation: Their data essentially follows a straight line as function of temperature, with the gas compressibility factor for CH 4 that roughly change linearly from around 0.9 at 280 K to approximately 0.85 at 282.5 K. Both end point presumably on the hydrate equilibrium pressure so that the corresponding pressures are 52 bar and 137 bars, respectively. Some other available data from open literature is plotted in as well.
Similar comparison for CO 2 is plotted in Figure 6b. Nasir and Lal [29] only plotted data up to 282.5 and such-the extrapolation is not fair above that temperature since they stopped before the density change for CO 2 . Clapeyron would not even be a fair approximation for higher temperature for several reasons; including the fact that SRK [22] is not good as calculating liquid densities (and compressibility factors) for CO 2 .
Their data essentially follows a straight line as function of temperature, with the gas compressibility factor for CH4 that roughly change linearly from around 0.9 at 280 K to approximately 0.85 at 282.5 K. Both end point presumably on the hydrate equilibrium pressure so that the corresponding pressures are 52 bar and 137 bars, respectively. Some other available data from open literature is plotted in as well.
Similar comparison for CO2 is plotted in Figure 6b. Nasir and Lal [29] only plotted data up to 282.5 and such -the extrapolation is not fair above that temperature since they stopped before the density change for CO2. Clapeyron would not even be a fair approximation for higher temperature for several reasons; including the fact that SRK [22] is not good as calculating liquid densities (and compressibility factors) for CO2.

Discussion
The advantage of the proposed method for calculating enthalpies of hydrate formation is that it is totally consistent with the same residual thermodynamic models which are utilized to calculate hydrate formation equilibrium and dynamics. This includes hydrate numbers and also a route for calculating heats of formation outside of the equilibrium curve. The latter may not be as straightforward in the Clapeyron method since that method is utilizing the gradients of the P-T equilibrium curve. Reported hydrate numbers vary substantially between different reported values for heat of hydrate dissociation. Sometimes only a fixed value (independent of T and P) is used and sometimes ideal value is used. Some care have to be taken in comparing values form different group, as seen in Figures 7a and b where comparison of values per mole hydrate is in fair agreement while our values for heat of formation per mole hydrate former is lower in the case of methane hydrate. The old data from Deaton and Frost [31] is hard to evaluate in terms of accuracy based on the original source from 1946.

Discussion
The advantage of the proposed method for calculating enthalpies of hydrate formation is that it is totally consistent with the same residual thermodynamic models which are utilized to calculate hydrate formation equilibrium and dynamics. This includes hydrate numbers and also a route for calculating heats of formation outside of the equilibrium curve. The latter may not be as straightforward in the Clapeyron method since that method is utilizing the gradients of the P-T equilibrium curve. Reported hydrate numbers vary substantially between different reported values for heat of hydrate dissociation. Sometimes only a fixed value (independent of T and P) is used and sometimes ideal value is used. Some care have to be taken in comparing values form different group, as seen in Figure 7a,b where comparison of values per mole hydrate is in fair agreement while our values for heat of formation per mole hydrate former is lower in the case of methane hydrate. The old data from Deaton and Frost [31] is hard to evaluate in terms of accuracy based on the original source from 1946.
In general any concept based on the Clapeyron method, or the simplified Clausius-Clapeyron method, using calculated equilibrium curves, will be affected by how the pressure temperature gradient is calculated. While we use residual thermodynamics also for the empty clathrate, so as to keep these values consistent with chemical potentials for ice and liquid water, most codes today use a fitting approach in which the difference between liquid (or ice) water and empty clathrate is treated as an empirical fitting parameter. In addition to fitting parameters in the interactions between water and hydrate formers filling the cavities the fitting of chemical potential differences often goes via fitting chemical difference between water phase and empty clathrate at a reference temperature. Then a corresponding enthalpy difference is fitted in order to obtain temperature dependency for the chemical potential difference according to a relationship similar to Equation (10), which also then requires a specific heat capacity difference. And also a volume difference is needed for the pressure difference. The sum of it all results in a variety of possibilities for fitted parameters towards data experimental data. In the final end it is highly unlikely that two different codes based on this chemical potential difference route should end up with same calculations of filling fractions of various cavities in hydrates. We have no possibility to examine how large the differences will be. Experimental equilibrium curves may of course also be used, although in many cases then without any information about hydrate numbers.  In general any concept based on the Clapeyron method, or the simplified Clausius-Clapeyron method, using calculated equilibrium curves, will be affected by how the pressure temperature gradient is calculated. While we use residual thermodynamics also for the empty clathrate, so as to keep these values consistent with chemical potentials for ice and liquid water, most codes today use a fitting approach in which the difference between liquid (or ice) water and empty clathrate is treated as an empirical fitting parameter. In addition to fitting parameters in the interactions between water and hydrate formers filling the cavities the fitting of chemical potential differences What is normally not discussed much is the limited stability window of hydrate. Pressure and Temperature, like in Figure 2a is just a projection of a stability window. So even if hydrate may be inside stability range in terms of temperature and pressure it is still unstable if the outside concentration of CH 4 dissolved in water is smaller than the values in Figure 4b. This can be seen clearly in Figure 7 as function of pressure. The solubility of CH 4 in water at 50 bars is very low compared to the two higher pressures. The corresponding window of possible hydrate growth from solution, as given by the distance in concentrations between 4a and 4b, narrows down when the pressures goes down to 50 bars. Similarly also the stability window for hydrate goes down with pressure. This can even be better illustrated in terms of thermodynamic driving forces (free energy change of the hydrate formation). In Figure 8 we therefore plot the values of molar (per mole hydrate) free energy changes according to Equation (23)  What is normally not discussed much is the limited stability window of hydrate. Pressure and Temperature, like in Figure 2a) is just a projection of a stability window. So even if hydrate may be inside stability range in terms of temperature and pressure it is still unstable if the outside concentration of CH4 dissolved in water is smaller than the values in Figure 4b). This can be seen clearly in Figure 7 as function of pressure. The solubility of CH4 in water at 50 bars is very low compared to the two higher pressures. The corresponding window of possible hydrate growth from solution, as given by the distance in concentrations between 4a) and 4b), narrows down when the pressures goes down to 50 bars. Similarly also the stability window for hydrate goes down with pressure. This can even be better illustrated in terms of thermodynamic driving forces (free energy change of the hydrate formation). In Figure 8 we therefore plot the values of molar (per mole hydrate) free energy changes according to Equation (23)  The example calculations plotted in Figure 8 illustrates the limited stability window for hydrates, which is rarely discussed in open literature, which is useful in various hydrate production schemes. But dissociation towards water which is under saturated (with reference to hydrate) will dissociate hydrate. This is the reason that new fractures through hydrate always will lead to new leakage fluxes of methane. And in geological settings formation which are in risk of developing geo mechanical instabilities and slides this route to hydrate dissociation have to be accounted for in addition to effects of temperature changes.
Another advantage is that the proposed method is applicable for all various situations of hydrate formation as long as the chemical potentials of water and hydrate formers are available. This is illustrated in this work through calculations of enthalpies of hydrate formation from aqueous solution. To our knowledge this has never been investigated theoretically before so we were not able to find neither theoretical values nor experimental data to compare with.
There is a need for a thorough analysis of available experimental data from various methods for measuring the heats of hydrate dissociation for hydrates. Both direct and indirect methods needs to be evaluated and also subjected to a sensitivity analysis so as to evaluate the relative impact of uncertainties related to unknown parameters and/or uncertainties in various key parameters. As discussed above these include hydrate number, degrees of super heating (calorimetry) and impact of applied pressure during dissociation of the hydrate. But the list can become longer as the screening proceeds. Work is in progress and expected to be submitted soon.
Yet another advantage lies in the possibility of studying hydrate nucleation towards solid surfaces. Earlier studies of water and hydrate formers [34][35][36][37][38] and references in these indicate that CO 2 may partially adsorb directly on some mineral surfaces, in competition with water, while CH 4 can adsorb in water layers that are structured by the interactions with the solid surface. Work is in progress on quantification of the adsorption thermodynamics so as to be able to investigate in more detail the associated kinetics of hydrate nucleation along the lines of the simplified classical theory [11] as well as modified version of our Phase Field Theory (PFT) model [7,8,[39][40][41]. Some examples of values for adsorption chemical potential of CO 2 on calcite can be found in reference [36]. But several other estimated adsorption chemical potential values are available in open literature for hydrate formers and various mineral surfaces, including data from our research group. Some precautions needs to be taken since various cuttings of model minerals will unique surfaces in terms of exposed atoms and associated atomic charges. Implementation is trivial since this chemical potential enters Equations (1) and (4). Water for hydrate formation can also come from the layers adsorbed on the mineral surface, although the two first water layers of water have too low chemical potential [35] to leave adsorption in favor of hydrate formation.
As we have discussed in earlier publications most situations of hydrates in nature and industry cannot reach thermodynamic equilibrium because there are too many active phases involved in hydrate phase transitions. Hydrate forming mixtures are not uniform in thermodynamics benefits of adsorption on liquid water. Heavier components drop out first and slightly polar components (CO 2 , H 2 S). The combined first and second laws of thermodynamics will therefore also shift mixtures over to creating several hydrates with different compositions, even for heterogeneous hydrate formation from a gas mixture and liquid water. For homogeneous hydrate formation from dissolved hydrate formers in liquid water every concentration of every dissolved hydrate former will result in a unique hydrate. The proposed method does not contain any limitations in the number of hydrate formers that can be included, although the method has been demonstrated for pure components due to available data from open literature to compare the results with.

Conclusions
Available experimental data range for hydrate dissociation range from direct measurements to indirect measurements of the equilibrium curve in pressure, temperature projection and conversion of this information over to the enthalpy change of the phase transition. Approximations and lack of detailed information on all the properties needed induces various levels of bias in the reported data for heats of dissociation. We propose an alternative method using residual thermodynamics. This method is consistent with equilibrium calculations using the same thermodynamic models. An advantage of this method is that the associated entropy changes for the phase transition are consistent with the changes in free energy and enthalpy. Yet another advantage is that the method is not limited to a system of hydrate formation from separate phases of hydrate formers and water. As we have demonstrated it can be applied for any phase transition of hydrate formation and dissociation as long as models for chemical potentials of water and hydrate formers are available. With modern tools for molecular dynamics simulations this is feasible also for studies of hydrate nucleation towards solid surfaces and other substances which may stimulate or affect hydrate nucleation. Although we have demonstrated the method for pure components there are no limits in the numbers of hydrate formers that enters hydrate.