版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

freewain

銅蟲 (初入文壇)

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

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

大家好，我在計算 ZnO 納米線的輸運(yùn)性質(zhì)。這是我的輸入文件：

#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

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

回復(fù)此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
296求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有13人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有15人回復(fù)
生物學(xué)調(diào)劑招人�。�！已經(jīng)有4人回復(fù)

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

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

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

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

查看全部散金貼

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

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

zhangguangping

木蟲 (著名寫手)

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

引用回帖:

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

#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
（這三項(xiàng)，尤其是LDOS，能帶，好像不能在算輸運(yùn)的時候算。因?yàn)樗皇且粋€周期體系了，怎么算能帶。LDOS只有方法是Diagon的時候能算的。）
還在你的原來那個目錄計算�？匆幌逻@次行不行。主要是再算一次�？吹侥愫孟裼玫氖且呀�(jīng)收斂好的*TSDE算的。如果是這樣，這個TSDE是怎么得到的呢？計算到這兒，已經(jīng)到了結(jié)尾部分了。

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

贊一下(2人)

回復(fù)此樓

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

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

zhangguangping

木蟲 (著名寫手)

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

引用回帖:

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

贊一下(2人)

回復(fù)此樓

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

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

freewain

銅蟲 (初入文壇)

應(yīng)助: 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，結(jié)果現(xiàn)在一個多小時了，輸出文件一直停在 Starting K-Point loop for spin: 1 那里。我查了一下官方的郵件列表，一共有兩個相同問題的郵件，你回答了其中一個，說這是正常情況，只要慢慢等就可以了�？墒菚r間也太久了。。。

贊一下

回復(fù)此樓

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

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

zhangguangping

木蟲 (著名寫手)

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+2):謝謝指教 2010-11-19 11:47:00

引用回帖:

Originally posted by freewain at 2010-11-18 16:30:31:

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

另外一個問題，算完 transiesta 之后，再算 tbtrans，結(jié)果現(xiàn)在一個多小時了，輸出文件一直停在 Starting K-Point loop for spin: 1 那里。我查了一下官方的 ...

就是這個樣子。如果你的體系是自旋極化的，還要有
Starting K-Point loop for spin: 2，自旋非極化的在這個地方節(jié)省一半的時間。
等等就好了。看你體系的大小和k點(diǎn)的多少，等待的時間不一。

贊一下(2人)

回復(fù)此樓

5樓2010-11-19 08:39:32

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

xiongxiong5712

新蟲 (小有名氣)

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

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

引用回帖:

30913樓: Originally posted by freewain at 2010-11-18 15:11:28
大家好，我在計算 ZnO 納米線的輸運(yùn)性質(zhì)。這是我的輸入文件：

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

NumberOfSpecies 3
NumberOfAtoms 50

%block ChemicalSpeciesLabel ...

麻煩請問一下，你的文件里這兩個參數(shù)是什么意思？能否解釋一下？我在手冊里面沒有看到過
TS.TBT.PDOSFrom 3
TS.TBT.PDOSTo 9

贊一下

回復(fù)此樓

6樓2012-11-17 17:17:39

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

xiongxiong5712

新蟲 (小有名氣)

應(yīng)助: 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點(diǎn)的設(shè)置是這樣的：
#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計算，是非周期結(jié)構(gòu)，這樣設(shè)置k點(diǎn)正確嗎？不是z方向?yàn)榉侵芷趩�，k點(diǎn)應(yīng)該為1*1*1吧？

贊一下

回復(fù)此樓

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

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

zhangguangping

木蟲 (著名寫手)

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

引用回帖:

1981532樓: Originally posted by xiongxiong5712 at 2012-11-19 09:23:44
你好，我注意到他的輸入文件里k點(diǎn)的設(shè)置是這樣的：
#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計算，是非 ...

z方向在做輸運(yùn)的時候，自動忽略掉k點(diǎn)的設(shè)置，程序會按照1來做的。
但是這個kz值在做開始的SIESTA的時候可能會用到。

贊一下

回復(fù)此樓

8樓2012-11-19 09:42:03

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

xiongxiong5712

新蟲 (小有名氣)

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

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

引用回帖:

1981535樓: Originally posted by zhangguangping at 2012-11-19 09:42:03
z方向在做輸運(yùn)的時候，自動忽略掉k點(diǎn)的設(shè)置，程序會按照1來做的。
但是這個kz值在做開始的SIESTA的時候可能會用到。...

原來這樣，謝謝guangping兄了，我還得多學(xué)習(xí)學(xué)習(xí)。

贊一下

回復(fù)此樓

9樓2012-11-19 10:05:33

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 freewain 的主題更新

返回列表

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

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

[考研] 307求調(diào)劑 +9	冷笙123 2026-03-17	9/450	2026-03-19 22:44 by 學(xué)員8dgXkO
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3	@taotao 2026-03-19	3/150	2026-03-19 14:07 by peike
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu

24小時熱門版塊排行榜

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

» 猜你喜歡

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

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

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

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

查看全部散金貼