计算化学公社

标题: 批量提取柔性扫描的能量及结构 [打印本页]

作者
Author: 123qwertybobo 时间: 2020-9-21 16:34
标题: 批量提取柔性扫描的能量及结构
批量提取柔性扫描的能量及结构，意义还是有的。但二维扫描的能量和坐标提取都比较麻烦。论坛见过有用perl和python的，但不适合我的使用习惯。个人喜欢自带的awk和sed，啥也不用装。基本原理：通过 Stationary point found 定位每个单点柔扫的结果，在此基础上通过Standard orientation:定位相应坐标。

脚本不长，没使用到特别高级的语句，除了为替换坐标符号用到了gsub。如果有些元素符号没有定义，可依脚本画葫芦。
使用tips：
需要一个/多个柔扫的log文件，把脚本放于当前目录，执行即可，不需要其他参数。使用both.sh，会批量提取能量且批量提取坐标。使用energy.sh可以监测当前扫至哪步。
两个脚本均已测试正确，若不幸存在，还望请告知。
希望能对初学者有所裨益，见笑于诸位old birds。

作者
Author: staq001 时间: 2020-10-6 07:25
谢谢分享

作者
Author: 模拟小菜 时间: 2023-5-9 17:30
感谢分享，十分好用

作者
Author: 霉弟top 时间: 2024-5-30 01:04
both.sh 脚本输出不了坐标呀，终端会报错：sed: -e expression #1, char 1: unknown command: `,'

作者
Author: 123qwertybobo 时间: 2024-5-30 16:20

霉弟top 发表于 2024-5-30 01:04
both.sh 脚本输出不了坐标呀，终端会报错：sed: -e expression #1, char 1: unknown command: `,'

可能是有的原子没有定义，可以按照如Chlorine原子对应的gsub(17,"Cl",$2)来补充定义。不行的话上传源文件。

作者
Author: 霉弟top 时间: 2024-5-31 01:44

123qwertybobo 发表于 2024-5-30 16:20
可能是有的原子没有定义，可以按照如Chlorine原子对应的gsub(17,"Cl",$2)来补充定义。不行的话上传源文件 ...

我把both脚本中Standard orientation 换成Input orientation就可以正常使用了

作者
Author: 123qwertybobo 时间: 2024-5-31 11:03

霉弟top 发表于 2024-5-31 01:44
我把both脚本中Standard orientation 换成Input orientation就可以正常使用了

你确定你算完了么？还是说算的不是柔性扫描的结果。选Input orientation一点都不正常，那是你输入的结构，你要搜索下Standard orientation有没有。

作者
Author: 霉弟top 时间: 2024-5-31 12:30
那烦请您帮我看看这个柔性扫描输出文件有问题吗？十分感谢

作者
Author: 123qwertybobo 时间: 2024-5-31 16:40

霉弟top 发表于 2024-5-31 12:30
那烦请您帮我看看这个柔性扫描输出文件有问题吗？十分感谢

由于添加nosymm，G09和G16均不进行坐标标准化，只会保留Input orientation。由于Input和Standard结构是一样的，所以改成Input orientation就没问题。

作者
Author: 霉弟top 时间: 2024-5-31 18:06

123qwertybobo 发表于 2024-5-31 16:40
由于添加nosymm，G09和G16均不进行坐标标准化，只会保留Input orientation。由于Input和Standard结构是一 ...

感谢分析解答！！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)