查看: 300  |  回复: 6
【悬赏金币】回答本帖问题,作者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
7楼2018-07-17 04:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 moulh 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 如何提高抗压能力 +8 1538933411 2018-07-16 13/650 2018-07-17 00:11 by 1538933411
[考博] 想读博 +5 徐徐杨洋 2018-07-16 10/500 2018-07-16 22:46 by 929368583
[考博] 联系了三个学校的博导一个回信都没有 +16 解梦大触 2018-07-16 17/850 2018-07-16 17:31 by 45℃天空下
[硕博家园] 对上个帖子的回复。 +7 xinkuanti 2018-07-16 10/500 2018-07-16 17:20 by suwen-dian
[基金申请] 前后时间戳中间夹着的这个信息感觉很关键 +9 hmdc 2018-07-16 9/450 2018-07-16 16:50 by Tombob999
[论文投稿] 请教,稿件经过初审,外审,退修,终审,待编现在又变成待复审了,什么情况? +3 owenlzy 2018-07-13 11/550 2018-07-16 15:16 by owenlzy
[基金申请] 未来面上基金会越来越难申请 +17 aloesnow 2018-07-16 18/900 2018-07-16 10:49 by agar
[基金申请] 地学口面上前后缀 +17 marost 2018-07-15 18/900 2018-07-16 09:16 by linjx2000
[基金申请] 无聊的文件名后缀,浮躁的科研圈小丑 +20 stupidgump 2018-07-15 21/1050 2018-07-16 07:11 by 6543yes
[电化学] 请问有想创办锂电池公司的朋友吗? +8 锂电池小毅 2018-07-12 8/400 2018-07-15 23:55 by ypnl
[硕博家园] 博士毕业慎重考虑地方所谓的人才引进 +60 liangyc111 2018-07-11 95/4750 2018-07-15 19:21 by chuguoliuxue
[基金申请] 有没有人成功提交重庆市科委基础与前沿项目任务书? 5+4 tudoutouha 2018-07-10 17/850 2018-07-14 22:20 by 公路小开
[考研] 找研友 +9 行若微尘 2018-07-13 13/650 2018-07-14 10:05 by 行若微尘
[高分子] GPC测多糖类高分子分子量色谱图求助 +3 tlzrx111 2018-07-13 3/150 2018-07-13 21:48 by xiaoLifeiD
[仿真模拟] 为什么ansys这样做只能求出前6阶的模态? +3 萌萌的太阳 2018-07-11 9/450 2018-07-13 20:11 by zhl1998
[考研] 求助王江涛作文~ +3 是老可爱 2018-07-12 3/150 2018-07-13 06:23 by zzu毛毛虫
[有机交流] 非常感谢大家能帮忙一下 +3 1786072632 2018-07-11 4/200 2018-07-12 18:58 by 1786072632
[基金申请] 源代码中没汉子是挂了吗 +29 Solar1986 2018-07-10 40/2000 2018-07-11 16:31 by gdszgd
[考博] 见博导 +12 bruiser000 2018-07-10 20/1000 2018-07-11 09:03 by 303461270
[考博] 工科想读博,但导师侧重实践做项目不太指导论文怎么办? +5 就慢慢来吧 2018-07-10 17/850 2018-07-10 20:39 by pigthree
信息提示
请填处理意见