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

Abstract

On the one hand, the excessive release of methane gas after decomposition will accelerate the greenhouse effect; on the other hand, decomposition will change the mechanical property of the hydrate settled layer, and the hydrate at solid state will be decomposed into the liquid water and methane gas, which will reduce the shear strength of the settled layer and weaken the stratum stability, thus inducing severe geological damage [14]. [...]it is of vital significance to investigate the related changes in the mechanical property and microscopic structure at the time of hydrate decomposition for its exploitation and comprehensive utilization. The Natural Gas Hydrate Research Center of Guangzhou Energy Institute, Chinese Academy of Sciences [36,37,38], had discovered that, the decomposition of heating and hydrate inhibitor on hydrate began on the surface, which was gradually transmitted layer by layer from outside in, and the high hole occupancy and low temperature were beneficial for the stabilization of the methane hydrate crystal structure. [...]this paper had utilized the MD method to study the microscopic structure and mechanical property of the defect-containing sI hydrate. 2.

1. Introduction

Methane hydrate [1,2] (hereinafter referred as “hydrate”), which is generally distributed in the abyssal deposit or the permafrost of land area, is the glacial crystalline material formed by methane molecule with water molecule under high-pressure and low-temperature conditions. As one of the future new energies, hydrate is associated with the superiorities of little pollution, great reserves, wide distribution, and high energy density. It is estimated by the related department that the organic carbon reserve in the global hydrate is 1 × 10¹⁶ kg [3], which is 2-times the sum of all coal, oil, and natural gas resources discovered on the earth so far. At present, over 130 kinds of compounds have been identified to form the hydrates with water molecule; typically, the existing hydrates can be classified according to the crystal configuration as sI [4], sII [5], and sH type [6]. The most commonly seen methane hydrate generally exists in sI type in nature, which was chosen to be the object of this study.

Currently, the hydrate technique has been extensively applied in fields such as gas storage and transport, mixed gas separation and abyssal storage of carbon dioxide [7,8,9,10]. Generally, the storage conditions of hydrate will be changed during the exploitation process, while pressure reduction and temperature increase will accelerate the decomposition of hydrate [11,12,13]. On the one hand, the excessive release of methane gas after decomposition will accelerate the greenhouse effect; on the other hand, decomposition will change the mechanical property of the hydrate settled layer, and the hydrate at solid state will be decomposed into the liquid water and methane gas, which will reduce the shear strength of the settled layer and weaken the stratum stability, thus inducing severe geological damage [14]. Consequently, it is of vital significance to investigate the related changes in the mechanical property and microscopic structure at the time of hydrate decomposition for its exploitation and comprehensive utilization. So far, several studies on mechanical properties of the methane hydrate have been done. Winters et al. and Waite et al. [15,16] had investigated the mechanical properties of water-bearing sediments in permafrost and marine areas, and discovered that the various pore fillers would affect the strength and pore pressure of the hydrate. Ruppel et al. [17,18,19] had adopted the tri-axial test to compare and analyze various mechanical properties of the experimentally synthesized and natural hydrates, including the bulk modulus, stress-strain curve, Poisson’s ratio, and shear strength. Kuniyuki et al. [20] discovered that the elastic modulus of the hydrate was sharply decreased with the decomposition of hydrate, which had produced great transient changes. Masayuki et al. [21] had analyzed the relationships of methane hydrate failure degree with temperature, pressure, and hydrate saturation through the tri-axial compression experiment; their research suggested that the shear strength and elastic modulus of the test specimen tended to be increased with the increase in hydrate saturation. Typically, the mechanical property of hydrate plays a decisive role in the mechanical property of the sediment [22].

