版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

暖秋

管理员

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

回复此楼

» 猜你喜欢

291求调剂已经有25人回复
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂已经有6人回复
一志愿北京科技，085601总分305求调剂已经有10人回复
化学工程专硕324分，一志愿中国矿业大学求调剂已经有6人回复
279求调剂已经有3人回复
26年申博已经有4人回复
材料工程322分已经有5人回复
318求调剂，计算材料方向已经有7人回复
292分，材料与化工，申请调剂已经有20人回复
085600 295分求调剂已经有18人回复

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

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

查看全部散金贴

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

赞一下