24小时热门版块排行榜    

查看: 2281  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

子虚乌有5388

木虫 (正式写手)

[交流] 【其他】关于GW(准粒子计算能带结构等)方法的工作 已有6人参与

想做关于准粒子计算能带,激发态,光谱等的工作,大家讨论下,这个算起来是不是很费劲?目前能算准粒子的软件可算多大体系?
另外这方面工作少的原因是不是和计算很耗时间有关系?

谢谢大家指点讨论!
回复此楼
房无一间,地无一亩,身无分文,胸无大志。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主


子虚乌有5388(金币+1): 2010-08-03 22:28:06
xiaohunhun(金币+1):xiexie 2010-08-05 20:57:30
体系稍大点都算不动!
6楼2010-08-03 21:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

hedaors

禁虫 (正式写手)

隐清堂主


子虚乌有5388(金币+1):我打算计算光吸收谱 要用到这个 看来只能罢手了 2010-08-02 23:25:50
qasd(金币+1):谢谢 2010-08-03 07:25:41
这个非常耗时。

并不是所有工作都必须要用 gw 近似的,由于计算量太大,所以一般非必须的话都不会算。
2楼2010-08-02 23:24:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sweetwind88

金虫 (著名写手)


子虚乌有5388(金币+1): 2010-08-03 12:45:35
xiaohunhun(金币+1):xiexie 2010-08-05 20:58:48
是非常耗时的事情。
change
3楼2010-08-03 00:09:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomagh

铁杆木虫 (正式写手)

★ ★
子虚乌有5388(金币+1): 2010-08-08 18:26:05
zzy870720z(金币+2):谢谢分享 2010-08-08 18:47:40
我看到过一个计算的例子 , 贴给大家看看。

-------------------------------------------------------------------------
GW Band Structure Plots
To obtain band structure plots from the GW method, first you have to perform a
DFT calculation. The INCAR file I used for DFT on Si looks like this:
ISTART = 0
ISYM = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
LORBIT = 11
NSW = 200
IBRION = 2
ISIF = 3
EDIFF = 0.0001
EDIFFG = -0.01

KPOINTS file for the run is:
Monkhorst Pack
0
Gamma
4 4 4
0 0 0

POSCAR file is:
Si
5.38936000000000
0.5000000000000000 0.5000000000000000 0.0000000000000000
0.0000000000000000 0.5000000000000000 0.5000000000000000
0.5000000000000000 0.0000000000000000 0.5000000000000000
2
Direct
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.2500000000000000 0.2500000000000000 0.2500000000000000
POTCAR file is supplied with VASP 5.2, remember to use paw pseudopotentials so that in the band structure calculation the Wigner-Seitz radius for Si need not be directly specified in INCAR.

After this run completes you are ready to perform the HF-type (HSE06)
calculation. Change the INCAR file to:
ISTART = 1
ICHARG = 2
EDIFF = 0.0001
EDIFFG = -0.01
ENCUT = 400
ENAUG = 800
LREAL = .FALSE.
LWAVE = .TRUE.
LCHARG = .TRUE.
NELM = 200
NSW = 0
IBRION = -1
LMAXMIX = 4
ISMEAR = 0
SIGMA = 0.1
NSIM = 4
IALGO = 48
ISYM = 0
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = Damped
TIME = 0.4
ENCUTFOCK = 0
NKRED = 2
AEXX = 0.25

Keep the other files constant for this run (e.g. KPOINTS, POSCAR, POTCAR need not be changed).
Once this run has completed, add to the bottom of the INCAR file from the HSE06 calculation:
LOPTICS = .TRUE.
This step calculates the dielectric matrix Once this run has completed, change ALGO = Damped to ALGO = Gwo and
ICHARG = 2 to ICHARG = 11 and add:
NOMEGA = 50
ENCUTGW = 100
LORBIT = 11
and change the KPOINTS file to:
k-points along high symmetry lines
20 ! 20 intersections
Line-mode
rec
0 0 0 ! gamma
0.5 0.0 0 ! X
0.5 0.0 0 ! X
0.5 0.5 0.0 ! M
0.5 0.5 0.0 ! M
0.0 0.0 0.0 ! gamma
0.0 0.0 0.0 ! gamma
0.5 0.5 0.0 ! M
0.5 0.5 0.0 ! M
0.5 0.5 0.5 ! R
0.5 0.5 0.5 ! R
0.0 0.0 0.0 ! gamma
These values are the points of high symmetry along the lines of high symmetry
within the first Brillouin Zone of the Si unit cell structure in reciprocal space. And that’s it! Once complete you can plot the PROCAR, EIGENVAL, and DOSCAR files from the output of this run as described in the previous posts. The band structure plot produced from the plotting utility is shown below.

---------------------------------------------------------------
7楼2010-08-08 12:35:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见