fix pimd/langevin: An Efficient Implementation of Path Integral Molecular Dynamics in LAMMPS

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Li, Yifan, Gomez, Axel, Cai, Kehan, Zhang, Chunyi, Fu, Li, Jia, Weile, Feldman, Yotam M. Y., Blumer, Ofir, Higer, Jacob, Hirshberg, Barak, Xu, Shenzhen, Kohlmeyer, Axel, Car, Roberto
Format: Preprint
Publié: 2026
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866911448003772416
author Li, Yifan
Gomez, Axel
Cai, Kehan
Zhang, Chunyi
Fu, Li
Jia, Weile
Feldman, Yotam M. Y.
Blumer, Ofir
Higer, Jacob
Hirshberg, Barak
Xu, Shenzhen
Kohlmeyer, Axel
Car, Roberto
author_facet Li, Yifan
Gomez, Axel
Cai, Kehan
Zhang, Chunyi
Fu, Li
Jia, Weile
Feldman, Yotam M. Y.
Blumer, Ofir
Higer, Jacob
Hirshberg, Barak
Xu, Shenzhen
Kohlmeyer, Axel
Car, Roberto
contents Path integral molecular dynamics (PIMD), which maps a quantum particle onto a fictitious classical system of ring polymers and propagates the "beads" of this extended classical system using molecular dynamics, is widely used to capture nuclear quantum effects (NQEs) in molecular simulations. Accurate PIMD calculations typically require a large number of beads and are therefore computationally demanding. While software packages such as i-PI offer comprehensive PIMD functionality, the high efficiency of simulations driven by machine learning interatomic potentials, such as Deep Potential (DP), calls for more efficient PIMD implementations that fully exploit modern massively parallel supercomputers. Here we present fix pimd/langevin, an efficient PIMD implementation in LAMMPS that supports commonly used features and leverages the Message Passing Interface architecture of LAMMPS to achieve high computational efficiency. We demonstrate the usage and validate the correctness of our code using liquid water as a representative example, and provide a comprehensive overview of the supported features. Then we discuss several important technical aspects of the implementation. Using DP simulations of water as a benchmark, we show that our implementation achieves several-fold acceleration compared to i-PI. Finally, we report strong and weak scaling results that demonstrate the favorable parallel performance of our code.
format Preprint
id arxiv_https___arxiv_org_abs_2602_13553
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle fix pimd/langevin: An Efficient Implementation of Path Integral Molecular Dynamics in LAMMPS
Li, Yifan
Gomez, Axel
Cai, Kehan
Zhang, Chunyi
Fu, Li
Jia, Weile
Feldman, Yotam M. Y.
Blumer, Ofir
Higer, Jacob
Hirshberg, Barak
Xu, Shenzhen
Kohlmeyer, Axel
Car, Roberto
Chemical Physics
Path integral molecular dynamics (PIMD), which maps a quantum particle onto a fictitious classical system of ring polymers and propagates the "beads" of this extended classical system using molecular dynamics, is widely used to capture nuclear quantum effects (NQEs) in molecular simulations. Accurate PIMD calculations typically require a large number of beads and are therefore computationally demanding. While software packages such as i-PI offer comprehensive PIMD functionality, the high efficiency of simulations driven by machine learning interatomic potentials, such as Deep Potential (DP), calls for more efficient PIMD implementations that fully exploit modern massively parallel supercomputers. Here we present fix pimd/langevin, an efficient PIMD implementation in LAMMPS that supports commonly used features and leverages the Message Passing Interface architecture of LAMMPS to achieve high computational efficiency. We demonstrate the usage and validate the correctness of our code using liquid water as a representative example, and provide a comprehensive overview of the supported features. Then we discuss several important technical aspects of the implementation. Using DP simulations of water as a benchmark, we show that our implementation achieves several-fold acceleration compared to i-PI. Finally, we report strong and weak scaling results that demonstrate the favorable parallel performance of our code.
title fix pimd/langevin: An Efficient Implementation of Path Integral Molecular Dynamics in LAMMPS
topic Chemical Physics
url https://arxiv.org/abs/2602.13553