Dynamic simulation of double-cased perforation in deepwater high temperature and high-pressure oil and gas wells

In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.

Molecular Dynamics Simulation of Shock Response of CL-20 Co-crystals Containing Void Defects

To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.

Unveiling the crystallization mechanism of cadmium selenide via molecular dynamics simulation with machine-learning-based deep potential

Cadmium selenide(CdSe)is an inorganic semiconductor with unique optical and electronic properties that make it useful in various applications,including solar cells,light-emitting diodes,and biofluorescent tagging.In order to synthesize high-quality crystals and subsequently integrate them into devices,it is crucial to understand the atomic scale crystallization mechanism of CdSe.Unfortunately,such studies are still absent in the literature.To overcome this limitation,we employed an enhanced sampling-accelerated active learning approach to construct a deep neural potential with ab initio accuracy for studying the crystallization of CdSe.Our brute-force molecular dynamics simulations revealed that a spherical-like nu-cleus formed spontaneously and stochastically,resulting in a stacking disordered structure where the competition between hexagonal wurtzite and cubic zinc blende polymorphs is temperature-dependent.We found that pure hexagonal crystal can only be obtained approximately above 1430 K,which is 35 K below its melting temperature.Furthermore,we observed that the solidification dynamics of Cd and Se atoms were distinct due to their different diffusion coefficients.The solidification process was initiated by lower mobile Se atoms forming tetrahedral frameworks,followed by Cd atoms occupying these tetra-hedral centers and settling down until the third-shell neighbor of Se atoms sited on their lattice posi-tions.Therefore,the medium-range ordering of Se atoms governs the crystallization process of CdSe.Our findings indicate that understanding the complex dynamical process is the key to comprehending the crystallization mechanism of compounds like CdSe,and can shed lights in the synthesis of high-quality crystals.

Study of the Relationship Between New Ionic Interaction Parameters and Salt Solubility in Electrolyte Solutions Based on Molecular Dynamics Simulation

Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is proposed by using molecular dynamics simulation,and the relationship between ion-to-ion interaction and salt solubility in a simulated seawater water-salt system is investigated.By analyzing the variation of distance and contact time between ions in an electrolyte solution,from both spatial and temporal perspectives,new parameters were proposed to describe the interaction between ions:interaction distance(ID),and interaction time ratio(ITR).The best correlation between characteristic time ratio and solubility was found for a molar ratio of salt-to-water of 10:100 with a correlation coefficient of 0.96.For the same salt,a positive correlation was found between CTR and the molar ratio of salt and water.For type 1-1,type 2-1,type 1-2,and type 2-2 salts,the correlation coefficients between CTR and solubility were 0.93,0.96,0.92,and 0.98 for a salt-to-water molar ratio of 10:100,respectively.The solubility of multiple salts was predicted by simulations and compared with experimental values,yielding an average relative deviation of 12.4%.The new ion-interaction parameters offer significant advantages in describing strongly correlated and strongly hydrated electrolyte solutions.

Unravelling biotoxicity of graphdiyne:Molecular dynamics simulation of the interaction between villin headpiece protein and graphdiyne

Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.

Molecular dynamics simulation of the flow mechanism of shear-thinning fluids in a microchannel

Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.

Temperature effect on nanotwinned Ni under nanoindentation using molecular dynamic simulation

Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.

Temperature-Induced Unfolding Pathway of Staphylococcal Enterotoxin B:Insights from Circular Dichroism and Molecular Dynamics Simulation

In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.

Investigation of Projectile Impact Behaviors of Graphene Aerogel Using Molecular Dynamics Simulations

Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.

Molecular Dynamics Numerical Simulation of Adsorption Characteristics and Exploitation Limits in Shale Oil Microscopic Pore Spaces

Microscopic pore structure in continental shale oil reservoirs is characterized by small pore throats and complex micro-structures.The adsorption behavior of hydrocarbons on the pore walls exhibits unique physical and chemical properties.Therefore,studying the adsorption morphology of hydrocarbon components in nanometer-sized pores and clarifying the exploitation limits of shale oil at the microscopic level are of great practical significance for the efficient development of continental shale oil.In this study,molecular dynamics simulations were employed to investigate the adsorption characteristics of various single-component shale oils in inorganic quartz fissures,and the influence of pore size and shale oil hydrocarbon composition on the adsorption properties in the pores was analyzed.The results show that different molecules have different adsorption capacities in shale oil pores,with lighter hydrocarbon components(C6H14)exhibiting stronger adsorption abilities.For the same adsorbed molecule,the adsorption amount linearly increases with the increase in pore diameter,but larger pores contribute more to shale oil adsorption.In shale pores,the thickness of the adsorption layer formed by shale oil molecules ranges from 0.4 to 0.5 nm,which is similar to the width of alkane molecules.Shale oil in the adsorbed state that is difficult to be exploited is mainly concentrated in the first adsorption layer.Among them,the volume fraction of adsorbed shale oil in 6 nm shale pores is 40.8%,while the volume fraction of shale oil that is difficult to be exploited is 16.2%.

