计算化学公社

标题: 如何给特定原子添加ma-方式的弥散函数 [打印本页]
作者
Author: cokie    时间: 2024-1-22 15:18
标题: 如何给特定原子添加ma-方式的弥散函数
各位老师好，我根据博文《给def2以ma-方式加弥散函数的Gaussian格式的基组定义文件（含所有def2支持的元素）》（http://sobereva.com/509）来调用“ma-def2-TZVP”。不过我想做的是给特定几个原子加弥散（水解的羧基），其他原子用def2-TZVP。遇到两个问题。我的输入文件如下：
#p PBE1PBE/gen em=gd3bj TD(singlets,NSTATES=15,ROOT=1) IOp(9/40=4)

0 1
……
空行
1 2 5 6 29 30 0
@C:\orcajob\HP\C3T-PBE0-diffuse-def2tzvp\ma-TZVP.txt
****
3-4 7-28 31-110 0
def2tzvp
****
空行
空行
我才用了sob老师提供的“file”里的.txt文件，但是提交任务后报错显示：定义的基组错误，从“-H”开始这一行就无法读取。之后我选择钯.txt里的“-C还有-O”的负号去了（加弥散的原子），可以正确读取并计算，但程序显示：
****
Warning:  center   9 has no basis functions!
Warning:  center  10 has no basis functions!
Warning:  center  23 has no basis functions!
Warning:  center  24 has no basis functions!
Warning:  center  35 has no basis functions!
Warning:  center  36 has no basis functions!
Warning:  center  37 has no basis functions!
Warning:  center  38 has no basis functions!
Warning:  center  42 has no basis functions!
Warning:  center  43 has no basis functions!
Warning:  center  45 has no basis functions!
Warning:  center  46 has no basis functions!
Warning:  center  48 has no basis functions!
Warning:  center  49 has no basis functions!
Warning:  center  51 has no basis functions!
Warning:  center  52 has no basis functions!
Warning:  center  54 has no basis functions!
Warning:  center  55 has no basis functions!
Warning:  center  57 has no basis functions!
Warning:  center  58 has no basis functions!
Warning:  center  60 has no basis functions!
Warning:  center  61 has no basis functions!
Warning:  center  63 has no basis functions!
Warning:  center  64 has no basis functions!

…………………………………………………………
是我不能用这种方式来给局部加弥散函数吗？正确的方式应该是什么呢？恳请老师们赐教！谢谢！
作者
Author: wzkchem5    时间: 2024-1-22 17:35
报错的这些center对应的原子分别是什么元素？
作者
Author: cokie    时间: 2024-1-22 17:47
本帖最后由 cokie 于 2024-1-22 17:48 编辑
wzkchem5 发表于 2024-1-22 17:35
报错的这些center对应的原子分别是什么元素？

C，H，O，N元素都有。谢谢老师，我目前已经能正确运行了（虽然不知道我解决的方案对不对）
我把自定义基组的地方改写作了这样：
  1. ……………………
  2. 空行
  3. 3-4 7-28 31-110 0
  4. def2tzvp
  5. ****
  6. 1 2 5 6 0
  7. S   6   1.00
  8.   27032.3826310              0.21726302465E-03
  9.    4052.3871392              0.16838662199E-02
  10.     922.32722710             0.87395616265E-02
  11.     261.24070989             0.35239968808E-01
  12.      85.354641351            0.11153519115
  13.      31.035035245            0.25588953961
  14. S   2   1.00
  15.      12.260860728            0.39768730901
  16.       4.9987076005           0.24627849430
  17. S   1   1.00
  18.       1.1703108158           1.0000000
  19. S   1   1.00
  20.       0.46474740994          1.0000000
  21. S   1   1.00
  22.       0.18504536357          1.0000000
  23. P   4   1.00
  24.      63.274954801            0.60685103418E-02
  25.      14.627049379            0.41912575824E-01
  26.       4.4501223456           0.16153841088
  27.       1.5275799647           0.35706951311
  28. P   1   1.00
  29.       0.52935117943          1.0000000
  30. P   1   1.00
  31.       0.17478421270          1.0000000
  32. D   1   1.00
  33.       2.31400000             1.0000000
  34. D   1   1.00
  35.       0.64500000             1.0000000
  36. F   1   1.00
  37.       1.42800000             1.0000000
  38. S   1   1.00
  39.   0.616817878567E-01   1.000
  40. P   1   1.00
  41.   0.582614042333E-01   1.000
  42. ****
  43. 29 30 0
  44. S   6   1.00
  45.   13575.3496820              0.22245814352E-03
  46.    2035.2333680              0.17232738252E-02
  47.     463.22562359             0.89255715314E-02
  48.     131.20019598             0.35727984502E-01
  49.      42.853015891            0.11076259931
  50.      15.584185766            0.24295627626
  51. S   2   1.00
  52.       6.2067138508           0.41440263448
  53.       2.5764896527           0.23744968655
  54. S   1   1.00
  55.       0.57696339419          1.0000000
  56. S   1   1.00
  57.       0.22972831358          1.0000000
  58. S   1   1.00
  59.       0.95164440028E-01            1.0000000
  60. P   4   1.00
  61.      34.697232244            0.53333657805E-02
  62.       7.9582622826           0.35864109092E-01
  63.       2.3780826883           0.14215873329
  64.       0.81433208183          0.34270471845
  65. P   1   1.00
  66.       0.28887547253          1.0000000
  67. P   1   1.00
  68.       0.10056823671          1.0000000
  69. D   1   1.00
  70.       1.09700000             1.0000000
  71. D   1   1.00
  72.       0.31800000             1.0000000
  73. F   1   1.00
  74.       0.76100000             1.0000000
  75. S   1   1.00
  76.   0.317214800093E-01   1.000
  77. P   1   1.00
  78.   0.335227455700E-01   1.000
  79.   ****
  80. 空行
  81. 空行
复制代码


就可以了，我猜测可能是因为sob老师的.txt文件是针对元素定义的，但是我加弥散的六个原子是C和O元素，和sob老师写的.txt里的基组定义冲突了。所以我就分别把里面对O和C的定义copy出来分开加的，现在就能正常跑了。（应该这种处理方式是对的吧）
作者
Author: sobereva    时间: 2024-1-22 23:00
cokie 发表于 2024-1-22 17:47
C，H，O，N元素都有。谢谢老师，我目前已经能正确运行了（虽然不知道我解决的方案对不对）
我把自定义基 ...

对



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3