| 查看: 1021 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
孤立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 |
回复此楼» 猜你喜欢
286求调剂
已经有7人回复
-
290求调剂
已经有5人回复
-
085600，320分求调剂
已经有6人回复
-
292分，材料与化工，申请调剂
已经有15人回复
-
0710生物学求调剂
已经有4人回复
-
一志愿北京科技，085601总分305求调剂
已经有5人回复
-
材料求调剂
已经有4人回复
-
一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业
已经有9人回复
-
材料专硕322分
已经有6人回复
-
070300化学求调剂
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
|
后来我成功的获得的文献的结果,三个峰,修改了 静态计算: 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
