| 查看: 768 | 回复: 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会发声相变，不知道这点有没有联系，还望大家帮我分析一下，非常感谢啊！十万火急啊！！！拜谢！！！  |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有24人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
2楼2013-07-25 17:51:28
clq5201234
铁杆木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 5701.9
- 散金: 155
- 帖子: 28
- 在线: 18.5小时
- 虫号: 1461285
- 注册: 2011-10-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2013-07-25 18:10:41