计算化学公社

标题: 求助：在PIMD下加压后平衡体积不准确 [打印本页]

作者
Author: Dragonmk3 时间: 2022-7-15 21:49
标题: 求助：在PIMD下加压后平衡体积不准确
使用的程序是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，但不能解决问题。

作者
Author: k64_cc 时间: 2022-7-16 09:39
不太懂，碳的NQE有影响吗……？

作者
Author: Dragonmk3 时间: 2022-7-16 12:57

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)