计算化学公社

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

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

  [复制链接 Copy URL]

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

46#
发表于 Post on 2020-3-29 21:01:59 | 只看该作者 Only view this author
本帖最后由 mutron 于 2020-3-30 01:47 编辑

编辑掉。。


犯了低级理解错误,需要“当发现文件名里面有crest字样的时候就会自动以crest_ensemble.xyz文件的格式来读取内容”,把含有crest的文件名写完整没问题了

68

帖子

4

威望

1566

eV
积分
1714

Level 5 (御坂)

大鸭鸭

47#
发表于 Post on 2020-4-6 22:16:12 | 只看该作者 Only view this author
想问一下社长用xtb做半经验动力学的时候能不能像Gaussian跑admp一样设置初速度的方向?我看手册上面好像没有这个选项
The leaping arc at your finger
is my faith and will eliminate never
Only the railgun will live forever

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

48#
 楼主 Author| 发表于 Post on 2020-4-7 06:59:11 | 只看该作者 Only view this author
linqiaosong 发表于 2020-4-6 22:16
想问一下社长用xtb做半经验动力学的时候能不能像Gaussian跑admp一样设置初速度的方向?我看手册上面好像没 ...

先跑一下,产生mdrestart文件,然后再篡改mdrestart文件内容,通过restart=true来读入初速度(PS:ORCA的AIMD目前也得用这种做法)

例如

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

32

帖子

0

威望

2884

eV
积分
2916

Level 5 (御坂)

49#
发表于 Post on 2020-6-6 15:41:44 | 只看该作者 Only view this author
请问,第一步,用xtb在GFN0-xTB下做动力学模拟,可不可以同时或者依次对多个化合物的xyz文件进行批量动力学模拟,分别形成多个结果文件?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

50#
 楼主 Author| 发表于 Post on 2020-6-6 18:43:14 | 只看该作者 Only view this author
Evanwill 发表于 2020-6-6 15:41
请问,第一步,用xtb在GFN0-xTB下做动力学模拟,可不可以同时或者依次对多个化合物的xyz文件进行批量动力学 ...

如果有明确目的,可以
北京科音自然科学研究中心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!

194

帖子

0

威望

3694

eV
积分
3888

Level 5 (御坂)

51#
发表于 Post on 2020-6-13 05:23:13 | 只看该作者 Only view this author
表示构象差别的overlay diagram是用什么软件做的?谢谢

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

52#
 楼主 Author| 发表于 Post on 2020-6-13 22:07:30 | 只看该作者 Only view this author
scf 发表于 2020-6-13 05:23
表示构象差别的overlay diagram是用什么软件做的?谢谢

在VMD中计算RMSD衡量两个结构间的差异以及叠合两个结构
http://sobereva.com/290http://bbs.keinsci.com/thread-1080-1-1.html
北京科音自然科学研究中心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!

194

帖子

0

威望

3694

eV
积分
3888

Level 5 (御坂)

53#
发表于 Post on 2020-6-14 04:10:02 | 只看该作者 Only view this author
thanks

141

帖子

4

威望

1406

eV
积分
1627

Level 5 (御坂)

科研泥瓦匠

54#
发表于 Post on 2020-6-25 21:43:15 | 只看该作者 Only view this author
sobereva 发表于 2020-6-13 22:07
在VMD中计算RMSD衡量两个结构间的差异以及叠合两个结构
http://sobereva.com/290(http://bbs.keinsci.c ...

老师,您好,在用XTB对14元大环柔性体系做构象搜索时,初始由chemdraw画的结构,弄成初始结构,进行构象搜索,得到的构象优化后与单晶差别很大;400~1000K中尝试不同温度,还是没有搜到与单晶相似的构象。
因为有多个14元的大环柔性化合物,仅有其中的1个有单晶数据,因此,以有单晶数据结构的化合物为参照,检测XTB做构象搜索方法的可靠性,尝试多种条件,终未搜索到与单晶相似的构象。

