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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于born effective charge
591/2返回列表12下一页
查看: 6194  |  回复: 58
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

spur

金虫 (正式写手)

[交流] 关于born effective charge

Born Effective Charge (BEC)张量 是铁电性质研究中一个很重要的物理量。根据公式，通过将BEC与相应原子的位移相乘可以得到材料的极化值大小。但在做乘积的时候，BEC用全张量还是只用对角元呢？BEC张量非对角元有没有什么意义？我此前对BiFeO3单胞做了一个测试，如果用BEC的全张量，得到的极化值为63.1 uC/cm^2;如果用BEC的对角元做乘积，得到的极化值为92.1 uC/cm^2，该值与以往发表的数值相吻合。请大家提提建议。谢谢!

各位，我把我计算BiFeO3的结果以及一篇BiFeO3的文章放在附件了，有兴趣的朋友帮我看看吧。[ Last edited by spur on 2012-5-8 at 20:49 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : L. Q. Chen-2011_Acta_BFO.pdf
    • 2012-05-08 20:29:37, 531.29 K
  • 附件 2 : BEC and displacement.doc
    • 2012-05-08 20:29:51, 37 K

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

VASP

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-05-08 14:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有9个 )

iamikaruk

木虫 (著名写手)
★
贺仪: 金币+1, 谢谢指教！ 2012-05-09 09:46:38
贺仪: 回帖置顶 2012-05-09 14:30:49
uuv2010: , 帖子编辑上面有latex可以让公式更好看，呵呵 2012-05-12 10:39:22
引用回帖:
20楼: Originally posted by spur at 2012-05-08 20:55:03:
谢谢你的讲解！不知道能不能请你帮忙看看我的计算结果或者另外一篇文章(首页附件)？如果这个结果是由于计算方法的问题，那我可能得用berry phase做。但如果用berry phase，怎么计算BEC的非对角元呢？谢谢！

用berry phase无外乎也是一次给一个原子小的位移 (i=1,2,3)，考虑对称性的话，你可以选取特定的位移使得可以同时计算不同分量上的BEC。用berry phase计算极化量 (i=1,2,3)，BEC电荷张量就是简单的了

[ Last edited by iamikaruk on 2012-5-12 at 10:43 ]
赞一下 回复此楼
21楼2012-05-08 21:04:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
贺仪: 回帖置顶 2012-05-09 14:30:57
引用回帖:
22楼: Originally posted by iamikaruk at 2012-05-08 21:12:06:
我大概看了一下你的计算结果，BEC没有太大问题，倒是你的位移量是怎么来的呢？保持晶格常数不变只弛豫内部离子坐标吗？

我是对晶格常数以及原子坐标都进行了优化(ISIF=3)来得到Final structure的，然后采用Final structure的晶格常数来构建Initial structure，让Initial structure的原子都处在高对称的立方相位置，比如0， 0.25， 0.5， 1.0这些位置。

如果BEC、位移量都没问题的话，为什么算出来的极化值跟以往别人用berry phase算得的结果差那么远呢？
赞一下 回复此楼
23楼2012-05-08 22:40:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★
贺仪: 金币+1, 谢谢指教！ 2012-05-09 09:46:55
贺仪: 回帖置顶 2012-05-09 14:31:01
引用回帖:
23楼: Originally posted by spur at 2012-05-08 22:40:13:
我是对晶格常数以及原子坐标都进行了优化(ISIF=3)来得到Final structure的，然后采用Final structure的晶格常数来构建Initial structure，让Initial structure的原子都处在高对称的立方相位置，比如0， 0.25， 0 ...

有一个可能就是你的原子位移量比较大，已经偏离了Born effective charge - 自发极化的线性范围了。
另外你的initial structure是Pm-3m还是R-3c？
赞一下 回复此楼
27楼2012-05-08 23:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
贺仪: 回帖置顶 2012-05-09 14:31:12
引用回帖:
27楼: Originally posted by iamikaruk at 2012-05-08 23:13:16:
有一个可能就是你的原子位移量比较大，已经偏离了Born effective charge - 自发极化的线性范围了。
另外你的initial structure是Pm-3m还是R-3c？

有一个可能就是你的原子位移量比较大，已经偏离了Born effective charge - 自发极化的线性范围了

这个没理解，能说得更具体吗？

