Abstract: Starting from an empirical force constant model of valence interactions and calculating by Ewald's method the electrostatic force constants, we derive the dynamical matrix for a monolayer and for multilayer systems of hexagonal boron nitride (h-BN). Solution of the secular problem leads to the corresponding phonon dispersion relations. The interplay between valence forces and Coulomb forces is discussed. A comparison with previous results on graphene and graphene multilayers is made. Our spectra on the h-BN monolayer are rather similar to previous ab initio theory results. Comparison is also made with Raman and infrared experimental results.

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: Born's long wave method is used to study the elastic properties of graphite and graphene. Starting from an empirical force constant model derived from full inplane phonon dispersions of graphite [Mohr et al., Phys. Rev. B 76, 035439 (2007)] we calculate the tension coefficients of graphene. Extending the model by interplanar interactions, we calculate the elastic constants of graphite. The agreement of our theoretical values with inelastic x-ray scattering results on elastic constants of graphite [Bosak et al., Phys. Rev. B 75, 153408 (2007)] is very satisfactory.

Abstract: For many decades the Boltzmann distribution function has been used to calculate the non-equilibrium properties of mobile particles undergoing the combined action of various scattering mechanisms and externally applied force fields. When the latter give rise to the occurrence of inhomogeneous potential profiles across the region through which the particles are moving, the numerical solution of the Boltzmann equation becomes a highly complicated task. In this work we highlight a particular algorithm that can be used to solve the time dependent Boltzmann equation as well as its quantum mechanical extension, the WignerBoltzmann equation. As an illustration, we show the calculated distribution function describing electrons propagating under the action of both a uniform and a pronouncedly non-uniform electric field.