The experiment to test the hydrate mechanical property is demanding on the used equipment, moreover, the synthesized hydrate in the experiment usually contains lattice defect or cavity and other impurities, which have resulted in great difficulties in investigating the mechanical property of hydrate through experiment. The molecular dynamic (MD) simulation method [23,24,25,26,27,28,29,30] can intuitively reappear the molecular trajectory, which has thus been extensively applied in studying the microscopic substance mechanism [10]. For instance, Rodger et al. [31,32,33,34,35] had adopted MD to investigate the hydrate dynamic inhibitor, and discovered that the inhibitor was amphoteric, among which, the hydrophilic sulfur-oxygen group could suppress hydrate production, while the hydrophobic methyl could promote hydrate formation. Additionally, the supersaturated CH₄ was to the benefit of hydrate synthesis, and the mechanisms of action of various hydrate catalysts were also different. The Natural Gas Hydrate Research Center of Guangzhou Energy Institute, Chinese Academy of Sciences [36,37,38], had discovered that, the decomposition of heating and hydrate inhibitor on hydrate began on the surface, which was gradually transmitted layer by layer from outside in, and the high hole occupancy and low temperature were beneficial for the stabilization of the methane hydrate crystal structure. Wu et al. [39] had discovered based on MD that, the yield limit of hydrate was greatly affected by temperature and grain size.

The actual hydrate system generally contains defects (such as cavity, polycrystal, and vacancy) and impurities (like sand sediment) [40], which will lead to certain differences in the performance and microscopic structure between the actual hydrate system and the complete hydrate system. The application of MD to investigate the mechanical property and microscopic structure evolution during the hydrate decomposition failure process can effectively avoid multiple problems in experimentally investigating the hydrate performance. Consequently, this paper had utilized the MD method to study the microscopic structure and mechanical property of the defect-containing sI hydrate.

2. Results and Discussion 2.1. Failure Process During Stretching

Herein, Figure 1 shows the tensile strain process of the model with the percentage of randomly deleted particles of 5.92%, where X and Y directions represent 2 cell units, and Z direction indicates 10 cell units. At the initial stage of strain, the model still maintained the crystal structure of methane hydrate and it was subject to elastic strain, which could partially increase the model length at the tensile direction. With the increase in strain, crack appeared in the model, the crystal structure near the crack was destroyed, and the entire fault was uniform. By contrast, the crystal structure of methane hydrate remained unchanged in the area far away from the crack. After the appearance of a crack in the model, the strain at the tensile direction in the material gradually decreased.

2.2. Strain–Stress Curves with Various Defect Ratios

In Figure 2, the curve of the complete hydrate structure (denoted as “Perfect”) was consistent with the simulation results by Wu et al. [39]. Also, Young’s modulus of the studied curves agree with the experimental data [10]. The stress–strain curve of each system experienced three stages. The first stage was the linear elastic strain stage, in which the stress showed linearity with the strain. At this stage, the stress–strain curve of each system superposed each other, and it could be figured out based on a rough calculation of Young’s modulus in the experiment that the Young’s modulus of each model was similar. The second stage was the plastic deformation stage, in which the stress–strain curve of each system displayed non-linear changes; when the stress had reached the maximum value, the stress value of each model rapidly decreased within the limited variation range of strain, suggesting that each system had displayed the characteristics of brittle material failure. The third stage was the stage after failure, in which the stress value of each system remained zero with the increase in strain.

It could be discovered when comparing the effects of various proportions of defects on the system mechanical property that there was no obvious rule regarding the effects of various defect proportions on the mechanical property with the reduction in total system particles. When the reduction of system particle number was <5.92%, the maximum system stress was reduced with the decrease in total particles, but it remained greater than the maximum stress of the 2 × 2 × 10 system, which indicated that the deletion of a certain proportion of particles in the system could partially enhance the mechanical property of hydrate. This might be because the microscopic structure energy of the particle-deleting part was optimized during the tensile process, and the motion of some water molecules had formed new hydrogen bonds, which could boost the mechanical property of the entire structure. When the system particle deletion proportion reached 9.02%, the stress-strain curve of the system was markedly reduced compared with that of the complete lattice, and its mechanical property was weakened. This was because excessive particles were deleted in the system, which had resulted in the incomplete structure of the methane hydrate system, along with the weak intermolecular forces on the tensile fracture surface.

