GPUMD: 高效分子动力学模拟的GPU加速包

GPUMD:高效分子动力学模拟的GPU加速包

GPUMD(Graphics Processing Units Molecular Dynamics)是一款由Zheyong Fan开发的开源分子动力学模拟软件包,专为在图形处理单元(GPU)上高效运行而设计。作为一个通用的分子动力学工具,GPUMD不仅能执行传统的分子动力学模拟,还支持先进的机器学习势能,为材料科学和凝聚态物理研究提供了强大的计算能力。

GPUMD的主要特性

GPU加速: GPUMD充分利用GPU的并行计算能力,大幅提升分子动力学模拟的速度。它要求使用计算能力不低于3.5的GPU卡和CUDA 9.0或更新版本的CUDA工具包。
跨平台支持: GPUMD可以在Linux(使用GCC编译器)和Windows(使用MSVC编译器)操作系统上运行,为研究人员提供了灵活的使用环境。
神经进化势能(NEP): GPUMD支持训练和使用一类称为神经进化势能(Neuroevolution Potentials, NEPs)的机器学习势能。这种势能结合了神经网络和进化算法的优势,能够准确高效地描述原子间相互作用。
多种经典势能: 除了支持NEP,GPUMD还实现了多种经典的原子间势能,如Tersoff势、Lennard-Jones势、嵌入原子方法(EAM)势等。
丰富的分析功能: GPUMD提供了多种分析工具,用于计算材料的热力学、动力学和结构性质,如热膨胀、态密度、热导率等。

安装和使用

安装GPUMD相对简单:

克隆或下载GPUMD源代码。
进入src目录,执行make命令进行编译。
编译完成后,src目录中会生成两个可执行文件:gpumd和nep。

使用GPUMD进行模拟时,只需进入包含输入文件的目录,然后运行相应的可执行文件即可:

path/to/gpumd

或

path/to/nep

GPUMD的应用领域

GPUMD在多个研究领域都有广泛应用:

热输运性质研究: GPUMD提供了多种方法计算材料的热导率,包括平衡分子动力学(EMD)、非平衡分子动力学(NEMD)和均匀非平衡分子动力学(HNEMD)方法。
机械性质模拟: 通过GPUMD,研究人员可以模拟材料在各种力学条件下的行为,如弹性变形、塑性变形等。
相变和结构转变: GPUMD能够模拟材料在不同温度和压力下的相变过程,帮助研究人员理解材料的结构演化。
纳米材料模拟: 对于纳米材料如二维材料、纳米线等,GPUMD提供了专门的分析工具,如声子色散关系计算等。
辐照损伤研究: 通过NEP势能,GPUMD能够模拟材料在辐照条件下的行为,这对核材料研究至关重要。

GPUMD Logo

NEP势能:GPUMD的核心创新

神经进化势能(NEP)是GPUMD的一大亮点。NEP结合了神经网络的表达能力和进化算法的全局优化能力,能够高效准确地描述复杂的原子间相互作用。相比传统的机器学习势能,NEP具有以下优势:

高精度: NEP能够准确捕捉原子间相互作用的复杂性,在多种材料系统中展现出优异的性能。
高效率: 得益于其简洁的网络结构和GPU加速,NEP在训练和预测阶段都表现出极高的计算效率。
可转移性: NEP模型具有良好的可转移性,能够处理训练数据中未出现的原子环境。
多功能性: 除了预测能量和力,NEP还可以用于计算偶极矩、极化率等物理量。

GPUMD的生态系统

围绕GPUMD,研究人员开发了多个Python包,进一步扩展了其功能:

calorine: 用于运行和分析GPUMD分子动力学模拟的Python包。
GPUMD-Wizard: 基于ASE(Atomic Simulation Environment)的材料结构处理软件。
gpyumd: GPUMD的Python3接口,用于生成输入文件和处理输出文件。
mdapy: 提供强大的原子轨迹分析工具。
PyNEP: NEP势能的Python接口。
SOMD: 一个与SIESTA DFT代码结合的从头算分子动力学包。

这些工具大大提高了GPUMD的易用性和功能多样性,使研究人员能够更方便地进行复杂的模拟和分析。

GPUMD的最新进展

GPUMD团队一直在不断改进和扩展软件的功能。最新版本(v3.9.5)引入了多项新特性和改进:

灵活的截断方案: 新增use_typewise_cutoff和use_typewise_cutoff_zbl关键字,允许更灵活的截断方案,提高了模拟性能。
电子阻止: 新增electron_stop关键字,用于模拟电子阻止效应,这对研究辐照损伤至关重要。
外电场模拟: 引入add_efield命令,允许基于给定固定电荷添加电场。
Monte Carlo分子动力学: 新增mc关键字,支持在正则系综、半大正则系综和方差约束半大正则系综中进行高效的MCMD模拟。
D3色散校正: 通过dftd3关键字,可以为NEP模型添加D3色散校正。