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

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

计算化学公社

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

[Molclus] 使用molclus程序做团簇构型搜索和分子构象搜索

  [复制链接 Copy URL]

68

帖子

0

威望

567

eV
积分
635

Level 4 (黑子)

466#
发表于 Post on 2025-4-2 17:46:08 | 只看该作者 Only view this author
sobereva 发表于 2025-4-2 17:33
1 无所谓,不是重点
2 实测试试,如果结构很不稳定,甚至酸根都跑飞了,也可以用更低温
3 不用。这跟化 ...

感谢卢老师耐心的回复,我现在就去跑一下试试

14

帖子

0

威望

151

eV
积分
165

Level 3 能力者

467#
发表于 Post on 2025-4-2 20:35:41 | 只看该作者 Only view this author
sobereva 发表于 2025-4-2 03:41
先手动用ORCA运行,确保ORCA能正常使用再说molclus调用的事
参看下文
量子化学程序ORCA的安装方法

谢谢老师,手动运行orca时,测试文件仍然出现错误:
--------------------------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 0 on node autodl-container-0c4b4cae70-aa64ddb9 exited on signal 7 (Bus error).
--------------------------------------------------------------------------
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run
我安装的是orca_6_0_1_linux_x86-64_shared_openmpi416.tar.xz,并且按照上面的说明安装了OpenMPI,版本是mpiexec (OpenRTE) 4.1.6

test.out前面看起来都正常,最后三行显示:
ORCA finished by error termination in Startup
Calling Command: mpirun -np 4  /root/autodl-tmp/orca_6_0_1/orca_startup_mpi test.int.tmp test
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

test.out

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

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

468#
 楼主 Author| 发表于 Post on 2025-4-2 20:49:11 | 只看该作者 Only view this author
caoyx0612 发表于 2025-4-2 20:35
谢谢老师,手动运行orca时,测试文件仍然出现错误:
------------------------------------------------ ...

确保运行ORCA时输入了完整的路径。博文里说了这点
北京科音自然科学研究中心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

4

帖子

0

威望

99

eV
积分
103

Level 2 能力者

469#
发表于 Post on 2025-4-3 17:34:51 | 只看该作者 Only view this author
sob老师太牛啦~

14

帖子

0

威望

151

eV
积分
165

Level 3 能力者

470#
发表于 Post on 2025-4-5 23:09:23 | 只看该作者 Only view this author
sobereva 发表于 2025-4-2 20:49
确保运行ORCA时输入了完整的路径。博文里说了这点

谢谢sob老师,我是按照《量子化学程序ORCA的安装方法》使用完整路径:
(base) root@autodl-container-0c4b4cae70-aa64ddb9:~/autodl-tmp# which orca
/root/autodl-tmp/orca_6_0_1/orca
单独测试时出现下面的输出:
(base) root@autodl-container-0c4b4cae70-aa64ddb9:~/autodl-tmp# /root/autodl-tmp/orca_6_0_1/orca test.inp > test.out
--------------------------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node autodl-container-0c4b4cae70-aa64ddb9 exited on signal 7 (Bus error).
--------------------------------------------------------------------------
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

test (1).out

15.81 KB, 下载次数 Times of downloads: 0

68

帖子

0

威望

567

eV
积分
635

Level 4 (黑子)

471#
发表于 Post on 2025-4-8 10:43:36 | 只看该作者 Only view this author
想请教一下老师,我目前在进行溶液中150个原子的有机分子笼和酸根阴离子体系的弱相互作用(C-H...O氢键)分析。我目前需要搜索酸根在有机笼空腔的结合位点,同时我也要搜索有机笼和酸根结合后最稳定的构象,我想通过xtb跑从头算动力学来一次达到这两个目的。我目前学习了您高级量化班和贴子里用xtb对瑞德西韦跑从头算动力学的相关教程,但是我发现那些例子里面都是不带电的中性体系。我看到您的讲义上也没有提到带电体系应该如何设置用xtb跑从头算动力学的输入文件。我想请教老师两个问题