2.3. F3 Order Parameters

To intensively explore the relationship between the hydrate system microscopic structure and its mechanical properties, the order parameter was introduced to analyze the hydrate microscopic structure. The order parameters of the water molecule are generally used to depict the arrangement regularity of the water-containing system microscopically. F3 order parameter [41] is the parameter constituted by three water molecule structures, which represents the degree of tetrahedral arrangement formed by any one oxygen atom with other surrounding oxygen atoms, as expressed below (Equation (1)):

F3={^{(cos_φjik|cos_φjik|+cos2(104.52°))2}}

where ϕ_jik is the included angle formed by the specified oxygen atom i with the two oxygen atoms j and k in the adjacent region (assigned as the region that is 3.2 Å from the oxygen atom). The H–O–H included angle in the TIP4P water molecule is 104.52°. The F3 order parameter of normal water is about 0.1, while the hydrate can form the saturated hydrogen bond with the oxygen atom in the ice system as well as the water molecule where the adjacent oxygen atom is located, and its F3 order parameter is 0. The program for computing F3 order parameter was verified with Li et al. [42,43] works.

The F3 order parameters of various defect systems were statistically analyzed, respectively, which focused on discussing the changes in F3 order parameter of the system before and after failure. Firstly, the 2 × 2 × 10 unit methane hydrate with complete lattice structure was analyzed, and the F3 order parameter output near the failure time step during the tensile process was shown in Figure 3.

The x-coordinate represented the size delamination of the system at the tensile direction, while the y-coordinate was the F3 order parameter. The changes in F3 order parameter of the system at the information output frames of 800, 900, 1000, and 2000 were analyzed (the failure occurred between 900–1000 frames). In the figure, the red arrow interval represented the failure location interval of the system microscopic structure in the Z-direction. As could be observed from Figure 3, at the output frames of 800 and 900, the F3 order parameter of the system remained almost unchanged, which fluctuated around 0, suggesting that the system had maintained the hydrate crystal structure. At the output frame of 1000, the F3 order parameter had fluctuated and peaked at 0.06, and it was close to that of normal water (0.1). Moreover, the Z-direction of the system with fluctuation had overlapped the fracture interval. At the output frame of 2000, the F3 order parameter of the system had little change in the Z-direction, which indicated that the part far away from the fracture maintained the hydrate crystal structure after failure. Upon comparison, the F3 order parameters of all the other structures (Supplementary Figure S1) had displayed a similar trend.

3. Materials and Methods 3.1. Simulation Model and Computational Method

The sI hydrate was employed to carry out the mechanical tests and MD simulation in this study. sI methane hydrate is a type of the common hydrate, a unit cell of body-centered cubic (bcc) structure containing 46 water molecules. The hydrate crystal structure used in this paper was modified based on our previous study [42,43], and the defect-free sI hydrate stable structure of the 2 × 2 × 10 (X × Y × Z) unit cell system was constructed as the initial model, as shown in Figure 4. Subsequently, the particles were randomly deleted from the structure, to obtain the hydrate structures with various defect rates. Further details about the stable structure were described in [42,43]. In this paper, the sI hydrate of the complete 2 × 2 × 10 unit cell system was deleted to obtain the hydrate systems with various defects.

In this study, the open-source LAMMPS software (15 Nov 2018 version, Sandia National Laboratories, Albuquerque, NM, USA) [44] is a large-scale atomic/molecular massively parallel simulator, which was adopted for simulation. The Lennard-Jones (LJ) form of hard sphere model [42,45] was utilized to describe the interaction between methane molecules, while the TIP4P model [46] was used to depict the physiochemical properties of water molecules. At the same time, the Lorentz–Berthelot compositing rules [47] were employed in this study to simulate the interaction among various particles in the system.

3.2. Computational Parameters

