Mechanical instability of monocrystalline and polycrystalline methane hydrates

Molecular insights into the interactions between PEG carriers and drug molecules from Celastrus hindsii: a multi-scale simulation study

Efficient drug delivery is crucial for the creation of effective pharmaceutical treatments, and polyethylene glycol (PEG) carriers have been emerged as promising candidates for this purpose due to their bio-compatibility, enhancement of drug solubility, and stability. In this study, we utilized molecular simulations to examine the interactions between PEG carriers and selected drug molecules extracted from Celastrus hindsii: Hindsiilactone A, Hindsiiquinoflavan B, Maytenfolone A, and Celasdin B. The simulations provided detailed insights into the binding affinity, stability, and structural properties of these drug molecules when complexed with PEG carriers. A multi-scale approach combining density functional theory (DFT), extended tight-binding (xTB), and molecular dynamics (MD) simulations was conducted to investigate both unbound and bound states of PEG/drug systems. The results from DFT and xTB calculations revealed that the unbound complex has an unfavorable binding free energy, primarily due to negative contributions of delta solvation free energy and entropy. The MD simulations provided more detailed insights into the interactions between PEG and drug molecules in water solutions. By integrating the findings from the multi-scale simulations, a comprehensive picture of the unbound and bound states of PEG and drug systems were obtained. This information is valuable for understanding the molecular mechanisms governing the binding of drugs in PEG-based delivery platforms, and it contributes to the rational design and optimization of these systems.

HTR+: a novel algorithm for identifying type and polycrystal of gas hydrates

In this work, the hierarchical topology ring (HTR+) algorithm, an extension of the HTR algorithm, was developed for identifying gas hydrate types, cage structures, and grain boundaries (GBs) within polycrystalline structures. Utilizing molecular dynamics trajectories of polycrystalline hydrates, the accuracy of the HTR+ algorithm is validated in identifying sI, sII and sH hydrate types, hydrate grains, and GBs in multi-hydrate polycrystals, as well as clathrate cages at GBs. Additionally, during the hydrate nucleation and growth processes, clathrate cages, hydrate type, hydrate grains and ice structures are accurately recognized. Significantly, this algorithm demonstrates high efficiency, particularly for large hydrate systems. HTR+ algorithm emerges a powerful tool for identifying micro/mesoscopic structures of gas hydrates, enabling an in-depth understanding of the formation mechanisms and properties of gas hydrates.

Insights into the mechanical stability of tetrahydrofuran hydrates from experimental, machine learning, and molecular dynamics perspectives

Natural gas hydrates (NGHs) hold immense potential as a future energy resource and for sustainable applications such as gas capture and storage. Due to the challenging formation conditions, however, their mechanical properties remain poorly understood. Herein, the mechanical characteristics of tetrahydrofuran (THF) hydrates, a proxy for methane hydrates, were investigated at different ice contents, strain rates, and temperatures using uniaxial compressive experiments. The results unveil a distinct behavior in the peak strength of THF hydrates with a varying ice content, strain rate and temperature, exhibiting an increase as the strain rate and temperature decrease, in contrast to the peak strength-strain rate relationship observed in polycrystalline ice. Based on the experimental data, four machine learning (ML) models including extreme gradient boosting (XGboost), multilayer perceptron (MLP), gradient boosting decision tree (GBDT) and decision tree (DT) were developed to predict the peak strength. The XGboost model demonstrates superior predictive performance, emphasizing the significant influence of ice content and temperature on the peak strength of hydrates. Furthermore, molecular dynamics (MD) simulations were employed to gain insights into the dissociation and formation processes of clathrate cages, as well as phase transitions and amorphization occurring at grain boundaries (GBs) involving diverse unconventional clathrate cages, including 51265, 4151062, 4151064, 425861 and 425862, with 425861 and 425862 cages being predominant. This study enhances our understanding of the mechanical properties and deformation mechanisms of hydrates and provides a ML-based predictive framework for estimating the compressive strength of hydrates under diverse coupling conditions. The findings have significant implications for stability assessments of NGHs and the exploitation of NGH resources.

Mechanical properties of amorphous CO2 hydrates: insights from molecular simulations

