计算化学公社

 找回密码 Forget password
 注册 Register
Views: 20522|回复 Reply: 20
打印 Print 上一主题 Last thread 下一主题 Next thread

[Molclus] molclus已支持Grimme的xtb程序

[复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

跳转到指定楼层 Go to specific reply
楼主
molclus程序(介绍见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=577,主页为http://www.keinsci.com/research/molclus.html)刚发布了1.4版,从此版本开始,支持了Grimme开发的xtb程序。这个程序主要是做GFB-xtb计算的,包括单点、优化、振动分析、动力学等。GFB-xtb具体介绍看J. Chem. Theory Comput., 2017, 13 (5), pp 1989–2009,简单介绍可以看《盘点Grimme迄今对理论化学的贡献》(http://sobereva.com/388)文中第10条。xtb预编译版程序可以向https://www.chemie.uni-bonn.de/pctc/mulliken-center/software/xtb/xtb里面提到的Grimme的邮箱发邮件免费索取,回复很快。

xtb使用很简单,先解压到比如/sob/xtb,然后在用户目录下的.bashrc中添加export XTBHOME=/sob/xtb,最好也添加export PATH=$PATH:/sob/xtb便于调用。之后重新进入终端即可使用。比如xtb test.xyz -sp就是对test.xyz里的坐标做GFN-xtb单点计算,把-sp改成-opt就是做优化计算(产生的xtbopt.coord是最后帧的xyz文件,xtbopt.log是含有优化过程每一帧的多帧xyz文件),等等。直接输入xtb可以看到所有运行选项,xtb目录下的.xtbrc文件含有各种任务用的细节参数。还有一大堆.开头的文件都是理论方法的参数文件。可以设置xtb并行计算核数和OpenMP并行的堆栈内存量,例:
export OMP_NUM_THREADS=4
export MKL_NUM_THREADS=4
export OMP_STACKSIZE=1000m

简单来说,在机子里按照如上方式把xtb的环境变量设好后,就可以在molclus的settings.ini中把iprog设为4,然后在molclus运行时调用xtb对traj.xyz里的结构进行单点计算或优化了。molclus结合xtb使用的时候不需要像调用Gaussian/ORCA/MOPAC那样写模板文件。

2019-Mar-20注:本文说的xtb安装方式已经过时。最新版本xtb的安装基和本使用看
将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421http://bbs.keinsci.com/thread-10106-1-1.html

评分 Rate

参与人数
Participants 4
eV +17 收起 理由
Reason
ggdh + 5 牛!
ccccn345 + 2
smutao + 5 谢谢
muxijiao + 5 好物!

查看全部评分 View all ratings

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

292

帖子

8

威望

1696

eV
积分
2148

Level 5 (御坂)

2#
发表于 Post on 2017-12-1 06:22:11 | 只看该作者 Only view this author
请问sob老师 这个xtb程序可以计算梯度和Hessian吗? 计算后能不能输出这两样信息?谢谢

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

3#
发表于 Post on 2017-12-1 06:48:40 | 只看该作者 Only view this author
请问老师,xtb程序可以做conformational search吗?

99

帖子

0

威望

1388

eV
积分
1487

Level 4 (黑子)

PIG

4#
发表于 Post on 2017-12-1 09:40:20 | 只看该作者 Only view this author
ccccn345 发表于 2017-12-1 06:48
请问老师,xtb程序可以做conformational search吗?

可以配合molclus做

99

帖子

0

威望

1388

eV
积分
1487

Level 4 (黑子)

PIG

5#
发表于 Post on 2017-12-1 09:58:12 | 只看该作者 Only view this author
昨晚mopac跑几个小时的东西xTB竟然十几秒就做完了,震惊震惊

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

6#
 楼主 Author| 发表于 Post on 2017-12-1 12:51:28 | 只看该作者 Only view this author
smutao 发表于 2017-12-1 06:22
请问sob老师 这个xtb程序可以计算梯度和Hessian吗? 计算后能不能输出这两样信息?谢谢

xtb带着-grad可以给出梯度,-hess可以给出Hessian
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

7#
 楼主 Author| 发表于 Post on 2017-12-1 12:53:50 | 只看该作者 Only view this author
muxijiao 发表于 2017-12-1 09:58
昨晚mopac跑几个小时的东西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

99

帖子

0

威望

1388

eV
积分
1487

Level 4 (黑子)

PIG

8#
发表于 Post on 2017-12-1 21:25:32 | 只看该作者 Only view this author
本帖最后由 muxijiao 于 2017-12-1 21:27 编辑
sobereva 发表于 2017-12-1 12:53
结果咋样?靠谱不?

我在搜索二硫化钼的团簇,感觉要比PM7和PM6-D3H4靠谱,这两很多构象都飞了我在想xTB是不是对于过渡族金属化合物要好于PM系列呢?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

9#
 楼主 Author| 发表于 Post on 2017-12-1 22:53:28 | 只看该作者 Only view this author
muxijiao 发表于 2017-12-1 21:25
我在搜索二硫化钼的团簇,感觉要比PM7和PM6-D3H4靠谱,这两很多构象都飞了我在想xTB是不是对于过渡族金属 ...

应该要好于。GFN-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

99

帖子

0

威望

1388

eV
积分
1487

Level 4 (黑子)

PIG

10#
发表于 Post on 2017-12-2 09:36:25 | 只看该作者 Only view this author
sobereva 发表于 2017-12-1 22:53
应该要好于。GFN-xtb普适性较强

export OMP_NUM_THREADS=20
export MKL_NUM_THREADS=20
export OMP_STACKSIZE=10000m
社长,并行设置的bashrc应该你写的少了等号了吧

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

11#
 楼主 Author| 发表于 Post on 2017-12-2 14:42:09 | 只看该作者 Only view this author
muxijiao 发表于 2017-12-2 09:36
export OMP_NUM_THREADS=20
export MKL_NUM_THREADS=20
export OMP_STACKSIZE=10000m

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

99

帖子

0

威望

1388

eV
积分
1487

Level 4 (黑子)

PIG

12#
发表于 Post on 2017-12-3 10:15:56 | 只看该作者 Only view this author

社长,突然发现GFN-xTB的HOWTO文件中发现了这么一段
  1. A better full QM alternative for conformational search is
  2. xtb coord -siman
  3. which runs MD/simulated annealing based conformational search (including SHAKE speeds up these
  4. calculations considerably by allowing the timestep to be increased to 3-4 fs).
  5. set mdtime, mddump, mdstep, shake, as well as nsiman and tend properly e.g.
  6. mdtime     50.0       # MD run time in ps
  7. mddump    1000.       # dump (=optimize) structure in siman every mddump fs      
  8. mdstep      5.0       # MD time step in fs (automatically determined if < 0), could be 7-8 fs with shake =2, md_hmass
  9. md_hmass      4       # increase hydrogen mass to this value in amu (at const. tot. mass) allowing large time steps (=0 off)
  10. shake         2       # shake on (=0: off which is default) for X-H bonds only (=1), or all bonds (=2)
  11. nsiman        3       # number of siman annealing blocks
  12. tend       1000       # highest siman annealing temperature (very system specific)
复制代码

是不是意味着可以代替gmx周期退火来产生初始构型呢?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

13#
 楼主 Author| 发表于 Post on 2017-12-3 12:39:33 | 只看该作者 Only view this author
muxijiao 发表于 2017-12-3 10:15
社长,突然发现GFN-xTB的HOWTO文件中发现了这么一段

是不是意味着可以代替gmx周期退火来产生初始构型 ...

没法代替。gmx是基于力场做的,GFN-xTB是基于半经验DFT做的,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

99

帖子

0

威望

1388

eV
积分
1487

Level 4 (黑子)

PIG

14#
发表于 Post on 2017-12-3 15:43:40 | 只看该作者 Only view this author
sobereva 发表于 2017-12-3 12:39
没法代替。gmx是基于力场做的,GFN-xTB是基于半经验DFT做的,xtb再怎么快也没法代替力场的价值

嗯嗯,速度上确实是的,但是精度上是不是要比力场准确呢?而且在易用程度上貌似要方便一些?对于过度金属团簇,gmx的力场是不是也不太好?(纯粹个人想法,对分子动力学不太熟悉)

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

15#
 楼主 Author| 发表于 Post on 2017-12-3 15:45:40 | 只看该作者 Only view this author
muxijiao 发表于 2017-12-3 15:43
嗯嗯,速度上确实是的,但是精度上是不是要比力场准确呢?而且在易用程度上貌似要方便一些?对于过度金属 ...

比一般力场准确
gmx不适合过渡金属团簇
北京科音自然科学研究中心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, 2025-8-16 18:51 , Processed in 0.169465 second(s), 21 queries , Gzip On.

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