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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » transiesta 计算错误
162/2返回列表上一页12
查看: 4078  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

天地逍遥

铁虫 (初入文坛)
引用回帖:
279761楼: Originally posted by zhangguangping at 2011-08-21 13:55:54
附带上一个套输入文件：
体系结构：
24/97/529624_1313905891_601.jpg
说明：左右电极相同，是一个3*4的Au(111)的ABC结构，三层，每层12个原子。中间是散射区，金连有点长而已。

电极的输入文件：
SystemNa ...

请问，再用这些数据之前，你是用CG＝100经过coordinate optimization的嘛？也是用这36个原子，输入文件也一样，就CG变了变是嘛？
赞一下 回复此楼
11楼2012-07-09 04:50:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★
fzx2008: 金币+1, 谢谢回帖！ 2012-07-09 08:36:56
引用回帖:
4624647楼: Originally posted by 天地逍遥 at 2012-07-09 04:50:49
请问，再用这些数据之前，你是用CG＝100经过coordinate optimization的嘛？也是用这36个原子，输入文件也一样，就CG变了变是嘛？...

这儿的CG=0不再进行优化了。优化结构是之前的事。
赞一下 回复此楼
弘德明志博学笃行
12楼2012-07-09 08:07:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天地逍遥

铁虫 (初入文坛)
引用回帖:
4624671楼: Originally posted by zhangguangping at 2012-07-09 08:07:17
这儿的CG=0不再进行优化了。优化结构是之前的事。...

对，这个结构就是优化完的结果吧，你之前的优化输入文件还有嘛，可否给我看下，我设定的CG＝100，然后Kpoints是10x10x20，16个原子优化了好久都没有优化完，而且原子坐标变化都好大。。。正发愁呢。。。

tangshengjie61@gmail.com
赞一下 回复此楼
13楼2012-07-09 08:55:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
4624693楼: Originally posted by 天地逍遥 at 2012-07-09 08:55:46
对，这个结构就是优化完的结果吧，你之前的优化输入文件还有嘛，可否给我看下，我设定的CG＝100，然后Kpoints是10x10x20，16个原子优化了好久都没有优化完，而且原子坐标变化都好大。。。正发愁呢。。。

tangsh ...

你的k点设置的有点大。如果带着部分Au优化的话，k设置一个6x6x1或者6x6x1就行，并且z方向不用使用周期性边界条件。如果你做的是固定电极的距离的优化的话，z方向使用周期性边界条件也行，不过会加大计算量的。
赞一下(1人) 回复此楼
弘德明志博学笃行
14楼2012-07-09 09:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天地逍遥

铁虫 (初入文坛)
引用回帖:
4624755楼: Originally posted by zhangguangping at 2012-07-09 09:53:27
你的k点设置的有点大。如果带着部分Au优化的话，k设置一个6x6x1或者6x6x1就行，并且z方向不用使用周期性边界条件。如果你做的是固定电极的距离的优化的话，z方向使用周期性边界条件也行，不过会加大计算量的。...

那个K点不是要从小到大不断测试，直到能量趋于稳定嘛？难道K点不用多次测试，去能量稳定的最少值，可以直接取嘛？做Au的时候我是6x6x1一直试到10x10x20都快把我算死了。
如何设置Z方向不是周期性边界条件呀，默认的是周期性的吧？
这是我做的Cu电极优化前的输入文件，请指教呀，哪里可以纠正一下呢？
SystemName  bulk
SystemLabel bulk

==================================================
==================================================
# SPECIES AND BASIS

# Number of species
NumberOfSpecies 1
%block ChemicalSpeciesLabel
  1  29 Cu
%endblock ChemicalSpeciesLabel

PAO.BasisType    split
PAO.BasisSize SZP
PAO.EnergyShift 50 meV

==================================================
==================================================
# K-points

