whw19850730

應(yīng)助: 1 (幼兒園)
金幣: 5875.8
紅花: 1
帖子: 64
在線: 121.3小時(shí)
蟲(chóng)號(hào): 1048864
注冊(cè): 2010-06-28
專業(yè): 物理學(xué)I

[求助] VASP 電子迭代發(fā)散詭異問(wèn)題

我做一個(gè)六角晶體LuMnO3，電子迭代發(fā)散，遇到不收斂問(wèn)題可以增加 K點(diǎn)，改變
電荷密度混合參數(shù)等很多辦法，可是我的不是迭代震蕩，而是迭代能量越來(lái)越大最后崩潰掉，我調(diào)了很久都解決不掉，希望高手大牛給些意見(jiàn)。
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.242426654415E+03 -0.24243E+03 -0.49552E+02 12636 0.152E+01
DAV: 2 -0.244734802042E+03 -0.23081E+01 -0.13019E+01  8820 0.106E+00
DAV: 3 -0.244911974268E+03 -0.17717E+00 -0.17362E+00 12996 0.458E-01
DAV: 4 -0.244917066602E+03 -0.50923E-02 -0.45943E-02  8136 0.601E-02
DAV: 5 -0.244917150324E+03 -0.83722E-04 -0.83640E-04  5904 0.106E-02
DAV: 6 -0.247540272456E+03 -0.26231E+01 -0.31613E+00  5976 0.408E-03
DAV: 7 -0.246111205273E+03 0.14291E+01 -0.36981E+01  5760 0.395E-01
DAV: 8 -0.260796035821E+03 -0.14685E+02 -0.25748E+01  5760 0.277E-01 0.148E+00
DAV: 9 -0.244897893351E+03 0.15898E+02 -0.23689E+01 12600 0.638E+00 0.122E+00
DAV:  10 -0.244879497047E+03 0.18396E-01 -0.12455E-01 12564 0.182E+01 0.761E-01
DAV:  11 -0.244883124169E+03 -0.36271E-02 -0.56369E-02 12564 0.110E+01 0.688E-01
DAV:  12 -0.244879821921E+03 0.33022E-02 -0.94784E-03 13140 0.449E+00 0.367E-01
DAV:  13 -0.245058699872E+03 -0.17888E+00 -0.16888E+00 13032 0.277E+00 0.294E-01
DAV:  14 -0.252709857733E+03 -0.76512E+01 -0.84904E+01 12960 0.465E+00 0.152E-01
DAV:  15 -0.245505026950E+03 0.72048E+01 -0.31400E+01 13068 0.124E+00 0.641E-02
DAV:  16 -0.245391127504E+03 0.11390E+00 -0.13041E+01 12888 0.127E+00 0.531E-02
DAV:  17 -0.244879335901E+03 0.51179E+00 -0.69883E+00 12348 0.371E-01 0.707E-02
DAV:  18 -0.267404373109E+03 -0.22525E+02 -0.22226E+02 12996 0.465E-01 0.757E-01
DAV:  19 -0.244846617464E+03 0.22558E+02 -0.92693E+01 10080 0.445E+00 0.452E-02
DAV:  20 -0.317293990069E+03 -0.72447E+02 -0.71922E+02 13284 0.191E-01 0.566E-01
DAV:  21 -0.520146965043E+05 -0.51697E+05 -0.51779E+05  9288 0.219E+00 0.187E+00
DAV:  22 -0.248042514172E+03 0.51767E+05 -0.54223E+02  8028 0.761E+00 0.460E-01
DAV:  23 -0.543114596657E+04 -0.51831E+04 -0.63755E+02 10764 0.247E+00 0.164E+00
DAV:  24 -0.589055666359E+03 0.48421E+04 -0.36204E+03 10296 0.657E+00 0.713E-01
DAV:  25 -0.308703180751E+06 -0.30811E+06 -0.24597E+04 10044 0.731E+00
1 F= -.30870318E+06 E0= -.30870318E+06  d E =-.308703E+06  mag=    0.0000

回復(fù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(fù)

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

關(guān)于electrostatic potential 已經(jīng)有5人回復(fù)
CFX 能量方程迭代發(fā)散已經(jīng)有5人回復(fù)
一個(gè)很菜的問(wèn)題。。。管道流量或速度的計(jì)算？已經(jīng)有4人回復(fù)
在做fluent例子的時(shí)候遇到了問(wèn)題，發(fā)現(xiàn)無(wú)法進(jìn)行迭代計(jì)算，求教大家已經(jīng)有9人回復(fù)
請(qǐng)問(wèn)~Fluent 中兩相流如何設(shè)置各自的壓力入口？另外，迭代后發(fā)現(xiàn)發(fā)散咋辦？已經(jīng)有4人回復(fù)

1樓 2013-08-17 02:47:00

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

oupengfei717

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

應(yīng)助: 3 (幼兒園)
金幣: 744.8
散金: 80
紅花: 3
帖子: 258
在線: 144小時(shí)
蟲(chóng)號(hào): 1435139
注冊(cè): 2011-10-10
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

感謝參與，應(yīng)助指數(shù) +1

樓主這個(gè)是靜態(tài)計(jì)算嗎？怎么只有一步離子迭代計(jì)算呢？我猜想可能是該結(jié)構(gòu)在初期全弛豫的時(shí)候沒(méi)有能找到平衡結(jié)構(gòu)？不知道樓主是否進(jìn)行了弛豫得到了平衡結(jié)構(gòu)？樓主的輸入文件是否都檢查過(guò)不存在問(wèn)題呢？

贊一下

回復(fù)此樓

2樓2013-08-17 10:31:52

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

whw19850730

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

應(yīng)助: 1 (幼兒園)
金幣: 5875.8
紅花: 1
帖子: 64
在線: 121.3小時(shí)
蟲(chóng)號(hào): 1048864
注冊(cè): 2010-06-28
專業(yè): 物理學(xué)I

引用回帖:

2樓: Originally posted by oupengfei717 at 2013-08-17 10:31:52
樓主這個(gè)是靜態(tài)計(jì)算嗎？怎么只有一步離子迭代計(jì)算呢？我猜想可能是該結(jié)構(gòu)在初期全弛豫的時(shí)候沒(méi)有能找到平衡結(jié)構(gòu)？不知道樓主是否進(jìn)行了弛豫得到了平衡結(jié)構(gòu)？樓主的輸入文件是否都檢查過(guò)不存在問(wèn)題呢？

我打算做幾何優(yōu)化的，第一步電子迭代就發(fā)散了，之所以迭代25步，是因?yàn)槲易鰷y(cè)試設(shè)定的NELM =25。

