| 查看: 818 | 回复: 0 | ||
[求助]
计算体模量时，晶格常数该如何变化？
|
|
若是从模拟实验的角度出发，体积的变化不应该超过1%，也就是晶格常数的取值以 0.002 A 增加。 但是，按照VASP手册中的说法，为了消除 Puley stress 的影响，体积的变化要大一些，有的教程推荐使用 0.1 A 的间隔。 在我的计算中，当使用 0.002 A 时，同一 K-point 取值下的一系列计算结果可以很好的拟合出 Birch-Murnaghan 方程。但是计算出的体模量并不能随着 K-point 取值的增加而稳定/收敛。 如果使用 0.1 A 时，则甚至连同一 K-point 取值下的一系列计算都不能很好的拟合出状态方程。 请问，按照大家的计算经验，应该如何选择？谢谢！ P.S. 我刚开始使用 VASP，现在主要进行一些小的计算上手。上面的结果来自于计算 NiFe 合金。 具体的输入文件如下： #! /bin/bash for i in `seq -w 3 1 15` do mkdir $i cd $i for j in `seq -w 3.0 0.1 4.0` do mkdir $j cd $j cat > KPOINTS << EOF automatic k-mesh generation 0 G $i $i $i 0 0 0 EOF cat > POSCAR << EOF fcc Ni_3_Fe_1 $j 1 0 0 0 1 0 0 0 1 3 1 c 0 0 0 0 0.5 0.5 0.5 0 0.5 0.5 0.5 0 EOF cat > INCAR << EOF SYSTEM = Ni_3_Fe_1 ISTART = 0 job: 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ENCUT = 350 eV ISMEAR = 2 SIGMA = 0.1 LREAL = .FALSE. PREC = Normal NGX = 16 NGY = 16 NGZ = 16 ISPIN = 2 MAGMOM = 1 1 1 5 NPAR = 2 IBRION = 2 NSW = 300 EOF |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
关于势函数的求助
已经有23人回复
- 求助castep计算中：晶胞总能量E 与晶胞体积V 的关系图如何得到，看内容
已经有24人回复
- 采用murn.f计算体积弹性模量不正确？
已经有17人回复
- 【求助】两种对物质加压的作法相同么
已经有9人回复
- 【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
- 【求助】一定压强下的弹性常数计算方法
已经有25人回复
- 【分享】各种仪器分析的基本原理及谱图表示方法大家一起共享！！！
已经有23人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有52人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
