计算化学公社

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

[Gaussian/gview] 给GaussView添加自定义基团模板

[复制链接 Copy URL]

178

帖子

3

威望

1916

eV
积分
2154

Level 5 (御坂)

GaussView自带一些常用的基团模板,主要分为ring fragment,R-group fragment和biological fragment三类,各自有对应的缩略图方便随时取用。但实际操作中还是缺少一些常用的片段,尤其是处理某一个特定的课题时,如果能添加一些特殊的基团可能会大大方便建模。
GaussView官方推荐的是Custom fragment,有帖子已经提到过,具体使用方式在GaussView的help中有。但Custom fragment只能添加模板,从下拉菜单中选择,却不能像内置模板那样有一个直观的缩略图,在选择的时候不那么方便。本帖介绍如何优雅地添加带缩略图的自定义基团。本帖基于Linux和MacOS系统的GaussView 6,对于Windows应该也是可行的。


1. GaussView内置基团模板文件在哪里?
这些文件在以下三个文件夹中,分别对应GaussView面板上的这三个按钮:
gv/data/rings & gv/data/rgoups & gv/data/biofrags


2. GaussView基团模板文件的格式是什么?
以gv/data/rings为例,其中包含bitmaps和fragments两个文件夹,以及rings.ftb。
rings.ftb中包含了所有36个内置的环状化合物模板,内容如下。
  1. # Rings Fragment Table
  2. #   Format: name pixmap fragment row col [atnum hotx hoty]
  3. #   Tabstops = 20, Indexes = 0 ... n-1
  4. nItems = 36 defaultItem = 0
  5. "benzene"  phenyl.xbm  phenyl.frg  0  0 # 基团名称 缩略图文件 结构文件 行号 列号
复制代码
以benzene为例,对应的缩略图存在bitmaps/phenyl.xbm,对应的结构文件存在fragments/phenyl.frg,"0 0"表示在GaussView面板中点开ring fragment后,苯环出现在第1行第1列(行列号从0开始计数)。如果我们要在ring fragment中添加一个酞菁铁的模板,首先修改nItems = 37,然后追加一行:
  1. "ironpc" ironpc.xbm ironpc.frg 6 0
复制代码
此行将会把酞菁铁添加到GaussView面板中的第7行第1列。

.xbm是一种古老的位图文件格式,人是读不懂的,不过我们可以将.png/.jpg/.tif等常用图片格式转换成.xbm格式。准备一个酞菁铁的结构图,并保存为ironpc.png。下载本帖提供的mkxbm.py并运行:
  1. python mkxbm.py ironpc.png
复制代码
将得到的ironpc.xbm放到gv/data/rings/bitmaps/ironpc.xbm

.frg是GaussView自己的片段文件格式,内容主要包含分子坐标和连接关系。首先用GaussView搭建一个酞菁铁的结构,保存为ironpc.mol2。下载本帖提供的mkfrg.py并运行:
  1. python mkfrg.py ironpc.mol2
复制代码
注意:标准的.mol2格式中芳香化合物中的共轭C-C键的键型记为ar,GaussView保存的.mol2文件则记为Ar,但.frg格式中记为4。
将得到的ironpc.frg放到gv/data/rings/fragments/ironpc.frg。
含有连接关系的文件格式很多,此处用mol2处理比较简洁,反正GaussView正好支持.mol2格式。

3. 结果
重启GaussView,这时候会看到ring fragment下多了一个酞菁铁的模板。实践表明清晰度感人,但勉强能辨识。这主要是转换成位图的时候的一些问题,如果需要更清晰的缩略图可能需要在一开始提供图片时就考虑一大堆问题,也可能需要借助更高级的转换工具,此处既然能辨识就不再折腾了。如果有大量模板,只需要写一个for循环处理一下就行。


mkfrg.py (2.08 KB, 下载次数 Times of downloads: 3)
mkxbm.py (2.83 KB, 下载次数 Times of downloads: 2)
ironpc.frg (4.13 KB, 下载次数 Times of downloads: 2)
ironpc.xbm (2.96 KB, 下载次数 Times of downloads: 1)

评分 Rate

参与人数
Participants 6
威望 +1 eV +23 收起 理由
Reason
LittlePupil + 3 GJ!
thou + 5 赞!
丁越 + 5 赞!
sobereva + 1
Kevin_Chen + 5 赞!
Uus/pMeC6H4-/キ + 5 好物!

查看全部评分 View all ratings

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

2#
发表于 Post on 2026-1-4 14:14:29 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2026-1-4 14:17 编辑

在我之前将数据集的金属配合物做成GaussView片段时也发现过扩展GaussView模板的方法、分析过.frg格式的特征,不过没把批量处理的脚本贴出来……

提醒一下mkfrg.py依赖numpy库,mkxbm.py依赖pillow库(PIL),没有的话需要pip install安装。

编辑:除了用mol2指定连接关系外,提供原子坐标再按gview自动成键判定的依据来生成连接关系或许也是个方法。
√546=23.36664289109

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 13:12 , Processed in 0.199128 second(s), 25 queries , Gzip On.

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