當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

fangyongxinxi

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 41
散金: 69
紅花: 3
帖子: 876
在線: 138.5小時(shí)
蟲號: 437418
注冊: 2007-08-26
性別: MM
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] VASP優(yōu)化結(jié)構(gòu)遇到的問題。

在用vasp 優(yōu)化一個(gè)體系，INCAR 文件如下。優(yōu)化完后，出現(xiàn)CONTCAR 文件，即優(yōu)化成功后的文件。重新建一個(gè)文件夾，把前面得到的CONTCAR文件改為POSCAR文件，用同樣的參數(shù)和方法再進(jìn)行一次優(yōu)化，得到CONTCAR_2文件。
問題：
1. 理論上這兩個(gè)CONTCAR文件中的值應(yīng)該相差很小，但晶體的a b c 的值卻相差不較大，而原子的分?jǐn)?shù)坐標(biāo)沒有差別，不知道為什么，或是哪里參數(shù)設(shè)置出了問題。
2. 自己這樣驗(yàn)證的目的是想看看VASP對結(jié)構(gòu)優(yōu)化的有效性，所以還有個(gè)問題，即如何正確的優(yōu)化一個(gè)體系。

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

回復(fù)此樓

» 猜你喜歡

11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) 已經(jīng)有6人回復(fù)
求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分已經(jīng)有8人回復(fù)
環(huán)境工程 085701，267求調(diào)劑已經(jīng)有7人回復(fù)
材料工程專碩求調(diào)劑已經(jīng)有8人回復(fù)
297求調(diào)劑已經(jīng)有3人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有7人回復(fù)
281求調(diào)劑已經(jīng)有7人回復(fù)
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行已經(jīng)有8人回復(fù)
求收留已經(jīng)有6人回復(fù)
085600 295分求調(diào)劑已經(jīng)有9人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

vasp 結(jié)構(gòu)優(yōu)化問題已經(jīng)有14人回復(fù)
vasp5.2 openmpi+intel+mkl編譯，能量正值問題已經(jīng)有7人回復(fù)
求助：vasp結(jié)構(gòu)優(yōu)化不能收斂已經(jīng)有18人回復(fù)
vasp計(jì)算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計(jì)算不收斂已經(jīng)有3人回復(fù)
vasp優(yōu)化結(jié)構(gòu)中出現(xiàn)的問題已經(jīng)有10人回復(fù)
【活動(dòng)or資源】vasp晶體結(jié)構(gòu)優(yōu)化已經(jīng)有142人回復(fù)
【求助】vasp 高壓相結(jié)構(gòu)優(yōu)化已經(jīng)有17人回復(fù)
【求助】vasp優(yōu)化結(jié)構(gòu)未結(jié)束，log文件出現(xiàn)下列問題停止計(jì)算，請高手指教已經(jīng)有6人回復(fù)
【求助】vasp靜態(tài)計(jì)算能量值錯(cuò)誤已經(jīng)有21人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化后，靜態(tài)計(jì)算竟然失��？已經(jīng)有8人回復(fù)
【討論】vasp結(jié)構(gòu)優(yōu)化之前的測試已經(jīng)有11人回復(fù)

1樓 2011-06-06 10:39:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kwanqsh

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1625.9
帖子: 84
在線: 133.5小時(shí)
蟲號: 2802650
注冊: 2013-11-15
專業(yè): 實(shí)驗(yàn)地球化學(xué)和計(jì)算地球化

請問問題解決了嗎我也遇到了誒

贊一下

回復(fù)此樓

9樓2014-03-17 15:45:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 10 個(gè)回答

fangyongxinxi

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 41
散金: 69
紅花: 3
帖子: 876
在線: 138.5小時(shí)
蟲號: 437418
注冊: 2007-08-26
性別: MM
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

Originally posted by fangyongxinxi at 2011-06-06 10:39:57:
在用vasp 優(yōu)化一個(gè)體系，INCAR 文件如下。優(yōu)化完后，出現(xiàn)CONTCAR 文件，即優(yōu)化成功后的文件。重新建一個(gè)文件夾，把前面得到的CONTCAR文件改為POSCAR文件，用同樣的參數(shù)和方法再進(jìn)行一次優(yōu)化，得到CONTCAR_2文件。 ...

這是待優(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

顯然兩個(gè)CONTCAR的結(jié)果差別很大：
5.6274121763945830    vs 5.5636005577964589
悲劇，不知道錯(cuò)在哪里。

贊一下

回復(fù)此樓

2樓2011-06-06 10:45:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

1ST強(qiáng)帖: 5
應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號: 995125
注冊: 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
fangyongxinxi(金幣+6): 2011-06-06 18:59:44
franch(金幣+2): 謝謝回帖交流 2011-06-06 20:05:57

體積比是88.9vs85.8，是有點(diǎn)大，看看兩者優(yōu)化的總能的大小，取小的一個(gè)再進(jìn)行優(yōu)化，有關(guān)優(yōu)化可以看http://www.gaoyang168.com/bbs/viewthread.php?tid=3234667&fpage=1，
不要直接用isif=3，這樣會(huì)很不準(zhǔn)確

贊一下(1人)

回復(fù)此樓

嘟啦啦嘟啦啦嘟

3樓2011-06-06 14:54:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fangyongxinxi

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 41
散金: 69
紅花: 3
帖子: 876
在線: 138.5小時(shí)
蟲號: 437418
注冊: 2007-08-26
性別: MM
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

經(jīng)過再一次驗(yàn)證，悲劇再次出現(xiàn)
將第一次優(yōu)化的結(jié)果，用作第二次優(yōu)化的結(jié)構(gòu)輸入?yún)?shù)，
發(fā)現(xiàn)這兩次優(yōu)化的能力和結(jié)構(gòu)參數(shù)差別很大。（如上面所說）

but
如果只優(yōu)化初始結(jié)構(gòu)，無論選什么參數(shù)，得到的能力和結(jié)構(gòu)參數(shù)都差不多。

是不是可以這樣說：
vasp的重復(fù)性好。
但如何解釋我遇到的問題？
蒼天，幫幫我。

贊一下

回復(fù)此樓

4樓2011-06-08 09:40:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +8	@taotao 2026-03-29	8/400	2026-03-30 20:48 by dophin1985
[考研] 085600 材料與化工 329分求調(diào)劑 +18	Mr. Z 2026-03-25	19/950	2026-03-30 19:01 by Wang200018
[考研] 105500藥學(xué)求調(diào)劑，一志愿山東大學(xué)藥學(xué)，348分 +3	gr哈哈哈 2026-03-28	3/150	2026-03-30 18:56 by 源_2020
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 材料專碩 085600求調(diào)劑 +7	BBQ233 2026-03-30	7/350	2026-03-30 17:44 by oooqiao
[考研] 0703化學(xué) +14	妮妮ninicgb 2026-03-27	14/700	2026-03-30 12:46 by fangnagu
[考研] 311求調(diào)劑 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 學(xué)碩274求調(diào)劑 +9	Li李魚 2026-03-26	9/450	2026-03-28 21:42 by bymhappy
[考研] 求調(diào)劑 +7	爭取九點(diǎn)睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點(diǎn)睡
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 085600，材料與化工321分，求調(diào)劑 +9	大饞小子 2026-03-27	9/450	2026-03-27 14:30 by mmm just
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike