含镧系原子的赝势基组输入问题

zhang1992

@sobereva 老师：        你好~

我这边想计算一个含有Yb的分子的能级，但是用以下输入文件后出现报错；如果把Yb换成Bi就可以计算，不知道哪里有问题·····
输入：
%chk=P-Yb-1.chk
%mem=20GB
%nprocshared=28
#p opt freq bp86/gen  pseudo=read

opt

3 1
Yb                -2.69400388   -0.21372361   -0.02989929
F                 -3.83400388    1.76081431   -0.02989929
F                 -3.83400388   -2.18826153   -0.02989929
F                 -0.41400388   -0.21372361   -0.02989929

F 0
6-31g(d)
****
Yb 0
MWB59
****

Yb 0
MWB59

输出：

····（前面的我省略了）
General basis read from cards:  (5D, 7F)
Centers:       2      3      4
6-31g(d)
****
Centers:       1
MWB59

***
======================================================================================================
                                       Pseudopotential Parameters
======================================================================================================
  Center     Atomic      Valence      Angular      Power
  Number     Number     Electrons     Momentum     of R      Exponent        Coefficient   SO-Coeffient
======================================================================================================
    1         70           11
                                      G and up
                                                     2        1.0000000        0.00000000    0.00000000
                                      S - G
                                                     2        7.1674000      140.92611800    0.00000000
                                                     2        3.5837000        1.45493500    0.00000000
                                      P - G
                                                     2        6.1460000       76.74522800    0.00000000
                                                     2        3.0730000        7.76620700    0.00000000
                                      D - G
                                                     2        4.2030000       46.06763600    0.00000000
                                                     2        2.1015000       -0.48238900    0.00000000
                                      F - G
                                                     2        8.1312000     -109.74769500    0.00000000
    2          9
                                   No pseudopotential on this center.——还有这里为啥读不出了？
    3          9
                                   No pseudopotential on this center.
    4          9
                                   No pseudopotential on this center.
======================================================================================================
Ernie: Thresh=  0.10000D-02 Tol=  0.10000D-05 Strict=F.
The combination of multiplicity 1 and    35 electrons is impossible.——？这里我也想不通

sobereva

直接写SDD，别写MWB59（赝化59个电子的大核赝势），要不然体系的电子数就成奇数了，显然自旋多重度不能为1

zhang1992

sobereva 发表于 2018-4-25 12:19
直接写SDD，别写MWB59（赝化59个电子的大核赝势），要不然体系的电子数就成奇数了，显然自旋多重度不能为1
...

我直接改成SDD，还是会出现F原子-No pseudopotential on this center，这个现象

ggdh

zhang1992 发表于 2018-4-25 12:54
我直接改成SDD，还是会出现F原子-No pseudopotential on this center，这个现象

因为F原子不需要赝势啊

zhang1992

ggdh 发表于 2018-4-25 16:42
因为F原子不需要赝势啊

我用的6-31g(d),没有用赝势