dyc_2008

應(yīng)助: 14 (小學(xué)生)
金幣: 124.2
帖子: 1185
在線: 244.3小時
蟲號: 1044421

[交流] 請大牛指教。關(guān)于截斷動能、原子數(shù)目，K點密度

基于周期性邊界條件的VASP，要用超胞堆疊整個3維空間，這樣的話，原子數(shù)目是無限的，而實際計算中，為什么原子數(shù)目多了就算不動？

長話短說，不提公式，只說原理：
1、類比于原子軌道。晶態(tài)材料中有‘晶體軌道’。也就是單電子函數(shù)。
2、原子中的電子，按泡利不相容原理，從低能級的軌道開始填充。 VS 晶體中的電子，也按泡利不相容原理從低能量的‘晶體軌道’開始填起。
3、怎樣得到體系的能量呢？（我這里說的不是密度泛函版，這個太難描述了。只說一個最簡單的原型）
a、求出體系的所有‘晶體軌道’，當(dāng)然個軌道所對能量也順便得出來了。
b、電子一個一個，從低能軌道排到高能軌道。各電子的能量，就是它所在軌道的能量。
c、各電子能量求和就OK了。

4、可見，關(guān)鍵在于求解各‘晶體軌道’�；贐loch地理，計算機是這樣處理的。
a、用超胞中各原子的原子軌道的線性疊加來表示‘晶體軌道’。
b、利用變分原理。把該問題轉(zhuǎn)化為求解一個齊次線性方程組。也就是對角化一個矩陣，就可求出各系數(shù)。
c、可見。如果超胞中的原子越多，則參與計算的原子軌道及越多。需要對角化的矩陣就越大。所有不能處理太大的超胞。
5、另外，每個原子的原子軌道是無限多條。晶體軌道也是無限多條。如果不做處理，那個矩陣是無窮維的�！疅o窮×無窮’的矩陣
a、引入截斷動能。能量高于這個值的原子軌道則不參與計算。
這樣一來‘無窮×無窮’簡化為 ‘無窮×有窮’
b、還有一個方向是無限維。該K空間踩點來處理了。你不是無窮嗎？我只踩有限個點來計算。
于是： ‘無窮×有窮’再次簡化為 ‘有窮×有窮’

總結(jié)：胞越大，截斷動能越大，K點越大。要對角化的矩陣就越大。

水平有限，必有疏漏。恭請指正。

回復(fù)此樓

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

VASP

» 猜你喜歡

11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) 已經(jīng)有6人回復(fù)
求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分已經(jīng)有8人回復(fù)
環(huán)境工程 085701，267求調(diào)劑已經(jīng)有7人回復(fù)
材料工程專碩求調(diào)劑已經(jīng)有8人回復(fù)
297求調(diào)劑已經(jīng)有3人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有7人回復(fù)
281求調(diào)劑已經(jīng)有7人回復(fù)
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行已經(jīng)有8人回復(fù)
求收留已經(jīng)有6人回復(fù)
085600 295分求調(diào)劑已經(jīng)有9人回復(fù)

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

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

