 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1303 | 回復(fù): 5 | |||||
redskywei木蟲 (著名寫手)
|
[交流]
【求助】transiesta關(guān)閉后可以接著原來的迭代繼續(xù)算么 已有4人參與
|
|
用transiesta計算時，由于掛載程序過多而內(nèi)存不足退出，都已經(jīng)計算到transiesta自洽迭代了： From SIESTA: Efermi= -4.805 qc0 : 0.000 qcn0 : 456.635 nou,NGL,NGR: 3731 832 832 TranSiesta: Qsol,Qtot: 483.480 957.000 1.000 transiesta: TSiscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) transiesta: 1 -25796.9060 -25694.0764 -25694.0764 85.4730 -4.8048 timer: Routine,Calls,Time,% = IterSCF 44 118292.741 80.25 elaps: Routine,Calls,Wall,% = IterSCF 44 59371.237 80.29 timer: Routine,Calls,Time,% = TS 1 28704.656 19.47 elaps: Routine,Calls,Wall,% = TS 1 14364.325 19.43 TranSiesta: Qsol,Qtot: 534.939 957.000 1.000 transiesta: 2 -88127.8734 12349.9694 12349.9694228.7030 -4.8048 TranSiesta: Qsol,Qtot: 488.859 957.000 1.000 transiesta: 3 -56270.7540 -21245.4598 -21245.4598 99.3472 -4.8048 TranSiesta: Qsol,Qtot: 488.432 957.000 1.000 transiesta: 4 -55294.0759 -21317.6327 -21317.6327 55.5982 -4.8048 TranSiesta: Qsol,Qtot: 479.158 957.000 1.000 transiesta: 5 -53117.3585 -21539.4531 -21539.4531 41.3231 -4.8048 TranSiesta: Qsol,Qtot: 470.471 957.000 1.000 transiesta: 6 -52472.9078 -21585.8679 -21585.8679 46.3186 -4.8048 TranSiesta: Qsol,Qtot: 515.301 957.000 1.000 transiesta: 7 -55305.2639 -18966.4175 -18966.4175 84.6345 -4.8048 TranSiesta: Qsol,Qtot: 463.658 957.000 1.000 transiesta: 8 -47134.7254 -21236.5423 -21236.5423 51.2251 -4.8048 TranSiesta: Qsol,Qtot: 456.971 957.000 1.000 transiesta: 9 -44427.7347 -21318.3401 -21318.3401 46.8924 -4.8048 TranSiesta: Qsol,Qtot: 444.560 957.000 1.000 transiesta: 10 -43802.3646 -21233.9120 -21233.9120166.2342 -4.8048 TranSiesta: Qsol,Qtot: 459.526 957.000 1.000 transiesta: 11 -43144.4210 -21205.0722 -21205.0722 40.8517 -4.8048 TranSiesta: Qsol,Qtot: 438.813 957.000 1.000 transiesta: 12 -41201.9274 -20950.2247 -20950.2247 56.0684 -4.8048 TranSiesta: Qsol,Qtot: 436.696 957.000 1.000 transiesta: 13 -37944.0731 -21239.0155 -21239.0155 87.3010 -4.8048 TranSiesta: Qsol,Qtot: 441.258 957.000 1.000 transiesta: 14 -36116.9398 -21165.2910 -21165.2910 55.2350 -4.8048 TranSiesta: Qsol,Qtot: 449.885 957.000 1.000 transiesta: 15 -33992.3977 -21101.5238 -21101.5238 85.3353 -4.8048 TranSiesta: Qsol,Qtot: 458.481 957.000 1.000 transiesta: 16 -33016.1655 -21078.3628 -21078.3628104.3443 -4.8048 TranSiesta: Qsol,Qtot: 468.988 957.000 1.000 transiesta: 17 -32348.3750 -21035.2013 -21035.2013101.3979 -4.8048 TranSiesta: Qsol,Qtot: 482.396 957.000 1.000 transiesta: 18 -30486.7028 -21057.5577 -21057.5577 59.6814 -4.8048 TranSiesta: Qsol,Qtot: 482.013 957.000 1.000 transiesta: 19 -29276.6606 -21011.9338 -21011.9338 67.5077 -4.8048 TranSiesta: Qsol,Qtot: 510.079 957.000 1.000 transiesta: 20 -26873.0004 -20920.4854 -20920.4854 59.2618 -4.8048 另外有個小問題想問，哪個Qsol對應(yīng)的電荷數(shù)是哪部分的電荷，我的這里是個分?jǐn)?shù)，是電子轉(zhuǎn)移到過渡層了么？ 還有一個關(guān)于電極原胞尺度問題： 我給的電極基矢為： 16.198320 0.000000 0.000000 0.000000 16.198320 0.000000 0.000000 0.000000 4.049580 但是不知道怎么到計算散射區(qū)時，這里卻有個原胞尺度為： Begin LEFT unit cell: 30.6104 0.0000 0.0000 0.0000 30.6104 0.0000 0.0000 0.0000 7.6526 Size: 692224 Efermi: 0.0000000000000000 GF: Spin number -> 1 Got LEFT Electrode GF Done Begin RIGHT unit cell: 30.6104 0.0000 0.0000 0.0000 30.6104 0.0000 0.0000 0.0000 7.6526 Size: 692224 Efermi: -0.0000000000000000 系統(tǒng)也沒報錯，奇怪。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 SIESTA | transport |
» 猜你喜歡
317求調(diào)劑
已經(jīng)有12人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有7人回復(fù)
-
324求調(diào)劑
已經(jīng)有5人回復(fù)
-
288求調(diào)劑
已經(jīng)有16人回復(fù)
-
317求調(diào)劑
已經(jīng)有9人回復(fù)
-
287求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有4人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有4人回復(fù)
-
085600材料與化工
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于electrostatic potential
已經(jīng)有5人回復(fù)
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
★ ★
zzy870720z(金幣+2):謝謝交流 2010-05-19 12:13:25
redskywei(金幣+8): 2010-05-19 14:52:09
zzy870720z(金幣+2):謝謝交流 2010-05-19 12:13:25
redskywei(金幣+8): 2010-05-19 14:52:09
|
可以接著自洽迭代啊！直接設(shè)置DM.UseSaveDM T就可以直接接著算，讀取上次的*.TSDE密度矩陣，（這個是transiesta的密度矩陣），如果已經(jīng)算出*.GF來了，重新計算的時候不想在計算這個了，因?yàn)檫@個花不少的時間。對同一個體系的話，沒必要再計算一次了。這個時候就設(shè)置一下TS.CalcGF F,這個時候就不會再計算了。你的第二個問題，可能是你對這個還不是很熟悉，或者以前沒有注意過。siesta和transiesta內(nèi)部執(zhí)行的時候總是用的原子單位。所以那是以原子單位輸出的原胞矢量。你是用的Ang做的單位吧。對于電荷的問題，我以前沒有關(guān)注過，你可以到mail-list里去問一下。你要是知道了，記得告訴我奧。呵呵呵，以后多交流！ [ Last edited by zhangguangping on 2010-5-19 at 04:37 ] |
2樓2010-05-19 11:35:55
redskywei
木蟲 (著名寫手)
- 應(yīng)助: 14 (小學(xué)生)
- 金幣: 2935.1
- 散金: 57
- 紅花: 4
- 帖子: 1413
- 在線: 300小時
- 蟲號: 222757
- 注冊: 2006-03-21
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-05-19 14:54:10
4樓2010-05-19 18:34:19
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
嗨 您好~~~ 看了這個貼我還有個問題。 我的output和他有點(diǎn)像。 qsol 和qtol 的值相差很大。 我不知道是什么原因 也不知道該怎樣改。 我查maillist 上說 這兩個值應(yīng)該比較接近的。 求指點(diǎn)。 多謝了。 我做的是一個nano電容的計算。 我把input貼在這里 供參考。 # SPECIES AND BASIS # Number of species NumberOfSpecies 3 %block ChemicalSpeciesLabel 1 12 Mg_LDA 2 8 O_LDA 3 79 Au_LDA %endblock ChemicalSpeciesLabel %block PAO.BasisSizes Mg_LDA SZP O_LDA SZP Au_LDA SZP %endblock PAO.BasisSizes PAO.BasisType split ================================================== ================================================== # K-points %block kgrid_Monkhorst_Pack ================================================== ================================================== # UNIT CELL LatticeConstant 4.31132 Ang %block LatticeVectors < lat_AuMgOAu.dat # Atomic coordinates NumberOfAtoms 36 AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies < coords_AuMgOAu.dat ================================================== ================================================== # General variables ElectronicTemperature 300 K OccupationFunction MP MeshCutoff 250. Ry xc.functional GGA # Exchange-correlation functional xc.authors PBE SpinPolarized .false. ================================================== ================================================== SolutionMethod Transiesta ================================================== ================================================== # SCF variables DM.MixingWeight 0.01 DM.NumberPulay 5 MaxSCFIterations 300 # Maximum number of SCF iter DM.Tolerance 1.d-4 # Tolerance in maximum difference DM.UseSaveDM true # to use continuation files Diag.DivideAndConquer no ================================================== ================================================== # MD variables MD.FinalTimeStep 1 MD.TypeOfRun CG MD.NumCGsteps 000 MD.UseSaveXV .true. ================================================= ================================================= # Output variables WriteMullikenPop 1 WriteBands .false. SaveRho .true. SaveDeltaRho .true. SaveHS .false. SaveElectrostaticPotential .true. SaveIonicCharge .true. SaveTotalCharge .true. WriteCoorXmol .true. WriteMDXmol .true. #animation created WriteMDhistory .false. WriteEigenvalues yes ================================================== ================================================== # Parallel variables Diag.ParallelOverK .true. ProcessorY 1 ================================================== ================================================== # Transiesta information # GF OPTIONS TS.ComplexContour.Emin -30.0 eV TS.ComplexContour.NPoles 06 TS.ComplexContour.NCircle 30 TS.ComplexContour.NLine 10 # BIAS OPTIONS TS.biasContour.NumPoints 15 TS.UpdateDMCROnly .false. # TS OPTIONS TS.Voltage 0.000000 eV #Write hamiltonian TS.SaveHS .true. # LEFE ELECTRODE TS.HSFileLeft 'elec.TSHS' TS.NumUsedAtomsLeft 4 # RIGHT ELECTROD TS.HSFileRight 'elec.TSHS' TS.NumUsedAtomsRight 4 不勝感激 不勝感激~~~ |
5樓2012-10-01 23:12:09
redskywei
木蟲 (著名寫手)
- 應(yīng)助: 14 (小學(xué)生)
- 金幣: 2935.1
- 散金: 57
- 紅花: 4
- 帖子: 1413
- 在線: 300小時
- 蟲號: 222757
- 注冊: 2006-03-21
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
6樓2012-10-04 10:11:24
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 redskywei 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +3 | 77chaselx 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +3 | shxchem 2026-03-20 | 4/200 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+7
|
.6lL 2026-03-18 | 7/350 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研]
11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|