| 查看: 2165 | 回复: 18 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
shala123456铁杆木虫 (知名作家)
|
[求助]
如何在vasp中TDOS费米能不归零 已有3人参与
|
||
|
用split 分DOSCAR后，分出的Tdos费米能级就在0ev处，比如费米能是7.1ev时，那么其他的能量是7.1的相对值。 现在我的问题是：如何在DOS图中Oev处不代表费米能，改为它的实际能量？ 谢谢大家啦！   5555555555.png |
回复此楼» 猜你喜欢
化学工程321分求调剂
已经有22人回复
-
322求调剂
已经有4人回复
-
考研调剂
已经有3人回复
-
326求调剂
已经有8人回复
-
333求调剂
已经有7人回复
-
一志愿东华大学控制学硕320求调剂
已经有3人回复
-
考研化学学硕调剂，一志愿985
已经有7人回复
-
0703化学调剂
已经有15人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有14人回复
-
326求调剂
已经有4人回复
小学英语本
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1280.8
- 红花: 1
- 帖子: 56
- 在线: 168.8小时
- 虫号: 3724327
- 注册: 2015-03-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
用写字板打开split_dos脚本，查找 $efermi 你应该能找到类似这样的几行代码 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 该行代码是由sed和awk两个处理文件的命令组合而成，其中 $1+(-1)*'$efermi' $1:表示指定行第一个变量,即DOSCAR中的第一列能量； $efermi:表示从OUTCAR中读取的费米能； 你把代码中所有$1+(-1)*'$efermi'改为$1即可。 |
7楼2016-05-22 15:15:18
wangpeng5226
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 72.8
- 帖子: 3
- 在线: 3.4小时
- 虫号: 3510868
- 注册: 2014-10-31
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
shala123456: 金币+6, ★★★很有帮助 2016-01-21 15:29:36
感谢参与，应助指数 +1
shala123456: 金币+6, ★★★很有帮助 2016-01-21 15:29:36
|
把数据导入origin中，能量整体加上费米能，然后以新的能量作图 发自小木虫IOS客户端 |
2楼2016-01-21 14:05:42
shala123456
铁杆木虫 (知名作家)
- 应助: 12 (小学生)
- 金币: 3842.2
- 散金: 40
- 红花: 9
- 沙发: 1
- 帖子: 5802
- 在线: 465.1小时
- 虫号: 2983895
- 注册: 2014-02-21
- 专业: 金属材料
3楼2016-01-21 15:29:12
小学英语本
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1280.8
- 红花: 1
- 帖子: 56
- 在线: 168.8小时
- 虫号: 3724327
- 注册: 2015-03-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shala123456: 金币+12, ★★★很有帮助, 感谢你的回复 2016-05-21 17:19:01
shala123456: 金币+6, ★★★★★最佳答案, 只能在这里加点金币了，谢谢你 2016-05-22 15:30:26
shala123456: 金币+12, ★★★很有帮助, 感谢你的回复 2016-05-21 17:19:01
shala123456: 金币+6, ★★★★★最佳答案, 只能在这里加点金币了，谢谢你 2016-05-22 15:30:26
|
费米能级处在0eV，是由于split_dos脚本导致的，VASP本身没有参数使费米能未0eV。 不知道你用的是哪个split_dos脚本，一般在脚本开始的注释处有以下注释： # Modified to also align the Fermi level from # the ground state to E=0.0 具体代码也做了使费米能出于0eV. 你可以把代码改了。 |
5楼2016-05-18 17:04:43
