| 8 | 1/1 | 返回列表 |
| 查看: 3467 | 回復(fù): 7 | |||
cindywei木蟲(chóng) (初入文壇)
|
[求助]
VASP計(jì)算穩(wěn)定體系的零點(diǎn)能為什么也會(huì)出現(xiàn)虛頻?
|
|
最近用VASP計(jì)算Gamma點(diǎn)的零點(diǎn)能做校正,原本熱力學(xué)及力學(xué)上都能穩(wěn)定存在的鈀氫化合物有些成分點(diǎn)出現(xiàn)了虛頻,是哪里出現(xiàn)問(wèn)題了嗎? INCAR設(shè)置 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 |
木蟲(chóng) (正式寫(xiě)手)
|
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. |
木蟲(chóng) (小有名氣)
|
cindywei您好,我想問(wèn)您一個(gè)問(wèn)題就是我要做零點(diǎn)能校正的話,按E=1/2hv公式,把所有頻率的能量相加,大約就是0.7eV,按理說(shuō),零點(diǎn)能校正的話,不會(huì)是這么大的值,也就是0.02eV,是不是這中計(jì)算方法有誤呢? 我如果要對(duì)氣態(tài)分子吸附態(tài)進(jìn)行零點(diǎn)能校正的話,比如H2O在Cu上的吸附,要計(jì)算頻率,我是不是應(yīng)該把H2O和Cu的第一層原子弛豫,把所有的頻率的能量相加,就算是零點(diǎn)能呢?是不是應(yīng)該把Cu原子弛豫的頻率去掉,只加H2O的頻率呢,但是那樣的話,我根本就不能從所有的頻率中區(qū)分H2O和Cu的頻率。 最近甚是苦惱啊,還望您能指點(diǎn)迷津,期待您的回復(fù),在此,萬(wàn)分感謝。。。 |

木蟲(chóng) (初入文壇)
|
計(jì)算零點(diǎn)能的方法有很多種吧,如果你只考慮輕質(zhì)的H及O的零點(diǎn)能的話,可以在POSCAR選擇H及O不固定而Cu固定,這種計(jì)算不完全等同于晶格優(yōu)化的弛豫,是每個(gè)不固定原子按特定幾個(gè)方向有限的弛豫來(lái)獲得其零點(diǎn)振動(dòng)能。我跟個(gè)INCAR的設(shè)置,你可以通過(guò)參數(shù)來(lái)理解一下,希望對(duì)你有用 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 |
木蟲(chóng) (小有名氣)

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

木蟲(chóng) (初入文壇)
|
是先把整個(gè)結(jié)構(gòu)弛豫優(yōu)化獲得最優(yōu)化結(jié)構(gòu),然后對(duì)優(yōu)化結(jié)構(gòu)選擇部分原子振動(dòng)的弛豫。零點(diǎn)能主要用來(lái)修正吧,計(jì)算公式你可以參考該文章及其引用文獻(xiàn),自己多搜一些吧,表面吸附考慮零點(diǎn)能的很多。http://www.sciencedirect.com/sci ... i/S0360319913024518 |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 0856材料專(zhuān)碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 材料專(zhuān)碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|