| 查看: 864 | 回复: 3 | ||
xiongxiong5712新虫 (小有名气)
|
[求助]
原子受力没有达到收敛标准计算就结束了 已有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 ] |
回复此楼» 猜你喜欢
0854求调剂
已经有24人回复
-
收到复试调剂但是去不了
已经有3人回复
-
22408 312求调剂
已经有25人回复
-
304求调剂
已经有4人回复
-
化学070300 求调剂
已经有27人回复
-
294求调剂
已经有14人回复
-
0831生医工第一轮调剂失败求助
已经有18人回复
-
求调剂
已经有6人回复
-
290求调剂
已经有25人回复
-
一志愿沪9，生物学326求调剂
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 燃烧计算，先收敛后发散
已经有12人回复
- 关于回流与计算结果收敛的讨论
已经有11人回复
- ansys接触计算不收敛问题
已经有16人回复
- pwscf计算integrated charge出错及不收敛
已经有13人回复
- 怎样判定fluent计算是否收敛
已经有4人回复
- NEB计算扩散路径的力收敛精度设置问题
已经有16人回复
- Castep计算无法收敛，请高人指教
已经有17人回复
- 加溶剂后计算，不收敛，请高手指点
已经有5人回复
- CASTEP计算掺杂不收敛
已经有5人回复
- 计算时候总是不能收敛如何解决？
已经有5人回复
- 静态计算能收敛，同样的参数弛豫不能收敛，这是为什么？
已经有4人回复
- ELK计算严重不收敛～
已经有10人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- IRC计算没收敛
已经有9人回复
- DMol3计算频率不收敛
已经有4人回复
- 【求助】用Orca计算一个片簇的单点能，不收敛。
已经有7人回复
- 【求助】计算不收敛怎么办？
已经有4人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
xiongxiong5712
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 543
- 散金: 12
- 红花: 2
- 帖子: 257
- 在线: 267.7小时
- 虫号: 706467
- 注册: 2009-02-23
- 专业: 半导体材料
|
自己顶一下。麻烦各位大侠百忙中帮我看看，或者提示一下也行。我先说说我的看法：是不是因为MD.TypeOfRun 不是CG，计算不会进行CGmove的循环，所以在Begin CG move = 1的时候就结束了，是否收敛也不是看force里的Max，而是有一个地方给出了每个原子的受力： siesta: Atomic forces (eV/Ang): siesta: 1 -0.368915 -0.160954 -0.000062 siesta: 2 0.323778 -0.238986 -0.000043 siesta: 3 0.045044 0.399972 0.000016 siesta: 4 0.368938 0.160918 0.000009 siesta: 5 -0.323802 0.238974 -0.000020 siesta: 6 -0.045044 -0.399962 -0.000040 siesta: 7 0.000000 -0.000020 -0.000041 siesta: 8 -0.000020 0.000009 0.000023 siesta: 9 0.450554 -0.573922 -0.000001 siesta: 10 0.271759 0.677123 0.000003 siesta: 11 -0.722307 -0.103215 -0.000009 siesta: 12 -0.450542 0.573917 0.000002 siesta: 13 -0.271775 -0.677130 0.000005 siesta: 14 0.722304 0.103217 0.000011 siesta: ---------------------------------------- siesta: Tot -0.000027 -0.000058 -0.000148 最后一行是总的受力，是不是看这个地方收敛即可？（CG方法没有给出这一结果）。不知道我猜测的对不对，请大家指正。 |
2楼2013-06-01 15:10:48
newton3915
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4742.3
- 散金: 40
- 红花: 4
- 帖子: 418
- 在线: 391.3小时
- 虫号: 609189
- 注册: 2008-09-22
- 专业: 凝聚态物性 II ：电子结构
3楼2014-12-15 04:22:57
4楼2016-05-10 22:10:45
