Equilibrium Molecular Dynamics Research Articles

Hydrodynamic slip characteristics of shear-driven water flow in nanoscale carbon slits.

This paper reports on the effects of shear rate and interface modeling parameters on the hydrodynamic slip length (LS) for water-graphite interfaces calculated using non-equilibrium molecular dynamics. Five distinct non-bonded solid-liquid interaction parameters were considered to assess their impact on LS. The interfacial force field derivations included sophisticated electronic structure calculation-informed and empirically determined parameters. All interface models exhibited a similar and bimodal LS response when varying the applied shear rate. LS in the low shear rate regime (LSR) is in good agreement with previous calculations obtained through equilibrium molecular dynamics. As the shear rate increases, LS sharply increases and asymptotes to a constant value in the high shear regime (HSR). It is noteworthy that LS in both the LSR and HSR can be characterized by the density depletion length, whereas solid-liquid adhesion metrics failed to do so. For all interface models, LHSR calculations were, on average, ∼28% greater than LLSR, and this slip jump was confirmed using the SPC/E and TIP4P/2005 water models. To address the LS transition from the LSR to the HSR, the viscosity of water and the interfacial friction coefficient were investigated. It was observed that in the LSR, the viscosity and friction coefficient decreased at a similar rate, while in the LSR-to-HSR transition, the friction coefficient decreased at a faster rate than the shear viscosity until they reached a new equilibrium, hence explaining the LS-bimodal behavior. This study provides valuable insights into the interplay between interface modeling parameters, shear rate, and rheological properties in understanding hydrodynamic slip behavior.

Simulation of interlayer coupling for electroactive covalent organic framework design.

Porous, stacked two-dimensional covalent organic frameworks (2D COFs) bearing semiconducting linkers can support directional charge transfer across adjacent layers of the COF. To better inform the current and possible future design rules for enhancing electron and hole transport in such materials, an understanding of how linker selection and functionalization affects interlayer electronic couplings is essential. We report electronic structure simulation and analysis of electronic couplings across adjacent linker units and to encapsulated species in functionalized electroactive 2D COFs. The detailed dependence of these electronic couplings on interlayer interactions is examined through scans along key interlayer degrees of freedom and through configurational sampling from equilibrium molecular dynamics on semiempirical potential energy surfaces. Beyond affirming the sensitivity of the electronic coupling to interlayer distance and orientation, these studies offer guidance toward linker functionalization strategies for enhancing charge carrier transport in electroactive 2D COFs.

Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators.

Long-range electrostatic effects are fundamental for describing chemical reactivity in the condensed phase. Here, we present the methodology of an efficient quantum mechanical/molecular mechanical (QM/MM) model in periodic boundary conditions (PBC) compatible with QM/MM boundaries at chemical bonds. The method combines electrostatic potential fitted charge operators and electrostatic potentials derived from the smooth particle-mesh Ewald (PME) sum approach. The total energy and its analytic first derivatives with respect to QM, MM, and lattice vectors allow QM/MM molecular dynamics (MD) in the most common thermodynamic ensembles. We demonstrate the robustness of the method by performing a QM/MM MD equilibration of methanol in water. We simulate the cis/trans isomerization free-energy profiles in water of proline amino acid and a proline-containing oligopeptide, showing a correct description of the reaction barrier. Our PBC-compatible QM/MM model can efficiently be used to study the chemical reactivity in the condensed phase and enzymatic catalysis.

Thermal conductivity of irregularly shaped nanoparticles from equilibrium molecular dynamics.

The computation of thermal conductivity for finite nanoparticulate systems, particularly those of irregular shapes, poses significant challenges. Much of the previous work on the thermal conductivity of nanoparticles has been carried out by applying traditional nonequilibrium molecular dynamics (NEMD) methods. One of our recent works [Physical Review B 103, 035417 (2021)] proposed that equilibrium molecular dynamics (EMD) methods can be used for the simulation of thermal conductivity of finite-scale systems and demonstrated their equivalence to NEMD methods. In this work, we explored the application of the EMD approach to calculating the thermal conductivity of zero-dimensional nanoparticles. In our initial step, we merged both methodologies to substantiate the equivalence in thermal conductivity calculation for cube and cylinder nanoparticles. Post data filtration, the usefulness of EMD for thermal conductivity evaluation of zero-dimensional materials was confirmed. The NEMD method struggles with accuracy in nanoparticle systems characterized by a cross-sectional area change in the direction of transport, whereas in the case where the volume can be determined, EMD methods can be used to calculate the thermal conductivity. In a subsequent study, we used the latest machine learning potential (NEP) to calculate the thermal conductivity of spherical nanoparticles and evaluated the results against the classical Tersoff potential.Ultimately, predictions were made for the thermal conductivity of nanoparticles with various unique shapes in all directions. Collectively, our findings underscore the simplicity and effectiveness of using EMD methods for thermal conductivity calculations in the context of nanoparticles. This paves the way for novel strategies in studying the thermal transport properties of finite nanoparticle systems as well as nanopowders.

