當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

rustboy

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 52.5
散金: 30
紅花: 1
帖子: 23
在線: 29.5小時(shí)
蟲號(hào): 1483024
注冊(cè): 2011-11-08
專業(yè): 原子和分子物理

[求助] vasp 計(jì)算單分子/原子時(shí)晶胞大小的選取

我覺得，原則上講應(yīng)該是使相鄰晶胞內(nèi)分子/原子距離足夠大，以消除計(jì)算過程中互相的影響。滿足這一條件后，在相同的精度下，所得到的單分子/原子的能量大致應(yīng)該是相當(dāng)?shù)摹?br />
但是本人在計(jì)算單個(gè)金屬Ni時(shí)發(fā)現(xiàn)，使用10*10*10的晶胞以及15*15*15晶胞，其它參數(shù)一樣。計(jì)算所得單原子的能量相差較大。

不知道大家平時(shí)計(jì)算類似的體系時(shí)，是如何設(shè)置晶胞大小的？

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

化學(xué)電源之鋰離子電池

» 猜你喜歡

一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)308分求調(diào)劑已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有11人回復(fù)
335分 | 材料與化工專碩 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有3人回復(fù)
一志愿武理材料工程348求調(diào)劑已經(jīng)有9人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有18人回復(fù)
0703化學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
材料292調(diào)劑已經(jīng)有3人回復(fù)
一志愿陜師大生物學(xué)071000，298分，求調(diào)劑已經(jīng)有3人回復(fù)
306求0703調(diào)劑一志愿華中師范已經(jīng)有7人回復(fù)

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

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

如何用VASP計(jì)算H原子的振動(dòng)頻率已經(jīng)有16人回復(fù)
態(tài)密度計(jì)算 vasp 4.6 計(jì)算Zn原子的PDOS 關(guān)于Zn（3d104s2）的4p對(duì)價(jià)帶的貢獻(xiàn) 已經(jīng)有4人回復(fù)
怎么用vasp的計(jì)算結(jié)果分析分子中各原子帶電多少情況已經(jīng)有6人回復(fù)
畫VASP計(jì)算的原子結(jié)構(gòu)圖的一個(gè)方法（自己瞎想的，如和他人方法雷同請(qǐng)見諒）已經(jīng)有4人回復(fù)
vasp計(jì)算界面結(jié)構(gòu)時(shí)原子弛豫的問題已經(jīng)有12人回復(fù)
vasp計(jì)算單個(gè)氧原子和分子的能量已經(jīng)有20人回復(fù)
VASP 四方晶系（TiO2 anatase ）單胞建模原子數(shù)的問題麻煩大家指導(dǎo) 已經(jīng)有19人回復(fù)
在vasp計(jì)算里，有沒有好的只優(yōu)化原子位置和晶胞體積的方法？已經(jīng)有8人回復(fù)
vasp可以固定單個(gè)原子的磁矩嗎？已經(jīng)有5人回復(fù)
如何用VASP計(jì)算體系中每一類原子的DOS？已經(jīng)有10人回復(fù)
vasp中如何構(gòu)建超晶胞已經(jīng)有3人回復(fù)
求助：vasp 計(jì)算氫原子能量是-1.16ev 而不是 -13.6ev，why ? 已經(jīng)有18人回復(fù)
【求助】如何顯示晶胞結(jié)構(gòu)中的原子編號(hào)與空間坐標(biāo)？已經(jīng)有10人回復(fù)
【求助完畢】vasp-slab-邊界附近原子頂位吸附小分子的結(jié)構(gòu)優(yōu)化很難收斂嗎？已經(jīng)有5人回復(fù)
【求助】vasp為什么采用周期性邊界條件（或超原胞模型）來處理原子分子表面體系等，有已經(jīng)有8人回復(fù)
【求助】能否發(fā)一份用vasp計(jì)算超晶胞某個(gè)表面的輸入文件給我已經(jīng)有7人回復(fù)
【求助】Castep中結(jié)構(gòu)優(yōu)化后晶胞參數(shù)差別較大 vasp高手也幫我看一下已經(jīng)有22人回復(fù)
【求助】關(guān)于鋰原子vasp計(jì)算后的1s2s2p分波態(tài)密度處理已經(jīng)有12人回復(fù)

1樓 2012-12-30 20:02:57

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

rustboy

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 52.5
散金: 30
紅花: 1
帖子: 23
在線: 29.5小時(shí)
蟲號(hào): 1483024
注冊(cè): 2011-11-08
專業(yè): 原子和分子物理

引用回帖:

2樓: Originally posted by emilyoyang at 2012-12-30 20:19:41
這個(gè)能量差別很多可能是你的cutoff 取得過小。。。不知道你取多少

