| 24小時(shí)熱門(mén)版塊排行榜 |
| 8 | 1/1 | 返回列表 |
| 查看: 972 | 回復(fù): 7 | ||||
xj20120828木蟲(chóng) (正式寫(xiě)手)
|
[求助]
頻率驗(yàn)證 已有1人參與
|
|
大家好,我用vasp計(jì)算完過(guò)渡態(tài)想驗(yàn)證一下頻率,設(shè)置 PREC=Normal ISTART=0 ICHARG=2 ISPIN=1 LREAL=AUTO LWAVE=.FALSE. LCHARG=.FALSE. ISYM=0 Electronic Relaxation ENCUT=500 NELM=60 EDIFF=0.1E-3 #NELMIN=4 GGA=PE VOSKOWN=1 Ionic relaxtion EDIFFG=-.1E-0 NSW=1 IBRION=5 ISIF=2 POTIM=0.2 DOS related values: ISMEAR=0 !Gaussian SIGMA=0.02 Magnetism #MAGMOM=8*0 18*0 18*0 3*2 LORBIT=10 #NUPDOWN=2 #For NEB calculation #IMAGES=4 #LCLIMB=.True. #SPRING=-5 #ICHAIN=0 #LSCALAPACK = .FALSE. 計(jì)算頻率,為什么設(shè)置NSW=1,還是會(huì)進(jìn)行許多步的離子位置優(yōu)化?不是應(yīng)該只有一步嗎? |
木蟲(chóng) (正式寫(xiě)手)
專家顧問(wèn) (著名寫(xiě)手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
For IBRION=5, all atoms are displaced in all three Cartesian directions, resulting in a significant computational effort even for moderately sized high symmetry systems. 因?yàn)槟汩_(kāi)了 IBRION=5 VASP 會(huì)對(duì)每一顆"可允許移動(dòng)"的原子 計(jì)算 x y z 三個(gè)方向移動(dòng) POTIM 距離後的受力 所以 你有 N 顆 "允許移動(dòng)" 的原子 就會(huì)有 3N 個(gè)離子步計(jì)算 (非優(yōu)化) |
木蟲(chóng) (正式寫(xiě)手)
專家顧問(wèn) (著名寫(xiě)手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
"為什么我用優(yōu)化好的結(jié)構(gòu)計(jì)算頻率也會(huì)算出來(lái)虛頻呢?" ㄟ~ 那就表示 你優(yōu)化的不夠好 做 PHONON 計(jì)算需要的收斂條件是相對(duì)嚴(yán)苛的 PHONON 計(jì)算事前優(yōu)化計(jì)算所需要的一個(gè)重要條件 就是 EDIFFG 最好能在 1E-4 或 1E-5 這參數(shù)對(duì)總能影響很小 但~對(duì)原子受力的精準(zhǔn)度有很大幫助 如果你只有一個(gè)虛頻 表示你計(jì)算的結(jié)構(gòu)可能在暫穩(wěn)態(tài) 已經(jīng)很接近 global minimum 如果是兩個(gè)以上的虛頻 ( 負(fù)值如果很小,可忽略不算) 那表示你結(jié)構(gòu)離最穩(wěn)態(tài)還有很長(zhǎng)一條路要走 |
木蟲(chóng) (正式寫(xiě)手)
|
EDIFFG要那么高的精度嗎?還是EDIFF? 我優(yōu)化的是EDIFFG=-0.01。我算出來(lái)四個(gè)虛頻, 171 f/i= 0.389539 THz 2.447547 2PiTHz 12.993627 cm-1 1.611005 meV 172 f/i= 1.438074 THz 9.035684 2PiTHz 47.968978 cm-1 5.947397 meV 173 f/i= 2.256721 THz 14.179397 2PiTHz 75.276111 cm-1 9.333052 meV 174 f/i= 22.875331 THz 143.729946 2PiTHz 763.038898 cm-1 94.604800 meV 68個(gè)原子。請(qǐng)問(wèn)大概多少頻率就算比較大的虛頻了? |
專家顧問(wèn) (著名寫(xiě)手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
??o?ε? PHONON ????????@?N??? EDIFFG EDIFFG=0.01 ????????????????????? ???????M??PHONON??????? ?????l????? ???е? 171, 172, 173 ????"????"???? 0 ??????????????l?????? 174 ?? (~94 meV) 68??????????????????????????????? > ??l??С?????????o?P ?f???????"?^??"??l...??~?????? ??????X 10 meV ????????????? 0 ???? ????l??С??l?g??????~??С????l"???"??????????l ????N??????????????l (????]??L??^??????????????????) ?z????l?? vibration mode (no. 174) ?? OUTCAR ??????????? frequency ???????? phonon mode ?? POSCAR????????l?? phonon mode ?????????????С??λ?? (???С??? ??? 0.02 ??????????????????????С.....) ??????????? ??r???? EDIFFG ? 1E-04 ????????????? phonon ??? ???ЙC????????l????? |
木蟲(chóng) (正式寫(xiě)手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 280求調(diào)劑 +7 | 咕嚕曉曉 2026-03-18 | 8/400 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|