Contributions of the individual domains of αIIbβ3 integrin to its extension: Insights from multiscale modeling.

The platelet integrin αIIbβ3 undergoes long-range conformational transitions between bent and extended conformations to regulate platelet aggregation during hemostasis and thrombosis. However, how exactly αIIbβ3 transitions between conformations remains largely elusive. Here, we studied how transitions across bent and extended-closed conformations of αIIbβ3 integrin are regulated by effective interactions between its functional domains. We first carried out μs-long equilibrium molecular dynamics (MD) simulations of full-length αIIbβ3 integrins in bent and intermediate conformations, the latter characterized by an extended headpiece and closed legs. Then, we built heterogeneous elastic network models, perturbed inter-domain interactions, and evaluated their relative contributions to the energy barriers between conformations. Results showed that integrin extension emerges from: (i) changes in interfaces between functional domains; (ii) allosteric coupling of the head and upper leg domains with flexible lower leg domains. Collectively, these results provide new insights into integrin conformational activation based on short- and long-range interactions between its functional domains and highlight the importance of the lower legs in the regulation of integrin allostery.

Insights into the mechanism of peptide fibril growth on gold surface

Understanding the formation of β-fibrils over the gold surface is of paramount interest in nano-bio-medicinal Chemistry. The intricate mechanism of self-assembly of neurofibrillogenic peptides and their growth over the gold surface remains elusive, as experiments are limited in unveiling the microscopic dynamic details, in particular, at the early stage of the peptide aggregation. In this work, we carried out equilibrium molecular dynamics and enhanced sampling simulations to elucidate the underlying mechanism of the growth of an amyloid-forming sequence of tau fragments over the gold surface. Our results disclose that the collective intermolecular interactions between the peptide chains and peptides with the gold surface facilitate the peptide adsorption, followed by integration, finally leading to the fibril formation.

Prediction of phonon properties of cubic boron nitride with vacancy defects and isotopic disorders by using a neural network potential

Cubic boron nitride (c-BN) is a promising ultra-wide bandgap semiconductor for high-power electronic devices. Its thermal conductivity can be substantially modified by controlling the isotope abundance and by the quality of a single crystal. Consequently, an understanding of the phonon transport in c-BN crystals, with both vacancy defects and isotopic disorders at near-ambient temperatures, is of practical importance. In the present study, a neural network potential (NNP) for c-BN has been developed, which has facilitated the investigation of phonon properties under these circumstances. As a result, the phonon dispersion and the three- and four-phonon scattering rates that were predicted with this NNP were in close agreement with those obtained from density-functional theory (DFT) calculations. The thermal conductivities of the c-BN crystals were also investigated, with boron (B) vacancies ranging from 0.0% to 0.6%, by using equilibrium molecular dynamics simulations based on the Green-Kubo formula. These simulations accurately capture vacancy-induced phonon softening, localized vibration modes, and phonon localization effects. As has previously been experimentally prepared, four isotope-modified c-BN samples were selected for analyses in the evaluation of the impact of isotopic disorders. The calculated thermal conductivities aligned well with the DFT benchmarks. In addition, the present study was extended to include a c-BN crystal with a natural abundance of B atoms, which also contained B vacancies. Reasonable thermal conductivities and vibrational characteristics, within the temperature range of 250–500 K, were then obtained.

Electrically regulated thermal conductivity of aramid polymer systems

