24小时热门版块排行榜    

查看: 447  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 12:13 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:02 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:30 by l0VvVHGBGRLv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:10 by l0VvVHGBGRLv
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 17/850 2026-08-30 10:28 by 孤独的英雄6
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见