24小时热门版块排行榜    

查看: 537  |  回复: 2

ywfsxh

金虫 (正式写手)


[交流] 【求助】关于MCSCF计算?

我想计算Xe2+的激发态能量,$DET里几个参量不会取,NCORE,NACT,NELS?哪位高手能指导一下啊?谢谢,
下面是我的输入文件,是不是SBKJC这个基在这里不合适啊?

$CONTRL SCFTYP=MCSCF RUNTYP=ENERGY ICHARG=2 MULT=3 PP=SBKJC $END
$SYSTEM TIMLIM=525600 MEMORY=1000000 $END
$BASIS GBASIS=SBKJC $END
$DATA
Xe 2+ MCSCF
DNH 2

Xe    54.0    0.00000     0.00000     0.00000
$END
$GUESS GUESS=MOREAD NORB=8 $END
$DET NCORE=2 NACT=2 NELS=2 NSTATE=1 $END
$VEC
1  1 9.08475190E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.07393411E-01
1  2 0.00000000E+00 0.00000000E+00 0.00000000E+00
2  1 0.00000000E+00 9.56512106E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00
2  2 6.08700309E-02 0.00000000E+00 0.00000000E+00
--------------------------------------------------------------------------------------------------------
8  1-1.58016604E+00-0.00000000E+00-0.00000000E+00-0.00000000E+00 1.81953800E+00
8  2-0.00000000E+00-0.00000000E+00-0.00000000E+00
$END

[ Last edited by ywfsxh on 2011-2-7 at 23:22 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ywfsxh

金虫 (正式写手)


引用回帖:
Originally posted by ywfsxh at 2011-02-07 23:10:35:
我想计算Xe2+的激发态能量,$DET里几个参量不会取,NCORE,NACT,NELS?哪位高手能指导一下啊?谢谢,
下面是我的输入文件,是不是SBKJC这个基在这里不合适啊?

$CONTRL SCFTYP=MCSCF RUNTYP=ENERGY ICHARG= ...

SBKJC可不可以用在$ORMAS计算里呢?
2楼2011-02-08 21:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwwolfchj

金虫 (正式写手)


Re

ywfsxh(金币+2): 2011-02-12 16:08:48
先做HF计算,将HF的轨道粘贴到MCSCF里面,NORB取多少看你HF分子的轨道有多少个,SBKJC基组没用过,但Xe2+那么大,仅有8个NORB吗?NCORE是冻结核轨道数,NACT是你要考虑的活性电子数, NELS是活性电子, NSTAT是要求的电子态能级,1好像默认是基态,如果考虑对称性还要加ISTSYM,如果要考虑某个态的N阶特性,还要加其它参数。你先参考自带例子8和9,搞清楚了对全电子就没问题了,再做ECP的。并且你的多重读是3,为开窍层结构,还要加单电子方面的轨道参数,可能还要考虑自选情况。具体也可看GAMESS的Reference里面也有例子。好运。
3楼2011-02-12 14:04:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ywfsxh 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见