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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问DFPT能算出来声子吗
323/4返回列表上一页1234下一页
查看: 2885  |  回复: 31
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的回帖

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的回帖

星星哥哥

金虫 (正式写手)
用免费的PWSCF软件计算更合适......
赞一下 回复此楼
心事浩茫连广宇,于无声处听惊雷.
24楼2014-11-30 09:55:58
 已阅   回复此楼   关注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的回帖

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的回帖

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射手座

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

多谢！
回复此楼
主做电弧增材制造方向
28楼2014-12-03 08:37:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjia

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

楼主，你好！我想问一下DFPT只能算Gamma点的力常数，我现在用phonopy能的到声子谱，那声子谱其他高对称点的数据是怎么来的？是Phonopy后续怎么处理的吗？主要是采用什么方法处理的？谢谢！
赞一下 回复此楼
29楼2015-10-21 11:02:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjia

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

楼主，你好！我想问一下DFPT只能算Gamma点的力常数，我现在用phonopy能的到声子谱，那声子谱其他高对称点的数据是怎么来的？是Phonopy后续怎么处理的吗？主要是采用什么方法处理的？谢谢！
赞一下 回复此楼
30楼2015-10-21 11:02:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
323/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