| 查看: 2270 | 回复: 6 | |||
[交流]
【其他】关于GW（准粒子计算能带结构等）方法的工作 已有6人参与
|
|
想做关于准粒子计算能带，激发态，光谱等的工作，大家讨论下，这个算起来是不是很费劲？目前能算准粒子的软件可算多大体系？ 另外这方面工作少的原因是不是和计算很耗时间有关系？ 谢谢大家指点讨论！ |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于女人和工作的问题
已经有27人回复
- HSE（GW类似）计算能带总结
已经有51人回复
- 关于考研和工作，大家给点意见
已经有9人回复
- vasp 中能用 GW算能带吗？
已经有7人回复
- 【求助】VASP计算光学性质时的GW方法（RPA, no local field effects)
已经有5人回复
- 请教大家关于项目申报或项目管理方面的问题
已经有10人回复
- 【求助】pwscf中的GW计算能隙的问题
已经有5人回复
- 关于国外硕博士论文搜索和下载的讨论
已经有81人回复
- 【求助】GW /GW0/G0W0 计算是不是一定要基于HSE 计算的结果？
已经有4人回复
hedaors
禁虫 (正式写手)
隐清堂主
- 应助: 11 (小学生)
- 贵宾: 0.225
- 金币: 3948.8
- 散金: 201
- 红花: 8
- 帖子: 855
- 在线: 271.1小时
- 虫号: 449622
- 注册: 2007-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-08-02 23:24:07
sweetwind88
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 2440.6
- 散金: 104
- 红花: 2
- 帖子: 1258
- 在线: 254.9小时
- 虫号: 1062935
- 注册: 2010-07-23
- 专业: 无机非金属基复合材料
3楼2010-08-03 00:09:22
 |
4楼2010-08-03 01:15:09
5楼2010-08-03 12:45:58
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
6楼2010-08-03 21:04:40
xiaomagh
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 6113.1
- 红花: 5
- 帖子: 652
- 在线: 150.8小时
- 虫号: 224710
- 注册: 2006-03-22
- 性别: GG
- 专业: 金属材料的合金相、相变及
★ ★
子虚乌有5388(金币+1): 2010-08-08 18:26:05
zzy870720z(金币+2):谢谢分享 2010-08-08 18:47:40
子虚乌有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
40