版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

tiger11212

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 597.6
帖子: 32
在線: 505.6小時(shí)
蟲號: 232030
注冊: 2006-03-28
性別: GG
專業(yè): 分子模擬

[求助] 孤立NH3分子的DOS計(jì)算與文獻(xiàn)不符

我首先按如下進(jìn)行:
1 優(yōu)化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點(diǎn)在tao點(diǎn)
得到的計(jì)算結(jié)果和文獻(xiàn)非常吻合.
2 靜態(tài)計(jì)算, 將1步的CONTCAR 復(fù)制為POSCAR, INCAR如下:
System =NH3
ISTART = 0 ; ICHARG = 2
EDIFF = 1E-5 ;
ENCUT  = 400
EDIFFG = -0.015
IBRION = -1
NSW = 0
PREC = Accurate
此處我沒有設(shè)置用了隱含的ISMEAR及sigma值等一些參數(shù)
3 計(jì)算DOS, 并且改變K點(diǎn)為 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分解得到總態(tài)密度DOS0,作圖,發(fā)現(xiàn),與文獻(xiàn)相差太大了,文獻(xiàn)有三個尖銳的峰,而我的都在一起,怎么回事?
另外,這是我的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

回復(fù)此樓

» 猜你喜歡

306求調(diào)劑已經(jīng)有5人回復(fù)
284求調(diào)劑已經(jīng)有4人回復(fù)
尋找調(diào)劑已經(jīng)有4人回復(fù)
一志愿中南化學(xué)（0703）總分337求調(diào)劑已經(jīng)有10人回復(fù)
324求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有7人回復(fù)
319求調(diào)劑已經(jīng)有4人回復(fù)
298求調(diào)劑一志愿211 已經(jīng)有3人回復(fù)
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級已過已經(jīng)有8人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tiger11212

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 597.6
帖子: 32
在線: 505.6小時(shí)
蟲號: 232030
注冊: 2006-03-28
性別: GG
專業(yè): 分子模擬

自己頂一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fengangelo

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hakuna

木蟲 (知名作家)

1ST強(qiáng)帖: 3
應(yīng)助: 82 (初中生)
金幣: 8735.9
散金: 158
紅花: 126
帖子: 5599
在線: 1367.2小時(shí)
蟲號: 167622
注冊: 2006-01-15
專業(yè): 考古理論

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+1, 謝謝交流 2013-04-13 07:32:57

凡分子計(jì)算，k點(diǎn)均設(shè)置為gama點(diǎn)，盒子要足夠大

贊一下

回復(fù)此樓

4樓2013-04-12 06:13:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fengangelo

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

5樓2013-04-12 19:54:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tiger11212

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 597.6
帖子: 32
在線: 505.6小時(shí)
蟲號: 232030
注冊: 2006-03-28
性別: GG
專業(yè): 分子模擬

后來我成功的獲得的文獻(xiàn)的結(jié)果,三個峰,修改了
靜態(tài)計(jì)算:
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計(jì)算:
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
但是得到的能量高的兩峰很尖銳,文獻(xiàn)有一些展寬,能量地的倒似乎正常.我現(xiàn)在很疑惑是我靜態(tài)計(jì)算修改起作用,還是dos計(jì)算里面的重要,需要測試.sigma越小,峰寬越窄,然后我把sigma分別改為0.1和0.2都得到相似的結(jié)果.

贊一下

回復(fù)此樓

6樓2013-04-13 23:45:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangkj

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 tiger11212 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9	niko- 2026-03-19	10/500	2026-03-22 16:08 by ColorlessPI
[考研] 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑 +6	2117205181 2026-03-21	7/350	2026-03-22 14:59 by allen-yin
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 293求調(diào)劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl