在线答疑:Wien2k的问题 - 第一性原理 - 第9页小木虫论坛-学术科研互动平台

…………

81楼2008-12-17 14:20:30

wax54

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 151.2
帖子: 86
在线: 28.4小时
虫号: 507923
注册: 2008-02-20
性别: GG
专业: 凝聚态物性 II ：电子结构

计算出来的DOS spin up部分与fp-lmto基本一致，自旋向下部分能带似乎整天想fermi能级方向移动了1eV多。导致原来导带底在fermi能级上，现在却变成了价带顶在fermi能级上了。是什么参数会和这个东西有关呢？

82楼2008-12-17 16:25:37

wuli8

荣誉版主 (知名作家)

…………

1ST强帖: 2
应助: 35 (小学生)
贵宾: 12.924
金币: 20188.4
散金: 15888
红花: 88
沙发: 4
帖子: 7840
在线: 1114.6小时
虫号: 465889
注册: 2007-11-23
专业: 物理学I
管辖: 计算模拟

★ ★
wuchenwf(金币+2,VIP+0):xiexie

引用回帖:

Originally posted by wax54 at 2008-12-17 16:25:
计算出来的DOS spin up部分与fp-lmto基本一致，自旋向下部分能带似乎整天想fermi能级方向移动了1eV多。导致原来导带底在fermi能级上，现在却变成了价带顶在fermi能级上了。是什么参数会和这个东西有关呢？

这个很正常。
和你的结构有关

赞一下(10人)

…………

83楼2008-12-17 16:49:00

wax54

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 151.2
帖子: 86
在线: 28.4小时
虫号: 507923
注册: 2008-02-20
性别: GG
专业: 凝聚态物性 II ：电子结构

我没明白，我是说我用两个软件算的同一种结构。并且，文献中和我用其他软件算的结果都是导带底在fermi能级上，价带顶在fermi能级以下，他们之间有大约1.5eV的能隙。现在用wien算出来的，spin dn部分相当于整体向高等级方向移动了大约1.5eV，spin up部分基本不变。

84楼2008-12-17 16:57:26

wuli8

荣誉版主 (知名作家)

…………

1ST强帖: 2
应助: 35 (小学生)
贵宾: 12.924
金币: 20188.4
散金: 15888
红花: 88
沙发: 4
帖子: 7840
在线: 1114.6小时
虫号: 465889
注册: 2007-11-23
专业: 物理学I
管辖: 计算模拟

引用回帖:

Originally posted by wax54 at 2008-12-17 16:57:
我没明白，我是说我用两个软件算的同一种结构。并且，文献中和我用其他软件算的结果都是导带底在fermi能级上，价带顶在fermi能级以下，他们之间有大约1.5eV的能隙。现在用wien算出来的，spin dn部分相当于整体向高 ...

能不能把两个图贴出来看看？

…………

85楼2008-12-17 17:07:34

wax54

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 151.2
帖子: 86
在线: 28.4小时
虫号: 507923
注册: 2008-02-20
性别: GG
专业: 凝聚态物性 II ：电子结构

附件1 wien2k GGA NSMO

附件2 lmtart GGA NSMO

还没学会怎么把wien2k的图导出来。看起来可能比较费劲。wien2k中我只给出了TDOS，因为版本比较老，要一个一个画。

lmtart中我给出了TDOS和PDOS

烦劳帮我分析一下，我还是怀疑是哪个地方的参数没搞定。

86楼2008-12-17 20:16:58

wax54

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 151.2
帖子: 86
在线: 28.4小时
虫号: 507923
注册: 2008-02-20
性别: GG
专业: 凝聚态物性 II ：电子结构

还有必要说明下，附件2中，实线表示自旋向上部分，点线表示自旋向下部分。
烦劳wuli8了。

87楼2008-12-17 20:19:52

wuli8

荣誉版主 (知名作家)

…………

1ST强帖: 2
应助: 35 (小学生)
贵宾: 12.924
金币: 20188.4
散金: 15888
红花: 88
沙发: 4
帖子: 7840
在线: 1114.6小时
虫号: 465889
注册: 2007-11-23
专业: 物理学I
管辖: 计算模拟

★ ★ ★
wuchenwf(金币+3,VIP+0):xiexie

赞一下(10人)

…………

88楼2008-12-17 20:48:21

wuli8

荣誉版主 (知名作家)

…………

1ST强帖: 2
应助: 35 (小学生)
贵宾: 12.924
金币: 20188.4
散金: 15888
红花: 88
沙发: 4
帖子: 7840
在线: 1114.6小时
虫号: 465889
注册: 2007-11-23
专业: 物理学I
管辖: 计算模拟

★ ★ ★
wuchenwf(金币+3,VIP+0):xiexie

你的结果和lmtart总体上是一致的体系整体均呈现金属性。
自旋向下的gap大小也是差不多的，只不过是gap出现的位置不一样。
你可以考虑尝试计算其他体系使用相近似的参数，看看是否还有这样的问题。

赞一下(10人)

…………

89楼2008-12-17 20:53:21

wax54

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 151.2
帖子: 86
在线: 28.4小时
虫号: 507923
注册: 2008-02-20
性别: GG
专业: 凝聚态物性 II ：电子结构

恩，总的来说都是一种半金属的能带结构。自旋向上部分无论是TDOS还是PDOS波形基本一致，只是自旋向下的部分有这个差别，貌似产生了一定程度的位移。lmtart程序（FPLMTO）计算出的DOS经过检验是没有问题的，与文献中的基本一致，我还用那个发了个小东西。现在看来可能还是wien2k在画DOS上存在一些问题（可能还没有掌握正确的使用方法）。

90楼2008-12-17 21:13:54