24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 5741  |  回复: 14

elliot.cc

铁虫 (初入文坛)

[交流] VASP计算极化强度时 铁电相、参考相和顺电相的设定 已有7人参与

用VASP计算极化强度时,分别算出铁电相和参考相的极化,两者按Barry phase方法中的公式算出结果后除以体积就是自发极化强度。但问题是参考相的设定,既然要除以体积,那么参考相和铁电相的体积应该相同才对。①如果已知材料的铁电相,想建立参考相的模型,是否就是在铁电相晶格参数的基础上仅仅修改原子坐标就行,这时铁电相和参考相的a,b,c分别相等,再对两者做等体积结构优化,然后算极化强度?②参考相和顺电相到底是不是同一个相?如果不是,那就按①中的步骤计算极化;如果是,那问题就来了,比如BaTiO3,铁电相是四方结构,顺电相是立方结构,晶格参数完全不同,这又该怎么计算极化呢?   
希望了解相关计算过程的伙伴能给讲解下!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

iamikaruk

木虫 (著名写手)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
ljw4010: 金币+1, 谢谢交流! 2014-12-05 20:14:51
1. 对铁电相做结构优化,然后从优化的结果反推参考相
2. 一般不是一个
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
2楼2014-06-14 00:45:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by jk射手座 at 2014-10-06 15:27:18
(路过看到这里)请问师兄,您所指的反推是什么意思呢,是怎么个反推?我就差参考相的坐标怎么取差一点点想明白。。。...

从非中心对称相强制变成中心对称性的相
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
4楼2014-10-08 10:17:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

jk射手座

新虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by iamikaruk at 2014-06-14 00:45:01
1. 对铁电相做结构优化,然后从优化的结果反推参考相
2. 一般不是一个

(路过看到这里)请问师兄,您所指的反推是什么意思呢,是怎么个反推?我就差参考相的坐标怎么取差一点点想明白。。。
主做电弧增材制造方向
3楼2014-10-06 15:27:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao123zou

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by iamikaruk at 2014-10-08 10:17:38
从非中心对称相强制变成中心对称性的相...

什么叫做强制变成对称中心呢?可否举个例子?拜托了,我也很苦恼
5楼2014-12-05 19:00:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by xiao123zou at 2014-12-05 19:00:19
什么叫做强制变成对称中心呢?可否举个例子?拜托了,我也很苦恼...

BaTiO3的四方相可以强制把Ti和O放在体心和面心位置
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
6楼2014-12-06 12:13:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyj1409

铜虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by iamikaruk at 2014-12-06 12:13:35
BaTiO3的四方相可以强制把Ti和O放在体心和面心位置...

你好   ,您看我理解的强制中心对称对不对?
这是优化后的四方相结构的CONTCAR:
New structure                           
   1.00000000000000     
     3.8468205145753176    0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.8468205145753176    0.0000000000000000
     0.0000000000000000   -0.0000000000000000    3.9284566005661996
   Pb   Ti   O
     1     1     3
Direct
-0.0000000000000000  0.0000000000000000  0.0365403318419011
  0.5000000000000000  0.5000000000000000  0.5601448334707568
  0.5000000000000000  0.5000000000000000  0.0980368213085622
  0.0000000000000000  0.5000000000000000  0.6141389656893873
  0.5000000000000000 -0.0000000000000000  0.6141389656893873
强制变为中心对称的参考相,如下,
New structure                           
   1.00000000000000     
     3.8468205145753176    0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.8468205145753176    0.0000000000000000
     0.0000000000000000   -0.0000000000000000    3.9284566005661996
   Pb   Ti   O
     1     1     3
Direct
-0.0000000000000000  0.0000000000000000  0.0
  0.5000000000000000  0.5000000000000000  0.5
  0.5000000000000000  0.5000000000000000  0.0
  0.0000000000000000  0.5000000000000000  0.5
  0.5000000000000000 -0.0000000000000000  0.5
然后直接用它作为极化计算的POSCAR。您看这么理解对吗
开心果
7楼2015-04-22 11:44:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by sunyj1409 at 2015-04-22 11:44:18
你好   ,您看我理解的强制中心对称对不对?
这是优化后的四方相结构的CONTCAR:
New structure                           
   1.00000000000000     
     3.8468205145753176    0.0000000000000000    0.0 ...

这样是可以的。唯一要注意的是是否满足绝缘的要求
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
8楼2015-04-23 07:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyj1409

铜虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by iamikaruk at 2015-04-23 07:53:56
这样是可以的。唯一要注意的是是否满足绝缘的要求...

多谢!根据这个思路,我昨天构造了PbVO3四方相的参考相试了下,发现这个参考相没有带隙。然后又试了下变形90%的过渡态,还是会出现错误,Error in subroutine BERRY: did not find all determinants。看来我选的参考相有问题。请问,我优化好的四方向应该怎么去选取它的参考相呢?
开心果
9楼2015-04-23 09:31:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by sunyj1409 at 2015-04-23 09:31:56
多谢!根据这个思路,我昨天构造了PbVO3四方相的参考相试了下,发现这个参考相没有带隙。然后又试了下变形90%的过渡态,还是会出现错误,Error in subroutine BERRY: did not find all determinants。看来我选的参 ...

试过选取其他变形吗?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
10楼2015-04-24 04:26:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 elliot.cc 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 一志愿安徽大学085601 +8 zyx上岸! 2026-04-04 8/400 2026-04-05 23:08 by lbsjt
[考研] 0855求调剂材料 +11 红桃灼灼 2026-04-04 11/550 2026-04-05 22:13 by 醉翁wl
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 材料0856 英一数二 323 求调剂 +14 袁sy 2026-04-01 14/700 2026-04-05 18:18 by cql1109
[考研] 一志愿南航,数一英一学硕317求调剂!! +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 293分求调剂,外语为俄语 +6 加一一九 2026-03-31 6/300 2026-04-04 14:57 by 聪明的大松鼠
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 321求调剂 一志愿 浙江工业大学生物医药 +5 嘿嘿HC 2026-04-01 6/300 2026-04-02 15:23 by sophie2180
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
信息提示
请填处理意见