| 查看: 487 | 回复: 1 | ||
mousekingadv金虫 (小有名气)
|
[求助]
为什么已经收敛的结构，打开LBERRY= .TURE.后不收敛
|
|
如题，已经完全计算收敛的结构，使用其CONTCAR作为POSCAR，仅增加以下参数计算： 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 而原本是收敛的结构。 两边参数完全相同，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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 被积函数很快收敛到0，为什么积分不收敛。
已经有3人回复
- siesta结构优化收敛标准
已经有5人回复
- fluent中continuity曲线不收敛
已经有8人回复
- siesta-3.1 计算石墨烯条带不收敛，不知是不是我的输入有问题。
已经有24人回复
- 用VASP优化结构收敛不了，求指教
已经有10人回复
- 求助，对于未relax的结构测试ecutwfc时不收敛
已经有13人回复
- 键结构没有收敛是什么原因呢？
已经有4人回复
- 用Dmol优化收敛的结构，用vasp不收敛？
已经有9人回复
- 计算时候总是不能收敛如何解决？
已经有5人回复
- fluent温度场分析结果不收敛
已经有6人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- 关于ABINIT计算四方结构弹性常数不收敛的问题
已经有5人回复
- scf自洽不收敛的解决办法
已经有7人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助】高斯计算过渡态不收敛，请问如何解决？
已经有8人回复
- 【求助】结构优化不收敛或是收敛太慢，怎么办？
已经有8人回复
- 【求助】计算不收敛怎么办？
已经有4人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
2楼2013-04-03 08:54:18
