应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】transiesta关闭后可以接着原来的迭代继续算么
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 1312  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

redskywei

木虫 (著名写手)

[交流] 【求助】transiesta关闭后可以接着原来的迭代继续算么 已有4人参与

用transiesta计算时，由于挂载程序过多而内存不足退出，都已经计算到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对应的电荷数是哪部分的电荷，我的这里是个分数，是电子转移到过渡层了么？


还有一个关于电极原胞尺度问题：
我给的电极基矢为：
16.198320    0.000000    0.000000
  0.000000   16.198320    0.000000
  0.000000    0.000000     4.049580

但是不知道怎么到计算散射区时，这里却有个原胞尺度为：
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     
系统也没报错，奇怪。
回复此楼

» 收录本帖的淘帖专辑推荐

SIESTA transport

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
坠入无边的网~~
1楼 2010-05-19 10:55:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★
zzy870720z(金币+2):谢谢交流 2010-05-19 12:13:25
redskywei(金币+8): 2010-05-19 14:52:09
引用回帖:
Originally posted by redskywei at 2010-05-19 03:55:37:
用transiesta计算时，由于挂载程序过多而内存不足退出，都已经计算到transiesta自洽迭代了：
From SIESTA: Efermi=  -4.805
qc0 :   0.000
qcn0 : 456.635
nou,NGL,NGR:        3731         832         83 ...

可以接着自洽迭代啊！直接设置DM.UseSaveDM           T就可以直接接着算，读取上次的*.TSDE密度矩阵，（这个是transiesta的密度矩阵），如果已经算出*.GF来了，重新计算的时候不想在计算这个了，因为这个花不少的时间。对同一个体系的话，没必要再计算一次了。这个时候就设置一下TS.CalcGF  F,这个时候就不会再计算了。你的第二个问题，可能是你对这个还不是很熟悉，或者以前没有注意过。siesta和transiesta内部执行的时候总是用的原子单位。所以那是以原子单位输出的原胞矢量。你是用的Ang做的单位吧。对于电荷的问题，我以前没有关注过，你可以到mail-list里去问一下。你要是知道了，记得告诉我奥。呵呵呵，以后多交流！

[ Last edited by zhangguangping on 2010-5-19 at 04:37 ]
赞一下(1人) 回复此楼
弘德明志博学笃行
2楼2010-05-19 11:35:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)
引用回帖:
Originally posted by zhangguangping at 2010-05-19 11:35:55:

可以接着自洽迭代啊！直接设置DM.UseSaveDM           T就可以直接接着算，读取上次的*.TSDE密度矩阵，（这个是transiesta的密度矩阵），如果已经算出*.GF来了，重新计算的时候不想在计算这个了，因为这个花不少 ...

呵呵，多谢了。金币不多，笑纳。
我问问看吧，有消息发出来大家都看看。
赞一下 回复此楼
坠入无边的网~~
3楼2010-05-19 14:54:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mayim2008

金虫 (正式写手)
redskywei(金币+2):小张哦，记住你的名字了。呵呵 2010-05-19 21:46:45
嘿嘿，SIESTA可以中途停掉是可以接着算的，好像好加个什么命令，师兄去问问王师兄，我的笔记放在家里了，命令都忘记了~TRSIESTA应该差不多吧~
赞一下 回复此楼
4楼2010-05-19 18:34:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bumubutai

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
1005721楼: Originally posted by zhangguangping at 2010-05-19 11:35:55
可以接着自洽迭代啊！直接设置DM.UseSaveDM           T就可以直接接着算，读取上次的*.TSDE密度矩阵，（这个是transiesta的密度矩阵），如果已经算出*.GF来了，重新计算的时候不想在计算这个了，因为这个花不少的 ...

嗨  您好~~~ 看了这个贴我还有个问题。  我的output和他有点像。 qsol 和qtol 的值相差很大。  我不知道是什么原因  也不知道该怎样改。 我查maillist 上说  这两个值应该比较接近的。  求指点。  多谢了。

我做的是一个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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)
引用回帖:
1273904楼: Originally posted by bumubutai at 2012-10-01 23:12:09
嗨  您好~~~ 看了这个贴我还有个问题。  我的output和他有点像。 qsol 和qtol 的值相差很大。  我不知道是什么原因  也不知道该怎样改。 我查maillist 上说  这两个值应该比较接近的。  求指点。  多谢了。

我做 ...

我很久没算这个了，电荷的问题主要是转移，你自己检查一下是否会在电极和中心之间发生转移，看这个输入我找不出啥问题。
赞一下 回复此楼
坠入无边的网~~
6楼2012-10-04 10:11:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 redskywei 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +3 种圣赐 2026-03-29 3/150 2026-03-29 12:06 by longlotian
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 318一志愿吉林大学生物与医药 求调剂 +4 笃行致远. 2026-03-28 4/200 2026-03-29 11:55 by longlotian
[考研] 一志愿双一流机械285分求调剂 +3 幸运的三木 2026-03-29 4/200 2026-03-29 10:46 by 想上岸了呜呜
[考研] 070300化学354求调剂 +5 101次希望 2026-03-28 5/250 2026-03-29 10:36 by 唐沐儿
[考研] 070305高分子化学与物理 304分求调剂 +8 c297914 2026-03-28 8/400 2026-03-29 10:23 by ms629
[考研] 0703化学调剂，求导师收 +9 天天好运来上岸� 2026-03-24 10/500 2026-03-28 22:17 by chemzp
[考研] 085600 286分 材料求调剂 +7 麻辣鱿鱼 2026-03-27 8/400 2026-03-28 12:17 by zllcz
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 07化学280分求调剂 +10 722865 2026-03-23 10/500 2026-03-27 15:51 by Plutoqq
[考研] 安徽大学专硕生物与医药专业(086000)324分，英语已过四六级，六级521，求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 085600材料与化工调剂 +10 A-哆啦Z梦 2026-03-23 16/800 2026-03-27 15:13 by caszguilin
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 321求调剂 +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
[考研] 中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培 +5 YangTyu1 2026-03-26 5/250 2026-03-26 18:27 by 猫咪猫咪呀
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 296求调剂 +4 汪！？！ 2026-03-25 7/350 2026-03-25 16:41 by 汪！？！
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
信息提示
关闭
请填处理意见
关闭