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.

Abstract: Neutron and high resolution X-ray diffraction investigations on perfect single crystals of RbH2PO4 (RDP), a hydrogen bonded ferroelectric of KDP type are reported. The results of crystal structure analysis from diffraction data, below and above the paraelectric – ferroelectric phase transition, support a disorder – order character Of [PO4H2](-)-groups. The tetragonal symmetry of the paraelectric phase with the double well potential of the hydrogen atoms obtained by diffraction, results simply from a time-space average of orthorhombic symmetry. According to the group – subgroup relation between the tetragonal space group 142d and the orthorhombic Fdd2 a short range order of ferroelectric clusters in the tetragonal phase is observed. With decreasing temperature the ferroelectric clusters increase and the long range interaction between their local polarisation vectors leads to the formation of lamellar ferroelectric domains with alternating polarisation directions at T-C = 147 K. From the high resolution X-ray data it is concluded that below T-C the ferroelastic strain in the (a,b)-plane leads to micro-angle grain boundaries at the domain walls. The tilt angle is enhanced by an applied electric field parallel to the ferroelectric axis. The resulting dislocations at the domain walls persist in the paraelectric phase leading to a memory effect for the arrangement of twin lamellae. With increased electric field the phase transition temperature T-C is decreased.

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: 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 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: 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: Born's long wave method is used to study the elastic properties of graphite and graphene. Starting from an empirical force constant model derived from full inplane phonon dispersions of graphite [Mohr et al., Phys. Rev. B 76, 035439 (2007)] we calculate the tension coefficients of graphene. Extending the model by interplanar interactions, we calculate the elastic constants of graphite. The agreement of our theoretical values with inelastic x-ray scattering results on elastic constants of graphite [Bosak et al., Phys. Rev. B 75, 153408 (2007)] is very satisfactory.

Abstract: From the theory of many-electron states in atoms, we know that there exists a strong Coulomb repulsion, which results in the electronic term structure of atoms and is responsible for Hunds rules. By expanding the Coulomb on-site repulsion into a multipolar series, we derive this interaction and show that it is also present in solids as a correlation effect, which means that the interaction requires a multideterminant version of the Hartree-Fock method. Of particular interest is the case where this interaction couples states of localized ( f) and delocalized ( s) electrons. We show that the interaction is bilinear in the creation/annihilation operators for localized electrons and bilinear in the operators for conduction electrons. To study the coupling, we consider a simple model in the framework of an effective limited configuration interaction method with one localized f-electron and one itinerant s-electron per crystal site. The on-site multipole interaction between the f- and s-electrons is explicitly taken into account. It is shown that depending on the low-lying excitation spectrum imposed by the crystal electric field, the model can lead not only to ferromagnetism but also to a nonmagnetic state. The model is relevant for solids with localized and itinerant electron states.

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.