Abstract: We have developed a microscopic theory which describes the orientational dynamics of C-70 molecules and its coupling to lattice displacements in the face-centered-cubic phase of C-70 fullerite. The single-molecule orientational density distribution in the disordered phase is calculated. The ferroelastic transition to the rhombohedral phase is investigated. The discontinuity of the orientational order parameter at the phase transition is calculated. It is found that the transition leads to a stretching of the primitive unit cell along a [111] cubic direction. A softening of the elastic constant c(44) at the transition is predicted.

Abstract: Starting from a model of microscopic interactions between C-70 molecules, we have developed a theory which describes the orientational dynamics and its coupling to lattice displacements in the rhombohedral phase of C-70 fullerite. The Landau free energy is calculated. We obtain a first-order phase transition to a monoclinic structure with the space group P2(1)/m. The transition is driven by the condensation of orientational quadrupoles at the F point of the Brillouin zone of the rhombohedral lattice. We find no evidence that the monoclinic structure is connected with the freezing in of orientations around the fivefold molecular axis. We calculate the lattice strains that are associated with the transition to the monoclinic structure. The theory is compared with a range of experimental data on the phase transition.

Abstract: The orientation of a C(70) fullerene molecule encapsulated in a single-walled carbon nanotube (SWCNT) depends on the tube radius. First we confirm that chirality effects do not affect the orientation as well by comparing discrete atomistic calculations with the results of a continuous tube approximation for a variety of SWCNTs. The molecular and the tube symmetry are exploited by using symmetry-adapted rotator functions. We accurately determine the optimal molecular orientation as a function of the tube radius; for low (less than or similar to 7 A) and high (greater than or similar to 7.2 A) tube radii, lying and standing molecular orientations are recovered, respectively. In between, we observe a transition regime. In addition, we consider off-axis positions. We perform a one-dimensional liquid description of a chain of on-axis C(70) molecules inside a SWCNT. All results agree well with recent x-ray diffraction experiments.

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: 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: The crystal field in the orientationally disordered phase of C90-fullerite is derived from an intermolecular potential model, which takes into account the geometric difference between double bonds and single bonds. The molecules are modelled as rigid bodies, atoms and single bonds are treated as single interaction centers, while double bonds are described by a distribution of interaction centers along the bond. The crystal field is expanded in terms of cubic rotator functions. The calculated expansion coefficients are compared with empirical values derived from diffraction data. The angular dependence of the crystal field, resulting from an anticlockwise rotation of the molecule around the [111] axis, exhibits an absolute and a secondary minimum at angles of 98-degrees and 38-degrees respectively. The self interaction of the molecule in a deformable lattice is investigated.

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