应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » molpro软件计算的组态信息解读
【调剂】北京石油化工学院2024年16个专业接受调剂
61/1返回列表
查看: 2735  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者18766315007将赠送您 50 个金币

18766315007

新虫 (著名写手)

[求助] molpro软件计算的组态信息解读已有1人参与

各位专家学者，如何理解CASSCF和CI计算的组态信息，里面的每个符号代表什么含义，如何和通常说的(1σg)2(1σu)2(2σg)2(2σu)(3σg)2(1πu)4(1πg)这种表达对应起来？

下图一个是CASSCF的结果，一个是CI 的结果。
希望大家帮忙，小弟不胜感激！

molpro软件计算的组态信息解读
CASSCF.png


molpro软件计算的组态信息解读-1
组态.png
回复此楼

» 猜你喜欢
1楼2018-04-06 10:45:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
18766315007: 金币+50, ★★★很有帮助 2018-04-22 11:51:27
要读懂这些需要跟MOs output结合起来阅读，否则不知道第几号是什么character的轨道。
加
print orbital
gprint,orbitals=5 （这个是多打印几个未占据轨道，一般选4，5就够了）
casscf还要加 natorb,****.2 把natural orbital打印出来。

然后你进output查看，在对应的casscf计算下有轨道打印出来了。然后你可以一个个查看轨道的character然后判断是哪个mo，
你会发现molpro默认的mo很难度，这里推荐一个script： https://github.com/RulingF/molpro_read_MOs
可以帮你大大减轻读mo的难度。这个script打印出来的大致是，以N2做例子:
Orb        Occ        Energy        Couls-En        Coefficients       
1.1        2        -15.6910 -43.5330           || 1 1s  :  1.00
2.1        2        -1.4789 -11.4890           || 1 2s  :  0.89  || 1 2pz :  0.49  || 1 4pz : -0.18  || 1 4s  : -0.20
3.1        2        -0.6352 -8.4685           || 1 2pz :  0.99  || 1 4pz : -0.15  || 1 4s  : -0.21  || 1 2s  : -0.34

1.2        2        -0.6175 -8.6822           || 1 2px :  0.89

1.3        2        -0.6175 -8.6822           || 1 2py :  0.89

1.5        2        -15.6874 -43.5303           || 1 1s  :  1.00
2.5        2        -0.7781 -8.4476           || 1 2pz :  0.46  || 1 4pz : -0.16  || 1 4s  : -0.31  || 1 2s  : -0.93
可自行读出MO类型。

好了以上是一个hf 轨道的例子，对于casscf来说，natural orbital会给一个轨道占据均值，所以你看到的占据会有非整数形式，对应的就是一个个ci vector所标识的行列式的占据均值。至于怎么跟传统mo联系，我只能告诉你，除了内层轨道之外，很多时候这些传统的概念已经没用了，这里就是很明白的第几个态由几个行列式由什么方式构成，至于你能从中获得什么信息你靠你自己的需求。
赞一下(1人) 回复此楼
2楼2018-04-19 05:20:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Lily0603

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 121099 at 2018-04-19 05:20:15
要读懂这些需要跟MOs output结合起来阅读，否则不知道第几号是什么character的轨道。
加
print orbital
gprint,orbitals=5 （这个是多打印几个未占据轨道，一般选4，5就够了）
casscf还要加 natorb,****.2 把na ...

您好，我看了这个script，确实很方便。不过我没用过python，我只能输出到桌面上，请问怎么输出到文件中?

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
3楼2018-05-06 16:57:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)
引用回帖:
3楼: Originally posted by Lily0603 at 2018-05-06 16:57:35
您好，我看了这个script，确实很方便。不过我没用过python，我只能输出到桌面上，请问怎么输出到文件中?
...

是需要打印到屏幕上么？
还是另外建一个文件
可以 molpro_read > userdefinedfilename.txt
赞一下(1人) 回复此楼
4楼2018-05-10 07:57:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo854639263
h

发自小木虫Android客户端
5楼2018-07-10 04:04:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cwllxy

新虫 (小有名气)
引用回帖:
2楼: Originally posted by 121099 at 2018-04-19 05:20:15
要读懂这些需要跟MOs output结合起来阅读，否则不知道第几号是什么character的轨道。
加
print orbital
gprint,orbitals=5 （这个是多打印几个未占据轨道，一般选4，5就够了）
casscf还要加 natorb,****.2 把na ...

您好，我使用了这个script，出现SyntaxError: Missing parentheses in call to 'print'. Did you mean print("%s\t"%item, end=" "?这个报错，应该怎么办呢
赞一下 回复此楼
6楼2022-08-16 15:53:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 18766315007 的主题更新
61/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 科研好难啊，想退学去打工。 +6 byywnyl 2024-04-27 7/350 2024-04-28 01:29 by 鱼翔浅底1
[论文投稿] AC审稿 +3 马winner 2024-04-26 7/350 2024-04-27 23:33 by topedit
[硕博家园] 聊天 +13 暮色恋伊人 2024-04-22 14/700 2024-04-27 23:26 by @tanzelin
[硕博家园] 失之东隅 收之桑榆 +5 半世迷离 2024-04-27 10/500 2024-04-27 22:19 by 半世迷离
[找工作] 普通院校药学硕士，做合成的，感觉找不到工作 +12 pom戴墨镜 2024-04-24 20/1000 2024-04-27 22:08 by saenpe
[教师之家] 博士论文被抄袭 +27 和尚敲小木鱼 2024-04-22 44/2200 2024-04-27 21:57 by ay_guobo
[硕博家园] 博士白读了 +49 Da_Meng_Zi 2024-04-21 54/2700 2024-04-27 20:12 by 慕容家族
[基金申请] 优青已经开始了么？ +3 majunge000 2024-04-27 3/150 2024-04-27 18:51 by spring_1987
[找工作] 江苏理工学院要慎选 +3 jjchenshui 2024-04-27 3/150 2024-04-27 17:45 by 阿司匹林888
[考博] 材料方向24博士申请/一作SCI三篇 +3 白天不碰 2024-04-24 5/250 2024-04-27 17:41 by 安塔瓦拉多
[考博] 真的好想读博！ +16 wangzhe_bs 2024-04-22 23/1150 2024-04-27 17:10 by @tanzelin
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +13 彩虹初见 2024-04-24 13/650 2024-04-27 15:24 by cuan118
[教师之家] 大学直属学院卸任的副院长退休后还享受副处级的养老待遇吗？ +5 苏东坡二世 2024-04-27 5/250 2024-04-27 14:54 by zpv2011
[考博] 25年博士申请 +7 Changzixuan 2024-04-25 13/650 2024-04-27 14:11 by 我属驴核动力驴
[考研] 学硕专硕 +5 小蜗牛* 2024-04-26 5/250 2024-04-26 16:43 by 鱼翔浅底1
[考研] 381求调剂 +4 小刺猬987654321 2024-04-25 6/300 2024-04-26 10:57 by czl12138
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[教师之家] 大家访学都是怎么找的啊？ +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
信息提示
关闭
请填处理意见
关闭