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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】transiesta关闭后可以接着原来的迭代继续算么
51/1返回列表
查看: 1308  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +4 CX 330 2026-03-21 4/200 2026-03-21 23:48 by ms629
[考研] 269专硕求调剂 +5 金恩贝 2026-03-21 5/250 2026-03-21 22:37 by zhyzzh
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 15/750 2026-03-21 10:16 by 脱颖而出
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿华中农业071010，总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
关闭
请填处理意见
关闭