(1)我的有机笼不带电,但是酸根带负电荷。当我用xtb在GFN0-xTB下给这个体系做做动力学模拟的时候,是否要设置体系所带的电荷呢?

(2)我目前计算了150个原子有机笼和酸根形成复合物时候的变形能,发现只有3.5KJ/mol。但是我的有机笼柔性很强,我感觉变形能应该不会这么小。我师兄说是当我把有机笼结构从复合物中抠出来单独优化时,它陷入了一个势能面上的极小点,无法跳出去。我目前准备将有机笼单独用xtb跑一次从头算动力学,找到有机笼单独存在时最稳定的构象,并计算它的单点能,然后再把它和复合物结构中的有机笼计算能量差,从而得到一个精确的变形能,请问您觉得我这样做可以吗?

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

472#
 楼主 Author| 发表于 Post on 2025-4-9 03:49:01 | 只看该作者 Only view this author
peter123 发表于 2025-4-8 10:43
想请教一下老师,我目前在进行溶液中150个原子的有机分子笼和酸根阴离子体系的弱相互作用(C-H...O氢键)分 ...

1 高级班里“DFTB、GFN系列方法与xtb程序的使用”那部分讲xtb使用的部分明确说了--chrg可以设置净电荷。实际情况该怎么设就怎么设,和Gaussian做计算一样

2 可以。注意GFN0-xTB做动力学可以,算变形能就太烂了,最起码也是用GFN2-xTB。另外,对很柔的体系,结构变化大不意味着能量变化也一定大。
北京科音自然科学研究中心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

68

帖子

0

威望

567

eV
积分
635

Level 4 (黑子)

473#
发表于 Post on 2025-4-9 21:40:26 | 只看该作者 Only view this author
本帖最后由 peter123 于 2025-4-9 21:44 编辑
sobereva 发表于 2025-4-9 03:49
1 高级班里“DFTB、GFN系列方法与xtb程序的使用”那部分讲xtb使用的部分明确说了--chrg可以设置净电荷。 ...

非常感谢卢老师耐心的指导,我还想请教一下您,我目前在进行溶液中150个原子的有机分子笼和酸根阴离子体系的弱相互作用(C-H...O氢键)分析。我目前需要搜索酸根在有机笼空腔的结合位点,同时我也要搜索有机笼和酸根结合后最稳定的构象。我参考您用molclus搜索瑞德西韦构象的例子,使用和您那个例子完全一样的计算级别,将温度设为300K(因为担心太高温酸根直接飞出去了),模拟的总时间设为50fs(第一次跑没跑太久),我按照您那个例子,在进行完“用Molclus调用xtb在GFN2-xTB结合隐式水模型做优化并用isostat进行去重和排序”这个步骤后得到了30个非重复结构,
但是如图所示,这30个非重复结构能量差居然非常非常小。
我后续用VMD观看了结构,发现前面17个非重复结构的分子笼核心部分(也就是容纳酸根的空腔)变化不大,酸根的位置和朝向也没有变化,变化大的只涉及到分子笼外侧不参与笼内弱相互作用的3个烯丙基的旋转。
想请教一下老师以下几个问题
1.请问这种情况是因为模拟温度和时间不够导致的吗?
2.请问您觉得我这种情况需要把前面的结构都拿去高斯做优化吗?因为我之前优化这个体系,在B3LYP-D3/6-31G**(酸根的氧加了弥散)级别下优化,都要算一天。而这里涉及到的结构太多了,而我还要研究有机笼和不同种类酸根的作用,我感觉计算量会非常大。
3.请问您觉得这么多个能量非常相近的结构,后续算自由能以及弱相互作用能、弱相互作用结合能应该怎么办呢?
4.请问这种情况是不是对我用波函数分析的手段研究弱相互作用影响并不大呢?因为像我前面提到的,前17个非重复结构涉及到弱相互作用的部分变化都非常小(也就是分子笼容纳酸根的空腔变化不大,酸根的位置和朝向也没有变化)。