Initial Structure中的原子处在立方相的位置，是Pm-3m对称性。
赞一下 回复此楼
38楼2012-05-09 08:20:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★
贺仪: 金币+2, 谢谢指教！ 2012-05-09 09:47:26
贺仪: 回帖置顶 2012-05-09 14:31:19
引用回帖:
38楼: Originally posted by spur at 2012-05-09 08:20:21:
有一个可能就是你的原子位移量比较大，已经偏离了Born effective charge - 自发极化的线性范围了

这个没理解，能说得更具体吗？

Initial Structure中的原子处在立方相的位置，是Pm-3m对称性。

用BEC来估算自发极化用的是线性近似，即\Delta P=Z*\delta，但是位移量很大的时候，这个线性近似就不成立了。
Pm-3m对称性的单胞从晶格常数上就跟R3c存在很大的区别，我建议你用一个跟R3c单胞大小完全相同的R-3c作为参考结构计算位移量
赞一下 回复此楼
44楼2012-05-09 09:01:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
贺仪: 回帖置顶 2012-05-09 14:31:21
引用回帖:
44楼: Originally posted by iamikaruk at 2012-05-09 09:01:21:
用BEC来估算自发极化用的是线性近似，即\Delta P=Z*\delta，但是位移量很大的时候，这个线性近似就不成立了。
Pm-3m对称性的单胞从晶格常数上就跟R3c存在很大的区别，我建议你用一个跟R3c单胞大小完全相同的R-3 ...

可能我此前没说清楚。我的Initial Structure的晶格常数跟R-3c是一样的，只是所有的原子处在立方结构的位置。这个结构可能就是你说的R3c对称性。
赞一下 回复此楼
45楼2012-05-09 09:19:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★
贺仪: 回帖置顶 2012-05-09 14:31:23
贺仪: 金币+1, 谢谢指教！ 2012-05-09 14:31:32
引用回帖:
45楼: Originally posted by spur at 2012-05-09 09:19:17:
可能我此前没说清楚。我的Initial Structure的晶格常数跟R-3c是一样的，只是所有的原子处在立方结构的位置。这个结构可能就是你说的R3c对称性。

明白你的意思了。
不过我看到你结果里面的位移有些大的吓人，应该有0.4埃？如果不介意的话，你能把R-3c和R3c结构的cif文件发给我吗？我研究看看
赞一下 回复此楼
46楼2012-05-09 09:55:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
贺仪: 回帖置顶 2012-05-09 17:36:58
引用回帖:
46楼: Originally posted by iamikaruk at 2012-05-09 09:55:20:
明白你的意思了。
不过我看到你结果里面的位移有些大的吓人，应该有0.4埃？如果不介意的话，你能把R-3c和R3c结构的cif文件发给我吗？我研究看看

不止呢，有些0.6埃。附件为两个文件，麻烦你帮忙看看。谢谢！
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : Initial-rBiFeO3.cif
    • 2012-05-09 10:32:41, 2.04 K
  • 附件 2 : Final-rBiFeO3.cif
    • 2012-05-09 10:32:49, 2.75 K
47楼2012-05-09 10:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★
贺仪: 回帖置顶 2012-05-09 17:37:02
贺仪: 金币+1, 谢谢指教！ 2012-05-09 17:37:11
引用回帖:
47楼: Originally posted by spur at 2012-05-09 10:33:01:
不止呢，有些0.6埃。附件为两个文件，麻烦你帮忙看看。谢谢！

大概算了一下，原子太多，好麻烦
不过我大概猜测一下，你给出的Born effective charge应该是笛卡尔坐标下的，而位移用的不是笛卡尔坐标，会不会你这里没有考虑这个问题？
赞一下 回复此楼
48楼2012-05-09 15:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

iamikaruk

木虫 (著名写手)
★ ★
spur(金币+1): 谢谢参与
御剑江湖: 金币+1, 谢谢～～～ 2012-05-08 19:49:33
必须是全张量

其中是原子的BEC，是原子沿j分量的位移
Neaton有篇文章不是说极化值在0~100都有可能，我了个囧……

[ Last edited by iamikaruk on 2012-5-12 at 10:45 ]
赞一下 回复此楼
2楼2012-05-08 15:31:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuanyeguhong

至尊木虫 (文坛精英)
★
spur(金币+1): 谢谢参与
谢谢楼主！
回复此楼
3楼2012-05-08 15:52:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
引用回帖:
2楼: Originally posted by iamikaruk at 2012-05-08 15:31:40:
必须是全张量
P_i=\sum_{\kappa,j} Z_{\kappa,ij}*\delta_{\kappa,j}
其中Z_{\kappa,ij}是原子\kappa的BEC，\delta_{\kappa,j}是原子\kappa沿j分量的位移
Neaton有篇文章不是说极化值在0~100都有可能，我了个囧 ...

谢谢你的回复。那能不能说一下非对角元的意义呢？另外，Neaton的文章是实验方面的吧？实验上各种因素都可能影响这个极化值。但理论计算如果差距那么大，貌似不太合理。
赞一下 回复此楼
4楼2012-05-08 16:02:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sesame_oil

荣誉版主 (文坛精英)
★
spur(金币+1): 谢谢参与
来学习了
回复此楼
5楼2012-05-08 16:15:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★
spur: 金币+5, 谢谢帮助！ 2012-05-08 19:46:49
fzx2008: 金币+2, 谢谢指教 2012-05-08 20:58:14
引用回帖:
4楼: Originally posted by spur at 2012-05-08 16:02:26:
谢谢你的回复。那能不能说一下非对角元的意义呢？另外，Neaton的文章是实验方面的吧？实验上各种因素都可能影响这个极化值。但理论计算如果差距那么大，貌似不太合理。

在高对称结构中，沿高对称性方向移动的原子对电子能带的影响也具有比较高的对称性。比如立方BaTiO3中Ti沿[001]分量方向产生位移，根据对称性O 2p电子是不会产生[100]和[010]分量方向上的极化，所以这部分非对角元为0。如果是低对称性结构，那么原子位移在其它分量上也将会产生微弱的电子极化，不过一般来说非对角元BEC比对角元BEC要小不少。所以你的变化那么大有些奇怪
Neaton那篇文章是理论计算，他解释成如果选取的顺电相结构不同的话，那么极化值也不同。我不太相信他这个解释。
赞一下 回复此楼
8楼2012-05-08 16:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaomeiyan

木虫 (正式写手)
★
spur(金币+1): 谢谢参与
祝福楼主
回复此楼
13楼2012-05-08 16:52:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

db20021212

铁虫 (小有名气)
★
spur(金币+1): 谢谢参与
學習了
回复此楼
19楼2012-05-08 18:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
引用回帖:
8楼: Originally posted by iamikaruk at 2012-05-08 16:34:20:
在高对称结构中，沿高对称性方向移动的原子对电子能带的影响也具有比较高的对称性。比如立方BaTiO3中Ti沿分量方向产生位移，根据对称性O 2p电子是不会产生和分量方向上的极化，所以这部分非对角元为0。如果是低对 ...

谢谢你的讲解！不知道能不能请你帮忙看看我的计算结果或者另外一篇文章(首页附件)？如果这个结果是由于计算方法的问题，那我可能得用berry phase做。但如果用berry phase，怎么计算BEC的非对角元呢？谢谢！
赞一下 回复此楼
20楼2012-05-08 20:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
我大概看了一下你的计算结果，BEC没有太大问题，倒是你的位移量是怎么来的呢？保持晶格常数不变只弛豫内部离子坐标吗？
赞一下 回复此楼
22楼2012-05-08 21:12:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)
★
spur(金币+1): 谢谢参与
皆张量~~~~~~~~
回复此楼
25楼2012-05-08 22:53:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啄虫鸟

新虫 (著名写手)
★
spur(金币+1): 谢谢参与
楼主大方！
回复此楼
29楼2012-05-08 23:49:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

捡金小儿

新虫 (著名写手)
★
spur(金币+1): 谢谢参与
楼主大方！
回复此楼
30楼2012-05-08 23:53:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拾金昧者

木虫 (著名写手)
★
spur(金币+1): 谢谢参与
谢谢楼主！
回复此楼
31楼2012-05-08 23:57:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
贺仪: , 呵呵 以后设置为手动发放金币吧 2012-05-09 14:32:12
太扯了！楼上这么多来捡金币的。
赞一下 回复此楼
36楼2012-05-09 07:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangzc

金虫 (正式写手)
★
spur(金币+1): 谢谢参与
受教啦
回复此楼
39楼2012-05-09 08:23:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangy1220

