| 查看: 2303 | 回复: 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 ] |
回复此楼» 猜你喜欢
一志愿北交大材料工程总分358求调剂
已经有4人回复
-
材料专硕322
已经有3人回复
-
336材料与化工085600求调剂
已经有5人回复
-
化学调剂
已经有16人回复
-
320分人工智能调剂
已经有8人回复
-
（调剂）一志愿报考哈尔滨工业大学0857资源与环境专业378分考生
已经有7人回复
-
331求调剂
已经有6人回复
-
材料334求调剂
已经有16人回复
-
308求调剂
已经有9人回复
-
332求调剂
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Click 反应的问题
已经有5人回复
- 提取DOS数据出现以下问题
已经有10人回复
- DOSCAR画态密度
已经有16人回复
- 【求助】怎么看分波态密度呢
已经有5人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
苏州科技大学化生学院蒋老师课题组招收研究生
+1/497
-
铟锡化合物半导体材料与器件团队[推介]
+3/145
-
上海交通大学化学化工学院膜分离-纳流控课题组诚聘博士后2-3名
+1/91
-
江苏大学，环境科学/环境工程专业，调剂
+1/85
-
大连海事大学轮机工程学院尚有少量调剂名额，如有兴趣的同学可以联系。
+2/46
-
2026年河北大学材料与化工专业研究生（专硕）招收调剂生
+1/39
-
河南工业大学化学化工学院 2026 年硕士调剂正在进行！
+1/38
-
哈尔滨工业大学化工学院中俄联合办学硕士生招生
+1/22
-
抗病毒性传染病创新药物全国重点实验室—接收化学、生物、药学、生物与医药调剂考生
+1/22
-
可调声波驱动仪
+1/20
-
大连工业大学 高分子材料 接收考研调剂
+1/19
-
招收调剂
+1/18
-
武汉工程大学化学与环境工程学院2026年硕士研究生调剂公告
+1/15
-
北华大学林产化学加工工程学硕和材料与化工专硕调剂招生
+1/12
-
广东工业大学生态环境与海洋学院张轩教授课题组招收分离膜领域2026年入学博士生1人
+1/9
-
欢迎调剂佛山大学环境工程（专硕）专业研究生
+1/8
-
重庆三峡科技大学光纤传感课题组招调剂硕士研究生（3.30）
+1/6
-
原中山大学池振国教授的五邑大学团队大量接收2026化学化工材料背景考研生调剂
+1/3
-
长江师范学院 材料与化工专业 招收调剂学生
+1/3
-
江西农业大学动科院“昆虫-共生菌协同进化”方向招收调剂研究生
+1/1
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
