24小时热门版块排行榜    

查看: 676  |  回复: 3

yangxuezhang

金虫 (正式写手)

[求助] siesta计算问题

我算了一个Ti3SiC2的机构,一个包而已,但是提交到机器后,算了很久,前面一次到了walltime掉下来了,后面改了walltime不过现在一直在算,不知道怎么解决,问了一个人,他说我的输出文件的dDmax这个参数不是一直保持很大,说明是收敛的
1.请问他说得对不对,这个收敛标准对不对,如果不对怎么判断?
2.如果收敛了,请问是什么问题导致这种问题,我的MD参数都是用的默认的
回复此楼

» 猜你喜欢

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

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

love5264

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你最好把输入文件贴出来 大家看看什么问题
谋定而动
2楼2012-12-29 17:21:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxuezhang

金虫 (正式写手)

引用回帖:
2楼: Originally posted by love5264 at 2012-12-29 17:21:11
你最好把输入文件贴出来 大家看看什么问题

麻烦大侠看一看
SystemLabel     Ti3SiC2
NumberOfAtoms    12
NumberOfSpecies  3

%block ChemicalSpeciesLabel
  1    22  Ti
  2    6   C
  3    14  Si
%endblock ChemicalSpeciesLabel

PAO.BasisSize     DZP
PAO.EnergyShift   0.02 Ry
PAO.SplitNorm     0.2

%block kgrid_Monkhorst_Pack
  7   0   0    0.5
  0   7   0    0.5
  0   0   5    0.5
%endblock kgrid_Monkhorst_Pack

MeshCutoff           150. Ry

SolutionMethod        diagon

MaxSCFIterations      60
DM.MixingWeight       0.1
DM.NumberPulay         4
ElectronicTemperature  5 K

#DM.UseSaveDM          true
#UseSaveData          true

SpinPolarized         false
xc.functional         GGA
xc.authors            PBE

# Molecular dynamics and relaxations
MD.TypeOfRun            CG          # Type of dynamics:
MD.NumCGsteps           5000
MD.MaxCGDispl 0.1 Ang # Maximum atomic displacement
MD.VariableCell   .true.


WriteMullikenPop       1
SaveDeltaRho                  .true.
SaveRho                       .true.
SaveElectrostaticPotential    .true.
SaveTotalPotential            .true.
DM.UseSaveDM     true
UseSaveData          true

LatticeConstant 3.082 Ang
%block LatticeVectors
  1.000  0.000  0.000
  0.500  0.866  0.000
  0.000  0.000  5.712
%endblock LatticeVectors

AtomicCoordinatesFormat Fractional
%block AtomicCoordinatesAndAtomicSpecies
      0.00000     0.00000     0.00000  1
      0.00000     0.00000     0.50000  1
      0.33300     0.66670     0.13550  2
      0.66670     0.33330     0.63550  2
      0.66670     0.33330     0.86450  2
      0.33330     0.66670     0.36450  2
      0.00000     0.00000     0.25000  3
      0.00000     0.00000     0.75000  3
      0.33330     0.66670     0.57220  1
      0.33330     0.66670     0.92780  1
      0.66670     0.33330     0.07220  1
      0.66670     0.33330     0.42780  1
%endblock AtomicCoordinatesAndAtomicSpecies

# Output options

WriteCoorInitial        .true.
WriteCoorStep           .true.
WriteForces             .true.
WriteKpoints            .true.
WriteEigenvalues        .true.
WriteDM                 .true.
WriteKbands             .true.
WriteBands              .true.
WriteWaveFunctions      .true.
WriteMullikenPop        1            # Write Mulliken Population Analysis
WriteCoorXmol           .true.
WriteCoorCerius         .true.
WriteMDCoorXmol         .true.
WriteMDhistory          .true.
WriteMDXmol             .true.
3楼2012-12-29 23:16:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

love5264

木虫 (正式写手)

【答案】应助回帖

我也在学,一起交流吧。
1. MeshCutoff           150. Ry
设大点试试 (大小别超过250)
2. ElectronicTemperature  5 K
你是研究这个温度吗?
这个温度可能是你不收敛的原因。或者改大点试试
谋定而动
4楼2012-12-31 05:42:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangxuezhang 的主题更新
信息提示
请填处理意见