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 lattice dynamics of low-frequency rigid-plane modes in metallic (graphene multilayers, GML) and in insulating (hexagonal boron-nitride multilayers, BNML) layered crystals is investigated. The frequencies of shearing and compression (stretching) modes depend on the layer number N and are presented in the form of fan diagrams. The results for GML and BNML are very similar. In both cases, only the interactions (van der Waals and Coulomb) between nearest-neighbor planes are effective, while the interactions between more distant planes are screened. A comparison with recent Raman scattering results on low-frequency shear modes in GML [Tan et al., Nat. Mater., in press, doi: 10.1038/nmat3245, (2012)] is made. Relations with the low-lying rigid-plane phonon dispersions in the bulk materials are established. Master curves, which connect the fan diagram frequencies for any given N, are derived. Static and dynamic thermal correlation functions for rigid-layer shear and compression modes are calculated. The results might be of use for the interpretation of friction force experiments on multilayer crystals.

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: 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: The recently discovered structural phase transition in Polymerized KC60 at about 50K leads to a doubling of the unit cell volume and is accompanied by a metal-insulator transition. Here, we show that the ((a) over right arrow + (c) over right arrow, (b) over right arrow, (a) over right arrow – (c) over right arrow) superstructure results from orientational charge density waves along the polymer chains and correlated displacements of the K+ ions. The presented model can also account for the metal-insulator transition. The effect is specific for the space group Pmnn of KC60 and is absent in both Rb- and CsC60 (space group 12/m), in agreement with the present experimental knowledge of these compounds.

Abstract: Piezo- and flexoelectricity are manifestations of electromechanical coupling in solids with potential applications in nanoscale materials. Naumov etal. [Phys. Rev. Lett. 102, 217601 (2009)] have shown by first principles calculations that a monolayer BN sheet becomes macroscopically polarized in-plane when in a corrugated state. Here, we investigate the interplay of layer corrugation and in-plane polarization by atomistic lattice dynamics. We treat the coupling between static flexural modes and in-plane atomic ion displacements as an anharmonic effect, similar to the membrane effect that is at the origin of negative thermal expansion in layered crystals. We have derived analytical expressions for the corrugation-induced static in-plane strains and the optical displacements with the resulting polarization response functions. Beyond h-BN, the theory applies to transition metal dichalcogenides and dioxides. Numerical calculations show that the effects are considerably stronger for 2D h-BN than for 2H-MoS2.

Abstract: Despite the general belief that hydrogen-bonded ferro- and antiferroelectrics undergo a pure order/disorder transition at the structural instability, new NMR data and a new theoretical concept yield convincing evidence that a pronounced displacive component is present in these systems, which modifies substantially the temperature dependencies of the tunnel and lattice mode frequencies. The experiments and their interpretation are presented. (C) 1999 Elsevier Science Ltd. All rights reserved.

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.