24小时热门版块排行榜    

查看: 438  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +4 a089 2026-08-07 4/200 2026-08-08 12:15 by arzu_hma
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 10:39 by 3OOjAIS77qg2
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] fileCode有新解读? +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见