计算化学在处理实际化学问题时,比如需要在某一化学平面的法向量上进行分子操作,这时最重要的是确定化学平面和求法向量,才能进行后续的操作(如下图所示),下面以高斯输入文件为例,用python代码实现该功能,包括以下两部分:
1. 确定化学平面:一般通过三个点确定,也就是三个原子的坐标确定一个化学平面;
2. 计算法向量和单位法向量
以下输入文件为例:
以下信息复制到新文件,命名为xx.gjf,具体如下:
%nprocshared=20
%mem=2000MB
%chk=xx.chk
# ub3lyp/6-31g(d) geom=connectivity
Title Card Required
0 2
C -0.07799900 -0.06595600 1.16020400
C 0.28486000 0.24116000 -0.05689700
Sc -0.57500900 -2.04502600 0.35332800
Sc 2.15465200 0.28669700 1.05119200
Sc -1.43262200 1.53516600 0.33122200
C 1.07300100 -3.58961700 -0.99271100
C 0.64679900 -2.84267000 -2.13946700
C 2.30603300 -3.23292800 -0.33821000
C 3.11036400 -2.22051000 -0.93080800
C 2.69014500 -1.51047400 -2.10623700
C 1.43348100 -1.78981800 -2.72004200
C -0.09411100 -4.12378500 -0.30887600
C