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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问DFPT能算出来声子吗
101/1返回列表
查看: 2880  |  回复: 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的回帖

future_wl

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jk射手座(ljw4010代发): 金币+2, 谢谢交流！ 2014-11-26 19:31:51
jk射手座: 金币+8, ★★★很有帮助 2014-11-27 09:08:26
VASP只能给出GAMMA点的声子振动频率，要计算整个布里渊区的声子谱需要外接程序如phonopy或phonon等。如果用PWSCF软件则可以直接计算声子谱！
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

tjia tjia
未来就是现在
6楼2014-11-26 08:48:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
8楼: Originally posted by jk射手座 at 2014-11-27 09:08:14
哦这样啊，谢谢你的回复！请问phonopy或phonon还有PWSCF软件哪个比较i容易使用一点呢，如果很费时间我只好舍掉了。...

还是学一下phonopy吧，这个是可以和VASP结合的，并且开源，phonon要收费
赞一下 回复此楼
未来就是现在
11楼2014-11-27 20:23:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
10楼: Originally posted by zyl雨田 at 2014-11-27 16:01:13
你好，
这句话是我在phononpy manual 里面看到的:VASP can calculate force constants using DFPT however only at Gamma-point.
而按照你的观点，通过pgonopy接口，可以得到整个不里渊区的色散谱。
感觉你的观 ...

读入的是vasp计算的所有键的力常数，phonopy根据这些力常数的方向和大小来计算声子谱，如果你用phonopy或phonon软件完整计算一遍声子谱你就明白了!
赞一下 回复此楼
未来就是现在
12楼2014-11-27 20:27:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
9楼: Originally posted by jk射手座 at 2014-11-27 14:55:36
对了，请问你计算过压电应变常数d33吗...

压电常数用DFTP就可以算啊，手册上有教程.
赞一下 回复此楼
未来就是现在
13楼2014-11-27 20:28:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
16楼: Originally posted by jk射手座 at 2014-11-28 09:21:55
计算出来的张量不是 压电应力常数eij 的吗？我是想得到 压电应变常数dij的，是不是需要公式换算的呢？另外，您说的手册是英文版的还是哪个版本的？不知道我看的是否和你看的一样，我在看一下，谢谢...

好像需要乘以弹性柔顺系数就可以转成实验的压电应变常数了，具体公式你搜一下，我忘记了
赞一下 回复此楼
未来就是现在
17楼2014-11-28 14:51:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
15楼: Originally posted by zyl雨田 at 2014-11-28 09:20:35
我算了 所以才更疑惑。为什么手册里面会说得到的仅仅是G点的力常数。...

VASP只能算G点的，但通过phonopy的冻声子法，将原子位置沿不同方向进行一定幅度的移动之后再算力常数，就可以拟合出布里渊区其他对称点的声子振动频率。
赞一下 回复此楼
未来就是现在
18楼2014-11-28 15:18:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
20楼: Originally posted by zyl雨田 at 2014-11-28 21:12:45
因为按照 我的理解VASP 计算声子有两种方法，一种是冷冻声子法也就是直接位移法IBRION=5/6 另一种就是DFPT IBRION=7/8。你在这里提到 通过phonopy的冷冻声子法可以得到整个布里区的其他高对称点的频率。难以理解。...

IBRION=5/6  不是冻声子法，是有限差分的方法
冻声子方法是计算声子谱的，VASP自身无法计算声子谱
赞一下 回复此楼
未来就是现在
21楼2014-11-29 15:42:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
23楼: Originally posted by jk射手座 at 2014-11-30 08:47:36
你好，我都到一个讲压电材料的PPT，这是里面截的两页图，里面好像有两种d的表达式，其中有你说的那个柔顺系数。  另外这个柔顺系数能算吗？还是特定的材料柔顺系数是固定的？

2.png

1.png
...

你好，柔顺系数就是弹性常数的倒数，你可以计算弹性矩阵，然后求逆，就得到了柔顺系数的矩阵。
赞一下 回复此楼
未来就是现在
25楼2014-12-02 09:40:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
26楼: Originally posted by jk射手座 at 2014-12-02 12:18:20
多谢回复。我找了下我的结果中的数据：
TOTAL ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------- ...

取第一个，总的
回复此楼
未来就是现在
27楼2014-12-02 16:19:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