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

登录注册

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

24小时热门版块排行榜

返回列表

identation

金虫 (著名写手)

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

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输出了这么多次？

很困惑，求高人解答

回复此楼

» 猜你喜欢

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

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

查看全部散金贴

1楼 2010-12-10 21:23:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★ ★ ★
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人)