Aramid polymers, renowned for their electronic insulation and thermal conductive properties, are widely adopted as thermal management materials in power electronics. However, the thermal conductivity of aramid polymers under electric field has not been thoroughly understood. In this study, we investigated the thermal conductivity of amorphous and aligned aramid polymer systems under electrical field utilizing equilibrium molecular dynamics simulations. Simulation results showed that the alignment of polymer chain can significantly enhance the thermal conductivity of aramid polymer systems, achieving up to 10.13 W/m-K. Moreover, polarization of aligned aramid polymer was observed when the applied electric field exceeded 14 V/nm. Interestingly, the thermal conductivity of aligned aramid polymer was selectively modulated by the applied electric field. To unravel the underlying phonon mechanism, the molecular orientation of polymer chains and phonon spectral information were analyzed. Our study provides guidance into understanding thermal transport mechanism and thermal conductivity modulation in polymers.

Calculating High-Pressure PAO4 Viscosity with Equilibrium Molecular Dynamics Simulations

The development of optimized lubricants is hindered by missing knowledge of fluid properties, in particular the viscosity, in the range of extreme pressures and temperatures relevant for application. Molecular dynamics simulations can be used to calculate viscosity, but the necessary computational effort imposes practical limits for high viscosities. In this study, the viscosity of PAO4 oil was extracted from equilibrium molecular dynamics simulations as a function of pressure and temperature reaching viscosities up to 20 Pas. Three calculation methods based on different microscopic expressions for the viscosity were used. The methods exhibit considerably different performance with respect to preciseness and computational efficiency. The highest viscosities were found to be calculated most efficiently via the Stokes–Einstein relation, by computing the diffusion coefficient from the velocity correlation function. This offers a new, more effective route to push viscosity calculations in equilibrium molecular dynamics simulations to higher pressure systems.Graphical

Shock Hugoniot Calculations of Newly Designed Thermoplastic Elastomers and Comparison with Classical Binder Estane.

Our purpose is to design excellent binder candidates used in polymer-bonded explosives (PBX) according to the calculated shock of Hugoniot. Here, we mainly examined the thermoplastic elastomer (TPE) binders commonly used in PBX formulations. Equilibrium molecular dynamics (MD) simulation and mixing rule methods were used to calculate the shock Hugoniot values of 180 newly designed TPEs. We focused on the influence of the polymerization degree, contents, and types of soft and hard segments composed of TPEs on the shock Hugoniot and compared them with the classic PBX binder, Estane. The results show that the hard segment has an effect on the Hugoniot curve, which gradually diminishes as the degree of polymerization increases. The underlying physical mechanism can be attributed to the presence of a large number of hydrogen bonds in hard segment domains. The shock Hugoniot of TPEs also depends on the type of soft segments. The volume compression rate of TPEs decreases with increasing content of hard segments under a given shock. By comparing with Estane, a TPE binder commonly used in PBX, we ultimately chose several new TPEs with the potential to serve as PBX binders in terms of shock performance.

Enhanced phonon transport in sI-type methane hydrate under uniaxial compression strain: Unveiling the opening of an advantageous channel

Unveiling the thermal transport mechanisms of natural gas hydrate under strain conditions is fundamental to the heat transfer of the hydrate reservoir during drilling and production. This work investigated the structural stability and phonon transport mechanisms of methane hydrate using equilibrium molecular dynamics simulations subjected to uniaxial compression strain. The results indicate that the stability of large cages within methane hydrates weakens under strain. However, the presence of methane molecules effectively stabilizes the cage structure, enabling it to remain intact even under significant deformation. Furthermore, the thermal conductivity of methane hydrate demonstrates non-linear behavior, characterized by an initial increase followed by a subsequent decrease with increasing uniaxial compression strain. Particularly, the thermal conductivity of methane hydrate attains a maximum value of 0.87±0.018 W⋅m−1⋅K−1 at a strain of -0.04. The results suggest that the elevated spectral energy density and the prolonged phonon lifetime contribute to the enhanced thermal conductivity of methane hydrate. Interestingly, it was observed that the overlapping of phonon vibration modes between water and methane molecules results in the suppression of the peaks and valleys when the strain reaches -0.04. This overlap eliminates the gap and opens up an advantageous channel for phonon transport, thereby enhancing the heat transport capability of methane hydrate. The opening of an advantageous channel for phonon transport results in an extended phonon lifetime and a subsequent increase in thermal conductivity. This study presents a novel perspective for comprehending the unusual thermal conductivity of methane hydrate under strain conditions.

Thermal conductivity and structural behavior of confined H2 from molecular dynamics simulation

