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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？
北京石油化工学院2026年研究生招生接收调剂公告
181/2返回列表12下一页
查看: 2547  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

donggaomin

金虫 (正式写手)

[求助] 请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？

请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？

说明书上只说了设置IBRION = 6

另外，怎么查看结果？论坛里提到是动力学矩阵的对角化，不是很明白。我知道动力学矩阵对角化后是频率的平方。

谢谢！

VASP5.2中是不是可以计算phonon frequency？
回复此楼

» 收录本帖的淘帖专辑推荐

声子谱

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
好好学习
1楼 2011-05-28 11:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
补充：这是在VASP论坛中找到的：
you can calculate the gamma-point phonons with vasp.4.6, using a frozen phonon approach
please use IBRION=5 in that case, and set POTIM = 0.015 (as the harmonic limit is assumed)
in case you do not need all modes, please use the Selective dymanics key in POSCAR and just allow the vibrations you are intrested in (by setting F/T ) in order to save CPU time. it is not possible to continue a broken phonon job (keeping the info) it has to be restarted from the beginning.
relativistic effects can be simulated by using spin-orbit splitting (please have a look at the LSORBIT tag, and the additional parameters you have to set.

这是第二个
To do \Gamma point phonon calculate

Firstly, fully relax your structure and save WAVECAR

Secondly, you can start a phonon calculation by set INCAR as

PREC= Accurate
ADDGRID = .TRUE.
ENCUT= 500.00 eV
ISTART=1
ICHARG=0
IBRION= 6
NFREE=4
POTIM=0.015 !displacement of atoms in Ang.
IALGO= 48
LREAL=Auto !needed for good convergence
NSIM = 4
LCHARG = .False.
LPLANE=.TRUE.
~


现在的问题是：是不是第二个的方式好点？如何查看结果？
赞一下 回复此楼
好好学习
2楼2011-05-28 16:27:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
请计算过Gamma点声子频率的高手出招！
赞一下 回复此楼
好好学习
3楼2011-05-28 20:42:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
下面是设置了IBRION = 6时的动力学矩阵本征值，按我的理解是本征值就是频率的平方，但...    请问怎么看下面的结果呢？

Eigenvectors and eigenvalues of the dynamical matrix
----------------------------------------------------


   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    -0.000519   -0.000028    0.081397  
      1.609613  0.929267  2.583248    -0.000519   -0.000028   -0.081397  
     -0.000113  1.858561  1.976701     0.001581    0.000057   -0.702404  
      1.609613  0.929267  4.559948     0.001581    0.000057    0.702404  

   2 f  =   13.134427 THz    82.526040 2PiTHz  438.117317 cm-1    54.319644 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000    -0.092238   -0.254007    0.000249  
      1.609613  0.929267  2.583248     0.092238    0.254007    0.000249  
     -0.000113  1.858561  1.976701    -0.222415   -0.614413   -0.000695  
      1.609613  0.929267  4.559948     0.222415    0.614413   -0.000695  

   3 f  =   13.116478 THz    82.413261 2PiTHz  437.518595 cm-1    54.245412 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.254559   -0.092022    0.000439  
      1.609613  0.929267  2.583248    -0.254559    0.092022    0.000439  
     -0.000113  1.858561  1.976701     0.614185   -0.222503   -0.000963  
      1.609613  0.929267  4.559948    -0.614185    0.222503   -0.000963  

   4 f  =   12.231166 THz    76.850683 2PiTHz  407.987772 cm-1    50.584055 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000    -0.128882   -0.287868   -0.000033  
      1.609613  0.929267  2.583248    -0.128882   -0.287868    0.000033  
     -0.000113  1.858561  1.976701     0.258437    0.577694    0.000731  
      1.609613  0.929267  4.559948     0.258437    0.577694   -0.000731  

   5 f  =   12.215895 THz    76.754733 2PiTHz  407.478389 cm-1    50.520900 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000    -0.288034    0.128850   -0.000113  
      1.609613  0.929267  2.583248    -0.288034    0.128850    0.000113  
     -0.000113  1.858561  1.976701     0.577610   -0.258454    0.001474  
      1.609613  0.929267  4.559948     0.577610   -0.258454   -0.001474  

   6 f  =   10.927284 THz    68.658148 2PiTHz  364.494939 cm-1    45.191629 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000    -0.000189    0.000636    0.313532  
      1.609613  0.929267  2.583248     0.000189   -0.000636    0.313532  
     -0.000113  1.858561  1.976701    -0.000728    0.000873   -0.633795  
      1.609613  0.929267  4.559948     0.000728   -0.000873   -0.633795  

   7 f  =    7.595291 THz    47.722620 2PiTHz  253.351625 cm-1    31.411609 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.000310    0.000022    0.702406  
      1.609613  0.929267  2.583248     0.000310    0.000022   -0.702406  
     -0.000113  1.858561  1.976701     0.000072   -0.000016    0.081397  
      1.609613  0.929267  4.559948     0.000072   -0.000016   -0.081397  

   8 f  =    2.947281 THz    18.518311 2PiTHz   98.310698 cm-1    12.188977 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.129670    0.640455   -0.000783  
      1.609613  0.929267  2.583248    -0.129670   -0.640455   -0.000783  
     -0.000113  1.858561  1.976701    -0.053713   -0.264796   -0.000086  
      1.609613  0.929267  4.559948     0.053713    0.264796   -0.000086  

   9 f  =    2.934466 THz    18.437793 2PiTHz   97.883244 cm-1    12.135980 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.640217   -0.129724    0.000119  
      1.609613  0.929267  2.583248    -0.640217    0.129724    0.000119  
     -0.000113  1.858561  1.976701    -0.265373    0.053582    0.000117  
      1.609613  0.929267  4.559948     0.265373   -0.053582    0.000117  

  10 f/i=    0.003591 THz     0.022565 2PiTHz    0.119792 cm-1     0.014852 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.000007   -0.000665   -0.633795  
      1.609613  0.929267  2.583248    -0.000007    0.000665   -0.633795  
     -0.000113  1.858561  1.976701    -0.000005    0.000373   -0.313533  
      1.609613  0.929267  4.559948     0.000005   -0.000373   -0.313533  

  11 f/i=    0.036145 THz     0.227107 2PiTHz    1.205674 cm-1     0.149485 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.063064    0.629725   -0.000046  
      1.609613  0.929267  2.583248     0.063064    0.629725    0.000046  
     -0.000113  1.858561  1.976701     0.031392    0.313821    0.000019  
      1.609613  0.929267  4.559948     0.031392    0.313821   -0.000019  

  12 f/i=    0.039877 THz     0.250556 2PiTHz    1.330159 cm-1     0.164919 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000     0.629641   -0.063053   -0.000333  
      1.609613  0.929267  2.583248     0.629641   -0.063053    0.000333  
     -0.000113  1.858561  1.976701     0.313989   -0.031415    0.000203  
      1.609613  0.929267  4.559948     0.313989   -0.031415   -0.000203
赞一下 回复此楼
好好学习
4楼2011-05-29 17:04:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★
donggaomin(金币+5): 谢谢您！ 2011-05-29 21:15:41
zzy870720z(金币+2): 谢谢指教 2011-05-29 23:32:08
NSW=1
IBRION=5
NFREE=2
POTIM=0.05 POTIM 的值取 0.04-0.1合理。
这是采用有限差分法计算的
赞一下(2人) 回复此楼
5楼2011-05-29 19:55:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★
donggaomin(金币+5): 谢谢！请问1f,2f 3f 4f 5f ... f/i 等分别表示什么？如何判断振动类型？ 2011-05-29 21:21:12
zzy870720z(金币+2): 谢谢指教 2011-05-29 23:32:18
引用回帖:
Originally posted by donggaomin at 2011-05-29 17:04:52:
下面是设置了IBRION = 6时的动力学矩阵本征值，按我的理解是本征值就是频率的平方，但...    请问怎么看下面的结果呢？

Eigenvectors and eigenvalues of the dynamical matrix
--------------------------- ...

举例
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    -0.000519   -0.000028    0.081397  
      1.609613  0.929267  2.583248    -0.000519   -0.000028   -0.081397  
     -0.000113  1.858561  1.976701     0.001581    0.000057   -0.702404  
      1.609613  0.929267  4.559948     0.001581    0.000057    0.702404  

第一行就是频率，不同的单位给出。
后面X Y Z是原子坐标
在后面的 dx dy dz是动力学矩阵的本征矢量，判断振动类型的时候用到。
赞一下(2人) 回复此楼
6楼2011-05-29 20:00:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★
donggaomin(金币+5): 谢谢dxcharlary!能不能举例说明一下关于强度和稳定性的计算？ 2011-05-29 21:26:56
zzy870720z(金币+2): 谢谢指教 2011-05-29 23:32:27
引用回帖:
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    ...

说的很好，补充几点：
5用的是差分，6用的是DFPT，感觉好像DFPT更靠谱一点。不过貌似VASP的计算结果没有PWSCF的好，而且只能算gamma点的。
dx,dy,dz是本征矢，也是震动方向，有时做一些强度或者是稳定性计算的时候可能会用到。
赞一下(1人) 回复此楼
7楼2011-05-29 20:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★
donggaomin(金币+5): 2011-05-29 21:34:28
zzy870720z(金币+1): 谢谢补充 2011-05-29 23:32:41
对了，在这些东西的上面会有你要的动力学矩阵。
赞一下(1人) 回复此楼
8楼2011-05-29 20:23:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★ ★
donggaomin(金币+5): 谢谢stractor!学习了！ 2011-05-29 21:30:19
zzy870720z(金币+3): 谢谢指教 2011-05-29 23:32:50
是的，有限差分的可信度取决于原子的位移
因为有限差分计算Hessen矩阵是
Hij=d2E/dxidxj=[E(ri0+dxi,rj0+dxj)+E(ri0-dxi,rj0-dxj)-2E(ri0,rj0)]/dxidxj
上面的d是微分，rjo就是j原子的平衡位置，E是能量。

上式反映出位移xi, xj不能太大。
太小的话位移后构型的能量和没位移的构型的能量没差别，也会导致误差。
还有一个误差，有些体系会遇到，就是VASP中计算频率是简谐近似，忽略了能量对位移泰勒展开的高阶项。
我测试过H2,CO分子，合理选择位移的话，计算值和实验值差别在10cm-1以内。
赞一下(1人) 回复此楼
9楼2011-05-29 21:18:00
已阅   回复此楼   关注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    ...

请问如何分类？如下图中的E（1TO），A1（1TO），F(1TO)是如何从数据中判断？

谢谢！
赞一下 回复此楼
好好学习
10楼2011-05-29 21:41:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 donggaomin 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +7 倘若起风了呢 2026-04-05 7/350 2026-04-05 20:15 by 南航~万老师
[考研] 285求调剂 +3 恶法大二的气味� 2026-04-05 4/200 2026-04-05 19:59 by 啵啵啵0119
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +16 披星河 2026-04-03 16/800 2026-04-05 09:00 by dick_runner
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 312求调剂 +4 赊月色 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[考研] 土木304求调剂 +4 兔突突突， 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突，
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 285求调剂 +14 AZMK 2026-04-02 14/700 2026-04-02 15:54 by 上九天揽月（好�
[考研] 材料专业求调剂 +10 月月鸟木 2026-04-01 10/500 2026-04-02 12:57 by wxiongid
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[考研] 085600，320分求调剂 +5 大馋小子 2026-04-01 6/300 2026-04-01 19:40 by 唐沐儿
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
信息提示
关闭
请填处理意见
关闭