24小时热门版块排行榜    

查看: 540  |  回复: 1

[交流] 【求助】请教声子计算,为什么会出现这样的问题呢

Representation #  1 mode #   1

     Self-consistent Calculation

     Pert. #  1: Fermi energy shift (Ry) =     0.2326E+80     0.5022E+60

      iter #   1 total cpu time :  1327.1 secs   av.it.:   8.9
      thresh= 0.100E-01 alpha_mix =  0.700 |ddv_scf|^2 =  0.464E-06

     Pert. #  1: Fermi energy shift (Ry) =    -0.1429E+79     0.5649E+59

      iter #   2 total cpu time :  2269.2 secs   av.it.:  20.3
      thresh= 0.681E-04 alpha_mix =  0.700 |ddv_scf|^2 =  0.617E-05

     Pert. #  1: Fermi energy shift (Ry) =     0.1217E+79    -0.2260E+60

      iter #   3 total cpu time :  3087.7 secs   av.it.:  17.4
      thresh= 0.248E-03 alpha_mix =  0.700 |ddv_scf|^2 =  0.142E-06

这里出现的Fermi energy shift为什么那么大?

     Mode symmetry, C_4h (4/m)  point group:

     omega(  1 -  1) =         40.3  [cm-1]   -->   ?
     omega(  1 -  1) =         40.3  [cm-1]   -->   0E_g             R  
     omega(  1 -  1) =         40.3  [cm-1]   -->   ?
     omega(  1 -  1) =         40.3  [cm-1]   -->   0E_g*            R  
     omega(  1 -  1) =         40.3  [cm-1]   --> A_u             I  
     omega(  1 -  1) =         40.3  [cm-1]   -->   ?
     omega(  1 -  1) =         40.3  [cm-1]   -->   0E_u             I  
     omega(  1 -  1) =         40.3  [cm-1]   -->   ?
     omega(  1 -  1) =         40.3  [cm-1]   -->   0E_u*            I  
     omega(  2 -  3) =         50.6  [cm-1]   --> E_u             I  
     omega(  2 -  3) =         50.6  [cm-1]   --> E_u*            I  
     omega(  4 -  5) =         52.2  [cm-1]   --> E_g             R  
     omega(  4 -  5) =         52.2  [cm-1]   --> E_g*            R  
     omega(  6 -  6) =         64.7  [cm-1]   --> A_g             R  
     omega(  6 -  6) =         64.7  [cm-1]   -->   ?
     omega(  6 -  6) =         64.7  [cm-1]   -->   0E_g             R  
     omega(  6 -  6) =         64.7  [cm-1]   -->   ?
     omega(  6 -  6) =         64.7  [cm-1]   -->   0E_g*            R  
     omega(  6 -  6) =         64.7  [cm-1]   -->   ?
     omega(  6 -  6) =         64.7  [cm-1]   -->   0E_u             I  
     omega(  6 -  6) =         64.7  [cm-1]   -->   ?
     omega(  6 -  6) =         64.7  [cm-1]   -->   0E_u*            I  
     omega(  7 -  8) =        133.9  [cm-1]   --> E_u             I  
     omega(  7 -  8) =        133.9  [cm-1]   --> E_u*            I  
     omega(  9 -  9) =        134.8  [cm-1]   -->   ?
     omega(  9 -  9) =        134.8  [cm-1]   -->   0E_g             R  
     omega(  9 -  9) =        134.8  [cm-1]   -->   ?
     omega(  9 -  9) =        134.8  [cm-1]   -->   0E_g*            R  
     omega(  9 -  9) =        134.8  [cm-1]   --> A_u             I  
     omega(  9 -  9) =        134.8  [cm-1]   -->   ?
     omega(  9 -  9) =        134.8  [cm-1]   -->   0E_u             I  
     omega(  9 -  9) =        134.8  [cm-1]   -->   ?
     omega(  9 -  9) =        134.8  [cm-1]   -->   0E_u*            I  
     omega( 10 - 11) =        149.3  [cm-1]   --> E_g             R  
     omega( 10 - 11) =        149.3  [cm-1]   --> E_g*            R  
     omega( 12 - 12) =        150.3  [cm-1]   --> B_g             R  
     omega( 12 - 12) =        150.3  [cm-1]   -->   ?
     omega( 12 - 12) =        150.3  [cm-1]   -->   0E_g             R  
     omega( 12 - 12) =        150.3  [cm-1]   -->   ?
     omega( 12 - 12) =        150.3  [cm-1]   -->   0E_g*            R  
     omega( 12 - 12) =        150.3  [cm-1]   -->   ?
     omega( 12 - 12) =        150.3  [cm-1]   -->   0E_u             I  
     omega( 12 - 12) =        150.3  [cm-1]   -->   ?
     omega( 12 - 12) =        150.3  [cm-1]   -->   0E_u*            I  
     omega( 13 - 13) =        171.7  [cm-1]   -->   ?
     omega( 13 - 13) =        171.7  [cm-1]   -->   0E_g             R  
     omega( 13 - 13) =        171.7  [cm-1]   -->   ?
     omega( 13 - 13) =        171.7  [cm-1]   -->   0E_g*            R  
     omega( 13 - 13) =        171.7  [cm-1]   --> A_u             I  
     omega( 13 - 13) =        171.7  [cm-1]   -->   ?
     omega( 13 - 13) =        171.7  [cm-1]   -->   0E_u             I  
     omega( 13 - 13) =        171.7  [cm-1]   -->   ?
     omega( 13 - 13) =        171.7  [cm-1]   -->   0E_u*            I  
     omega( 14 - 15) =        178.1  [cm-1]   --> E_u             I  
     omega( 14 - 15) =        178.1  [cm-1]   --> E_u*            I  
     omega( 16 - 17) =        275.2  [cm-1]   --> E_u             I  
     omega( 16 - 17) =        275.2  [cm-1]   --> E_u*            I  
     omega( 18 - 18) =        275.8  [cm-1]   -->   ?
     omega( 18 - 18) =        275.8  [cm-1]   -->   0E_g             R  
     omega( 18 - 18) =        275.8  [cm-1]   -->   ?
     omega( 18 - 18) =        275.8  [cm-1]   -->   0E_g*            R  
     omega( 18 - 18) =        275.8  [cm-1]   --> B_u               
     omega( 18 - 18) =        275.8  [cm-1]   -->   ?
     omega( 18 - 18) =        275.8  [cm-1]   -->   0E_u             I  
     omega( 18 - 18) =        275.8  [cm-1]   -->   ?
     omega( 18 - 18) =        275.8  [cm-1]   -->   0E_u*            I  
     omega( 19 - 19) =        399.5  [cm-1]   -->   ?
     omega( 19 - 19) =        399.5  [cm-1]   -->   0E_g             R  
     omega( 19 - 19) =        399.5  [cm-1]   -->   ?
     omega( 19 - 19) =        399.5  [cm-1]   -->   0E_g*            R  
     omega( 19 - 19) =        399.5  [cm-1]   --> A_u             I  
     omega( 19 - 19) =        399.5  [cm-1]   -->   ?
     omega( 19 - 19) =        399.5  [cm-1]   -->   0E_u             I  
     omega( 19 - 19) =        399.5  [cm-1]   -->   ?
     omega( 19 - 19) =        399.5  [cm-1]   -->   0E_u*            I  
     omega( 20 - 21) =        400.6  [cm-1]   --> E_u             I  
     omega( 20 - 21) =        400.6  [cm-1]   --> E_u*            I  
输出的结果中为什么有这么多问好,然后同一个omega输出了这么多次?

很困惑,求高人解答
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★
zhang668(金币+2):多谢交流 2010-12-12 23:02:36
mazuju028(金币+2):谢谢交流 2010-12-12 23:07:40
identation(金币+10): 2011-04-25 21:31:57
这个应该和你的体系对称性有关系。我没碰到过这个问题,但是你应该看看你的结构是不是已经最优化了,或者里面有什么新的物理,超导什么的。或者你的体系根本有什么现象没有在这里的DFT里面考虑
2楼2010-12-12 22:53:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见