求助使用CP2K计算熔盐黏度方法

KaZikl

各位老师，学生使用CP2K，GFN-XTB方法结合Green-Kubo公式进行熔盐黏度的计算，参考了Reliable Viscosity Calculation from Equilibrium Molecular Dynamics Simulations: A Time Decomposition Method这篇文章，模拟运行时长100ps,得出的结果仍然不理想，是否是模拟时长太短了，还是我模拟过程有问题，
（1）请问各位老师，有没有使用G-K公式进行黏度计算的实例或者数据处理的方法，阅读了很多使用第一性原理分子动力学进行熔盐黏度计算的文章，大多只是给出计算的结果和理论，并没有给出计算的过程和吹得过程；
（2）是否应该换用Gromacs,进行CMD计算并结合较长时间模拟来进行，因为GFN-XTB方法也并不算特别精确的方法；

sobereva

GFN-xTB描述熔盐本身的精度靠不靠谱是个关键问题。诸如你可以取一段轨迹，对比像样的泛函结合像样的基组跑的能量曲线和GFN-xTB跑的能量曲线（或者对比原子受力），看看是什么情况，势能面描述的误差如果都很大，就根本没有后面分析的事了。

不要在没有定义的情况下用可能造成混淆的词CMD。如果指基于经典力场的MD，和GFN-xTB的模拟精度比完全取决于力场的好坏，既可以更好也可以更差。但无论如何一个关键优点是在简单形式的力场下MD能轻易跑长两个数量级以上的时间，不需要顾虑采样充分程度方面的问题。

KaZikl

sobereva 发表于 2026-5-10 16:35
GFN-xTB描述熔盐本身的精度靠不靠谱是个关键问题。诸如你可以取一段轨迹，对比像样的泛函结合像样的基组跑 ...

好的卢老师，学生受教了