24小时热门版块排行榜    

查看: 2014  |  回复: 15

菲顶顶

铜虫 (正式写手)

[求助] berry phase相的计算时,POSCAR是如何改变的?是否有一些原则

看了手册上berry phase方法的例子NaF:
1.关于POSCAR
原来的是
NaF
4.5102
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 1
Direct
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.5000000000000000 0.5000000000000000 0.5000000000000000
位移后的是
NaF
4.5102
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 1
Direct
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.5100000000000000 0.5100000000000000 0.4900000000000000
不太明白POSCAR为什么要这样改动,之前看过一个帖子:http://muchong.com/html/201506/9084252.html
但是用它说的好像和这个例子对不上,原子位置的变化是有什么规则和原则的吗?希望大家帮忙解答一下
2.另外,手册上的例子给出的赝势是LDA,Na_sv和F PAW,请问这种算法在赝势的选取上是有什么要求吗?谢谢大家
3.手册上给出最后的结果是Z*是什么意思?不是计算了3个方向上的吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

顶起~~~大家帮帮忙吧

发自小木虫Android客户端
2楼2015-12-03 21:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

3楼2015-12-07 16:28:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

原则就是保证计算精度的前提下越小越好
当然,我是不建议用vasp的,如果可以的话请尽量使用abinit
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
4楼2015-12-07 17:49:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

引用回帖:
4楼: Originally posted by iamikaruk at 2015-12-07 17:49:29
原则就是保证计算精度的前提下越小越好
当然,我是不建议用vasp的,如果可以的话请尽量使用abinit

是说原子位置移动越小越好吗? 那为什么x和y方向是变大了0.01而z方向减小了0.01.而且看到有人说什么找他的参考相,中心对称想之类之类的 ,不太明白这其中的含义,您方便解答一下吗
5楼2015-12-07 20:19:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

引用回帖:
4楼: Originally posted by iamikaruk at 2015-12-07 17:49:29
原则就是保证计算精度的前提下越小越好
当然,我是不建议用vasp的,如果可以的话请尽量使用abinit

而且为什么要说只移动了F原子
6楼2015-12-07 20:20:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
6楼: Originally posted by 菲顶顶 at 2015-12-07 20:20:18
而且为什么要说只移动了F原子...

因为要算的是F原子的Born effective charge。
移动的话可以0.51 0.51 0.51的。
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
7楼2015-12-07 20:40:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

引用回帖:
7楼: Originally posted by iamikaruk at 2015-12-07 20:40:13
因为要算的是F原子的Born effective charge。
移动的话可以0.51 0.51 0.51的。...

哦哦 好的 谢谢,那大家平时所说的寻找参考相是随便这样移动一下得到的吗~~~就是如果平时算铁电体的极化率的话, 移动哪个原子,还是整体移动,这是改怎么考虑,看文献的话有些看不太明白,因为我刚接触这方面,麻烦了
8楼2015-12-07 20:46:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
8楼: Originally posted by 菲顶顶 at 2015-12-07 20:46:07
哦哦 好的 谢谢,那大家平时所说的寻找参考相是随便这样移动一下得到的吗~~~就是如果平时算铁电体的极化率的话, 移动哪个原子,还是整体移动,这是改怎么考虑,看文献的话有些看不太明白,因为我刚接触这方面,麻 ...

参考相你要找中心对称相。
算铁电极化必须用弛豫的铁电相的原子位置进行计算出formal polarization,然后再用它减去参考相的formal polarization。

» 本帖已获得的红花(最新10朵)

[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
9楼2015-12-07 21:54:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

送红花一朵
引用回帖:
9楼: Originally posted by iamikaruk at 2015-12-07 21:54:41
参考相你要找中心对称相。
算铁电极化必须用弛豫的铁电相的原子位置进行计算出formal polarization,然后再用它减去参考相的formal polarization。...

恩恩 只个知道的,就是说这个中心对称相是怎么得到的,该怎么移动原子位置什么的,是该怎么办呢
10楼2015-12-07 22:27:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 菲顶顶 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 16:45 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:04 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 14:55 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 14:35 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 12:53 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 12:40 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见