24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 1212  |  回复: 1
【悬赏金币】回答本帖问题,作者大力2010将赠送您 10 个金币

大力2010

金虫 (著名写手)

[求助] 请教高手,我在用vasp进行结构优化. 优化1周多了,还没结果,怎样解决?已有1人参与

请教高手,我在用vasp进行结构优化. 优化1周多了,还没结果,怎样解决?  我用的18核进行。 体系一共132个原子,我固定了中部的一些原子(X  Y Z三个方向固定), INCAR  和POSCAR如下所示:

                                Global Parameters
ISTART =  1            (Read existing wavefunction; if there)
# ISPIN =  2           (Spin polarised DFT)
# ICHARG =  11         (Non-self-consistent: GGA/LDA band structures)
LREAL  = .FALSE.       (Projection operators: automatic)
# ENCUT  =  400        (Cut-off energy for plane wave basis set, in eV)
PREC   =  Normal       (Precision level)
LWAVE  = .TRUE.        (Write WAVECAR or not)
LCHARG = .TRUE.        (Write CHGCAR or not)
ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)
# LVTOT  = .TRUE.      (Write total electrostatic potential into LOCPOT or not)
# LVHAR  = .TRUE.      (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =             (No. of electrons: charged cells; be careful)
# LPLANE = .TRUE.      (Real space distribution; supercells)
# NPAR   = 4           (Max is no. nodes; don't set for hybrids)
# NWRITE = 2           (Medium-level output)
# KPAR   = 2           (Divides k-grid into separate groups)
# NGX    = 500         (FFT grid mesh density for nice charge/potential plots)
# NGY    = 500         (FFT grid mesh density for nice charge/potential plots)
# NGZ    = 500         (FFT grid mesh density for nice charge/potential plots)

Electronic Relaxation
ISMEAR =  0            (Gaussian smearing; metals:1)
SIGMA  =  0.05         (Smearing value in eV; metals:0.2)
NELM   =  60           (Max electronic SCF steps)
NELMIN =  6            (Min electronic SCF steps)
EDIFF  =  1E-08        (SCF energy convergence; in eV)
# GGA  =  PS           (PBEsol exchange-correlation)

Ionic Relaxation
NSW    =  100          (Max ionic steps)
IBRION =  2            (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
ISIF   =  2            (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)
EDIFFG = -2E-02        (Ionic convergence; eV/AA)
# ISYM =  2            (Symmetry: 0=none; 2=GGA; 3=hybrids)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

感谢参与,应助指数 +1
自洽精度太高了,EDIFF  =  1E-08,我一般就是1E-05
18核心算132原子也是压力太大了,一般来说1个核心1个原子效率可以达到最大化,最高也不要超过1个核心2原子

How to choose the number of cores
Rule of thumb:

1 atom per core = Good

0.5 atom per core = Could work (but bad efficiency and time wasted)

< 0.5 atom per core = Don’t do it

https://www.pdc.kth.se/software/ ... la/index_using.html
路漫漫其修远兮
2楼2021-08-09 22:03:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大力2010 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见