24小时热门版块排行榜    

查看: 923  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

newton3915

铁杆木虫 (正式写手)

[求助] siesta 优化石墨结构

请教各位,小弟初学siesta,尝试着优化石墨的结构,发现程序运行到如下所示没多少就停了,不知什么问题
siesta: ******************** Simulation parameters ****************************
siesta:
siesta: The following are some of the parameters of the simulation.
siesta: A complete list of the parameters used, including default values,
siesta: can be found in file out.fdf
siesta:
redata: Non-Collinear-spin run           =     F
redata: SpinPolarized (Up/Down) run      =     F
redata: Number of spin components        =     1
redata: Long output                      =     F
redata: Number of Atomic Species         =        1
redata: Charge density info will appear in .RHO file
---------------------------------------------------------------------------------------------------------------------------
下面是我的fdf文件,不知道哪里出错了,希望大家点拨一下。不胜感激!
# General System Descriptors

SystemName                graphite geoopt
SystemLabel                graphite_opt

NumberOfAtoms                4
NumberOfSpecies                1
%block ChemicalSpeciesLabel
1  6  C
%endblock ChemicalSpeciesLabel

PAO.BasisSize                DZP
PAO.EnergyShift                0.02 Ry

# Structure and K-sampling

LatticeConstant                2.46 Ang
%block LatticeVectors
  1.000  0.000  0.000
-0.500  0.866  0.000
  0.000  0.000  2.725
%endblock LatticeVectors

AtomicCoordinatesFormat        Fractional
AtomicCoorFormatOut        Ang
%block AtomicCoordinatesAndAtomicSpecies
  0.000  0.000  0.000  1
  0.333  0.667  0.000  1
  0.000  0.000  0.500  1
  0.667  0.333  0.500  1
%endblock AtomicCoordinatesAndAtomicSpecies

%block kgrid_Monkhorst_Pack
  12  0  0  0.5
   0 12  0  0.5
   0  0  4  0.5
%endblock kgrid_Monkhorst_Pack

# DFT

XC.functional                LDA
XC.authors                CA
SpinPolarized                false
MaxSCFIterations        50     
DM.MixingWeight                0.25
DM.NumberPulay                3
DM.Tolerance                1.d-4
MeshCutoff                100 Ry
SolutionMethod                diagon
ElectronicTemperature        300.0 K

# MD and Relaxations

MD.TypeOfRun                cg
MD.VariableCell                true
MD.MaxForceTol                0.04 eV/Ang
MD.MaxStressTol                1.0 GPa
MD.NumCGsteps                200
MD.MaxCGDispl                0.2 Bohr

# Output options

WriteCoorInitial        true
WriteCoorStep                true
WriteForces                true
WriteKpoints                false
WriteEigenvalues        false
Writekbands                false
WriteBands                false
WriteMullikenPop        false
WriteWaveFunction        false
WriteCoorXmol                false
WriteCoorCerius                false
WriteMDCoorXmol                false
WriteMDhistory                false
WriteMDXmol                false
WriteDM                        true

# Options for saving/reading information

DM.UseSaveDM                false
MD.UseSaveXV                false
MD.UseSaveCG                false
SaveRho                        false
SaveDeltaRho                false
SaveElectrostaticPotential        false
SaveTotalPotential        false
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

newton3915

铁杆木虫 (正式写手)

问题已解决,WriteMullikenPop 是integer不是logical。
3楼2014-11-24 22:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

newton3915

铁杆木虫 (正式写手)

顶一下

[ 发自小木虫客户端 ]
2楼2014-11-24 15:28:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 7/350 2026-09-03 13:20 by T0rGB46095mJ
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 13:20 by T0rGB46095mJ
[基金申请] 国社科又开始会评了,不知道这次命运如何 +10 雨打竹帘 2026-08-30 14/700 2026-09-03 13:17 by qsd10086
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:05 by T0rGB46095mJ
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:03 by T0rGB46095mJ
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
请填处理意见