swx0789

應(yīng)助: 0 (幼兒園)
金幣: 12.1
帖子: 124
在線: 512小時(shí)
蟲(chóng)號(hào): 648879

[交流] 【求助】結(jié)構(gòu)優(yōu)化

各位前輩，小弟優(yōu)化MoO3晶體結(jié)構(gòu)（已有原始數(shù)據(jù)a=3.9628，b=13.8550，c=3.6964；alpha=beta=gama=90度)。我按照候博士的步驟：
1，先優(yōu)化得到晶胞最佳體積；
2，然后再在這個(gè)最佳體積下用ISIF=4;NSW=120;IBRION=1進(jìn)行結(jié)構(gòu)馳豫；
3，用第二步的到的CONTCAR作為POSCAR做靜態(tài)計(jì)算，發(fā)現(xiàn)OUTCAR里面的：
  external pressure =    11.18 kB  Pullay stress =       0.00 kB;
4，因?yàn)榈?步一直計(jì)算到120步才結(jié)束，我再在第2步的基礎(chǔ)上從新用ISIF=4;NSW=120;IBRION=1進(jìn)行結(jié)構(gòu)馳豫；
5，用第4步的到的CONTCAR作為POSCAR做靜態(tài)計(jì)算，發(fā)現(xiàn)OUTCAR里面的：
  external pressure =       0.12 kB  Pullay stress =       0.00 kB

下面我有幾個(gè)問(wèn)題想問(wèn)各位前輩：
1，一個(gè)結(jié)構(gòu)優(yōu)化好是不是我第3，5步的“external pressure =    0 kB”為標(biāo)志，或者接近0kB？
2，你們覺(jué)得我這個(gè)結(jié)構(gòu)算優(yōu)化好了嗎？
3，我在第2步進(jìn)行馳豫并沒(méi)有收斂，它一直計(jì)算到120步，我在第5步再繼續(xù)計(jì)算120步還是一樣不收斂。但我直接還是用CONTACR作為POSCAR做靜態(tài)計(jì)算，你們覺(jué)得這里面有問(wèn)題嗎？

[ Last edited by swx0789 on 2011-2-17 at 09:46 ]

回復(fù)此樓

» 猜你喜歡

317求調(diào)劑已經(jīng)有12人回復(fù)
一志愿武漢理工材料工程專(zhuān)碩調(diào)劑已經(jīng)有9人回復(fù)
【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑已經(jīng)有7人回復(fù)
324求調(diào)劑已經(jīng)有5人回復(fù)
288求調(diào)劑已經(jīng)有16人回復(fù)
317求調(diào)劑已經(jīng)有9人回復(fù)
287求調(diào)劑已經(jīng)有8人回復(fù)
一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。已經(jīng)有4人回復(fù)
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分已經(jīng)有4人回復(fù)
085600材料與化工已經(jīng)有8人回復(fù)

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

結(jié)構(gòu)優(yōu)化 a=b ----> a不等于b 已經(jīng)有3人回復(fù)
【課件】結(jié)構(gòu)優(yōu)化設(shè)計(jì)（崔昌禹，哈爾濱工業(yè)大學(xué)）已經(jīng)有65人回復(fù)
MS 結(jié)構(gòu)優(yōu)化的結(jié)果不能理解求高人指點(diǎn) 已經(jīng)有10人回復(fù)
vasp 結(jié)構(gòu)優(yōu)化問(wèn)題已經(jīng)有14人回復(fù)
請(qǐng)教初始結(jié)構(gòu)優(yōu)化問(wèn)題已經(jīng)有10人回復(fù)
分子結(jié)構(gòu)優(yōu)化已經(jīng)有8人回復(fù)
結(jié)構(gòu)優(yōu)化中問(wèn)題已經(jīng)有6人回復(fù)
【活動(dòng)or資源】vasp晶體結(jié)構(gòu)優(yōu)化已經(jīng)有142人回復(fù)
【求助結(jié)束】結(jié)構(gòu)優(yōu)化很離奇，請(qǐng)問(wèn)可能是什么原因？--初始構(gòu)型不合理已經(jīng)有18人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化不收斂已經(jīng)有9人回復(fù)
【原創(chuàng)】關(guān)于結(jié)構(gòu)優(yōu)化的總結(jié)，體會(huì) 已經(jīng)有33人回復(fù)
【求助】Dmol3結(jié)構(gòu)優(yōu)化已經(jīng)有8人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-02-17 09:37:53

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

215215

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
帖子: 125
在線: 53.5小時(shí)
蟲(chóng)號(hào): 860062

關(guān)注

回復(fù)此樓

2樓2011-02-17 13:03:31

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

tider

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時(shí)
蟲(chóng)號(hào): 1190614

★
sunyang1988(金幣+1): 謝謝交流 2011-04-10 00:02:14

問(wèn)題問(wèn)的好啊�。№敚。�！
我也不明白external pressure是什么意義呢，只知道Pullay force。記得manual上講external pressure<0 就好。但不曉得為什么。
有沒(méi)有高手給點(diǎn)見(jiàn)解啊

贊一下(1人)

回復(fù)此樓

3樓2011-04-09 22:50:41

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

yangprc

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2439
帖子: 103
在線: 179.5小時(shí)
蟲(chóng)號(hào): 435768

關(guān)注

回復(fù)此樓

4樓2011-04-10 16:38:51

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

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
帖子: 565
在線: 244.3小時(shí)
蟲(chóng)號(hào): 853446

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

樓主解決這個(gè)問(wèn)題了嗎？我也想知道

贊一下(1人)

回復(fù)此樓

5樓2011-04-24 22:18:38

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

swx0789

銅蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 12.1
帖子: 124
在線: 512小時(shí)
蟲(chóng)號(hào): 648879

