| 查看: 924 | 回复: 2 | ||
clq5201234铁杆木虫 (初入文坛)
|
[求助]
求助：用vasp算B4结构的E－V线的时候，出现莫名拐点，困扰了很久，急求解决！！
|
|
计算过程如下： INCAR.relax: SYSTEM=HfN ENCUT=600 ISTART=0 ICHARG=2 ISMEAR=1;SIGMA=0.18 IBRION=2 NSW=100 POTIM=0.2 ISIF=4 EDIFF=1E-5 EDIFFG=-1E-3 PREC=Accurate INCAR.static: SYSTEM=HfN ENCUT=600 ISTART=0 ICHARG=2 ISMEAR=-5 PREC=Accurate K点取19*19*15，以Gamma点为中心。 运行脚本为： for i in `seq -w -53 -1 -70`; do cat > POSCAR <<! hcp $i 0.5000000000000000 -0.86600000000000000 0.000000000000000 0.5000000000000 0.86600000000000000 0.000000000000000 0.000000000000000 0.000000000000000 1.65715000000000000 2 2 Direct 0.33333333 0.66666666 0.00000000 0.66666666 0.33333333 0.50000000 0.33333333 0.66666666 0.37500000 0.66666666 0.33333333 0.87500000 ! cd /public/users/jy002/chenlongqing/clq/B4/Vother cp INCAR.relax INCAR vasp（运行命令） cp CONTCAR POSCAR cp INCAR.static INCAR vasp（运行命令） 也就是，根据参考文献取一系列体积，然后对每个体积进行基矢和原子坐标优化，初始的c/a和u因子都是从参考文献查找确定。 但是结果出现了很头痛的问题： 体积从57到70的优化和能量都很不错很平滑，但是体积56却是一个突变点，能量开始骤降，得到的EV图如下  我查看了体积为56的优化过程，发现vasp把它优化成了正交晶系(base centered orthorhombic)，并且c/a也突变到了1.32（本来应该是1.65左右）。我用了很多方法(比如加入SYMPREC=1E-4，或者调小EDIFFG)，但是都不能让vasp明白我计算的是六方晶系，这点大家能给点意见吗？ 另外，我怀疑是不是在体积为56时对应压强下发生了相变，所以会出现这样的情况，所以我看了下我的相图(相图我是只取体积为57之后的正常的EV线拟合得出的)，发现当B4结构体积为56时，对应的压强为23.2GPa，而相变图上在18GPa的时候B1和B4会发声相变，不知道这点有没有联系，还望大家帮我分析一下，非常感谢啊！十万火急啊！！！拜谢！！！  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
2楼2013-07-25 17:51:28
clq5201234
铁杆木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 5679.9
- 散金: 175
- 帖子: 29
- 在线: 18.5小时
- 虫号: 1461285
- 注册: 2011-10-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2013-07-25 18:10:41
28