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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？
182/2返回列表上一页12
查看: 2522  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

donggaomin

金虫 (正式写手)
引用回帖:
Originally posted by dxcharlary at 2011-05-29 20:21:14:
说的很好，补充几点：
5用的是差分，6用的是DFPT，感觉好像DFPT更靠谱一点。不过貌似VASP的计算结果没有PWSCF的好，而且只能算gamma点的。
dx,dy,dz是本征矢，也是震动方向，有时做一些强度或者是稳定性计算 ...

您好！

请问：1f,2f,3f...等是不是类似于能带图，1f表示其中的一条声频色散带在gamma点上的值？2f表示另一条声频色散带在gamma点上的值？...


我们在文章中常看到的“Gamma点A1(1TO)模式声子最低频率的平方“是不是这些”1f,2f,3f,...“中最小的一个值取平方得到？


四楼中给出的数据：最后的三个"10 f/i, 11f/i, 12f/i "中的i 好像是坐标轴，请问这三个是什么？

非常感谢！
赞一下 回复此楼
好好学习
11楼2011-05-29 21:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
引用回帖:
Originally posted by dxcharlary at 2011-05-29 20:23:02:
对了，在这些东西的上面会有你要的动力学矩阵。

我看了一下，这些东西上面是“SECOND DERIVATIVES (NOT SYMMETRIZED)”，（在IBRION = 6 时得到的）。


我现在的理解是“力常数矩阵就是HESSIAN矩阵，矩阵元是能量对位移的二阶倒数。动力学矩阵是力常数矩阵的傅里叶变换”

请问如何利用这个动力学矩阵？

第二，动力学矩阵的本征值应该是频率，而不是频率的平方是吗？

谢谢您的帮助！
赞一下 回复此楼
好好学习
12楼2011-05-29 22:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
引用回帖:
Originally posted by stractor at 2011-05-29 20:00:30:
举例
1 f  =   15.937321 THz   100.137142 2PiTHz  531.611795 cm-1    65.911486 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000    ...

请问stractor,
     E（1TO），A1（1TO），F(1TO)  我现在的理解是TO表示横波光学模声子，而1TO中的1表示能带一的，不知道我的理解对不对？

另外，E，A1，F 分别表示什么？从数据中怎么判断？


谢谢！
赞一下 回复此楼
好好学习
13楼2011-05-29 22:57:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

donggaomin(金币+3): 谢谢stractor! 2011-05-30 16:37:33
引用回帖:
Originally posted by donggaomin at 2011-05-29 21:53:15:
您好！

请问：1f,2f,3f...等是不是类似于能带图，1f表示其中的一条声频色散带在gamma点上的值？2f表示另一条声频色散带在gamma点上的值？...


我们在文章中常看到的“Gamma点A1(1TO)模式声子最低频率的 ...

1f就是gamma点的第一支频率，以此类推，f/i表示虚频，但数据太小，是计算误差造成的，真实只应该是0.
赞一下 回复此楼
14楼2011-05-30 04:36:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
引用回帖:
Originally posted by donggaomin at 2011-05-29 22:57:23:
请问stractor,
     E（1TO），A1（1TO），F(1TO)  我现在的理解是TO表示横波光学模声子，而1TO中的1表示能带一的，不知道我的理解对不对？

另外，E，A1，F 分别表示什么？从数据中怎么判断？


谢谢！

TO表示横波光学模声子,E, A1, F是按群操作对振动模式的分类。我不搞晶格动力学，对这个不懂。给你一个连接，你可以摸索一些东西出来：
http://www.cryst.ehu.es/rep/sam.html
赞一下 回复此楼
15楼2011-05-30 04:54:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

donggaomin(金币+2): 谢谢dxcharlary!!! 2011-05-30 16:39:59
引用回帖:
Originally posted by donggaomin at 2011-05-29 22:01:55:
我看了一下，这些东西上面是“SECOND DERIVATIVES (NOT SYMMETRIZED)”，（在IBRION = 6 时得到的）。


我现在的理解是“力常数矩阵就是HESSIAN矩阵，矩阵元是能量对位移的二阶倒数。动力学矩阵是力常数矩 ...

那个东西是hessian矩阵，k=0的时候，exp（-kr）=1.且单胞内的原子相互作用要远大于跨元胞之间的相互作用，所以动力学矩阵退化为hessian矩阵除以质量的平方根。
另外不好意思，告诉你错了。。。5，6是差分，7，8是DFPT，5，7是不考虑对称性，6，8是考虑对称性的。
赞一下 回复此楼
16楼2011-05-30 14:28:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)
关于那些模我一般都是看人家的文献的
赞一下 回复此楼
paper我要好多好的的paper
17楼2011-05-31 10:31:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ufo991

铁虫 (正式写手)
引用回帖:
475878楼: Originally posted by stractor at 2011-05-29 21:18:00
是的，有限差分的可信度取决于原子的位移
因为有限差分计算Hessen矩阵是
Hij=d2E/dxidxj=/dxidxj
上面的d是微分，rjo就是j原子的平衡位置，E是能量。

上式反映出位移xi, xj不能太大。
太小的话位移后构型的能 ...

你好，请问你哪里还有计算H2或者CO的输入文件嘛？能否给出一份？谢谢ufo991@163.com
赞一下 回复此楼
18楼2012-05-25 15:05:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 donggaomin 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 3/150 2026-03-15 17:32 by 小物理化学
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 085600材料与化工 求调剂 +11 enenenhui 2026-03-13 12/600 2026-03-14 21:48 by laoshidan
[考研] 化学工程321分求调剂（南京工业，浙江工业） +3 大米饭！ 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 293求调剂，一志愿陕师大生物学 +3 ??????.?.??? 2026-03-09 3/150 2026-03-11 10:02 by 学员8dgXkO
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
关闭
请填处理意见
关闭