24小时热门版块排行榜

返回列表

暖秋

新虫 (小有名气)

应助: 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

回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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小时热门版块排行榜

暖秋

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

» 猜你喜欢

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

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

暖秋

输出

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

查看全部散金贴