Parameters such as the total particle number, number of methane molecules, and number of water molecules are shown in Figure 5, where the x-coordinate represents the deletion percentage of system particles, and the y-coordinate indicates the total number of various particles.

The periodic boundary conditions were applied in the X, Y, and Z directions of the simulator, the Nose–Hoover algorithm was used to control the system temperature and pressure, the Velocity–Verlet algorithm [43] was utilized to solve the equation of motion of the particles, and the Ewald algorithm [48] was adopted to calculate the remote coulomb force of H₂O, with the cutoff radius of 10 Å.

The systems after deletion of various percentages of particles were performed stable relaxation; firstly, the isothermal–isochoric (NVT) ensemble [49] (number of particles (N), volume (V), and temperature (T)) was relaxed for 1,000,000 steps at 200 K at the time step of 0.2 fs; subsequently, the isothermal–isobaric (NPT) ensemble [50] (number of particles (N), pressure (P), and temperature (T)) was relaxed under 200 K and 10 MPa conditions for 1,000,000, with the time step of 0.2 fs.

Afterward, the NPT ensemble was adopted for all test results at the constant temperature of 200 K; then, the test specimens were applied tensile strain along the Z-direction until failure at the time step of 0.2 fs and the strain rate during the tension of 1 × 10⁷/s. The system particle information was output at every 1000 steps, and the corresponding data were counted for analysis.

4. Conclusions

The study on molecular dynamics of methane hydrate involves the exploitation and application of methane hydrate technology, which attracted increasing worldwide attention. This study has utilized the MD method to investigate the influence of defect on the microscopic structure and mechanical property of the sI methane hydrate system based on the complete lattice hydrate. Additionally, the F3 order parameters are adopted to examine the microscopic structure of the sI methane hydrate before and after failure. Firstly, the 2 × 2 × 10 unit sI methane hydrate is selected as the initial model to construct the defect-containing hydrate structure. In addition, particle deletion at various ratios is performed on the initial model structure, and then the models with various defect proportions are subject to tensile strain, to analyze the evolution of microscopic structures and discover the mechanical property for each model system to maintain the brittle materials.

Subsequently, the stress–strain relationship of each system is calculated and analyzed, and it is discovered that the mechanical property of the system will be reduced at the system particle deletion rate of 9.02%. The maximum stress value of this system is lower than that of the initial model and will develop brittle failure prior to the initial structure; moreover, the mechanical properties of models at other defect ratios are partially enhanced compared with the initial model. Additionally, the microscopic structures of the sI methane hydrate before and after failure were analyzed by means of the F3 order parameter value method, the F3 order parameters near the crack will be subject to great fluctuations at the time of failure of the sI methane hydrate structure. Since the actual methane hydrate usually contains lattice defects or cavity and other impurities, the phenomenon and conclusions drawn in this study provide a basis for the study of the microscopic structure, mechanical characteristics, and practical exploitation of methane hydrate. Moreover, the structural properties and microscopic mechanism of the defect-containing hydrate will be further studied, and the experimental verification will be conducted in our future work.

Figure 1. Snapshots of 2 × 2 × 10 (deletion percentage of 5.92%) system during stretching. (The arrow is the place where hydrate cracked).

Figure 2. The strain-stress curves in the Z-axis of different deletion percentage of hydrate systems during stretching.

Figure 3. F3 order parameter for 2 × 2 × 10 unit cells of sI methane hydrate during stretching. (The arrows are the place where hydrate cracked).

Figure 4. Simulation model for 2 × 2 × 10 unit cells of sI methane hydrate.

Figure 5. Deletion percentage of system particles for 2 × 2 × 10 unit cells of sI methane hydrate.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/9/2305/s1.

Author Contributions

S.C. and Q.T. participated in the design of the work, methodology, data interpretation, and analysis for the work; carried out the statistical analyses; and drafted the manuscript. S.T. and Y.L. designed the study; participated in data interpretation, analysis for the work, and methodology; carried out the statistical analyses; and drafted the manuscript. X.G. participated in the design of the work, data interpretation, and analysis for the work.

