应《网络安全法》要求，自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的回帖

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的回帖
查看全部 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:22 by o3EMxHUAQDtS
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 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
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
关闭
请填处理意见
关闭