版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

tiger11212

金蟲 (初入文壇)

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

[求助] 孤立NH3分子的DOS計算與文獻不符

我首先按如下進行:
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點在tao點
得到的計算結(jié)果和文獻非常吻合.
2 靜態(tài)計算, 將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 計算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分解得到總態(tài)密度DOS0,作圖,發(fā)現(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ù)此樓

» 猜你喜歡

本人考085602 化學工程專碩已經(jīng)有23人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學070300 學碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(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小時
蟲號: 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強帖: 3
應(yīng)助: 82 (初中生)
金幣: 8735.9
散金: 158
紅花: 126
帖子: 5599
在線: 1367.2小時
蟲號: 167622
注冊: 2006-01-15
專業(yè): 考古理論

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

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

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

贊一下

回復(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小時
蟲號: 232030
注冊: 2006-03-28
性別: GG
專業(yè): 分子模擬

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

[考研] 一志愿211，0703化學310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-21 22:21 by peike
[考研] 化學工程321分求調(diào)劑 +18	大米飯！ 2026-03-15	22/1100	2026-03-21 20:20 by HH領(lǐng)袖
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 0703化學調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學員8dgXkO
[考研] 生物學一志愿985，分數(shù)349求調(diào)劑 +3	zxts12 2026-03-21	3/150	2026-03-21 16:34 by 33來了真來了
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 一志愿華中科技大學071000，求調(diào)劑 +3	沿岸有貝殼6 2026-03-21	3/150	2026-03-21 10:35 by 暮云清寒
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 0856調(diào)劑，是學校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 0703化學336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 334求調(diào)劑 +3	志存高遠意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015