Funding

This research is supported by Chongqing Research Program of Basic Research and Frontier Technology (cstc2018jcyjAX0522), the National Natural Science Foundation of China (51874055, 21878147), the Program for Changjiang Scholars and Innovative Research Team in University (IRT_17R112), and Innovation Support Program for Chongqing Overseas Returnees (cx2018071).

Acknowledgments

The authors would like to acknowledge Qibin Li (Chongqing University) for providing their structures to verify the program of the F3 order parameter, and the colleagues from the State Key Laboratory of Coal Mine Disaster Dynamics and Control for their perspectives and suggestions related to experimental study and data analyses.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

References

1. Sloan, E.D., Jr.; Carolyn, K. Clathrate Hydrates of Natural Gases, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2007.

2. Feng, Y.; Chen, L.; Suzuki, A.; Kogawa, T.; Okajima, J.; Komiya, A.; Maruyama, S. Enhancement of gas production from methane hydrate reservoirs by the combination of hydraulic fracturing and depressurization method. Energy Conv. Manag. 2019, 184, 194-204.

3. Kvenvolden, K.A. Methane hydrate-A major reservoir of carbon in the shallow geosphere? Chem. Geol. 1988, 71, 41-51.

4. McMullan, R.K.; Jeffrey, G.A. Polyhedral Clathrate Hydrates. IX. Structure of Ethylene Oxide Hydrate. J. Chem. Phys. 1965, 42, 2725-2732.

5. Mak, T.C.W.; McMullan, R.K. Polyhedral Clathrate Hydrates. X. Structure of the Double Hydrate of Tetrahydrofuran and Hydrogen Sulfide. J. Chem. Phys. 1965, 42, 2732-2737.

6. Ripmeester, J.A.; Tse, J.S.; Ratcliffe, C.I.; Powell, B.M. A new clathrate hydrate structure. Nature 1987, 325, 135-136.

7. Kvenvolden, K.A. Gas Hydrates-Geological Perspective and Global Change. Rev. Geophys. 1993, 32, 173-187.

8. Englezos, P.; Lee, J. Gas hydrates: A cleaner source of energy and opportunity for innovative technologies. Korean J. Chem. Eng. 2005, 22, 671-681.

9. Linga, P.; Kumar, R.; Englezos, P. Gas hydrate formation from hydrogen/carbon dioxide and nitrogen/carbon dioxide gas mixtures. Chem. Eng. Sci. 2007, 62, 4268-4276.

10. Ning, F.; Yu, Y.; Kjelstrup, S.; Vlugt, T.J.H.; Glavatskiy, K. Mechanical properties of clathrate hydrates: Status and perspectives. Energy Environ. Sci. 2012, 5, 6779-6795.

11. Gladis, A.; Gundersen, M.T.; Neerup, R.; Fosbol, P.L.; Woodley, J.M.; Sloms, N.V. CO2 mass transfer model for carbonic anhydrase-enhanced aqueous MDEA solutions. Chem. Eng. J. 2018, 335, 197-208.

12. Luis, D.P.; García-González, A.; Saint-Martin, H. A Theoretical Study of the Hydration of Methane, from the Aqueous Solution to the sI Hydrate-Liquid Water-Gas Coexistence. Int. J. Mol. Sci. 2016, 17, 378.

13. Kobayashi, N.; Minami, T.; Tani, A.; Nakagoshi, M.; Sugahara, T.; Takeya, K.; Ohgaki, K. Intermolecular Hydrogen Transfer in Isobutane Hydrate. Energies 2012, 5, 1705-1712.

14. Vysniauskas, A.; Bishnoi, P.R. A Kinetic-Study of Methane Hydrate Formation. Chem. Eng. Sci. 1983, 38, 1061-1072.

15. Waite, W.F.; Kneafsey, T.J.; Winters, W.J.; Mason, D.H. Physical property changes in hydrate-bearing sediment due to depressurization and subsequent repressurization. J. Geophys. Res. Sol. Ea. 2008, 113.

