wzkchem5 发表于 2024-1-19 20:43
（1）把你的分子结构存成xyz、gjf等格式，用记事本打开，把所有16个碳、氮原子对应的行都挪到最前面
（2） ...

好的，谢谢老师

sobereva 发表于 2024-1-19 20:19
自己写程序/脚本实现，参考下文
基于OpenBabel批量产生特定基团以任意方式接到苯上的结构的方法
http://s ...

好的，谢谢老师

（1）把你的分子结构存成xyz、gjf等格式，用记事本打开，把所有16个碳、氮原子对应的行都挪到最前面
（2）随便用某个脚本语言，按以下伪代码写一个脚本

1. for i=1,...,16
   
2. for j=i+1,...,16
   
3. 把你的结构文件拷一份
   
4. 把拷的这一份的分子结构信息部分里，第i行和第j行的C替换为N
   
5. end
   
6. end

复制代码

例如用bash写循环体结构，用sed实现替换，就行。当然用python等语言也能写
（3）批量提交计算即可。
以上批量把两个C替换成N。不难举一反三把三个、四个C替换成N

自己写程序/脚本实现，参考下文
基于OpenBabel批量产生特定基团以任意方式接到苯上的结构的方法
http://sobereva.com/440（http://bbs.keinsci.com/thread-10824-1-1.html）