24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2116  |  回复: 10

wanfei

金虫 (小有名气)

[交流] 【求助】分析HOMO LUMO已有6人参与

我初学MS希望大家帮我看下下面的HOMO LUMO值怎么取,谢谢
Fermi Energy  -0.18133Ha    -4.934eV   -kTS_e= -0.054975Ha
  HOMO is orbital number   478
  LUMO is orbital number   479

    state                         eigenvalue        occupation
                              (au)          (ev)

    1   -   1   a         -9.969555      -271.286     1.000
    2   +   1   a         -9.969555      -271.286     1.000
    3   -   2   a         -9.969525      -271.285     1.000
    4   +   2   a         -9.969525      -271.285     1.000
    5   +   3   a         -9.969184      -271.275     1.000
    6   -   3   a         -9.969184      -271.275     1.000
    7   +   4   a         -9.968916      -271.268     1.000
    8   -   4   a         -9.968916      -271.268     1.000
    9   -   5   a         -9.968785      -271.265     1.000
   10   +   5   a         -9.968785      -271.265     1.000
   11   -   6   a         -9.968658      -271.261     1.000
   12   +   6   a         -9.968658      -271.261     1.000
   13   +   7   a         -9.968179      -271.248     1.000
   14   -   7   a         -9.968179      -271.248     1.000
   29   +  15   a         -9.964444      -271.146     1.000
  
   
  407   - 204   a         -0.331496        -9.020     1.000
  408   + 204   a         -0.331496        -9.020     1.000
  409   - 205   a         -0.325938        -8.869     1.000
  410   + 205   a         -0.325938        -8.869     1.000
  411   + 206   a         -0.316792        -8.620     1.000
  412   - 206   a         -0.316792        -8.620     1.000
  413   + 207   a         -0.316144        -8.603     1.000
  414   - 207   a         -0.316144        -8.603     1.000
  415   - 208   a         -0.311408        -8.474     1.000
  416   + 208   a         -0.311408        -8.474     1.000
  417   - 209   a         -0.303482        -8.258     1.000
  418   + 209   a         -0.303482        -8.258     1.000
  419   - 210   a         -0.298791        -8.131     1.000
  420   + 210   a         -0.298791        -8.131     1.000
  421   + 211   a         -0.296928        -8.080     1.000
  422   - 211   a         -0.296928        -8.080     1.000
  423   - 212   a         -0.293587        -7.989     1.000
  424   + 212   a         -0.293587        -7.989     1.000
  425   + 213   a         -0.287264        -7.817     1.000
  426   - 213   a         -0.287264        -7.817     1.000
  427   + 214   a         -0.287097        -7.812     1.000
  428   - 214   a         -0.287097        -7.812     1.000
  429   + 215   a         -0.276340        -7.520     1.000
  430   - 215   a         -0.276340        -7.520     1.000
  431   + 216   a         -0.275441        -7.495     1.000
  432   - 216   a         -0.275441        -7.495     1.000
  433   + 217   a         -0.274374        -7.466     1.000
  434   - 217   a         -0.274374        -7.466     1.000
  435   + 218   a         -0.270211        -7.353     1.000
  436   - 218   a         -0.270211        -7.353     1.000
  437   - 219   a         -0.269856        -7.343     1.000
  438   + 219   a         -0.269856        -7.343     1.000
  439   + 220   a         -0.268701        -7.312     1.000
  440   - 220   a         -0.268701        -7.312     1.000
  441   + 221   a         -0.267891        -7.290     1.000
  442   - 221   a         -0.267891        -7.290     1.000
  443   + 222   a         -0.266600        -7.255     1.000
  444   - 222   a         -0.266600        -7.255     1.000
  445   - 223   a         -0.265851        -7.234     1.000
  446   + 223   a         -0.265851        -7.234     1.000
  447   - 224   a         -0.265107        -7.214     1.000
  448   + 224   a         -0.265107        -7.214     1.000
  449   + 225   a         -0.262779        -7.151     1.000
  450   - 225   a         -0.262779        -7.151     1.000
  451   + 226   a         -0.261559        -7.117     1.000
  452   - 226   a         -0.261559        -7.117     1.000
  453   + 227   a         -0.260397        -7.086     1.000
  454   - 227   a         -0.260397        -7.086     1.000
  455   - 228   a         -0.257449        -7.006     1.000
  456   + 228   a         -0.257449        -7.006     1.000
  457   + 229   a         -0.256691        -6.985     0.999
  458   - 229   a         -0.256691        -6.985     0.999
  459   + 230   a         -0.254453        -6.924     0.999
  460   - 230   a         -0.254453        -6.924     0.999
  461   - 231   a         -0.251578        -6.846     0.999
  462   + 231   a         -0.251578        -6.846     0.999
  463   - 232   a         -0.250646        -6.820     0.999
  464   + 232   a         -0.250646        -6.820     0.999
  465   + 233   a         -0.242699        -6.604     0.998
  466   - 233   a         -0.242699        -6.604     0.998
  467   - 234   a         -0.238706        -6.496     0.997
  468   + 234   a         -0.238706        -6.496     0.997
  469   - 235   a         -0.232908        -6.338     0.994
  470   + 235   a         -0.232908        -6.338     0.994
  471   + 236   a         -0.230564        -6.274     0.993
  472   - 236   a         -0.230564        -6.274     0.993
  473   + 237   a         -0.226456        -6.162     0.989
  474   - 237   a         -0.226456        -6.162     0.989
  475   - 238   a         -0.211887        -5.766     0.955
  476   + 238   a         -0.211887        -5.766     0.955
  477   - 239   a         -0.203470        -5.537     0.902
  478   + 239   a         -0.203470        -5.537     0.902
  479   - 240   a         -0.177873        -4.840     0.415
  480   + 240   a         -0.177873        -4.840     0.415
  481   + 241   a         -0.154393        -4.201     0.063
  482   - 241   a         -0.154393        -4.201     0.063
  483   + 242   a         -0.153586        -4.179     0.059
  484   - 242   a         -0.153586        -4.179     0.059
  485   - 243   a         -0.153538        -4.178     0.058
  486   + 243   a         -0.153538        -4.178     0.058
  487   + 244   a         -0.150324        -4.091     0.043
  488   - 244   a         -0.150324        -4.091     0.043
  489   - 245   a         -0.147892        -4.024     0.034
  490   + 245   a         -0.147892        -4.024     0.034
  491   - 246   a         -0.141360        -3.847     0.018
  492   + 246   a         -0.141360        -3.847     0.018
  493   - 247   a         -0.135381        -3.684     0.010
  494   + 247   a         -0.135381        -3.684     0.010
  495   - 248   a         -0.135092        -3.676     0.010
  496   + 248   a         -0.135092        -3.676     0.010
  497   + 249   a         -0.130042        -3.539     0.006
  498   - 249   a         -0.130042        -3.539     0.006
  499   + 250   a         -0.129341        -3.520     0.005
  500   - 250   a         -0.129341        -3.520     0.005
  501   + 251   a         -0.121936        -3.318     0.003
  502   - 251   a         -0.121936        -3.318     0.003
  503   - 252   a         -0.117864        -3.207     0.002
  504   + 252   a         -0.117864        -3.207     0.002
  505   - 253   a         -0.109848        -2.989     0.001
  506   + 253   a         -0.109848        -2.989     0.001
  507   + 254   a         -0.106967        -2.911     0.001
  508   - 254   a         -0.106967        -2.911     0.001
  509   + 255   a         -0.103894        -2.827     0.000
  510   - 255   a         -0.103894        -2.827     0.000
  511   + 256   a         -0.098862        -2.690     0.000
  512   - 256   a         -0.098862        -2.690     0.000
  513   - 257   a         -0.098511        -2.681     0.000
  514   + 257   a         -0.098511        -2.681     0.000
  515   + 258   a         -0.093410        -2.542     0.000
  516   - 258   a         -0.093410        -2.542     0.000
  517   - 259   a         -0.092371        -2.514     0.000
  518   + 259   a         -0.092371        -2.514     0.000
  519   - 260   a         -0.083665        -2.277     0.000
  520   + 260   a         -0.083665        -2.277     0.000
  521   - 261   a         -0.081790        -2.226     0.000
  522   + 261   a         -0.081790        -2.226     0.000
  523   - 262   a         -0.079663        -2.168     0.000
  524   + 262   a         -0.079663        -2.168     0.000
  525   + 263   a         -0.078240        -2.129     0.000
  526   - 263   a         -0.078240        -2.129     0.000
  527   - 264   a         -0.071906        -1.957     0.000
  528   + 264   a         -0.071906        -1.957     0.000
  529   + 265   a         -0.067904        -1.848     0.000
  530   - 265   a         -0.067904        -1.848     0.000
  531   - 266   a         -0.060448        -1.645     0.000
  532   + 266   a         -0.060448        -1.645     0.000
  533   + 267   a         -0.042190        -1.148     0.000
  534   - 267   a         -0.042190        -1.148     0.000
  535   - 268   a         -0.039995        -1.088     0.000
  536   + 268   a         -0.039995        -1.088     0.000
  537   + 269   a         -0.031573        -0.859     0.000
  538   - 269   a         -0.031573        -0.859     0.000
  539   - 270   a         -0.029615        -0.806     0.000
  540   + 270   a         -0.029615        -0.806     0.000
  541   - 271   a         -0.020366        -0.554     0.000
  542   + 271   a         -0.020366        -0.554     0.000
  543   + 272   a         -0.018596        -0.506     0.000
  544   - 272   a         -0.018596        -0.506     0.000
  545   - 273   a          0.002193         0.060     0.000
  546   + 273   a          0.002193         0.060     0.000
  547   + 274   a          0.005364         0.146     0.000
  550   + 275   a          0.006215         0.169     0.000
  551   + 276   a          0.010034         0.273     0.000
  553   + 277   a          0.013892         0.378     0.000
  556   + 278   a          0.030601         0.833     0.000
  557   + 279   a          0.033528         0.912     0.000
  560   + 280   a          0.038694         1.053     0.000
  562   + 281   a          0.039440         1.073     0.000
  563   + 282   a          0.052195         1.420     0.000

  Orbital occupation is:
  232  a  ( 1,) alpha and 232  a  ( 1) beta, and
     7.56 alpha elec in  27 orbitals plus    7.56  beta elec in  27 orbitals
  Total number electrons:             479.1


        Total E (au)     Binding E (au)     Convergence     Time (m)   Iter
