24小时热门版块排行榜    

查看: 449  |  回复: 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文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 17:34 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 16:34 by xiangy672
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 4/200 2026-08-30 16:03 by sundiv
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 13:54 by GkCVyODlx4mR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:02 by l0VvVHGBGRLv
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见