16. Winters, W.J.; Waite, W.F.; Mason, D.H.; Gilbert, L.Y.; Pecher, I.A. Methane gas hydrate effect on sediment acoustic and strength properties. J. Pet. Sci. Eng. 2007, 56, 127-135.

17. Yun, T.S.; Francisca, F.M.; Santamarina, J.C.; Ruppel, C. Compressional and shear wave velocities in uncemented sediment containing gas hydrate. Geophys. Res. Lett. 2005, 32.

18. Francisca, F.; Yun, T.S.; Ruppel, C.; Santamarina, J.C. Geophysical and geotechnical properties of near-seafloor sediments in the northern Gulf of Mexico gas hydrate province. Earth Planet. Sci. Lett. 2005, 237, 924-939.

19. Lee, J.Y.; Santamarina, J.C.; Ruppel, C. Mechanical and electromagnetic properties of northern Gulf of Mexico sediments with and without THF hydrates. Mar. Pet. Geol. 2008, 25, 884-895.

20. Miyazaki, K.; Masui, A.; Aoki, K.; Sakamoto, Y.; Yamaguchi, T.; Okubo, S. Strain-Rate Dependence of Triaxial Compressive Strength of Artificial Methane-Hydrate-Bearing Sediment. Int. J. Offshore Polar Eng. 2010, 20, 9.

21. Hyodo, M.; Wu, Y.; Nakashima, K.; Kajiyama, S.; Nakata, Y. Influence of Fines Content on the Mechanical Behavior of Methane Hydrate-Bearing Sediments. J. Geophys. Res. Sol. Ea. 2017, 122, 7511-7524.

22. Rehder, G.; Kirby, S.H.; Durham, W.B.; Stern, L.A.; Peltzer, E.T.; Pinkston, J.; Brewer, P.G. Dissolution rates of pure methane hydrate and carbon-dioxide hydrate in undersaturated seawater at 1000-m depth. Geochim. Cosmochim. Acta 2004, 68, 285-292.

23. Hu, J.; Liu, C.; Li, Q.; Shi, X. Molecular simulation of thermal energy storage of mixed CO2/IRMOF-1 nanoparticle nanofluid. Int. J. Heat Mass Transf. 2018, 125, 1345-1348.

24. Li, Q.; Liu, C. Molecular dynamics simulation of heat transfer with effects of fluid-lattice interactions. Int. J. Heat Mass Transf. 2012, 55, 8088-8092.

25. Zhan, H.; Gu, Y. A fundamental numerical and theoretical study for the vibrational properties of nanowires. J. Appl. Phys. 2012, 111, 124303.

26. Li, Q.; Liu, C.; Chen, X. Molecular dynamics simulation of sulphur nucleation in S-H2S system. Mol. Phys. 2014, 112, 947-955.

27. Peng, T.; Li, Q.; Liu, C. Accelerated aqueous nano-film rupture and evaporation induced by electric field: A molecular dynamics approach. Int. J. Heat Mass Transf. 2016, 94, 39-48.

28. Li, Q.; Xiao, Y.; Shi, X.; Song, S. Rapid evaporation of water on graphene/graphene-oxide: A molecular dynamics study. Nanomaterials 2017, 7, 265.

29. Hu, J.; Liu, C.; Liu, L.; Li, Q. Thermal energy storage of R1234yf, R1234ze, R134a and R32/MOF-74 nanofluids: A molecular simulation study. Materials 2018, 11, 1164.

30. Li, Q.; Yu, Y.; Liu, Y.; Liu, C.; Lin, L. Thermal properties of the mixed n-octadecane/Cu nanoparticle nanofluids during phase transition: A molecular dynamics study. Materials 2017, 10, 38.

31. Rodger, P.M.; Forester, T.R.; Smith, W. Simulations of the methane hydrate/methane gas interface near hydrate forming conditions. Fluid Phase Equilib. 1996, 116, 326-332.