In this work, we perform equilibrium molecular dynamics simulation to study the thermal conductivity of hydrogen molecules (H2) under extreme confinement within graphene nanochannel. We analyze the structural behavior of H2 molecules inside the nanochannel and also examine the effect of nanochannel height, the number of H2 molecules, and temperature of the system on the thermal conductivity. Our results reveal that H2 molecules exhibit a strong propensity for absorption onto the nanochannel wall, consequently forming a dense packed layer in close to the wall. This phenomenon significantly impacts the thermal conductivity of the confined system. We made a significant discovery, revealing a strong correlation between the mass density near the nanochannel wall and the thermal conductivity. This finding highlights the crucial role played by the density near the wall in determining the thermal conductivity behavior. Surprisingly, the average thermal conductivity for nanochannels with a height (h) less than 27 Å exhibited an astonishing increase of over 12 times when compared to the bulk. Moreover, we observe that increasing the nanochannel height, while the number of H2 molecules fixed, leads to a notable decrease in thermal conductivity. Furthermore, we investigate the influence of temperature on thermal conductivity. Our simulations demonstrate that higher temperature enhance the thermal conductivity due to increased phonon activity and energy states, facilitating more efficient heat transfer and higher thermal conductivity. To gain deeper insights into the factors affecting thermal conductivity, we explored the phonon density of states. Studying the behavior of hydrogen in confined environments can offer valuable insights into its transport properties and its potential for industrial applications.

Effects of the crystallinity of fluoropolymer binders components in polymer-bonded explosives on shock Hugoniots: A computational study

Fluoropolymers play a crucial role as binders in polymer-bonded explosive (PBX) formulations. However, there is a lack of clear understanding of the effects of increased fluoropolymer crystallinity on the shock response of PBXs in the service environment. This study investigated the shock Hugoniots of two widely applied fluoropolymer binders: (1) F2314 from China—a copolymer with a molar ratio of vinylidene fluoride (VDF) to chlorotrifluoroethylene (CTFE) of 1:4 and (2) F2313 from the United States, also known as Kel F-800, with a VDF to CTFE molar ratio of 1:3. The Hugoniot curves of both fluoropolymers were calculated based on equilibrium molecular dynamics (MD) and a mixing rule. Furthermore, the corresponding P–V curves were obtained through fitting using the Tait equation of state (EOS). Their calculated parameters, including zero-pressure bulk modulus (κo) and sound velocity (co), agreed well with experimental data. The results reveal that the Hugoniots of amorphous F2314 and F2313 exhibited negligible differences. However, increasing crystallinity significantly impacted the Hugoniot curves of both fluoropolymers, especially for F2314 with high crystallinity. The obtained macroscopic characteristic parameters, namely κo and co, exhibited an exponential dependence on crystallinity. Physically, this phenomenon can be attributed to a reduction in the compressible free volume of the fluoropolymers due to a more orderly chain arrangement. Additionally, under the same compression ratio, the shock temperature of the fluoropolymers increased with the crystallinity, posing potential safety risks to explosives. These findings establish a correlation between the crystallinity of fluoropolymers and the shock properties of PBXs, providing a theoretical reference for the formulation design of fluoropolymer-based PBXs.

Thermal transport characterization of monolayer GaN nanoribbon doped with group IV materials: An equilibrium molecular dynamics study

Bulk gallium nitride (GaN) has been considered as one of the prominent materials for a long time in the application of laser, LED and transistors. Due to its superior electrical characteristics, its low-dimensional hexagonal counterpart is currently drawing attention like graphene. In this study, we have investigated the effect of impurity atoms on the heat energy conduction within this honeycomb nanostructure of GaN to explore its scope of applications. The room temperature thermal conductivity of zigzag GaN nanoribbon (ZGaNNR) is found to be ∼ 11 Wm−1K−1 using the molecular dynamics approach with the Green-Kubo formulation. Various group IV elements namely- carbon, silicon, germanium and tin, have been incorporated within the pristine structure in four distinct doping patterns; like- point (Type-I and Type-II), pair and delta doping, to explore the thermal transport in chemically doped nanostructures. A gradual enhancement in the thermal conductivity is observed for the increasing concentration of carbon and silicon as dopants. On the other hand, selection of germanium and tin as impurity atoms suppresses the heat conduction within the crystal due to the increased low frequency phonon scattering. In all cases, Type-I point doping affects the thermal transport of ZGaNNR significantly. The narrower ribbons impede the phonon transportation due to the strong boundary effects resulting in a lower thermal conductivity. Moreover, phonon transportation is largely suppressed, and heat energy is observed to flow at a slower rate as the temperature of the ribbon rises. We have also studied these doped structures under the influence of point, edge and pair vacancies and observed a decreasing trend in thermal conductivity for each type of vacancy defect patterns introduced in the crystal. Such a study would promote the rising interest in properly characterizing the heat transport of low dimensional nanostructures and it can be an effective technique to tune the thermal conductivity of ZGaNNR for improving thermoelectric performance and thermal management in nanoscale devices.

