版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (3545)
>文献求助 (298)
>导师招生 (292)
>虫友互识 (272)
>休闲灌水 (94)
>论文投稿 (80)
>硕博家园 (68)
>招聘信息布告栏 (66)
>考博 (62)
>博后之家 (57)
>公派出国 (52)
>催化 (47)
>基金申请 (41)
>论文道贺祈福 (40)
>考研 (37)
>教师之家 (34)

返回列表

tiger11212

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 597.6
帖子: 32
在线: 505.6小时
虫号: 232030
注册: 2006-03-28
性别: GG
专业: 分子模拟

[求助] 孤立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

回复此楼

» 猜你喜欢

求助:我三月中下旬出站，青基依托单位怎么办？已经有9人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2013-04-10 17:23:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tiger11212

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 597.6
帖子: 32
在线: 505.6小时
虫号: 232030
注册: 2006-03-28
性别: GG
专业: 分子模拟

自己顶一下

回复此楼

2楼2013-04-11 10:50:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fengangelo

禁虫 (小有名气)

本帖内容被屏蔽

3楼2013-04-11 22:51:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hakuna

木虫 (知名作家)

1ST强帖: 3
应助: 82 (初中生)
金币: 8735.9
散金: 158
红花: 126
帖子: 5599
在线: 1367.2小时
虫号: 167622
注册: 2006-01-15
专业: 考古理论

【答案】应助回帖

★
感谢参与，应助指数 +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

金虫 (初入文坛)

应助: 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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wangkj

7楼2013-04-14 00:11:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 tiger11212 的主题更新

返回列表