计算化学公社

标题: 能带计算相关问题 [打印本页]

作者
Author:
小螃蟹    时间: 2019-3-19 11:07
标题: 能带计算相关问题
sobereva大神,我想请教一下有关能带计算的一些问题。我计算了一个具有504个原子的体系,下面图中是我进行能带计算的输入:
(, 下载次数 Times of downloads: 74)        (, 下载次数 Times of downloads: 71)
首先我已经对构型进行了优化,在这基础上图一显示的是我静态自洽的输入,然后,图二是进行能带分析的输入,我加入了NPAR=1这个关键词,因为看手册中提到可以减小作业内存,第一步静态自洽可以计算,但第二步进行能带分析时,作业会跳出,是否是因为我所计算的体系太大导致的,还是说我的输入的参数有错误?

作者
Author:
Aristotler    时间: 2019-3-19 11:12
K点在算能带的时候是什么? 原因可能是内存溢出了,体系太大是一方面,能带计算本身就很消耗内存。NPAR=1没见过这么用的,我建议你改成NPAR=2再试一下,你也可以设置其他的,因为我也不知道你的配置是什么,NPAR和你用了多少核心有关。

作者
Author:
小螃蟹    时间: 2019-3-19 11:21
Aristotler 发表于 2019-3-19 11:12
K点在算能带的时候是什么? 原因可能是内存溢出了,体系太大是一方面,能带计算本身就很消耗内存。NPAR=1没 ...

谢谢你,我的K点设置的是:
k-points along high symmetry lines               
10 ! 10 intersections                           
Line-mode                                       
rec                                             
0      0       0   ! gamma                     
0      0.5     0   ! X                          
                                                
0      0.5     0   ! X                          
0      0.5   0.5   ! R                          
                                                
0      0.5   0.5   ! R                          
0      0     0.5   ! Z                          
                                                
0      0     0.5   ! Z                          
0      0     0     ! gamma                     
                                                
0      0     0   ! gamma                        
0.5    0.5   0   ! M                           
                                                
0.5    0.5   0   ! M                           
0.5    0.5   0.5 ! A                           
                                                
0.5   0.5  0.5   ! A                           
0     0    0.5   ! Z                           
                                                

作者
Author:
Marisa    时间: 2019-3-21 20:53
输出日志可以贴下
作者
Author:
hpfamily    时间: 2019-3-22 09:04
建议用小一点的体系计算
作者
Author:
小螃蟹    时间: 2019-3-22 15:55
hpfamily 发表于 2019-3-22 09:04
建议用小一点的体系计算

嗯 我之前算过小体系,可以正常算出,但因为现在讨论的是大体系,所以就出现了此问题。
作者
Author:
小螃蟹    时间: 2019-3-22 16:28
Marisa 发表于 2019-3-21 20:53
输出日志可以贴下

嗯,我因为在计算能带时作业跳出,OUTCAR文件中显示的内容如下:
maximum and minimum number of plane-waves per node :    484596   484381

maximum number of plane-waves:    484596
maximum index in each direction:
   IXMAX=   33   IYMAX=   43   IZMAX=   78
   IXMIN=  -34   IYMIN=  -44   IZMIN=  -78

WARNING: aliasing errors must be expected set NGX to 136 to avoid them
WARNING: aliasing errors must be expected set NGY to 176 to avoid them
WARNING: aliasing errors must be expected set NGZ to 314 to avoid them
aliasing errors are usually negligible using standard VASP settings
and one can safely disregard these warnings

serial   3D FFT for wavefunctions
parallel 3D FFT for charge:
    minimum data exchange during FFTs selected (reduces bandwidth)
还请多多指教

作者
Author:
啦啦黑还黑    时间: 2019-3-24 14:04
报错写的很明白NGX 这些参数需要你设置。。
不要用PREC = medium,设导致你的NGX和其他参数不匹配  那是vasp4时代的产物了,手册早就说明不要用Low medium high.
可以试一下改成PREC = Normal或者Accurate
作者
Author:
小螃蟹    时间: 2019-3-25 09:16
啦啦黑还黑 发表于 2019-3-24 14:04
报错写的很明白NGX 这些参数需要你设置。。
不要用PREC = medium,设导致你的NGX和其他参数不匹配  那是va ...

谢谢你,但是我还有一个问题想请教您。就是错误信息中NGX提示的参数设置我不明白该如何设置?
作者
Author:
啦啦黑还黑    时间: 2019-3-25 11:35
小螃蟹 发表于 2019-3-25 09:16
谢谢你,但是我还有一个问题想请教您。就是错误信息中NGX提示的参数设置我不明白该如何设置?

看看按手册的说法,一般都不用设置。只要通过PREC=N或A,控制这个参数即可。特殊需要要增加网格精度再手动设置。
作者
Author:
小螃蟹    时间: 2019-3-25 17:40
哦 好的 谢谢你




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3