24小时热门版块排行榜    

查看: 369  |  回复: 0

fwang2011

金虫 (小有名气)

[求助] 有熟悉cclib method的吗

系统信息:
TOTAL NUMBER OF BASIS SET SHELLS             =  937
NUMBER OF CARTESIAN GAUSSIAN BASIS FUNCTIONS = 2343
NUMBER OF ELECTRONS                          = 1008
CHARGE OF MOLECULE                           =    0
SPIN MULTIPLICITY                            =    1
NUMBER OF OCCUPIED ORBITALS (ALPHA)          =  504
NUMBER OF OCCUPIED ORBITALS (BETA )          =  504
TOTAL NUMBER OF ATOMS                        =  340

从Python调用cclib:
>>> from cclib.parser import GAMESS
>>> from cclib.method import CSPA
>>> s=GAMESS("new_b_00.log"
>>> ps=s.parse()
[GAMESS new_b_00.log INFO] Creating attribute atomcoords[]
[GAMESS new_b_00.log INFO] Creating attribute atomnos[]
[GAMESS new_b_00.log INFO] Creating attribute gbasis[]
[GAMESS new_b_00.log INFO] Creating attribute nbasis: 2343
[GAMESS new_b_00.log INFO] Creating attribute charge: 0
[GAMESS new_b_00.log INFO] Creating attribute mult: 1
[GAMESS new_b_00.log INFO] Creating attribute homos[]
[GAMESS new_b_00.log INFO] Creating attribute natom: 340
[GAMESS new_b_00.log INFO] Creating attribute aooverlaps[]
[GAMESS new_b_00.log INFO] Creating attribute scftargets[]
[GAMESS new_b_00.log INFO] Creating attribute scfvalues[]
[GAMESS new_b_00.log INFO] Creating attribute scfenergies[]
[GAMESS new_b_00.log INFO] Creating attribute moenergies[]
[GAMESS new_b_00.log INFO] Creating attribute mosyms[]
[GAMESS new_b_00.log INFO] Creating attribute nmo: 2343
[GAMESS new_b_00.log INFO] Creating attribute mocoeffs[]
[GAMESS new_b_00.log INFO] Creating attribute atombasis[]
[GAMESS new_b_00.log INFO] Creating attribute aonames[]
[GAMESS new_b_00.log INFO] Creating attribute coreelectrons[]
>>> m=CSPA(ps)
>>>>>> len(m.fragresults[0])
2343
>>> len(m.fragresults[0][0])
204
为何这里的结果不是340 (natom)呢?

谢谢

有可能帮忙解读一下这个<http://sourceforge.net/apps/medi ... php?title=C_squared>:
fragresults: a Numeric array[3] with spin, molecular orbital, and atoms as the axes (atomresults[1][23][4] gives the contribution of the 5th atomic/fragment orbital to the 24th beta molecular orbital)
这里的描述似乎比较乱,请达人解读一下。

Happy New Year everyone!

[ Last edited by fwang2011 on 2011-12-27 at 17:38 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 fwang2011 的主题更新
信息提示
请填处理意见