swx0789

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 12.1
帖子: 124
在線: 512小時
蟲號: 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，然后再在這個最佳體積下用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

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

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

回復(fù)此樓

» 猜你喜歡

290求調(diào)劑已經(jīng)有3人回復(fù)
一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有4人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有8人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有8人回復(fù)
317求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有5人回復(fù)
274求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有9人回復(fù)
招收調(diào)劑碩士已經(jīng)有11人回復(fù)

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

結(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)化問題已經(jīng)有14人回復(fù)
請教初始結(jié)構(gòu)優(yōu)化問題已經(jīng)有10人回復(fù)
分子結(jié)構(gòu)優(yōu)化已經(jīng)有8人回復(fù)
結(jié)構(gòu)優(yōu)化中問題已經(jīng)有6人回復(fù)
【活動or資源】vasp晶體結(jié)構(gòu)優(yōu)化已經(jīng)有142人回復(fù)
【求助結(jié)束】結(jié)構(gòu)優(yōu)化很離奇，請問可能是什么原因？--初始構(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é)，體會已經(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

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
帖子: 125
在線: 53.5小時
蟲號: 860062

關(guān)注

回復(fù)此樓

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

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

tider

木蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時
蟲號: 1190614

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

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

贊一下(1人)

回復(fù)此樓

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

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

yangprc

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2439
帖子: 103
在線: 179.5小時
蟲號: 435768

關(guān)注

回復(fù)此樓

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

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
帖子: 565
在線: 244.3小時
蟲號: 853446

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

樓主解決這個問題了嗎？我也想知道

贊一下(1人)

回復(fù)此樓

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

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

swx0789

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 12.1
帖子: 124
在線: 512小時
蟲號: 648879

引用回帖:

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

你是多external嗎

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

贊一下

回復(fù)此樓

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

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
帖子: 565
在線: 244.3小時
蟲號: 853446

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
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這一行前三個數(shù)也有幾個KB;

偶爾會有大的，如一次全優(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

這時，是不是說沒有全優(yōu)化好�。�

贊一下(2人)

回復(fù)此樓

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

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

swx0789

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 12.1
帖子: 124
在線: 512小時
蟲號: 648879

引用回帖:

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

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

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

贊一下

回復(fù)此樓

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

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
帖子: 565
在線: 244.3小時
蟲號: 853446

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

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

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

贊一下(1人)

回復(fù)此樓

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

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

tongji888

禁蟲 (小有名氣)

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

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

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

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

tider

木蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時
蟲號: 1190614

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

Originally posted by tongji888 at 2011-05-07 13:07:49:
可以問下你對Pullay force的理解嗎？謝謝。。。。

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

贊一下(1人)

回復(fù)此樓

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

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

tider

木蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時
蟲號: 1190614

我問過自己的導(dǎo)師了。他說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

木蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 1950.4
帖子: 313
在線: 87.2小時
蟲號: 1190614

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

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，然后再在這個最佳體積下用ISIF=4;NSW=120;IBRION ...

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

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	8/400	2026-03-20 11:24 by kkcoco25
[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	8/400	2026-03-20 11:20 by 申子申申
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +8	步川酷紫123 2026-03-13	8/400	2026-03-20 08:47 by xingguangj
[考研] 0703化學(xué)調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學(xué)員8dgXkO
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky