Analysis of dissociation process for gas hydrates by molecular dynamics simulation

The dissociation processes of methane and carbon dioxide hydrates were investigated by molecular dynamics simulation. The simulations were performed with 368 water molecules and 64 gas molecules using NPT ensembles. The TraPPE (single-site) and 5-site models were adopted for methane molecules. The EPM2 (3-site) and SPC/E models were used for carbon dioxide and water molecules, respectively. The simulations were carried out at 270 K and 5.0 MPa for hydrate stabilisation. Then, temperature was increased up to 370 K. The temperature increasing rates were 0.1–20 TK/s. The gas hydrates dissociated during increasing temperature or at 370 K. The potential models of methane molecule did not much influence the dissociation process of methane hydrate. The mechanisms of dissociation process were analysed with the coordination numbers and mean square displacements. It was found that the water cages break down first, then the gas molecules escape from the water cages. The methane hydrate was more stable than the carbon dioxide hydrate at the calculated conditions.     

Microhydration of caesium compounds: Cs,CsOH, CsI and Cs2I2 complexes with one to three H2O molecules of nuclear safety interest

The structure and properties of natural gas hydrates containing hydrocarbons, CO₂, and N₂ molecules were studied by using computational quantum chemistry methods via the density functional theory approach. All host cages involved in I, II, and H types structures where filled with hydrocarbons up to pentanes, CO₂ and N₂ molecules, depending on their size, and the structures of these host–guest systems optimized. Structural properties, vibrational spectra, and density of states were analyzed together with results from atoms-in-a-molecule and natural bond orbitals methods. The inclusion of dispersion terms in the used functional plays a vital role for obtaining reliable information, and thus, B97D functional was shown to be useful for these systems. Results showed remarkable interaction energies, not strongly affected by the type of host cage, with molecules tending to be placed at the center of the cavities when host cages and guest molecules cavities are of similar size, but with molecules approaching hexagonal faces for larger cages. Vibrational properties show remarkable features in certain regions, with shiftings rising from host-guest interactions, and useful patterns in the terahertz region rising from water surface vibrations strongly coupled with guest molecules. Likewise, calculations on crystal systems for the I and H types were carried out using a pseudopotential approach combined with Grimme’s method to take account of dispersion.

The thermal properties of binary structure sI clathrate hydrate from molecular dynamics simulation

In this work, we present temperature dependence of lattice parameter and normalised lattice parameter in the atmospheric pressure and 120 bar and also pressure dependence of unit cell volume and normalised unit cell volume at 150 and 250?K for variety guests with different size, polarity and guest–host hydrogen bonding capability such as trimethylene oxide (TMO), ethylene oxide (EO), formaldehyde (FA), cyclobutane (CB), cyclopropane (CP) and ethane (Et) in the large cages with CH₄ in small cages of sI clathrate hydrates by molecular dynamics simulations. The obtained values of lattice parameters for the guest species are compatible with the experimental values. These clathrate hydrates are simulated with TIP4P/ice four-site water potential. Herein, isobaric thermal expansivity and isothermal compressibility are calculated at a temperature range of 50–250?K and a wide pressure range. These structural properties have been compared for guests which they are isoelectronic and have similar masses but with different size and polarity. We use molecular dynamics simulations to relate microscopic guest properties, like guest–host hydrogen bonding to macroscopic sI clathrate hydrate properties. The temperature dependence of thermodynamic properties such as constant-volume and constant-pressure heat capacity is presented in the atmospheric pressure for these guest species.     

Diffusion of N2, O2, H2S and SO2 in MFI and 4A zeolites by molecular dynamics simulations

