Abstract: A microscopic theory for the coupling of molecular orientations with acoustic lattice displacements is proposed. The relevant interaction is quadratic in the orientational-order-parameter variables and linear in longitudinal strains. The coupling is evaluated for a complex molecular structure. The intermolecular potential is obtained from interaction centers placed at nuclei and at the centers of interatomic bonds. The free energy is derived and the experimental consequences of the theory are studied. The Clausius-Clapeyron equation for orientational melting is obtained. The theory explains the decrease of the lattice constant at the first-order phase transition and the increase of the transition temperature with applied pressure.

Abstract: The recently discovered structural transition in polymerized KC60 at about 50 K results in a doubling of the unit cell volume and accompanies the 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 small orientational charge density waves along the polymer chains and concomitant displacements of the surrounding K+ ions. The effect is specific for the space group Pmnn of KC60 and is absent in RbC60 and CsC60 (space group I2/m). The mechanism is relevant for the metal-insulator transition.

Abstract: Neutron and high resolution X-ray diffraction investigations on perfect single crystals of RbH2PO4 (RDP), a hydrogen bonded ferroelectric of KDP type are reported. The results of crystal structure analysis from diffraction data, below and above the paraelectric – ferroelectric phase transition, support a disorder – order character Of [PO4H2](-)-groups. The tetragonal symmetry of the paraelectric phase with the double well potential of the hydrogen atoms obtained by diffraction, results simply from a time-space average of orthorhombic symmetry. According to the group – subgroup relation between the tetragonal space group 142d and the orthorhombic Fdd2 a short range order of ferroelectric clusters in the tetragonal phase is observed. With decreasing temperature the ferroelectric clusters increase and the long range interaction between their local polarisation vectors leads to the formation of lamellar ferroelectric domains with alternating polarisation directions at T-C = 147 K. From the high resolution X-ray data it is concluded that below T-C the ferroelastic strain in the (a,b)-plane leads to micro-angle grain boundaries at the domain walls. The tilt angle is enhanced by an applied electric field parallel to the ferroelectric axis. The resulting dislocations at the domain walls persist in the paraelectric phase leading to a memory effect for the arrangement of twin lamellae. With increased electric field the phase transition temperature T-C is decreased.

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

Abstract: The transition Fm (3) over bar m -->Pa (3) over bar in solid C-60 is driven by the condensation of orientational modes belonging to X-5(+) irreducible representations (irreps) of Fm (3) over bar m. Taking into account irreps up to the manifold l = 12, we have studied the primary and secondary orientational order parameters loops). We have numerically solved the coupled molecular field equations for these oops and calculated the temperature dependence of Bragg reflection intensities. (C) 1998 Elsevier Science B.V. All rights reserved.

Abstract: Differential cross sections for neutron and x-ray scattering have been derived for the orientationally disordered phase of solid C60. Interaction centres are placed at nuclei and at the centres of interatomic bonds. Bragg and diffuse scattering cross sections, for single crystals and for powders, are formulated using symmetry-adapted rotator functions. Thermal averages are calculated taking account of crystal field effects. Thermally averaged orientational distribution functions have also been calculated.

Abstract: From the theory of many-electron states in atoms, we know that there exists a strong Coulomb repulsion, which results in the electronic term structure of atoms and is responsible for Hunds rules. By expanding the Coulomb on-site repulsion into a multipolar series, we derive this interaction and show that it is also present in solids as a correlation effect, which means that the interaction requires a multideterminant version of the Hartree-Fock method. Of particular interest is the case where this interaction couples states of localized ( f) and delocalized ( s) electrons. We show that the interaction is bilinear in the creation/annihilation operators for localized electrons and bilinear in the operators for conduction electrons. To study the coupling, we consider a simple model in the framework of an effective limited configuration interaction method with one localized f-electron and one itinerant s-electron per crystal site. The on-site multipole interaction between the f- and s-electrons is explicitly taken into account. It is shown that depending on the low-lying excitation spectrum imposed by the crystal electric field, the model can lead not only to ferromagnetism but also to a nonmagnetic state. The model is relevant for solids with localized and itinerant electron states.

Abstract: Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semiempirical model for anharmonic couplings of third and fourth orders, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Gruneisen parameters, the thermal tension and the thermal expansion coefficients as a function of temperature and crystal size. From the resonances of the displacement-displacement correlation functions, we obtain the renormalization and decay rate of in-plane and flexural phonons as a function of temperature, wave vector, and crystal size in the classical and in the quantum regime. Quantitative results are presented for graphene. There, we find that the transition temperature T-alpha from negative to positive thermal expansion is lowered with smaller system size. Renormalization of the flexural mode has the opposite effect and leads to values of T-alpha approximate to 300 K for systems of macroscopic size. Extensive numerical analysis throughout the Brillouin zone explores various decay and scattering channels. The relative importance of normal and umklapp processes is investigated. The work is complementary to crystalline membrane theory and computational studies of anharmonic effects in two-dimensional crystals.

Abstract: Using time-differential perturbed angular correlation spectroscopy we have measured the electric field gradient (EFG) at 111Cd probe nuclei in solid Ce in a pressure range up to 8 GPa. Covering various allotropic phases of Ce, we find that the value of the EFG in the cubic α phase is almost four times larger than in the cubic γ phase and close to values in the noncubic phases α′ and α″. These results together with the differences in time modulation of the spectra are interpreted as evidence for quadrupolar electronic charge-density ordering and symmetry lowering at the γ→α transition while the lattice remains face-centered cubic

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.