| 8 | 1/1 | 返回列表 |
| 查看: 3466 | 回復: 7 | ||
cindywei木蟲 (初入文壇)
|
[求助]
VASP計算穩(wěn)定體系的零點能為什么也會出現(xiàn)虛頻?
|
|
最近用VASP計算Gamma點的零點能做校正,原本熱力學及力學上都能穩(wěn)定存在的鈀氫化合物有些成分點出現(xiàn)了虛頻,是哪里出現(xiàn)問題了嗎? INCAR設置 ISTART = 0 ICHARG = 2 ENCUT = 450 ISMEAR = 0 PREC = A EDIFF = 1E-6 EDIFFG = -0.001 NSW = 1 IBRION = 5 POTIM = 0.015 NFREE = 2 OUTCAR 1 f = 20.970164 THz 131.759427 2PiTHz 699.489360 cm-1 86.725659 meV 2 f/i= 15.743741 THz 98.920841 2PiTHz 525.154652 cm-1 65.110902 meV 3 f/i= 15.775374 THz 99.119597 2PiTHz 526.209815 cm-1 65.241726 meV |
木蟲 (正式寫手)
|
http://wolverton.northwestern.ed ... phonon-calculations One drawback of the frozen-phonon method is that large supercells are needed to accurately calculate the force constant matrix. The periodic boundary conditions used in DFT calculations cause problems for the frozen phonon method. Displacing one atom in a small unit cell creates forces on all the atoms in the same unit cell, but also on the periodic images of these atoms. |
木蟲 (小有名氣)
|
cindywei您好,我想問您一個問題就是我要做零點能校正的話,按E=1/2hv公式,把所有頻率的能量相加,大約就是0.7eV,按理說,零點能校正的話,不會是這么大的值,也就是0.02eV,是不是這中計算方法有誤呢? 我如果要對氣態(tài)分子吸附態(tài)進行零點能校正的話,比如H2O在Cu上的吸附,要計算頻率,我是不是應該把H2O和Cu的第一層原子弛豫,把所有的頻率的能量相加,就算是零點能呢?是不是應該把Cu原子弛豫的頻率去掉,只加H2O的頻率呢,但是那樣的話,我根本就不能從所有的頻率中區(qū)分H2O和Cu的頻率? 最近甚是苦惱啊,還望您能指點迷津,期待您的回復,在此,萬分感謝。。。 |

木蟲 (初入文壇)
|
計算零點能的方法有很多種吧,如果你只考慮輕質的H及O的零點能的話,可以在POSCAR選擇H及O不固定而Cu固定,這種計算不完全等同于晶格優(yōu)化的弛豫,是每個不固定原子按特定幾個方向有限的弛豫來獲得其零點振動能。我跟個INCAR的設置,你可以通過參數(shù)來理解一下,希望對你有用 ISTART = 0 ICHARG = 2 ENCUT = 375 ISMEAR = 0 PREC = A EDIFF = 1E-7 EDIFFG = -0.001 NSW = 1 IBRION = 5 POTIM = 0.015 NFREE = 2 |
木蟲 (小有名氣)

木蟲 (初入文壇)
木蟲 (小有名氣)

木蟲 (初入文壇)
|
是先把整個結構弛豫優(yōu)化獲得最優(yōu)化結構,然后對優(yōu)化結構選擇部分原子振動的弛豫。零點能主要用來修正吧,計算公式你可以參考該文章及其引用文獻,自己多搜一些吧,表面吸附考慮零點能的很多。http://www.sciencedirect.com/sci ... i/S0360319913024518 |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學070300 學碩336求調劑 +4 | vv迷 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 生物學071000 329分求調劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|
[考研] 一志愿華中科技大學071000,求調劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 269專碩求調劑 +5 | 金恩貝 2026-03-21 | 5/250 |
|
|
[考研] 311求調劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|
[考研] 354求調劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|
[考研] 0703化學調劑 ,六級已過,有科研經(jīng)歷 +14 | 曦熙兮 2026-03-15 | 14/700 |
|
|
[考研] 0805 316求調劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 材料工程專碩 348分求調劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 296求調劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 0856材料專碩353求調劑 +3 | NIFFFfff 2026-03-20 | 3/150 |
|
|
[考研] 南昌大學材料專碩311分求調劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 求調劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 材料 336 求調劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩,總分308求調劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 工科材料085601 279求調劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 293求調劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|