应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 弹性常数计算中用原包体积还是晶包体积？
101/1返回列表
查看: 2332  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yxz476440843

铁虫 (小有名气)

[求助] 弹性常数计算中用原包体积还是晶包体积？

大家好，我是用原包去计算弹性常数的，侯老师计算弹性常数的公式(见下图）中的体积，我是用原包的体积还是晶包的体积去做计算啊？谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-06-09 22:33:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ufo991

铁虫 (正式写手)
这个能量与体积是不是对应的呀?
赞一下 回复此楼
2楼2012-06-11 09:18:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yxz476440843

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by ufo991 at 2012-06-11 09:18:40
这个能量与体积是不是对应的呀?

不是的，能量是加应变之后的能量与未加应变时的能量之差，体积是平衡体积，但不知是用原包还是晶包的体积？各位大侠帮忙啊
赞一下 回复此楼
3楼2012-06-11 09:26:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-11 10:11:44
引用回帖:
3楼: Originally posted by yxz476440843 at 2012-06-11 09:26:40
不是的，能量是加应变之后的能量与未加应变时的能量之差，体积是平衡体积，但不知是用原包还是晶包的体积？各位大侠帮忙啊...

既然分子是vasp算出的能量，分母当然应该是你用来计算的晶胞体积了。
赞一下 回复此楼
4楼2012-06-11 09:34:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yxz476440843

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by mousekingadv at 2012-06-11 09:34:01
既然分子是vasp算出的能量，分母当然应该是你用来计算的晶胞体积了。...

谢谢你的回复，但我没明白你的意思，我用来计算的是原包而不是晶包，你的意思是采用原包还是晶包体积啊？有什么依据吗？
赞一下 回复此楼
5楼2012-06-12 11:31:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

引用回帖:
5楼: Originally posted by yxz476440843 at 2012-06-12 11:31:56
谢谢你的回复，但我没明白你的意思，我用来计算的是原包而不是晶包，你的意思是采用原包还是晶包体积啊？有什么依据吗？...

如果我没理解错的话你计算用的晶胞就是这个结构的“原胞”（即最小重复单元）
好像不举例解释不清...贴POSCAR吧...
赞一下 回复此楼
6楼2012-06-12 14:02:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yxz476440843

铁虫 (小有名气)
引用回帖:
6楼: Originally posted by mousekingadv at 2012-06-12 14:02:27
如果我没理解错的话你计算用的晶胞就是这个结构的“原胞”（即最小重复单元）
好像不举例解释不清...贴POSCAR吧......

是的，呵呵，我用的就是原包，下面是POSCAR，谢谢！
xxxx
$i  
-1.75100000000000     1.75100000000000    5.27926500000000
1.75100000000000    -1.75100000000000    5.27926500000000
1.75100000000000     1.75100000000000   -5.27926500000000
2 1 2 1
Direct
0.3596999943256378    0.3596999943256378   -0.0000000000000000
-0.3596999943256378   -0.3596999943256378   0.0000000000000000
0.5000000000000001    0.5000000000000000    1.0000000000000000
0.2500000000000000    0.7500000000000000    0.5000000000000000
0.7500000000000000    0.2500000000000000    0.5000000000000000
0.0000000000000000    0.0000000000000000    0.0000000000000000
赞一下 回复此楼
7楼2012-06-13 12:55:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
WDD880227: 金币+1, 多谢您的指导 2012-06-14 08:40:06
yxz476440843: 金币+10, ★★★很有帮助, 懂了，非常感谢您热心的帮助！ 2012-06-15 09:33:19
引用回帖:
7楼: Originally posted by yxz476440843 at 2012-06-13 12:55:41
是的，呵呵，我用的就是原包，下面是POSCAR，谢谢！
xxxx
$i  
-1.75100000000000     1.75100000000000    5.27926500000000
1.75100000000000    -1.75100000000000    5.27926500000000
1.7510000000000 ...

类似这个就是这个体系的最简晶胞，称为原胞。有时候为了精细地控制控制掺杂比例会人为扩大计算用晶胞的大小。如果（假设）把x扩大一倍,原子数目也从6个变成12个，这样的晶胞就被称为超胞。

你的计算用晶胞（和原胞大小相同）体积是a点乘(b叉乘c)，或者可以在OUTCAR中读取…………
赞一下 回复此楼
8楼2012-06-13 15:28:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgf9f18

金虫 (小有名气)
是啊，我也碰到了这样的问题，在计算Lavas NbCr2 中的C15结构时，选择晶胞还是元胞都无法确定，用386的电脑天天在测试，CUtoff和Kpoint都没有得到准确值，郁闷。
赞一下 回复此楼
因未知而快乐
9楼2012-06-13 16:24:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)
引用回帖:
9楼: Originally posted by lgf9f18 at 2012-06-13 16:24:37
是啊，我也碰到了这样的问题，在计算Lavas NbCr2 中的C15结构时，选择晶胞还是元胞都无法确定，用386的电脑天天在测试，CUtoff和Kpoint都没有得到准确值，郁闷。

ENCUT取为ENMAX的1.3倍就非常好了。KPOINTS对机时和结果精度都影响巨大，建议在可承受的时间内尽量多取（不过11*11*11以上就意义不大了）
赞一下(1人) 回复此楼
10楼2012-06-15 20:24:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yxz476440843 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