新虫 (职业作家)
★
spur(金币+1): 谢谢参与
祝福楼主
回复此楼
40楼2012-05-09 08:26:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangy1220

新虫 (职业作家)
祝福楼主
回复此楼
42楼2012-05-09 08:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
引用回帖:
48楼: Originally posted by iamikaruk at 2012-05-09 15:11:58:
大概算了一下，原子太多，好麻烦
不过我大概猜测一下，你给出的Born effective charge应该是笛卡尔坐标下的，而位移用的不是笛卡尔坐标，会不会你这里没有考虑这个问题？

你好。我的位移就是采用的笛卡尔坐标，见首页附件2。
我还确实有点担心BEC究竟是什么坐标。我后来也验算了一下，如果BEC是内坐标，那位移肯定是需要用内坐标乘以晶格常数了，但得到的极化值更小，只有40多。
赞一下 回复此楼
49楼2012-05-09 15:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
49楼: Originally posted by spur at 2012-05-09 15:41:10:
你好。我的位移就是采用的笛卡尔坐标，见首页附件2。
我还确实有点担心BEC究竟是什么坐标。我后来也验算了一下，如果BEC是内坐标，那位移肯定是需要用内坐标乘以晶格常数了，但得到的极化值更小，只有40多。

我用的是六方坐标，这样的话x-y平面上的极化全部相互抵消，只剩下沿z轴的极化。
算出来的Bi沿z轴位移为-0.00088，Fe沿z轴位移为-0.02904，O沿z轴位移为-0.04436，算出来才40~50 \mu C/cm^2而已
赞一下 回复此楼
50楼2012-05-09 18:45:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
人间风月6楼
2012-05-08 16:17   回复  
spur(金币+1): 谢谢参与
jackycj7楼
2012-05-08 16:20   回复  
spur(金币+1): 谢谢参与
jybai9楼
2012-05-08 16:34   回复  
spur(金币+1): 谢谢参与
xdkevin10楼
2012-05-08 16:39   回复  
spur(金币+1): 谢谢参与
yangli89021311楼
2012-05-08 16:42   回复  
spur(金币+1): 谢谢参与
zhang731112楼
2012-05-08 16:48   回复  
spur(金币+1): 谢谢参与
祝福
tangbohejin14楼
2012-05-08 17:07   回复  
spur(金币+1): 谢谢参与
nanmann15楼
2012-05-08 17:45   回复  
spur(金币+1): 谢谢参与
qingheyouyou16楼
2012-05-08 17:48   回复  
spur(金币+1): 谢谢参与
mooncakeaa17楼
2012-05-08 17:48   回复  
spur(金币+1): 谢谢参与
jackycj18楼
2012-05-08 18:01   回复  
86739294924楼
2012-05-08 22:50   回复  
spur(金币+1): 谢谢参与
肖函26楼
2012-05-08 22:57   回复  
spur(金币+1): 谢谢参与
Easony28楼
2012-05-08 23:19   回复  
spur(金币+1): 谢谢参与
Scifinder232楼
2012-05-09 02:13   回复  
spur(金币+1): 谢谢参与
fengmeiyi33楼
2012-05-09 02:22   回复  
spur(金币+1): 谢谢参与
jdh92534楼
2012-05-09 07:11   回复  
spur(金币+1): 谢谢参与
紫芳枝叶35楼
2012-05-09 07:29   回复  
spur(金币+1): 谢谢参与
应用人生37楼
2012-05-09 08:03   回复  
spur(金币+1): 谢谢参与
icebergwu41楼
2012-05-09 08:35   回复  
spur(金币+1): 谢谢参与
105647872643楼
2012-05-09 08:50   回复  
spur(金币+1): 谢谢参与
v [ 发自手机版 http://muchong.com/3g ]
相关版块跳转 我要订阅楼主 spur 的主题更新
591/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 v7Aev8ORQjin 2026-08-22 3/150 2026-08-23 20:58 by 0p9y1isJp2jM
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 19:55 by lXgU0RvNTDC9
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:43 by lXgU0RvNTDC9
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:34 by lXgU0RvNTDC9
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:22 by lXgU0RvNTDC9
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 7/350 2026-08-23 19:10 by lXgU0RvNTDC9
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 18:55 by lXgU0RvNTDC9
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