Structural and thermal analyses in semiconducting and metallic zigzag single-walled carbon nanotubes using molecular dynamics simulations.

Equilibrium molecular dynamics (EMD) simulations have been performed to investigate the structural analysis and thermal conductivity (λ) of semiconducting (8,0) and metallic (12,0) zigzag single-walled carbon nanotubes (SWCNTs) for varying ±γ(%) strains. For the first time, the present outcomes provide valuable insights into the relationship between the structural properties of zigzag SWCNTs and corresponding thermal behavior, which is essential for the development of high-performance nanocomposites. The radial distribution function (RDF) has been employed to assess the buckling and deformation understandings of the (8,0) and (12,0) SWCNTs for a wide range of temperature T(K) and varying ±γ(%) strains. The visualization of SWCNTs shows that the earlier buckling and deformation processes are observed for semiconducting SWCNTs as compared to metallic SWCNTs for high T(K) and it also evident through an abrupt increase in RDF peaks. The RDF and visualization analyses demonstrate that the (8,0) SWCNTs can more tunable under compressive than tensile strains, however, the (12,0) zigzag SWCNTs indicate an opposite trend and may tolerate more tensile than compressive strains. Investigations show that the tunable domain of ±γ(%) strains decreases from (-10%≤ γ ≤+19%) to (-5%≤ γ ≤+10%) for (8,0) SWCNTs and the buckling process shifts to lower ±γ(%) for (12,0) SWCNTs with increasing T(K). For intermediate-high T(K), the λ(T) of (12,0) SWCNTs is high but the (8,0) SWCNTs show certainly high λ(T) for low T(K). The present λ(T, ±γ) data are in reasonable agreement with parts of previous NEMD, GK-HNEMD data and experimental investigations with simulation results generally under predicting the λ(T, ±γ) by the ∼1% to ∼20%, regardless of the ±γ(%) strains, depending on T(K). Our simulation data significantly expand the strain range to -10% ≤ γ ≤ +19% for both zigzag SWCNTs, depending on temperature T(K). This extension of the range aims to establish a tunable regime and delve into the intrinsic characteristics of zigzag SWCNTs, building upon previous work.

Salt Effects on the Mechanical Properties of Ionic Conductive Polymer: A Molecular Dynamics Study.

Functional polymers can be used as electrolyte and binder materials in solid-state batteries. This often requires performance targets in terms of both the transport and mechanical properties. In this work, a model ionic conductive polymer system, i.e., poly(ethylene oxide)-LiTFSI, was used to study the impact of salt concentrations on mechanical properties, including different types of elastic moduli and the viscoelasticity with both nonequilibrium and equilibrium molecular dynamics simulations. We found an encouragingly good agreement between experiments and simulations regarding Young's modulus, bulk modulus, and viscosity. In addition, we identified an intermediate salt concentration at which the system shows high ionic conductivity, high Young's modulus, and short elastic restoration time. Therefore, this study laid the groundwork for investigating ionic conductive polymer binders with self-healing functionality from molecular dynamics simulations.

Glass Transition Temperatures and Thermal Conductivities of Polybutadiene Crosslinked with Randomly Distributed Sulfur Chains Using Molecular Dynamic Simulation.

