wzkchem5 发表于 2023-8-3 16:42 有道理,谢谢wzk老师! |
sobereva 发表于 2023-8-3 17:37 那还是很有意思的,谢谢社长指点 |
dantevinsky 发表于 2023-8-3 08:51 Molclus的isostat是根据原子间距离序列来去重的,确保不会因为原子顺序问题影响判断。同时结合能量判断阈值,可以保证去重效果的合理。光通过能量在很多情况下还没法充分去重。 |
dantevinsky 发表于 2023-8-3 01:51 RMSD默认不考虑原子的交换对称性,但也存在考虑原子的交换对称性的RMSD算法,比如https://github.com/charnley/rmsd这个脚本就可以。但是考虑原子的交换对称性的RMSD算法计算量比较大,所以可以用能量、偶极矩、转动常数等初步去重,但是最后验证两个结构到底是不是一样,还是得做考虑原子的交换对称性的RMSD计算,只有这样判断才严格,不然总有一定可能有两个不同结构的能量、偶极矩、转动常数碰巧一致。不严格的判据越多,把不同结构误判断为一致的概率就越低,但永远无法降到0 |
dantevinsky 发表于 2023-8-3 08:51 好嘞,谢谢您 |
zjxitcc 发表于 2023-8-2 23:22 好的,谢谢老师! |
sobereva 发表于 2023-8-2 23:09 收到,谢谢老师! |
本帖最后由 zjxitcc 于 2023-8-3 09:49 编辑 dantevinsky 发表于 2023-8-3 08:51 这是因为常规RMSD没考虑同种元素原子可交换的特点,默认两个结构的原子1-1, 2-2, 3-3一一对应。这个主要在一些中小分子上体现为没法完全去重。它只是没法保证100%剔除重复结构,而不是会去剔除有用的结构,因此并非不能用。 MOKIT中已有RMSD功能,很快将会加入原子序号无关的RMSD功能,届时将不存在这种问题。 |
前排提醒,尽量不要采用RMSD去重,因为存在原子的交换对称性,可以采用如:能量、偶极矩等方式去重 |
凡是构象搜索程序肯定有RMSD去重功能的。。。你要不要看看程序输出、程序手册 |
用molclus(http://www.keinsci.com/research/molclus.html)构象搜索,isostat自动就会根据能量和结构相似性去重,不需要自己顾虑 参看下文的例子 ,写得都极其清楚,遇到问题随时可以问我 使用molclus程序做团簇构型搜索和分子构象搜索 http://bbs.keinsci.com/thread-577-1-1.html gentor:扫描方式做分子构象搜索的便捷工具 http://bbs.keinsci.com/thread-2388-1-1.html genmer:生成团簇初始构型结合molclus做团簇结构搜索的超便捷工具 http://bbs.keinsci.com/thread-2369-1-1.html 将Confab或Frog2与Molclus联用对有机体系做构象搜索 http://bbs.keinsci.com/thread-20063-1-1.html 使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索 http://bbs.keinsci.com/thread-16255-1-1.html |
wzkchem5 发表于 2023-8-2 22:42 我是用crest搜索出很多构象的gjf文件,然后优化 |
先说清楚你用的是什么构象搜索软件,构象搜索软件又不唯一 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-26 23:35 , Processed in 0.212734 second(s), 26 queries , Gzip On.