计算化学公社

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

[新手求助] 求助:左氧氟沙星优化分子结构一直运算也没有结果,求原因以及解决方法。

[复制链接 Copy URL]

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 SamShen 于 2025-12-3 15:17 编辑

求助:纯新手,左氧氟沙星优化分子结构,使用了# opt b3lyp/6-31+g(d,p) geom=connectivity,结果运算到如图所示就没有后续了,没有报错也没有其他反应,就是一直运算也没有结果,请问是什么原因?应该怎么处理呢?(程序和图片见附件)(已经按照大家的指导改了建模,甚至还用大佬给出的建模来算优化,但还是出现了同样的问题)



错误优化.png (80.46 KB, 下载次数 Times of downloads: 2)

我用大佬的进行优化,但计算又进行到这里没有后续了

我用大佬的进行优化,但计算又进行到这里没有后续了

错误频率.png (349 KB, 下载次数 Times of downloads: 0)

我用大佬的直接进行振动频率的计算,结果到这里又运行没有结果了

我用大佬的直接进行振动频率的计算,结果到这里又运行没有结果了

错误优化(自己).png (299.18 KB, 下载次数 Times of downloads: 0)

我自己改了建模进行的优化,运算到这里没后续了

我自己改了建模进行的优化,运算到这里没后续了

zuoyangfushaxing-opt.gjf

3.37 KB, 下载次数 Times of downloads: 1

自己最新搭建的,N原子问题改了

zuoyangfushaxing-opt-taren.gjf

3.38 KB, 下载次数 Times of downloads: 3

下面大佬给出的回复,我拿来接着做的优化

683

帖子

12

威望

2818

eV
积分
3741

Level 5 (御坂)

鸩羽

2#
发表于 Post on 2025-12-2 22:22:00 | 只看该作者 Only view this author
这建模是灾难 需要检查

202512022221402518..png (87.03 KB, 下载次数 Times of downloads: 1)

202512022221402518..png

评分 Rate

参与人数
Participants 2
eV +7 收起 理由
Reason
北大-陶豫 + 5 我很赞同
SharkYYX2025 + 2 233333

查看全部评分 View all ratings

某不知名实验组从苞米地里长出来的计算选手

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-12-2 22:38:36 | 只看该作者 Only view this author
wal 发表于 2025-12-2 22:22
这建模是灾难 需要检查

啊。。。感谢感谢,这个是我修改好多次优化后弄得,有好几次算一半给关了,可能弄乱了,我改一改,感谢!

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2025-12-2 23:21:30 | 只看该作者 Only view this author
我改完了,但还是如图运作到这一步没有后续了,求助!(下面是我改完的程序)
%chk=D:\SCX\zuoyangfushaxing-opt.chk
# opt b3lyp/6-31+g(d,p) geom=connectivity

zuoyangfushaxing-opt

0 2
C                 -0.25123270   -0.62278659    0.38744395
C                  1.08962902   -0.20370509    0.45919949
C                  1.39904239    1.12829477    0.15559686
C                  0.38463467    2.08115882    0.00925992
C                 -0.95469769    1.68361954    0.07364100
C                 -1.27136455    0.32910503    0.23666598
H                  0.63260808    3.10891330   -0.15526845
F                 -1.93917584    2.60146035   -0.03012503
O                 -0.60454140   -2.01027962    0.45016192
N                  2.14581630   -1.16013708    0.88145824
C                  2.85805293    1.51847550   -0.02685887
O                  3.20062847    2.71571180   -0.20829248
C                  3.48782649   -0.86850360    0.31062355
H                  4.17209213   -1.67761306    0.16235851
C                  3.88321600    0.39079390   -0.00943535
C                  5.35540056    0.67578762   -0.36012121
O                  5.76554046    2.00293468   -0.69980440
O                  6.19746083   -0.25926912   -0.34773589
C                 -3.51546643    0.98529970    0.78104693
C                 -3.10304962   -0.41379361   -1.12537518
C                 -4.98891564    0.53775537    0.79087730
H                 -3.41083408    1.85440259    0.16573834
H                 -3.20595537    1.21550214    1.77907206
C                 -4.57647808   -0.86120938   -1.11573273
H                 -2.99822645    0.45539869   -1.74052136
H                 -2.49351072   -1.20214952   -1.51501350
H                 -5.59854848    1.32603701    1.18056986
H                 -5.09355641   -0.33141178    1.40609230
H                 -4.88590086   -1.09116247   -2.11384246
H                 -4.68114968   -1.73041427   -0.50061524
N                 -2.67797731   -0.09780448    0.24583975
N                 -5.41406073    0.22169029   -0.58036396
H                 -5.31626939    1.03406926   -1.15526017
C                 -6.82052247   -0.20569847   -0.57116345
H                 -6.92507382   -1.07490713    0.04399005
H                 -7.43028944    0.58249181   -0.18151319
H                 -7.12993563   -0.43580640   -1.56927900
H                  6.68021395    2.13194933   -0.43847009
C                  0.37759506   -2.82709906    1.09455684
H                  0.12646463   -3.85405602    0.92976482
H                  0.39784551   -2.62922763    2.14591141
C                  1.76270601   -2.52410815    0.49994387
H                  2.48512757   -3.20687393    0.89601410
C                  1.71036442   -2.66904204   -1.03235609
H                  0.99123308   -1.98432274   -1.43102794
H                  2.67423570   -2.45382322   -1.44413765
H                  1.42924415   -3.66967466   -1.28655489

