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

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

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

回復(fù)此樓

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

VASP

» 猜你喜歡

食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
311（085601）求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有4人回復(fù)
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑已經(jīng)有7人回復(fù)
一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑已經(jīng)有12人回復(fù)
289求調(diào)劑已經(jīng)有11人回復(fù)
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 已經(jīng)有6人回復(fù)
321求調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(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的回帖

uuv2010

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

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

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

贊一下

回復(fù)此樓

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

已閱回復(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ìn)去。

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

謝謝參與

贊一下

回復(fù)此樓

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

已閱回復(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ā)表

[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +4	前行必有光 2026-03-28	5/250	2026-03-28 11:27 by xuxiang
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>
[考研] 0856材料化工調(diào)劑總分330 +10	zhubinhao 2026-03-27	10/500	2026-03-28 03:34 by fmesaito
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 283求調(diào)劑（080500） +4	A child 2026-03-27	4/200	2026-03-27 15:34 by XPU李慶
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 340求調(diào)劑 +3	Amber00 2026-03-26	3/150	2026-03-26 18:57 by 不吃魚的貓
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚。 m
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪！？！
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike