计算化学公社

 找回密码 Forget password
 注册 Register
Views: 13541|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助,高斯D.01版modredundant似乎有bug

[复制链接 Copy URL]

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

各位老师,高斯似乎有个很严重的bug,我用高斯D.01版本时,mod redundant会发生非常奇怪的问题。
输入文件如下:
#p opt=modredundant hf/3-21g geom=connectivity

Title Card Required

0 1
H                  0.00000000    0.00000000    0.00000000
H                  0.00000000    0.00000000    1.20000000

1
2

1 2  1.41 F
一输入高斯中,modredundant立刻会被读错。
我希望高斯将键长固定住(freeze),但是无论如何都会被读成Build。
输出文件的开头是这个样子的:
----------------------------------------------
#p opt=modredundant hf/3-21g geom=connectivity
----------------------------------------------
1/18=120,19=15,38=1,57=2/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=5,11=9,16=1,25=1,30=1,71=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7//1,2,3,16;
1/18=20,19=15/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=5,11=9,16=1,25=1,30=1,71=1/1,2,3;
4/5=5,16=3,69=1/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/18=20,19=15/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
Leave Link    1 at Sat Jan 03 16:34:04 2015, MaxMem=           0 cpu:         0.0
(Enter C:\G09W\l101.exe)
-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 1
H                     0.        0.        0.
H                     0.        0.        1.2

Add virtual bond connecting atoms H2         and H1         Dist= 2.27D+00.
The following ModRedundant input section has been read:
B       1       2 1.4100 B                                                   
不知道为什么,不论如何改,modredundant一直都是 build 模式,没办法用freeze,也没有办法relaxed PES scan。

请问,大家有什么好的解决方法吗? 比如 IOP命令什么的?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2015-1-4 08:17:00 | 只看该作者 Only view this author
modredundant是D.01大家经常提起的bug。用之前版本或者等下一个版本吧(估计快出了)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2015-1-4 09:31:14 | 只看该作者 Only view this author
谢谢sob,原来如此。 目前准备来B.01版的文件覆盖掉一部分D.01文件,看看能否解决。

9

帖子

0

威望

108

eV
积分
117

Level 2 能力者

4#
发表于 Post on 2015-1-5 10:16:17 | 只看该作者 Only view this author
以前用D01也有类似的情况,不过如果先手动调整键长,下面直接写1 2 F,或者直接去掉1.41估计是可以的

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2015-1-5 23:32:18 | 只看该作者 Only view this author
lou_cheung 发表于 2015-1-5 10:16
以前用D01也有类似的情况,不过如果先手动调整键长,下面直接写1 2 F,或者直接去掉1.41估计是可以的

感谢!这样写确实能行。

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

6#
发表于 Post on 2015-1-6 21:15:31 来自手机 | 只看该作者 Only view this author
楼主的版本还能用,我连装都装不上,覆盖老版本都不行

230

帖子

0

威望

940

eV
积分
1170

Level 4 (黑子)

7#
发表于 Post on 2015-11-6 05:01:34 | 只看该作者 Only view this author
参看高斯的更新说明,这个据说不是Gaussian 09 D.01的bug,而是“特性”。因为指定坐标可能导致几何关系被破坏,F,S关键字都不允许显式制定冗余坐标的值了。以后可能连"B"都取消掉。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

16

帖子

0

威望

117

eV
积分
133

Level 2 能力者

8#
发表于 Post on 2017-5-9 15:52:27 | 只看该作者 Only view this author
stecue 发表于 2015-11-6 05:01
参看高斯的更新说明,这个据说不是Gaussian 09 D.01的bug,而是“特性”。因为指定坐标可能导致几何关系被 ...

这个是真的吗?

4

帖子

0

威望

4184

eV
积分
4188

Level 6 (一方通行)

9#
发表于 Post on 2017-5-9 16:28:27 | 只看该作者 Only view this author
D.01版本更新的内容之一啊。。。
Opt=ModRedundant 不再支持指定一个特定的初始值。取而代之的是输入结构需给出所期望的初始值,并且坐标定义中不需要输入此值。

16

帖子

0

威望

117

eV
积分
133

Level 2 能力者

10#
发表于 Post on 2017-5-9 19:23:34 | 只看该作者 Only view this author
kulaomega 发表于 2017-5-9 16:28
D.01版本更新的内容之一啊。。。
Opt=ModRedundant 不再支持指定一个特定的初始值。取而代之的是输入结构 ...

非常感谢!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 05:30 , Processed in 0.236623 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list