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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助DFTP算出来的压电系数
81/1返回列表
查看: 2989  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jk射手座

新虫 (著名写手)

[求助] 求助DFTP算出来的压电系数 已有1人参与

我用的VASP5.2 下的LEPSILON-tag算的压电系数和BEC等东西，之前发过类似的贴，不过这次想一起贴上来在好好求助一下：
我的INCAR
SYSTEM=PLPK
ENCUT=400
PREC=Normal
ISMEAR=0
SIGMA=0.05
EDIFF=1E-3
ISTART=0
ICHARG=2
IBRION=8
POTIM=0.1
LEPSILON=.TRUE.
LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 -1 3 -1 -1 -1 -1 -1 -1
LDAUU=0 0 6 0 0 0 0 0 0
LDAUJ=0 0 0 0 0 0 0 0 0
LDAUPRINT=2

算的是一个20个原子的四方相结构，含极性不连续的。主要算它的极化强度和压电特性在应变条件下的变化，现在只是OUTCAR中的结果还差那么点儿分析，想多请教下算过相关内容的前辈和同行，希望不吝赐教。
比如0应变的结构，我从OUTCAR中找到的有用内容如下：（目前提取出来了弹性模量，BEC和压电张量矩阵，还有个介电张量）

MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects in DFT)
---------------------------------------------
           7.562     0.000     0.000
           0.000     7.562     0.000
           0.000     0.000     6.072
---------------------------------------------


PIEZOELECTRIC TENSOR  for field in x, y, z        (e  Angst)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x    -0.57655     0.53337     1.41549     0.00000     0.00000    -2.48765
  y     0.74595    -0.80028     0.02920     0.00000    -2.29471     0.00000
  z    -2.13003    -2.13003    -1.71919     0.00000     0.00000     0.00000

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.16001
  y     0.00000     0.00000     0.00000     0.00000    -0.16005     0.00000
  z    -0.14256    -0.14256    -0.11506     0.00000     0.00000     0.00000
BORN EFFECTIVE CHARGES (in e, cummulative output)
-------------------------------------------------
ion    1
    1     4.16337     0.00000    -0.35555
    2     0.00000     4.25218    -0.00273
    3     0.00000     0.00000     3.08453
ion    2
    1     5.65192     0.00000    -1.45143
    2     0.00000     5.62322     0.02908
    3     0.00000     0.00000     6.09140
ion    3
    1     4.59437     0.00000    -0.77293
    2     0.00000     4.57861     0.02654
    3     0.00000     0.00000     4.69899
ion    4
    1     2.97250     0.00000    -0.46660
    2     0.00000     2.96029     0.01771
    3     0.00000     0.00000     4.07500
ion    5
    1     4.17827     0.00000    -0.05032
    2     0.00000     4.15128     0.00691
    3     0.00000     0.00000     3.68609
ion    6
    1     6.14756     0.00000     0.98037
    2     0.00000     6.12975    -0.01380
    3     0.00000     0.00000     4.99445
ion    7
    1     1.18633     0.00000     0.29818
    2     0.00000     1.18496    -0.00823
    3     0.00000     0.00000     1.34329
ion    8
    1     8.15684     0.00000     0.91563
    2     0.00000     7.88065    -0.02060
    3     0.00000     0.00000     6.62523
ion    9
    1    -2.67478     0.00000     0.42944
    2     0.00000    -4.21338     0.03375
    3     0.00000     0.00000    -2.01709
ion   10
    1    -4.10024     0.00000     0.31148
    2     0.00000    -2.70871    -0.04404
    3     0.00000     0.00000    -2.01709
ion   11
    1    -3.29422     0.00000     0.52821
    2     0.00000    -3.19993    -0.00621
    3     0.00000     0.00000    -4.46866
ion   12
    1    -3.03210     0.00000     0.35551
    2     0.00000    -2.55061     0.00564
    3     0.00000     0.00000    -3.02006
