| 查看: 3529 | 回復(fù): 11 | ||||
yashumimi銀蟲 (小有名氣)
|
[交流]
【求助】二維材料k點(diǎn)的設(shè)置 已有1人參與
|
|||
|
今天在算k點(diǎn)不同對transiesta計(jì)算結(jié)果的影響 k點(diǎn)設(shè)置為 1 0 0 0 1 0 0 0 i 其中i為變量 我將其設(shè)為2-20 因?yàn)殚_始用1算了 出錯(cuò) 是個(gè)二維的結(jié)構(gòu) x方向有一定的寬度 沒有周期性 加了真空層 y方向不變 z方向具有周期性 請問k點(diǎn)應(yīng)該怎么設(shè)置才合理呢? 另外我更改k點(diǎn)后用tbtrans計(jì)算了I-V 同樣的操作 兩次 電流數(shù)量級有小幅度變化 趨勢變化規(guī)律也不一樣 這說明計(jì)算的結(jié)果不穩(wěn)定么? |
transport |
木蟲 (著名寫手)
|
對于scattering region一般是通過優(yōu)化得到中間結(jié)構(gòu)的。但是一旦做輸運(yùn)計(jì)算的時(shí)候,就不在考慮偏壓下的幾何結(jié)構(gòu)變化了。因?yàn)檫@個(gè)問題Transiesta是處理不了的。再說即使處理得了,這個(gè)時(shí)間也是負(fù)擔(dān)不起的。所以這個(gè)地方就有一個(gè)近似了。至于transiesta計(jì)算中輸出的力,已經(jīng)不再具有參考價(jià)值。不過聽說atk-2010.8可以在偏壓下進(jìn)行幾何結(jié)構(gòu)的優(yōu)化。具體的執(zhí)行方法我不是很清楚,但是明確的告訴你,transiesta不能勝任這個(gè)工作。 所以在做輸運(yùn)的時(shí)候,要將MD.NumCGsteps 設(shè)置為0,也就是說所有的關(guān)于MD的選項(xiàng)不再使用了。 關(guān)于這個(gè)你可以到siesta自帶的\Tests\transiesta目錄下,查看一下所有的散射區(qū)輸運(yùn)計(jì)算的輸入文件都是這樣設(shè)置的。 關(guān)于這個(gè)我就知道這么寫,如果有疑問,可以到siesta的mail_list上提問。 |

木蟲 (著名寫手)

木蟲 (著名寫手)
|
你把你的輸入文件貼出來,高手可以用你的輸入文件具體查看一下,有時(shí)候你只是這么一說很難排查。我也想知道這個(gè)問題的答案。 我一直感覺對于碳納米管這種一維體系,k點(diǎn)設(shè)置如下就行: %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 6 0.0 %endblock kgrid_Monkhorst_Pack 其中第一二行都設(shè)置為1,因?yàn)檫@兩個(gè)方向是沒有周期性的。(碳管沿著z方向)。第三行是為了計(jì)算電極部分的碳管用的,對于碳管這種體系,可能這個(gè)地方需要的k點(diǎn)數(shù)目比較多。需要做一個(gè)收斂測試。但是對于輸運(yùn)計(jì)算的時(shí)候,第三行就沒設(shè)么用了。只看第一二行。 |

銀蟲 (小有名氣)
銀蟲 (小有名氣)
銀蟲 (小有名氣)
|
SystemName SystemLabel %block ChemicalSpeciesLabel 1 1 H 2 6 C %endblock ChemicalSpeciesLabel PAO.EnergyShift 50 meV PAO.BasisType split PAO.BasisSize DZP %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 20 0.0 %endblock kgrid_Monkhorst_Pack xc.functional GGA # Exchange-correlation functional xc.authors PBE # Exchange-correlation version SpinPolarized F # Logical parameters are: yes or no FixSpin F MeshCutoff 300. Ry # Equivalent plane wave cutoff for the grid # SCF options DM.Tolerance 1.d-4 # Tolerance in maximum difference between input and output DM MaxSCFIterations 1000 # Maximum number of SCF iter DM.UseSaveDM T # to use continuation files DM.MixingWeight 0.02 # New DM amount for next SCF cycle DM.NumberPulay 6 DM.MixSCF1 F DM.PulayOnFile F # Store in memory ('F') or in files ('T') # NeglNonOverlapInt T # 'F'=do not neglect SolutionMethod Diagon # OrderN or Diagon ElectronicTemperature 300 K # Temp. for Fermi smearing # MD options MD.TypeOfRun CG # Type of dynamics: MD.VariableCell F MD.NumCGsteps 100 # Number of CG steps for coordinate optimization %include POSITION.fdf 輸入文件差不多就是這個(gè)樣子的 有個(gè)問題用一個(gè)小的架構(gòu)優(yōu)化后supercell 力還是很大 只能再優(yōu)化么 有沒有好的辦法提高優(yōu)化速度 另外用siesta優(yōu)化和transiesta優(yōu)化的區(qū)別在哪里? |
木蟲 (著名寫手)

銀蟲 (小有名氣)
木蟲 (著名寫手)

銀蟲 (小有名氣)
|
Dear Transiesta user: I have one question about the calculation of Transiesta. I am using Transiesta to calculate the transmission in nanowires. I am not sure about the optimization of the scattering region using siesta or transiesta. Since if I use siesta to optimize the scattering region, the results are not correct and if I use the Transiesta to optimize the scattering region, it can not converge . Could anybody help me? Thanks in advance. Lihong 我覺得我和他遇到的問題是一樣的 就是加上偏壓后力變得很大 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|