Understanding physicochemical properties of amorphous gas hydrate systems is of great significance to reveal structural stabilities of polycrystalline gas hydrate systems. Furthermore, amorphous gas hydrates can occur ordinarily in the nucleation events of gas hydrate systems. Herein, the mechanical properties of amorphous carbon dioxide hydrates are examined by means of all-atom classical molecular dynamic simulations. Our molecular simulation results reveal that mechanical strengths of amorphous carbon dioxide hydrates are evidently governed by temperatures, confining pressures, and ratios of water to carbon dioxide molecules. Notably, under compressive loads, amorphous carbon dioxide hydrates firstly exhibit monotonic strain hardening, followed by an interesting distinct phenomenon characterized by a steady flow stress at further large deformation strains. Furthermore, structural evolutions of amorphous carbon dioxide hydrates are analyzed on the basis of the N-Hbond DOP order parameter. These important findings can not only contribute to our understanding of the structural stabilities of amorphous gas hydrate systems, but also help to develop fundamental understandings about grain boundaries of gas hydrate systems.

Atomistic-geometric simulations to investigate the mechanical stability of monocrystalline sI methane hydrates under pressure

Gas hydrate mechanical stability under pressure is critically important in energy supply, global warming, and carbon-neutral technologies. The stability of these polyhedral guest-host crystals under increasing pressure is affected by host cage type and face connectivity as well as guest gas occupancy. The geometry-imposed cage connectivity generates crystal lattices that include inclusion-matrix material composite structures. In this paper, we integrate Density Functional Theory simulations with a polyhedral-inspired composite material model that quantifies stability limits, failure modes, and the impact of the type of cage occupancy. DFT reveals the existence of two failure mechanisms under increasing pressure: (i) a multistep lattice breakdown under total occupancy and under only large cage occupancy and (ii) a single-step breakdown under zero occupancy as well as with only small cage occupancy. The DFT-composite model predicts optimal occupancy pathways to generate strength and critical occupancy pathways to promote decomposition.

Interfacial mechanical properties of tetrahydrofuran hydrate-solid surfaces: Implications for hydrate management

Understanding the interfacial mechanical properties between hydrate and solids is vital to designing and fabricating surfaces for hydrate management. Herein, the role of the surface wettability, the type of solid substrate and temperature on the interfacial adhesion properties of tetrahydrofuran (THF) hydrate and ice were examined by force analysis based shearing measurements and molecular dynamics (MD) simulations. The results showed that the adhesion strength of THF hydrate and ice on silica varies with the compositions of coating, and the adhesion strength of ice is larger than that of THF hydrate for all investigated solid substrates. Particularly, in contrast to a linear relationship between 1 + cosθr and hydrate adhesion on organic silanes/thiols/polymer surfaces, the hydrate adhesion on the coated inorganic glass surfaces is enhanced as a function of 1 + cosθr, in which θr is the receding contact angle. MD simulations uncovered that the adhesion strength of ice on solid substrates is dominated by the quasi-liquid water layer, however, that of hydrate is governed not only by the quasi-liquid layer but also newly formed unconventional clathrate cages. This study provides new insights and perspectives into the hydrate adhesion on solid surfaces, which is of help to develop hydrate-phobic coatings for advanced hydrate management.

"Ladetes"-A novel device to test deformation behaviors of hydrate-bearing sediments

Natural gas hydrate (NGH) exploitation is severely restricted by geotechnical problems. Deformation behaviors of the hydrate-bearing strata (HBS) control the occurrence and evolution of geotechnical problems during extracting natural gas from HBS. In this paper, a novel approach named Ladetes is introduced to evaluate the lateral deformation behaviors of the near-wellbore and fracture-filling regions of the HBS. The pressuremeter test and the flat dilatometer test are designed to simulate the inner boundaries of an NGH-producing well and an artificial stimulation fracture for the first time. The device can realize the in situ hydrate formation prior to the experiment and axial loading application throughout the experiment. Both the strain control mode and the stress control mode can be achieved to estimate the deformation characteristics of HBS under different downhole conditions. Prime experiments proved their adaptability and reliability. The Ladetes provides an effective and alternative way of obtaining geotechnical parameters for HBS.

Role of mechanical deformation in the thermal transport of sI-type methane hydrate