我是輸入文件

SYSTEM = LuMnO3 calculation
ISTART = 0
ICHARG = 2
ISPIN = 2
LREAL = F
#PREC= Accurate
#ADDGRID=.TRUE.

Electronic Relaxation
ENCUT = 500
EDIFF = 1E-07
GGA=PS
NELM =40
NELMDL=8
TIME = 0.1

Ionic relaxation
NSW= 2
IBRION = 1
ISIF = 3
EDIFFG = -0.001
MAGMOM= 6*0 4 4 4 -4 -4 -4 18*0

LDAU= .TRUE.
LDAUTYPE= 2
LDAUL= -1  2  -1
LDAUU=  0  4  0
LDAUJ=  0  1  0
LMAXMIX= 4

AMIX=0.1
BMIX=0.0001
AMIX_MAG=0.4
BMIX_MAG=0.0001

DOS related parameters
ISMEAR = 0
SIGMA = 0.05

OUTPUT OPTION
LWAVE =  True
LCHARG = True

贊一下

回復(fù)此樓

3樓2013-08-17 11:18:27

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

xiaoluoj

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

應(yīng)助: 12 (小學(xué)生)
金幣: 690.1
散金: 5
紅花: 17
帖子: 487
在線: 504.1小時(shí)
蟲(chóng)號(hào): 1526088
注冊(cè): 2011-12-06
專業(yè): 通信理論與系統(tǒng)

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

感謝參與，應(yīng)助指數(shù) +1

NELM確定的是一個(gè)離子步中電子步迭代的最大步數(shù)，25設(shè)定的太小了

贊一下

回復(fù)此樓

4樓2013-08-17 12:47:07

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

whw19850730

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

應(yīng)助: 1 (幼兒園)
金幣: 5875.8
紅花: 1
帖子: 64
在線: 121.3小時(shí)
蟲(chóng)號(hào): 1048864
注冊(cè): 2010-06-28
專業(yè): 物理學(xué)I

引用回帖:

4樓: Originally posted by xiaoluoj at 2013-08-17 12:47:07
NELM確定的是一個(gè)離子步中電子步迭代的最大步數(shù)，25設(shè)定的太小了

25步已經(jīng)發(fā)現(xiàn)-0.30811E+06了在增加迭代步數(shù)沒(méi)有作用，我已經(jīng)試過(guò)了，最后會(huì)越來(lái)越大。

贊一下

回復(fù)此樓

5樓2013-08-18 02:28:01

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

oupengfei717

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

應(yīng)助: 3 (幼兒園)
金幣: 744.8
散金: 80
紅花: 3
帖子: 258
在線: 144小時(shí)
蟲(chóng)號(hào): 1435139
注冊(cè): 2011-10-10
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

樓主將IBRION改為2試試，感覺(jué)LuMnO3的晶體結(jié)構(gòu)應(yīng)該比較復(fù)雜吧？是多少個(gè)原子的晶胞？如果開(kāi)始的晶體結(jié)構(gòu)不是很確定的話，建議IBRION改為2。另外，你的EDIFFG，EDIFF和SIGMA值都設(shè)置的太小了，建議將這些參數(shù)的數(shù)值增大，比如EDIFFG=-0.02，EDIFF=1E-6，SIGMA=0.2應(yīng)該已經(jīng)足以，有些參數(shù)我在計(jì)算的時(shí)候也不常用，所以沒(méi)有太多的經(jīng)驗(yàn)。
如果這樣還是不能收斂的話，第一步試著將ISIF改為7試試，只弛豫晶胞體積，然后第二步改為5，弛豫晶胞形狀，最后再全部弛豫找到平衡結(jié)構(gòu)。

贊一下

回復(fù)此樓

6樓2013-08-18 11:10:56

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

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

[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 336求調(diào)劑 +4	收到VS 2026-03-20	4/200	2026-03-23 19:02 by macy2011
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +3	我愛(ài)學(xué)電池 2026-03-23	3/150	2026-03-23 17:16 by AZMK
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑 +3	Grand777 2026-03-21	3/150	2026-03-21 19:23 by 簡(jiǎn)之-
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