| 查看: 2305 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料专硕322
已经有4人回复
-
308求调剂
已经有11人回复
-
机械专硕274求调剂，不挑专业学校
已经有7人回复
-
336材料与化工085600求调剂
已经有6人回复
-
一志愿郑州大学085600求调剂
已经有18人回复
-
302分求调剂 一志愿安徽大学085601
已经有8人回复
-
294求调剂
已经有5人回复
-
生物与医药求调剂
已经有6人回复
-
求调剂
已经有8人回复
-
一志愿北京化工085600 310分求调剂
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Click 反应的问题
已经有5人回复
- 提取DOS数据出现以下问题
已经有10人回复
- DOSCAR画态密度
已经有16人回复
- 【求助】怎么看分波态密度呢
已经有5人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
长江大学化学与环境工程学院2026年硕士研究生调剂公告（理学、工学均可）
+5/845
-
谷岩课题组招收调剂同学
+1/187
-
双一流南京医科大学招计算机、生物信息等方向研究生调剂
+1/91
-
欢迎调剂北京石油化工学院聚合物功能材料组
+1/84
-
2026年宁波工程学院杜俊毅副研究员招收硕士研究生（含调剂）
+1/49
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
+2/46
-
材料、化学类硕士研究生最后调剂机会，学硕专硕均可，有意请联系，非诚勿扰！
+2/44
-
2026年河北大学无机化学专业招收调剂生
+1/40
-
企业出国访问交流求助
+1/37
-
湖南工大李广利教授课题组接受化学、材料、生物、生医、计算机相关专业调剂生8-10名
+1/17
-
长江师范学院化学化工学院化学工程方向接收大量调剂考生
+1/11
-
湖北省重点大学，化工领域
+1/9
-
【博士招生】天津理工大学国家杰青王铁课题组招收2026年博士研究生
+1/8
-
广东工业大学机电工程学院高伟教授招收博士生
+1/8
-
武汉纺织大学招收化学化工类和材料类研究生调剂生若干名
+1/8
-
渤海大学储能材料课题组招收调剂研究生
+1/8
-
一志愿211电子信息347求调剂
+1/8
-
哈尔滨工程大学动力学院赵建辉团队招收2026博士研究生
+1/3
-
招收【★计算化学、AI药物设计★】方向的硕士研究生（代码：1055/1007/0703）
+1/3
-
武汉工程大学化学与环境工程学院 2026 年硕士研究生调剂公告
+1/2
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
