 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3401 | 回復(fù): 12 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
zhangguangping木蟲 (著名寫手)
|
[求助]
K點(diǎn)在TranSIESTA計(jì)算中的作用
|
|||
|
在計(jì)算TranSIESTA輸運(yùn)的時(shí)候，如果橫向(x,y方向)使用了周期性邊界條件(這兒指的是真正的電極周期性，而不是程序包默認(rèn)的周期性邊界條件)，那么就需要在對(duì)應(yīng)的二維布里淵區(qū)要取k點(diǎn)。 首先，根據(jù)固體物理我們知道，對(duì)于周期體系，所有的物理量都是k點(diǎn)的函數(shù)。所以要想計(jì)算一個(gè)物理量，需要對(duì)k點(diǎn)求平均。原則上k點(diǎn)需要取無窮多個(gè)才能計(jì)算準(zhǔn)確。但是，考慮到數(shù)值計(jì)算時(shí)間和精度的容忍度，我們可以選擇有限的k點(diǎn)來近似計(jì)算這些物理量。因此，k點(diǎn)越多，計(jì)算結(jié)果越接近準(zhǔn)確值。 在TranSIESTA中，體系的哈密頓對(duì)k點(diǎn)收斂很快，而透射譜的計(jì)算對(duì)k點(diǎn)收斂卻很慢。如圖FIG1是分別采用4*4 k點(diǎn)和20*20 k點(diǎn)的結(jié)果。從計(jì)算結(jié)果可以看出，在費(fèi)米能級(jí)上方附近，4*4的k點(diǎn)有很多尖銳的透射峰。當(dāng)取20*20之后，這個(gè)情況得到緩解。但是如果對(duì)20*20的再繼續(xù)放大，如FIG2,就會(huì)發(fā)現(xiàn)，這種鋸齒狀的波動(dòng)還是存在。 問題1：這種鋸齒狀的情況以前從來沒有碰到過，并且對(duì)于以前的體系，采用較少的4*4的k點(diǎn)就很平滑了。想請(qǐng)教大家，這種鋸齒狀的透射譜，能說明什么呢？為什么隨著k點(diǎn)的收斂這么慢？這個(gè)能給出一些什么信息呢？ 問題2：對(duì)于計(jì)算透射譜的時(shí)候，這個(gè)k點(diǎn)什么含義呢？為什么要對(duì)k點(diǎn)求平均呢？ 對(duì)于某一個(gè)能量E，要對(duì)所有取的k點(diǎn)求平均。這個(gè)k點(diǎn)對(duì)應(yīng)電子的能量嗎？如果對(duì)應(yīng)的話，是否對(duì)應(yīng)一個(gè)k點(diǎn)（k_x,k_y）有E=(k_x*k_x+k_y*k_y)/2m嗎？那k_z怎么處理呢？ 真心請(qǐng)大牛出來指教. |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 Siesta（Transiesta、ATK） | 計(jì)算-siesta |
» 猜你喜歡
342求調(diào)劑
已經(jīng)有3人回復(fù)
-
070300化學(xué)專業(yè)279調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有20人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿a區(qū)211，085601-307分求調(diào)劑
已經(jīng)有23人回復(fù)
-
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑
已經(jīng)有8人回復(fù)
-
材料調(diào)劑
已經(jīng)有10人回復(fù)
-
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有13人回復(fù)
-
材料科學(xué)與工程339求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué)308分調(diào)劑
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 討論：在ATK和TranSIESTA輸運(yùn)中，大家對(duì)k-sampling的選取問題
已經(jīng)有16人回復(fù)
- 【求助】siesta中關(guān)于幾何結(jié)構(gòu)優(yōu)化的問題【完結(jié)】
已經(jīng)有17人回復(fù)
|
我用的是我導(dǎo)師自己開發(fā)的軟件（PRB 70.085410），類似于ATK的那種，我確實(shí)問過我導(dǎo)師，為什么計(jì)算輸運(yùn)前是進(jìn)行單點(diǎn)的DFT計(jì)算（不是優(yōu)化），而不是取很多的K點(diǎn)，老師的回答只是說跟K點(diǎn)沒有關(guān)系了，可是我明明發(fā)現(xiàn)Qup-Qdown 不一樣。 如果你有興趣讀我導(dǎo)師的博士后，當(dāng)然你可以直接問我導(dǎo)師的（同濟(jì)大學(xué)柯三黃）,我相信他在電子輸運(yùn)這塊應(yīng)該是國內(nèi)頂尖的， |
4樓2013-06-09 17:00:21
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+2, 謝謝交流 2013-06-09 12:42:25
zhangguangping: 金幣+10, ★有幫助 2013-06-09 17:26:44
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+2, 謝謝交流 2013-06-09 12:42:25
zhangguangping: 金幣+10, ★有幫助 2013-06-09 17:26:44
|
對(duì)于輸運(yùn)圖，計(jì)算的時(shí)候應(yīng)該是單點(diǎn)計(jì)算，也就是跟K點(diǎn)沒有關(guān)系，所以一般我計(jì)算輸運(yùn)的時(shí)候都是取單個(gè)的K點(diǎn)，，具體步驟是先用很多K點(diǎn)優(yōu)化，然后進(jìn)行單點(diǎn)的DFT計(jì)算，最后進(jìn)行輸運(yùn)計(jì)算，關(guān)于為什么取單個(gè)K點(diǎn)，我問過我導(dǎo)師，他說跟K沒關(guān)系，不過好像Datta的第一本關(guān)于輸運(yùn)的里面涉及了 ，最后你這個(gè)輸運(yùn)圖好奇怪,不知道你算得是什么的輸運(yùn)？ |
2樓2013-06-09 12:35:14
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
3樓2013-06-09 16:15:30
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
建議讀讀有關(guān)波函數(shù)匹配的方法來計(jì)算電子輸運(yùn)的文獻(xiàn)，就可以比較清楚地理解transimission spectra等計(jì)算中的波矢k的含義。 比如：PHYSICAL REVIEW B 70, 045417 (2004) 有關(guān)k點(diǎn)取樣對(duì)transimission spectra計(jì)算的影響： M. Stilling , K. Stokbro & K. Flensberg (2007): Crystalline Magnetotunnel Junctions: Fe-mgo-fe, Fe- feomgo-fe And Fe-aumgoau-fe, Molecular Simulation, 33:7, 557-561 http://www.tandfonline.com/doi/a ... 930508#.UbZstpzRL0c |
6樓2013-06-11 08:18:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +6 | 小小機(jī)靈蟲 2026-03-29 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +20 | Mr. Z 2026-03-25 | 21/1050 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +10 | 傾____蕭 2026-03-27 | 11/550 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +15 | Micky11223 2026-03-25 | 16/800 |
|
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|