查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 突然失去了斗志 +12 17332179410 2018-09-22 14/700 2018-09-23 01:00 by 邯郸赵先生
[教师之家] 从教感悟。 +11 苏小姐的猫 2018-09-21 32/1600 2018-09-22 16:06 by 苏小姐的猫
[硕博家园] 迷茫 +4 丹之布丁 2018-09-22 4/200 2018-09-22 10:50 by 木元雨
[论文投稿] ACS nano审稿意见,拒稿重投 +12 火木丁 2018-09-17 17/850 2018-09-21 22:33 by lugia2
[硕博家园] 有没有天津大学的博士? +9 miss聪花 2018-09-18 14/700 2018-09-21 15:49 by 小邹小邹105
[找工作] 现在想明白了,为啥学校要搞群众路线考核 +3 医学院渣渣 2018-09-20 7/350 2018-09-21 15:37 by 医学院渣渣
[文学芳草园] 一个人的词不达意 +3 zhouyanan91 2018-09-21 3/150 2018-09-21 15:02 by tianliang218
[考研] 有机化学问题 +4 wrgty 2018-09-21 9/450 2018-09-21 10:31 by linyu哒哒哒
[考研] sci对考研有帮助吗 +8 为祖国献石油 2018-09-18 12/600 2018-09-21 10:18 by 前面有樱桃
[硕博家园] 看到身边的人过得很幸福~ (金币+2) +4 yyzx2017 2018-09-17 4/200 2018-09-20 20:41 by 我是小老板
[硕博家园] 计算机研究生学生物信息学方向……一头雾水 +8 虎奇科马 2018-09-20 9/450 2018-09-20 10:37 by 沈晓晨4251
[硕博家园] 机械制造专业女教师在职博士该选什么研究方向? +5 Shmily_dmn 2018-09-17 13/650 2018-09-19 15:21 by skyking20
[有机交流] 聚四氟乙烯盖子怎么会融化了? 5+3 北方sdy 2018-09-19 5/250 2018-09-19 11:14 by 天涯不归路
[论文投稿] 想出文章,创新点找不到 +5 Nstudent 2018-09-19 7/350 2018-09-19 08:55 by Nstudent
[硕博家园] 无人会登临意 +7 A8A8888 2018-09-18 10/500 2018-09-19 07:41 by A8A8888
[基金申请] 基金计划书状态 +9 shyangyang1 2018-09-16 11/550 2018-09-18 21:57 by 独步青云
[论文投稿] 小修变大修,审稿人抓住问题不放怎么办?关键是问题真的很棘手 +5 tanzhijian 2018-09-17 5/250 2018-09-18 07:49 by hsigma
[硕博家园] 很烦躁!!! +18 18838517524 2018-09-16 20/1000 2018-09-17 21:53 by 白宇青田
[考博] 纠结 +14 TIT TAT 2018-09-17 27/1350 2018-09-17 20:42 by TIT TAT
[考博] 请问各位考上博的前辈,研究生必须具有SCI才能考博吗 +5 ttqk211314 2018-09-16 5/250 2018-09-17 08:18 by 2441430010
信息提示
请填处理意见