Python实现:已知化学分子的输入文件坐标(高斯计算输入文件为例),求其中任意三个原子确定的平面的法向量和单位法向量

        计算化学在处理实际化学问题时,比如需要在某一化学平面的法向量上进行分子操作,这时最重要的是确定化学平面和求法向量,才能进行后续的操作(如下图所示),下面以高斯输入文件为例,用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              
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值