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: 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: A theoretical model is presented for the description of the metal-insulator transition which accompanies the structural phase transition at T approximate to 50 K in polymerized KC60. The model involves orientational charge density waves (along the C-60 polymer chains) which were introduced previously for a description of the structural phase transition. A satisfactory qualitative and quantitative understanding is obtained when the three-dimensionality of the crystal and the presence of the K+ counterions is properly taken into account.

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: 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 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: Superexchange interactions in alkali fullerides AC(60) are derived for C-60 molecular ions separated by interstitial alkali-metal ions. We use a multiconfiguration approach which comprises the lowest molecular orbital states of the C-60 molecule and the excited s and d states of the alkali-metal atom A. Interactions are described by the valence bond (Heitler-London) method for a complex (C-60 – A – C-60) – with two valence electrons. The electronic charge transfer between the alkali-metal atom and a neighboring C-60 molecule is not complete. The occupation probability of excited d and s states of the alkali atom is not negligible. In correspondence with the relative positions of the C-60 molecules and A atoms in the polymer crystal, we consider 180degrees and 90degrees (angle) superexchange pathways. For the former case the ground state is found to be a spin singlet separated from a triplet at similar to20 K. For T < 20 K there appear strong spin correlations for the 180degrees superexchange pathway. The results are related to spin lattice relaxation experiments on CsC60 in the polymerized and in the quenched cubic phase. (C) 2005 American Institute of Physics.

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: 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: The influence of triaxial in-plane strain on the electronic properties of a hexagonal boron-nitride sheet is investigated using density functional theory. Different from graphene, the triaxial strain localizes the molecular orbitals of the boron-nitride flake in its center depending on the direction of the applied strain. The proposed technique for localizing the molecular orbitals that are close to the Fermi level in the center of boron nitride flakes can be used to actualize engineered nanosensors, for instance, to selectively detect gas molecules. We show that the central part of the strained flake adsorbs polar molecules more strongly as compared with an unstrained sheet.

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: 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: The giant isotope effect on the ferro- and antiferroelectric transition temperature upon deuteration of hydrogen-bonded systems is well known experimentally since various decades. Yet, theoretically only recently a microscopic understanding of this effect has been achieved which, specifically, took into account the geometry of the O ... H ... O bond. The implications of this modeling are multiple as numerous hydrogen-bonded organic systems show the same effects as ferro- and antiferroelectrics, i.e., cooperative proton tunneling at a well-defined temperature.