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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 孤立NH3分子的DOS计算与文献不符
71/1返回列表
查看: 1323  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2013-04-10 17:23:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiger11212

金虫 (初入文坛)
自己顶一下
回复此楼
2楼2013-04-11 10:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengangelo

禁虫 (小有名气)
本帖内容被屏蔽
3楼2013-04-11 22:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2013-04-13 07:32:57
凡分子计算，k点均设置为gama点，盒子要足够大
赞一下 回复此楼
4楼2013-04-12 06:13:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengangelo

禁虫 (小有名气)
本帖内容被屏蔽
5楼2013-04-12 19:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiger11212

金虫 (初入文坛)
后来我成功的获得的文献的结果,三个峰,修改了
静态计算:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangkj
7楼2013-04-14 00:11:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tiger11212 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 03:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 20/1000 2026-11-01 03:33 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:28 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 00:54 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-10-31 23:18 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