Dissociation mechanism of gas hydrates (I,II, H) of alkane molecules: a comparative molecular dynamics simulation

Employing NPT molecular dynamics method with consistent valence force field, the dissociation processes of sI, sII and sH gas hydrates are simulated at different temperatures and at a constant pressure of 100 MPa. The dissociation mechanisms of gas hydrates are revealed by analysing the structural snapshots, radial distribution functions and diffusion coefficients at different temperatures. As temperature increases, the diffusion rates of water molecules and guest molecules increase; thus the clathrate skeleton formed by water molecules with hydrogen bonds distorts and breaks down; meanwhile the guest molecules encapsulated in the water cavities are released. The size of guest molecules affects the dissociation behaviour of gas hydrate. In addition, the dissociation behaviour also relies on the structural phase of gas hydrates.     

Ab initio DFT study of structural and mechanical properties of methane and carbon dioxide hydrates

The structural and mechanical properties of methane and carbon dioxide hydrates were investigated using density functional theory simulations. Well-established equations of state of solids and exchange-correlation functionals were used for fitting the unit lattice total energy as a function of volume, and the full second-order elastic constants of these two gas hydrates were determined by energy–strain analyses. The polycrystalline elastic properties were also calculated from the unit lattice results. The final results for methane hydrate agree well with available experimental data and with other theoretical results. The two gas hydrates were found to be highly elastically isotropic, but they differed significantly in shear properties. The presented results for carbon dioxide hydrates are the first complete set reported so far. The results are a significant contribution to the ab initio material characterisation of gas hydrates required for ongoing fundamental studies and technological applications.     

Quantized Water Clusters Around Apolar Molecules

Abstract

In order to understand the mechanism of gas hydrate kinetics and to explore the existance of other new cavities in the hydrate structure, we have used Molecular Dynamics (MD) simulation to study a system comprising two Lennard-Jones particles and 214 water molecules. Equilibrium structure and properties of twelve cases have been investigated. Our findings were as follows: ? Apolar molecules promote spherical liquid water clusters in a hydrate-like labile cavity.

? The size of the cavity and the coordination number is dependent upon the size of the apolar molecule.

? The coordination number of water molecules is quantized in jumps of four.

? Similarities are observed between the labile cavities and cavities in solid hydrates and in other chemical structures such as Buckminsterfullerene.

? Such a simulation procedure suggests the possibility of other clusters which may exist in yet-to-be-found hydrates. A separate question involves whether such suggested cavities can be combined with other cavities into a space-filling crystal.

     

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.

Molecular simulations as a tool for selection of kinetic hydrate inhibitors

Natural gas hydrates are ice-like structures in which water molecules form a cage around gas molecules. They have been a problem in the petroleum industry. The heavy cost of alcohol and glycol injections needed to suppress the formation of hydrates has spurred an interest in so-called “kinetic inhibitors”, able to slow down the hydrate formation rather than prevent it. An earlier work (Kvamme, B. et al. 1997, Mol. Phys., 90, p. 979) proposed a simulation-based scheme to assess the comparative performance of prospective inhibitors and select the best candidates for experimental testing. In this work, we employed molecular dynamics simulations to test several kinetic inhibitors in a multiphase water–hydrate system with rigid hydrate interface. In addition, a long-scale run was implemented for a system where the hydrate was free to melt and reform. Our conclusion that PVCap inhibitor will outperform PVP as a kinetic hydrate inhibitor is supported by experimental data. We demonstrate that numerical experiments can be a valuable tool for selecting kinetic inhibitors as well as provide insight into mechanisms of kinetic inhibition and hydrate melting and reformation.     

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.     

First-principles study of the structural,elastic and electronic properties of SbXI (X=S,Se, Te) crystals

The structural, elastic, elastic anisotropy and electronic properties of ferroelectric SbSI and paraelectric SbSI, SbSeI and SbTeI crystals were computed using the local density approximation with first-principle calculations, based on density functional theory. The independent elastic constants of SbXI compounds were computed and the results reveal that they are mechanically stable. Some polycrystalline quantities such as bulk modulus, shear modulus, acoustic velocities, Young’s modulus, Poisson’s ratio, elastic anisotropy and elastic Debye temperatures of these compounds were derived from computed elastic constants. Energy band structures show that these compounds have an indirect band gap. The electronic charge distribution and partial density of states of SbXI compounds indicate that the Sb-X bond is typically covalent with a strong hybridization as well as Sb-I compounds that have strong ionic character. The results obtained were compared with experimentally measured values and other theoretical data.     

Monte Carlo Studies on Water in the dCpG/Proflavin Crystal Hydrate

Abstract