Natural gas hydrates (NGHs) are rising as an unconventional energy resource. The fundamental thermal characteristics of NGHs are of importance for natural gas exploitation from permafrost and oceanic sediments that are geomechanically deformed. Here, utilizing classic molecular dynamics simulations with all-atom (AA) and coarse-grained (CG) models of the methane guest molecule, the effects of mechanical strain on the thermal conductivity of sI-type methane hydrate are for the first time examined. Upon triaxial tension and compression, methane hydrate exhibits strong asymmetry in the stress responses. As the triaxial loads go from compression to tension, a reduction trend in the thermal conductivity is revealed for methane hydrate with both AA and CG models of methane, within a maximum reduction of over 44%. This reduction is because triaxial strain from compression to tension softens the phonon modes. Interestingly, there is a sudden rise in thermal conductivity at critical triaxial strain of 0.06, originating from that, at which, the phonon modes are hardened and the peaks of radial distribution functions are shifted back. This study provides important information on the thermal conductivity of methane hydrate, which is helpful for the practical production of natural gas from geo-deformed NGH-bearing sediments via a heating technique as well as evaluating their stability.

Unraveling Adhesion Strength between Gas Hydrate and Solid Surfaces

Natural gas hydrate is a promising future energy source, but it also poses a huge threat to oil and gas production due to its ability to deposit within and block pipelines. Understanding the atomistic mechanisms of adhesion between the hydrate and solid surfaces and elucidating its underlying key determining factors can shed light on the fundamentals of novel antihydrate materials design. In this study, large-scale molecular simulations are employed to investigate the hydrate adhesion on solid surfaces, especially with focuses on the atomistic structures of intermediate layer and their influences on the adhesion. The results show that the structure of the intermediate layer formed between hydrate and solid surface is a competitive equilibrium of induced growth from both sides, and is regulated by the content of guest molecules. By comparing the fracture behaviors of the hydrate-solid surface system with different intermediate structures, it is found that both the lattice areal density of water structure and the adsorption of guest molecules on the interface together determine the adhesion strength. Based on the analysis of the adhesion strength distribution, we have also revealed the origins of the drastic difference in adhesion among different water structures such as ice and hydrate. Our simulation indicates that ice-adhesion strength is approximately five times that of lowest hydrate adhesion strength. This finding is surprisingly consistent with the available experimental results.

Mechanical Instability of Methane Hydrate-Mineral Interface Systems

Massive methane hydrates occur on sediment matrices in nature. Therefore, sediment-based methane hydrate systems play an essential role in the society and hydrate community, including energy resources, global climate changes, and geohazards. However, a fundamental understanding of mechanical properties of methane hydrate-mineral interface systems is largely limited due to insufficient experimental techniques. Herein, by using large-scale molecular simulations, we show that the mechanical properties of methane hydrate-mineral (silica, kaolinite, and Wyoming-type montmorillonite) interface systems are strongly dictated by the chemical components of sedimentary minerals that determine interfacial microstructures between methane hydrates and minerals. The tensile strengths of hydrate-mineral systems are found to decrease following the order of Wyoming-type montmorillonite- > silica- > kaolinite-based methane hydrate systems, all of which show a brittle failure at the interface between methane hydrates and minerals under tension. In contrast, upon compression, methane hydrates decompose into water and methane molecules, resulting from a large strain-induced mechanical instability. In particular, the failure of Wyoming-type montmorillonite-based methane hydrate systems under compression is characterized by a sudden decrease in the compressive stress at a strain of around 0.23, distinguishing it from those of silica- and kaolinite-based methane hydrate systems under compression. Our findings thus provide a molecular insight into the potential mechanisms of mechanical instability of gas hydrate-bearing sediment systems on Earth.

Molecular Origins of Deformation in Amorphous Methane Hydrates

Water and methane can stay together under low temperature and high pressure in the forms of liquid solutions and crystalline solids. From liquid and gaseous states to crystalline solids or the contrary processes, amorphous methane hydrates can occur in these evolution scenarios. Herein, mechanical properties of amorphous methane hydrates are explored for the first time to bridge the gap between mechanical responses of monocrystalline and polycrystalline methane hydrates. Our results demonstrate that mechanical properties of amorphous methane hydrates are strongly governed by our original proposed order parameter, namely, normalized hydrogen-bond directional order parameter. Followed by this important achievement, a multistep deformation mechanism core is proposed to explain mechanical properties of amorphous methane hydrates. Through an extensive detailed analysis of amorphous methane hydrates, our simulation results not only greatly enlarge our fundamental understanding for mechanical responses of amorphous methane hydrates in geological systems but also offer a fresh perspective in structure-property topics of solid materials in future science and technology.

