Atomistic Insights into Nucleation and Formation of Hexagonal Boron Nitride on Nickel from First-Principles-Based Reactive Molecular Dynamics Simulations

Atomistic-scale insights into the growth of a continuous, atomically thin hexagonal boron nitride (hBN) lattice from elemental boron and nitrogen on Ni substrates were obtained from multiscale modeling combining density functional theory (DFT) and reactive molecular dynamics. The quantum mechanical calculations focused on the adsorption and reaction energetics for the hBN building-block species, i.e., atomic B, N, B x N y (x, y = 1, 2), on Ni(111) and Ni(211), and the diffusion pathways of elemental B and N on these slab model surfaces and in the sublayer. B can diffuse competitively on both the surface and in the sublayer, while N diffuses strictly on the substrate surface. The DFT data were then used to generate a classical description of the Ni–B and Ni–N pair interactions within the formulation of the reactive force field, ReaxFF. Using the potential developed from this work, the elementary nucleation and growth process of an hBN monolayer structure from elemental B and N is shown at the atomistic scale. The nucleation initiates from the growth of linear BN chains, which evolve into branched and then hexagonal lattices. Subsequent DFT calculations confirmed the structure evolution energetically and validate the self-consistency of this multiscale modeling framework. On the basis of this framework, the fundamental aspects regarding crystal quality and the role of temperature and substrates used during hBN growth can also be understood.