版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

freewain

銅蟲 (初入文壇)

應助: 3 (幼兒園)
金幣: 293.9
帖子: 39
在線: 35.4小時
蟲號: 987873

[交流] 【求助】Transiesta 計算出錯

大家好，我在計算 ZnO 納米線的輸運性質。這是我的輸入文件：

#General system descriptors
SystemLabel ZnO
SystemName Au-ZnO-Au

NumberOfSpecies 3
NumberOfAtoms 50

%block ChemicalSpeciesLabel
1 30 Zn
2 8 O
3 79 Au
%endblock ChemicalSpeciesLabel

#Basis set and KB projectors
PAO.BasisSize DZP
PAO.EnergyShift 0.02 Ry

#Structural information
LatticeConstant 1.00 Ang
MD.UseStructFile true
WriteCoorXmol true
WriteMDXmol true

#k-point sampling
%block kgrid_Monkhorst_Pack
1  0  0 0.0
0  1  0 0.0
0  0 30 0.0
%endblock kgrid_Monkhorst_Pack

#Exchange-correlation functionals
XC.functional GGA
XC.authors PBE

#SCF loop
MaxSCFIterations 150
DM.MixingWeight 0.05
DM.NumberPulay 5

#Other
MeshCutoff 400 Ry
SaveHS true
SolutionMethod transiesta

#Band structure analysis
BandLinesScale ReciprocalLatticeVectors
%block BandLines
1 0.0  0.0  0.0  \Gamma
20 0.0  0.0  0.5  X
%endblock BandLines

#DOS
%block ProjectedDensityOfStates
-10.0 5.0 0.10 151 eV
%endblock ProjectedDensityOfStates

%block LocalDensityOfStates
-6.00 2.00 eV
%endblock LocalDensityOfStates

COOP.Write true

#charge densities and potentials on the grid
SaveRho true
SaveDeltaRho true
SaveElectrostaticPotential true
SaveTotalPotential true

WriteDenchar true

#Relaxation
MD.TypeOfRun CG
MD.NumCGSteps 0

UseSaveData true

WriteMDhistory true

%block GeometryConstraints
  position from  1 to  19
  position from -1 to -19
%endblock GeometryConstraints

TS.HSFileLeft /home/niu/new/elec/Au.TSHS
TS.NumUsedAtomsLeft 19
TS.HSFileRight /home/niu/new/elec/Au.TSHS
TS.NumUsedAtomsRight 19

TS.ComplexContour.Emin -50.0 eV
TS.ComplexContour.NPoles 03
TS.ComplexContour.NCircle 30
TS.ComplexContour.NLine 10
TS.BiasContour.NumPoints 20
WriteEigenValues true

TS.SaveHS true
TS.MixH true

TS.TBT.NPoints 400
TS.TBT.Emin -1.00 eV
TS.TBT.Emax    1.00 eV
TS.TBT.OutputRegionData False
TS.TBT.NEigen 5
TS.TBT.PDOSFrom 3
TS.TBT.PDOSTo 9
TS.Voltage -1.00 eV

但是計算中間出現(xiàn)錯誤，輸入文件的后半部分是這樣的：

Begin LEFT
unit cell:
37.7945  0.0000  0.0000
  0.0000 37.7945  0.0000
  0.0000  0.0000 13.0819
Size:       81225
Efermi: -3.6748994224532205E-002
GF: Spin number ->          1
Got LEFT Electrode GF
Done
Begin RIGHT
unit cell:
37.7945  0.0000  0.0000
  0.0000 37.7945  0.0000
  0.0000  0.0000 13.0819
Size:       81225
Efermi: 3.6748994224532205E-002
GF: Spin number ->          1
Got RIGHT  Electrode GF
Done
L-ELEC: lastoL:          0          15          30          45
         60          75          90       105       120       135
      150       165       180       195       210       225
      240       255       270       285
R-ELEC: lastoR:          0          15          30          45
         60          75          90       105       120       135
      150       165       180       195       210       225
      240       255       270       285
TRANSIESTA: Initializing lists
From SIESTA: Efermi=  -3.171
qc0 : 0.000
qcn0 : 419.530
nou,NGL,NGR:       738       285       285
TranSiesta: Qsol,Qtot:  526.113  526.000 1.000

* Maximum dynamic memory allocated = 423 MB

siesta: Program's energy decomposition (eV):
siesta: Ebs    =    -3712.819471
siesta: Eions =    56290.310725
siesta: Ena    =    2049.915289
siesta: Ekin =    35787.957698
siesta: Enl    = -18060.566835
siesta: DEna =    301.485341
siesta: DUscf =       36.453196
siesta: DUext =       0.000000
siesta: Exc    = -131533.815024
siesta: eta*DQ  =       0.000000
siesta: Emadel  =       0.000000
siesta: Emeta =       0.000000
siesta: Emolmec =       0.000000
siesta: Ekinion =       0.000000
siesta: Eharris = -167708.882424
siesta: Etot = -167708.881061
siesta: FreeEng = -167708.881061