Mechanical creep instability of nanocrystalline methane hydrates

Mechanical creep behaviors of natural gas hydrates are of importance for understanding the mechanical instability of gas hydrate-bearing sediments on Earth. Limited by the experimental challenges, intrinsic creep mechanisms of nanocrystalline methane hydrates remain largely unknown yet at the molecular scale. Herein, using large-scale molecular dynamics simulations, mechanical creep behaviors of nanocrystalline methane hydrates are investigated. It is revealed that mechanical creep responses are greatly dictated by internal microstructures of crystalline grain size and external conditions of temperature and static stress. Interestingly, a long steady-state creep is observed in nanocrystalline methane hydrates, which can be described by a modified constitutive Bird-Dorn-Mukherjee model. Microstructural analysis shows that deformations of crystalline grains, grain boundary diffusion and grain boundary sliding collectively govern the mechanical creep behaviors of nanocrystalline methane hydrates. Furthermore, structural transformation also appears to be important in their mechanical creep behaviors. This study provides new insights into understanding the mechanical creep scenarios of gas hydrates.

Crustal fingering facilitates free-gas methane migration through the hydrate stability zone

Widespread seafloor methane venting has been reported in many regions of the world oceans in the past decade. Identifying and quantifying where and how much methane is being released into the ocean remains a major challenge and a critical gap in assessing the global carbon budget and predicting future climate [C. Ruppel, J. D. Kessler. Rev. Geophys. 55, 126-168 (2017)]. Methane hydrate ([Formula: see text]) is an ice-like solid that forms from methane-water mixture under elevated-pressure and low-temperature conditions typical of the deep marine settings (>600-m depth), often referred to as the hydrate stability zone (HSZ). Wide-ranging field evidence indicates that methane seepage often coexists with hydrate-bearing sediments within the HSZ, suggesting that hydrate formation may play an important role during the gas-migration process. At a depth that is too shallow for hydrate formation, existing theories suggest that gas migration occurs via capillary invasion and/or initiation and propagation of fractures (Fig. 1). Within the HSZ, however, a theoretical mechanism that addresses the way in which hydrate formation participates in the gas-percolation process is missing. Here, we study, experimentally and computationally, the mechanics of gas percolation under hydrate-forming conditions. We uncover a phenomenon-crustal fingering-and demonstrate how it may control methane-gas migration in ocean sediments within the HSZ.

Unraveling nucleation pathway in methane clathrate formation

Methane clathrates are widespread on the ocean floor of the Earth. A better understanding of methane clathrate formation has important implications for natural-gas exploitation, storage, and transportation. A key step toward understanding clathrate formation is hydrate nucleation, which has been suggested to involve multiple evolution pathways. Herein, a unique nucleation/growth pathway for methane clathrate formation has been identified by analyzing the trajectories of large-scale molecular dynamics (MD) simulations. In particular, ternary water-ring aggregations (TWRAs) have been identified as fundamental structures for characterizing the nucleation pathway. Based on this nucleation pathway, the critical nucleus size and nucleation timescale can be quantitatively determined. Specifically, a methane hydration layer compression/shedding process is observed to be the critical step in (and driving) the nucleation/growth pathway, which is manifested through overlapping/compression of the surrounding hydration layers of the methane molecules, followed by detachment (shedding) of the hydration layer. As such, an effective way to control methane hydrate nucleation is to alter the hydration layer compression/shedding process during the course of nucleation.

Contactless probing of polycrystalline methane hydrate at pore scale suggests weaker tensile properties than thought

Methane hydrate is widely distributed in the pores of marine sediments or permafrost soils, contributing to their mechanical properties. Yet the tensile properties of the hydrate at pore scales remain almost completely unknown, notably the influence of grain size on its own cohesion. Here we grow thin films of the hydrate in glass capillaries. Using a novel, contactless thermal method to apply stress, and video microscopy to observe the strain, we estimate the tensile elastic modulus and strength. Ductile and brittle characteristics are both found, dependent on sample thickness and texture, which are controlled by supercooling with respect to the dissociation temperature and by ageing. Relating the data to the literature suggests the cohesive strength of methane hydrate was so far significantly overestimated.

