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: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: The structural phase transitions with electric-quadrupole long-range order in NpO2 (Fm (3) over barm-->Pn (3) over barm) and UO2 (Fm (3) over barm-->Pa (3) over bar) are analyzed from a group theoretical point of view. In both cases, the symmetry lowering involves three quadrupolar components belonging to the irreducible representation T-2g (Gamma(5)) of O-h and condensing in a triple-q structure at the X point of the Brillouin zone. The Pa (3) over bar structure is close to Pn (3) over barm, but allows for oxygen displacements. The Pa (3) over bar ordering leads to an effective electrostatic attraction between electronic quadrupoles while the Pn (3) over barm ordering results in a repulsion between them. It is concluded that the Pn (3) over barm structure can be stabilized only through some additional process such as strengthening of the chemical bonding between Np and O. We also derive the relevant structure-factor amplitudes for Pn (3) over barm and Pa (3) over bar, and the effect of domains on resonant x-ray scattering experiments.

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: 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: The electronic properties of boron nitride nanoribbons (BNNRs) doped with a line of carbon atoms are investigated using density functional calculations. By replacing a line of alternating B and N atoms with carbons, three different configurations are possible depending on the type of the atoms which bond to the carbons. We found very different electronic properties for these configurations: (i) the NCB arrangement is strongly polarized with a large dipole moment having an unexpected direction, (ii) the BCB and NCN arrangements are nonpolar with zero dipole moment, (iii) the doping by a carbon line reduces the band gap regardless of the local arrangement of the borons and the nitrogens around the carbon line, and (iv) the polarization and energy gap of the carbon-doped BNNRs can be tuned by an electric field applied parallel to the carbon line. Similar effects were found when either an armchair or zigzag line of carbon was introduced.

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

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: Polymerized KC60 undergoes a structural phase transition accompanied by a metal-insulator transition around 50 K. To explain the structural aspect, a mechanism involving small orientational deviations of the valence electron density on every C-60 monomer orientational charge density waves (OCDWs) – has already been proposed earlier. In the present work, we address the metal-insulator transition using the OCDW concept. We are inspired by the analogy between a polymer chain exhibiting an OCDW and a linear atomic chain undergoing a static lattice deformation doubling the unit cell: such a deformation implies a band gap at the zone boundary, yielding an insulating state (Peierls instability). Within our view, a similar mechanism occurs in polymerized KC60; the OCDW plays the role of the lattice deformation. We present tight-binding band structure calculations and conclude that the metal-insulator transition can indeed be explained using OCDWs, but that the threedimensionality of the crystal plays an unexpected key role.