你似乎在做钙钛矿太阳能电池相关的计算？这恰好是我的领域。
文章中提到的结合能其实是binding free energy，而out文件中直接给出的这个值是interaction energy。获取binding free energy的方法可以参考《使用Shermo结合量子化学程序方便地计算分子的各种热力学数据》

顺带一提，看到你对Pb用了LANL2DZ基组，这是极不推荐的。这个基组的精度如今看来已经明显过差，对Pb应当至少用Lanl2TZ。（使用混合基组时赝势基组通常需要比全电子基组更好，且当前就是在讨论Pb和的结合能，Pb的基组精度很重要。）详见：《谈谈赝势基组的选用》。当前也没有恰当定义赝势，所以会算不出来。定义方法见：《详解Gaussian中混合基组、自定义基组和赝势基组的输入》：

6 使用赝势基组
赝势基组需要和对应的赝势一起使用，如果不熟悉赝势或想深入了解内在细节，可参见《赝势的函数形式以及在量子化学程序中定义的方式》（http://sobereva.com/188）。使用赝势基组需要用genecp关键词（等价于gen pseudo=read），代表先从分子坐标后面读取赝势基组定义，再读取赝势定义。例如Cu(CO)+：
# B3LYP/genecp
[空行]
Cu(CO)+
[空行]
1 1
Cu            
C                  1            B1
O                  2            B2    1            180.
[空行]
   B1             1.94000000
   B2             1.11540000
[空行]
Cu 0
Lanl2DZ    <-- Cu用Lanl2DZ赝势基组
****
C O 0
6-31G*   <-- C O用的基组为6-31G*
****
[空行]
Cu 0
Lanl2   <-- Cu用的赝势为Lanl2，这种赝势是Lanl2DZ赝势基组对应的赝势。如果不知道赝势名是什么，此处直接写赝势基组名就可以。比如这里也可以写成Lanl2DZ，和写Lanl2是等价的
[空行]
[空行]