| 查看: 4557 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
léon_WANG木虫 (小有名气)
|
[求助]
为什么vasp DFT与GW0计算后的wannier插值能带完全一样 已有1人参与
|
||
|
如题，照官网SrVO3算例弄的Bulk Si能带结构计算。从DFT获得的gap为0.62eV，G0W0获得的gap为1.16eV，但两者的wannier插值能带完全一样。求解。 脚本如下： $ECHO " The real work starts here ...\\n" # DFT calculation $ECHO " step 1: scf calculation ...\\n" cat > KPOINTS.8 << EOF Automatic 0 Gamma 8 8 8 0 0 0 EOF cp KPOINTS.8 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 NBANDS = 8 LWANNIER90_RUN = .TRUE. EDIFF = 1.E-8 ! required tight tolerance for ground state orbitals EOF # Wannier90.win for bandstructure cat > wannier90_NONE.win << EOF num_wann=8 num_bands=8 # for GW uncomment # exclude_bands: 9-24 Begin Projections Si:sp3 End Projections dis_froz_max=9 dis_num_iter=1000 guiding_centres=true # Bandstructure plot bands_plot = true begin kpoint_path L 0.50000 0.50000 0.5000 G 0.00000 0.00000 0.0000 G 0.00000 0.00000 0.0000 X 0.50000 0.00000 0.5000 X 0.50000 0.00000 0.5000 K 0.37500 -0.37500 0.0000 K 0.37500 -0.37500 0.0000 G 0.00000 0.00000 0.0000 end kpoint_path bands_num_points 40 bands_plot_format: gnuplot EOF $ECHO " running the DFT calculation for WAVECAR of Si ...\\n" cp INCAR.DFT INCAR cp wannier90_NONE.win wannier90.win $VASP_COMMAND check_failure $? $ECHO " done" cp OUTCAR OUTCAR.DFT cp vasprun.xml vasprun.DFT.xml cp wannier90_band.dat wannier90_NONE_band.dat $ECHO " step 2: obtain virtual obitals ...\\n" 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. ! Derivative of the orbitals w.r.t. the k-point OMEGAMAX = 10 NEDOS = 2000 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 # GW0 calculation $ECHO " step 3: GW0 calculation ...\\n" 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. ! use screened WXXXX.tmp files instead of (screened) HF kernel 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 ！determines how many QP energies are calculated and updated in GW calculations. integer twice the number of occupied states LWAVE = .TRUE. PRECFOCK = Normal LOPTICS=.TRUE. EOF $ECHO " running the GW0 calculation for le ...\\n" cp INCAR.GW0 INCAR $VASP_COMMAND check_failure $? $ECHO " done" cp OUTCAR OUTCAR.GW0 cp vasprun.xml vasprun.GW0.xml # Wannier90.win for bandstructure cat > wannier90_GW0.win << EOF num_wann=8 num_bands=8 # for GW uncomment exclude_bands: 9-96 Begin Projections Si:sp3 End Projections dis_froz_max=9 dis_num_iter=1000 guiding_centres=true # Bandstructure plot bands_plot = true begin kpoint_path L 0.50000 0.50000 0.5000 G 0.00000 0.00000 0.0000 G 0.00000 0.00000 0.0000 X 0.50000 0.00000 0.5000 X 0.50000 0.00000 0.5000 K 0.37500 -0.37500 0.0000 K 0.37500 -0.37500 0.0000 G 0.00000 0.00000 0.0000 end kpoint_path bands_num_points 40 bands_plot_format: gnuplot EOF cat >INCAR.BANDS << EOF #Si_gw System = Si PREC = Normal ; ENCUT = 250.0 ALGO = Nothing ; NELM = 1 ISMEAR = 0 ; SIGMA = 0.01 KPAR = 6 NBANDS = 96 LWAVE = .FALSE. LOPTICS = .TRUE. ; LPEAD = .TRUE. OMEGAMAX = 10 LWANNIER_RUN =.TRUE. EOF cp wannier90_GW0.win wannier90.win cp INCAR.BANDS INCAR $ECHO " running for bands with GW0 ...\\n" $VASP_COMMAND |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 |
» 猜你喜欢
面上合作单位盖章
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
newton3915
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4742.3
- 散金: 40
- 红花: 4
- 帖子: 418
- 在线: 391.3小时
- 虫号: 609189
- 注册: 2008-09-22
- 专业: 凝聚态物性 II ：电子结构
5楼2017-03-18 10:41:17
newton3915
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4742.3
- 散金: 40
- 红花: 4
- 帖子: 418
- 在线: 391.3小时
- 虫号: 609189
- 注册: 2008-09-22
- 专业: 凝聚态物性 II ：电子结构
2楼2017-03-18 00:26:22
léon_WANG
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 2577.5
- 散金: 20
- 红花: 4
- 帖子: 156
- 在线: 226.1小时
- 虫号: 2458360
- 注册: 2013-05-10
- 专业: 凝聚态物性I:结构、力学和
|
感觉最后一步vasp计算就是第一步重新来一遍，wannier插值用的都是DFT的值而不是中间GW0计算的IP值，官网上这个算例http://cms.mpi.univie.ac.at/wiki ... ture_of_SrVO3_in_GW不知你跑过没有。照理它没理由算两个一样的能带前后对比。难道要在GW0计算中加LWANNIER_RUN =.TRUE.才行 |
3楼2017-03-18 10:09:42
léon_WANG
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 2577.5
- 散金: 20
- 红花: 4
- 帖子: 156
- 在线: 226.1小时
- 虫号: 2458360
- 注册: 2013-05-10
- 专业: 凝聚态物性I:结构、力学和
4楼2017-03-18 10:23:11
