24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 1047  |  回复: 4

逆流的船123

新虫 (初入文坛)

[求助] molpro用ccsd(t)/aVDZ方法计算水团簇的单点能的问题。求助!

molpro用ccsd(t)/aVDZ方法计算水团簇的单点能,文献上说计算相关能计算(corerlation caculation)中,对基组要frozen氧原子的1s轨道,请问aVDZ基组是默认冻结氧原子的1s轨道还是要另外加卡指定?加的话,怎么写卡?
memory,1000,m
symmetry,nosym,noorient
basis=avdz
geomtyp=xyz
geometry={
18           ! number of atoms
This is an example of geometry input for water with an XYZ file
H,4.14240300, 0.82298200,-1.17574900
H,1.31553400, 1.42775900,-0.16018200
H,1.44169300,-0.12802500,-0.69386900
H,4.33058500, 1.25417600, 0.28721000
H,0.05125400, 0.43083200,-0.23695400
H,0.38312300,-3.69823800, 1.03877500
H,2.16938400, 2.77736400, 1.49885300
H,2.81196100,-1.86388100,-1.22868200
H,1.29933900,-0.99768300, 0.56110400
H,1.63939800,-2.26335700,-0.29024900
H,2.15057100, 0.25560300, 0.12409900
H,2.27722900, 3.46647700, 0.12244100
O,0.85952600, 0.69680700,-0.67592500
O,2.19058200,-1.51624800,-0.57481300
O,1.25799600,-0.03988400, 0.41900900
O,2.12749800, 2.61077000, 0.54744000
O,0.15381800,-2.83792400, 0.63545300
O,3.67803200, 0.90097900,-0.33201400
}
{rhf;wf,59,1,1}
uccsd(t)
回复此楼
不抛弃不放弃
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

逆流的船123

新虫 (初入文坛)

不抛弃不放弃
2楼2016-04-29 09:54:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

逆流的船123

新虫 (初入文坛)

3楼2016-04-29 16:13:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清新俊逸

木虫 (正式写手)

这么写默认采用frozen core 近似

发自小木虫Android客户端
4楼2016-04-29 23:18:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清新俊逸

木虫 (正式写手)

你的表述也有些问题。冻结1S轨道跟基函数没关系,是指计算相关能的时候不包括core-valence correlation,这是所谓的冻芯近似。当然,不同软件的具体处理方法不一样,比如,gaussian和molpro的ccsd算的结果就有微小的差别。如果想计算芯电子的影响,就在cc下面加选项core,0。这时候用avdz基函数计算也是可以的,但用aug-cc-pcvdz这样的基函数会更好。

发自小木虫Android客户端
5楼2016-04-29 23:25:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 逆流的船123 的主题更新
信息提示
请填处理意见