24小时热门版块排行榜

返回列表

pariah

铜虫 (小有名气)

应助: 6 (幼儿园)
金币: 41.4
散金: 5
红花: 4
帖子: 135
在线: 171.3小时
虫号: 1320067
注册: 2011-06-10
专业: 光学

[求助] 对BaTiO3的非线性光学系数计算求助

最近在重复tutorial里2005年的PRB的非线性系数计算的结果，按tutorial可以重复AlAs，可是自己在想重复BaTiO3时就不行了。
我的输入是
acell    3*7.14
rprim    0.5 0.0 0.0
            0.0 0.5 0.0
            0.0 0.0 0.5

#Definition of atomic positions and types
#*************************************
natom    5
ntypat    3
znucl    56 22 8
typat    1 2 3 3 3
xred       0.0  0.0  0.0
            0.5  0.5  0.5
            0.5  0.5  0.0
            0.0  0.5  0.5
            0.5  0.0  0.5
ixc       3

#Parameters of the SCF cycles
#****************************
iscf       5
nstep    100

#Plane wave basis and k-point grid
#*********************************
ecut       45
ecutsm    0.5
dilatmx    1.05
ngkpt    10 10 10
nshiftk    1
shiftk    0.5 0.5 0.5
以上是跟AlAs不一样的输入，其它都一样。可是现在会出现错误
occeig : ERROR -
  In a non-metallic case (occopt<3), for a RF calculation,
  if the eigenvalues are degenerate, the occupation numbers must also be degenerate.
  However, the following pair of states gave :
  k -state, band number  20, occ= 2.000000E+00, eigenvalue= 4.134323E+00,
  kq-state, band number  21, occ= 0.000000E+00, eigenvalue= 4.134323E+00.
  Action : change occopt, consistently, in GS and RF calculations.

请大家指点

回复此楼