Boundary conditions for transition-metal dichalcogenide monolayers in the continuum model

We derive the boundary conditions for MoS$_2$ and similar transition-metal dichalcogenide honeycomb (2H polytype) monolayers with the same type of $\mathbf{k}\!\cdot\!\mathbf{p}$ Hamiltonian within the continuum model around the K points. [1] In an effective 2-band description, the electron-hole symmetry breaking quadratic terms are also taken into account. We model the effect of the edges with a linear edge constraint method that has been applied previously to graphene. Focusing mainly on zigzag edges, we find that different reconstruction geometries with different edge-atoms can generally be described with one scalar parameter varying between 0 and $2\pi$. We analyze the edge states and their dispersion relation in MoS$_2$ in particular, and we find good agreement with the results of previous density functional theory calculations for various edge types. \\ {[1]} Cs.~G.~P\'eterfalvi, A.~Korm\'anyos, G.~Burkard, arXiv:1509.00184 (2015).