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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算极化强度时 铁电相、参考相和顺电相的设定
151/2返回列表12下一页
查看: 6129  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

elliot.cc

铁虫 (初入文坛)

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

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

» 收录本帖的淘帖专辑推荐

磁性计算 soft study

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-13 15:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

iamikaruk

木虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
ljw4010: 金币+1, 谢谢交流！ 2014-12-05 20:14:51
1. 对铁电相做结构优化，然后从优化的结果反推参考相
2. 一般不是一个
赞一下(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
（路过看到这里）请问师兄，您所指的反推是什么意思呢，是怎么个反推？我就差参考相的坐标怎么取差一点点想明白。。。...

从非中心对称相强制变成中心对称性的相
赞一下(1人) 回复此楼
[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 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了，不知道这次命运如何 +3 雨打竹帘 2026-08-30 3/150 2026-08-30 09:43 by KNIGHT2011
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:43 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 06:13 by ZPa0EcMwuECS
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:51 by ZPa0EcMwuECS
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