This paper presents diffusion data of N₂, O₂, H₂S and SO₂ in MFI and 4A zeolites obtained by molecular dynamics simulations, especially its dependence on temperature and loading. At high loadings and temperatures, the order of self-diffusivity of guests in two zeolites is O₂ > N₂ > H₂S>SO₂. The diffusion behaviour is different in different zeolites at lower loadings, reflecting different influences from straight channels (MFI) and α-cages (4A). Furthermore, with increasing loading, the self-diffusivity of guest molecules decreases in MFI but generally increases in 4A. The centre of mass (COM) probability densities and diffusion trajectories of guests give insight into molecular-level diffusion process. The simulation results reveal that with increasing loading, the diffusion mechanism would change from the inter-pore to intra-pore diffusion in MFI. However, in 4A, the intra-pore diffusion is predominant at low and high loadings, but inter-pore diffusion is more important at moderate loadings.     

Molecular dynamics simulation of diffusion behaviour of gas molecules within oil–paper insulation system

The diffusion behaviour of hydrogen, carbon monoxide, carbon dioxide, methane, acetylene, ethylene and ethane in oil and paper medium was examined using molecular dynamics to reveal the diffusion mechanism of gas molecules in transformer oil–paper insulation system at the microscopic level. These compounds are commonly used in the dissolved gas analysis of power transformers and produced during the ageing process of oil–paper composite insulating material. Two groups of models were constructed using molecular dynamics simulation software to simulate the diffusion behaviour of the aforementioned seven types of small gas molecules in oil and paper. The diffusion coefficients, displacement features, free volume characteristics and interaction energies of the gas molecules were investigated. In particular, the diffusion micro-mechanism of the gas molecules was observed. The differences in diffusion features among the gas molecules were discussed, and the factors influencing the diffusion of the gas molecules were compared. Simulation results indicate that the diffusion coefficients of gas molecules in cellulose is an order of magnitude lower than that in oil, and the diffusion coefficients of these gas molecules in the two types of insulation media have different orders. Free volume of gas molecules is the main factor that influences the diffusion behaviour in oil, whereas intermolecular interaction is the main influencing factor of diffusion behaviour in cellulose.     

Molecular dynamics study on thermal dehydration process of epsomite (MgSO4·7H2O)

Water vapour sorption in salt hydrates is one of the most promising means of compact, low loss and long-term solar heat storage in the built environment. Among all, epsomite (MgSO₄·7H₂O) excels for its high-energy storage density and vast availability. However, in practical applications, the slow kinetics and evident structural changes during hydration and dehydration significantly jeopardise the heat storage/recovery rate. A molecular dynamics (MD) study is carried out to investigate the thermal properties and structural changes in the thermal dehydration process of the epsomite. The MD simulation is carried out at 450 K and a vapour pressure of 20 mbar, in accordance with experimental heat storage conditions. The study identifies the dehydration as multiple stages from the initial quick water loss and collapse of the crystal framework to the adsorption of water molecules, which inhibits complete dehydration. Further, the anisotropic diffusion behaviour supports the important role of the porous matrix structure in the heat and mass transfer process. The enthalpy changes, partial densities, mass diffusion coefficients of water and radial distribution functions are calculated and compared with corresponding experimental data to support the conclusions.     

Inclusion of carvone enantiomers in cyclomaltoheptaose (beta-cyclodextrin): thermal behaviour and H-->D and D-->H exchange

The inclusion compounds of carvone enantiomers in cylcomaltoheptaose (beta-cyclodextrin, betaCD) are studied at defined temperatures above room temperature and in relation to H-->D and D-->H exchanges. Loss of water molecules and release of carvone molecules from the betaCD cavity are caused by increase of temperature above room temperature and are measured by the integrated intensities of the O-H and C-H Raman stretching bands, respectively. In turn, H-->D and D-->H exchanges are monitored by the integrated intensities of the O-H and O-D Raman stretching bands, respectively. All of these processes were followed in real time with a Raman spectrometer equipped with CCD detection. The results indicate that distinct carvone enantiomers lead to the formation of different betaCD inclusion hydrates that have different water content and hydration structures. In particular, the results suggest that SCarv-betaCD has a greater water content, dehydrates strongly for temperatures above room temperature, and exchanges protons faster than the RCarv-betaCD complex.     

Microbial Stabilization and Kinetic Enhancement of Marine Methane Hydrates

Abstract

In clathrate hydrates, a water host lattice encages small guest molecules in cavities. Methane hydrates are the most widespread in-situ clathrate in the permafrost and continental-shelf ocean regions, constituting a significant energy resource, and prompting recent marine-hydrate gas-production trials. Despite exciting engineering advances and a few marine-mimicking laboratory studies of methane-hydrate kinetics and stabilization, from microbial perspectives, little is known about a potential microbial origin of marine hydrates, nor their possible formation kinetics or potential stabilization by microbial sources. Here, for the first time, we show that an exported, extra-cytoplasmic porin – produced by a marine methylotrophic bacterium culture – provides the basis for kinetic enhancement and stabilization of methane hydrates under conditions simulating the seabed environment. We then identify the key protein at play, and we therefore suggest microbe-based stabilization of marine hydrates is evidently a property likely to be found in many marine bacteria. Our research opens the possibility of managing marine-hydrate deposits using microbiological strategies for environmental and societal benefit.     

Structural properties of sH hydrate: a DFT study of anisotropy and equation of state

ABSTRACT

Structure-H (sH) hydrate is one of the canonical gas hydrates with significant potential applications and scarce characterised material properties despite the wide knowledge available on other gas hydrates. In this work we characterise some of the important physical properties of this hydrate at the atomistic level using Density Functional Theory. Two exchange-correlation functionals (revPBE and DRSLL) were used to simulate six sH hydrate systems encapsulating neohexane and different help gas molecules. The important role of dispersion forces is quantified. The density and isothermal bulk modulus of sH hydrate are higher when dispersion interactions are considered. The presence of those interactions imposes a direct relationship between the hydrate density and its bulk modulus, while their absence reveals the bulk modulus dependency on hydrogen bond density. Anisotropy is a distinguishing feature of this hydrate in distinction to nearly isotropic sI and sII hydrates. Structure-H hydrate experiences a compressional anisotropy in which the a-lattice and the c-lattice constants respond differently to applied pressure showing less compressibility along the c-axis. This compressional anisotropy was found dependant on the chemistry of help gas molecules. Taken together, these property characterisation results and analysis are a significant and novel contribution to the material physics of sH hydrates.     

Diffusional behavior and guest conformational analysis of hexadecane-1,16-diol and hexadecane in urea crystal model via molecular dynamics simulation approach

Diffusion at the atomic or molecular level is a source of many physical, chemical, and biological processes taking place in plentiful materials. This work is an endeavor toward investigating the diffusional behavior of two different type of guests, hexadecane-1,16-diol and hexadecane enclathration in urea tunnel architecture, whereby the correlation of the diffusion mechanism with the guest’s structural and conformational properties is explored. To carry out this study, molecular dynamics simulation approach is adopted. It is found that hexadecane-1,16-diol exhibit slower diffusion with an average diffusion coefficient value \( \sim 1.83\times {10}^{-10} \), where hexadecane diffuse more rapidly with an average diffusion coefficient value \( \sim 2.58\times {10}^{-9} \). It is also observed that the structural properties influence the guest’s travel distance and torsion angle distribution of the trans and gauche conformational proportion. Furthermore, the observed high energy barrier accounted for hexadecane-1,16-diol and low energy barrier for hexadecane along urea tunnel systems was analyzed. The comparison of our obtained results are in close agreement with the available experimental measurements, i.e., gauche proportion properties between two different guest molecules correlate well with Raman spectroscopy investigation on α,ω-dihalogenoalkane/urea inclusion compounds. Our calculations also successfully endorse the structure-property relation between the two systems.     

Molecular dynamics simulation on the decomposition of type SII hydrogen hydrate and the performance of tetrahydrofuran as a stabiliser