-----------------------------------------------------------------------------
Ef     -3034.2278824      -21.3353035         0.66E-06      878.033     14

binding energy -21.335304Ha -580.56339eV -13388.372kcal/mol  -56016.950kJ/mol
-----------------------------------------------------------------------------

~~~~~~~~~~~~~~~~~~~~~ End Computing SCF Energy/Gradient ~~

因为太长,中间occupation为1.00的扇区了部分。另外能量值是取Ef     -3034.2278824吗?



我算是纳米管,想知道CNT(5,5) CNT(5,0) CNT(8,0)这些有什么区别?谢谢大家了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
zzy870720z(金币+2):感谢专家指导 2010-05-31 18:41:32
开头的部分已经很清楚的告诉你了:
Fermi Energy  -0.18133Ha    -4.934eV   -kTS_e= -0.054975Ha
  HOMO is orbital number   478
  LUMO is orbital number   479

因为你采用smearing的方法来确定费米能级的,所以第479个轨道的占有数为小于1接近0的小数。

CNT具有一维周期,看Fermi level,画出能带图来,找出导带底和价带顶这样更有意义。

[ Last edited by valenhou001 on 2010-5-31 at 17:05 ]
2楼2010-05-31 17:02:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanfei

金虫 (小有名气)

478   + 239   a         -0.203470        -5.537     0.902
  479   - 240   a         -0.177873        -4.840     0.415