The extensive water network identified in the crystallographic studies of the dCpG/Proflavin hydrate by Neidle, Berman and Shieh (Nature 288, 129, 1980) forms an ideal test case for a) assessing the accuracy of theoretical calculations on nucleic acid—water systems based on statistical thermodynamic computer simulation, and b) the possible use of computer simulation in predicting the water positions in crystal hydrates for use in the further refinement and interpretation of diffraction data. Monte Carlo studies have been carried out on water molecules in the unit cell of dCpG/proflavin, with the nucleic acid complex fixed and the condensed phase environment of the system treated by means of periodic boundary conditions. Intermolecular interactions are described by potential functions representative of quantum mechanical calculations developed by Clementi and coworkers, and widely used in recent studies of the aqueous hydration of various forms of DNA fragments. The results are analyzed in terms of hydrogen bond topology, hydrogen bond distances and energies, mean water positions, and water crystal probability density maps. Detailed comparison of calculated and experimentally observed results are given, and the sensitivity of results to choice of potential is determined by comparison with simulation results based on a set of empirical potentials.     

The Stability of Clathrate Hydrates: Temperature Dependence of Dissociation Pressure in Xe and Ar Hydrate

Abstract

The thermodynamic stability of clathrate hydrates I and II encaging xenon or argon has been investigated by examining the temperature dependence of the dissociation pressure. The evaluation of the stability is made based on the generalized van der Waals and Platteeuw theory developed by Tanaka and Kiyohara [J. Chem. Phys. 98, 4098 (1993)]. In the new treatment, the free energy of formation of hydrates in equilibrium with ice is calculated by taking the coupling of the host lattice vibrations with guests into consideration. The predicted dissociation pressures of Xe and Ar hydrates agree well with experiments in higher temperature range. A poor agreement between experiment and calculation for Ar clathrate hydrate at low temperature is improved by the use of a quantum mechanical partition function for a harmonic oscillator in evaluating the free energy difference between ice and empty hydrate.     

Van der Waals interactions involving proteins.

Van der Waals (dispersion) forces contribute to interactions of proteins with other molecules or with surfaces, but because of the structural complexity of protein molecules, the magnitude of these effects is usually estimated based on idealized models of the molecular geometry, e.g., spheres or spheroids. The calculations reported here seek to account for both the geometric irregularity of protein molecules and the material properties of the interacting media. Whereas the latter are found to fall in the generally accepted range, the molecular shape is shown to cause the magnitudes of the interactions to differ significantly from those calculated using idealized models, with important consequences. First, the roughness of the molecular surface leads to much lower average interaction energies for both protein-protein and protein-surface cases relative to calculations in which the protein molecule is approximated as a sphere. These results indicate that a form of steric stabilization may be an important effect in protein solutions. Underlying this behavior is appreciable orientational dependence, one reflection of which is that molecules of complementary shape are found to exhibit very strong attractive dispersion interactions. Although this has been widely discussed previously in the context of molecular recognition processes, the broader implications of these phenomena may also be important at larger molecular separations, e.g., in the dynamics of aggregation, precipitation, and crystal growth.     

Activities and distribution of methanogenic and methane‐oxidizing microbes in marine sediments from the Cascadia Margin

We investigated methane production and oxidation and the depth distribution and phylogenetic affiliation of a functional gene for methanogenesis, methyl coenzyme M reductase subunit A (mcrA), at two sites of the Integrated Ocean Drilling Program Expedition 311. These sites, U1327 and U1329, are respectively inside and outside the area of gas hydrate distribution on the Cascadia Margin. Radiotracer experiments using ¹⁴C‐labelled substrates indicated high potential methane production rates in hydrate‐bearing sediments [128–223 m below seafloor (mbsf)] at U1327 and in sediments between 70 and 140 mbsf at U1329. Tracer‐free experiments indicated high cumulative methane production in sediments within and below the gas hydrate layer at U1327 and in sediments below 70 mbsf at U1329. Stable tracer experiments using ¹³C‐labelled methane showed high potential methane oxidation rates in near‐surface sediments and in sediments deeper than 100 mbsf at both sites. Results of polymerase chain reaction amplification of mcrA in DNA were mostly consistent with methane production: relatively strong mcrA amplification was detected in the gas hydrate‐bearing sediments at U1327, whereas at U1329, it was mainly detected in sediments from around the bottom‐simulating reflector (126 mbsf). Phylogenetic analysis of mcrA separated it into four phylotype clusters: two clusters of methanogens, Methanosarcinales and Methanobacteriales, and two clusters of anaerobic methanotrophic archaea, ANME‐I and ANME‐II groups, supporting the activity measurement results. These results reveal that in situ methanogenesis in deep sediments probably contributes to gas hydrate formation and are inconsistent with the geochemical model that microbial methane currently being generated in shallow sediments migrates downward and contributes to the hydrate formation. At Site U1327, gas hydrates occurred in turbidite sediments, which were absent at Site U1329, suggesting that a geological setting suitable for a gas hydrate reservoir is more important for the accumulation of gas hydrate than microbiological properties.     

Pressure effects in confined nanophases

In this article, we review how pressure effects in pores affect both the physics of the confined fluid and the properties of the host porous material. Molecular simulations in which high-pressure effects were observed are first discussed; we will see how the strong dependence on bulk phase pressure of the freezing temperature of a fluid confined in nanopores can be explained by important variations of the pressure within the pore. We then discuss recent works in which direct calculations of the pressure tensor of fluids confined in pores provide evidence for large pressure enhancements. Finally, practical applications of these pressure effects in which gas adsorption in microporous solids (pore size < 2 nm) was found to enhance their mechanical properties by increasing the elastic modulus by a factor 4 are discussed.     