Molecular dynamics simulation is used to study the decomposition and stability of SII hydrogen and hydrogen/tetrahydrofuran (THF) hydrates at 150 K, 220 K and 100 bar. The modelling of the microscopic decomposition process of hydrogen hydrate indicates that the decomposition of hydrogen hydrate is led by the diffusive behaviour of H₂ molecules. The hydrogen/THF hydrate presents higher stability, by comparing the distributions of the tetrahedral angle of H₂O molecules, radial distribution functions of H₂O molecules and mean square displacements or diffusion coefficients of H₂O and H₂ molecules in hydrogen hydrate with those in hydrogen/THF hydrate. It is also found that the resistance of the diffusion behaviour of H₂O and H₂ molecules can be enhanced by encaging THF molecules in the (5¹²6⁴) cavities. Additionally, the motion of THF molecules is restricted due to its high interaction energy barrier. Accordingly, THF, as a stabiliser, is helpful in increasing the stability of hydrogen hydrate.     

Theoretical study of conformational disorder and selective adsorption of small organic molecules in the flexible metal-organic framework material MIL-53-Fe

Many of the potential applications for metal-organic framework (MOF) materials require molecular level understanding of their adsorption of small organic molecules, which are not readily accessible from experiment. Through high-level van der Waals corrected, hybrid density functional theory calculations, we elucidate the adsorption configurations of several, representative small organic guest molecules in an archetypal flexible MOF material, MIL-53-Fe. The predicted relative energies between low-energy adsorption configurations of 1,4-benzoquinone in MIL-53-Fe are in very good agreement with the thermal transition temperatures observed experimentally and suggest that thermodynamic factors govern the precise arrangements and loading of guests in the MOF host. Experimentally observed conformational disorder of small organic molecules in MIL-53-Fe is explained by predicted multiple low-energy adsorption configurations that are comparable with the thermal energy of the guests, kT. Finally, we show that the previously observed selective adsorption of pyridine and 2,6-lutidine molecules in water by MIL-53-Fe, can be rationalised through a careful analysis of the host–guest and guest–guest interactions and is controlled by thermodynamic factors.     

Structure and dynamics of parallel beta-sheets, hydrophobic core, and loops in Alzheimer's A beta fibrils

We explore the relative contributions of different structural elements to the stability of Abeta fibrils by molecular-dynamics simulations performed over a broad range of temperatures (298 K to 398 K). Our fibril structures are based on solid-state nuclear magnetic resonance experiments of Abeta(1-40) peptides, with sheets of parallel beta-strands connected by loops and stabilized by interior salt bridges. We consider models with different interpeptide interfaces, and different staggering of the N- and C-terminal beta-strands along the fibril axis. Multiple 10-20 ns molecular-dynamics simulations show that fibril segments with 12 peptides are stable at ambient temperature. The different models converge toward an interdigitated side-chain packing, and present water channels solvating the interior D23/K28 salt bridges. At elevated temperatures, we observe the early phases of fibril dissociation as a loss of order in the hydrophilic loops connecting the two beta-strands, and in the solvent-exposed N-terminal beta-sheets. As the most dramatic structural change, we observe collective sliding of the N- and C-terminal beta-sheets on top of each other. The interior C-terminal beta-sheets in the hydrophobic core remain largely intact, indicating that their formation and stability is crucial to the dissociation/elongation and stability of Abeta fibrils.     

Vibrational spectra of deuterated methane and water molecules in structure I clathrate hydrate from ab initio MD simulation

The deuteration of the lattice molecules in clathrate hydrates is a widely used experimental technique to clearly separate the vibrational modes. However, the effect of the deuteration on the vibrational spectra and molecular motions is not fully understood. Since the guest–host coupling may change the vibrational spectra, a detailed analysis of the vibrational spectra of deuterated clathrate hydrate is significant in the understanding of the mechanism of the vibrational shift. In this study, the vibrational spectra of the deuterated methane hydrates were calculated by ab initio molecular dynamics simulation. The intramolecular vibrational frequency of the methane in D₂O lattice and deuterated methane in H₂O lattice was calculated and compared with the pure methane hydrate. The bending, rocking and overtone of the bending mode was also reported. The effect of coupling of the rattling motions of guest and host molecules on the vibrational spectra was revealed.     