Dynamic simulation of differential accumulation history of deep marine oil and gas in superimposed basin:A case study of Lower Paleozoic petroleum system of Tahe Oilfield,Tarim Basin,NW China

According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.

Molecular dynamics simulation study of nitrogen vacancy color centers prepared by carbon ion implantation into diamond

Nitrogen vacancy(NV)color centers in diamond have useful applications in quantum sensing andfluorescent marking.They can be gen-erated experimentally by ion implantation,femtosecond lasers,and chemical vapor deposition.However,there is a lack of studies of the yield of NV color centers at the atomic scale.In the molecular dynamics simulations described in this paper,NV color centers are pre-pared by ion implantation in diamond with pre-doped nitrogen and subsequent annealing.The differences between the yields of NV color centers produced by implantation of carbon(C)and nitrogen(N)ions,respectively,are investigated.It is found that C-ion implantation gives a greater yield of NV color centers and superior location accuracy.The effects of different pre-doping concentrations(400–1500 ppm)and implantation energies(1.0–3.0 keV)on the NV color center yield are analyzed,and it is shown that a pre-doping concentra-tion of 1000 ppm with 2 keV C-ion implantation can produce a 13%yield of NV color centers after 1600 K annealing for 7.4 ns.Finally,a brief comparison of the NV color center identification methods is presented,and it is found that the error rate of an analysis utiliz-ing the identify diamond structure coordination analysis method is reduced by about 7%compared with conventional identification+methods.

Exploring the molecular mechanism of action of curcumin for the treatment of diabetic retinopathy,using network pharmacology,molecular docking,and molecular dynamics simulation

Background:Based on network pharmacology and molecular docking,the present study investigated the mechanism of curcumin(CUR)in diabetic retinopathy treatment.Methods:Based on the DisGeNET,Swiss TargetPrediction,GeneCards,Online Mendelian Inheritance in Man,Gene Expression Omnibus,and Comparative Toxicogenomics Database,the intersection core targets of CUR and diabetic retinopathy were identified.The intersection target was imported into the STRING database to obtain the protein-protein interaction map.According to the Database for Annotation,Visualization and Integrated Discovery database,the intersected targets were enriched in Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes pathways.Then Cytoscape 3.9.1 is used to make the drug-target-disease-pathway network.The mechanism of CUR and diabetic retinopathy was further verified by molecular docking and molecular dynamics simulation.Results:There were 203 intersecting targets of CUR and diabetic retinopathy identified.1320 GO entries were enriched for GO functions,which were primarily involved in the composition of cells such as identical protein binding,protein binding,enzyme binding,etc.It was found that 175 pathways were enriched using Kyoto Encyclopedia of Genes and Genomes pathway enrichment methods,which were mainly included in the lipid and atherosclerosis,AGE-RAGE signaling pathway in diabetic complications,pathways in cancer,etc.In the molecular docking analysis,CUR was found to have a good ability to bind to the core targets of albumin,IL-1B,and IL-6.The binding of albumin to CUR was further verified by molecular dynamics simulation.Conclusion:As a result of this study,CUR may exert a role in the treatment of diabetic retinopathy through multi-target and multi-pathway regulation,which indicates a possible direction of future research.

Molecular simulation study of the microstructures and properties of pyridinium ionic liquid[HPy][BF_(4)]mixed with acetonitrile

The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.

Dynamics simulation of tertiary amines adsorbing on kaolinite(001) plane

The collecting power of tertiary amines（DRN,DEN and DPN） on kaolinite follows the order of DENDPNDRN.After reacting with DRN,DEN and DPN,the surface potentials of kaolinite increase remarkably,and the recruitments caused by collectors also follow the order of DENDPNDRN.The results of dynamics simulation show that the geometries of substituent groups bonding to N are deflected and twisted,and some of bond angles are changed when tertiary amines cations adsorb on kaolinite（001） surface.Based on the results of dynamics simulations and quantum chemistry calculations,the electrostatic forces between three tertiary amines cations and 4×4×3（001） plane of kaolinite are 1.38×10？7 N（DRN12H＋）,1.44×10-6 N（DEN12H＋）,1.383×10-6 N（DPN12H＋）,respectively.

Molecular dynamics simulation of relationship between local structure and dynamics during glass transition of Mg_7Zn_3 alloy

The rapid solidification process of Mg7Zn3 alloy was simulated by the molecular dynamics method. The relationship between the local structure and the dynamics during the liquid-glass transition was deeply investigated. It was found that the Mg-centered FK polyhedron and the Zn-centered icosahedron play a critical role in the formation of Mg7Zn3 metallic glass. The self-diffusion coefficients of Mg and Zn atoms deviate from the Arrhenius law near the melting temperature and then satisfy the power law. According to the time correlation functions of mean-square displacement, incoherent intermediate scattering function and non-Gaussian parameter, it was found that the β-relaxation in Mg7Zn3 supercooled liquid becomes more and more evident with decreasing temperature, and the α-relaxation time rapidly increases in the VFT law. Moreover, the smaller Zn atom has a faster relaxation behavior than the Mg atom. Some local atomic structures with short-range order have lower mobility, and they play a critical role in the appearance of cage effect in theβ-relaxation regime. The dynamics deviates from the Arrhenius law just at the temperature as the number of local atomic structures begins to rapidly increase. The dynamic glass transition temperature （Tc） is close to the glass transition point in structure （TgStr）.

