Abstract: Using the tight-binding formalism we show that in the recently experimentally realized AA-stacked graphite in essence two types of massless relativistic Dirac particles are present with a different effective speed of light. We also investigate how the electronic structure evolves from a single graphene sheet into AA-stacked graphite. It is shown that in contrast to AB-stacked graphene layers, the spectrum of AA-stacked graphene layers can be considered as a superposition of single-layer spectra and only particles with a linear spectrum at the Fermi energy around the K point are present. From the evolution of the band overlap we show that 6 multilayers of AA-stacked graphene already behave as AA-stacked graphite. The evolution of the effective speeds of light of the Dirac particles to their bulk values shows exactly the same behavior. The tight-binding parameters we use to describe AA-stacked graphite and multilayers of graphene are obtained by ab initio calculations.

Abstract: Thr structural, electronic, and magnetic properties of the highly mismatched perovskite oxides, Th(0.35)A(0.65)MnO(3), where Ais for the alkaline earth divalent cations (Ca, Ba, Sr), which are all characterized by the same large tolerance factor (t=0.934), have been investigated by using electron microscopy, electrical resistivity, magnetic susceptibility, and magnetization. It is clearly established that a transition from ferromagnetic metallic towards spin-glass insulator samples is induced as the A-site cationic size mismatch is increased. Moreover, the magnetoresistance (MR) properties of these manganites are strongly reduced for the spin-glass insulators, demonstrating that the A-sire cationic disorder is detrimental for the colossal MR properties. Based on these results, a new electronic and magnetic diagram is established that shows that the A-site disorder, rather than the A-site average cationic size (or t) is the relevant factor for generating spin-glass insulating manganites. [S0163-1829(99)01746-4].

Abstract: The hydrogen interstitial and the substitutional AlZn, GaZn, and InZn are all shallow donors in ZnO and lead to n-type conductivity. Although shallow donors are expected to repel each other, we show by first-principles calculations that in ZnO these shallow donor impurities attract and form a complex, leading to a donor level deep in the band gap. This puts a limit on the n-type conductivity of (Al,Ga,In)-doped ZnO in the presence of hydrogen.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

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.