| 查看: 534 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求siesta的mesh_cutoff和 kgrid测试脚本【需要继续讨论】
|
||
|
求siesta的mesh_cutoff和 kgrid测试脚本 [ Last edited by 贺仪 on 2012-5-11 at 08:48 ] |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 从outcar中提取原子坐标的脚本
已经有14人回复
- 请问南京或常州哪里有四探针电阻测试仪----测试薄膜
已经有11人回复
- 原子力显微镜的样品用什么拿到测试的地方去？自己实验室和那个地方离得不近。
已经有7人回复
- 采用PBS脚本第一次提交任务，算了1步就被杀死，有错误提示但不知道原因？
已经有19人回复
- 如何测试材料的脆性
已经有3人回复
- 收敛测试，这种情况我该继续测试么,简直就是跳跃式的嘛
已经有24人回复
- siesta-3.1编译成功，但有几个问题求助
已经有7人回复
- IE 总是弹出这个烦人的脚本错误
已经有5人回复
- 【求助】K-point或 energy cut-off 测试
已经有12人回复
- 【求助】SIESTA 中ATOM生成赝势测试后如何看赝势的好坏
已经有4人回复
- 【求助】请教MeshCutOff，PAO.EnergyShift和KgridCutoff等具体物理意义和取值依据
已经有6人回复
- 【求助】VMD 中cutoff的设置
已经有3人回复
3楼2012-05-11 07:47:17
guohuazhong
至尊木虫 (职业作家)
- 1ST强帖: 1
- 应助: 52 (初中生)
- 贵宾: 0.032
- 金币: 35131.5
- 散金: 5
- 红花: 19
- 沙发: 1
- 帖子: 3924
- 在线: 885.6小时
- 虫号: 550098
- 注册: 2008-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
love5264: 金币+30, ★★★★★最佳答案, 学习了！谢谢！！！ 2012-05-11 07:41:47
感谢参与，应助指数 +1
love5264: 金币+30, ★★★★★最佳答案, 学习了！谢谢！！！ 2012-05-11 07:41:47
|
看看下面这个算I-V曲线的脚本就明白怎么弄了 #!/bin/sh # #BSUB -a intelmpi # Bu satiri silmeyin #BSUB # LSF stdout dosyasi ismi .%J.out #BSUB # LSF stderr dosyasi ismi .%J.err #BSUB -q salim.ciraci.q # kuyruk ismi #BSUB -m karadeniz_e1 # -m unam_bladez OR -m n_anadolu #BSUB -n 8 # islemci sayisi 4 8 16 32 #BSUB -R "span[ptile=8]" # #BSUB -J a8-b2-yb # LSF is ismi ncpu=8 ################################################################################ cp ../step-1-electrodes/elec.TSHS . # copy TSHS file from step-1 mkdir cont # read the comment at the end of this script. for i in `seq -w 0.0 0.1 1.8` do cp -r cont $i cd $i cp ../*.psf . cp ../elec.TSHS . cp ../POSITIONS.fdf . cat > INPUT.fdf < SystemLabel scat %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 3 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 1000 #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 50 meV PAO.SplitNorm 0.1500000000 PAO.BasisType split PAO.BasisSize SZP WriteMullikenPop 0 WriteBands F SaveRho F SaveElectrostaticPotential F SaveTotalPotential F WriteCoorXmol T %block ChemicalSpeciesLabel 1 6 C-GGA %endblock ChemicalSpeciesLabel %include POSITIONS.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 $i 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 ./elec.TSHS TS.NumUsedAtomsLeft 16 TS.BufferAtomsLeft 16 # RIGHT ELECTRODE TS.HSFileRight ./elec.TSHS TS.NumUsedAtomsRight 16 TS.BufferAtomsRight 16 ================================================== ================================================== EOF mpirun.lsf -np $ncpu /XXX/transiesta < INPUT.fdf > OUTPUT.fdf cd .. rm -rf cont mkdir cont cp ./$i/scat.TSDE ./$i/scat.TSHS ./$i/scat.DM cont # copy these files for continuation of the next bias step. done |
2楼2012-05-10 21:11:11
25