计算化学公社

标题: 含镧系原子的赝势基组输入问题 [打印本页]

作者
Author: zhang1992 时间: 2018-4-25 11:32
标题: 含镧系原子的赝势基组输入问题

@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.——？这里我也想不通

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

作者
Author: zhang1992 时间: 2018-4-25 12:54

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

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

作者
Author: ggdh 时间: 2018-4-25 16:42

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

因为F原子不需要赝势啊

作者
Author: zhang1992 时间: 2018-4-25 16:49

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)