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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » transiesta 计算错误
162/2返回列表上一页12
查看: 3618  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 焦虑 +8 水冰月月野兔 2026-03-13 12/600 2026-03-18 15:27 by 咪呜喵呜
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 0703化学调剂 ，六级已过，有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么！ 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 285化工学硕求调剂（081700） +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
关闭
请填处理意见
关闭