你試著計(jì)算一下9.5*10*10.5 和14.5*15*15.5這兩個(gè)超胞中單個(gè)Ni原子的能量差是不是也很大？也有可能是對(duì)稱性造成的

cutoff我取的是400，具體INCAR參數(shù)請(qǐng)見三樓
我新算了一個(gè)Co原子。
結(jié)果是10*10*10的超胞能量1.480eV，9.5*10*10.5的1.785eV，15*15*15的1.931, 14.5*15*15.5的1.944。貌似沒有規(guī)律性。
這些能量差距，請(qǐng)問前輩是何原因？

另外，本人菜鳥，前輩提到的 “對(duì)稱性造成的” 是什么意思？

贊一下

回復(fù)此樓

5樓2013-01-01 10:51:27

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

查看全部 7 個(gè)回答

emilyoyang

木蟲 (正式寫手)

應(yīng)助: 232 (大學(xué)生)
金幣: 4340.4
紅花: 42
帖子: 667
在線: 189.9小時(shí)
蟲號(hào): 1630291
注冊(cè): 2012-02-20
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

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

★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+1, 謝謝交流 2013-01-03 12:09:43

這個(gè)能量差別很多可能是你的cutoff 取得過小。。。不知道你取多少

你試著計(jì)算一下9.5*10*10.5 和14.5*15*15.5這兩個(gè)超胞中單個(gè)Ni原子的能量差是不是也很大？也有可能是對(duì)稱性造成的

贊一下

回復(fù)此樓

2樓2012-12-30 20:19:41

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

hn19870519

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 45 (小學(xué)生)
金幣: 1225.2
散金: 1861
紅花: 11
帖子: 379
在線: 203.1小時(shí)
蟲號(hào): 867945
注冊(cè): 2009-10-10
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+1, 謝謝提示 2013-01-03 12:09:52

你說的能量相差很大是多大？最好把你計(jì)算的參數(shù)貼出來，這樣大家也好幫你。

贊一下

回復(fù)此樓

3樓2012-12-31 22:17:09

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

rustboy

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 52.5
散金: 30
紅花: 1
帖子: 23
在線: 29.5小時(shí)
蟲號(hào): 1483024
注冊(cè): 2011-11-08
專業(yè): 原子和分子物理

引用回帖:

3樓: Originally posted by hn19870519 at 2012-12-31 22:17:09
你說的能量相差很大是多大？最好把你計(jì)算的參數(shù)貼出來，這樣大家也好幫你。

如下是我的INCAR參數(shù)：
  NWRITE =    2; LPETIM=F write-flag & timer
ISTART =    0 job : 0-new  1-cont  2-samecut
ICHARG =    2 charge: 1-file 2-atom 10-const
INIWAV =    1 electr: 0-lowe 1-rand  2-diag

Electronic Relaxation 1
PREC = Accurate ! standard precision
ENCUT =    400.000
NELM =    500; NELMIN= 2; NELMDL= 0    # of ELM steps
ISPIN  =    2 spin polarized calculation?
EDIFF  =    0.1E-04

Ionic Relaxation
EDIFFG = -0.03
NSW =    500  number of steps for IOM
NBLOCK =    1;
IBRION =    2
ISIF =    2
ISYM =    2 0-nonsym 1-usesym
LCORR  =    T Harris-correction to forces
POTIM  =    0.5 time-step for ion-motion

Electronic Relaxation 2
LPLANE = .TRUE.
ALGO = Fast
LDIAG  = T    sub-space diagonalisation
LWAVE =  F    write WAVECAR
LCHARG  = F    write CHGCAR and CHG
LVDW = .TRUE.
LREAL  = Auto    real-space projection
NSIM = 4
昨天又算了一個(gè)Co原子，情況也是類似，10*10*10的能量是1.480eV,15的是1.93，20的是1.476。
謝前輩幫忙！

贊一下

回復(fù)此樓

4樓2013-01-01 09:51:57

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

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

[考研] 一志愿武理材料工程348求調(diào)劑 +6	￣^￣゜汗 2026-03-19	9/450	2026-03-23 19:53 by pswait
[考研] 336求調(diào)劑 +4	收到VS 2026-03-20	4/200	2026-03-23 19:02 by macy2011
[考研] 07化學(xué)280分求調(diào)劑 +3	722865 2026-03-23	3/150	2026-03-23 18:48 by macy2011
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +4	石能偉 2026-03-17	4/200	2026-03-23 17:51 by 17862566385
[考研] 333求調(diào)劑 +6	87639 2026-03-21	10/500	2026-03-23 10:41 by Iveryant
[考研] 291求調(diào)劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong