Abstract:The theory developed in our earlier papers is extended to predict dynamical and thermodynamic properties of clathrate structures by accounting for the possibility of multiple filling of cavities by guest molecules. The method is applied to the thermodynamic properties of argon and krypton hydrates, considering both structures I (sI) and II (sII), in which the small cages can be singly occupied and large cages of sII can be singly or doubly occupied. It was confirmed that the structure of the clathrate hydrate … Show more
“…3 The dissociation curve of the Kr-water clathrate systems has been measured by Dyadin et al 16,46 up to 1.5 GPa, and they have observed transitions from sII to sI and to sH structures. Such transitions have been also found from X-ray diffraction studies 14 15 Further, in a recent theoretical study by Subbotin et al 21 where the corresponding transitions occur.…”
Classical thermodynamic properties of a size-specific Kr clathrate-like cluster, modeled by an ab initio and a semiempirical Kr-water potential, were computed from parallel-tempering Monte Carlo simulations. The temperature and pressure dependence of the cluster's heat capacity was studied, and phase diagrams were constructed using a multiple-histogram method. By associating the heat capacity (maxima) and the Pearson correlation coefficient (minima/maxima) values with the phase transitions, attempts were made to identify such changes to particular cluster structures. Various isomers were computed by local optimizations of the interaction enthalpy at a set of randomly selected configurations at each temperature-pressure grid point. Their energy distributions and relative abundances were then employed to assign the observed phase transitions of the cluster. It was found that the structural changes at high pressures are related to Kr clathrate-like cages, such as those of the sI, sII and sH hydrates, as well as a new one at higher pressures, formed by tetragons and pentagons. Such transitions, at low temperatures and as pressure increases, are related to the topology of the intermolecular interactions, that are getting accessible by the sampling in the MC simulations, through the employed volume model.
“…3 The dissociation curve of the Kr-water clathrate systems has been measured by Dyadin et al 16,46 up to 1.5 GPa, and they have observed transitions from sII to sI and to sH structures. Such transitions have been also found from X-ray diffraction studies 14 15 Further, in a recent theoretical study by Subbotin et al 21 where the corresponding transitions occur.…”
Classical thermodynamic properties of a size-specific Kr clathrate-like cluster, modeled by an ab initio and a semiempirical Kr-water potential, were computed from parallel-tempering Monte Carlo simulations. The temperature and pressure dependence of the cluster's heat capacity was studied, and phase diagrams were constructed using a multiple-histogram method. By associating the heat capacity (maxima) and the Pearson correlation coefficient (minima/maxima) values with the phase transitions, attempts were made to identify such changes to particular cluster structures. Various isomers were computed by local optimizations of the interaction enthalpy at a set of randomly selected configurations at each temperature-pressure grid point. Their energy distributions and relative abundances were then employed to assign the observed phase transitions of the cluster. It was found that the structural changes at high pressures are related to Kr clathrate-like cages, such as those of the sI, sII and sH hydrates, as well as a new one at higher pressures, formed by tetragons and pentagons. Such transitions, at low temperatures and as pressure increases, are related to the topology of the intermolecular interactions, that are getting accessible by the sampling in the MC simulations, through the employed volume model.
“…This modification of the SPC/E water-water interaction potential significantly improves the agreement of the unit cell volume calculated in our model for ice I h with the experimental values as described in our previous works. 27,28 The long-range electrostatic interaction has been considered as the interaction between charges on both the hydrogen and oxygen atoms of different water molecules and has been computed by the Ewald method. 35 The charges on hydrogen (q H ¼ þ0.4238jej) and oxygen (q O ¼ À0.8476jej) within the SPCE model were not changed.…”
Section: Theory and Computational Detailsmentioning
confidence: 99%
“…At low methane concentration, the second factor begins to play a more important role because in this case the chemical potential of the guests in the gas phase is significantly increased according to Eq. (28). Thus, the application of methane as a second guest may be preferable especially for stabilization of the CS-II structure due to weaker guest-host interactions in comparison with the commonly used THF molecules.…”
“…It was chosen because it qualitatively describes the kinetics of dissociation [12,23,24]. The modified version was used according to good description of thermodynamics for these phases in our previous studies [43,44]. The short-range interaction was described by Lennard-Jones potential with parameter σ = 3.1556 Å and the energy parameter ε = 0.65063 kJ mol −1 .…”
The object of this study is the kinetic process of solid-liquid first-order phase transitionmelting of carbon dioxide CS-I hydrate with various cavity occupation ratios. The work was done within a framework of study on the local structure of water molecules. These include the time depending change of the short-range order at temperatures close to the melting point and comparison with hexagonal ice structure. Using molecular dynamics method, dependencies of the internal energy of the studied systems on the time of heating were found. Jumps in the internal energy of solids in the range at 275-300 K indicate a phase transition. The study of oxygen-oxygen radial distribution and hydrogen-oxygen-oxygen mutual orientation angles between molecules detached at no more than 3.2 Å allowed to find the H-bond coordination number of all molecules and full number of H-bonds and showed the instant (less than 1 ns) reorganization of short-range order of all molecules. The structure analysis of every neighbor water molecules pairs showed the ~10-15% decrease of H-bond number after the melting whereas all molecules form single long-range hydrogen bond network. The analysis of hydrogen bond network showed the minor changes in the H-bond interaction energy at solid-liquid phase transition.
KeywordsGas hydrate, phase transition, molecular dynamics, hydrogen bond network, short-range order. Recently, the molecular dynamic (MD) studies of gas hydrate described the dissociation processes by
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.