| 查看: 1029 | 回复: 6 | ||
tiger11212金虫 (初入文坛)
|
[求助]
孤立NH3分子的DOS计算与文献不符
|
|
我首先按如下进行: 1 优化NH3,INCAR如下: general: System = NH3 PREC=Accurate ISTART = 0 ; ICHARG = 2 ; EDIFF = 1E-5 ; ENCUT = 400 ISMEAR = 1; SIGMA = 0.2 Ion Relaxation: EDIFFG = -0.015 NSW = 500 IBRION = 2 POTIM = 0.2 NPAR = 2 K点在tao点 得到的计算结果和文献非常吻合. 2 静态计算, 将1步的CONTCAR 复制为POSCAR, INCAR如下: System =NH3 ISTART = 0 ; ICHARG = 2 EDIFF = 1E-5 ; ENCUT = 400 EDIFFG = -0.015 IBRION = -1 NSW = 0 PREC = Accurate 此处我没有设置用了隐含的ISMEAR及sigma值等一些参数 3 计算DOS, 并且改变K点为 5 * 5 * 5 INCAR如下: System =NH3 ISTART = 1 ; ICHARG = 11 ISMEAR = -5; LORBIT = 11 EDIFF = 1E-5 ; ENCUT = 400 EDIFFG = -0.015 IBRION = -1 NSW = 0 PREC = Accurate 然后用split_dos分解得到总态密度DOS0,作图,发现,与文献相差太大了,文献有三个尖锐的峰,而我的都在一起,怎么回事? 另外,这是我的POSCAR: NH3 1.00000000000000 15.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 15.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 15.0000000000000000 N H 1 3 Selective dynamics Direct 0.5229799999999969 0.4556699999999978 0.4820899999999995 F F F 0.4747572401218676 0.4089732894885714 0.4943562314884602 T T T 0.5828934127766203 0.4248546175374699 0.4929506476151957 T T T 0.5166913251309211 0.5031757768430012 0.5306787459577365 T T T 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 谢谢了  Graph1.jpg |
回复此楼» 猜你喜欢
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
-
一志愿北京科技，085601总分305求调剂
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
tiger11212
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 597.6
- 帖子: 32
- 在线: 505.6小时
- 虫号: 232030
- 注册: 2006-03-28
- 性别: GG
- 专业: 分子模拟
2楼2013-04-11 10:50:18
fengangelo
禁虫 (小有名气)
|
本帖内容被屏蔽 |
3楼2013-04-11 22:51:43
hakuna
木虫 (知名作家)
- 1ST强帖: 3
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
4楼2013-04-12 06:13:32
fengangelo
禁虫 (小有名气)
|
本帖内容被屏蔽 |
5楼2013-04-12 19:54:02
tiger11212
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 597.6
- 帖子: 32
- 在线: 505.6小时
- 虫号: 232030
- 注册: 2006-03-28
- 性别: GG
- 专业: 分子模拟
|
后来我成功的获得的文献的结果,三个峰,修改了 静态计算: System =NH3 ISTART = 0 ICHARG = 2 EDIFF = 1E-5 ; ENCUT = 400 EDIFFG = -0.015 ISMEAR = 0 SIGMA = 0.2 LWAVE = .FALSE. LCHARG = .TRUE. IBRION = -1 NSW = 0 POTIM = 0.2 PREC = Accurate DOS计算: System =NH3 ISTART = 0 ; ICHARG = 11 ISMEAR = -5; LORBIT = 11 EDIFF = 1E-5 ; ENCUT = 400 EDIFFG = -0.015 POTIM = 0.2 IBRION = -1 NSW = 0 NPAR = 1 NEDOS = 3500 PREC = Accurate 但是得到的能量高的两峰很尖锐,文献有一些展宽,能量地的倒似乎正常.我现在很疑惑是我静态计算修改起作用,还是dos计算里面的重要,需要测试.sigma越小,峰宽越窄,然后我把sigma分别改为0.1和0.2都得到相似的结果. |
6楼2013-04-13 23:45:41
 |
7楼2013-04-14 00:11:19