1 2 1.5 6 1.5 9 1.0
2 3 1.5 10 1.0
3 4 1.5 11 1.0
4 5 1.5 7 1.0
5 6 1.5 8 1.0
6 31 1.0
7
8
9 39 1.0
10 13 1.0 42 1.0
11 12 2.0 15 1.0
12
13 14 1.0 15 2.0
14
15 16 1.0
16 17 1.0 18 2.0
17 38 1.0
18
19 21 1.0 22 1.0 23 1.0 31 1.0
20 24 1.0 25 1.0 26 1.0 31 1.0
21 27 1.0 28 1.0 32 1.0
22
23
24 29 1.0 30 1.0 32 1.0
25
26
27
28
29
30
31
32 33 1.0 34 1.0
33
34 35 1.0 36 1.0 37 1.0
35
36
37
38
39 40 1.0 41 1.0 42 1.0
40
41
42 43 1.0 44 1.0
43
44 45 1.0 46 1.0 47 1.0
45
46
47

53

帖子

0

威望

705

eV
积分
758

Level 4 (黑子)

5#
发表于 Post on 2025-12-2 23:40:23 | 只看该作者 Only view this author
SamShen 发表于 2025-12-2 23:21
我改完了,但还是如图运作到这一步没有后续了,求助!(下面是我改完的程序)
%chk=D:\SCX\zuoyangfushaxi ...

32号的N有问题 另外可以检查一下自旋多重度 一般正常的中性体系建模出来为1,不然就是加氢少氢 或者其他原子的成键方式有问题

93

帖子

0

威望

453

eV
积分
546

Level 4 (黑子)

开心小鲨鱼 ORCA的饲养员/食物

6#
发表于 Post on 2025-12-3 00:37:56 | 只看该作者 Only view this author
现在32号N上头多一个H。如果不是特意做成这样的话,我给你删了。这是xtb优化过的,已经转成fch和gjf,可以拿gview打开另存为,然后继续优化
这种复杂结构可以考虑拿chemdraw画,画错了会提醒你,然后保存mol再拿gview打开
opt.gjf (252.51 KB, 下载次数 Times of downloads: 3) opted.fch (392.29 KB, 下载次数 Times of downloads: 1)

评分 Rate

参与人数
Participants 2
eV +10 收起 理由
Reason
北大-陶豫 + 5 赞!
agone997 + 5 你太可爱

查看全部评分 View all ratings

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2025-12-3 12:19:10 | 只看该作者 Only view this author
啊,太感谢了,我查了半天眼花了没看出来。。。感谢大家

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2025-12-3 13:30:56 | 只看该作者 Only view this author
wcy鱼饵 发表于 2025-12-2 23:40
32号的N有问题 另外可以检查一下自旋多重度 一般正常的中性体系建模出来为1,不然就是加氢少氢 或者其他 ...

感谢指导!

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2025-12-3 15:08:24 | 只看该作者 Only view this author
求助:已经按照大家的指导改了建模,甚至还用大佬给出的建模来算优化,但还是出现了同样的问题。【大哭】(用的程序和运行的结果图我放到一楼更新了)

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2025-12-3 15:08:59 | 只看该作者 Only view this author
SharkYYX2025 发表于 2025-12-3 00:37
现在32号N上头多一个H。如果不是特意做成这样的话,我给你删了。这是xtb优化过的,已经转成fch和gjf,可以 ...