我的差分電荷密度正常嗎？能不能給個提示。謝謝啦！已經(jīng)有19人回復(fù)
TEM 原子序數(shù)大的原子對比度更大嗎已經(jīng)有8人回復(fù)
求大牛指教：怎么根據(jù)DTA-TG曲線制定燒結(jié)制度~ 已經(jīng)有4人回復(fù)
兩個原子的局域態(tài)密度在費米處重疊已經(jīng)有4人回復(fù)
新手，請高手幫忙，計算小分子中某個碳原子上的電子云密度已經(jīng)有15人回復(fù)
關(guān)于K點選擇與原子間的力已經(jīng)有9人回復(fù)
【求助完畢】這篇JACS的communicaton是如何從電荷密度圖看出共價鍵來的？已經(jīng)有12人回復(fù)
【求助】請教幾個有關(guān)表面活性劑的問題(標(biāo)題要長，金幣要多，虛心求指教，額外金幣奉?/a> 已經(jīng)有13人回復(fù)

【請教】向大俠請教，已知空間群，并且已經(jīng)有粉末衍射數(shù)據(jù)，怎么確定原子位置？已經(jīng)有4人回復(fù)
【求助】誰做ZSM-5分子篩，請高手指教已經(jīng)有24人回復(fù)
【求助】電子密度圖中背景電子密度的概念？已經(jīng)有12人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2012-07-06 11:20:48

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

dyc_2008

捐助貴賓 (著名寫手)

應(yīng)助: 14 (小學(xué)生)
金幣: 124.2
帖子: 1185
在線: 244.3小時
蟲號: 1044421

引用回帖:

2樓: Originally posted by uuv2010 at 2012-07-06 12:05:27
只是說其中一點，也就是3c所說的“c、各電子能量求和就OK了�！�
恐怕是不可以的，即簡單地將各個電子的能量求和并不能描述體系的真實的能量狀態(tài)，至少交換能沒有考慮進去。

我已經(jīng)說了這是最基本的模型
不是DFT
沒有考慮交互作用
僅僅就是單電子方程

謝謝參與

贊一下

回復(fù)此樓

3樓2012-07-06 12:33:33

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

查看全部 5 個回答

uuv2010

榮譽版主 (職業(yè)作家)

★ ★
dyc_2008(金幣+1): 謝謝參與
dyc_2008: 金幣+1 2012-09-11 09:10:26

只是說其中一點，也就是3c所說的“c、各電子能量求和就OK了。”
恐怕是不可以的，即簡單地將各個電子的能量求和并不能描述體系的真實的能量狀態(tài)，至少交換能沒有考慮進去。

贊一下

回復(fù)此樓

2樓2012-07-06 12:05:27

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

uuv2010

榮譽版主 (職業(yè)作家)

★ ★ ★
dyc_2008: 金幣+3 2012-09-11 09:10:35

無論DFT還是HF，都是需要考慮交換能的，呵呵
如果是想計算氫原子的能級的話(這個可能是LZ所說的最簡單模型)，沒有必要考慮原子個數(shù)和多k點計算

另外，關(guān)于“為什么原子數(shù)目多了就算不動”的問題，是因為DFT的計算量是隨著原子數(shù)目的立方的關(guān)系增加的

贊一下

回復(fù)此樓

4樓2012-07-06 15:50:41

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

dyc_2008

捐助貴賓 (著名寫手)

應(yīng)助: 14 (小學(xué)生)
金幣: 124.2
帖子: 1185
在線: 244.3小時
蟲號: 1044421

引用回帖:

4樓: Originally posted by uuv2010 at 2012-07-06 15:50:41
無論DFT還是HF，都是需要考慮交換能的，呵呵
如果是想計算氫原子的能級的話(這個可能是LZ所說的最簡單模型)，沒有必要考慮原子個數(shù)和多k點計算

另外，關(guān)于“為什么原子數(shù)目多了就算不動”的問題，是因為DFT的計 ...

感謝您老的熱情關(guān)注。（老，不是真的老，而是一種尊稱）
我就是想問問，VASP的骨架是不是這樣的？

交互作用什么的，就是血肉了。不可忽略。

贊一下

回復(fù)此樓

5樓2012-07-06 16:27:30

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

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

[考研] 環(huán)境工程 085701，267求調(diào)劑 +7	minht 2026-03-29	7/350	2026-03-30 20:46 by dophin1985
[考研] 303求調(diào)劑 +5	DLkz1314. 2026-03-30	5/250	2026-03-30 20:06 by 無際的草原
[考研] 085602化工求調(diào)劑（331分） +7	111@127 2026-03-30	7/350	2026-03-30 20:02 by 無際的草原
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +10	舉個栗子oi 2026-03-24	11/550	2026-03-30 19:31 by michael2011
[考研] 求調(diào)劑 +10	家佳佳佳佳佳 2026-03-29	10/500	2026-03-30 18:34 by 544594351
[考研] 291求調(diào)劑 +8	HanBeiNingZC 2026-03-24	8/400	2026-03-30 17:01 by 無際的草原
[碩博家園] 求調(diào)劑有機化學(xué)考研356分 +10	Nadiums 2026-03-25	11/550	2026-03-30 11:56 by yjolah
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 生物技術(shù)與工程 +7	1294608413 2026-03-25	8/400	2026-03-30 11:36 by 唐沐兒
[考研] 求化學(xué)調(diào)劑 +11	wulanna 2026-03-28	11/550	2026-03-30 10:59 by 探123
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8	mikasa的圍巾 2026-03-28	8/400	2026-03-29 12:48 by 無際的草原
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8	AZMK 2026-03-25	12/600	2026-03-28 18:15 by AZMK
[考博] 26申博 +3	加油沖�。� 2026-03-26	3/150	2026-03-27 15:38 by cls512
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 296求調(diào)劑 +4	汪��？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper