计算化学公社

标题: Gaussian输入文件关键词替换脚本、单点能数据提取后处理脚本 [打印本页]

作者
Author: gaogaoamm 时间: 2025-5-7 14:00
标题: Gaussian输入文件关键词替换脚本、单点能数据提取后处理脚本
本帖最后由 gaogaoamm 于 2025-5-7 14:05 编辑

（一）问题背景（适用计算场景）

问题一：

笔者最近在做单点能计算，因想临时更换计算精度，需要对Gaussian输入文件更换关键词，需要批量对输入文件进行编辑，因此写了一个脚本辅助处理。

问题二：

在单点能计算过程中，笔者只需要输出文件中的MP2和QCISD(T)数据，所以也写了一个python脚本进行批处理。

（二）代码介绍

2.1 关键词更换脚本（命名为replace.py放在附件中）
只需将此py文件和待处理的输入文件放在同一文件夹下，运行即可生成。
2.2 单点能后处理脚本（命名为LINE.py和MP2.py放在附件中）
只需将此py文件和待处理的输入文件放在同一文件夹下，运行即可生成。

思路说明：
（1）代码实现逻辑为直接读取文件中“MP2=”和“QCISD(T)=”以及后面带的数据，并输出到excel表中。
（2）基于逻辑与笔者代码水平，故先用LINE.py文件读取输出文件的后50行，排除了文件中其他“MP2=”和“QCISD(T)=”关键词的干扰，保证MP2.py读取的为最终计算出的单点能数据。
（3）局限性：因为有些单点能数据恰好呈现下图上图中的样式（即要读取的数据不在同一行）：所以读取会读取出空白如下图。实测下来，这种格式发生率在2%左右，可以自行打开文件输入。

(, 下载次数 Times of downloads: 123) (, 下载次数 Times of downloads: 123)

附件：
(, 下载次数 Times of downloads: 3)
(, 下载次数 Times of downloads: 4)
(, 下载次数 Times of downloads: 6)

作者
Author: imasen 时间: 2025-5-7 16:48
以前写过一个提取单点能的shell脚本，可以避免行末换行的问题。使用时把脚本里的HF全部替换成MP2就可以了。 (, 下载次数 Times of downloads: 7)

作者
Author: gaogaoamm 时间: 2025-5-7 16:52

imasen 发表于 2025-5-7 16:48
以前写过一个提取单点能的shell脚本，可以避免行末换行的问题。使用时把脚本里的HF全部替换成MP2就可以了。

谢谢！已经用成功了！

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