|
本帖最后由 冰释之川 于 2020-5-9 09:02 编辑
最近手头有一批复合物体系需要利用xtb软件计算其吸附质与吸附剂之间的相互作用能,这就需要复合物、吸附质与吸附剂三者的坐标文件(*.xyz)。为了能够快速的从复合物坐标中分别提取出吸附质与吸附剂坐标,靠手动一个个复制粘贴是一件很枯燥的事情,因而我写了一个从复合物xyz坐标文件中提取片段坐标的py脚本,从而避免不必要的重复劳作,大幅提升工作效率。
值得注意的是,本脚本中的“r”选项可以选择所有残余的原子标号,这对拆分相同原子数的吸附剂但是吸附质原子数各异的一大批复合体系来说方便不少。
举个例子:
1.xyz文件中1-50号原子属于吸附剂,51-60号原子属于吸附质;
2.xyz文件中1-50号原子属于吸附剂,51-80号原子属于吸附质;
为了能够一键提取这两个体系的吸附质与吸附剂坐标,第一个片段我们照常输入"1-50"选择所有的吸附剂原子;
而对于第二片段,我们可以直接输入"r",让脚本自动计算剩余未选择的原子的编号并将其归为吸附质原子。
下载地址见附件:
XYZ_Coordinate_Extraction.7z
(1.63 KB, 下载次数 Times of downloads: 33)
|
评分 Rate
-
查看全部评分 View all ratings
|