版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

mousekingadv

金蟲 (小有名氣)

應(yīng)助: 51 (初中生)
金幣: 1141.8
紅花: 1
帖子: 140
在線: 28.7小時
蟲號: 1721002
注冊: 2012-03-27
專業(yè): 理論和計算化學(xué)

[求助] 為什么已經(jīng)收斂的結(jié)構(gòu)，打開LBERRY= .TURE.后不收斂

如題，已經(jīng)完全計算收斂的結(jié)構(gòu)，使用其CONTCAR作為POSCAR，僅增加以下參數(shù)計算：
LBERRY = .TRUE.
IGPAR = 1
NPPSTR = 4
DIPOL = 0.5 0.5 0.5

此后OUTCAR中受力：
POSITION                                     TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
   0.00001    -0.00016    13.42408       0.002270    -0.001368    0.075811
   0.00000    -0.00007    6.60260       0.007010    0.008876    0.065903
   0.00004    -0.00005    3.63661       0.001253    0.013916    0.138615
   0.00004    0.00000    10.45852       -0.006488    -0.010749    0.127437
   1.11037    1.73639    7.23323       0.102562    0.494922    -0.068024
   -2.05901    0.09334    7.23325       -0.489949    -0.155533    -0.070947
   0.94862    -1.82991    7.23327       0.375283    -0.337733    -0.065236
   -0.64638    -1.46858    4.95926       0.494762    -0.153782    -0.062439
   -0.94860    1.29402    4.95925       -0.377880    -0.346448    -0.070688
   1.59499    0.17437    4.95921       -0.108822    0.487900    -0.070432
-----------------------------------------------------------------------------------
total drift:                            -0.007562    -0.010924    -0.080503
同時，free  energy TOTEN  =    -49.840167 eV也比沒有使用LBERRY開關(guān)時增加了3eV之多（原來是52eV左右）。
同時，  external pressure =    210.78 kB  Pullay stress =       0.00 kB
而原本是收斂的結(jié)構(gòu)。

兩邊參數(shù)完全相同，KPOINTS均為Monkhorst Pack 6*6*6，但vasp分析得出的對稱性卻不一樣，生成的IBZKPT也不一樣，不知道是否和這個有關(guān)。
附上POSCAR供參考：
PbNiO
5.43451220000000
   0.4978366681631410 -0.2874349764565433 0.8368208902611921
-0.0000038192337877 0.5748425284672941 0.8368146314927309
-0.4978298375016811 -0.2874311648782443 0.8368270884266549
Pb Ni O
   2    2    6
Direct
  0.9839462595309030  0.9839369543350216  0.9839474441821920
  0.4839490595662141  0.4839479234216597  0.4839530644653465
  0.2665608034278364  0.2665479575295109  0.2665481933520226
  0.7665717645148789  0.7665976760798827  0.7665623625575790
  0.5501051149618448  0.9007352224456184  0.1396860815992866
  0.1396820797655920  0.5501050231659832  0.9007355125994536
  0.9007285823991609  0.1396877102424341  0.5501095422365694
  0.4007286164638728  0.0501122480270979  0.6396502452535557
  0.0501124915364741  0.6396510728241935  0.4007308319153040
  0.6396542146109852  0.4007171987063387  0.0501157086164170

INCAR也貼上吧：
iSYSTEM= TEMP

Start parameter for this run
PREC= Accurate
ENCUT= 400
ALGO= F
LREAL= Auto
GGA= PE
SYMPREC=1E-4
Electronic Relaxation 1
ISTART= 2
ICHARG= 11
ISMEAR= 0
EDIFF= 1E-5
ISPIN= 2
MAGMOM= 2*0 2 -2 6*0

LBERRY = .TRUE.
IGPAR = 1
NPPSTR = 4
DIPOL = 0.5 0.5 0.5

LDAU=T
LDAUL= 2 2 2
LDAUU= 0 4.6 0
LDAUJ= 0 0 0

#DOS related values
#RWIGS= 0.8  1.3
Lorbit= 11
EMIN= -5
EMAX= 15
NEDOS= 2000
#POTIM=0.3
NELM=60

盼望回答！在線等！不勝感激！

[ Last edited by mousekingadv on 2013-3-31 at 17:05 ]

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有6人回復(fù)
招08考數(shù)學(xué) 已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有6人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
285求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有7人回復(fù)
318求調(diào)劑已經(jīng)有6人回復(fù)

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

被積函數(shù)很快收斂到0，為什么積分不收斂。已經(jīng)有3人回復(fù)
siesta結(jié)構(gòu)優(yōu)化收斂標準已經(jīng)有5人回復(fù)
fluent中continuity曲線不收斂已經(jīng)有8人回復(fù)
siesta-3.1 計算石墨烯條帶不收斂，不知是不是我的輸入有問題。已經(jīng)有24人回復(fù)
用VASP優(yōu)化結(jié)構(gòu)收斂不了，求指教已經(jīng)有10人回復(fù)
求助，對于未relax的結(jié)構(gòu)測試ecutwfc時不收斂已經(jīng)有13人回復(fù)
鍵結(jié)構(gòu)沒有收斂是什么原因呢？已經(jīng)有4人回復(fù)
用Dmol優(yōu)化收斂的結(jié)構(gòu)，用vasp不收斂？已經(jīng)有9人回復(fù)
計算時候總是不能收斂如何解決？已經(jīng)有5人回復(fù)
fluent溫度場分析結(jié)果不收斂已經(jīng)有6人回復(fù)
求助：vasp結(jié)構(gòu)優(yōu)化不能收斂已經(jīng)有18人回復(fù)
關(guān)于ABINIT計算四方結(jié)構(gòu)彈性常數(shù)不收斂的問題已經(jīng)有5人回復(fù)
scf自洽不收斂的解決辦法已經(jīng)有7人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化不收斂已經(jīng)有9人回復(fù)
【求助】高斯計算過渡態(tài)不收斂，請問如何解決？已經(jīng)有8人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化不收斂或是收斂太慢，怎么辦？已經(jīng)有8人回復(fù)
【求助】計算不收斂怎么辦？已經(jīng)有4人回復(fù)
【求助】請教PWSCF晶體建模和自洽計算，不收斂已經(jīng)有21人回復(fù)

1樓 2013-03-31 17:03:01

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

identation

金蟲 (著名寫手)

1ST強帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
mousekingadv(franch代發(fā)): 金幣+2, 專家考核, 謝謝回帖交流，， 2013-04-08 21:52:57

ISYM=-1試試看

回復(fù)此樓

2樓2013-04-03 08:54:18