Molecular origins of osmotic second virial coefficients of proteins.

The thermodynamic properties of protein solutions are determined by the molecular interactions involving both solvent and solute molecules. A quantitative understanding of the relationship would facilitate more systematic procedures for manipulating the properties in a process environment. In this work the molecular basis for the osmotic second virial coefficient, B22, is studied; osmotic effects are critical in membrane transport, and the value of B22 has also been shown to correlate with protein crystallization behavior. The calculations here account for steric, electrostatic, and short-range interactions, with the structural and functional anisotropy of the protein molecules explicitly accounted for. The orientational dependence of the protein interactions is seen to have a pronounced effect on the calculations; in particular, the relatively few protein-protein configurations in which the apposing surfaces display geometric complementarity contribute disproportionately strongly to B22. The importance of electrostatic interactions is also amplified in these high-complementarity configurations. The significance of molecular recognition in determining B22 can explain the correlation with crystallization behavior, and it suggests that alteration of local molecular geometry can help in manipulating protein solution behavior. The results also have implications for the role of protein interactions in biological self-organization.     

Modeling and simulation of the mechanical response from nanoindentation test of DNA-filled viral capsids

Viruses can be described as biological objects composed mainly of two parts: a stiff protein shell called a capsid, and a core inside the capsid containing the nucleic acid and liquid. In many double-stranded DNA bacterial viruses (aka phage), the volume ratio between the liquid and the encapsidated DNA is approximately 1:1. Due to the dominant DNA hydration force, water strongly mediates the interaction between the packaged DNA strands. Therefore, water that hydrates the DNA plays an important role in nanoindentation experiments of DNA-filled viral capsids. Nanoindentation measurements allow us to gain further insight into the nature of the hydration and electrostatic interactions between the DNA strands. With this motivation, a continuum-based numerical model for simulating the nanoindentation response of DNA-filled viral capsids is proposed here. The viral capsid is modeled as large- strain isotropic hyper-elastic material, whereas porous elasticity is adopted to capture the mechanical response of the filled viral capsid. The voids inside the viral capsid are assumed to be filled with liquid, which is modeled as a homogenous incompressible fluid. The motion of a fluid flowing through the porous medium upon capsid indentation is modeled using Darcy’s law, describing the flow of fluid through a porous medium. The nanoindentation response is simulated using three-dimensional finite element analysis and the simulations are performed using the finite element code Abaqus. Force-indentation curves for empty, partially and completely DNA-filled capsids are directly compared to the experimental data for bacteriophage λ. Material parameters such as Young’s modulus, shear modulus, and bulk modulus are determined by comparing computed force-indentation curves to the data from the atomic force microscopy (AFM) experiments. Predictions are made for pressure distribution inside the capsid, as well as the fluid volume ratio variation during the indentation test.     

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.     

Ultraviolet-induced thymine hydrates in DNA are excised by bacterial and human DNA glycosylase activities

Ultraviolet irradiation of DNA results in various pyrimidine modifications. We studied the excision of an ultraviolet thymine photoproduct by Escherichia coli endonuclease III and by a preparation of human WI-38 cells. These enzymes cleave UV-irradiated DNA at apyrimidinic sites formed by glycosylic removal of the photoproduct. Poly(dA-[3H]dT).poly(dA-[3H]dT) was UV irradiated and incubated with purified E. coli endonuclease III. 3H-Containing material was released in a manner consistent with Michaelis-Menten kinetics. This 3H-labeled material was determined to be a mixture of thymine hydrates (6-hydroxy-5,6-dihydrothymine), separable from unmodified thymine by chromatography in three independent systems. Both cis-thymine hydrate and trans-thymine hydrate were chemically and photochemically synthesized. These coeluted with the enzyme-released 3H-containing material. No thymine glycol was released from the UV-irradiated polymer. Similar results were obtained with extracts of WI-38 cells as the enzyme source. The release of thymine hydrates by both glycosylase activities was directly proportional to the amount of enzyme and the irradiation dose to the DNA substrate. These results demonstrate the modified thymine residues recognized and excised by endonuclease III and the human enzyme to be a mixture of cis-thymine hydrate and trans-thymine hydrate. The reparability of these thymine hydrates suggests that they are stable in DNA and therefore potentially genotoxic.     

Ab initio and DFT study of the aromaticity of some Fulvalenes derived from Methylidenecyclopropabenzene

Methylidencyclopropabenzene (MCPB) 1 and Fulvalenes 2–4 are molecules of special interest due to the relation between structure and aromaticity. The aim of this work was to analyze this relation and to quantify the aromaticity in 1–4 using different methods. Magnetic properties are directly related with aromaticity; here we studied the magnetic susceptibility and the anisotropy of the magnetic susceptibility. Nucleus indepedent chemical shift (NICS) and the anisotropy of the induced current density (ACID) were also employed. Tools of very different nature, geometric indexes HOMA and Bird, were determinated too for 1–4. All of these measures were found to be in agreement. Figure Both spatial NICS and ACID plot allow to show the aromaticity/antiaromaticity of a ring