The authors here report tensile properties of polycrystalline methane hydrate at the micron scale by applying a contactless, thermos-induced stress to a tenuous shell of hydrate grown in a thin glass capillary. The results suggest that the cohesive strength of methane hydrate in marine settings may be an order of magnitude less than currently thought.

Mechanical Response of Nanocrystalline Ice-Contained Methane Hydrates: Key Role of Water Ice

Water ice and gas hydrates can coexist in the permafrost and polar regions on Earth and in the universe. However, the role of ice in the mechanical response of ice-contained methane hydrates is still unclear. Here, we conduct direct million-atom molecular simulations of ice-contained polycrystalline methane hydrates and identify a crossover in the tensile strength and average compressive flow stress due to the presence of ice. The average mechanical shear strengths of hydrate-hydrate bicrystals are about three times as large as those of hydrate-ice bicrystals. The ice content, especially below 70%, shows a significant effect on the mechanical strengths of the polycrystals, which is mainly governed by the proportions of the hydrate-hydrate grain boundaries (HHGBs), the hydrate-ice grain boundaries (HIGBs), and the ice-ice grain boundaries (IIGBs). Quantitative analysis of the microstructure of the water cages in the polycrystals reveals the dissociation and reformation of various water cages due to mechanical deformation. These findings provide molecular insights into the mechanical behavior and microscopic deformation mechanisms of ice-contained methane hydrate systems on Earth and in the universe.

Molecular dynamics simulation of sI methane hydrate under compression and tension

Molecular dynamics (MD) analysis of methane hydrate is important for the application of methane hydrate technology. This study investigated the microstructure changes of sI methane hydrate and the laws of stress–strain evolution under the condition of compression and tension by using MD simulation. This study further explored the mechanical property and stability of sI methane hydrate under different stress states. Results showed that tensile and compressive failures produced an obvious size effect under a certain condition. At low temperature and high pressure, most of the clathrate hydrate maintained a stable structure in the tensile fracture process, during which only a small amount of unstable methane broke the structure, thereby, presenting a free-motion state. The methane hydrate cracked when the system reached the maximum stress in the loading process, in which the maximum compressive stress is larger than the tensile stress under the same experimental condition. This study provides a basis for understanding the microscopic stress characteristics of methane hydrate.

Gas hydrates in sustainable chemistry

This review includes the current state of the art understanding and advances in technical developments about various fields of gas hydrates, which are combined with expert perspectives and analyses.

Molecular Simulation Study on the Microscopic Structure and Mechanical Property of Defect-Containing sI Methane Hydrate

The study of changes in the related mechanical property and microscopic structure of methane hydrate during the decomposition process are of vital significance to its exploitation and comprehensive utilization. This paper had employed the molecular dynamics (MD) method to investigate the influence of defects on the microscopic structure and mechanical property of the sI methane hydrate system, and to discover the mechanical property for the defect-containing hydrate system to maintain its brittle materials. Moreover, the stress-strain curve of each system was analyzed, and it was discovered that the presence of certain defects in the methane hydrate could promote its mechanical property; however, the system mechanical property would be reduced when the defects had reached a certain degree (particle deletion rate of 9.02% in this study). Besides, the microscopic structures of the sI methane hydrate before and after failure were analyzed using the F3 order parameter value method, and it was found that the F3 order parameters near the crack would be subject to great fluctuations at the time of failure of the hydrate structure. The phenomenon and conclusions drawn in this study provide a basis for the study of the microscopic structure and mechanical characteristics of methane hydrate.

Energy and Environmental Analysis of Membrane-Based CH4-CO2 Replacement Processes in Natural Gas Hydrates

Natural gas hydrates are the largest reservoir of natural gas worldwide. This paper proposes and analyzes the CH4-CO2 replacement in the hydrate phase and pure methane collection through the use of membrane-based separation. The investigation uses a 1 L lab reactor, in which the CH4 hydrates are formed in a quartz sand matrix partially saturated with water. CH4 is subsequently dissociated with a CO2 stream supplied within the sediment inside the reactor. An energy and environmental analysis was carried out to prove the sustainability of the process. Results show that the process energy consumption constitutes 4.75% of the energy stored in the recovered methane. The carbon footprint of the CH4-CO2 exchange process is calculated as a balance of the CO2 produced in the process and the CO2 stored in system. Results provide an estimated negative value, equal to 0.004 moles sequestrated, thus proving the environmental benefit of the exchange process.