求助:已经按照您的指导改了建模,甚至还用您给出的建模来算优化和振动频率,但还是出现了同样的问题。【大哭】(用的程序和运行的结果图我放到一楼更新了)

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2025-12-3 15:09:33 | 只看该作者 Only view this author
wcy鱼饵 发表于 2025-12-2 23:40
32号的N有问题 另外可以检查一下自旋多重度 一般正常的中性体系建模出来为1,不然就是加氢少氢 或者其他 ...

求助:已经按照您的指导改了建模,甚至还用楼下帮助者给出的建模来算优化和振动频率,但还是出现了同样的问题。【大哭】(用的程序和运行的结果图我放到一楼更新了)

53

帖子

0

威望

705

eV
积分
758

Level 4 (黑子)

12#
发表于 Post on 2025-12-3 15:12:20 | 只看该作者 Only view this author
SamShen 发表于 2025-12-3 15:09
求助:已经按照您的指导改了建模,甚至还用楼下帮助者给出的建模来算优化和振动频率,但还是出现了同样的 ...

可以给一下你的报错信息和输入文件吗

19

帖子

0

威望

111

eV
积分
130

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2025-12-3 15:19:54 | 只看该作者 Only view this author
wcy鱼饵 发表于 2025-12-3 15:12
可以给一下你的报错信息和输入文件吗

问题在于没有报错信息就是在这一步一直运行没有后续,输入文件和程序运行的最后显示部分我更新在一楼了,回复这里我发不了

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

14#
发表于 Post on 2025-12-3 16:22:14 | 只看该作者 Only view this author
SamShen 发表于 2025-12-3 15:19
问题在于没有报错信息就是在这一步一直运行没有后续,输入文件和程序运行的最后显示部分我更新在一楼了, ...

程序仍然在运行吗(也就是任务管理器能看到计算进程吗)?注意不是g09w.exe,是l***.exe,***是某个数字。程序的cpu和内存占用情况如何?
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

93

帖子

0

威望

453

eV
积分
546

Level 4 (黑子)

开心小鲨鱼 ORCA的饲养员/食物

15#
发表于 Post on 2025-12-3 16:23:03 | 只看该作者 Only view this author
本帖最后由 SharkYYX2025 于 2025-12-3 16:31 编辑

仔细看此文解决方法:Gaussian的安装方法及运行时的相关问题  http://sobereva.com/439
谈谈弥散函数和“月份”基组 - 思想家公社的门口:量子化学·分子模拟·二次元  http://sobereva.com/119

如果你是卡在L502这一步,可能是算的太慢。是win版吗?可以描述一下你用的核数和电脑性能。win版gaussian速度非常不给力,和玩具一样。我试了下(拿性能差的轻薄办公本,4核心,64W-G16RevA.03),在下一步没算完之前,显示框里头就是半天一动不动,看起来像是卡住了,其实未必,就是算太慢了。如果出错会报错停止,而不是卡着不动。我看这个原子数拿win版算已经有点吃力了。我把你这个结构弥散函数去掉然后做opt,拿笔记本电脑计算还是极慢的。可以考虑换Linux版,配上服务器或者超算平台 new.out (62.94 KB, 下载次数 Times of downloads: 0) new.gjf (2.8 KB, 下载次数 Times of downloads: 0)
也有可能是这个09版本太老,有什么未知bug,可以考虑更新一下其他版本。
可以看看out文件,有没有写在哪一步终止。没写的话估计不是报错,就是算太慢了。
可以把计算关键词保留,算个甲烷或者苯,看看是这个分子的问题还是电脑的问题。或者,把分子结构保留,拿半经验方法算下,看看能不能正常运行结束。这样排除一下是哪方面问题。下面,半经验这个我拿win版算了下,几分钟结束,如果这个也不能正确算出来,那是电脑的原因
  1. #p opt freq pm6d3
复制代码

然后建议目录不要有中文和空格,建议不要装到C盘。(虽然也不是不能用,我电脑就一个盘,只能装在C盘了)

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 18:46 , Processed in 0.191967 second(s), 24 queries , Gzip On.

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