以Al为例计算平衡晶格常数和结合能

——作者水平有限，仅供学习交流使用，如有不足还请指出

本文以Al为例分子模拟计算了面心立方(fcc)晶格结构中铝的平衡晶格常数和相应的结合能。

首先明确晶格常数和结合能的概念，对于晶格常数，我相信各位早已烂熟于心，这里不再过多赘述。对于结合能，我简单去介绍一下其定义：若某个物体由两部分及以上组成时，各部分必然存在一定吸引力使其结合在一起，要将各部分分开至无穷远必然对其做功，所做的功的大小即表示他们的结合能。单位为能量单位(eV)，对于平均结合能其计算公式为体系总能量/粒子总数。

运行结束后会在屏幕以及log文件中输出。