“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[Molclus] 使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索

  [复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

331#
 楼主 Author| 发表于 Post on 2024-12-3 01:13:59 | 只看该作者 Only view this author
hujj29 发表于 2024-12-2 22:52
请教sob老师,寻找最低能量的过渡态构象有什么有效、常用的方法或流程吗

同一个问题别在两处问
http://bbs.keinsci.com/thread-50180-1-1.html我已经回答了
北京科音自然科学研究中心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

72

帖子

0

威望

389

eV
积分
461

Level 3 能力者

332#
发表于 Post on 2024-12-17 18:44:55 | 只看该作者 Only view this author
用isostat时出现错误,我复合物的原子数目598个,1000个结构,是什么原因?

202412171843326639..png (54.76 KB, 下载次数 Times of downloads: 62)

202412171843326639..png

72

帖子

0

威望

389

eV
积分
461

Level 3 能力者

333#
发表于 Post on 2024-12-17 19:03:00 | 只看该作者 Only view this author
wilely 发表于 2024-12-17 18:44
用isostat时出现错误,我复合物的原子数目598个,1000个结构,是什么原因?

解决了

14

帖子

0

威望

151

eV
积分
165

Level 3 能力者

334#
发表于 Post on 2025-4-23 13:52:47 | 只看该作者 Only view this author
在重复第四步:用Molclus调用Gaussian和ORCA得到水环境下的高精度自由能》时:
完全按照说明设置,出现以下错误:
Error: This Gaussian task did not perform geometry optimization at all! Please check if you omitted "opt" keyword or there is any problem in the template file
Press ENTER button to exit program
使用的版本是G16

已上传运行附件


gau.gjf

1.99 KB, 下载次数 Times of downloads: 4

gau.out

6.59 KB, 下载次数 Times of downloads: 4

template.gjf

80 Bytes, 下载次数 Times of downloads: 1

1296

帖子

6

威望

2761

eV
积分
4177

Level 6 (一方通行)

傻傻的木瓜

335#
发表于 Post on 2025-4-23 14:56:45 | 只看该作者 Only view this author
caoyx0612 发表于 2025-4-23 13:52
在重复《第四步:用Molclus调用Gaussian和ORCA得到水环境下的高精度自由能》时:[ ...

在提交计算前请先检查化学合理性,不能默认从某些数据库下载的结构文件可以直接用。这个gjf文件里的分子缺失氢原子且存在O9-O35的异常接触、O14-C38-C39-C40的诡异直线等等,是从哪来的、叫什么名?
√546=23.36664289109

14

帖子

0

威望

151

eV
积分
165

Level 3 能力者

336#
发表于 Post on 2025-4-23 15:08:04 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-4-23 14:56
在提交计算前请先检查化学合理性,不能默认从某些数据库下载的结构文件可以直接用。这个gjf文件里的分子 ...

明白了,感谢您的解答。我是初学者,按照教程复现案例。因为用的高斯版本是新版的,原来的案例参数可能要修改,但是我对这些参数还不会调整,只能慢慢学习了

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

337#
 楼主 Author| 发表于 Post on 2025-4-23 15:09:04 | 只看该作者 Only view this author
caoyx0612 发表于 2025-4-23 13:52
在重复《第四步:用Molclus调用Gaussian和ORCA得到水环境下的高精度自由能》时:[ ...

由gau.out可见Gaussian完全没开始计算,8成是运行环境或者Default.Rou配置有问题
先确保能手动用Gaussian成功运行gau.gjf,再说molclus调用的事
北京科音自然科学研究中心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

14

帖子

0

威望

151

eV
积分
165

Level 3 能力者

338#
发表于 Post on 2025-4-24 08:44:14 | 只看该作者 Only view this author
sobereva 发表于 2025-4-23 15:09
由gau.out可见Gaussian完全没开始计算,8成是运行环境或者Default.Rou配置有问题
先确保能手动用Gaussia ...

谢谢老师

14

帖子

0

威望

151

eV
积分
165

Level 3 能力者

339#
发表于 Post on 2025-4-24 09:22:30 | 只看该作者 Only view this author
sobereva 发表于 2025-4-23 15:09
由gau.out可见Gaussian完全没开始计算,8成是运行环境或者Default.Rou配置有问题
先确保能手动用Gaussia ...

卢老师,根据您的建议,我在配置G16时完全按照《Gaussian的安装方法及运行时的相关问题》来安装:
安装版本为G16-A03-AVX2.tbz
cpu配置如下:

flags显示是支持avx2的
在手动运行时显示如下错误。
Error: illegal instruction, illegal opcode
   rax 00000000014bd760, rbx 00007ffe0437bfb0, rcx 00007ffe0437bfb0
   rdx 00007ffe0437bf80, rsp 00007ffe0437bed8, rbp 00007ffe0437bee0
   rsi 00007ffe0437bf88, rdi 00007ffe0437bf18, r8  0000000000000080
   r9  00000000014a3d00, r10 0000000000000032, r11 0000000000000035
   r12 0000000002cf5ad0, r13 0000144f18f47250, r14 0000144f18f10778
   r15 0000144f18f05b98
  /lib/x86_64-linux-gnu/libc.so.6(+0x42520) [0x145e0dc0d520]
  /root/autodl-tmp/g16/l101.exe() [0x14bd80a]
Aborted (core dumped)





gau.gjf

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

gau.out

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

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

340#
 楼主 Author| 发表于 Post on 2025-4-24 18:28:31 | 只看该作者 Only view this author
caoyx0612 发表于 2025-4-24 09:22
卢老师,根据您的建议,我在配置G16时完全按照《Gaussian的安装方法及运行时的相关问题》来安装:
安装 ...

和文中说的“运行时提示Error: illegal instruction , illegal opcode”情况一样处理
北京科音自然科学研究中心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

8

帖子

0

威望

213

eV
积分
221

Level 3 能力者

341#
发表于 Post on 2025-8-18 15:41:46 | 只看该作者 Only view this author
各位老师好,我用该方法对一个66原子的分子进行构象搜索时,遇到了两个问题
1.用Molclus调用xtb在GFN2-xTB结合溶剂模型下做批量优化这一步,由于GBSA模型中包含的溶剂种类有限,我实验所使用的溶剂不在其中,请问我这一步不加GBSA溶剂模型可以吗?也正是因为我在计算时尝试不加GBSA溶剂模型,所以这一步保留相对能量<=8kcal/mol的结构,不知道这样能否减轻不带溶剂模型带来的误差影响
2.我在用gaussian进行优化和振动分析任务前,保留了能量最低的114个结构,但是在计算到第85个时出现了震荡,该任务运行到步数上限就停止了,并且也没有继续算剩余的结构。请问有没有办法在进行多个优化和振动分析任务时,尽管中间对某一个结构的任务因为SCF不收敛而停止,但是仍然会继续算下一个,能够最后统一处理没能收敛的结构呢。

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

342#
 楼主 Author| 发表于 Post on 2025-8-19 07:29:48 | 只看该作者 Only view this author
谢却人间事 发表于 2025-8-18 15:41
各位老师好,我用该方法对一个66原子的分子进行构象搜索时,遇到了两个问题
1.用Molclus调用xtb在GFN2-xTB ...

1 找一个介电常数与你实际用的溶剂接近的,也比完全不考虑溶剂效应强

2 Molclus调用量子化学程序批量计算的过程中本来就不会因为某个任务失败而导致整个molclus任务停止。如果就是molclus调用的还这样,把molclus输出的末尾信息截图说明
北京科音自然科学研究中心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

8

帖子

0

威望

213

eV
积分
221

Level 3 能力者

343#
发表于 Post on 2025-8-19 10:07:23 | 只看该作者 Only view this author
sobereva 发表于 2025-8-19 07:29
1 找一个介电常数与你实际用的溶剂接近的,也比完全不考虑溶剂效应强

2 Molclus调用量子化学程序批量 ...

明白了,谢谢卢老师,我会加上介电常数接近的溶剂进行计算的。
由于我的molclus调用界面没能截图保存就被关掉了,之后遇到问题一定会及时截图保存,只能从gaussian进行优化和振动分析任务后的out文件中的
Optimization stopped.
    -- Number of steps exceeded,  NStep=  373
    -- Flag reset to prevent archiving.
和Error termination request processed by link 9999.
Error termination via Lnk1e in d:\study\G09W\l9999.exe at Wed Aug 06 18:33:37 2025.
判断是几何收敛失败,guassview里看出震荡。因为之前遇到震荡我都及时手动停止了任务,还没有让程序一直算超然后自动停止的经历。因为我在输出文件里没有设置步数上限,我觉得这个默认的373这个上限有点太高了,所以想再请教一下您,我能否在输入文件里这样写:opt=(calcfc,maxcyc=N)或者opt=(recalc=n,maxcyc=N)在使用更好的Hessian矩阵的同时限制一下步数上限。

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

344#
 楼主 Author| 发表于 Post on 2025-8-20 01:09:33 | 只看该作者 Only view this author
谢却人间事 发表于 2025-8-19 10:07
明白了,谢谢卢老师,我会加上介电常数接近的溶剂进行计算的。
由于我的molclus调用界面没能截图保存就 ...

应该给gjf里设置优化步数上限,比如80-120,要不然震荡不收敛的话会一直傻算着白浪费时间。并且注意下面贴子里说的template2.gjf的用法,可以用于当template.gjf由于震荡导致难收敛、达到预设的步数上限时,切换到帮助收敛的关键词上
使用molclus程序做团簇构型搜索和分子构象搜索
http://bbs.keinsci.com/thread-577-1-1.html

北京科音自然科学研究中心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

8

帖子

0

威望

213

eV
积分
221

Level 3 能力者

345#
发表于 Post on 2025-8-20 09:48:52 | 只看该作者 Only view this author
sobereva 发表于 2025-8-20 01:09
应该给gjf里设置优化步数上限,比如80-120,要不然震荡不收敛的话会一直傻算着白浪费时间。并且注意下面 ...

明白了,谢谢卢老师!

本版积分规则 Credits rule

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

GMT+8, 2026-4-15 07:09 , Processed in 0.360704 second(s), 22 queries , Gzip On.

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