应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【其他】问，DOSCAR的格式？
81/1返回列表
查看: 1103  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gleerat

木虫 (正式写手)

[交流] 【其他】问，DOSCAR的格式？已有3人参与

DOSCAR中的分态密度的一部分，第一列是能量，后边的九列是什么？
引用回帖:
15.00000000    -10.00000000 5000      5.46520431      1.00000000
    -10.000  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.995  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.990  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.985  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.980  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.975  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.970  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.965  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.960  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.955  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.950  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     -9.945  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
…… …… ……
…… …… ……

说明书上是这样说的：
引用回帖:
If RWIGS or LORBIT (Wigner Seitz radii, see section 6.326.33) is set in the INCAR file, a lm- and site-projected DOS is calculated and also written to the file DOSCAR.
…… ……
Information on the individual spin components is available only for the site projected density of states, which has the format:
energy s-DOS(total) s-DOS(mx) s-DOS(my) s-DOS(mz) p-DOS(total) p-DOS(mx) ...
In this case, the (site projected) total density of states (total) and the (site projected) energy resolved magnetization density in the   (mx),  (my) and  (mz) direction are available.

在下看不懂。这九列（不算能量）怎么回事？
引用回帖:
energy s-DOS(total) s-DOS(mx) s-DOS(my) s-DOS(mz) p-DOS(total) p-DOS(mx) ...

要是按照这个格式来，到就是从s(tot)到d(tot)，d轨道就没有分量了，这样是怎么回事？如果不是按照这个来，那么按照哪个来？谢谢。
我的INCAR：
引用回帖:
System = CaFe2As2 ThCr2Si2 DOS
PREC = Accurate
ISPIN = 1
ISMEAR = -5 #for DOS,tethraedron method with Blochl corrections
SGIMA = 0.02
ISTART = 1 #contiune from wavecar
EDIFF = 10E-5
ICHARG = 11 #non SC-calculation,to obtain egenvalues
NSW = 0 #for DOS, the steps of ions relaxtions
IBRION = -1 #for DOS, ions no update, not moved

LREAL = Auto #for supercell or largecell
LORBIT = 11 #use the PAW
EMIN = -10
EMAX = 15
NEDOS = 5000
ENLMDL = -5 #no update charge for 3steps
NELMIN = 2 #min of two steps
NELM = 60

LPLANE = .TRUE. #for linux cluster, as follow
NPAR = 1 #for site-projected DOS
LSCALU = .FALSE.
NSIM = 4

[ Last edited by gleerat on 2010-6-5 at 15:08 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
明察、慎思、笃行
1楼2010-06-05 14:37:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anionxt

铁杆木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+2):鼓励新虫!呵呵 2010-06-05 20:59:52
顺序应该是
S
Py
Pz
Px
Dxy
Dyz
Dz2
Dxz
Dx2
赞一下(2人) 回复此楼
2楼2010-06-05 15:55:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ioribest

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
帮顶，我也不懂
energy s-DOS(total) s-DOS(mx) s-DOS(my) s-DOS(mz) p-DOS(total) p-DOS(mx) ...
前面的total跟后面的分量到底是什么关系呢，我的结果里不是我想象的叠加关系啊。
赞一下(1人) 回复此楼
3楼2010-06-05 15:55:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ioribest

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
2楼的顺序是哪里看到的啊？有相关资料么？我想看看~
赞一下(1人) 回复此楼
4楼2010-06-05 15:58:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by anionxt at 2010-06-05 15:55:22:
顺序应该是
S
Py
Pz
Px
Dxy
Dyz
Dz2
Dxz
Dx2

这东西的原则是怎么样的？我在《固体与表面》上也看到了像“z2,xy,xz,yz”等等的东西，说的是耦合么？！
麻烦请给出相关参考资料吧，“授之于鱼不如授之于渔”啊！再次谢谢了

[ Last edited by gleerat on 2010-6-5 at 16:19 ]
赞一下 回复此楼
明察、慎思、笃行
5楼2010-06-05 16:17:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anionxt

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by ioribest at 2010-06-05 15:58:52:
2楼的顺序是哪里看到的啊？有相关资料么？我想看看~

顺序可以看看计算完之后的vasprun.xml，和DOSCAR里面数据对比一下就知道了
赞一下(2人) 回复此楼
6楼2010-06-05 21:32:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anionxt

铁杆木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复.!周末愉快! 2010-06-05 22:07:16
引用回帖:
Originally posted by gleerat at 2010-06-05 16:17:15:

这东西的原则是怎么样的？我在《固体与表面》上也看到了像“z2,xy,xz,yz”等等的东西，说的是耦合么？！
麻烦请给出相关参考资料吧，“授之于鱼不如授之于渔”啊！再次谢谢了

[ Last edited by gleerat on ...

我觉得严格的来说，应该是不同角动量本征波函数D-2 D-1 D0 D+1 D+2 的组合。看看这里
http://hi.baidu.com/sobereva/blo ... 42c43b5c6008cf.html
赞一下(2人) 回复此楼
7楼2010-06-05 21:44:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by anionxt at 2010-06-05 21:44:24:


我觉得严格的来说，应该是不同角动量本征波函数D-2 D-1 D0 D+1 D+2 的组合。看看这里
http://hi.baidu.com/sobereva/blo ... 42c43b5c6008cf.html

好的，谢谢了。学习学习。
赞一下 回复此楼
明察、慎思、笃行
8楼2010-06-06 10:25:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