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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？
52/1返回列表
查看: 2909  |  回复: 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的回帖

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的回帖
查看全部 18 个回答

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的回帖
查看全部 18 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 面上函评意见出来了，像什么等级？ 20+3 Tsingking1 2026-08-27 14/700 2026-09-01 08:29 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