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

计算化学公社

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

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

  [复制链接]

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-1-15 08:03:12 | 显示全部楼层
tianyongpan 发表于 2021-1-13 18:43
卢老师,Molclus文件夹里的settings.ini设定了freeze冻结的原子序数;在后续计算过程中采用xtb计算,并且结 ...
直接用molclus调用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!

17

帖子

0

威望

206

eV
积分
223

Level 3 能力者

发表于 2021-1-15 10:08:54 | 显示全部楼层
sobereva 发表于 2021-1-15 08:03
直接用molclus调用xtb才起效

也就是说./crest -mdopt traj.xyz -gfn0 -opt -chrg 1 -uhf 0,这种调用方式是没有办法固定某些原子不动的,这么理解可对?

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-1-16 01:30:30 | 显示全部楼层
tianyongpan 发表于 2021-1-15 10:08
也就是说./crest -mdopt traj.xyz -gfn0 -opt -chrg 1 -uhf 0,这种调用方式是没有办法固定某些原子不动 ...

显然不起效,根本没有任何选项和设置在这里体现了要冻结原子
北京科音自然科学研究中心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!

17

帖子

0

威望

206

eV
积分
223

Level 3 能力者

发表于 2021-1-24 17:42:21 | 显示全部楼层
sobereva 发表于 2021-1-16 01:30
显然不起效,根本没有任何选项和设置在这里体现了要冻结原子

./crest -mdopt traj.xyz md.inp -gfn0 -opt -chrg 1 -uhf 0 这个命令可否实现的固定序号为1的原子。md.inp的内容是$fix atoms:1 force constant=0 $end

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-1-26 19:27:41 | 显示全部楼层
tianyongpan 发表于 2021-1-24 17:42
./crest -mdopt traj.xyz md.inp -gfn0 -opt -chrg 1 -uhf 0 这个命令可否实现的固定序号为1的原子。md.i ...

xtb目前的MD不支持冻结
北京科音自然科学研究中心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!

1

帖子

0

威望

33

eV
积分
34

Level 2 能力者

发表于 2021-1-26 20:35:08 | 显示全部楼层
谢谢分享

77

帖子

2

威望

984

eV
积分
1101

Level 4 (黑子)

发表于 2021-3-18 12:43:18 | 显示全部楼层
如果我想固定某些原子来跑动力学, 不知道 xtb 能否做到? (因为想对过渡态的构型取样, 但是又不是单纯的键转动, 可能有环反转等等) 或者最简单的方法是什么?

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-3-19 18:24:43 | 显示全部楼层
elvisng 发表于 2021-3-18 12:43
如果我想固定某些原子来跑动力学, 不知道 xtb 能否做到? (因为想对过渡态的构型取样, 但是又不是单纯的键转 ...

xtb的MD不支持固定原子(至少对6.3版而言。6.4没测试过)
北京科音自然科学研究中心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!

13

帖子

0

威望

77

eV
积分
90

Level 2 能力者

发表于 2021-3-29 09:09:24 | 显示全部楼层
sob老师您好  我看到第三步是用Molclus调用xtb在GFN2-xTB结合隐式水模型下做批量优化    请问如果我是做气相状态下的话   可不可以不加溶剂   直接用这个计算crest -mdopt traj.xyz -gfn2 -opt normal -niceprint

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-3-29 18:05:57 | 显示全部楼层
wangchunchun 发表于 2021-3-29 09:09
sob老师您好  我看到第三步是用Molclus调用xtb在GFN2-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!

14

帖子

0

威望

445

eV
积分
459

Level 3 能力者

发表于 2021-4-1 16:02:41 | 显示全部楼层
老师,提取前3低能量结构去Shermo里计算后得到的分布比例与文章给出的isostat里排序后得到的构象分布不同;做这样的构象搜索后,若要增加平均光谱的绘制,选择哪一种分布比例合适些?

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-4-1 16:19:01 | 显示全部楼层
jingsi 发表于 2021-4-1 16:02
老师,提取前3低能量结构去Shermo里计算后得到的分布比例与文章给出的isostat里排序后得到的构象分布不同; ...

如果当前isostat是基于电子能量算的分布概率,明显不如通过Shermo算出来的自由能得到的波尔兹曼分布概率合理。用电子能量算的分布概率只是一个很糙的近似,只是便于筛选过程而已。
北京科音自然科学研究中心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!

14

帖子

0

威望

445

eV
积分
459

Level 3 能力者

发表于 2021-4-1 20:02:19 | 显示全部楼层
sobereva 发表于 2021-4-1 16:19
如果当前isostat是基于电子能量算的分布概率,明显不如通过Shermo算出来的自由能得到的波尔兹曼分布概率 ...

谢谢,这个isostat的分布比例是本篇博文里第四步高精度自由能计算后的

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2021-4-2 03:30:13 | 显示全部楼层
jingsi 发表于 2021-4-1 20:02
谢谢,这个isostat的分布比例是本篇博文里第四步高精度自由能计算后的

按照第四节isostat给出的也能用,不过靠Shermo可以用更好的模型(quasi-RRHO),而且溶解自由能最佳的计算级别是M05-2X/6-31G*,这些因素都会带来差异
北京科音自然科学研究中心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!

15

帖子

0

威望

77

eV
积分
92

Level 2 能力者

发表于 2021-4-2 10:21:12 | 显示全部楼层
我用win10 wsl2计算,按照步骤算到第四步的时候,总内存32GB,如此设置6核4GB内存,但还是提示内存不足
! PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX grid4 gridx4 tightSCF noautostart miniprint nopop
%maxcore 4000
%pal nprocs 6 end
报错
[file orca_tools/Tool-CI/qccitrafo2.cpp, line 10958, Process 2]: Error (ORCA/RI-TRAFO): Not enough memory in RI-Trafo (need 2877.82568359375 MB)

[file orca_tools/Tool-CI/qccitrafo2.cpp, line 10958, Process 3]: Error (ORCA/RI-TRAFO): Not enough memory in RI-Trafo (need 2877.82568359375 MB)

[file orca_tools/Tool-CI/qccitrafo2.cpp, line 10958, Process 4]: Error (ORCA/RI-TRAFO): Not enough memory in RI-Trafo (need 2877.82568359375 MB)

[file orca_tools/Tool-CI/qccitrafo2.cpp, line 10958, Process 5]: Error (ORCA/RI-TRAFO): Not enough memory in RI-Trafo (need 2877.82568359375 MB)

[file orca_tools/Tool-CI/qccitrafo2.cpp, line 10958, Process 0]: Error (ORCA/RI-TRAFO): Not enough memory in RI-Trafo (need 2877.82568359375 MB)

[file orca_tools/Tool-CI/qccitrafo2.cpp, line 10958, Process 1]: Error (ORCA/RI-TRAFO): Not enough memory in RI-Trafo (need 2877.82568359375 MB)

难道说是wsl2虚拟的linux实际内存仍不够么?是不是wsl2总内存加上win10的内存直接占满了,虽然分配了4GB,但是实际已经满了,需要减小核数或增大物理内存
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2021-4-11 17:40 , Processed in 0.190148 second(s), 23 queries .

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