Understanding the Pathway of Gas Hydrate Formation with Porous Materials for Enhanced Gas Separation

The reason that the stoichiometry of gas to water in artificial gas hydrates formed on porous materials is much higher than that in nature is still ambiguous. Fortunately, based on our experimental thermodynamic and kinetic study on the gas hydrate formation behavior with classic ordered mesoporous carbon CMK-3 and irregular porous activated carbon combined with density functional theory calculations, we discover a microscopic pathway of the gas hydrate formation on porous materials. Two interesting processes including (I) the replacement of water adsorbed on the carbon surface by gas and (II) further replacement of water in the pore by gas accompanied with the gas condensation in the pore and growth of gas hydrate crystals out of the pore were deduced. As a result, a great enhancement of the selectivity and regeneration for gas separation was achieved by controlling the gas hydrate formation behavior accurately.

Molecular-scale thermally activated fractures in methane hydrates: a molecular dynamics study

Cage by cage slow fracture initiation in methane hydrates is consistent with a thermal activation model.

Pressure-Induced Densification of Ice Ih under Triaxial Mechanical Compression: Dissociation versus Retention of Crystallinity for Intermediate States in Atomistic and Coarse-Grained Water Models

Molecular-dynamics (MD) simulation of triaxially pressurized ice I_h up to 30 kbar at 240 K (with sudden mechanical pressurization from its ambient-pressure structure) has been carried out with both the single-particle mW and atomistic TIP4P-Ice water potentials on systems of up to ∼1 million molecules, for times of the order of 100 ns. It was found that the TIP4P-Ice systems adopted a high-density liquid state above ∼7 kbar, while densification of the mW systems retained essentially crystalline order, owing to a failure for the tetrahedral network to break down appreciably from its ice I_h lattice structure. Both are intermediate states adopted along the path toward respective thermodynamically stable states (and with pressure removal show reversion to I_h for mW and to supercooled liquid for TIP4P-Ice), similar to recent ice electro-freezing simulations in "No Man's Land". Densification kinetics showed faster mW-system adaptation.

Grain-Size-Controlled Mechanical Properties of Polycrystalline Monolayer MoS2

Pristine monocrystalline molybdenum disulfide (MoS₂) possesses high mechanical strength comparable to that of stainless steel. Large-area chemical-vapor-deposited monolayer MoS₂ tends to be polycrystalline with intrinsic grain boundaries (GBs). Topological defects and grain size skillfully alter its physical properties in a variety of materials; however, the polycrystallinity and its role played in the mechanical performance of the emerging single-layer MoS₂ remain largely unknown. Here, using large-scale atomistic simulations, GB structures and mechanical characteristics of realistic single-layered polycrystalline MoS₂ of varying grain size prepared by confinement-quenched method are investigated. Depending on misorientation angle, structural energetics of polar-GBs in polycrystals favor diverse dislocation cores, consistent with experimental observations. Polycrystals exhibit grain-size-dependent thermally induced global out-of-plane deformation, although defective GBs in MoS₂ show planar structures that are in contrast to the graphene. Tensile tests show that presence of cohesive GBs pronouncedly deteriorates the in-plane mechanical properties of MoS₂. Both stiffness and strength follow an inverse pseudo Hall-Petch relation to grain size, which is shown to be governed by the weakest link mechanism. Under uniaxial tension, transgranular crack propagates with small deflection, whereas upon biaxial stretching, the crack grows in a kinked manner with large deflection. These findings shed new light in GB-based engineering and control of mechanical properties of MoS₂ crystals toward real-world applications in flexible electronics and nanoelectromechanical systems.

Influence of temperature on methane hydrate formation

During gas hydrate formation process, a phase transition of liquid water exists naturally, implying that temperature has an important influence on hydrate formation. In this study, methane hydrate was formed within the same media. The experimental system was kept at 1.45, 6.49, and 12.91 °C respectively, and then different pressurization modes were applied in steps. We proposed a new indicator, namely the slope of the gas flow rates against time (dν g /dt), to represent the intrinsic driving force for hydrate formation. The driving force was calculated as a fixed value at the different stages of formation, including initial nucleation/growth, secondary nucleation/growth, and decay. The amounts of gas consumed at each stage were also calculated. The results show that the driving force during each stage follows an inverse relation with temperature, whereas the amount of consumed gas is proportional to temperature. This opposite trend indicates that the influences of temperature on the specific formation processes and final amounts of gas contained in hydrate should be considered separately. Our results also suggest that the specific ambient temperature under which hydrate is formed should be taken into consideration, when explaining the formation of different configurations and saturations of gas hydrates in natural reservoirs.

