玩了一下新潮的计算软件ORCA(话说这主页有点蹩脚),等它算的过程中突然有点不确了:我的输入文件用的单位是不是Å?懒得去对了,我直接把输入文件里相邻两个碳原子间的距离算出来,看和键长是否一样不就知道了。

嗯,现在问题转化为计算两个点的欧氏距离了,因为有一个打开的Octave,工具就是Octave了。我首先写成这个样子:

1
C1 = [2.72700 4.54400 6.53900]; C2 = [2.00400 4.25100 5.43400]; sum((C1 - C2).^2)^0.5

太蹩脚。搜一下,发现有人问过了,别人给答案是:

1
C1 = [2.72700 4.54400 6.53900]; C2 = [2.00400 4.25100 5.43400]; norm(C1 - C2, 2)

嗯,不错。这个问题下面有人回复“STFW”,于是我也搜了一下,发现Matlab网站上给的几种写法:

1
2
3
4
5
6
7
8
9
10
11
C1 = [2.72700 4.54400 6.53900]; C2 = [2.00400 4.25100 5.43400];

#1
sqrt(sum((C1 - C2) .^ 2))

#2
V = C1 - C2;
sqrt(V * V')

#3
norm(C1 - C2)

第二种是新的,没学过线代不知道原理。

😢写了这么多才算了一轮,DLPNO-CCSD还是好慢啊。