版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (2818)
>考研 (176)
>導師招生 (168)
>蟲友互識 (93)
>基金申請 (65)
>休閑灌水 (37)
>教師之家 (26)
>文獻求助 (24)
>博后之家 (20)
>碩博家園 (19)
>考博 (16)
>微米和納米 (13)
>論文投稿 (12)
>找工作 (11)
>公派出國 (8)
>論文道賀祈福 (7)

返回列表

mousekingadv

金蟲 (小有名氣)

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

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

如題，已經(jīng)完全計算收斂的結構，使用其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開關時增加了3eV之多（原來是52eV左右）。
同時，  external pressure =    210.78 kB  Pullay stress =       0.00 kB
而原本是收斂的結構。

兩邊參數(shù)完全相同，KPOINTS均為Monkhorst Pack 6*6*6，但vasp分析得出的對稱性卻不一樣，生成的IBZKPT也不一樣，不知道是否和這個有關。
附上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 ]

回復此樓

» 猜你喜歡

山東省面上項目限額評審已經(jīng)有4人回復
一志愿北京化工大學070300 學碩336求調(diào)劑已經(jīng)有5人回復
生物學071000 329分求調(diào)劑已經(jīng)有4人回復
一志愿華中科技大學071000，求調(diào)劑已經(jīng)有4人回復
生物學一志愿985，分數(shù)349求調(diào)劑已經(jīng)有4人回復
生物學調(diào)劑已經(jīng)有4人回復
求調(diào)劑院校信息已經(jīng)有4人回復
085600材料與化工306 已經(jīng)有4人回復
286求調(diào)劑已經(jīng)有10人回復
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有9人回復

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

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

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

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

identation

金蟲 (著名寫手)

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

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

ISYM=-1試試看

回復此樓

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