| 10 | 1/1 | 返回列表 |
| 查看: 3628 | 回復(fù): 9 | |||
fangyongxinxi新蟲 (正式寫手)
|
[求助]
VASP優(yōu)化結(jié)構(gòu)遇到的問題。
|
|
在用vasp 優(yōu)化一個體系,INCAR 文件如下。優(yōu)化完后,出現(xiàn)CONTCAR 文件,即優(yōu)化成功后的文件。重新建一個文件夾,把前面得到的CONTCAR文件改為POSCAR文件,用同樣的參數(shù)和方法再進行一次優(yōu)化,得到CONTCAR_2文件。 問題: 1. 理論上這兩個CONTCAR文件中的值應(yīng)該相差很小,但晶體的a b c 的值卻相差不較大,而原子的分?jǐn)?shù)坐標(biāo)沒有差別,不知道為什么,或是哪里參數(shù)設(shè)置出了問題。 2. 自己這樣驗證的目的是想看看VASP對結(jié)構(gòu)優(yōu)化的有效性,所以還有個問題,即如何正確的優(yōu)化一個體系。 system ISTART = 0 ICHARG = 2 ISPIN =2 VOSKOWN = 1 ENCUT = 302.00 eV ISMEAR = 0 SIGMA = 0.1 EDIFF = 0.1E-05 NSW = 100 IBRION = 2 ISIF = 3 LORBIT = 11 |
新蟲 (正式寫手)
|
這是待優(yōu)化的結(jié)構(gòu): TiNi B2 1.00000000000000 5.7516 0.0000000000000000 0.0000000000000000 0.0000000000000000 4.06762 0.0000000000000000 0.0000000000000000 0.0000000000000000 4.06762 Ni Ti 4 4 Direct 0.0000000000000000 0.7500000000000000 0.2500000000000000 0.5 0.75 0.25 0.0000000000000000 0.2500000000000000 0.7500000000000000 0.5 0.25 0.75 0.25000000000000000 0.2500000000000000 0.2500000000000000 0.75 0.25 0.25 0.25000000000000000 0.7500000000000000 0.7500000000000000 0.75 0.75 0.75 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 這是用上面INCAR文件參數(shù)后,得到的CONTCAR NiTi B2 1.00000000000000 5.6274121763945830 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.9757576568911093 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.9757576568911093 Ni Ti 4 4 Direct 0.0000000000000000 0.7500000000000000 0.2500000000000000 0.5000000000000000 0.7500000000000000 0.2500000000000000 0.0000000000000000 0.2500000000000000 0.7500000000000000 0.5000000000000000 0.2500000000000000 0.7500000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.7500000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 將上面第一次優(yōu)化的CONTCAR 改名為 poscar ,用同樣的INCAR, 再次優(yōu)化得到contcar_2 NiTi B2 1.00000000000000 5.5636005577964589 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.9265497075565161 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.9265497075565161 Ni Ti 4 4 Direct 0.0000000000000000 0.7500000000000000 0.2500000000000000 0.5000000000000000 0.7500000000000000 0.2500000000000000 0.0000000000000000 0.2500000000000000 0.7500000000000000 0.5000000000000000 0.2500000000000000 0.7500000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.7500000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 顯然兩個CONTCAR的結(jié)果差別很大: 5.6274121763945830 vs 5.5636005577964589 悲劇,不知道錯在哪里。 |
|
體積比是88.9vs85.8,是有點大,看看兩者優(yōu)化的總能的大小,取小的一個再進行優(yōu)化,有關(guān)優(yōu)化可以看http://www.gaoyang168.com/bbs/viewthread.php?tid=3234667&fpage=1, 不要直接用isif=3,這樣會很不準(zhǔn)確 |

新蟲 (正式寫手)

鐵蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
木蟲 (正式寫手)

金蟲 (小有名氣)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
孅華 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 336化工調(diào)劑 +3 | 王大坦1 2026-03-23 | 4/200 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|