请问老师,xtb有没有适合大环柔性体系构象搜索的条件?

TIM图片20200625213007.jpg (94.21 KB, 下载次数 Times of downloads: 96)

TIM图片20200625213007.jpg
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

1187

帖子

5

威望

2841

eV
积分
4129

Level 6 (一方通行)

55#
发表于 Post on 2020-6-25 21:48:29 | 只看该作者 Only view this author
欢乐多 发表于 2020-6-25 21:43
老师,您好,在用XTB对14元大环柔性体系做构象搜索时,初始由chemdraw画的结构,弄成初始结构,进行构象 ...

上面整个流程计算的本身就不是凝聚相的结构,但是单晶是凝聚相。你看看有没有可能是这个原因。

141

帖子

4

威望

1406

eV
积分
1627

Level 5 (御坂)

科研泥瓦匠

56#
发表于 Post on 2020-6-25 23:41:34 | 只看该作者 Only view this author
snljty 发表于 2020-6-25 21:48
上面整个流程计算的本身就不是凝聚相的结构,但是单晶是凝聚相。你看看有没有可能是这个原因。

凝聚相是怎么回事?高斯可以体现凝聚相吗?
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

1632

帖子

4

威望

4088

eV
积分
5800

Level 6 (一方通行)

喵星人

57#
发表于 Post on 2020-6-27 00:36:25 | 只看该作者 Only view this author
欢乐多 发表于 2020-6-25 21:43
老师,您好,在用XTB对14元大环柔性体系做构象搜索时,初始由chemdraw画的结构,弄成初始结构,进行构象 ...

无论用Gaussian还是xtb等优化构型的时候都是游离的分子,即在真空/溶剂中进行的,而在晶体中,分子构象明显会受到相邻其他分子的相互作用,可能包括氢键,pi-pi堆积等等,所以晶体中的构象在游离态中未必是稳定的。搜索不到晶体的构象是很正常的。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
欢乐多 + 5 正解

查看全部评分 View all ratings

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

58#
 楼主 Author| 发表于 Post on 2020-6-27 01:10:52 | 只看该作者 Only view this author
欢乐多 发表于 2020-6-25 21:43
老师,您好,在用XTB对14元大环柔性体系做构象搜索时,初始由chemdraw画的结构,弄成初始结构,进行构象 ...

晶体环境对于这种高度柔性的分子构象影响是相当大的
真空下模拟不可能指望得到和晶体中完全一样的结构

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
欢乐多 + 5 谢谢

查看全部评分 View all ratings

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

141

帖子

4

威望

1406

eV
积分
1627

Level 5 (御坂)

科研泥瓦匠

59#
发表于 Post on 2020-6-27 09:33:44 | 只看该作者 Only view this author
本帖最后由 欢乐多 于 2020-6-27 09:45 编辑
喵星大佬 发表于 2020-6-27 00:36
无论用Gaussian还是xtb等优化构型的时候都是游离的分子,即在真空/溶剂中进行的,而在晶体中,分子构象明 ...

好的,老师,感谢回复!老师,如果用真空下进行团簇搜索,比如放进去2~4个分子,搜到的稳定构象与单晶团簇进行对比,是否有意义?
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

1632

帖子

4

威望

4088

eV
积分
5800

Level 6 (一方通行)

喵星人

60#
发表于 Post on 2020-6-27 11:08:21 | 只看该作者 Only view this author
欢乐多 发表于 2020-6-27 09:33
好的,老师,感谢回复!老师,如果用真空下进行团簇搜索,比如放进去2~4个分子,搜到的稳定构象与单晶团 ...

没有,单晶结构就是只有在无限重复的情况下成立的。只有几个分子连完整的晶胞都构建不出来。
一般来说只有像蛋白,MOF(按团簇处理)之类分子内相互作用很强的,远远强于分子间作用的物种,在溶液中的结构会和晶体基本一致。

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 09:56 , Processed in 0.180838 second(s), 24 queries , Gzip On.

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