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

24小时热门版块排行榜    

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

金虫 (正式写手)
下面是设置了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 个回答

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的回帖

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的回帖
查看全部 18 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[硕博家园] 极限拉扯，如何抉择 +5 Isak0320 2026-10-12 7/350 2026-10-12 18:03 by 有人要害朕
[教师之家] 课题组会帮你通过考核么？ +3 好运来～9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