| 查看: 3009 | 回復(fù): 17 | ||||
xirainbow木蟲 (正式寫手)
|
[交流]
【原創(chuàng)】The meaning of nband 已有1人參與
|
|
MS版上發(fā)了個(gè)貼,討論nband的意義 PWscf也牽涉到這個(gè)問題。 以下為我自己的理解,其中最后一段有待商榷 原帖網(wǎng)址:http://blog.sina.com.cn/s/blog_5f15ead20100e1l5.html 由于上傳圖片比較麻煩 所以原帖中的很多公式我沒有上傳 有興趣的可以去博客上看;) This is my second try. I spend almost two days to get some progress in understanding of Davidson diagonalization. I am not interested in Davidson method itself. But without knowing it, you can not understand the real meaning of parameter "nband". 在計(jì)算中指定了截止動(dòng)能,對(duì)于一個(gè)固定的k點(diǎn),也就意味著指定了有多少個(gè)G(=2pi/a的整數(shù)倍)參與傅立葉變換(假設(shè)為M個(gè))。而G的個(gè)數(shù)就是相應(yīng)哈密頓量的維度,M*M維。通常M是非常大的,而我們只需要知道較低的幾個(gè)被占據(jù)的能級(jí),即可求出總能等一些列剛興趣的物理量。因此直接對(duì)M*M維的哈密頓量進(jìn)行對(duì)角化是不合適的。所以需要采用迭代的方法來求解幾個(gè)低能的本征值和本征能量即可。 其實(shí)迭代方法在矩陣對(duì)角化中是非常常用的數(shù)值計(jì)算方法,他可以通過迭代的方式只求出某一個(gè)或者幾個(gè)本征值。Davidson是其中一個(gè)非常出名的,因?yàn)樗诹孔恿W(xué)中取得了很大的成功,這個(gè)成功來源于量子力學(xué)的哈密頓矩陣通常都是稀疏的,就是非對(duì)角元很多都是0,或者接近0。 當(dāng)nband被指定后,本征方程的就被給定。。本征方程的維度N應(yīng)該是等于nband的。此時(shí)N遠(yuǎn)小于有ecutwfc決定的G的個(gè)數(shù),即H的維度M*M。而波函數(shù)psi則是由M個(gè)平面波的疊加。 然后求得對(duì)初始猜測波函數(shù)的修正。這樣的修正波函數(shù)作為獨(dú)立的基矢,加入到初始猜測的N個(gè)波函數(shù)上,構(gòu)成N+1維的子空間。如果對(duì)N個(gè)初始猜測波函數(shù)均采用這種辦法,則得到2*N維的子空間。然后,在2*N維的子空間中求得H矩陣,并對(duì)其對(duì)角化,比較前N個(gè)本征值是否有變化,判定是否收斂。 (本征方程的維度是否在不斷擴(kuò)大,有待商榷) 納米盤下載: http://d.namipan.com/sd/1558396 [ Last edited by wuli8 on 2009-12-22 at 19:13 ] |
第一性原理計(jì)算 |
木蟲 (著名寫手)
榮譽(yù)版主 (著名寫手)

木蟲 (正式寫手)
木蟲 (正式寫手)
鐵蟲 (初入文壇)
木蟲 (正式寫手)
金蟲 (著名寫手)
木蟲 (正式寫手)
金蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 308 085701 四六級(jí)已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|