计算化学公社

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

[Molclus] 使用molclus程序做团簇构型搜索和分子构象搜索

  [复制链接 Copy URL]

147

帖子

4

威望

1567

eV
积分
1794

Level 5 (御坂)

科研泥瓦匠

211#
发表于 Post on 2019-9-10 19:29:22 | 只看该作者 Only view this author
sobereva 发表于 2018-8-29 21:48
如果你当前是用的isostat默认的归簇标准,可以把两个归簇标准都设成0.5,然后把前10个取出来用于之后DFT ...

sob老师,两个归簇标准设置有很大的随意性,怎么设置才能计算的合理,时间合适,主要的构象也不会漏掉?比如下面是我的结构,有50个原子,在初步gentor的构象搜索及简单的半经验优化后,准备下一步用B97-3C进行优化,优化时间控制在8小时内。

isomers.xyz

132.51 KB, 下载次数 Times of downloads: 1

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

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

212#
 楼主 Author| 发表于 Post on 2019-9-10 19:34:18 | 只看该作者 Only view this author
欢乐多 发表于 2019-9-10 19:29
sob老师,两个归簇标准设置有很大的随意性,怎么设置才能计算的合理,时间合适,主要的构象也不会漏掉? ...

先用默认的标准。如果你用VMD观看cluster.xyz,发现里头有一些结构你肉眼看上去觉得太接近、优化完了肯定也会收敛到相同结构去,可以尝试更宽的阈值重新归簇
北京科音自然科学研究中心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

147

帖子

4

威望

1567

eV
积分
1794

Level 5 (御坂)

科研泥瓦匠

213#
发表于 Post on 2019-9-19 23:18:32 | 只看该作者 Only view this author
sobereva 发表于 2019-9-10 19:34
先用默认的标准。如果你用VMD观看cluster.xyz,发现里头有一些结构你肉眼看上去觉得太接近、优化完了肯定 ...

老师,在对我的一个分子构象搜索后选出了前十四个构象进行了优化及振动分析,发现10和14有虚频,进行结构微调,虚频消除,那么10与14微调后优化的构象会不会与其他构象重复,使得构象搜索不全面?
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

214#
 楼主 Author| 发表于 Post on 2019-9-21 00:43:59 | 只看该作者 Only view this author
欢乐多 发表于 2019-9-19 23:18
老师,在对我的一个分子构象搜索后选出了前十四个构象进行了优化及振动分析,发现10和14有虚频,进行结构 ...

有可能优化后变成了其它现有的结构。这和全面不全面没有必然关系,只能说有效的样本数从14个变少了点
北京科音自然科学研究中心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

147

帖子

4

威望

1567

eV
积分
1794

Level 5 (御坂)

科研泥瓦匠

215#
发表于 Post on 2019-9-21 09:23:20 | 只看该作者 Only view this author
sobereva 发表于 2019-9-21 00:43
有可能优化后变成了其它现有的结构。这和全面不全面没有必然关系,只能说有效的样本数从14个变少了点

好的,老师
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

147

帖子

4

威望

1567

eV
积分
1794

Level 5 (御坂)

科研泥瓦匠

216#
发表于 Post on 2019-9-22 09:54:32 | 只看该作者 Only view this author
sobereva 发表于 2019-9-21 00:43
有可能优化后变成了其它现有的结构。这和全面不全面没有必然关系,只能说有效的样本数从14个变少了点

老师,博文中提到做振动分析,需要将itask改为1(single point energy), 在设置文件中itask为0:optimization (frequency analysis keyword can also be presented),是不是振动分析itask也可以为0?这个怎么理解?
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

217#
 楼主 Author| 发表于 Post on 2019-9-23 07:53:18 | 只看该作者 Only view this author
欢乐多 发表于 2019-9-22 09:54
老师,博文中提到做振动分析,需要将itask改为1(single point energy), 在设置文件中itask为0:optimiz ...

itask=0的时候可以同时写上振动分析的关键词,此时molclus也不会报错,不过molclus不会去从中提取额外的和振动分析有关的数据,需要用户自己提取

下一个版本可能加上itask=-1,可以专门做振动分析任务,自动判断有无虚频、提取热力学校正量等
北京科音自然科学研究中心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

147

帖子

4

威望

1567

eV
积分
1794

Level 5 (御坂)

科研泥瓦匠

218#
发表于 Post on 2019-9-24 00:32:00 | 只看该作者 Only view this author
sobereva 发表于 2019-9-23 07:53
itask=0的时候可以同时写上振动分析的关键词,此时molclus也不会报错,不过molclus不会去从中提取额外的 ...

好的,老师,我会时刻关注molclus版本的更新
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

70

帖子

0

威望

4882

eV
积分
4952

Level 6 (一方通行)

219#
发表于 Post on 2019-9-25 22:48:20 | 只看该作者 Only view this author
社长,您好,想请教您一个关于构象归簇的问题,我计算的对象主要是天然产物小分子,:
在对一批构象进行优化之后,一般会根据构象两两之间的能量差异和几何差异(基于RMSD或距离矩阵)对构象进行归簇,然后再对归簇后的构象进行下一步计算。
事实上,molclus中的isostat就提供了这样的功能。但问题在于,我一般是将opt+freq一起做的,在此之后用更高级别做单点能计算,与吉布斯自由能的校正能相加获得吉布斯自由能。但如果在opt+freq与单点能计算之间使用isostat,将使得前后的数值(校正能与单点能)无法对应。
事实上,我目前采用的方法是将优化后的构象在一些可视化软件(VMD、Pymol、Macromodel等)中进行align或手动叠加后计算RMSD,但是在构象较多的时候还是非常费时。
综上,想请教您,是否有更快捷的方法,对于一小批优化后的构象,输出构象两两之间的几何差异。或者,可否令isostat在运行时,输出构象两两之间差异的矩阵(一般进行DFT优化的构象不会太多,所以矩阵大小应该可以接受)
非常感谢您!

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

220#
 楼主 Author| 发表于 Post on 2019-9-26 04:20:31 | 只看该作者 Only view this author
happyknighthawk 发表于 2019-9-25 22:48
社长,您好,想请教您一个关于构象归簇的问题,我计算的对象主要是天然产物小分子,:
在对一批构象进行优 ...

你的意思是做完opt freq之后,用isostat筛选出能量较低的结构之后,对这些结构你再做的高精度单点计算?
如果是这样,我不建议opt freq一起做,不如只做opt,待isostat筛选后,把选出的那批结构再分别做freq和单点,这样算完了自由能热校正量和单点能直接就能相加。而且本身freq也是比较耗时的任务,最好对筛选出来的能量比较低的那些做。

你可以用VMD把轨迹保存为pdb格式,用gromacs里的cluster命令处理,这个命令是做归簇的,计算过程中也会把每一帧之间的RMSD值一次性给出来
北京科音自然科学研究中心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

70

帖子

0

威望

4882

eV
积分
4952

Level 6 (一方通行)

221#
发表于 Post on 2019-9-26 09:26:07 | 只看该作者 Only view this author
sobereva 发表于 2019-9-26 04:20
你的意思是做完opt freq之后,用isostat筛选出能量较低的结构之后,对这些结构你再做的高精度单点计算?
...

好的,非常感谢社长的建议,我一会儿就去试试Gromacs的Cluster命令,再次感谢!

70

帖子

0

威望

4882

eV
积分
4952

Level 6 (一方通行)

222#
发表于 Post on 2019-9-26 17:32:48 | 只看该作者 Only view this author
社长,我还有一个关于归簇的两个标准的问题,一般在进行完构象搜索后,会采用半经验方法对构象进行预优化,优化完生成的isomer.xyz会用isostat进行归簇,我的问题是,在这一阶段进行归簇时使用构象间能量差异作为归簇标准的意义有多大?是否可以仅仅按照构象间的几何差异来进行归簇,而完全放弃能量差异(比如区分重复结构的能量阈值直接设置为99)?谢谢社长

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

223#
 楼主 Author| 发表于 Post on 2019-9-27 03:23:43 | 只看该作者 Only view this author
happyknighthawk 发表于 2019-9-26 17:32
社长,我还有一个关于归簇的两个标准的问题,一般在进行完构象搜索后,会采用半经验方法对构象进行预优化, ...

只用结构作为标准应当也可以,设个能量阈值算是增加个双保险,以免一些特殊情况
北京科音自然科学研究中心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

70

帖子

0

威望

4882

eV
积分
4952

Level 6 (一方通行)

224#
发表于 Post on 2019-9-27 10:42:15 | 只看该作者 Only view this author
sobereva 发表于 2019-9-27 03:23
只用结构作为标准应当也可以,设个能量阈值算是增加个双保险,以免一些特殊情况

好的,谢谢社长!

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

225#
 楼主 Author| 发表于 Post on 2019-9-28 13:51:56 | 只看该作者 Only view this author
今日molclus发布了1.8.9版,更新内容如下

• 更新了isostat。改进在于将某isomer纳入已有的簇的时候,如果这个isomer的能量比那个簇的能量更低,则这个簇的结构和能量将会被这个isomer的所代替。另外,从本版本开始的Windows版中,isostat.exe对应64bit版,而isostat_32bit.exe对应32bit版。
• 当ibkout>0而MOPAC运行没有正常结束(没生成arc)的时候,原先版本无法备份文件,而当前版本会将out进行备份。
• gentor里的critskip参数的默认值从原来的0.5改大为了0.8,避免得到的初始结构里有一些原子因为离得太近在优化后变成诡异结构。
• 修正了itask=0的时候与最新版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

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 05:43 , Processed in 0.221562 second(s), 23 queries , Gzip On.

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