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

24小时热门版块排行榜    

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

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的回帖
查看全部 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 重要来源：本周末出结果 +9 瞬息宇宙 2026-08-12 9/450 2026-08-13 10:12 by Vivilian
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
关闭
请填处理意见
关闭