transiesta: TSiscf Eharris(eV)    E_KS(eV) FreeEng(eV) dDmax  Ef(eV)
transiesta: 1  -167708.8824  -167708.8811  -167708.8811  0.0001 -3.1707
timer: Routine,Calls,Time,% = IterSCF       1    0.000 0.00
elaps: Routine,Calls,Wall,% = IterSCF       1    0.000 0.00
timer: Routine,Calls,Time,% = TS          1 3611.599  89.88
elaps: Routine,Calls,Wall,% = TS          1    902.900  89.88
timer: Routine,Calls,Time,% = IterSCF       1 3992.787  99.36
elaps: Routine,Calls,Wall,% = IterSCF       1    998.197  99.36

siesta: E_KS(eV) =       -167708.8811

siesta: E_KS - E_eggbox = -167708.8811
rank 0 in job 5  amd_34368 caused collective abort of all ranks
  exit status of rank 0: killed by signal 9

請問這是什么原因造成的呢？謝謝解答：）

回復此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有20人回復
285求調(diào)劑已經(jīng)有8人回復
材料與化工（0856）304求B區(qū)調(diào)劑已經(jīng)有8人回復
370求調(diào)劑已經(jīng)有3人回復
070300化學279求調(diào)劑已經(jīng)有5人回復
福建理工大學材料學院先進合金團隊招收考研調(diào)劑學生已經(jīng)有4人回復
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑已經(jīng)有5人回復
083000環(huán)境科學與工程調(diào)劑，總分281 已經(jīng)有4人回復
求調(diào)劑已經(jīng)有3人回復
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有10人回復

» 本主題相關商家推薦: (我也要在這里推廣)

» 本主題相關價值貼推薦，對您同樣有幫助:

VASP并行計算出錯已經(jīng)有9人回復
MS discover 計算出錯，請高手指點一下已經(jīng)有4人回復
求助：gaussian計算出錯。已經(jīng)有3人回復
transiesta計算運行出錯已經(jīng)有4人回復
討論：在ATK和TranSIESTA輸運中，大家對k-sampling的選取問題已經(jīng)有16人回復
關于Dictra的計算問題（主要是關于計算過程出錯的問題）已經(jīng)有8人回復
【求助】BSSE計算出錯已經(jīng)有16人回復
【求助】siesta3.0 安裝，錯誤一籮筐已經(jīng)有11人回復
【求助】高斯計算出錯已經(jīng)有11人回復
【求助】siesta中transiesta模塊的收斂問題已經(jīng)有13人回復
【求助】BSSE計算出錯已經(jīng)有19人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-11-18 15:11:28

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiongxiong5712

新蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 543
帖子: 257
在線: 267.7小時
蟲號: 706467

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

32744樓: Originally posted by zhangguangping at 2010-11-18 16:01:29
你把你的輸入文件里面加上一行
TS.CalcGF  F
并且把下面的東西注釋掉
#Band structure analysis
BandLinesScale             ReciprocalLatticeVectors
%block BandLines
1 0.0  0.0  0.0  \Gamma
20 ...

你好，我注意到他的輸入文件里k點的設置是這樣的：
#k-point sampling
%block kgrid_Monkhorst_Pack
1  0  0 0.0
0  1  0 0.0
0  0 30 0.0
%endblock kgrid_Monkhorst_Pack

用transiesta計算，是非周期結構，這樣設置k點正確嗎？不是z方向為非周期嗎，k點應該為1*1*1吧？

贊一下

回復此樓

7樓2012-11-19 09:23:44

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個回答

zhangguangping

木蟲 (著名寫手)

★ ★ ★ ★
freewain(金幣+1):謝謝參與
freewain(金幣+9):謝謝廣平師兄！問題解決了。我是測好電極距離后直接拿過來用了，也沒刪掉態(tài)密度的內(nèi)容。上面的運算是出錯之后又計算的，TSDE是第一遍就算出來的。不過之后的計算，每次會使用上一次的TSDE和DM文件。 2010-11-18 17:58:08
sunyang1988(金幣+3):謝謝指點 2010-11-18 22:02:15

引用回帖:

Originally posted by freewain at 2010-11-18 08:11:28:
大家好，我在計算 ZnO 納米線的輸運性質。這是我的輸入文件：

#General system descriptors
SystemLabel ZnO
SystemName Au-ZnO-Au

NumberOfSpecies 3
NumberOfAtoms 50

