版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

byin

铜虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 245
散金: 332
红花: 3
帖子: 395
在线: 174.4小时
虫号: 948127
注册: 2010-01-25
性别: GG
专业: 固体力学

[求助] DOSCAR的问题，输出结果很诡异

计算开了自旋，并且计算原子DOS。输出结果如下：
++++++++++++++++++++++++++++++++++
  50  50 1 0
  0.1203034E+02  0.3944000E-09  0.3944000E-09  0.3867000E-08  0.1000000E-15
1.0000000000000000E-004
  CAR
F
   20.32379460 -31.70330242  301    5.58077025    1.00000000
-31.703  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
-31.530  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
……
   20.150  0.0000E+00  0.0000E+00  0.3200E+03  0.3200E+03
   20.324  0.0000E+00  0.0000E+00  0.3200E+03  0.3200E+03

   20.32379460 -31.70330242  301    5.58077025    1.00000000
-31.703  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  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
-31.530  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  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
-31.356  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  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
……
++++++++++++++++++++++++++++++++
开始总的DOS还好，按照
energy up down    int_up    int_down
格式输出。
但是后面每个原子的DOS不知道怎么回事？还请高手相救！

INCAR：
+++++++++++++++++++++++++
SYSTEM=F
#------------------------
PREC=Accurate
ISTART=0
ICHARG=2
ISPIN=2
MAGMOM=19*1 6*-1 20*3 5*1
#------------------------
ENCUT=600
NELM=200
NELMIN=8
EDIFF=1E-5
#------------------------
EDIFFG=-1E-2
IBRION=2
ISIF=2
ISYM=0
POTIM=0.1
#------------------------
LORBIT=12
#------------------------
ISMEAR=0
SIGMA=0.1
#------------------------
+++++++++++++++++++++++++

[ Last edited by byin on 2011-10-27 at 22:34 ]

回复此楼

» 猜你喜欢

1楼 2011-10-27 22:32:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

byin

铜虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 245
散金: 332
红花: 3
帖子: 395
在线: 174.4小时
虫号: 948127
注册: 2010-01-25
性别: GG
专业: 固体力学

引用回帖:

8楼: Originally posted by cenwanglai at 2011-10-31 14:48:47:
分量相加和总的差别多少？

及Px+Py+Px 和P相差多少？

分别给出上旋和下旋的结果。

我不知道这时候的P，只知道分量。怎么比较？

赞一下

回复此楼

11楼2011-11-03 15:43:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答

windring

禁虫 (小有名气)

★
uuv2010(金币+1): 多谢提示 2011-10-28 17:19:52
byin(金币+10): 2011-10-31 10:42:53

本帖内容被屏蔽

2楼2011-10-28 05:32:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 1469.1
散金: 10
红花: 3
帖子: 235
在线: 40.2小时
虫号: 911880
注册: 2009-11-24
性别: MM
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-10-28 17:19:57
byin(金币+10): 2011-10-31 10:42:58

你的procar里面有相应的排列顺序的应该是你设置了LORBIT>10之后会输出的每一个原子的每一个轨道上的态密度

赞一下(1人)

回复此楼

paper我要好多好的的paper

3楼2011-10-28 11:57:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

byin

铜虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 245
散金: 332
红花: 3
帖子: 395
在线: 174.4小时
虫号: 948127
注册: 2010-01-25
性别: GG
专业: 固体力学

引用回帖:

2楼: Originally posted by windring at 2011-10-28 05:32:49:
可以推断，LORBIT=11.
那么对于每个原子给出partial dos，分别是：
s_up s_dn px_up px_dn py_up py_dn pz_up pz_dn dx2-y2_up dx2-y2_dn ...

两个问题：

1，电子轨道的Z向怎么取？和晶体结构的XYZ什么联系？

2，是不是可以认为：Px+Py+Pz=P，以此类推？

赞一下

回复此楼

4楼2011-10-31 10:53:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答