边界的碳原子肯定在优化时得固定住
在Gaussian中做限制性优化的方法
http://sobereva.com/404（http://bbs.keinsci.com/thread-9022-1-1.html）

用团簇模型可以做，但是石墨烯的边界需要特殊处理，参见http://sobereva.com/540
另外提醒一下，从共振式可以预期，石墨烯的简单碳空位极其不稳定，实验上估计做不出这种缺陷，做出来也会马上重构。务必确认你研究的缺陷是可以稳定存在的，而且之前文献没算过，再去研究。