%block ChemicalSpeciesLab ...

你把你的輸入文件里面加上一行
TS.CalcGF  F
并且把下面的東西注釋掉
#Band structure analysis
BandLinesScale             ReciprocalLatticeVectors
%block BandLines
1 0.0  0.0  0.0  \Gamma
20 0.0  0.0  0.5  X
%endblock BandLines

#DOS
%block ProjectedDensityOfStates
-10.0 5.0 0.10 151 eV
%endblock ProjectedDensityOfStates

%block LocalDensityOfStates
-6.00 2.00 eV
%endblock LocalDensityOfStates
（這三項，尤其是LDOS，能帶，好像不能在算輸運的時候算。因為他不是一個周期體系了，怎么算能帶。LDOS只有方法是Diagon的時候能算的。）
還在你的原來那個目錄計算�？匆幌逻@次行不行。主要是再算一次。看到你好像用的是已經(jīng)收斂好的*TSDE算的。如果是這樣，這個TSDE是怎么得到的呢？計算到這兒，已經(jīng)到了結尾部分了。

[ Last edited by zhangguangping on 2010-11-18 at 09:11 ]

贊一下(2人)

回復此樓

2樓2010-11-18 16:01:29

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangguangping

木蟲 (著名寫手)

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
sunyang1988(金幣+2):謝謝指點 2010-11-18 22:02:34

引用回帖:

找到出錯原因了？原因是什么��？可以分享一下。是我猜測的原因嗎？以后想用上一個計算的密度，直接用*TSDE文件就行。如果沒有*TSDE，還想加速的話，就是用*DM。程序會先檢查*TSDE文件，如果沒有才去找*DM文件的。如果兩個都沒有，那么就從頭計算，先產(chǎn)生*DM,再去產(chǎn)生*TSDE。

贊一下(2人)

回復此樓

3樓2010-11-18 21:48:34

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

freewain

銅蟲 (初入文壇)

應助: 3 (幼兒園)
金幣: 293.9
帖子: 39
在線: 35.4小時
蟲號: 987873

引用回帖:

Originally posted by zhangguangping at 2010-11-18 21:48:34:

找到出錯原因了？原因是什么��？可以分享一下。是我猜測的原因嗎？以后想用上一個計算的密度，直接用*TSDE文件就行。如果沒有*TSDE，還想加速的話，就是用*DM。程序會先檢查*TSDE文件，如果沒有才去找*DM文件的 ...

就是廣平師兄你說的原因，把那些注釋掉就好了。

另外一個問題，算完 transiesta 之后，再算 tbtrans，結果現(xiàn)在一個多小時了，輸出文件一直停在 Starting K-Point loop for spin: 1 那里。我查了一下官方的郵件列表，一共有兩個相同問題的郵件，你回答了其中一個，說這是正常情況，只要慢慢等就可以了�？墒菚r間也太久了。。。

贊一下

回復此樓

4樓2010-11-18 23:30:31

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個回答

普通表情龍兔虎貓高級回復 (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 291求調(diào)劑 +8	Y-cap 2026-03-29	11/550	2026-03-30 21:03 by Y-cap
[考研] 各位老師好，我的一志愿為北京科技大學085601材料專碩 +10	Koxui 2026-03-28	10/500	2026-03-30 19:33 by 源_2020
[考研] 327求調(diào)劑 +5	小卡不卡. 2026-03-29	5/250	2026-03-30 19:30 by Wang200018
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +3	lmnlzy 2026-03-30	3/150	2026-03-30 17:49 by wangjy2002
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +12	芋泥小鈴鐺 2026-03-28	12/600	2026-03-30 17:06 by wangjy2002
[考研] 求調(diào)劑323材料與化工 +10	1124361 2026-03-24	10/500	2026-03-30 16:26 by 690616278
[考研] 求調(diào)劑 +7	青春裁為三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導師 +7	傾____蕭 2026-03-27	8/400	2026-03-30 09:37 by longlotian
[考研] 一志愿北京工業(yè)大學，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +7	腳滑的守法公民 2026-03-27	7/350	2026-03-29 20:21 by dophin1985
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無際的草原
[考研] 0703化學調(diào)劑，求導師收 +9	天天好運來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 調(diào)劑 +3	好好讀書。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 0703化學求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[有機交流] 高溫高壓反應求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 085601 材料工程 313分求調(diào)劑 +5	Ong3 2026-03-27	5/250	2026-03-27 12:24 by goldfish51
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

freewain

[交流] 【求助】Transiesta 計算出錯

» 猜你喜歡

» 本主題相關商家推薦: (我也要在這里推廣)

» 本主題相關價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

xiongxiong5712

zhangguangping

zhangguangping

freewain

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 本主題相關價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