%block kgrid_Monkhorst_Pack
10   0   0   0.0
0   10   0   0.0
0   0   20  0.5
%endblock kgrid_Monkhorst_Pack
==================================================
==================================================
# UNIT CELL
LatticeConstant       3.615 Ang
%block LatticeVectors
1   0  0
0  1  0
0   0  1
%endblock LatticeVectors


# Atomic coordinates
NumberOfAtoms 16
AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
  2.5562  0.0000  0.0000       1
  1.2781  2.2137  0.0000       1
  3.8343  2.2137  0.0000       1
  0.0000  0.0000  0.0000       1
  1.2781  0.7379  2.0871       1
  3.8343  0.7379  2.0871       1
  0.0000  4.4275  0.0000       1
  2.5562  4.4275  0.0000       1
  2.5562  2.9516  2.0871       1
  1.2781  5.1654  2.0871       1
  3.8343  5.1654  2.0871       1
  0.0000  2.9516  2.0871       1
  0.0000  1.4758  4.1742       1
  1.2781  3.6895  4.1742       1
  3.8343  3.6895  4.1742       1
  2.5562  1.4758  4.1742       1
%endblock AtomicCoordinatesAndAtomicSpecies

==================================================
==================================================
# General variables

ElectronicTemperature  300 K
OccupationFunction MP
MeshCutoff           350. Ry
xc.functional         LDA           # Exchange-correlation functional
xc.authors            CA
SpinPolarized .false.

==================================================
==================================================
SolutionMethod Diagon


==================================================
==================================================
# SCF variables

DM.MixSCF1   T
MaxSCFIterations      300           # Maximum number of SCF iter
DM.MixingWeight       0.01          # New DM amount for next SCF cycle
DM.Tolerance          1.d-4         # Tolerance in maximum difference
DM.UseSaveDM          true          # to use continuation files
DM.NumberPulay         3
DM.NumberKick         50
DM.KickMixingWeight   0.05
Diag.DivideAndConquer     no

==================================================
==================================================
# MD variables

MD.FinalTimeStep         1
MD.TypeOfRun                 CG
MD.NumCGsteps             100
MD.UseSaveXV              .true.
MD.VariableCell                .true.

=================================================
=================================================


==================================================
==================================================
# Output variables

WriteMullikenPop                1
WriteBands                      .false.
#SaveRho                         .true.
#SaveDeltaRho                    .true.
SaveHS                          .false.
#SaveElectrostaticPotential      .true.
#SaveTotalPotential              .true.
#SaveIonicCharge                        .true.
#SaveTotalCharge                        .true.
WriteCoorXmol                   .true.
WriteMDXmol                     .true.  #animation created
WriteMDhistory                  .false.
#WriteEigenvalues                yes

==================================================
==================================================
# Parallel variables

Diag.ParallelOverK      true

==================================================
==================================================
我看到我的xc.functional和xc.authors与你的设置不太一样，这个有关系嘛？
赞一下 回复此楼
15楼2012-07-09 12:25:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天地逍遥

铁虫 (初入文坛)
引用回帖:
4624755楼: Originally posted by zhangguangping at 2012-07-09 09:53:27
你的k点设置的有点大。如果带着部分Au优化的话，k设置一个6x6x1或者6x6x1就行，并且z方向不用使用周期性边界条件。如果你做的是固定电极的距离的优化的话，z方向使用周期性边界条件也行，不过会加大计算量的。...

您好，请问，就是如何设置非周期性边界呀，默认的是周期性的呀，我找了半天的命令都没找到。我是想先优化Cu电极，然后算电极，再优化整个体系（电极距离固定），然后再算整个体系的transiesta，我这样可以z方向非周期，然后6x6x1嘛？还有就是怎么设置呀。。。唉，现在非常迷惑中。
赞一下 回复此楼
16楼2012-07-10 00:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 guohuazhon 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:38 by G6APbkg8SA6w
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:28 by G6APbkg8SA6w
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