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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】参数对计算结果的影响
61/1返回列表
查看: 1047  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yashumimi

银虫 (小有名气)

[交流] 【求助】参数对计算结果的影响 已有1人参与

发现改变计算的一个参数，对结果都有很大的影响，不太明白。

比如改变MaxSCFIterations值，计算都是收敛的，在同一偏压下，电流值增加了4个数量级。
比如PAO.BasisSize类型DZP为SZP，在同一偏压下，电流值就增加了5个数量级。
比如，我将TS.BufferatomsLeft/Right设置为和左右电极数一样，增加原子数，散射区未改变。在同一偏压下，电流值就增加了5个数量级。

这种情况正常么？
以下是我的输入文件：

SystemName                                                                                                                     
SystemLabel   

%block kgrid_Monkhorst_Pack
  20  0   0    0.0
  0   1   0    0.0
  0   0   20    0.0
%endblock kgrid_Monkhorst_Pack


xc.functional           GGA                     
xc.authors              PBE                     

MeshCutoff              300.00000000 Ry
SolutionMethod Transiesta                       
OccupationFunction      MP                     
ElectronicTemperature   300 K                  

SpinPolarized           F                       
FixSpin                 F                       

MaxSCFIterations       100            
#DM.Tolerance           1.d-5                  
DM.NumberPulay          6                       
#DM.NumberBroyden        0                       
DM.MixingWeight         0.1000000000            
DM.OccupancyTolerance    0.1000000000E-11        
#DM.NumberKick           0                       
#DM.KickMixingWeight     0.5000000000            
DM.Tolerance            0.1000000000E-03        

UseSaveData             T
                  
MD.NumCGsteps           0                        
MD.TypeOfRun            CG                       
MD.VariableCell         F                        
MD.MaxCGDispl           0.2000000000  Bohr      
MD.MaxForceTol          0.05 eV/Ang              
#MD.MaxStressTol        0.0001 eV/Ang**3        


Diag.ParallelOverK      F                                                
PAO.EnergyShift         20 meV                    
PAO.SplitNorm           0.1500000000              
PAO.BasisType           split                     
        
                          

WriteMullikenPop                0
WriteBands                      F
SaveRho                         F
SaveElectrostaticPotential      F
SaveTotalPotential              F
WriteCoorXmol                   T


                                                                                                                              
#%block ChemicalSpeciesLabel
#    1  6 C-GGA
#%endblock ChemicalSpeciesLabel

%include POSITION.fdf


# Transiesta information
SolutionMethod Transiesta
# GENGF OPTIONS
TS.ComplexContour.Emin        -28 eV
TS.ComplexContour.NPoles       16
TS.ComplexContour.NCircle      16
TS.ComplexContour.NLine        10
# BIAS OPTIONS
TS.biasContour.NumPoints       10


# TS OPTIONS
TS.Voltage 1.7  eV

# TBT OPTIONS
TS.TBT.Emin  -3 eV
TS.TBT.Emax  +3 eV
TS.TBT.NPoints  500
TS.TBT.NEigen  3
TS.TBT.Eta        0.000001 Ry

# Write electrode hamiltonian
TS.SaveHS          .true.
TS.SaveLead        .true.



# LEFT ELECTRODE
TS.HSFileLeft  ./g-e.TSHS
TS.NumUsedAtomsLeft   16
TS.BufferAtomsLeft    16

# RIGHT ELECTRODE
TS.HSFileRight  ./g-e.TSHS
TS.NumUsedAtomsRight  16
TS.BufferAtomsRight   16
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-11 15:50:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

隆长江

银虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
楼主，你那个问题现在解决了没，我也遇到相似的问题了。可否交流一下这是什么原因呢？
赞一下(1人) 回复此楼
新手来报到，请君多关照！
2楼2011-05-24 17:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yashumimi

银虫 (小有名气)
引用回帖:
Originally posted by 隆长江 at 2011-05-24 17:18:29:
楼主，你那个问题现在解决了没，我也遇到相似的问题了。可否交流一下这是什么原因呢？

你的是什么问题。。。。好久了我都忘记了 我应该是改变结构了吧 大致参数变动的不多。。。
赞一下 回复此楼
3楼2011-05-25 18:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

隆长江

银虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
我的是在0偏压时，算出来的透射谱在10的负5次方左右。我没有去改变电压和左右电极，我改变PAO.BasisSize时没有什么变化。电流还没有算过，估计也会很小，我想你可能也遇到的是差不多的情况，所以问一下。
赞一下(1人) 回复此楼
新手来报到，请君多关照！
4楼2011-05-25 19:15:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yashumimi

银虫 (小有名气)
引用回帖:
Originally posted by 隆长江 at 2011-05-25 19:15:50:
我的是在0偏压时，算出来的透射谱在10的负5次方左右。我没有去改变电压和左右电极，我改变PAO.BasisSize时没有什么变化。电流还没有算过，估计也会很小，我想你可能也遇到的是差不多的情况，所以问一下。

我的问题貌似已经放在一边没有管了。。。0电压的电流不就是0么
赞一下 回复此楼
5楼2011-06-17 20:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

隆长江

银虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
电压是零，我是说在.AVTRANS中的T(E)这一列数字。
赞一下(1人) 回复此楼
新手来报到，请君多关照！
6楼2011-09-08 20:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yashumimi 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 如何实现卤原子转化 10+3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 3/150 2026-05-16 12:49 by vcdazktkjx
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 12:24 by h3oerqvkv9
[硕博家园] 考博自荐 +3 科研狗111 2026-05-13 4/200 2026-05-16 11:45 by 科研狗111
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 11:36 by h3oerqvkv9
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 cjf4bx70cj 2026-05-14 6/300 2026-05-16 11:16 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[文学芳草园] 裁员滚滚，退居二线 +4 J_wei 2026-05-10 4/200 2026-05-16 10:52 by zh10246
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +11 dd921ww 2026-05-12 13/650 2026-05-16 09:16 by Howard28
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 04:44 by k37jurhrau
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-16 01:49 by k37jurhrau
[教师之家] 上海大学实验技术岗位非升即走 +5 嘻嘻哈哈乐呵呵 2026-05-15 5/250 2026-05-16 00:17 by caiyun
[基金申请] 青C资助名额大幅增加！ +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
[考博] 现在不知道怎么办，感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
关闭
请填处理意见
关闭