24小时热门版块排行榜    

查看: 423  |  回复: 4
当前主题已经存档。

xli6

新虫 (初入文坛)

[交流] 【求助】how to calculate system energy?

Guys, when I calculate a 2-D system energy, I just use the following options in siesta fdf file, but after running, the structure has been changed, do you know why? how can I fix it. when I set the lattice vector I make it larger in one direction(y-direction) since I want to leave some vacuum space in that direction in order to avoid the interaction. Thanks!
LatticeConstant  1.0 Ang
%block LatticeVectors
  5.3340      0.0000    0.0000
  0.0000      16.1592    0.0000
  0.0000      0.0000    15.000
%endblock LatticeVectors
%block kgrid_Monkhorst_Pack
    10   0  0     0.0
     0  10  0     0.0
     0   0  1     0.0
MD.TypeOfRun            CG
MD.NumCGsteps           100
MaxSCFIterations        200
MeshCutoff                150 Ry
ElectronicTemperature   300 K
DM.NumberPulay          5
DM.MixingWeight               0.01
SolutionMethod                Diagon
DM.Tolerance            2.d-5  
MD.MaxForceTol                0.02 eV/Ang
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ice_rain(金币+3):谢谢交流讨论 2010-02-03 20:23
(1)如果只算个能量,结构怎么会变呢?
(2)如果只算个能量,请问您设置参数“MD.MaxForceTol  0.02 eV/Ang”是有什么意图呢?是不是程序有特殊之处?
(3)参数的问题,应该在说明书里挨个仔细对照,楼主说明书读了乎?
先进使役材料计算(https://aidme.nimte.ac.cn)
2楼2010-02-03 18:42:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhsh

木虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
MD.NumCGsteps           0
3楼2010-02-03 20:48:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xli6

新虫 (初入文坛)

引用回帖:
Originally posted by xli6 at 2010-02-03 12:52:50:
Guys, when I calculate a 2-D system energy, I just use the following options in siesta fdf file, but after running, the structure has been changed, do you know why? how can I fix it. when I set the ...

I set this option just because I want to relax my model and at the same time get the minimum energy.
4楼2010-02-05 02:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xli6

新虫 (初入文坛)

引用回帖:
Originally posted by lfhuang at 2010-02-03 18:42:34:
(1)如果只算个能量,结构怎么会变呢?
(2)如果只算个能量,请问您设置参数“MD.MaxForceTol  0.02 eV/Ang”是有什么意图呢?是不是程序有特殊之处?
(3)参数的问题,应该在说明书里挨个仔细对照,楼主说明 ...

I set this option just because I want to relax this model and calculate the minimum energy.
5楼2010-02-05 02:15:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xli6 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:50 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:22 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:11 by Kddn1pB6pAzX
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[博后之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:02 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 10/500 2026-08-01 17:51 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:51 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:43 by Kddn1pB6pAzX
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:40 by Kddn1pB6pAzX
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +16 低垂的野花 2026-07-31 22/1100 2026-08-01 16:39 by Cathyau
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见