Influence of moisture on the crystal forms of niclosamide obtained from acetone and ethyl acetate

The purpose of this study was to elucidate the formation of crystal hydrates of niclosamide and to delineate the effect of relative humidity on the crystal forms obtained from acetone and ethyl acetate. Recrystallization of niclosamide was performed in the presence and absence of moisture. Two hydrates and their corresponding anhydrates were isolated. The hydrates obtained by the process of recrystallization from acetone (Form I) and that obtained from ethyl acetate (Form II) were classified based on differences in their dehydration profile, crystal structure, shape, and morphology. Crystals obtained in the absence of moisture were unstable, and when exposed to the laboratory atmosphere transformed to their corresponding hydrates. Differential scanning calorimetry thermograms indicate that Form I changes to an anhydrate at temperatures below 100°C, while Form II dehydrates in a stepwise manner above, 140°C. This finding was further confirmed by thermogravimetric analysis. Dehydration of Form II was accompanied by a loss of structural integrity, demonstrating that water molecules play an important role in maintaining its crystal structure. Form I, Form II, and the anhydrate of Form II showed no significant moisture sorption over the entire range of relative humidity. Although the anhydrate of Form I did not show any moisture uptake at low humidity, it converted to the monohydrate at elevated relative humidity (>95%). All forms could be interconverted depending on the solvent and humidity conditions.     

NMR diffusion as a novel tool for measuring the association constant between cyclodextrin and guest molecules

In this paper we introduce the use of diffusion measurements by nuclear magnetic resonance (NMR) spectroscopy for determining association constants of weak and very weak interactions between cyclodextrin and guest molecules, as long as both the free and complexed guest molecules are soluble to an extent that allows good sensitivity in the NMR experiment. The experimental setup and data analysis is discussed for three different guest molecules: L-phenylalanine, L-leucine and L-valine, representing different strengths of interaction. The underlying assumptions are discussed and the scope of the method (range of K(a) values, requirements to the guest molecule) are discussed. The method's main advantage is its general applicability independent of chromogenic or electrochemical properties of the guest molecule. Whereas calorimetric methods that exhibit a similar generality, are applicable mainly to strong interactions, NMR diffusion measurements are applicable to weaker interactions down to the theoretical limit of 1 M(-1), the upper limit for K(a) values to be determined by it is approximately 200. A further advantage of the method is the low amount of sample needed. The method is in principle applicable to any case of molecular recognition between a host and guest molecule leading to weak interactions.     

Temperature dependence of protein-hydration hydrodynamics by molecular dynamics simulations

The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.     

Microstructures,interactions and dynamics properties studies of N-methyldiethanolamine + guanidinium triflate ionic liquid + water tertiary system at the standard temperature

Molecular dynamics simulations with an all-atom force field have been carried out in order to understand the phase equilibrium behaviour of ternary aqueous mixtures containing guanidinium triflate ionic liquid [gua][OTf] and water mixed with N-methyldiethanolamine (MDEA) in different function composition at the standard temperature of 298.15 K. A very good numerical agreement has been obtained for the prediction of the mixture densities. The analysis of structural and dynamic properties showed that the molecular level of ternary mixtures is slightly affected by the presence of MDEA and [gua][OTf] molar fractions. For MDEA–water interactions in [gua][OTf] media, we found that MDEA prefers to be surrounded by water molecules rather than by MDEA molecules even at a high MDEA molar fraction. While for [gua][OTf]–water interaction in MDEA media, as [gua][OTf] molar fraction increases, water molecules replace counterions in the coordination shell of both ions, thus weakening their interaction. On the other hand, for MDEA–[gua][OTf] interactions in water media, we have found that as the molar fraction of [gua][OTf] increases, a sulfonate group from anion appears to have a stronger association by making hydrogen bonding with MDEA molecules. The chemical process using ionic liquids (ILs) as solvents is commonly limited by their high viscosity. Based on their physical properties such as viscosities, these ternary solvents can be applied in natural gas industry, such as removing carbon dioxide using aqueous MDEA and IL at high pressure.