The thermal conductivities and glass transition temperatures of polybutadiene crosslinked with randomly distributed sulfur chains having different lengths from mono-sulfur (S1) to octa-sulfur (S8) were investigated. The thermal conductivities of the related models as a function of the heat flux autocorrelation function, applying an equilibrium molecular dynamic (EMD) simulation and the Green-Kubo method, were studied for a wide range of temperatures. The influence of the length of sulfur chains, degree of crosslinking, and molar mass of the crosslinker on the glass transition temperature and final values of thermal conductivities were studied. First, the degree of crosslinking is considered constant for the eight simulation models, from mono-sulfur (S1) to octa-sulfur (S8), while the molar mass of the sulfur is increases. The results show that the thermal conductivities of the crosslinked structure decrease with increasing temperature for each model. Moreover, by increasing the lengths of the sulfur chains and the molar weight of the crosslinker, thermal conductivity increases at a constant temperature. The MD simulation demonstrates that the glass transition temperature and density of the crosslinked structure enhance as the length of the sulfur chains and molar mass of the sulfur increase. Second, the molar weight of sulfur is considered constant in these eight models; therefore, the degree of crosslinking decreases with the increase in the lengths of the sulfur chains. The results show that the thermal conductivities of the crosslinked structure decrease with the increase in the temperature for each model. Moreover, by increasing the lengths of sulfur chains and thus decreasing the degree of crosslinking, the trend in changes in thermal conductivities are almost the same for all of these models, so thermal conductivity is constant for a specific temperature. In addition, the glass transition temperature and density of the crosslinked structure decrease.

Theoretical insights into He/CH4 separation by MXene nanopore

Natural gas is currently the main industrial source of helium, thus He/CH4 separation has become increasing important. Nanopores based on two-dimensional materials have been considered as very promising for He/CH4 separation, but the transport mechanism through them is yet unclear. Here, both non-equilibrium molecular dynamics (NEMD) and equilibrium MD (EMD) were employed to study the transport of He/CH4 through nanopores of different MXenes with various pore diameters (d), pore density and functional groups, which showed two main transport mechanisms, molecular sieving and Knudsen diffusion. Small d was found to be crucial for the molecular sieving and consequently, high He/CH4 selectivity (SHe/CH4). The density profiles of gas molecules in EMD simulations were used to predict SHe/CH4 and the permeance for gas passing through different nanopores, yielding results well agreeing with NEMD simulations. This work provides theoretical insights for gas transport through MXene nanopore and may help designing nanopores in the future.

Efficient single-run implementation of generalized Einstein relation to compute transport coefficients: A binary-based time sampling.

A robust and simple implementation of the generalized Einstein formulation using single equilibrium molecular dynamics simulation is introduced to compute diffusion and shear viscosity. The unique features underlying this framework are as follows: (1) The use of a simple binary-based method to sample time-dependent transport coefficients results in a uniform distribution of data on a logarithmic time scale. Although we sample "on-the-fly," the algorithm is readily applicable for post-processing analysis. Overlapping same-length segments are not sampled as they indicate strong correlations. (2) Transport coefficients are estimated using a power law fitting function, a generalization of the standard linear relation, that accurately describes the long-time plateau. (3) The use of a generalized least squares (GLS) fitting estimator to explicitly consider correlations between fitted data points results in a reliable estimate of the statistical uncertainties in a single run. (4) The covariance matrix for the GLS method is estimated analytically using the Wiener process statistics and computed variances. (5) We provide a Python script to perform the fits and automate the procedure to determine the optimal fitting domain. The framework is applied to two fluids, binary hard sphere and a Lennard-Jones near the triple point, and the validity of the single-run estimates is verified against multiple independent runs. The approach should be applicable to other transport coefficients since the diffusive limit is universal to all of them. Given its rigor and simplicity, this methodology can be readily incorporated into standard molecular dynamics packages using on-the-fly or post-processing analysis.

Thermal conductivity of GeTe crystals based on machine learning potentials

GeTe has attracted extensive research interest for thermoelectric applications. In this paper, we first train a neuro-evolution potential (NEP) based on a dataset constructed by ab initio molecular dynamics, with the Gaussian approximation potential (GAP) as a reference. The phonon density of states is then calculated by two machine learning potentials and compared with density functional theory results, with the GAP potential having higher accuracy. Next, the thermal conductivity of a GeTe crystal at 300 K is calculated by the equilibrium molecular dynamics method using both machine learning potentials, and both of them are in good agreement with the experimental results; however, the calculation speed when using the NEP potential is about 500 times faster than when using the GAP potential. Finally, the lattice thermal conductivity in the range of 300 K–600 K is calculated using the NEP potential. The lattice thermal conductivity decreases as the temperature increases due to the phonon anharmonic effect. This study provides a theoretical tool for the study of the thermal conductivity of GeTe.