| 查看: 2548 | 回复: 14 | |||
[交流]
【求助成功】split_dos分割DOSCAR的问题
|
|
各位大侠，有熟悉split_dos的帮帮忙啊，我看了一下split_dos的程序，对于DOS1、DOS2。。。貌似输出的是ene S Py Pz ，侯老师的教程里说第一列为能量，第二三四列是s、p、d态的分波态值。小弟糊涂了。。。下面是部分程序： Atomic projected DOS if [ $nspin -eq 2 ]; then sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, -1*$3, $4, -1*$5, $6, -1*$7}' >> DOS$i else sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, $3, $4 }' >> DOS$i fi else 哪位有好用的split_dos，请给我链接或传我一份，多谢！！！ [ Last edited by wujiandongde on 2011-3-24 at 08:56 ] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
怎么查啊
已经有6人回复
-
基金系统什么内容也没有
已经有8人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Click 反应的问题
已经有5人回复
- 提取DOS数据出现以下问题
已经有10人回复
- DOSCAR画态密度
已经有16人回复
- 【求助】怎么看分波态密度呢
已经有5人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
2楼2011-03-23 16:46:03
3楼2011-03-23 17:00:20
4楼2011-03-23 17:49:44
★ ★
wujiandongde(金币+2): 2011-03-23 19:10:34
zhang668(金币+2): 多谢指导 2011-03-23 20:34:34
wujiandongde(金币+2): 2011-03-23 19:10:34
zhang668(金币+2): 多谢指导 2011-03-23 20:34:34
|
你应该仔细看看源代码和对应的文件格式。原本的split_dos（我不知道你的从哪里来的，原版的是从vtstscript来的）是没有办法处理Px Py Pz的，只有 s p d f. $2,$3就是s, $4,$5就是p， $6,$7就是d 你可以在版面搜索split_dos，我以前修改过一个版本可以处理分波的 [ Last edited by bingmou on 2011-3-23 at 20:14 ] |
5楼2011-03-23 18:27:13
★
zzy870720z(金币+1): 谢谢提示 2011-03-24 08:22:43
zzy870720z(金币+1): 谢谢提示 2011-03-24 08:22:43
6楼2011-03-23 18:33:15
7楼2011-03-23 19:12:20
8楼2011-03-23 21:06:44
9楼2011-03-23 21:46:12
10楼2011-03-23 21:48:05
|
您好！最近计算DOS又遇到了问题，想请您帮忙了，计算DOS时出现了警告WARNING: stress and forces are not correct 我的INCAR如下： SYSTEM = PREC = Normal ISTART=1 ICHARG=11 ENMAX = 500 ! maximum of ENMAX for B and N ALGO = Fast ! electronic relaxation ISMEAR = -5 ; SIGMA = 0.2 ! Gaussian smearing, -5 for tetrahedron method IBRION=-1 ISPIN = 2 !magnetic system LDAU = .TURE. LDATYPE = 2 LDAUL =-1 2 -1 -1 -1 LDAUU = 0 3 0 0 0 LDAUJ = 0 1 0 0 0 LDAUPRINT = 0 LREAL = Auto ! real space projection for a large system ! remove this for an accurate calculation LORBIT=11 EMIN=-15 EMAX=15 NEDOS=3000 不知道哪里设置的有问题，多谢！ |
11楼2011-04-11 12:18:20
12楼2011-04-11 17:55:15
13楼2011-04-11 19:07:26
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
franch(金币+2): 谢谢指点 2011-04-11 20:19:38
小木虫(金币+0.5):给个红包，谢谢回帖交流
franch(金币+2): 谢谢指点 2011-04-11 20:19:38
|
为什么我搜第一个就给出这个呢： It makes no sense to keep the charge density constant (ICHARG-11) and relax the ions (ISIF=2, IBRION=1) at one time. Actually, as you set NSW=0, the positions are not updated anyway, and hence you can ignore that warning. In your case I would prepare 2 separate INCAR files, 1) for the ionic relaxations (ICHARG not set) 2) for the post-processing after the structure is relaxed (ICHARG=11, NSW=0, but IBRION and ISIF not set ) from：http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?4.714 |
14楼2011-04-11 20:03:41
15楼2011-04-11 21:19:38
投票:
30