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 elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

Abstract: Starting from our previous work in which we obtained a system of coupled integrodifferential equations for acoustic sound waves and phonon density fluctuations in two-dimensional (2D) crystals, we derive here the corresponding hydrodynamic equations, and we study their consequences as a function of temperature and frequency. These phenomena encompass propagation and damping of acoustic sound waves, diffusive heat conduction, second sound, and Poiseuille heat flow, all of which are characterized by specific transport coefficients. We calculate these coefficients by means of correlation functions without using the concept of relaxation time. Numerical calculations are performed as well in order to show the temperature dependence of the transport coefficients and of the thermal conductivity. As a consequence of thermal tension, mechanical and thermal phenomena are coupled. We calculate the dynamic susceptibilities for displacement and temperature fluctuations and study their resonances. Due to the thermomechanical coupling, the thermal resonances such as the Landau-Placzek peak and the second-sound doublet appear in the displacement susceptibility, and conversely the acoustic sound wave doublet appears in the temperature susceptibility, Our analytical results not only apply to graphene, but they are also valid for arbitrary 2D crystals with hexagonal symmetry, such as 2D hexagonal boron nitride, 2H-transition-metal dichalcogenides, and oxides.

Abstract: We have developed a microscopic theory which describes the orientational dynamics of C60 molecules in the face-centered-cubic phase of C60 fullerite. The molecular interaction potential and the crystal-field potential are formulated in terms of symmetry-adapted rotator functions. The phase transition to the Pa3BAR structure is driven by an active multipolar mode of T2g symmetry belonging to the l = 10 manifold. The Birman criterion is satisfied. The transition is found to be of first order.