请选择 进入手机版 | 继续访问电脑版
第9届北京科音分子动力学与GROMACS培训班将于4月17~20日于北京举办,请点击此链接查看培训详情,欢迎参加和相互转告!

计算化学公社

 找回密码
 现在注册!
楼主: sobereva

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

  [复制链接]

161

帖子

0

威望

546

eV
积分
707

Level 4 (黑子)

蓝卫兵

发表于 2020-2-10 08:38:06 | 显示全部楼层
exity 发表于 2020-2-9 15:29
XTB果然万能,社长有机会能出一个crest的帖子吗?

你也不想想 社长花这么多时间写这个是为了推销Grimme的crest程序吗? 当然不是
你居然还在这边问crest的事情。。。
虽然crest真的很香 亲测。。。
B样条插值
个人专栏https://zhuanlan.zhihu.com/p/21936803

161

帖子

0

威望

546

eV
积分
707

Level 4 (黑子)

蓝卫兵

发表于 2020-2-10 08:42:02 | 显示全部楼层
本帖最后由 pyscf 于 2020-2-10 08:46 编辑
sobereva 发表于 2020-2-9 16:53
文中说了,用crest没意义,所以不可能写帖子。本文的做法就已经完美了。

需要注意的是 你的protocol第一步就用了Grimme的xtb 殊不知crest就是原生基于xtb的
而且所有步骤一步到位(Grimme弄成了全自动的,可能“更完美”)
另外 由于其metadynamics的CV是rmsd 因此就不存在(克服了)一定需要升温至x00度才能“有机会”access到一个构象的问题
B样条插值
个人专栏https://zhuanlan.zhihu.com/p/21936803

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2020-2-11 12:09:48 | 显示全部楼层
pyscf 发表于 2020-2-10 08:42
需要注意的是 你的protocol第一步就用了Grimme的xtb 殊不知crest就是原生基于xtb的
而且所有步骤一步到 ...

文中已然明确写了
名为crest的利用xtb跑metadynamics的构象搜索工具
正因为实际用过,知道效果,所以在对比之下才说没有必要用。

没有什么情况是“一定升温到xxx”才能翻越的,这和模拟时间长度直接相关,想加速采样效率提升到更高温度就完了。为什么没必要用crest在帖子和回复里我都已经说了。

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口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!

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2020-2-11 17:03:26 | 显示全部楼层
更新了本文,主要是在文末加入了借用crest大幅降低xtb批量优化耗时的说明,还增加了DFT与GFN-xTB优化结构的对比等内容。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口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!

302

帖子

2

威望

1156

eV
积分
1498

Level 4 (黑子)

发表于 2020-2-11 17:21:15 | 显示全部楼层
对小分子来说,确定质子化状态还是一个很蛋疼的问题,crest这方面做的还挺方便的感觉。

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2020-2-11 20:12:37 | 显示全部楼层
这里说一下crest的构象搜索的结果(输入结构和本文中高温MD初始结构相同):
crest结合GFN2-xTB真空下搜索,两个半小时结束,给出的能量最低的至少前50个结构里没有任何一个和真实的能量最低结构有相似性
crest结合GFN2-xTB在GBSA表现的水溶剂下搜索,两个小时结束,给出的能量最低的至少前50个结构里没有任何一个和真实的能量最低结构有相似性
crest结合GFN0-xTB真空下做搜索,45分钟结束,所有给出的516个结构没有任何一个结构和真实的能量最低结构有相似性
可见对于当前体系大失败。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口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!

77

帖子

2

威望

984

eV
积分
1101

Level 4 (黑子)

发表于 2020-2-11 22:08:09 | 显示全部楼层
sobereva 发表于 2020-2-11 20:12
这里说一下crest的构象搜索的结果(输入结构和本文中高温MD初始结构相同):
crest结合GFN2-xTB真空下搜索 ...

可否分析一下是什么可能的原因让 crest 的结果这么差吗?

2

帖子

0

威望

167

eV
积分
169

Level 3 能力者

发表于 2020-2-12 09:50:25 | 显示全部楼层
运行第二部优化出现如下错误:
                        *** Configuration     1  ***
Current date: 2020-02-11   Time: 20:40:32
Loading geometry         1 from the inputted geometry file
Generating xtb.xyz file...
Running xtb: xtb xtb.xyz  --opt --gfn 0 --chrg 0 --uhf 0 > xtb.out
normal termination of xtb
Extracting energy and geometry from xtbopt.xyz
forrtl: severe (59): list-directed I/O syntax error, unit 10, file /home/gauss/molclus_1.9.2_Linux/xtbopt.xyz
Image              PC                Routine            Line        Source            
molclus            0000000000419DDB  Unknown               Unknown  Unknown
molclus            0000000000446D0B  Unknown               Unknown  Unknown
molclus            000000000040A36B  Unknown               Unknown  Unknown
molclus            0000000000404422  Unknown               Unknown  Unknown
libc-2.17.so       00007F8EDB51E505  __libc_start_main     Unknown  Unknown
molclus            0000000000404329  Unknown               Unknown  Unknown

689

帖子

2

威望

1615

eV
积分
2344

Level 5 (御坂)

