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.

Abstract: The thermal expansion alpha(T) in layered crystals is of fundamental and technological interest. As suggested by I. M. Lifshitz in 1952, in thin solid films (crystalline membranes) a negative contribution to alpha(T) is due to anharmonic couplings between in-plane stretching modes and out-of-plane bending (flexural modes). Genuine in-plane anharmonicities give a positive contribution to alpha(T). The competition between these two effects can lead to a change of sign (crossover) from a negative value of alpha(T) in a temperature (T) range T <= T-alpha to a positive value of alpha(T) for T > T-alpha in layered crystals. Here, we present an analytical lattice dynamical theory of these phenomena for a two-dimensional (2D) hexagonal crystal. We start from a Hamiltonian that comprises anharmonic terms of third and fourth order in the lattice displacements. The in-plane and out-of-plane contributions to the thermal expansion are studied as functions of T for crystals of different sizes. Besides, renormalization of the flexural mode frequencies plays a crucial role in determining the crossover temperature T-alpha. Numerical examples are given for graphene where the anharmonic couplings are determined from experiments. The theory is applicable to other layer crystals wherever the anharmonic couplings are known. (C) 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

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: A unified theory of phonon dispersions and piezoelectricity in bulk and multilayers of hexagonal boron nitride (h-BN) is derived. The dynamical matrix is calculated on the basis of an empirical force constant model of intralayer valence and interlayer van der Waals interactions. Coulomb interactions are calculated by Ewalds method, adapted for the three-dimensional (3D) and the multilayer case. The deformation of the ionic charge distribution with long-wave lattice displacements is taken into account. Special attention is devoted to the nonanalytic long-range Coulomb contribution to the dynamical matrix which is different for the 3D crystal and the multilayer case. Consequently there is a splitting of the transverse optical (TO) and longitudinal optical (LO) phonon branches of E1u symmetry and a discontinuity of the A2u branch at the Γ point in 3D h-BN. No such splitting and discontinuity at Γ are present in multilayer crystals with a finite number N of layers. There a diverging bundle of N overbending optical phonon branches emerges from Γ. Borns long-wave theory is applied and extended for the study of piezoelectricity in layered crystals. While 3D h-BN and h-BN multilayers with an even number of layers (symmetry D6h) are not piezoelectric, multilayers with an uneven number of Nu layers (symmetry D3h) are piezoelectric; the piezoelectric coefficient e1,11 is inversely proportional to Nu.

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: 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.

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: The confinement of a C-60 molecule encapsulated in a cylindrical nanotube depends on the tube radius. In small tubes with radius R-T less than or similar to 7 A, a fivefold axis of the molecule coincides with the tube axis. The interaction between C-60 molecules in the nanotube is then described by a O-2-rotor model on a 1D liquid chain with coupling between orientational and displacive correlations. This coupling leads to chain contraction. The structure factor of the 1D liquid is derived. In tubes with a larger radius the molecular centers of mass are displaced off the tube axis. The distinction of two groups of peapods with on- and off-axis molecules suggests an explanation of the apparent splitting of A(g) modes of C-60 in nanotubes measured by resonant Raman scattering.