Interaction of boron with graphite: A van der Waals density functional study

Boron doping has been widely investigated to improve oxidation resistance of graphite. In this work the interaction of boron with graphite is investigated by a van der Waals density-functional approach (vdW-DF). The traditional density-functional theory (DFT) is well accounted for the binding in boron-substituted graphite. However, to investigate the boron atom on graphite surface and the interstitial impurities require use of a description of graphite interlayer binding. Traditional DFT cannot describe the vdW physics, for instance, GGA calculations show no relevant binding between graphite sheets. LDA shows some binding, but they fail to provide an accurate account of vdW forces. In this paper, we compare the calculation results of graphite lattice constant and cohesive energy by several functionals, it shows that vdW-DF such as two optimized functionals optB88-vdW and optB86b-vdW give much improved results than traditional DFT. The vdW-DF approach is then applied to study the interaction of boron with graphite. Boron adsorption, substitution, and intercalation are discussed in terms of structural parameters and electronic structures. When adsorbing on graphite surface, boron behaves as π electron acceptor. The π electron approaches boron atom because of more electropositive of boron than carbon. For substitution situation, the hole introduced by boron mainly concentrates on boron and the nearest three carbon atoms. The B-doped graphite system with the hole has less ability to offer electrons to oxygen, ultimately resulted in the inhibition of carbon oxidation. For interstitial doping, vdW-DFs show more accurate formation energy than LDA. PBE functional cannot describe the interstitial boron in graphite reasonably because of the ignoring binding of graphite sheets. The investigation of electron structures of boron doped graphite will play an important role in understanding the oxidation mechanism in further study.