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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » CASSCF计算如何计算CCSD(T)？
61/1返回列表
查看: 2552  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tanghao86

金虫 (初入文坛)

[求助] CASSCF计算如何计算CCSD(T)？已有2人参与

对于一个金属化合物，不知道该化合物的基态是闭壳层单重态，还是开窍成单重态，或三重态。对这个化合物做了CASSCF（10，6）计算，请教大家，我因为不知道有该化合物在单重态上会有多少个激发态，那么在计算时 NROOT 该等于几呢？
另外，好像CASSCF算完一个组态会出现一个HF能量，但是我想要CCSD(T)的能量，我该怎么做呢？是先进行CASSCF计算，再把计算出来的每个组态进行CCSD(T)单点计算 还是 可以在计算CASSCF时就加上ccsd(t)关键词？

请大家指教！万分感谢！
回复此楼

» 收录本帖的淘帖专辑推荐

经验

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-10-28 09:42:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
1,先进行基态计算，为什么会不能确定稳定基态？?
2,需要多少个态取决于你感兴趣的跃迁是否出现,看跃迁能量或者轨道分析
3,不知你用的什么软件，如果是molcas的话，ccsdt还必须得先做rasscf，具体你查看手册，有详细例子。
赞一下 回复此楼
2楼2014-10-29 06:50:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangning198511

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
2楼: Originally posted by gmy1990 at 2014-10-29 06:50:10
1,先进行基态计算，为什么会不能确定稳定基态？?
2,需要多少个态取决于你感兴趣的跃迁是否出现,看跃迁能量或者轨道分析
3,不知你用的什么软件，如果是molcas的话，ccsdt还必须得先做rasscf，具体你查看手册，有详 ...

1.对于单参考方法，如DFT，CCSD等，开窍层单重态不能描述（需要两个斯莱特行列式），所以就要使用多参考CASSCF等方法来计算。（注意活性空间选取）

2.这个可以多取几个态，因为你只是要寻求最低能量的态。对于阿贝尔群体系而言，
CAS（10，6）总共有：２s+1 : 1,3,5,7  states Number :4950,6930,2310,210
计算公式可以在molcas公司提供的材料里找到。可以看出有非常多的态。

3.CCSD(t)是以ＨＦ为参考的，以ＣＡＳＳＣＦ为参考那就是ＭＲＣＣ，这个目前都没有什么成熟的软件（有个程序叫ＭＲＣＣ，里面有ＭＲＣＣ模块，目前只能计算非常小的体系）。所以建议你使用ＣＡＳＰＴ２　，具有可以参考软件手册。
赞一下 回复此楼
3楼2014-10-29 08:56:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangning198511

木虫 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by gmy1990 at 2014-10-29 06:50:10
1,先进行基态计算，为什么会不能确定稳定基态？?
2,需要多少个态取决于你感兴趣的跃迁是否出现,看跃迁能量或者轨道分析
3,不知你用的什么软件，如果是molcas的话，ccsdt还必须得先做rasscf，具体你查看手册，有详 ...

3,不知你用的什么软件，如果是molcas的话，ccsdt还必须得先做rasscf，具体你查看手册，有详细例子。

ＭＯＬＣＡＳ是利用ＲＡＳＳＣＦ模块计算ＨＦ能量，它的活性空间只是１，所以就是ＲＯＨＦ方法，并不是利用ＣＡＳＳＣＦ参考态计算ＣＣＳＤ（t）. 在ＣＣ计算方面，比较推荐ＭＯＬＰＲＯ和ＣＦＯＵＲ。
赞一下 回复此楼
4楼2014-10-29 09:04:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tanghao86

金虫 (初入文坛)
非常感谢你们的提示。我现在用高斯做CASSCF, 用高斯做ccsd(t)单点，但是单重态正常结束，而在进行三重态ccsd(t)单点计算时，出现了l913错误，错误信息如下，请问有什么办法解决吗？

我本想用molpro做ccsd(T) 但是我不知道该如何写基组的输入文件，还想请问 如果我想对金属Ru 定义SDD基组和SDD赝势; 对P S原子定义SDD基组; 对C H定义6-31g基组，我该怎么写molpro输入文件中基组部分呢？
**********************
DD1Dir will call FoFMem   1 times, MxPair=      4294
NAB=  1443 NAA=   741 NBB=   666.
Norm of the A-vectors is  7.8132901D-05 conv= 1.00D-05.
RLE energy=       -1.4618693533
DE(Corr)=  -1.4618698     E(CORR)=     -2543.9889259     Delta=-1.57D-07
NORM(A)=   0.15023207D+01
*************
*MAX. CYCLES*
*************
Dominant configurations:
***********************
Spin Case        I    J    A    B          Value
    BB           70        78           -0.112989D+00
    BB           71        78            0.257942D+00
    BB           72        78           -0.160335D+00
    BB           73        78           -0.162565D+00
    BB           76        79            0.288029D+00
    BB           77        78            0.108302D+00
Largest amplitude= 2.88D-01
Error termination via Lnk1e in /usr/local/g09_D01/g09/l913.exe at Fri Oct 31 07:56:46 2014.
Job cpu time:       5 days  5 hours  3 minutes 33.8 seconds.
File lengths (MBytes):  RWF=  18105 Int=      0 D2E=      0 Chk=     10 Scr=      1
赞一下 回复此楼
5楼2014-11-01 00:51:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tanghao86

金虫 (初入文坛)
引用回帖:
2楼: Originally posted by gmy1990 at 2014-10-29 06:50:10
1,先进行基态计算，为什么会不能确定稳定基态？?
2,需要多少个态取决于你感兴趣的跃迁是否出现,看跃迁能量或者轨道分析
3,不知你用的什么软件，如果是molcas的话，ccsdt还必须得先做rasscf，具体你查看手册，有详 ...

我用dft方法计算是闭壳层单重态为基态，但是实验室报道是开窍层单重态为基态，我现在不能与实验一致，所以要用CASSCF和CCSD(T)计算。
赞一下 回复此楼
6楼2014-11-01 00:53:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tanghao86 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