We present our latest advancements of machine-learned potentials (MLPs) based on the neuroevolution potential (NEP) framework introduced in [Fan et al., Phys. Rev. B 104, 104309 (2021)] and their implementation in the open-source package GPUMD.We increase the accuracy of NEP models both by improving the radial functions in the atomic-environment descriptor using a linear combination of Chebyshev basis functions and by extending the angular descriptor with some four-body and five-body contributions as in the atomic cluster expansion approach.We also detail our efficient implementation of the NEP approach in graphics processing units as well as our workflow for the construction of NEP models, and we demonstrate their application in large-scale atomistic simulations.By comparing to state-of-the-art MLPs, we show that the NEP approach not only achieves above-average accuracy but also is far more computationally efficient.These results demonstrate that the GPUMD package is a promising tool for solving challenging problems requiring highly accurate, large-scale atomistic simulations.To enable the construction of MLPs using a minimal training set, we propose an active-learning scheme based on the latent space of a pre-trained NEP model.Finally, we introduce three separate Python packages, GPYUMD, CALORINE, and PYNEP, which enable the integration of GPUMD into Python workflows.
Thermal properties of molybdenum disulfide (MoS2) have recently attracted attention related to fundamentals of heat propagation in strongly anisotropic materials, and in the context of potential applications to optoelectronics and thermoelectrics. Multiple empirical potentials have been developed for classical molecular dynamics (MD) simulations of this material, but it has been unclear which provides the most realistic results. Here, we calculate lattice thermal conductivity of singleand multi-layer pristine MoS2 by employing three different thermal transport MD methods: equilibrium, nonequilibrium, and homogeneous nonequilibrium ones. These methods allow us to verify the consistency of our results and also facilitate comparisons with previous works, where different schemes have been adopted. Our results using variants of the Stillinger-Weber potential are at odds with some previous ones and we analyze the possible origins of the discrepancies in detail. We show that, among the potentials considered here, the reactive empirical bond order (REBO) potential gives the most reasonable predictions of thermal transport properties as compared to experimental data. With the REBO potential, we further find that isotope scattering has only a small effect on thermal conduction in MoS2 and the in-plane thermal conductivity decreases with increasing layer number and saturates beyond about three layers. We identify the REBO potential as a transferable empirical potential for MD simulations of MoS2 which can be used to study thermal transport properties in more complicated situations such as in systems containing defects or engineered nanoscale features. This work establishes a firm foundation for understanding heat transport properties of MoS2 using MD simulations.
Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H 2 ), such as underground H 2 storage and H 2 O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H 2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H 2 O, (b) self-diffusivities of infinitely diluted H 2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/ kg H 2 O, and (c) solubilities of H 2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H 2 O. The force fields used are the TIP4P/2005 for H 2 O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H 2 . Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.