求助: 不能做出和文献工作相似的结果 (cp2k进行aimd测试)

Ron3319

小弟尝试重复文献的工作来检查自己的cp2k计算是否正确｡  主要看在差不用的设定下，我的RDF是否和文献工作差不多。文献是"Structure, dynamics and vibrational spectrum of supercritical CO2/H2O mixtures from ab initio molecular dynamics as a function of water cluster formation" Structure, dynamics and vibrational spectrum of supercritical CO2/H2O mixtures from ab initio molecular dynamics as a function of water cluster formation - Physical Chemistry Chemical Physics (RSC Publishing)

做32个CO2在318.15K时的计算。使用基组为DZVP-MOLOPT-SR-GTH，CUTOFF为300Ry，REL_CUTOFF为50 Ry

Pair Potential:
Type DFTD3
REFERENCE FUNCTIONAL PBE

但我的计算结果和文献的RDF有明显差异。虽然文献的工作是跑了42或52ps，而我的只跑了约22ps，但这差异应该不会在20ps有大改善。

另一个问题是我的RDF没有归一，32CO2在14.29 x 14.29 x 14.29A的立方体中，但RDF的R距离只到7.06A。

可否请大家告诉我那个地方设定有误

文件已附在这帖中。


抱歉out和xyz文件过大，压缩后也未能上传

sobereva

设置没明显问题
跑的时间太短，rdf只应统计平衡后的轨迹
弄清楚文献里用的是DZVP-GTH还是DZVP-MOLOPT-SR-GTH，差异极其悬殊

诸如Test-1.ener这样尺寸很大的文本文件上传前记得压缩