引用回帖:

Originally posted by donggaomin at 2011-04-24 22:18:38:
樓主解決這個(gè)問(wèn)題了嗎？我也想知道

你是多external嗎

平衡結(jié)夠的external為0，我一般也有幾kb吧

贊一下

回復(fù)此樓

6樓2011-04-27 14:01:08

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

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
帖子: 565
在線: 244.3小時(shí)
蟲(chóng)號(hào): 853446

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖
bingmou(金幣+1): 感謝交流！ 2011-05-06 22:53:23

引用回帖:

Originally posted by swx0789 at 2011-04-27 14:01:08:
你是多external嗎

平衡結(jié)夠的external為0，我一般也有幾kb吧

是external pressure, 一般有負(fù)的幾KB，另外，in kB這一行前三個(gè)數(shù)也有幾個(gè)KB;

偶爾會(huì)有大的，如一次全優(yōu)化晶格常數(shù)和位置：
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.2
NSW = 60; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 3
POTIM = 0.2
PREC = Accurate
LWAVE = .FALSE.
LCHARG =.FALSE.

OUTCAR:
Total    -2.27547 -2.27547 -1.52158    0.00000    0.00000    0.000
00
  in kB    -57.62827 -57.62827 -38.53531    0.00000    0.00000    0.000
00
  external pressure =    -51.26 kB  Pullay stress =       0.00 kB

這時(shí)，是不是說(shuō)沒(méi)有全優(yōu)化好��？

贊一下(2人)

回復(fù)此樓

7樓2011-05-02 15:48:28

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

swx0789

銅蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 12.1
帖子: 124
在線: 512小時(shí)
蟲(chóng)號(hào): 648879

引用回帖:

Originally posted by donggaomin at 2011-05-02 15:48:28:
是external pressure, 一般有負(fù)的幾KB，另外，in kB這一行前三個(gè)數(shù)也有幾個(gè)KB;

偶爾會(huì)有大的，如一次全優(yōu)化晶格常數(shù)和位置：
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEA ...

你的這個(gè)收斂了嗎？你設(shè)置的NSW = 60，到60步收斂了嗎？我想你的這個(gè)結(jié)構(gòu)不一定收斂了，我看你設(shè)的其他參數(shù)比較的苛刻。

贊一下

回復(fù)此樓

8樓2011-05-04 08:37:01

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

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
帖子: 565
在線: 244.3小時(shí)
蟲(chóng)號(hào): 853446

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

Originally posted by swx0789 at 2011-05-04 08:37:01:
你的這個(gè)收斂了嗎？你設(shè)置的NSW = 60，到60步收斂了嗎？我想你的這個(gè)結(jié)構(gòu)不一定收斂了，我看你設(shè)的其他參數(shù)比較的苛刻。

對(duì)的，確實(shí)NSW = 60可能不夠。另外，您說(shuō)“其他參數(shù)比較的苛刻”比如哪一個(gè)？請(qǐng)您指點(diǎn)，因?yàn)槲易约哼@樣用習(xí)慣了，呵呵

贊一下(1人)

回復(fù)此樓

9樓2011-05-06 22:48:43

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

tongji888

禁蟲(chóng) (小有名氣)

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

本帖內(nèi)容被屏蔽

10樓2011-05-07 13:07:49

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

tider

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時(shí)
蟲(chóng)號(hào): 1190614

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

Originally posted by tongji888 at 2011-05-07 13:07:49:
可以問(wèn)下你對(duì)Pullay force的理解嗎？謝謝。。。。

嗯，manual后面有講到的，你找找kan

贊一下(1人)

回復(fù)此樓

11樓2011-05-07 22:51:48

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

tider

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時(shí)
蟲(chóng)號(hào): 1190614

我問(wèn)過(guò)自己的導(dǎo)師了。他說(shuō)extential energy 就是晶體結(jié)構(gòu)中的內(nèi)應(yīng)力。

贊一下

回復(fù)此樓

12樓2011-05-07 22:53:14

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

tider

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時(shí)
蟲(chóng)號(hào): 1190614

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

Originally posted by swx0789 at 2011-02-17 09:37:53:
各位前輩，小弟優(yōu)化MoO3晶體結(jié)構(gòu)（已有原始數(shù)據(jù)a=3.9628，b=13.8550，c=3.6964；alpha=beta=gama=90度)。我按照候博士的步驟：
1，先優(yōu)化得到晶胞最佳體積；
2，然后再在這個(gè)最佳體積下用ISIF=4;NSW=120;IBRION ...

聽(tīng)學(xué)長(zhǎng)說(shuō)一般不會(huì)馳豫到100多步，你可能是晶格常數(shù)的初始值設(shè)置得偏離平衡位置太遠(yuǎn)了吧。不過(guò)NSW設(shè)置多些應(yīng)該總沒(méi)壞處。還有我的導(dǎo)師建議我馳豫時(shí)ＩＢＲＩＯＮ取２，這樣可能馳豫速度要快些。當(dāng)然我不曉得為什么。。算法的問(wèn)題。。。

贊一下(1人)

回復(fù)此樓

13樓2011-05-07 22:57:19

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

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

[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽(yáng)
[考研] 350求調(diào)劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿武理材料工程348求調(diào)劑 +3	￣^￣゜汗 2026-03-19	4/200	2026-03-20 21:01 by zhukairuo
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 279分求調(diào)劑一志愿211 +10	chaojifeixia 2026-03-19	11/550	2026-03-20 20:32 by JourneyLucky
[考研] 材料與化工專(zhuān)碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無(wú)懈可擊111
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 0703化學(xué)調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學(xué)員8dgXkO
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 化學(xué)求調(diào)劑 +3	臨澤境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無(wú)懈可擊111
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染