24小时热门版块排行榜     【意得辑11月】润色翻译满减最高600元

查看: 350  |  回复: 7
【悬赏金币】回答本帖问题,作者moulh将赠送您 6 个金币

moulh

新虫 (小有名气)

[求助] 用molrp算RCCSD(T)单点能

各位大神们,想请教一个问题:你们用molpro算RCCSD(T)单点能时是怎么个流程呢?手册上都是hf+rccsd(t),这对小分子可以,含过渡金属的就不适用了吧,一些文献中对含过渡金属体系算得是hf+casscf+rccsd(t),问题是ccsd不能读取cas的轨道,不知道他是怎么操作的,希望大神指点一下,说说一般都怎么算RCCSD(T),不胜感激!
ps:之前有前辈说cas后面加一个hf,是个很好的方法!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moulh

新虫 (小有名气)

顶一下免得沉…

发自小木虫IOS客户端
2楼2018-07-11 18:04:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)

请贴一下文献
3楼2018-07-12 03:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moulh

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 121099 at 2018-07-12 03:03:47
请贴一下文献

Chem. Eur.J .2015, 21,17748–17756, DOI :10.1002/chem.201503060
计算描述请见附图,只是文章里没说怎么用casscf波函算的RCCSD(T)
用molrp算RCCSD(T)单点能
1.PNG

4楼2018-07-12 23:27:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)

引用回帖:
4楼: Originally posted by moulh at 2018-07-12 23:27:24
Chem. Eur.J .2015, 21,17748–17756, DOI :10.1002/chem.201503060
计算描述请见附图,只是文章里没说怎么用casscf波函算的RCCSD(T)

1.PNG
...

他们可能用casscf+1 iteration HF 来获得ccsd需要的单参考轨道。
如果直接casscf+ccsd的话轨道不是nicely defined最后的绝对能量可能会较高一点。

单从这篇文章上看,没有提到任何为何要用casscf来得到reference的理由,我能想到的只有降低T1 diagnostic,或者依赖molpro multi的良好收敛性。

你如果想要完全复制这篇文章的结果 可能需要问问作者 他们具体到底是怎么做的。

如果没有这个想法,完全可以用 hf+ccsd另行计算。

» 本帖已获得的红花(最新10朵)

5楼2018-07-14 03:42:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moulh

新虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by 121099 at 2018-07-14 03:42:09
他们可能用casscf+1 iteration HF 来获得ccsd需要的单参考轨道。
如果直接casscf+ccsd的话轨道不是nicely defined最后的绝对能量可能会较高一点。

单从这篇文章上看,没有提到任何为何要用casscf来得到referen ...

casscf+1 iteration 什么意思呢?
6楼2018-07-14 07:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)

引用回帖:
6楼: Originally posted by moulh at 2018-07-14 07:46:58
casscf+1 iteration 什么意思呢?...

操作上就是casscf之后加一个 maxit,1的hf,记得读取casscf的natural orbital

» 本帖已获得的红花(最新10朵)

7楼2018-07-17 04:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moulh

新虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by 121099 at 2018-07-17 04:30:09
操作上就是casscf之后加一个 maxit,1的hf,记得读取casscf的natural orbital...

谢谢!读取cas的自然轨道就是您之前说的打印的CI vector的方法吧?

发自小木虫IOS客户端
8楼2018-07-17 14:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 moulh 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 有人一起学雅思吗 +14 mingnanqi 2018-11-15 19/950 2018-11-15 11:37 by 柠檬小咔咔
[环境] 重金属浓度检测中反应后马上测和隔天测有区别吗 +4 zhaoningzn 2018-11-13 4/200 2018-11-15 09:25 by wilddonkey
[考博] 万万没想到被快递员坑了 +9 ?快乐人生 2018-11-14 20/1000 2018-11-15 08:35 by ?快乐人生
[教师之家] 教育部发文在高校开展“唯论文、唯帽子、唯职称、唯学历、唯奖项”(“五唯”)清理 +23 傲我华科 2018-11-13 26/1300 2018-11-15 05:56 by Jessy_CN
[找工作] 高校三选一,恳求各位的意见与想法 10+3 wuxl515 2018-11-12 17/850 2018-11-14 23:33 by einstain666
[硕博家园] 在职研究生论文已经两次外审不通过 +6 i6lan 2018-11-14 14/700 2018-11-14 19:54 by robertfairy
[公派出国] 法国中央科学院(CNRS)-LHFA实验室,Nicolas MéZAILLES 课题组长期招收CSC公派留学生 (金币+5) +3 franklfls 2018-11-13 5/250 2018-11-14 18:24 by franklfls
[晶体] 单晶,小分子,晶体结构是核磁解的结构的一半,如何处理 2+6 冰山网子1.0 2018-11-14 14/700 2018-11-14 17:33 by 冰山网子1.0
[论文投稿] 杂志submit系统只能指定一个通讯作者,而我这篇需要写两个 +8 love_korla 2018-11-09 9/450 2018-11-13 14:54 by baichen0310
[考研] 考研 +6 Seniray 2018-11-12 6/300 2018-11-13 14:19 by YL170
[硕博家园] 三本给排水 硕士双非,如果有机会读同济大学博士,如何抉择 +3 l1214l 2018-11-12 5/250 2018-11-13 13:43 by l1214l
[论文投稿] Required reviews completed,又15天了。 +4 孤独前行007 2018-11-09 13/650 2018-11-13 10:08 by 孤独前行007
[论文投稿] 可以给编辑发邮件并要求重新考虑文章吗?(目前是reject状态) +5 kaen0221 2018-11-13 7/350 2018-11-13 10:00 by Group尼克
[找工作] 药企咨询 +5 化工小女卒 2018-11-10 9/450 2018-11-12 21:57 by 13961985205
[海外博后] 洪堡基金申请时间求助 5+3 爱有机 2018-11-09 5/250 2018-11-12 17:49 by meichao1989
[考博] 申博研究计划 +9 碧涧流泉 2018-11-11 13/650 2018-11-12 10:16 by zcwaxc
[考研] 考研 +3 1303285013 2018-11-08 3/150 2018-11-12 00:34 by 东城旧事
[硕博家园] 期刊多久可以online +7 qinxueyang 2018-11-10 8/400 2018-11-12 00:14 by mengnan216
[论文投稿] 超级尴尬的问题 +8 小蜗牛高速 2018-11-11 8/400 2018-11-11 23:42 by yangjinian
[硕博家园] 如何抉择? +12 92林夕一梦 2018-11-09 12/600 2018-11-10 09:52 by greaterwei
信息提示
请填处理意见