| 查看: 2338 | 回复: 17 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
donggaomin金虫 (正式写手)
|
[求助]
请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？
|
||
|
请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？ 说明书上只说了设置IBRION = 6 另外，怎么查看结果？论坛里提到是动力学矩阵的对角化，不是很明白。我知道动力学矩阵对角化后是频率的平方。 谢谢！ VASP5.2中是不是可以计算phonon frequency？ |
回复此楼» 收录本帖的淘帖专辑推荐
 声子谱 |
» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有13人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
急问 现在哪里的计算中心有vasp5.2.12？
已经有6人回复- time-on-stream什么意思啊
已经有10人回复
- 全国第二届理论与计算化学软件Molpro Molcas培训班
已经有19人回复
- 很菜的问题：不设置参比电极不行么？
已经有9人回复
- 请教有机分子中O原子和H原子电荷数怎么计算，O-H键键级怎么计算！谢谢！
已经有4人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 计划书，word2007,宏设置问题
已经有5人回复
- ThinkPad SL410 2842EKC如何恢复出厂设置？
已经有9人回复
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 770.6
- 散金: 805
- 红花: 9
- 帖子: 564
- 在线: 242.7小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
下面是设置了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
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 770.6
- 散金: 805
- 红花: 9
- 帖子: 564
- 在线: 242.7小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
补充：这是在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
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 770.6
- 散金: 805
- 红花: 9
- 帖子: 564
- 在线: 242.7小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-05-28 20:42:54
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
5楼2011-05-29 19:55:42