24小时热门版块排行榜     新春论文润色,领200元京东卡

查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] ieee trans. cyber 被拒后申诉疑惑 +5 15690326186 2019-02-15 16/800 2019-02-16 07:50 by 15690326186
[硕博家园] 相处五年无奈分手 +90 王士冠1 2019-02-14 104/5200 2019-02-16 00:16 by abc876524408
[考研] 总分374求调剂 学前教育学 +6 大可爱吖 2019-02-15 10/500 2019-02-15 22:42 by echo_simon
[考研] 281求调剂 +4 花花花? 2019-02-15 6/300 2019-02-15 22:37 by 花花花?
[论文投稿] SCIE四区好中吗? +6 学者学者 2019-02-14 6/300 2019-02-15 21:17 by 文刃丶
[考研] 19年临床医学考研295,求调剂到B区院校 +3 水瓶座女生 2019-02-15 5/250 2019-02-15 20:13 by 我是个聪明蛋
[硕博家园] SCI论文里的流程图一般用什么软件比较方便做出来 +13 snowinghrb 2019-02-12 15/750 2019-02-15 17:40 by 火舞凤凰
[硕博家园] 转行 +6 C1化学 2019-02-15 6/300 2019-02-15 15:31 by 天蝎水金
[考博] 考博运气成份 +4 军10060112 2019-02-15 9/450 2019-02-15 11:27 by 军10060112
[考研] 等成绩的过程才是最痛苦的 +8 xiuxiu杨 2019-02-14 14/700 2019-02-14 23:04 by 你清醒清醒
[考研] 210求调剂 +4 我是葡萄干 2019-02-14 4/200 2019-02-14 21:43 by 吏逃远远的
[硕博家园] 嗨起来 +7 17332179410 2019-02-09 10/500 2019-02-14 15:58 by 兰慧龙
[考研] 是选择南航,还是选择719研究所 +7 磊?悦520 2019-02-14 11/550 2019-02-14 10:40 by 213070312
[教师之家] 在全国范围内全额拨款的事业单位,总体看包括教师等待遇比公务员低太多了 +18 文化小革命 2019-02-11 24/1200 2019-02-14 09:15 by 儿子毅
[论文投稿] 文章见刊,但是莫名其妙被挂了个不认识的基金号 +18 weiguang7910 2019-02-11 24/1200 2019-02-13 15:02 by weiguang7910
[版块工场] 【早起签到贴】2019年2月13日(q q 48626260 进群申请:小木虫——您个人ID) +57 8475 2019-02-13 57/2850 2019-02-13 13:03 by 携手人生路
[论文投稿] 请问各位 影响因子大于6的天然产物方向杂志有哪些? +10 Wlinger 2019-02-11 10/500 2019-02-12 15:47 by 左波小小虫
[硕博家园] 最近在看老版天龙八部,发现人物性格塑造的很好。 +5 玲珑文哲 2019-02-11 6/300 2019-02-12 09:03 by 露雨夹雪A
[考研] 关于考研调剂 +9 含含嘻嘻 2019-02-09 10/500 2019-02-11 20:49 by 实验迷在此
[硕博家园] 突然想到马上收假了 (金币+2) +10 DEhiyo 2019-02-10 22/1100 2019-02-10 16:56 by 漫步云端LWG
信息提示
请填处理意见