202504092140157480..png (142.13 KB, 下载次数 Times of downloads: 37)

202504092140157480..png

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

474#
 楼主 Author| 发表于 Post on 2025-4-10 00:49:32 | 只看该作者 Only view this author
peter123 发表于 2025-4-9 21:40
非常感谢卢老师耐心的指导,我还想请教一下您,我目前在进行溶液中150个原子的有机分子笼和酸根阴离子体 ...

1 是

2 不用。你肉眼都看到了很多结构差异很小。这也体现出可以考虑把去重的阈值放得更宽。并且最好把xtb做几何优化的收敛限设得更严

3 把几何优化收敛限设严、isostat去重阈值放宽后,不会有那么多高度重复的结构。真正有明显差异性的能量较低的结构不会有太多

4 语义不明
北京科音自然科学研究中心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

68

帖子

0

威望

567

eV
积分
635

Level 4 (黑子)

475#
发表于 Post on 2025-4-16 20:06:50 | 只看该作者 Only view this author
sobereva 发表于 2025-4-10 00:49
1 是

2 不用。你肉眼都看到了很多结构差异很小。这也体现出可以考虑把去重的阈值放得更宽。并且最好把 ...

非常感谢卢老师的耐心指导。我还有一些问题想请教一下您,我目前在进行溶液中150个原子的有机分子笼和酸根阴离子体系的弱相互作用(C-H...O氢键)分析。我对有机笼做了振动分析,发现有40个低于100波数的低频,并且全都是有机笼的整体晃动和摆动。因此我感觉酸根和有机笼结合时,有机笼必然会发生整体的变形。我目前需要搜索酸根在有机笼空腔的结合位点,同时我也要搜索有机笼和酸根结合后最稳定的构象,我想通过xtb跑从头算动力学来一次达到这两个目的。我目前学习了您高级量化班和贴子里用xtb跑从头算动力学的相关教程,我参考您那个瑞德西韦的教程来跑我这个结构的动力学。目前有三个问题想请教一下您。


1.您之前指导我说温度设置要通过实测来判断,如果结构很不稳定,或者酸根都跑飞了,就要降低温度。我有点不太明白如何判断结构是否很不稳定?我目前将温度设为350K,发现有机笼的三条侧链都发生了不小的变形和扭曲(我的有机笼的腔体是通过3条较长的侧链构成的),有机笼的腔体也发生了不小的变化,但是有机笼整体的结构并没有散架,没有发生特别大的破坏整体框架的变化。请问您觉得我这种情况是否属于结构很不稳定的情况呢?


2.还想请教一下老师,当我将温度设为350K时,发生跑动力学的时候温度一直350K左右。请问这是否会因为温度较高从而导致整个复合物结构一直处于能量较高的不稳定状态,无法采样到室温下的稳定状态呢?但是我看您那个瑞德西韦的例子里面似乎也没有涉及到降温的操作,因此想请教一下您。


3.请问第一步用GFN0-xTB跑从头算动力学,还有第二步用GFN0-xTB批量优化的时候是不是不用带上溶剂模型?我看到您那个瑞德西韦的例子里面一直到后续用GFN2-xTB进行批量优化的时候才带上溶剂模型。

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

476#
 楼主 Author| 发表于 Post on 2025-4-17 15:11:32 | 只看该作者 Only view this author
peter123 发表于 2025-4-16 20:06
非常感谢卢老师的耐心指导。我还有一些问题想请教一下您,我目前在进行溶液中150个原子的有机分子笼和酸 ...

1 根据你的描述属于正常情况。酸根近乎跑出去或笼的整体框架畸变严重才算是很不稳定。很长的链普遍都有显著柔性,这种温度下变形、扭曲很正常

