[论文解读] Accelerating glassy dynamics using graphics processing units

Peter H. Colberg, Felix Höfling|arXiv (Cornell University)|Dec 20, 2009

Theoretical and Computational Physics参考文献 1被引用 2

一句话总结

本文提出了一种使用CUDA实现的GPU加速分子动力学（MD）模拟，可模拟最多一百万个粒子，与串行CPU代码相比，速度最高提升了80倍。通过采用双单精度浮点数模拟，解决了单精度浮点数运算导致的数值不稳定性问题，确保在10⁸个MD步长内严格保持能量守恒，并防止过冷Lennard-Jones混合体系中出现物理上不正确的结果。

ABSTRACT

Modern graphics processing units (GPUs) provide impressive computing resources, which can be accessed conveniently through the CUDA programming interface. We describe how GPUs can be used to considerably speed up molecular dynamics (MD) simulations for system sizes ranging up to about 1 million particles. Particular emphasis is put on the numerical long-time stability in terms of energy and momentum conservation, and caveats on limited floating-point precision are issued. Strict energy conservation over 10^8 MD steps is obtained by double-single emulation of the floating-point arithmetic in accuracy-critical parts of the algorithm. For the slow dynamics of a supercooled binary Lennard-Jones mixture, we demonstrate that the use of single-floating point precision may result in quantitatively and even physically wrong results. For simulations of a Lennard-Jones fluid, the described implementation shows speedup factors of up to 80 compared to a serial implementation for the CPU, and a single GPU was found to compare with a parallelised MD simulation using 64 distributed cores.

研究动机与目标

• 利用现代GPU加速玻璃态体系的分子动力学模拟。
• 尽管GPU浮点数精度有限，仍确保能量和动量守恒的长期数值稳定性。
• 证明单精度算术在过冷液体慢动态模拟中可能产生物理上不正确的结果。
• 仅使用单个GPU即可实现与64核分布式CPU集群相当的性能。

提出的方法

• 利用CUDA编程接口，将MD模拟移植至GPU架构。
• 在对精度要求高的部分实现双单精度算术，以模拟更高精度的浮点运算。
• 采用二元Lennard-Jones混合体系作为测试系统，评估长时间稳定性和物理准确性。
• 与串行CPU代码及64核并行化CPU实现进行性能基准测试。
• 应用标准MD积分方法，并对力计算进行优化，以适配GPU的内存层次结构与并行性。

实验结果

研究问题

• RQ1 GPU加速能否在玻璃形成体系的分子动力学模拟中实现显著的速度提升？
• RQ2 单精度浮点数算术是否会在过冷液体的长时间MD模拟中引入物理上不正确的行为？
• RQ3 双单精度模拟能否在GPU上实现10⁸个MD步长内的严格能量守恒？
• RQ4 GPU性能与高度并行化的CPU集群在MD模拟中的表现如何比较？

主要发现

• 与串行CPU实现相比，GPU实现最高实现了80倍的速度提升。
• 单个GPU的性能可与64核分布式CPU并行化实现相媲美。
• 双单精度模拟实现了10⁸个MD步长内的严格能量守恒，确保了数值稳定性。
• 在过冷二元Lennard-Jones混合体系的模拟中，单精度算术导致了定性和物理上不正确的结果。
• 该方法保持了动量守恒，并在长时间模拟时间尺度上保持稳定。