发表于 2020-2-12 09:54:44 | 显示全部楼层
本帖最后由 snljty 于 2020-2-12 09:56 编辑
szzjd 发表于 2020-2-12 09:50
运行第二部优化出现如下错误:
                        *** Configuration     1  ***
Current date:  ...

十有八九是你的xtb版本太老,Molclus版本"太新"。
Molclus的开发日志有:
2020-Jan-19:1.9.1版发布。
• 兼容了2019年12月18日发布的xtb 6.3版预览版,对老版本xtb不再兼容。

到git-hub下载最新版xtb。

2

帖子

0

威望

167

eV
积分
169

Level 3 能力者

发表于 2020-2-12 10:11:43 | 显示全部楼层
snljty 发表于 2020-2-12 09:54
十有八九是你的xtb版本太老,Molclus版本"太新"。
Molclus的开发日志有:
2020-Jan-19:1.9.1版发布。
...

是的,问题解决,谢谢!

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2020-2-12 13:52:25 | 显示全部楼层
elvisng 发表于 2020-2-11 22:08
可否分析一下是什么可能的原因让 crest 的结果这么差吗?

从原理以及实际结果的分析,有两个可能
(1)crest默认用的iMTD-GC方法在其流程的设计上,目的就是去找能量最低构象(也包括利用MD在能量最低构象附近进行采样),因此采样并不全面、充分。然而,GFN-xTB系列方法中哪怕精度最高的GFN2-xTB,其精度还是很有限,其势能面上能量最低构象并非是真实的能量最低构象(而且可能结构相差很大),这就导致crest搜索过程采样区域高度侧重于GFN-xTB势能面的能量最低区域,这反倒往往把真正能量最低结构区域漏过去了(这点在GFN0-xTB上尤为显著,给出来的能量最低构象凭化学直觉就知道其实是高能构象才对)。
(2)crest搜索过程默认的参数对于当前体系可能还不合理,但这个因素相对于(1)来说较为次要。
另外,crest的搜索算法有明显的历史依赖性,这导致最终结果可能对初始结构的选取敏感。

相对来说,crest的搜索方式,对于那种势垒很高(直接做高温MD都不好翻越的),而且GFN-xTB的精度又能正确区分不同构象/构型能量高低次序的时候有用。而对于瑞德西韦这种靠高温MD就足矣充分采样、遍历各种构象区域,而且GFN-xTB的精度又显得明显不足的柔性大体系,用crest不仅没意义,反倒结果可能严重误导人。

评分

参与人数 1eV +5 收起 理由
snljty + 5 GJ!

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口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!

29

帖子

0

威望

502

eV
积分
531

Level 4 (黑子)

发表于 2020-2-12 15:12:45 | 显示全部楼层
szzjd 发表于 2020-2-12 09:50
运行第二部优化出现如下错误:
                        *** Configuration     1  ***
Current date:  ...

前两天我遇到了同样的问题。谢谢你汇报这个错误。

29

帖子

0

威望

502

eV
积分
531

Level 4 (黑子)

发表于 2020-2-12 15:13:22 | 显示全部楼层
snljty 发表于 2020-2-12 09:54
十有八九是你的xtb版本太老,Molclus版本"太新"。
Molclus的开发日志有:
2020-Jan-19:1.9.1版发布。
...

厉害!

77

帖子

2

威望

984

eV
积分
1101

Level 4 (黑子)

发表于 2020-2-12 21:06:17 | 显示全部楼层
本帖最后由 elvisng 于 2020-2-12 21:47 编辑

我尝试用 crest 做批次优化,但几秒就完成,出来的 crest_ensemble.xyz 是空白的。请问有人遇到这问题吗?(换成 gfn1 就没问题, 只是 gfn0 有问题)

以下是輸出內容:



> crest -mdopt traj.xyz -gfn0 -opt normal -niceprint

  -gfn0 : Use of GFN0-xTB requested.
  -opt 0
=============================
  # threads =          36
=============================

------------------------------
GFNn-xTB Geometry Optimization
------------------------------
Geometry successfully optimized.

     =======================================
     |               M D O P T             |
     =======================================
     Optimization along trajectory (or ensemble).

Input file: <traj.xyz>
Containing 2000 structures.
writing TMPCONF* Dirs from file "traj.xyz" ... done.
Starting optimization of generated structures
2000 jobs to do.
[##################################################]  100.00% finished.
done.
Now appending opt.xyz file with new structures

Optimized ensemble on file <crest_ensemble.xyz>

-----------------
Wall Time Summary
-----------------
               MDOPT wall time :         0h : 0m : 3s
--------------------
Overall wall time  : 0h : 0m : 3s

CREST terminated normally.

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2020-2-13 17:36:28 | 显示全部楼层
elvisng 发表于 2020-2-12 21:06
我尝试用 crest 做批次优化,但几秒就完成,出来的 crest_ensemble.xyz 是空白的。请问有人遇到这问题吗?( ...

没遇见过
先看看xtb直接优化此体系能否正常完成
另外,mdopt运行中途会在当前目录下产生临时目录,里面有自动产生的xtb输入输出信息,看看里面内容
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口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!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )

GMT+8, 2021-4-11 18:24 , Processed in 0.281536 second(s), 24 queries .

快速回复 返回顶部 返回列表