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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问DFPT能算出来声子吗
323/4返回列表上一页1234下一页
查看: 3411  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中，共勉 +5 eulota 2026-08-26 5/250 2026-08-26 23:33 by young90
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 国自然面上复盘~欢迎讨论 +4 晴天加油 2026-08-26 5/250 2026-08-26 21:27 by 木风9012
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
关闭
请填处理意见
关闭