应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 金属氢氧化物优化出现错误
61/1返回列表
查看: 2175  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

薛定谔的糖果

木虫 (小有名气)

[求助] 金属氢氧化物优化出现错误 已有1人参与

Rotational constants (GHZ):      0.3554359      0.2987240      0.2117731
General basis read from cards:  (5D, 7F)
Warning:  center   1 has no basis functions!
Warning:  center   3 has no basis functions!
Warning:  center   4 has no basis functions!
Warning:  center   5 has no basis functions!
Warning:  center   7 has no basis functions!
Warning:  center   8 has no basis functions!
Warning:  center   9 has no basis functions!
Warning:  center  11 has no basis functions!
Warning:  center  13 has no basis functions!
Warning:  center  14 has no basis functions!
Warning:  center  15 has no basis functions!
Warning:  center  16 has no basis functions!
Warning:  center  17 has no basis functions!
Warning:  center  18 has no basis functions!
Warning:  center  19 has no basis functions!
Warning:  center  20 has no basis functions!
Warning:  center  21 has no basis functions!
Warning:  center  22 has no basis functions!
Warning:  center  23 has no basis functions!
Warning:  center  24 has no basis functions!
There are    88 symmetry adapted basis functions of A   symmetry.
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
There are  96 occupied orbitals but only  88 basis functions!
Error termination via Lnk1e in /opt/gaussian/g09C01/g09/l301.exe at Wed Aug  1 17:24:29 2018.
Job cpu time:  0 days  0 hours  0 minutes  1.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

错误信息见上
输入文件见下
%nprocshared=4
%mem=10GB
%chk=Ti.chk
# opt freq b3lyp/gen pop=npa geom=connectivity

****

0 1
O                 -4.37960000   -1.09060000    1.91760000
Ti                -6.19020000   -0.86820000    1.21900000
O                 -5.18760000   -2.06090000    0.00280000
O                 -6.59260000   -0.15770000   -0.56220000
O                 -5.85290000    1.07090000    1.35700000
Ti                -4.38140000    2.38800000    1.23670000
O                 -4.96810000    1.92130000   -0.57810000
O                 -3.68150000    1.18400000    2.60610000
O                 -2.81490000    1.92220000    0.13770000
Ti                -3.76880000   -2.04330000   -1.36930000
O                 -4.25700000   -0.17890000   -1.69090000
Ti                -1.62240000    0.37080000    0.19150000
O                 -1.35900000    1.15590000   -1.60450000
O                  0.23970000    0.92970000    0.00530000
O                 -2.98150000   -2.69440000    0.30690000
O                 -2.10280000   -1.02010000   -1.10440000
H                 -4.07500000   -1.97290000    1.81190000
H                 -6.63910000    0.78110000   -0.54730000
H                 -4.26000000    1.72430000   -1.16080000
H                 -3.11370000    1.64530000    3.19770000
H                 -5.15860000   -0.08500000   -1.93050000
H                 -2.03280000    0.85420000   -2.18710000
H                  0.38370000    1.30560000   -0.84570000
H                 -3.47110000   -3.43600000    0.61810000

1 2 1.0 17 1.0
2 3 1.0 4 1.0 5 1.0
3 10 1.0
4 18 1.0
5 6 1.0
6 7 1.0 8 1.0 9 1.0
7 19 1.0
8 20 1.0
9 12 1.0
10 11 1.0 15 1.0 16 1.0
11 21 1.0
12 13 1.0 14 1.0 16 1.0
13 22 1.0
14 23 1.0
15 24 1.0
16
17
18
19
20
21
22
23
24

Ti 0
Lanl2dz
****

O 0
6-31+g(d,p)
****

H 0
6-31+g(d,p)
****
回复此楼

» 猜你喜欢
1楼 2018-08-02 11:10:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ABetaCare

捐助贵宾 (著名写手)
http://sobereva.com/60
把赝势基组输入方法写对
赞一下(1人) 回复此楼
懂一点点弱相互作用分析的电脑小白、EQ贫困户、高SQ烧脑症、自然科学忠实信徒、日常CSI下饭族。
2楼2018-08-02 11:36:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
薛定谔的糖果: 金币+30, ★★★★★最佳答案 2018-08-03 15:59:05
基组部分写法有问题
1. 非赝势部分基组不同原子之间不应加入空行 加入空行高斯会默认空行之后是赝势部分，这就是为什么错误提示说的是 O 和 H 的基组找不到， 这就是 在Ti 基组末尾的 **** 之后加入了空行导致的。

2. LanL2DZ是赝势基组，需要把赝势部分单独写出来。 做法就是在写完非赝势部分的 最后一个**** 之后加入一个空行 之后加上赝势部分

3 另外 如果O 和 H用同样基组 是能够写在一起的。

input的大概写法应该是
%nprocshared=4
%mem=10GB
%chk=Ti.chk
# opt freq b3lyp/gen pop=npa