Elasticity and Stability of Clathrate Hydrate: Role of Guest Molecule Motions

Molecular dynamic simulations were performed to determine the elastic constants of carbon dioxide (CO₂) and methane (CH₄) hydrates at one hundred pressure–temperature data points, respectively. The conditions represent marine sediments and permafrost zones where gas hydrates occur. The shear modulus and Young’s modulus of the CO₂ hydrate increase anomalously with increasing temperature, whereas those of the CH₄ hydrate decrease regularly with increase in temperature. We ascribe this anomaly to the kinetic behavior of the linear CO₂ molecule, especially those in the small cages. The cavity space of the cage limits free rotational motion of the CO₂ molecule at low temperature. With increase in temperature, the CO₂ molecule can rotate easily, and enhance the stability and rigidity of the CO₂ hydrate. Our work provides a key database for the elastic properties of gas hydrates, and molecular insights into stability changes of CO₂ hydrate from high temperature of ~5 °C to low decomposition temperature of ~−150 °C.

Mechanical properties of monocrystalline and polycrystalline monolayer black phosphorus

The mechanical properties of monocrystalline and polycrystalline monolayer black phosphorus (MBP) are systematically investigated using classic molecular dynamic simulations. For monocrystalline MBP, it is found that the shear strain rate, sample dimensions, temperature, atomic vacancies and applied statistical ensemble affect the shear behaviour. The wrinkled morphology is closely connected with the direction of the in-plane shear, dimensions of the samples, and applied ensembles. Particularly, small samples subjected to loading/unloading of the shear deformation along the armchair direction demonstrate a clear mechanical hysteresis loop. For polycrystalline MBP, the maximum shear stress as a function of the average grain size follows an inverse pseudo Hall-Petch type relationship under an isothermal-isobaric (NPT) ensemble, whereas under a canonical (NVT) ensemble, the maximum shear stress of polycrystalline MBP exhibits a 'flipped' behaviour. Furthermore, polycrystalline MBP subjected to uniaxial tension also exhibits a strongly grain size-dependent mechanical response, and it can fail by brittle intergranular and transgranular fractures because of its weaker grain boundary structures and the direction-dependent edge energy, respectively. These findings provide useful insight into the mechanical design of BP for nanoelectronic devices.

Inhibition of methane and natural gas hydrate formation by altering the structure of water with amino acids

Natural gas hydrates are solid hydrogen-bonded water crystals containing small molecular gases. The amount of natural gas stored as hydrates in permafrost and ocean sediments is twice that of all other fossil fuels combined. However, hydrate blockages also hinder oil/gas pipeline transportation, and, despite their huge potential as energy sources, our insufficient understanding of hydrates has limited their extraction. Here, we report how the presence of amino acids in water induces changes in its structure and thus interrupts the formation of methane and natural gas hydrates. The perturbation of the structure of water by amino acids and the resulting selective inhibition of hydrate cage formation were observed directly. A strong correlation was found between the inhibition efficiencies of amino acids and their physicochemical properties, which demonstrates the importance of their direct interactions with water and the resulting dissolution environment. The inhibition of methane and natural gas hydrate formation by amino acids has the potential to be highly beneficial in practical applications such as hydrate exploitation, oil/gas transportation, and flow assurance. Further, the interactions between amino acids and water are essential to the equilibria and dynamics of many physical, chemical, biological, and environmental processes.

Nanoscale deicing by molecular dynamics simulation

Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice adhesion strength by an aqueous water layer, and provide atomistic details that support previous experimental studies. Our results contribute quantitative comparison of nanoscale adhesion strength of ice on hydrophobic and hydrophilic surfaces, and supply for the first time theoretical references for understanding the mechanics at the atomistic origins of macroscale ice adhesion.

Hydraulic fracturing in methane-hydrate-bearing sand

Hydraulic fracturing in methane-hydrate-bearing sand can generate laminar fractures and increase the effective permeability.