 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1102 | 回復(fù): 3 | ||
zhangguangping木蟲(chóng) (著名寫(xiě)手)
|
[求助]
VASP WAVECAR初始化加速收斂的問(wèn)題
|
|
我在算一個(gè)金屬的slab體系的電子結(jié)構(gòu)，想做偶極矯正。如果直接進(jìn)行偶極矯正，算不出來(lái)。如果不做偶極矯正，很快就收斂。 因此，我采用先做一個(gè)沒(méi)有偶極矯正的計(jì)算，得到WAVECAR，然后用這個(gè)WAVECAR去初始化偶極矯正的體系，很快就收斂了。 如果我用偶極矯正之后收斂的WAVECAR作為初始波函數(shù)，為什么就不收斂了？按照我的理解應(yīng)該更快的收斂才對(duì)。因?yàn)檫@個(gè)WAVECAR本身就是對(duì)偶極矯正收斂了的。 事實(shí)是不收斂。不理解，求高手來(lái)解答。 注：我在初始話的時(shí)候，只用了WAVECAR文件。 （1）不做偶極矯正，收斂情況，得到WAVECAR-1 N E dE d eps ncg rms rms(c) RMM: 1 0.215703574514E+05 0.21570E+05 -0.50682E+05 1040 0.103E+03 RMM: 2 0.738803812819E+04 -0.14182E+05 -0.14976E+05 1040 0.266E+02 RMM: 3 0.281239919696E+04 -0.45756E+04 -0.53362E+04 1040 0.141E+02 RMM: 4 0.803056736557E+03 -0.20093E+04 -0.18779E+04 1040 0.955E+01 RMM: 5 -0.222460740543E+02 -0.82530E+03 -0.66398E+03 1040 0.524E+01 RMM: 6 -0.380721447783E+03 -0.35848E+03 -0.28158E+03 1040 0.334E+01 RMM: 7 -0.530123352098E+03 -0.14940E+03 -0.12798E+03 1040 0.214E+01 RMM: 8 -0.605434301842E+03 -0.75311E+02 -0.60532E+02 1040 0.148E+01 RMM: 9 -0.669664938241E+03 -0.64231E+02 -0.58458E+02 2766 0.995E+00 RMM: 10 -0.673551569208E+03 -0.38866E+01 -0.28471E+01 3158 0.239E+00 RMM: 11 -0.673685124635E+03 -0.13356E+00 -0.11441E+00 3164 0.501E-01 RMM: 12 -0.673697301134E+03 -0.12176E-01 -0.11015E-01 3207 0.153E-01 0.926E+01 RMM: 13 -0.636115298836E+03 0.37582E+02 -0.15835E+02 2279 0.119E+01 0.655E+01 RMM: 14 -0.683509423723E+03 -0.47394E+02 -0.10110E+02 2358 0.556E+00 0.549E+01 RMM: 15 -0.668155052094E+03 0.15354E+02 -0.16404E+02 2338 0.847E+00 0.417E+01 RMM: 16 -0.630484188137E+03 0.37671E+02 -0.28445E+01 2259 0.405E+00 0.215E+01 RMM: 17 -0.685947392736E+03 -0.55463E+02 -0.96447E+00 2315 0.210E+00 0.320E+01 RMM: 18 -0.621473850574E+03 0.64474E+02 -0.52893E+00 2261 0.167E+00 0.108E+01 RMM: 19 -0.618837841720E+03 0.26360E+01 -0.17932E+00 2333 0.102E+00 0.702E+00 RMM: 20 -0.618173509222E+03 0.66433E+00 -0.47228E-01 2343 0.467E-01 0.540E+00 RMM: 21 -0.617482369950E+03 0.69114E+00 -0.83465E-01 2298 0.642E-01 0.262E+00 RMM: 22 -0.617425133468E+03 0.57236E-01 -0.21972E-01 2341 0.345E-01 0.219E+00 RMM: 23 -0.617532998741E+03 -0.10787E+00 -0.55507E-02 2316 0.168E-01 0.289E+00 RMM: 24 -0.617290090512E+03 0.24291E+00 -0.29929E-02 2314 0.125E-01 0.649E-01 RMM: 25 -0.617277515774E+03 0.12575E-01 -0.85941E-03 2315 0.691E-02 0.327E-01 RMM: 26 -0.617273761862E+03 0.37539E-02 -0.40333E-03 2318 0.402E-02 0.219E-01 RMM: 27 -0.617272786321E+03 0.97554E-03 -0.32609E-03 2277 0.318E-02 0.154E-01 RMM: 28 -0.617273058184E+03 -0.27186E-03 -0.10491E-03 1803 0.186E-02 0.805E-02 RMM: 29 -0.617274746763E+03 -0.16886E-02 -0.53841E-04 1540 0.132E-02 0.969E-02 RMM: 30 -0.617276157491E+03 -0.14107E-02 -0.27372E-04 1459 0.100E-02 0.705E-02 RMM: 31 -0.617276943507E+03 -0.78602E-03 -0.93009E-05 1409 0.583E-03 0.327E-02 RMM: 32 -0.617278362306E+03 -0.14188E-02 -0.13615E-04 1399 0.615E-03 0.197E-02 RMM: 33 -0.617279466711E+03 -0.11044E-02 -0.97225E-05 1305 0.512E-03 0.165E-02 RMM: 34 -0.617280160088E+03 -0.69338E-03 -0.32354E-05 1019 0.340E-03 0.859E-03 RMM: 35 -0.617280713817E+03 -0.55373E-03 -0.19710E-05 911 0.285E-03 0.876E-03 RMM: 36 -0.617281138459E+03 -0.42464E-03 -0.12955E-05 830 0.254E-03 0.618E-03 RMM: 37 -0.617281411872E+03 -0.27341E-03 -0.10039E-05 737 0.228E-03 0.370E-03 RMM: 38 -0.617281541513E+03 -0.12964E-03 -0.53457E-06 570 0.201E-03 0.236E-03 RMM: 39 -0.617281599275E+03 -0.57762E-04 -0.38373E-06 547 0.199E-03 （2）我用WAVECAR-1來(lái)初始化做偶極矯正的體系，得到WAVECAR-2 N E dE d eps ncg rms rms(c) RMM: 1 -0.619344694650E+03 -0.61934E+03 -0.14347E-02 2269 0.661E-02 0.718E+00 RMM: 2 -0.350538948397E+04 -0.28860E+04 -0.11684E+03 3099 0.102E+01 0.178E+02 RMM: 3 -0.801088918231E+03 0.27043E+04 -0.51753E+01 2283 0.404E+00 0.711E+01 RMM: 4 -0.655700042436E+03 0.14539E+03 -0.78366E+00 2305 0.184E+00 0.138E+01 RMM: 5 -0.620810364391E+03 0.34890E+02 -0.18897E+00 2335 0.903E-01 0.744E+00 RMM: 6 -0.617856946183E+03 0.29534E+01 -0.37157E-01 2321 0.453E-01 0.359E+00 RMM: 7 -0.618270265827E+03 -0.41332E+00 -0.98294E-02 2295 0.244E-01 0.466E+00 RMM: 8 -0.617375153783E+03 0.89511E+00 -0.25550E-02 2324 0.115E-01 0.163E+00 RMM: 9 -0.617370621743E+03 0.45320E-02 -0.47104E-03 2332 0.533E-02 0.160E+00 RMM: 10 -0.617337743937E+03 0.32878E-01 -0.69003E-03 2273 0.634E-02 0.127E+00 RMM: 11 -0.617593438930E+03 -0.25569E+00 -0.23585E-02 2233 0.117E-01 0.307E+00 RMM: 12 -0.617315259950E+03 0.27818E+00 -0.19643E-02 2243 0.104E-01 0.107E+00 RMM: 13 -0.617289636662E+03 0.25623E-01 -0.34083E-03 2304 0.421E-02 0.510E-01 RMM: 14 -0.617298933662E+03 -0.92970E-02 -0.16248E-03 2232 0.314E-02 0.777E-01 RMM: 15 -0.617286114518E+03 0.12819E-01 -0.10389E-03 2219 0.243E-02 0.378E-01 RMM: 16 -0.617282658705E+03 0.34558E-02 -0.11504E-03 2169 0.271E-02 0.173E-01 RMM: 17 -0.617282005095E+03 0.65361E-03 -0.21119E-04 1457 0.116E-02 0.127E-01 RMM: 18 -0.617281818120E+03 0.18698E-03 -0.16814E-04 1415 0.100E-02 0.133E-01 RMM: 19 -0.617281470952E+03 0.34717E-03 -0.55315E-05 1356 0.521E-03 0.970E-02 RMM: 20 -0.617281402108E+03 0.68844E-04 -0.19169E-05 925 0.309E-03 （3）我用WAVECAR-2同樣來(lái)初始化（2）計(jì)算中用的體系 N E dE d eps ncg rms rms(c) RMM: 1 -0.635408979780E+03 -0.63541E+03 -0.46663E-01 2270 0.349E-01 0.136E+01 RMM: 2 -0.587150192989E+04 -0.52361E+04 -0.21556E+03 3100 0.167E+01 0.190E+02 RMM: 3 -0.802576284254E+03 0.50689E+04 -0.24057E+02 2284 0.750E+00 0.686E+01 RMM: 4 -0.189846693072E+06 -0.18904E+06 -0.18002E+05 3333 0.164E+03 0.927E+02 RMM: 5 -0.295410388026E+04 0.18689E+06 -0.42425E+04 3490 0.135E+02 0.223E+02 RMM: 6 -0.257017901718E+06 -0.25406E+06 -0.31875E+05 3129 0.223E+03 0.143E+03 RMM: 7 0.309556257787E+04 0.26011E+06 -0.44927E+04 3775 0.179E+02 0.279E+02 RMM: 8 -0.807832907454E+06 -0.81093E+06 -0.23071E+06 3093 0.127E+04 0.196E+03 RMM: 9 0.501436454573E+04 0.81285E+06 -0.14483E+05 4160 0.372E+02 0.882E+02 RMM: 10 -0.139064714325E+06 -0.14408E+06 -0.10275E+06 3429 0.110E+04 0.223E+03 RMM: 11 0.129838979838E+04 0.14036E+06 -0.14362E+05 4160 0.368E+02 0.145E+03 RMM: 12 -0.144032234536E+06 -0.14533E+06 -0.11086E+06 3412 0.261E+03 0.257E+03 RMM: 13 0.881743429912E+04 0.15285E+06 -0.62072E+04 4156 0.210E+02 0.960E+02 RMM: 14 -0.460974531041E+05 -0.54915E+05 -0.22373E+05 3607 0.559E+02 0.217E+03 RMM: 15 -0.668075561127E+06 -0.62198E+06 0.81775E+05 3292 0.137E+04 0.342E+03 RMM: 16 0.219569463762E+05 0.69003E+06 -0.10915E+05 4160 0.343E+02 0.130E+03 RMM: 17 -0.342406208787E+06 -0.36436E+06 -0.72604E+05 3648 0.537E+03 0.191E+03 RMM: 18 0.353335346143E+03 0.34276E+06 -0.18256E+05 4160 0.439E+02 0.163E+03 RMM: 19 0.141289018611E+05 0.13776E+05 -0.57101E+04 3919 0.201E+02 0.807E+02 RMM: 20 0.163815991815E+05 0.22527E+04 -0.57281E+04 3854 0.206E+02 0.873E+02 RMM: 21 0.114491016167E+05 -0.49325E+04 -0.19919E+04 4151 0.816E+01 0.629E+02 RMM: 22 0.105640996311E+05 -0.88500E+03 -0.94818E+03 4146 0.508E+01 0.568E+02 RMM: 23 0.890036399913E+04 -0.16637E+04 -0.96320E+03 4127 0.459E+01 0.569E+02 RMM: 24 0.101107966556E+05 0.12104E+04 -0.10282E+04 4135 0.523E+01 0.535E+02 RMM: 25 0.945038614460E+04 -0.66041E+03 -0.36972E+03 4129 0.334E+01 0.493E+02 RMM: 26 0.911927818592E+04 -0.33111E+03 -0.63625E+03 4143 0.365E+01 0.495E+02 RMM: 27 0.877241016239E+04 -0.34687E+03 -0.33588E+03 4139 0.290E+01 0.470E+02 RMM: 28 0.685759686077E+04 -0.19148E+04 -0.11458E+04 4107 0.476E+01 0.558E+02 RMM: 29 0.753315758343E+04 0.67556E+03 -0.29950E+03 4109 0.263E+01 0.469E+02 RMM: 30 0.734875761516E+04 -0.18440E+03 -0.22085E+03 4083 0.225E+01 0.455E+02 RMM: 31 0.801156500329E+04 0.66281E+03 -0.22746E+03 4042 0.224E+01 0.412E+02 RMM: 32 0.830584883622E+04 0.29428E+03 -0.40099E+03 4140 0.285E+01 0.421E+02 RMM: 33 0.746204161945E+04 -0.84381E+03 -0.10945E+04 4158 0.515E+01 0.517E+02 RMM: 34 0.676862265457E+04 -0.69342E+03 -0.72662E+03 4145 0.420E+01 0.574E+02 RMM: 35 0.593350742310E+04 -0.83512E+03 -0.61020E+03 4160 0.365E+01 0.618E+02 RMM: 36 0.558985295988E+04 -0.34365E+03 -0.45523E+03 4098 0.325E+01 0.657E+02 RMM: 37 0.528138589350E+04 -0.30847E+03 -0.35727E+03 4160 0.282E+01 0.685E+02 RMM: 38 0.487823298447E+04 -0.40315E+03 -0.30134E+03 4105 0.262E+01 0.708E+02 RMM: 39 0.461322301400E+04 -0.26501E+03 -0.24204E+03 4153 0.241E+01 0.727E+02 .... 高手來(lái)解答吧。這個(gè)問(wèn)題的出現(xiàn)肯定是兩次做的偶極矯正不一樣導(dǎo)致的。第一次的偶極矯正得到的收斂波函數(shù)對(duì)第二次矯正來(lái)說(shuō)不是好的初始猜測(cè)了。 |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
g_wcai
木蟲(chóng) (正式寫(xiě)手)
開(kāi)心
- 應(yīng)助: 35 (小學(xué)生)
- 金幣: 4696.6
- 紅花: 8
- 帖子: 460
- 在線: 163.6小時(shí)
- 蟲(chóng)號(hào): 915747
- 注冊(cè): 2009-11-29
- 性別: GG
- 專(zhuān)業(yè): 半導(dǎo)體科學(xué)與信息器件
★
蛋蛋小童鞋: 金幣+1, 謝謝參與 2016-08-01 14:58:49
蛋蛋小童鞋: 金幣+1, 謝謝參與 2016-08-01 14:58:49
|
是不是先搞清楚加偶標(biāo)的實(shí)質(zhì)？！這樣會(huì)有助于理解問(wèn)題，根據(jù)題主的一子描述，我認(rèn)為你是在加了一個(gè)偶極子收斂后再加了一個(gè)偶極子，這樣會(huì)不會(huì)導(dǎo)致不收斂！也有可能是金屬材料的緣故！拋個(gè)磚，大牛接上！ 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
2樓2016-03-31 08:14:42
zhangjing233
金蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 552.8
- 帖子: 88
- 在線: 77.2小時(shí)
- 蟲(chóng)號(hào): 3599487
- 注冊(cè): 2014-12-17
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2016-07-31 16:02:43
4樓2016-07-31 16:20:57
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhangguangping 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 材料專(zhuān)碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|