24小时热门版块排行榜

返回列表

forum_user

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 26017.5
红花: 5
帖子: 2967
在线: 65小时
虫号: 1120335
注册: 2010-10-12
专业: 环境化工

[求助] 计算动力学矩阵遇到的收敛问题

各位虫友：

我在用Ph.x计算如下体系的动力学矩阵时，遇到了如下问题:

   Electric Fields Calculation
   ik 1 ibnd  83 linter: root not converged  2.076E-09
   ik 1 ibnd  83 linter: root not converged  9.963E-10
   ik 1 ibnd  83 linter: root not converged  3.695E-09
   ik 2 ibnd  83 linter: root not converged  1.564E-09
   ik 2 ibnd  83 linter: root not converged  1.586E-09
   ik 2 ibnd  83 linter: root not converged  1.386E-09
   ik 3 ibnd  83 linter: root not converged  2.491E-09
   ik 3 ibnd  83 linter: root not converged  2.811E-09
   ik 3 ibnd  83 linter: root not converged  1.902E-09
   ik 4 ibnd  83 linter: root not converged  1.426E-09
   ik 4 ibnd  83 linter: root not converged  1.709E-09
   ik 4 ibnd  83 linter: root not converged  1.195E-09
   ik 5 ibnd  83 linter: root not converged  2.458E-09

最后并没有输出所要的动力学矩阵。

计算的输入文件如下:

title_line
&inputph
  tr2_ph=1.0d-12,
  prefix='mgalo',
  fildvscf='mgfeodv',
  amass(1)=24.30,
  amass(2)=27.0,
  amass(3)=16.00,
  outdir='./tmp',
  fildyn='mgfeo.dyn',
  ldisp=.true.
  nq1=2, nq2=2, nq3=2
/

请各位指点迷津，奖励金币10个，谢谢！

回复此楼