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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问DFPT能算出来声子吗
111/2返回列表12下一页
查看: 2886  |  回复: 31
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jk射手座

新虫 (著名写手)

[求助] 请问DFPT能算出来声子吗已有3人参与

如题，密度泛函微扰理论理论能算出来压电超晶格的声子频率吗？？对声子的计算不实很了解，只知道大多是拿ABINIT算压电材料的声子。不过今天重翻某项目书，发现居然写着“原子有效电荷、压电系数以及声子频率可以通过密度泛函微扰理论进行计算”，  BEC 和压电系数张量我知道能算，OUTCAR中我都能找到，但是声子频率我没找到呀，我是VASP5.2的，请问INCAR中是不是还要加什么东西才能算声子呢？比较急，求路过的大神同僚们帮忙解答下，如果不能算也请留个肯定的答复！谢谢！！！！

PS：另外想问一个压电张量结果的问题：
PIEZOELECTRIC TENSOR  for field in x, y, z        (e  Angst)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x    -0.62281     0.65564    -0.03830     0.00000     0.00000     0.84200
  y     0.70177    -0.70921     0.00903     0.00000     0.81327     0.00000
  z     0.19100     0.19100    -0.57731     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.11620
  y     0.00000     0.00000     0.00000     0.00000     0.11659     0.00000
  z     0.02687     0.02687    -0.08122     0.00000     0.00000     0.00000

这两个在OUTCAR中挨着的结果，除了单位不一样，各自代表什么意思呢？应该取哪个结果来分析呢？
回复此楼

» 猜你喜欢

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

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

jk射手座

新虫 (著名写手)
引用回帖:
2楼: Originally posted by KalaShayminS at 2014-11-25 17:33:39
声子的计算量非常大，需要单独做，手册里有。
后面两个是单位不同。

额。。好吧，对于声子还需酌情考虑。

那么这两组结果好像不太一样啊，请问一般应该取哪个单位的呢
为啥后者的XX YY为零了呢。。。
赞一下 回复此楼
主做电弧增材制造方向
4楼2014-11-25 20:22:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
内容已删除
回复此楼
主做电弧增材制造方向
5楼2014-11-25 20:24:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
6楼: Originally posted by future_wl at 2014-11-26 08:48:33
VASP只能给出GAMMA点的声子振动频率，要计算整个布里渊区的声子谱需要外接程序如phonopy或phonon等。如果用PWSCF软件则可以直接计算声子谱！

哦这样啊，谢谢你的回复！请问phonopy或phonon还有PWSCF软件哪个比较i容易使用一点呢，如果很费时间我只好舍掉了。
赞一下 回复此楼
主做电弧增材制造方向
8楼2014-11-27 09:08:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
6楼: Originally posted by future_wl at 2014-11-26 08:48:33
VASP只能给出GAMMA点的声子振动频率，要计算整个布里渊区的声子谱需要外接程序如phonopy或phonon等。如果用PWSCF软件则可以直接计算声子谱！

对了，请问你计算过压电应变常数d33吗
赞一下 回复此楼
主做电弧增材制造方向
9楼2014-11-27 14:55:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
13楼: Originally posted by future_wl at 2014-11-27 20:28:12
压电常数用DFTP就可以算啊，手册上有教程....

计算出来的张量不是 压电应力常数eij 的吗？我是想得到 压电应变常数dij的，是不是需要公式换算的呢？另外，您说的手册是英文版的还是哪个版本的？不知道我看的是否和你看的一样，我在看一下，谢谢
赞一下 回复此楼
主做电弧增材制造方向
16楼2014-11-28 09:21:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
21楼: Originally posted by future_wl at 2014-11-29 15:42:42
IBRION=5/6  不是冻声子法，是有限差分的方法
冻声子方法是计算声子谱的，VASP自身无法计算声子谱...

请问对于DFPT下 IBRION=8计算出的压电矩阵和介电矩阵，您了解吗，可否看下我近期求助的一个问题帖子
赞一下 回复此楼
主做电弧增材制造方向
22楼2014-11-29 15:48:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
17楼: Originally posted by future_wl at 2014-11-28 14:51:49
好像需要乘以弹性柔顺系数就可以转成实验的压电应变常数了，具体公式你搜一下，我忘记了...

你好，我都到一个讲压电材料的PPT，这是里面截的两页图，里面好像有两种d的表达式，其中有你说的那个柔顺系数。  另外这个柔顺系数能算吗？还是特定的材料柔顺系数是固定的？
请问DFPT能算出来声子吗
2.png


请问DFPT能算出来声子吗-1
1.png
赞一下 回复此楼
主做电弧增材制造方向
23楼2014-11-30 08:47:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
25楼: Originally posted by future_wl at 2014-12-02 09:40:24
你好，柔顺系数就是弹性常数的倒数，你可以计算弹性矩阵，然后求逆，就得到了柔顺系数的矩阵。...

多谢回复。我找了下我的结果中的数据：
TOTAL ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX        3403.2062    481.0583   2439.5217      0.0000      0.0000      0.3415
YY         481.0583   3566.4192   2139.4804      0.0000      0.0000      0.0574
ZZ        2439.5217   2139.4804 -12184.7143      0.0000      0.0000     -4.1046
XY           0.0000      0.0000      0.0000   1089.6760      0.0000      0.0000
YZ           0.0000      0.0000      0.0000      0.0000   1056.7884      0.0000
ZX           0.3415      0.0574     -4.1046      0.0000      0.0000   1056.7830
--------------------------------------------------------------------------------
这个应该是总的，还有离子贡献的弹性模量，请问是对哪个矩阵取逆矩阵呢？
SYMMETRIZED ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX        4387.6503   1222.8740   1195.7458      0.0000      0.0000      0.0000
YY        1222.8740   4387.6503   1195.7458      0.0000      0.0000      0.0000
ZZ        1195.7458   1195.7458   4488.2980      0.0000      0.0000      0.0000
XY           0.0000      0.0000      0.0000   1089.6760      0.0000      0.0000
YZ           0.0000      0.0000      0.0000      0.0000   1056.7884      0.0000
ZX           0.0000      0.0000      0.0000      0.0000      0.0000   1056.7885
--------------------------------------------------------------------------------

ELASTIC MODULI CONTR FROM IONIC RELAXATION (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX        -984.4441   -741.8156   1243.7758      0.0000      0.0000      0.3415
YY        -741.8156   -821.2311    943.7346      0.0000      0.0000      0.0574
ZZ        1243.7758    943.7346 -16673.0123      0.0000      0.0000     -4.1046
XY           0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
YZ           0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
ZX           0.3415      0.0574     -4.1046      0.0000      0.0000     -0.0055
赞一下 回复此楼
主做电弧增材制造方向
26楼2014-12-02 12:18:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
27楼: Originally posted by future_wl at 2014-12-02 16:19:56
取第一个，总的...

多谢！
回复此楼
主做电弧增材制造方向
28楼2014-12-03 08:37:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