关于MD运行速度问题

HTZ

在进行TDDFT测试的时候，我发现TDAP的运行速度比ABACUS要快许多。TDAP是基于Siesta做的扩展，所以我对比了Siesta 4.0.2版本和ABACUS 3.1.0版本MD的运行速度。以下是两者完成任务的时间：


Siesta用时：1381s

ABACUS用时：4508s
可以看出有较大的速度差异。测试体系是乙炔分子(C2H2),两者输入的初始结构一致，均采用verlet算法，PBE泛函，进行1000步的md模拟，只用单核计算，详细的输入输出可以参照上传的附件。
由于我不了解siesta，也只对齐了一部分全局的精度参数，所以我无法断言是算法的差异，还是默认输入的差异。但在结果上Siesta确实也能给出相对合理的结果。由于现在TDDFT测试进度普遍较慢，如果是算法或者优化的问题，ABACUS能否有相应的改进？如果是默认精度的问题，要如何调整输入文件来使ABACUS接近Siesta的速度呢?

zdy-AISI

你可以设置ecutwfc=60, gamma_only=1 来优化计算速度，优化后的时间约为SIESTA的两倍（~2561s）。

zdy-AISI

SIESTA的能量截断应该是电荷密度的吧，我们是波函数的。有个4倍的差别