A high-efficiency neuroevolution potential for tobermorite and calcium silicate hydrate systems with ab initio accuracy

Tobermorite and Calcium Silicate Hydrate (C-S-H) systems are indispensable cement materials but still lack a satisfactory interatomic potential with both high accuracy and high computational efficiency for better understanding their mechanical performance. Here, we develop a Neuroevolution Machine Learning Potential (NEP) with Ziegler-Biersack-Littmark hybrid framework for tobermorite and C-S-H systems, which conveys unprecedented efficiency in molecular dynamics simulations with substantially reduced training datasets. Our NEP model achieves prediction accuracy comparable to DFT calculations using just around 300 training structures, significantly fewer than other existing machine learning potentials trained for tobermorite. Critically, the GPU-accelerated NEP computations enable scalable simulations of large tobermorite systems, reaching several thousand atoms per GPU card with high efficiency. We demonstrate the NEP's versatility by accurately predicting mechanical properties, phonon density of states, and thermal conductivity of tobermorite. Furthermore, we extend the NEP application to large-scale simulations of amorphous C-S-H, highlighting its potential for comprehensive analysis of structural and mechanical behaviors under various realistic conditions.