| 5 | 1/1 | 返回列表 |
| 查看: 2326 | 回復(fù): 13 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
redskywei木蟲 (著名寫手)
|
[交流]
【求助】siesta中transiesta模塊的收斂問題 已有3人參與
|
|||
|
最近在做transieta的計算,發(fā)現(xiàn)收斂困難,看maillist里主要考慮增加屏蔽層和設(shè)置TS.UpdateDMCROnly,我在設(shè)置TS.UpdateDMCROnly時,一直不明白這個東西到底是作用在整個散射區(qū)還是僅僅與電極結(jié)構(gòu)相重復(fù)的散射區(qū)中左右單元的過度層?這個coupling terms該怎么理解? 還有期中有個人提到TS.UpdateDMCROnly + extra elec layers,又該怎么理解,是說對外電極設(shè)置這個參數(shù)么? 我用3層Al電極做左右屏蔽層,TS.UpdateDMCROnly設(shè)置為 T了,還是收斂非差差。近300步才收斂。也不知道哪里問題哦。我的輸入?yún)?shù)如下,請大家指教: NumberOfSpecies 3 NumberOfAtoms 287 %block ChemicalSpeciesLabel 1 6 C 2 7 N 3 13 Al %endblock ChemicalSpeciesLabel PAO.EnergyShift 200 meV PAO.BasisType split PAO.BasisType split PAO.BasisSize DZP LatticeConstant 1.0 Ang %block LatticeVectors 16.19832 0.00000000 0.00000000 0.00000000 16.19832 0.00000000 0.00000000 0.00000000 26.80535368 %endblock LatticeVectors #kgrid_cutoff 20. Ang %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 1 0.5 %endblock kgrid_Monkhorst_Pack # SPIN options xc.functional GGA # Exchange-correlation functional xc.authors PBE # Exchange-correlation version SpinPolarized F # Logical parameters are: yes or no FixSpin F TotalSpin 0.0 NonCollinearSpin F # 'T', 'F' MeshCutoff 150. 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 300 # Maximum number of SCF iter DM.UseSaveDM F # to use continuation files DM.MixingWeight 0.03 # New DM amount for next SCF cycle DM.NumberPulay 5 DM.MixSCF1 F DM.PulayOnFile F # Store in memory ('F') or in files ('T') # NeglNonOverlapInt T # 'F'=do not neglect SolutionMethod transiesta # OrderN or Diagon or transiesta ElectronicTemperature 300 K # Temp. for Fermi smearing # MD options MD.TypeOfRun CG # Type of dynamics: MD.VariableCell F MD.NumCGsteps 000 # Number of CG steps for coordinate optimization ZM.UnitsLength Ang #the units of length used during Z-matrix input ZM.UnitsAngle deg #the units of angles used during Z-matrix input ZM.ForceTolLength 0.02 eV/Ang #controls the convergence with respect to forces on Z-matrix lengths ZM.ForceTolAngle 0.00356549 Ry/rad #controls the convergence with respect to forces on Z-matrix angles ZM.MaxDisplLength 0.1 Ang # controls the maximum change in a Zmatrix length during an optimisation step ZM.MaxDisplAngle 0.003 rad # controls the maximum change in a Z-matrix angle during an optimisation step ZM.CalcAllForces T # Default value [ Last edited by redskywei on 2010-5-29 at 18:21 ] |
SIESTA |

木蟲 (著名寫手)

木蟲 (著名寫手)

金蟲 (著名寫手)
實習(xí)版主
|
這個收斂參數(shù)沒有問題,200meV。 MeshCutoff=150.Ry 你在看一下其他參數(shù)。 還有可能就是體系的問題,同樣結(jié)構(gòu)的BN比Graphene就會要難收斂好多,體系完全一樣,把原子替換成BN之后 計算量是Graphene的好幾倍 |

木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境285分,過六級,求調(diào)劑 +5 | xhr12 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +13 | jjjjjjjjjjl 2026-03-27 | 13/650 |
|
|
[考研] 314求調(diào)劑 +11 | 1xiaojun23 2026-03-31 | 12/600 |
|
|
[考研] 求調(diào)劑推薦 +3 | 南山南@ 2026-04-01 | 3/150 |
|
|
[考研] 08開頭看過來!。 +4 | wwwwffffff 2026-03-31 | 6/300 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +13 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 385分 生物學(xué)(071000)求調(diào)劑 +6 | qf626 2026-04-01 | 6/300 |
|
|
[考研] 324分 085600材料與化工 +20 | 呆鵝oor 2026-03-27 | 20/1000 |
|
|
[考研] 367求調(diào)劑 +8 | 芋泥啵! 2026-03-28 | 8/400 |
|
|
[考研] 材料調(diào)劑 +12 | 一樣YWY 2026-04-01 | 12/600 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +7 | 烤麥芽 2026-03-27 | 10/500 |
|
|
[考研] 材料專碩調(diào)劑 +17 | 椰椰。 2026-03-29 | 17/850 |
|
|
[考研] 291求調(diào)劑 +3 | 迷蒙木木 2026-04-01 | 4/200 |
|
|
[考研] 299求調(diào)劑 +8 | 嗯嗯嗯嗯2 2026-03-27 | 8/400 |
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
[考研] 一志愿大連理工大學(xué),機械工程學(xué)碩,341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|