| 查看: 637 | 回复: 1 | |||
[交流]
【求助】晶格优化崩溃？
|
|
各位朋友。本人驰豫原子位置的时候，坐标一直变大？偏离实验数据跟文献上的计算数据很大，麻烦大家给我看看，我的输入文件有什么不妥的地方，万分感激。。。 以下是我的输入文件内容： &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人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有6人回复
-
我面上完蛋了
已经有13人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于优化晶格常数
已经有4人回复
- 晶格常数优化
已经有7人回复
- 【求助】请高人指点一下掺氮TiO2的制备
已经有19人回复
- 【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】简单问题----关于vasp晶格优化
已经有7人回复
- 【求助成功】晶格优化的结果很诡异
已经有13人回复
- 【求助成功】晶格优化中，murn.f的使用方法
已经有7人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【原创】晶格优化中ISIF取值的探究（六方MgO为例）
已经有9人回复
输出
|
以下是输出： 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
20