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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » siesta计算问题
41/1返回列表
查看: 678  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yangxuezhang

金虫 (正式写手)

[求助] siesta计算问题

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-12-28 14:20:30
已阅   回复此楼   关注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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