2 室温下模拟会采样到的结构,高温下虽然采样这些结构的比率占整个轨迹更小,但只要跑的时间够长,终究还是会多多少少采样到

3 如果显明显离子性的区域暴露在外、能和溶剂接触,带上溶剂模型原理上会更好,但也未必对动力学采样这个目的有多少改进,如果实测发现用溶剂模型造成耗时显著增加或者结果异常,那就去掉溶剂模型。

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

68

帖子

0

威望

567

eV
积分
635

Level 4 (黑子)

477#
发表于 Post on 2025-4-17 15:34:45 | 只看该作者 Only view this author
sobereva 发表于 2025-4-17 15:11
1 根据你的描述属于正常情况。酸根近乎跑出去或笼的整体框架畸变严重才算是很不稳定。很长的链普遍都有显 ...

非常感谢卢老师的耐心指导,我还想请教您一下。我之前在本帖466楼请教过您,我当时说我做完GFN2-xTB优化并用isostat去重以后得到的30个非重复结构能量差非常小,您在本帖477楼给我的指导是最好把xtb做几何优化的收敛限设置的更严。我看到您高级量化班里面也提到用xtb优化弱相互作用的复合物时最好把几何优化的收敛限设置成tight。

我现在想请教一个问题,我是参考您瑞德西韦的那个例子来进行构象搜索的,您那个例子先用GFN0-xTB做了一次优化并去重,然后再用GFN2-xTB进一步优化并去重。我想请教一下我能否在GFN2-xTB那一步再去把几何优化的收敛限设置成tight,而GFN0-xTB那一步保持默认收敛限呢?因为收敛限设严可能会导致几何优化的时间变长,我目前跑动力学一共得到了6000帧结构,GFN0-xTB默认收敛限下一分钟只能优化4个结构,我感觉第一步就把收敛限设严格会导致时间变长很多。请问我这样做可以吗?我想听一下您的指导

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

478#
 楼主 Author| 发表于 Post on 2025-4-17 15:47:55 | 只看该作者 Only view this author
peter123 发表于 2025-4-17 15:34
非常感谢卢老师的耐心指导,我还想请教您一下。我之前在本帖466楼请教过您,我当时说我做完GFN2-xTB优化 ...

如果最开始最粗糙的批量优化的收敛限较松,虽然看似当前优化的耗时降低了,由于收敛得不够充分,可能导致去重后留下的构象更多,之后进一步优化时花更多时间。
GFN0-xTB优化也可以改成耗时更低的GFN-FF。我写下文的时候当时xtb还不支持GFN-FF,要不然当时我很可能就用GFN-FF了。GFN-FF描述构象能量差、相互作用能的能力并不比很糙的GFN0-xTB更差
使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索
http://bbs.keinsci.com/thread-16255-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

68

帖子

0

威望

567

eV
积分
635

Level 4 (黑子)

479#
发表于 Post on 2025-4-17 16:31:22 | 只看该作者 Only view this author
sobereva 发表于 2025-4-17 15:47
如果最开始最粗糙的批量优化的收敛限较松,虽然看似当前优化的耗时降低了,由于收敛得不够充分,可能导致 ...

非常感谢卢老师耐心的指导,我刚才去复习了一下高级量化班GFN-FF的相关部分,请问我用GFN-FF代替GFN0-xTB做第一步弱相互作用复合物的几何优化,同样可以将收敛限设置成tight是吧?我看了您的指导,决定还是第一步优化就把收敛限设置成tight

6万

帖子

99

威望

6万

eV
积分
126179

管理员

公社社长

480#
 楼主 Author| 发表于 Post on 2025-4-17 16:35:25 | 只看该作者 Only view this author
peter123 发表于 2025-4-17 16:31
非常感谢卢老师耐心的指导,我刚才去复习了一下高级量化班GFN-FF的相关部分,请问我用GFN-FF代替GFN0-xTB ...

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

本版积分规则 Credits rule

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

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

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