24小时热门版块排行榜    

查看: 1107  |  回复: 3

skydark

木虫 (正式写手)

[求助] siesta 有關xcrysden作圖一問

???????D??xcrysden??siesta????????????????
???u?????????? ????siesta??fdf?n???e?^????SaveTotalCharge   .true.
???????selfconsistent???  ???.TOCH ?? .RHO?n??
????rho2xsf???D?n?????.TOCH ?? .RHO?D??.xsf?n??
???????xcrysden?_??xsf?n??

???????n?????_??   xcrysden?@?????????????? ??N???]??  ?????????N?????????????  ??N?|?????]??

???????linux??terminal??vim??????.xsf?n???????_
?l?F?e?^??ж?????:
Specify  SystemLabel (or 'siesta' if none):   Error opening file H^@^@^@U??6?^E?a??6?^E@U??6?^E@U?XV as old formatted

???????? ????????e???e???  ??????N????????????шD?????????


?????????fdf input?n??:

SystemName  Fe
SystemLabel Fe
         
NumberOfSpecies 1
%block ChemicalSpeciesLabel
  1  26   Fe
%endblock ChemicalSpeciesLabel

UserBasis True
      
# K-points

%block kgrid_Monkhorst_Pack
20   0     0     0.0
0    20    0     0.0
0    0     20    0.0
%endblock kgrid_Monkhorst_Pack
   
   
%block ProjectedDensityOfStates
-10.00 10.00 0.200 1000 eV
%endblock ProjectedDensityOfStates

%block BandLines
1  0.000 0.000  0.000 \Gamma
30 0.500 -0.500 0.500 H
30 0.000 0.000  0.500 N
30 0.250 0.250  0.250 P
30 0.000 0.000  0.000 \Gamma
30 0.000 0.000  0.500 N
%endblock BandLines


# UNIT CELL AND ATOMIC POSITIONS

# UNIT CELL
LatticeConstant       2.87 Ang
%block LatticeVectors
-0.50000000    0.5000000000   0.500000000
0.50000000    -0.5000000000   0.500000000
0.50000000     0.5000000000   -0.500000000
%endblock LatticeVectors

# Atomic coordinates
NumberOfAtoms 1
AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
0                    0              0             1
%endblock AtomicCoordinatesAndAtomicSpecies


# General variables

ElectronicTemperature  50 K
MeshCutoff             40. Ry  # 350. Ry
xc.functional            LDA         # Exchange-correlation functional
#xc.authors               PBE
SpinPolarized .true.
SolutionMethod diagon

# MD variables

MD.TypeOfRun        CG
MD.VariableCell    .false.
MD.NumCGsteps       1000
MD.MaxCGDispl       0.05 Ang
MD.MaxForceTol      0.02 eV/Ang
DM.NumberPulay      3
DM.MixingWeight     0.05
DM.Tolerance        1.d-5
MaxSCFIterations    1000


# Output variables

WriteMullikenPop                1
WriteKBands                      .true.
WriteBands                      .true.
SaveRho                         .true.
SaveDeltaRho                    .true.
SaveHS                          .true.
SaveElectrostaticPotential      True
SaveTotalPotential              True
WriteCoorXmol                   .true.
WriteMDXmol                     .true.
WriteMDhistory                  .false.
WriteEigenvalues                .true.
COOP.Write                      .true.
SaveTotalCharge                 .true.
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)

【答案】应助回帖

★
sunyang1988: 金币+1, 谢谢交流 2013-03-07 09:52:46
这个跟你的输入文件没什么关系,如果你的siesta正常的运行结束,那输入就是没有问题的,你可以看一下你的*.TOCH 或者*.RHO的大小,如果里面有数字,说明是没有问题的,而且这两个文件应该不是很小
所以最有可能是你用rho2xsf 转格式的时候出的问题,你没有转好罢了
这只是一场游戏
2楼2013-03-06 22:56:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)

【答案】应助回帖

★
sunyang1988: 金币+1, 谢谢交流 2013-03-07 09:52:54
最重要的一点忘了说了
siesta必须是ifort编译的
要不然rho2xsf 不能用,你的问题八成是这个问题
这只是一场游戏
3楼2013-03-06 23:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo服diyi

新虫 (正式写手)

引用回帖:
3楼: Originally posted by lu_wei07 at 2013-03-06 23:01:17
最重要的一点忘了说了
siesta必须是ifort编译的
要不然rho2xsf 不能用,你的问题八成是这个问题

不用吧

发自小木虫Android客户端
4楼2017-03-09 16:06:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 skydark 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:45 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 11:18 by Equinoxhua
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 11:09 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 10:53 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
请填处理意见