24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiongxiong5712

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 543
散金: 12
红花: 2
帖子: 257
在线: 267.7小时
虫号: 706467
注册: 2009-02-23
专业: 半导体材料

[求助] 原子受力没有达到收敛标准计算就结束了已有1人参与

用siesta进行结构优化，由于是六方晶系，首先通过拟合得到了能量最低对应的体积，MD.TypeOfRun设置为CG，然后再保持这个体积不变，优化晶格参数，MD.TypeOfRun 设置为Broyden，具体的fdf文件如下：
SystemName       bulk
SystemLabel       bulk

NumberOfAtoms    14
NumberOfSpecies    2

%block ChemicalSpeciesLabel
1 14  Si
2 7 N
%endblock ChemicalSpeciesLabel

LatticeConstant    1 Ang

%block LatticeVectors
   6.571864087507342    -3.794267500000001    0.000000000000000
   0.000000000000000    7.588534999999999    0.000000000000000
   0.000000000000000    0.000000000000000    2.906249000000000
%endblock LatticeVectors

%block kgrid_Monkhorst_Pack
5 0 0 0
0 5 0 0
0 0 12 0
%endblock kgrid_Monkhorst_Pack

AtomicCoordinatesFormat    Fractional
%block AtomicCoordinatesAndAtomicSpecies
0.3323000000000000 0.0314000000000000 0.2500000000000000 2
-0.0314000000000000 0.3009000000000000 0.2500000000000000 2
-0.3009000000000000 -0.3323000000000000 0.2500000000000000 2
-0.3323000000000000 -0.0314000000000000 0.7500000000000000 2
0.0314000000000000 -0.3009000000000000 0.7500000000000000 2
0.3009000000000000 0.3323000000000000 0.7500000000000000 2
0.3333333333333330 0.6666666666666670 0.2500000000000000 2
-0.3333333333333330 -0.6666666666666670 0.7500000000000000 2
0.1733000000000000 -0.2306000000000000 0.2500000000000000 1
0.2306000000000000 0.4039000000000000 0.2500000000000000 1
-0.4039000000000000 -0.1733000000000000 0.2500000000000000 1
-0.1733000000000000 0.2306000000000000 0.7500000000000000 1
-0.2306000000000000 -0.4039000000000000 0.7500000000000000 1
0.4039000000000000 0.1733000000000000 0.7500000000000000 1
%endblock AtomicCoordinatesAndAtomicSpecies

xc.functional       LDA
xc.authors          CA
MeshCutoff          100 Ry

PAO.BasisSize       DZP
PAO.EnergyShift    75 mev

MD.MaxForceTol    0.05 eV/Ang

DM.MixingWeight    0.01
DM.NumberPulay    10
MaxSCFIterations    200

SolutionMethod    diagon

MD.TypeOfRun       Broyden
MD.NumCGsteps       500

MD.VariableCell       .true.
MD.ConstantVolume    .true.
MD.RelaxCellOnly    .true.

优化很快， Begin CG move =  1的时候就结束了，dMax达到了收敛标准，但是原子受力没有收敛，结果是：
siesta: Atomic forces (eV/Ang):
----------------------------------------
Tot -0.000014 -0.000041 -0.000144
----------------------------------------
Max 0.722244
Res 0.314921 sqrt( Sum f_i^2 / 3N )
----------------------------------------
Max 0.722244 constrained

可见原子受力没有达到0.05以下，而计算就正常结束了。请教大家，这是什么原因引起的啊？是不是我MD.TypeOfRun 设置为Broyden的原因呢？用这种方法优化的时候应该注意什么？有劳各位帮帮忙，多谢了。

[ Last edited by xiongxiong5712 on 2013-5-31 at 21:46 ]

回复此楼

» 猜你喜欢

自荐读博已经有9人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

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

1楼 2013-05-31 16:23:13

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

newton3915

铁杆木虫 (正式写手)

应助: 4 (幼儿园)
金币: 4742.3
散金: 40
红花: 4
帖子: 418
在线: 391.3小时
虫号: 609189
注册: 2008-09-22
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by xiongxiong5712 at 2013-06-01 15:10:48
自己顶一下。麻烦各位大侠百忙中帮我看看，或者提示一下也行。我先说说我的看法：是不是因为MD.TypeOfRun 不是CG，计算不会进行CGmove的循环，所以在Begin CG move = 1的时候就结束了，是否收敛也不是看force里的M ...

LZ问题解决了吗，我也是遇到同样的问题？

赞一下