****

0 1
O                 -4.37960000   -1.09060000    1.91760000
Ti                -6.19020000   -0.86820000    1.21900000
O                 -5.18760000   -2.06090000    0.00280000
O                 -6.59260000   -0.15770000   -0.56220000
O                 -5.85290000    1.07090000    1.35700000
Ti                -4.38140000    2.38800000    1.23670000
O                 -4.96810000    1.92130000   -0.57810000
O                 -3.68150000    1.18400000    2.60610000
O                 -2.81490000    1.92220000    0.13770000
Ti                -3.76880000   -2.04330000   -1.36930000
O                 -4.25700000   -0.17890000   -1.69090000
Ti                -1.62240000    0.37080000    0.19150000
O                 -1.35900000    1.15590000   -1.60450000
O                  0.23970000    0.92970000    0.00530000
O                 -2.98150000   -2.69440000    0.30690000
O                 -2.10280000   -1.02010000   -1.10440000
H                 -4.07500000   -1.97290000    1.81190000
H                 -6.63910000    0.78110000   -0.54730000
H                 -4.26000000    1.72430000   -1.16080000
H                 -3.11370000    1.64530000    3.19770000
H                 -5.15860000   -0.08500000   -1.93050000
H                 -2.03280000    0.85420000   -2.18710000
H                  0.38370000    1.30560000   -0.84570000
H                 -3.47110000   -3.43600000    0.61810000

Ti 0
Lanl2dz
****
O H 0
6-31+g(d,p)
****

Ti 0
Lanl2dz
赞一下 回复此楼
3楼2018-08-03 00:12:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

薛定谔的糖果

木虫 (小有名气)
引用回帖:
3楼: Originally posted by paramecium86 at 2018-08-03 00:12:05
基组部分写法有问题
1. 非赝势部分基组不同原子之间不应加入空行 加入空行高斯会默认空行之后是赝势部分，这就是为什么错误提示说的是 O 和 H 的基组找不到， 这就是 在Ti 基组末尾的 **** 之后加入了空行导致的。 ...

你好，我按照您写的您写的内容运行了一下，该开始的时候没有报错，但我今天查看结果的时候又提示了错误：Error
Error termination via Lnk1e in /opt/gaussian/g09C01/g09/l607
能再帮我解答一下吗？感谢。
赞一下 回复此楼
4楼2018-08-09 14:16:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
4楼: Originally posted by 薛定谔的糖果 at 2018-08-09 01:16:14
你好，我按照您写的您写的内容运行了一下，该开始的时候没有报错，但我今天查看结果的时候又提示了错误：Error
Error termination via Lnk1e in /opt/gaussian/g09C01/g09/l607
能再帮我解答一下吗？感谢。...

这个错误一般是因为 g09自带的NBO3.1软件有 bug导致。 如果能用编译进高级版本NBO5.9 的高斯软件，可能这问题就解决了。如果没法换软件那么换换基组也是可以的
这个error提示完整的是如下这样。

Subroutine NAOANL could not find a S-type valence orbital on atom Ti    2.
IVAL :  1  1  1  0     M :  1     LA :  1
Error
Error termination via Lnk1e in /home/yc2698v4/g16/l607.exe

那么既然问题出在Ti的lanl2dz 赝势基组上。那么尝试换一个基组给ti。 我试了试 def2tzvp 计算就能正常结束。
我把我用的input和算出来的output都放附件了
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : test_08022018_redo3.log
    • 2018-08-10 03:44:31, 2.87 M
  • 附件 2 : test_08022018_redo3.com
    • 2018-08-10 03:44:32, 1.56 K

» 本帖已获得的红花（最新10朵）

薛定谔的糖果
5楼2018-08-10 03:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

薛定谔的糖果

木虫 (小有名气)
送红花一朵
引用回帖:
5楼: Originally posted by paramecium86 at 2018-08-10 03:45:56
这个错误一般是因为 g09自带的NBO3.1软件有 bug导致。 如果能用编译进高级版本NBO5.9 的高斯软件，可能这问题就解决了。如果没法换软件那么换换基组也是可以的
这个error提示完整的是如下这样。

Subroutine N ...

非常感谢！

发自小木虫IOS客户端
回复此楼
6楼2018-08-10 05:03:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 薛定谔的糖果 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 8rmuugja8q 2026-02-22 6/300 2026-02-23 06:39 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:34 by w4l55oybr1
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 10/500 2026-02-23 04:58 by 5jlh3qtdvx
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 6/300 2026-02-23 02:08 by 5jlh3qtdvx
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +5 successhe 2026-02-16 6/300 2026-02-22 21:38 by 山西悬空寺空悬�
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[教师之家] 为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 5/250 2026-02-22 21:15 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +9 瞬息宇宙 2026-02-20 12/600 2026-02-21 10:39 by 欢乐颂叶蓁
信息提示
关闭
请填处理意见
关闭