那我的这个能隙就是-5.537-(-4.840)=0.697ev?
3楼2010-05-31 17:12:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)


cenwanglai(金币+1):谢谢参与 2010-06-03 16:30:06
对   不过感觉这种不太准  不过要是算非周期性的话 也只好用这个了
周期性的话还是能带来的好些
4楼2010-05-31 20:38:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanfei

金虫 (小有名气)

我算非金属c纳米管,能带计算选择不了啊,没办法算。能告诉我纳米管计算怎么让HOMO LUMO值更准确吗,另外怎么设置就是周期性还是非周期性啊,刚刚学,请大家指点下,谢谢
5楼2010-06-02 10:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ym23

铁杆木虫 (著名写手)


zzy870720z(金币+1):谢谢交流 2010-06-03 08:18:39
引用回帖:
Originally posted by wanfei at 2010-06-02 10:48:50:
我算非金属c纳米管,能带计算选择不了啊,没办法算。能告诉我纳米管计算怎么让HOMO LUMO值更准确吗,另外怎么设置就是周期性还是非周期性啊,刚刚学,请大家指点下,谢谢

是不是周期性依据的是你的模型,而不是计算参数。模型可以分成周期性模型和非周期性模型。因为你选择的是非周期性模型,所以不能计算能带。

» 本帖已获得的红花(最新10朵)

博观约取、厚积薄发---壁立千仞、无欲则刚
6楼2010-06-02 22:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanfei

金虫 (小有名气)

谢谢各位了
7楼2010-06-03 09:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanfei

金虫 (小有名气)

谢谢各位了
8楼2010-06-03 10:02:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolianshun

金虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by wanfei at 2010-05-31 17:12:16:
478   + 239   a         -0.203470        -5.537     0.902
  479   - 240   a         -0.177873        -4.840     0.415
那我的这个能隙就是-5.537-(-4.840)=0.697ev?

学习路过
9楼2010-07-16 08:59:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolianshun

金虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖
引用回帖:
Originally posted by wanfei at 2010-06-02 10:48:50:
我算非金属c纳米管,能带计算选择不了啊,没办法算。能告诉我纳米管计算怎么让HOMO LUMO值更准确吗,另外怎么设置就是周期性还是非周期性啊,刚刚学,请大家指点下,谢谢

HOMO,LUMO,的轨道图怎么做啊
10楼2011-07-17 18:00:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wanfei 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见