 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1316 | 回復(fù): 5 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||
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 |
» 猜你喜歡
求調(diào)劑
已經(jīng)有6人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有4人回復(fù)
-
327求調(diào)劑
已經(jīng)有4人回復(fù)
-
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于electrostatic potential
已經(jīng)有5人回復(fù)
4樓2010-05-19 18:34:19
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 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來了，重新計算的時候不想在計算這個了，因為這個花不少的時間。對同一個體系的話，沒必要再計算一次了。這個時候就設(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
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
嗨 您好~~~ 看了這個貼我還有個問題。 我的output和他有點像。 qsol 和qtol 的值相差很大。 我不知道是什么原因 也不知道該怎樣改。 我查maillist 上說 這兩個值應(yīng)該比較接近的。 求指點。 多謝了。 我做的是一個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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工328分調(diào)劑 +7 | 。，。，。，。i 2026-03-23 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9 | 天天好運來上岸?/a> 2026-03-24 | 10/500 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|