ion   13
    1    -2.56321     0.00000     0.32517
    2     0.00000    -2.93166    -0.02865
    3     0.00000     0.00000    -3.02006
ion   14
    1    -2.59899     0.00000     0.66267
    2     0.00000    -2.64079    -0.01770
    3     0.00000     0.00000    -3.98642
ion   15
    1    -2.45927     0.00000    -0.27148
    2     0.00000    -4.36881     0.03943
    3     0.00000     0.00000    -1.74718
ion   16
    1    -4.39566     0.00000    -0.19074
    2     0.00000    -2.44264    -0.03256
    3     0.00000     0.00000    -1.74718
ion   17
    1    -3.00654     0.00000    -0.27400
    2     0.00000    -3.01064    -0.00666
    3     0.00000     0.00000    -3.64142
ion   18
    1    -2.20096     0.00000     0.01679
    2     0.00000    -5.89128     0.02448
    3     0.00000     0.00000    -1.75444
ion   19
    1    -6.05680     0.00000    -0.01867
    2     0.00000    -2.16964    -0.03114
    3     0.00000     0.00000    -1.75444
ion   20
    1    -1.57589     0.00000    -0.92688
    2     0.00000    -1.56198     0.02343
    3     0.00000     0.00000    -4.92454
另外还有个离子贡献的介电张量：
MACROSCOPIC STATIC DIELECTRIC TENSOR IONIC CONTRIBUTION
---------------------------------------------
           0.155    -0.004    -0.114
          -0.004     0.000     0.001
          -0.114     0.001     7.979
---------------------------------------------
那么问题来了：
1、首先是两个压电张量，为什么除了单位不一样，它们输出的格式也不一样，后者好多项是0，请问那是eij 的矩阵吗？可不可以换算成dij 的矩阵？
2、离子的介电张量贡献有用吗？或者说这部分贡献分析什么东西有用？
目前就这俩问题，请大神们帮看看吧！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
主做电弧增材制造方向
1楼 2014-11-27 15:33:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
自顶一下，顺便说下，又试了算了下，加不加LRPA=.FALSE. ，提取的介电矩阵和压电系数的结果是一样 的。。。这正常吗？
赞一下 回复此楼
主做电弧增材制造方向
2楼2014-11-29 22:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
晕，不要沉啊
回复此楼
主做电弧增材制造方向
3楼2014-12-02 09:24:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
拜托，难道没人能帮说一下介电张量的问题吗？
赞一下 回复此楼
主做电弧增材制造方向
4楼2014-12-04 20:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主
★ ★
ljw4010: 金币+2, 谢谢指导解答！ 2014-12-05 20:16:51
引用回帖:
2楼: Originally posted by jk射手座 at 2014-11-29 22:02:02
自顶一下，顺便说下，又试了算了下，加不加LRPA=.FALSE. ，提取的介电矩阵和压电系数的结果是一样 的。。。这正常吗？

因为LRPA的默认值就是FALSE
赞一下 回复此楼
5楼2014-12-05 19:50:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
5楼: Originally posted by fzx2008 at 2014-12-05 19:50:57
因为LRPA的默认值就是FALSE...

谢谢回复。请问您算过介电常数吗，那个介电张量矩阵没太搞懂怎么处理。。。
赞一下 回复此楼
主做电弧增材制造方向
6楼2014-12-06 18:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

obaica

禁言 (著名写手)
本帖内容被屏蔽
7楼2016-10-27 15:32:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小蓝天空

新虫 (正式写手)

【答案】应助回帖

楼主你好，我尝试用过这个参数 这个好像是密度泛函线性响应的方法来计算压电矩阵 vasp可直接计算 只用于三维体系，我也是初学者 刚接触 说的不对的地方希望多指正多讨论。
赞一下 回复此楼
8楼2017-10-23 11:01:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了，能看出什么问题？ +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息，中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
关闭
请填处理意见
关闭