求助：在PIMD下加压后平衡体积不准确

Dragonmk3

使用的程序是JAEA开发的PIMD code。
我想用路径积分方法计算弹性模量之类的东西。
我的体系是一个2x2x2的碳12金刚石，64个原子。用VASP计算作用势，赝势用paw-pbe，NTT系综，MNHC控温，三个主应力加10GPa。温度控制在300K，珠子的初始温度设在500K，珠子数设为4，步长试过1fs和0.25fs，跑5000步。
简单的说就是同一个程序在加压前后的两次结构优化以及计算普通的aimd时得到的结果都算靠谱，但是当珠子数大于1后，体系明显变硬了很多。加压前结构优化结果与平衡后的体积差别不算太大；在加压后的结构优化时盒子体积平衡在2432 bohr左右，但是有了珠子数以后加压体积平衡在2475 bohr左右，导致计算的体积模量基本在一千以上，而正常应该在四五百左右。

试过把珠子数设到16，但不能解决问题。

k64_cc

不太懂，碳的NQE有影响吗……？

Dragonmk3

k64_cc 发表于 2022-7-16 09:39
不太懂，碳的NQE有影响吗……？

哪怕是没影响，这个结果也不对啊