24小时热门版块排行榜

返回列表

freewain

铜虫 (初入文坛)

应助: 3 (幼儿园)
金币: 293.9
散金: 59
红花: 3
帖子: 39
在线: 35.4小时
虫号: 987873
注册: 2010-04-02

[求助] 计算体模量时，晶格常数该如何变化？

若是从模拟实验的角度出发，体积的变化不应该超过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

回复此楼

» 猜你喜欢

1楼 2011-12-06 03:24:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

关于势函数的求助已经有23人回复
求助castep计算中：晶胞总能量E 与晶胞体积V 的关系图如何得到，看内容已经有24人回复
采用murn.f计算体积弹性模量不正确？已经有17人回复
【求助】两种对物质加压的作法相同么已经有9人回复
【求助】CASTEP计算弹性常数得到应力应变曲线已经有17人回复
【求助】一定压强下的弹性常数计算方法已经有25人回复
【分享】各种仪器分析的基本原理及谱图表示方法大家一起共享！！！已经有23人回复
【讨论】用vasp计算过弹性常数的请进！！更新中！！已经有51人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 freewain 的主题更新

返回列表