24小时热门版块排行榜    

查看: 4833  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

junfeiyan

新虫 (小有名气)

[求助] 求助二维材料如MoS2用berry phase 方法计算极化强度, 已有1人参与

如题,本身一直从事铁电材料的计算,最近看到一些文献,对关于二维材料的压电效应(应力/应变诱导极化)比较感兴趣,想算算。我对块体材料的极化计算比较擅长。

对于二维材料,如MoS2 在应变作用下,应该会有极化,那么该怎么计算呢? 对于块体材料极化可定义为单位体积的偶极矩,这个很确定。
那么二维材料的极化呢, 单位体积的偶极矩? 体积怎么确定,a, b , c 之积吗, 还是a, b, 乘以上下两个S原子的距离, 感觉都不靠谱。
看了一篇PRB 算MoS2 压电系数的文章, 极化定义为:However, in the case of a 2D system it is more meaningful to define a 2D polarization per unit area, P2D
这样看来,应该是偶极矩除以表面积,得到极化强度, 小弟试了一下,跟文献值差了个2倍的关系(也许是巧合)。所以还是不确定,该如何利用berry phase 的结果得到MoS2的极化强度。

另外本人跑了DFPT 计算了波恩有效电荷(两个Mo, 4个S), 想利用有效电荷法计算,极化强度,一个很奇怪的现象出现了, Mo的有效电荷是负的,S 的有效电荷是正的。有效电荷如下
BORN EFFECTIVE CHARGES (in e, cummulative output)
-------------------------------------------------
ion    1
    1     0.05108     0.00000     0.00000
    2     0.00007     0.53235     0.00000
    3     0.00000     0.00001     0.53261
ion    2
    1     0.05108     0.00001     0.00000
    2     0.00000     0.53259     0.00000
    3     0.00000     0.00000     0.53236
ion    3
    1     0.05109     0.00000     0.00000
    2     0.00005     0.53241     0.00000
    3     0.00000    -0.00001     0.53263
ion    4
    1     0.05109     0.00001     0.00000
    2     0.00001     0.53265     0.00000
    3     0.00000    -0.00001     0.53239
ion    5
    1    -0.10264    -0.00001     0.00000
    2    -0.00006    -1.06469     0.00000
    3     0.00000     0.00000    -1.06461
ion    6
    1    -0.10263    -0.00001     0.00000
    2    -0.00007    -1.06455     0.00000
    3     0.00000     0.00000    -1.06463
压电系数如下
PIEZOELECTRIC TENSOR  for field in x, y, z        (e  Angst)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x    -0.00001    -0.00008     0.00008     0.00003     0.00000     0.00000
  y    -0.00006    -3.33781     3.33756     0.00067     0.00000     0.00000
  z     0.00001     0.00000     0.00001     0.00000     3.33792    -0.00067

PIEZOELECTRIC TENSOR  for field in x, y, z        (C/m^2)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  y     0.00000    -0.19052     0.19051     0.00004     0.00000     0.00000
  z     0.00000     0.00000     0.00000     0.00000     0.19053    -0.00004

PIEZOELECTRIC TENSOR IONIC CONTR  for field in x, y, z        (C/m^2)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  y    -0.00060    -0.03825     0.03881    -0.00004     0.00000     0.00000
  z     0.00000     0.00000     0.00000     0.00000     0.04669     0.00005

小弟不是做二维材料的,对其性质不是特别了解,纯粹兴趣使然,想看看二维材料的极化怎么算, 求助各位大神,小弟就剩这么多金币了,在此谢过!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

physicsshh

新虫 (初入文坛)

小弟最近在学习DFPT计算波恩有效电荷。 OUTCAR里面给出的是每个原子的伯恩有效电荷张量,是个3*3矩阵,可以指导一下需要如何处理来得到每个原子的波恩有效电荷吗?
8楼2020-10-05 17:33:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

iamikaruk

木虫 (著名写手)

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。。 2016-04-13 02:07:00
对于有真空的体系,用DFPT可能不太合适?
感觉文献中对二维材料的极化定义都不妥
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
3楼2016-04-10 22:21:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蜀道难

铁虫 (小有名气)

请问dfpt算压电时初始结构是超胞还是周期原胞呢

发自小木虫IOS客户端
4楼2016-12-21 00:39:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

small_grass

银虫 (小有名气)

引用回帖:
2楼: Originally posted by junfeiyan at 2016-04-10 21:28:16
自己顶一下

你好,你会算二维材料的极化么?请教一下
温柔要有,但不妥协,要在安静中,不慌不忙的刚强
5楼2018-01-16 21:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:37 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:34 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:23 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 19:07 by xiangy672
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 G6APbkg8SA6w 2026-08-29 6/300 2026-08-30 19:02 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:07 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 17:34 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +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
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见