Abstract: Piezo- and flexoelectricity are manifestations of electromechanical coupling in solids with potential applications in nanoscale materials. Naumov etal. [Phys. Rev. Lett. 102, 217601 (2009)] have shown by first principles calculations that a monolayer BN sheet becomes macroscopically polarized in-plane when in a corrugated state. Here, we investigate the interplay of layer corrugation and in-plane polarization by atomistic lattice dynamics. We treat the coupling between static flexural modes and in-plane atomic ion displacements as an anharmonic effect, similar to the membrane effect that is at the origin of negative thermal expansion in layered crystals. We have derived analytical expressions for the corrugation-induced static in-plane strains and the optical displacements with the resulting polarization response functions. Beyond h-BN, the theory applies to transition metal dichalcogenides and dioxides. Numerical calculations show that the effects are considerably stronger for 2D h-BN than for 2H-MoS2.

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: 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: 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: 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: 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: 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 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: Starting from our previous work in which we obtained a system of coupled integrodifferential equations for acoustic sound waves and phonon density fluctuations in two-dimensional (2D) crystals, we derive here the corresponding hydrodynamic equations, and we study their consequences as a function of temperature and frequency. These phenomena encompass propagation and damping of acoustic sound waves, diffusive heat conduction, second sound, and Poiseuille heat flow, all of which are characterized by specific transport coefficients. We calculate these coefficients by means of correlation functions without using the concept of relaxation time. Numerical calculations are performed as well in order to show the temperature dependence of the transport coefficients and of the thermal conductivity. As a consequence of thermal tension, mechanical and thermal phenomena are coupled. We calculate the dynamic susceptibilities for displacement and temperature fluctuations and study their resonances. Due to the thermomechanical coupling, the thermal resonances such as the Landau-Placzek peak and the second-sound doublet appear in the displacement susceptibility, and conversely the acoustic sound wave doublet appears in the temperature susceptibility, Our analytical results not only apply to graphene, but they are also valid for arbitrary 2D crystals with hexagonal symmetry, such as 2D hexagonal boron nitride, 2H-transition-metal dichalcogenides, and oxides.