Abstract: Starting from an empirical force constant model of valence interactions and calculating by Ewalds method the ion-ion force constants, we derive the dynamical matrix for a monolayer crystal of hexagonal boron nitride (h-BN). The phonon dispersion relations are calculated. The interplay between valence and Coulomb forces is discussed. It is shown by analytical methods that the longitudinal and the transverse optical (LO and TO) phonon branches for in-plane motion are degenerate at the Γ point of the Brillouin zone. Away from Γ, the LO branch exhibits pronounced overbending. It is found that the nonanalytic Coulomb contribution to the dynamical matrix causes a linear increase of the LO branch with increasing wave vector starting at Γ. This effect is general for two-dimensional (2D) ionic crystals. Performing a long-wavelength expansion of the dynamical matrix, we use Borns perturbation method to calculate the elastic constants (tension coefficients). Since the crystal is noncentrosymmetric, internal displacements due to relative shifts between the two sublattices (B and N) contribute to the elastic constants. These internal displacements are responsible for piezoelectric and dielectric phenomena. The piezoelectric stress constant and the dielectric susceptibility of 2D h-BN are calculated.

Abstract: We present a unified theory of the phonon dispersions and elastic properties of graphene, graphite, and graphene multilayer systems. Starting from a fifth-nearest-neighbor force-constant model derived from full in-plane phonon dispersions of graphite [Mohr et al., Phys. Rev. B 76, 035439 (2007)], we use Born's long-wave method to calculate the tension and bending coefficients of graphene. Extending the model by interplanar interactions, we study the phonon dispersions and the elastic constants of graphite, and the phonon spectra of graphene multilayers. We find that the inner displacement terms due to sublattice shifts between inequivalent C atoms are quantitatively important in determining the elastomechanical properties of graphene and of graphite. The overall agreement between theory and experiment is very satisfactory. We investigate the evolution from graphene to graphite by studying the increase in the rigid plane optical mode as a function of the number of layers N. At N=10 the graphite value B2g1127 cm−1 is attained within a few percent.

Abstract: The appearance of massless Dirac fermions in graphene requires two equivalent carbon sublattices of trigonal shape. While the generation of an effective mass and a band gap at the Dirac point remains an unresolved problem for freestanding extended graphene, it is well established by breaking translational symmetry by confinement and by breaking sublattice symmetry by interaction with a substrate. One of the strongest sublattice-symmetry-breaking interactions with predicted and measured band gaps ranging from 400 meV to more than 3 eV has been attributed to the interfaces of graphene with Ni and Co, which are also promising spin-filter interfaces. Here, we apply angle-resolved photoemission to epitaxial graphene on Ni (111) and Co(0001) to show the presence of intact Dirac cones 2.8 eV below the Fermi level. Our results challenge the common belief that the breaking of sublattice symmetry by a substrate and the opening of the band gap at the Dirac energy are in a straightforward relation. A simple effective model of a biased bilayer structure composed of graphene and a sublattice-symmetry-broken layer, corroborated by density-functional-theory calculations, demonstrates the general validity of our conclusions.

Abstract: The determination of the layer number ${\cal N}$ in nanoscale thin layered crystals is a challenging problem of technological relevance. In addition to innovative experimental techniques, a thorough knowledge of the underlying lattice dynamics is required. Starting from phenomenological atomic interaction potentials we have carried out an analytical study of the low-frequency optical phonon dispersions in layered crystals. At the gamma point of the two-dimensional Brillouin zone the optical phonon frequencies correspond to rigid-plane shearing and compression modes. We have investigated graphene multilayers (GML) and hexagonal boron-nitride multilayers (BNML). The frequencies show a characteristic dependence on ${\cal N}$. The results which are represented in the form of fan diagrams are very similar for both materials. Due to charge neutrality within layers Coulomb forces play no role, only van der Waals forces between nearest neighbor layers are relevant. The theoretical results agree with recent low-frequency Raman results on rigid-layer modes [Tan et al., Nature Mater. 11, 294 (2012)] in GML and double-resonant Raman scattering data on rigid-layer compression modes [Herziger et al., Phys. Rev. B 85, 235447 (2012)] in GML. (C) 2012 WILEY-VCH Verlag GmbH Co. KGaA, Weinheim

Abstract: Taking into account covalent, Coulomb and van der Waals interactions, we construct the dynamical matrix and calculate the phonon dispersion relations for multilayer crystals of hexagonal boron-nitride. Coulomb interactions account for a strong overbending of optical phonons. Applying and extending Born's long-wave theory to the case of multilayer crystals, we calculate the piezoelectric stress constant equation image as a function of the number of layers equation image. In agreement with group theory, we find that equation image for equation image even; for an uneven number equation image of layers we obtain equation image, i.e. the piezoelectric constant decreases as equation image.