24小时热门版块排行榜    

查看: 1719  |  回复: 5

daofengwalan

金虫 (小有名气)

[交流] 关于pwscf进行优化的问题 已有4人参与

请问我进行relax计算的时候,原子的坐标完全没有变化,这正常吗?
输出文件(部分结果):

     the Fermi energy is    -3.2199 ev

!    total energy              =     -77.04604446 Ry
     Harris-Foulkes estimate   =     -77.04604446 Ry
     estimated scf accuracy    <          4.5E-13 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =     -58.02346949 Ry
     hartree contribution      =      32.67217156 Ry
     xc contribution           =     -16.21517278 Ry
     ewald contribution        =     -35.45905453 Ry
     smearing contrib. (-TS)   =      -0.02051921 Ry

     convergence has been achieved in  36 iterations

     Forces acting on atoms (Ry/au):

     atom    1 type  1   force =     0.00000038    0.00000000    0.00000000
     atom    2 type  2   force =    -0.00000533    0.00000000    0.00000000
     atom    3 type  3   force =     0.00000516    0.00000000    0.00000000
     atom    4 type  3   force =    -0.00000031    0.00000000    0.00000000
     atom    5 type  3   force =     0.00000010    0.00000000    0.00000000

     Total force =     0.000007     Total SCF correction =     0.000000


     entering subroutine stress ...

          total   stress  (Ry/bohr**3)                   (kbar)     P=  -10.24
  -0.00008839   0.00000000   0.00000000        -13.00      0.00      0.00
   0.00000000  -0.00003109   0.00000000          0.00     -4.57      0.00
   0.00000000   0.00000000  -0.00008925          0.00      0.00    -13.13


     BFGS Geometry Optimization

     bfgs converged in   1 scf cycles and   0 bfgs steps
     (criteria: energy <  1.0E-04, force <  1.0E-03)

     End of BFGS Geometry Optimization

     Final energy   =     -77.0460444599 Ry
Begin final coordinates
回复此楼

» 猜你喜欢

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

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

1244464635

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这说明你的初始结构就已经很好了啊,很正常的。
2楼2014-05-21 08:20:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
relax是只考虑力对原子位置的影响。从计算中可以看到,原子受力非常小,自然原子位置也不会发生变化了。
另外,使用relax时,不需要&cell,相信输出文件的最开始也会提到discard &cell这样的语句吧~
思想重于技巧,内涵重于表象
3楼2014-05-21 13:14:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daofengwalan

金虫 (小有名气)

谢谢楼上的回复,我还有两个问题想请教。
做relax的时候,看到输入结果中也做了 Self-consistent Calculation,这是不是意味着,relax结束之后,就不需要再进行单独的SCF计算了呢?
vc-relax(variable-cell)和relax的区别是什么呢?确切的说,是在什么情况下需要进行vc-relax?
Many thanks in advance!
4楼2014-05-22 08:38:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
vc-relax的话貌似是会变胞,晶格参数什么的,一般优化要选这个吧;relax的话,是不是只是优化原子位置?!
5楼2014-05-22 21:32:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

padian

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by daofengwalan at 2014-05-22 08:38:05
谢谢楼上的回复,我还有两个问题想请教。
做relax的时候,看到输入结果中也做了 Self-consistent Calculation,这是不是意味着,relax结束之后,就不需要再进行单独的SCF计算了呢?
vc-relax(variable-cell)和rel ...

vc-relax是要优化cell的,而relax只是优化原子的位置

我觉得如果要优化结构获得晶格常数之类的,需要用vc-relax
6楼2014-06-03 00:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daofengwalan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] E0414, 我的本子有没有希望? +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
请填处理意见