24小时热门版块排行榜

>论坛更新日志 (3064)
>虫友互识 (235)
>导师招生 (197)
>文献求助 (188)
>论文投稿 (69)
>硕博家园 (56)
>休闲灌水 (55)
>公派出国 (49)
>考博 (42)
>招聘信息布告栏 (37)
>论文道贺祈福 (32)
>催化 (29)
>教师之家 (27)
>找工作 (27)
>考研 (25)
>博后之家 (23)

返回列表

暖秋

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.5
帖子: 70
在线: 36小时
虫号: 797547

[交流] 【求助】晶格优化崩溃？

各位朋友。本人驰豫原子位置的时候，坐标一直变大？偏离实验数据跟文献上的计算数据很大，麻烦大家给我看看，我的输入文件有什么不妥的地方，万分感激。。。
以下是我的输入文件内容：

&control
calculation='vc-relax',
restart_mode='from_scratch',
prefix='Ba'
tstress = .true.
tprnfor = .true.
pseudo_dir = './',
outdir='./'
nstep =  105  ,
etot_conv_thr = 1.0E-6  ,
forc_conv_thr = 1.0D-5 ,
dt = 50 ,
/
&system
                        ibrav = 0,
                     nat = 5,
                     ntyp = 3,
                  ecutwfc = 40.0 ,
                  ecutrho = 400.0 ,
                     nosym = .false. ,
                     nbnd = 80,
               occupations = 'smearing' ,
                  degauss = 0.05,
                  smearing = 'marzari-vanderbilt' ,
                     nspin = 1 ,
/

&electrons
mixing_beta = 0.7
conv_thr =  1.0d-8
/
&IONS
/
&CELL
cell_dynamics = 'bfgs'
press = 1000
/
。。。
。。。
。。。

K_POINTS automatic
  8 8 8 0 0 0

回复此楼

» 猜你喜欢

求助:我三月中下旬出站，青基依托单位怎么办？已经有8人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-11 19:23:58

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

暖秋

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.5
帖子: 70
在线: 36小时
虫号: 797547

输出

以下是输出：

Forces acting on atoms (Ry/au):

   atom 1 type  1 force =    0.00000000 0.00000000 0.00000000
   atom 2 type  2 force = -0.00076425 0.00076425 0.00000000
   atom 3 type  2 force =    0.00076425 -0.00076425 0.00000000
   atom 4 type  3 force = -0.00079442 -0.00079442 0.43102097
   atom 5 type  3 force =    0.00079442 0.00079442 -0.43102097

   Total force =    0.000000    Total SCF correction =    0.000000

   entering subroutine stress ...

      total stress  (Ry/bohr**3)                (kbar)    P=  989.12
0.00670743 0.00000026  -0.00001451       986.70    0.04    -2.13
0.00000026 0.00670743  -0.00001451       0.04 986.70    -2.13
  -0.00001451  -0.00001451 0.00675680       -2.13    -2.13 993.96

   number of scf cycles =  11
   number of bfgs steps =  10

   enthalpy old          = -666.0630713336 Ry
   enthalpy new          = -666.0638799665 Ry

   CASE: enthalpy_new < enthalpy_old

   new trust radius       =    0.0033684758 bohr
   new conv_thr          =    0.0000000010 Ry

new unit-cell volume = 379.17935 a.u.^3 ( 56.18858 Ang^3 )

CELL_PARAMETERS (alat=  3.96250000)
1.627610836  -0.000000404 0.000086123
  -0.000000404 1.627610836 0.000086123
0.813717998 0.813717998 2.300657444

赞一下

回复此楼

2楼2010-12-11 19:27:24

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

相关版块跳转我要订阅楼主暖秋的主题更新

返回列表

24小时热门版块排行榜

暖秋

[交流] 【求助】晶格优化崩溃？

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到: 查看全部散金贴

暖秋

输出

» 抢金币啦！回帖就可以得到:

查看全部散金贴