求助molcas基组怎么解收缩

banxia

我想使用解收缩的ANO-RCC基组，应该怎么写输出文件呢？求助各位老师，谢谢

zjxitcc

基组数据库网站上Get Basis Set右边有一个Advanced，在里面选择Format，选择Molcas或Molcas basis_library，底下有3种Uncontract选项，根据你了解到的“解收缩”的程度或者说方案，选择一种，比如谁告诉你的 你就去问谁，或者你从哪篇文献看来的 你就问文献作者。我不知道哪种是最常用的，因为我没用过。

-HMDS-

我这里放一段Molcas的ANO-RCC基组数据：

/Zr.ANO-rcc.Roos.21s18p13d6f4g2h.10s9p8d5f4g2h.
B.O. Roos, R. Lindh, P.-A. Malmqvist, V. Veryazov, P.-O. Widmark. J. Phys. Chem. A 109 (2005) 6575-6579. doi:10.1021/jp0581126
ZIRCONIUM (21s,18p,13d,6f,4g,2h) -> [10s,9p,8d,5f,4g,2h]
* Contraction range: 5s4p2d - 10s9p8d5f4g2h                            *
* Results(eV)                                                          *
* Contraction         CASPT2(Zr,3F)  5F(d3s1) Zr+ 1S(d2s1) Zr- 4F(d3s2)*
* MB: 5s4p2d          -3595.88282151   1.361         6.437      -2.263     *
* DZP:6s5p3d1f        -3596.06195671   0.604         6.536      -0.428     *
* TZP:7s6p4d2f1g      -3596.14301533   0.472         6.593       0.083     *
* QZP:8s7p5d3f2g      -3596.16291192   0.492         6.642       0.213     *
* large:10s9p8d5f4g   -3596.17083268   0.483         6.661       0.397     *
* Primitive:          -3596.17129385   0.484         6.663       0.410     *


你是可以根据他的基组定义直接调用然后选择任意想要的波函数基组大小的
比如你不想要h轨道的QZP就可以用
Zr.ANO-RCC...8S7P5D3F...

8S7P5D3F这里你可以自己选择放进多少，解收缩的我没试过，但是你可以试试这样自定义一下。