32. Rodger, P.M. Methane Hydrate: Melting and Memory. Ann. N. Y. Acad. Sci. 2000, 912, 474-482.

33. Moon, C.; Taylor, P.C.; Rodger, P.M. Molecular dynamics study of gas hydrate formation. J. Am. Chem. Soc. 2003, 125, 4706-4707.

34. Storr, M.T.; Taylor, P.C.; Monfort, J.; Roger, P.M. Kinetic Inhibitor of Hydrate Crystallization. J. Am. Chem. Soc. 2004, 126, 1569-1576.

35. Moon, C.; Hawtin, R.W.; Rodger, P.M. Nucleation and control of clathrate hydrates: Insights from simulation. Faraday Discuss. 2007, 136, 367-382.

36. Yan, K.; Li, X.; Chen, Z.; Li, G.; Tang, G.; Fan, S. Molecular dynamics simulation of methane hydrate dissociation by thermal stimulation. Acta Phys. Sin. 2007, 56, 4994-5002.

37. Yan, K.; Li, X.; Chen, Z.; Li, G.; Li, Z. Molecular dynamics simulation of methane hydrate dissociation by thermal stimulation in conjunction with chemical injection method. Acta Phys. Sin. 2007, 56, 6727-6735.

38. Wan, L.; Yan, K.; Li, X.; Fan, S. Molecular Dynamics Simulation of Methane Hydrate Dissociation Process in the Presence of Thermodynamic Inhibitor. Acta Phys. Chim. Sin. 2009, 25, 486-494.

39. Wu, J.; Ning, F.; Trinh, T.T.; Kjelstrup, S.; Vlugt, T.J.H.; He, J.; Skallerud, B.H.; Zhang, Z. Mechanical instability of monocrystalline and polycrystalline methane hydrates. Nat. Commun. 2015, 6, 8743.

40. Wan, L.; Liang, D.; Wu, N.; Guan, J. Molecular dynamics simulations of the mechanisms of thermal conduction in methane hydrates. Sci. China Chem. 2012, 55, 167-174.

41. BÁEz, L.A.; Clancy, P. Computer Simulation of the Crystal Growth and Dissolution of Natural Gas Hydrates. Ann. N. Y. Acad. Sci. 1994, 715, 177-186.

42. Li, Q.; Liu, C.; Chen, X. Molecular characteristics of dissociated water with memory effect from methane hydrates. Int. J. Mod. Phys. B 2014, 28, 1450062.

43. Li, Q.; Tang, Q.; Peng, T.; Zhang, X.; Liu, C.; Shi, X. Molecular characteristics of H2O in hydrate/ice/liquid water mixture. Int. J. Mod. Phys. B 2015, 29, 1550185.

44. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1-19.

45. Delhommelle, J.; MilliÉ, P. Inadequacy of the Lorentz-Berthelot combining rules for accurate predictions of equilibrium properties by molecular simulation. Mol. Phys. 2001, 99, 619-625.

46. Sandberg, T.; Eklund, P.; Hotokka, M. Conformational Solvation Studies of LIGNOLs with Molecular Dynamics and Conductor-Like Screening Model. Int. J. Mol. Sci. 2012, 13, 9845-9863.

47. Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications, 2nd ed.; Academic Press: Cambridge, MA, USA, 2002.

48. Ewald, P.P. Die berechnung optischer und elektrostatischer gitterpotentiale. Ann. Phys. 1921, 369, 253-287.

49. Zhang, L.; Tian, S.; Peng, T. Molecular Simulations of Sputtering Preparation and Transformation of Surface Properties of Au/Cu Alloy Coatings Under Different Incident Energies. Metals 2019, 9, 259.

50. Peng, T.; Li, Q.; Xu, L.; He, C.; Luo, L. Surface interaction of nanoscale water film with SDS from computational simulation and film thermodynamics. Entropy 2017, 19, 620.