计算化学公社

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

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

  [复制链接 Copy URL]

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

106#
 楼主 Author| 发表于 Post on 2020-11-11 22:46:46 | 只看该作者 Only view this author
CYQ@ 发表于 2020-11-11 10:27
社长,我想请教一下,我用xtb跑动力学然后优化的构相,在用gaussian优化和算单点能之前,能够通过vmd自己 ...

如果你相信自己的直觉,在molclus的settings.ini里设置ngeom,直接设置想要考虑的帧的序号就完了,非常灵活
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

42

帖子

0

威望

530

eV
积分
572

Level 4 (黑子)

107#
发表于 Post on 2020-11-12 09:39:45 | 只看该作者 Only view this author
sobereva 发表于 2020-11-11 22:46
如果你相信自己的直觉,在molclus的settings.ini里设置ngeom,直接设置想要考虑的帧的序号就完了,非常灵 ...

好的,谢谢sob老师

42

帖子

0

威望

530

eV
积分
572

Level 4 (黑子)

108#
发表于 Post on 2020-11-14 10:14:19 | 只看该作者 Only view this author
sobereva 发表于 2020-11-11 22:46
如果你相信自己的直觉,在molclus的settings.ini里设置ngeom,直接设置想要考虑的帧的序号就完了,非常灵 ...

sob老师,我想请教一下,因为我下载的是xtb6.3.3版本。那它调用gfn2优化时加入的隐式水溶剂命令,由--gbsa h2o是要换成--alpb water吗?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

109#
 楼主 Author| 发表于 Post on 2020-11-14 23:42:49 | 只看该作者 Only view this author
CYQ@ 发表于 2020-11-14 10:14
sob老师,我想请教一下,因为我下载的是xtb6.3.3版本。那它调用gfn2优化时加入的隐式水溶剂命令,由--gbs ...

我对alpb持保留意见。虽然原理上比GBSA更好,但xtb里面的实现效果如何、有没有bug我不清楚。最好Grimme发表一篇对于GFN-xTB用的溶剂模型的专门测试文章讨论这个事情
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

12

帖子

0

威望

485

eV
积分
497

Level 3 能力者

110#
发表于 Post on 2020-12-12 17:07:56 | 只看该作者 Only view this author
请教社长:xtb成功安装后,如何再调用crest模块呢?
解压crest.tar后,得到了一个名为crest的文件,执行example-0报错“Cannot find crest binary. Exit.”
不知在.bashrc中添加crest需要输入什么命令。

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

111#
 楼主 Author| 发表于 Post on 2020-12-13 04:39:50 | 只看该作者 Only view this author
CCmaters 发表于 2020-12-12 17:07
请教社长:xtb成功安装后,如何再调用crest模块呢?
解压crest.tar后,得到了一个名为crest的文件,执行ex ...

我不知道你说的example-0指什么
该说的帖子里都说了,根本不需要安装,直接写crest可执行文件路径就能跑
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

25

帖子

0

威望

638

eV
积分
663

Level 4 (黑子)

112#
发表于 Post on 2020-12-15 14:17:41 | 只看该作者 Only view this author
本帖最后由 compXB 于 2020-12-15 15:20 编辑

社长你好,用xtb进行md时用参数 --gfn 0会有如下报错,换成1或2就不会有报错

########################################################################
[ERROR] Program stopped due to fatal error
-3- Could not setup parameterisation
-2- main_setup_newCalculator: Could not construct new calculator
-1- main_setup_newXTBCalculator: Parameter file param_gfn0-xtb.txt not found!
########################################################################
abnormal termination of xtb


这个问题已在另一个帖子中解决了:
[xtb] 求助xtb的安装和使用
http://bbs.keinsci.com/thread-19068-1-1.html

25

帖子

0

威望

638

eV
积分
663

Level 4 (黑子)

113#
发表于 Post on 2020-12-15 15:25:45 | 只看该作者 Only view this author
请教社长:没有gaussian,第四步中涉及到gaussian的怎么用orca代替?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

114#
 楼主 Author| 发表于 Post on 2020-12-16 01:52:29 | 只看该作者 Only view this author
compXB 发表于 2020-12-15 15:25
请教社长:没有gaussian,第四步中涉及到gaussian的怎么用orca代替?

看此文了解ORCA基本使用常识,知道怎么写基本的关键词,自然就知道了
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490http://bbs.keinsci.com/thread-13560-1-1.html
基于ORCA量子化学程序对分子做优化、振动分析、观看红外光谱、观看轨道的简单演示
https://www.bilibili.com/video/av59599938
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

25

帖子

0

威望

638

eV
积分
663

Level 4 (黑子)

115#
发表于 Post on 2020-12-16 10:31:26 | 只看该作者 Only view this author
sobereva 发表于 2020-12-16 01:52
看此文了解ORCA基本使用常识,知道怎么写基本的关键词,自然就知道了
详谈Multiwfn产生ORCA量子化学程序 ...

谢谢社长,我去研读一下这两篇教程

273

帖子

0

威望

4345

eV
积分
4618

Level 6 (一方通行)

116#
发表于 Post on 2020-12-24 22:15:07 | 只看该作者 Only view this author
本帖最后由 函数与激情 于 2020-12-24 22:18 编辑

Sob老师您好,在此请教一个xtb跑BOMD的关键词,在本文中采取了shake算法约束了与氢相关的化学键,且shake约束是默认开启的。之前没有跑过MD,对此关键词不太了解,对于我的体系(近100原子过渡金属体系)在BOMD过程中关闭shake约束结果更好,我看您的另外一篇orca相关的AIMD贴中也没有涉及到shake算法,所以我想知道在跑MD中关闭此关键词是否可以呢?另外,开启或关闭shake约束哪种情况更加接近真实的环境呢?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

117#
 楼主 Author| 发表于 Post on 2020-12-25 04:52:31 | 只看该作者 Only view this author
函数与激情 发表于 2020-12-24 22:15
Sob老师您好,在此请教一个xtb跑BOMD的关键词,在本文中采取了shake算法约束了与氢相关的化学键,且shake约 ...

看你用的步长。常温下模拟,通常>=2 fs才需要开SHAKE避免氢的振动导致动力学不稳定
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

273

帖子

0

威望

4345

eV
积分
4618

Level 6 (一方通行)

118#
发表于 Post on 2020-12-25 07:44:02 | 只看该作者 Only view this author
sobereva 发表于 2020-12-25 04:52
看你用的步长。常温下模拟,通常>=2 fs才需要开SHAKE避免氢的振动导致动力学不稳定

我用的步长是1fs,十分感谢sob老师。

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

119#
 楼主 Author| 发表于 Post on 2020-12-26 05:33:20 | 只看该作者 Only view this author
函数与激情 发表于 2020-12-25 07:44
我用的步长是1fs,十分感谢sob老师。

完全没必要SHAKE,除非是很高温时原子运动很快
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

23

帖子

0

威望

276

eV
积分
299

Level 3 能力者

120#
发表于 Post on 2021-1-13 18:43:29 | 只看该作者 Only view this author
卢老师,Molclus文件夹里的settings.ini设定了freeze冻结的原子序数;在后续计算过程中采用xtb计算,并且结合crest,那么在优化过程中freeze设定的原子是否是冻结了呢

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 10:03 , Processed in 0.305588 second(s), 21 queries , Gzip On.

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