Local structure of calcium silicate melts from classical molecular dynamics simulation and a newly constructed thermodynamic model

The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of five kinds of Si-O tetrahedra Qi from these two methods was compared with each other and also with the experimental Raman spectra, an excellent agreement was achieved. These not only displayed the panorama distribution of microstructural units in the whole composition range, but also proved that the thermodynamic model is suitable for the utilization as the subsequent application model of spectral experiments for the thermodynamic calculation. Meanwhile, the five refined regions mastered by different disproportionating reactions were obtained. Finally, the distributions of two kinds of connections between Qi were obtained, denoted as Qi-Ca-Qj and Qi-[Ob]-Qj, from the thermodynamic model, and a theoretical verification was given that the dominant connections for any composition are equivalent connections.

Molecular Dynamic Simulation on the Absorbing Process of Isolating and Coating of α-olefin Drag Reducing Polymer

The absorbing process in isolating and coating process of α-olefin drag reducing polymer was studied by molecular dynamic simulation method, on basis of coating theory of α-olefin drag reducing polymer particles with polyurethane as coating material. The distributions of sodium laurate, sodium dodeeyl sulfate, and sodium dodeeyl benzene sulfonate on the surface of α-olefin drag reducing polymer particles were almost the same, but the bending degrees of them were obviously different. The bending degree of SLA molecules was greater than those of the other two surfactant molecules. Simulation results of absorbing and accumulating structure showed that, though hydrophobie properties of surfactant molecules were almost the same, water density around long chain sulfonate sodium was bigger than that around alkyl sulfate sodium. This property goes against useful absorbing and accumulating on the surface of α-olefin drag reducing polymer particles; simulation results of interactions of different surfactant and multiple hydroxyl compounds on surface of particles showed that, interactions of different surfaetant and one kind of multiple hydroxyl compound were similar to those of one kind of surfaetant and different multiple hydroxyl compounds. These two contrast types of interactions also exhibited the differences of absorbing distribution and closing degrees to surface of particles. The sequence of closing degrees was derived from simulation; control step of addition polymerization interaction in coating process was absorbing mass transfer process, so the more closed to surface of particle the multiple hydroxyl compounds were, the easier interactions With isoeyanate were. Simulation results represented the compatibility relationship between surfactant and multiple hydroxyl compounds. The isolating and coating processes of α-olefin drag reducing polymer were further understood on molecule and atom level through above simulation research, and based on the simulation, a referenced theoretical basis was provided for practical optimal selection and experimental preparation of α-olefin drag reducing polymer particles suspension isolation agent.

Research on the Dynamic Model of Spindle Rotation Induced Error Compensation System of Boring and its Simulation

In this paper we address the dynamics of compensation cutting process from both Laplace s frequency domain and the time domain of the first time, using the two computer aided analyzing softwares: MATLAB and SIMULINK. Theoretical analysis and simulation experiments firstly show that not only the systematical stiffness of workpiece, spindle and tools, but also the regenerated coefficient affects the compensation displacement effect. The results show that the SREC is practicable in reality to decease the spindle induced errors in many engineering applications such as hard boring through simulation and the preliminary experiment results.

Molecular Dynamic Simulation of Melting Points of Trans-1,4,5,8- tetranitro-1,4,5,8-tetraazadacalin （TNAD） with Some Propellants

Molecular dynamic simulation was employed to predict the melting points Tm of TNAD/HMX, TNAD/RDX, TNAD/DINA, and TNAD/DNP systems （tans-1,4,5,8- tetranitro-1,4,5,8-tetraazadacalin （TNAD）, dinitropiperazine （DNP）, cyclotetramethylenetetranitroamine （HMX）, cyclotrimethylenetrinitramine （RDX）, and N-nitrodihydroxyethylaminedinitrate （DINA））. Tm was determined from the inflexion point on the curve of mean specific volume vs. temperature. The result shows that the Tm values of TNAD/HMX, TNAD/RDX, and TNAD/DINA systems are 500, 536, and 488 K, respectively. The TNAD/DNP system has no obvious Tm value, which shows the system is insoluble. Using Tm, the solubility of the four systems was analyzed. The radial distribution functions of the four systems were analyzed and the main intermolecular forces between TNAD and other energetic components are short-range interactions. The better the solubility is, the stronger the intermoleenlar interaction is. In addition, the force field energy at different temperature was also analyzed to predict Tm of the four systems.