应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » GW0计算出错
51/1返回列表
查看: 4278  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

léon_WANG

木虫 (小有名气)

[求助] GW0计算出错 已有1人参与

vasp5.3.5 跑bulk Si GW0 算例；前两步没问题，GW0计算时出现错误：
reading now W0001.tmp
N1 differs from file     32767       372

-----------------------------------------------------------------------------
|                                                                             |
|     EEEEEEE  RRRRRR   RRRRRR   OOOOOOO  RRRRRR      ###     ###     ###     |
|     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
|     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
|     EEEEE    RRRRRR   RRRRRR   O     O  RRRRRR       #       #       #      |
|     E        R   R    R   R    O     O  R   R                               |
|     E        R    R   R    R   O     O  R    R      ###     ###     ###     |
|     EEEEEEE  R     R  R     R  OOOOOOO  R     R     ###     ###     ###     |
|                                                                             |
|      The file WXXXX.tmp for XXXX=     1   can not be read.                  |
|      Rerun the GW calculations with LRPA=  .TRUE., and do not change        |
|      ENCUTGW or NBANDS between runs.                                        |
|                                                                             |
-----------------------------------------------------------------------------
按提示设置也无济于事，请问有知道怎么回事的没？

计算脚本：
$ECHO "  The real work starts here ...\n"
# DFT calculation
cat > KPOINTS.6 << EOF
Automatic
0
Gamma
6 6 6
0 0 0
EOF
cp KPOINTS.6 KPOINTS
#rm WAVECAR* WAVEDER*

cat > POSCAR << EOF
Si
5.4300
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
2
cart
0.00 0.00 0.00
0.25 0.25 0.25
EOF


cat > INCAR.DFT << EOF
#Si_dft_exc
System  = Si
PREC = Normal ; ENCUT = 250.0
ISMEAR = 0 ; SIGMA = 0.01            ! small sigma is required to avoid partial occupancies
KPAR = 6
EDIFF = 1.E-8                        ! required tight tolerance for ground state orbitals
EOF

$ECHO "  running the DFT calculation for WAVECAR of Si ...\n"
cp INCAR.DFT INCAR
#$VASP_COMMAND
check_failure $?
$ECHO " done"
cp OUTCAR OUTCAR.DFT
cp vasprun.xml vasprun.DFT.xml

cat > INCAR.DIAG << EOF
#Si_dft_exc
System  = Si

PREC = Normal ; ENCUT = 250.0

ALGO = EXACT                          ! use exact diagonalization of the Hamiltonian
NELM = 1                              ! since we are already converged stop after first step
ISMEAR = 0 ; SIGMA = 0.01             ! small sigma is required to avoid partialoccupancies
KPAR = 6                              ! number of cores per compute node or sqrt of number of cores
NBANDS = 96
LOPTICS = .TRUE.                      ! to write WAVEDER the differ of wavefunction
LPEAD = .TRUE.
OMEGAMAX = 40
EOF

$ECHO "  running the DFT calculation for WAVECAR of Si ...\n"
cp INCAR.DIAG INCAR
$VASP_COMMAND
check_failure $?
$ECHO " done"
cp OUTCAR OUTCAR.DIAG
cp vasprun.xml vasprun.DIAG.xml
./extract_optics.sh
./extract_chi.sh
mv optics.dat optics.DFT.dat
mv chi.dat chi.DFT.dat

# GW0 calculation
cat > INCAR.GW0 << EOF
#Si_gw
System  = Si
PREC = Normal ; ENCUT = 250.0
ALGO = GW0  
LSPECTRAL = .TRUE.                          ! use the spectral method:This reduces the computational work load by almost a factor NOMEGA
ISMEAR = 0
SIGMA = 0.01
ENCUTGW = 150
NELM = 1
NOMEGA =  48                                ! specifies the number of frequency grid points
OMEGATL = 280                               ! maximum frequency for coarse part of frequency grid)
KPAR = 6
#NBANDSO=4 ; NBANDSV=8                      ! number of bands for BSE around Fermi level
LADDER=.TRUE.
LUSEW=.TRUE.
LRPA =.FALSE.                               ! local field effects on the Hartree level (RPA); the dielectric matrix is written after the lines:  MACROSCOPIC STATIC DIELECTRIC TENSOR
LWANNIER90=.TRUE.
NBANDS = 96
NBANDSGW = 12                               ! integer twice the number of occupied states
LWAVE = .TRUE.
PRECFOCK = Normal

EOF

$ECHO "  running the GW0 calculation for Si ...\n"
cp INCAR.GW0 INCAR
$VASP_COMMAND
check_failure $?
$ECHO " done"
回复此楼

» 猜你喜欢
1楼 2017-03-14 19:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

léon_WANG

木虫 (小有名气)
引用回帖:
2楼: Originally posted by newton3915 at 2017-03-14 22:29:36
去掉LUSEW这个参数，这个参数有问题。
https://cms.mpi.univie.ac.at/vasp-forum/viewtopic.php?t=16224

果然是这问题。LUSEW=T (use screened WXXXX.tmp files instead of (screened) HF kernel), 看资料WXXXX.tmp files 是G0W0计算才产生的吧？那这参数只对后面GW计算有用是吧
赞一下 回复此楼
3楼2017-03-15 10:16:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

newton3915

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
léon_WANG: 金币+20, ★★★★★最佳答案, 问题搞定 2017-03-15 10:17:39
去掉LUSEW这个参数，这个参数有问题。
https://cms.mpi.univie.ac.at/vasp-forum/viewtopic.php?t=16224
赞一下(1人) 回复此楼
2楼2017-03-14 22:29:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newton3915

铁杆木虫 (正式写手)
引用回帖:
3楼: Originally posted by léon_WANG at 2017-03-15 10:16:55
果然是这问题。LUSEW=T (use screened WXXXX.tmp files instead of (screened) HF kernel), 看资料WXXXX.tmp files 是G0W0计算才产生的吧？那这参数只对后面GW计算有用是吧...

目前版本还没用到。
赞一下 回复此楼
4楼2017-03-15 11:11:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