| 查看: 226 | 回复: 2 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】优化Mg2Nii结构发生变化
|
|||
|
最近想算下Mg2Ni的声子谱，要求能量和力的收敛标准分别为1e-8 eV和1e-6，但在VASP中优化结构时发现，随着收敛判据的提高，Mg2Ni空间群发生变化，由 P6222 (SG Number: 180 )变成了C222 (SG Number: 21 ). 请问下，这种情况算声子时应该采用哪种结构？用C222的话，收敛要求达到了，但结构发生了变化；采用原来的结构P6222，计算声子的精度又达不到要求？这样的情况应该如何解决？恳请知道的赐教！ 附上INCAR和KPOINTS Automatic mesh > 0 > Monkhorst Pack > 8 8 8 > 0 0 0 > INCAR > SYSTEM = Mg2Ni-180 > PREC = Accurate > ENCUT = 350.000 > IBRION = 1 > NSW = 150 > ISIF = 3 > EDIFF = 1.0E-08 > EDIFFG = 1.0E-06 > ISMEAR = 0 > SIGMA = 0.05 > IALGO = 38 > LREAL = .FALSE. > ADDGRID = .TRUE. |
回复此楼» 猜你喜欢
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
多组分精馏求助
已经有6人回复
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):谢谢详细指点 12-24 15:58
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):谢谢详细指点 12-24 15:58
|
看你的INCAR你做的case应该是0GPa的，一种情况可能是因为你的初始结构远离0GPa，从而需要向0GPa优化，如果压力差的较大，优化时，离子和格矢都变化，这样结构有可能跑了。 你可以先做一下静态的scf，然后再看一下三轴的stress是否相等。 你在那个external pressure的压力下，现优化一下，用isif=2来做，如果找到的空间群是180就没问题。 然后你再一点一点向0压优化，结构应该就没问题了。 也有可能是你找空间群的程序有点问题。 再用,如果你的结构a=b=c的话,你完全就可以用isif来优化, 然后用拟和物态方程的方法来得到你想做的压力点的晶格常数. 你必须要收敛测试,k点和encut都得测试! |
2楼2009-12-24 15:52:09
3楼2009-12-25 13:46:21