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: 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: 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: 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 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: The elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

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 superexchange interactions between the fullerenes arise as a result of the electron transfer from the C-60 molecule to the alkali atom and back. We present a scheme, which is a configuration interaction approach based on the valence bond (Heitler-London) method. The effect of superexchange is described together with chemical bonding by constructing and solving a secular equation, rather than by using a perturbation treatment. We have considered 180degrees and 90degrees superexchange for the C-60 Cs-C-60 pathways. The calculations account for unusual electronic properties of polymer orthorhombic and quenched cubic phases of CsC60: two lines in nuclear magnetic resonance experiments, the development of a spin-singlet ground state and a decrease of magnetic susceptibility as T-->0.

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: 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: We report single-crystal X-ray diffuse scattering measurements on C-60.C8H8 fullerene-cubane showing that the C-60 molecules are orientationally disordered at 300 and 150K and get ordered at low temperatures. Monte Carlo simulations provide further insight in the orientational behavior of both C-60 and C8H8 molecules; low-temperature molecular orientations are predicted.