Abstract: We study the confinement of a C-60 molecule encapsulated in a cylindrical nanotube as a function of the tube radius. Drawing the Mercator maps of the potential, we find two distinct molecular orientations; for tubes with small radii, R-T less than or similar to 7 angstrom, a fivefold axis of the molecule coincides with the tube long axis, for larger radii, R-T less than or similar to 8 angstrom, a threefold axis of the molecule coincides with the tube long axis. These different orientations are caused by the relative importance of the repulsive and the attractive parts of the van der Waals potentials of the molecule with the tube wall for small and large tubes respectively. Experimental evidence is provided by the apparent splitting of A(g) modes of the C-60 molecule in resonant Raman scattering.

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: 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: Starting from a model of rigid interacting C-60 polymer chains on an orthorhombic lattice, we study the mutual orientation of the chains and the stability of the crystalline structures Pmnn and I2/m. We take into account (i) van der Waals interactions and electric quadrupole interactions between C-60 monomers on different chains as well as (ii) interactions of the monomers with the surrounding alkali atoms. The direct interactions (i) always lead to an antiferrorotational structure Pmnn with alternate orientation of the C-60 chains in planes (001). The interactions (ii) with the alkalis consist of two parts: translation-rotation (TR) coupling where the orientations of the chains interact with displacements of the alkalis, and quadrupolar electronic polarizability (ep) coupling, where the electric quadrupoles on the C-60 monomers interact with induced quadrupoles due to excited electronic d-states of the alkalis. Both interactions (ii) lead to an effective orientation-orientation interaction between the C-60 chains and always favor the ferrorotational structure I2/m, where C-60 chains have a same orientation. The structures Pmnn for KC60 and I2/m for Rb- and CsC60 are the result of a competition between the direct interaction (i) and the alkali-mediated interactions (ii). In Rb- and CsC60 the latter are found to be dominant, the preponderant role being played by the quadrupolar electronic polarizability of the alkali ions. (C) 2002 American Institute of Physics.

Abstract: C-60 fullerenes confined inside single-walled carbon nanotubes form an archetypal one-dimensional system. X-ray diffraction experiments, from room temperature to 1073 K, reveal an increasing melting phenomenon. Detailed analysis of the sawtooth peak characteristic of the fullerene organization allows the quantitative determination of fluctuations in intermolecular distances. The present results validate the predictions of one-dimensional statistical models.

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