Anharmonic properties from a generalized third order ab~initio approach: theory and applications to graphite and graphene

We have implemented a generic method, based on the 2n+1 theorem within density functional perturbation theory, to calculate the anharmonic scattering coefficients among three phonons with arbitrary wavevectors. The method is used to study the phonon broadening in graphite and graphene mono- and bi-layer. The broadening of the high-energy optical branches is highly nonuniform and presents a series of sudden steps and spikes. At finite temperature, the two linearly dispersive acoustic branches TA and LA of graphene have nonzero broadening for small wavevectors. The broadening in graphite and bi-layer graphene is, overall, very similar to the graphene one, the most remarkable feature being the broadening of the quasi acoustical ZO branch. Finally, we study the intrinsic anharmonic contribution to the thermal conductivity of the three systems, within the single mode relaxation time approximation. We find the conductance to be in good agreement with experimental data for the out-of-plane direction but to underestimate it by a factor 2 in-plane.

Ab initio variational approach for evaluating lattice thermal conductivity

We present a first-principles theoretical approach for evaluating the lattice thermal conductivity based on the exact solution of the Boltzmann transport equation. We use the variational principle and the conjugate gradient scheme, which provide us with an algorithm faster than the one previously used in literature and able to always converge to the exact solution. Three-phonon normal and umklapp collision, isotope scattering and border effects are rigorously treated in the calculation. Good agreement with experimental data for diamond is found. Moreover we show that by growing more enriched diamond samples it is possible to achieve values of thermal conductivity up to three times larger than the commonly observed in isotopically enriched diamond samples with 99.93% C12 and 0.07 C13.

Zener tunneling in the electrical transport of quasi-metallic carbon nanotubes

We study theoretically the impact of Zener tunneling on the charge-transport properties of quasi-metallic (Qm) carbon nanotubes (characterized by forbidden band gaps of few tens of meV). We also analyze the interplay between Zener tunneling and elastic scattering on defects. To this purpose we use a model based on the master equation for the density matrix, that takes into account the inter-band Zener transitions induced by the electric field (a quantum mechanical effect), the electron-defect scattering and the electron-phonon scattering. In presence of Zener tunnelling the Qm tubes support an electrical current even when the Fermi energy lies in the forbidden band gap. In absence of elastic scattering (in high quality samples), the small size of the band gap of Qm tubes enables Zener tunnelling for realistic values of the the electric field (above $sim$ 1 V/mu m). The presence of a strong elastic scattering (in low quality samples) further decreases the values of the field required to observe Zener tunnelling. Indeed, for elastic-scattering lengths of the order of 50 nm, Zener tunnelling affects the current/voltage characteristic already in the linear regime. In other words, in quasi-metallic tubes, Zener tunneling is made more visible by defects.

Hidden polymorphs drive the vitrification in B2O3

Understanding the conditions which favor crystallisation or vitrification of liquids has been a long-standing scientific problem. Another connected, and not yet well understood question is the relationship between the glassy and the various possible crystalline forms a system may adopt. In this context, B2O3 is a puzzling case of study since i) it is one of the best glass-forming systems despite an apparent lack of low-pressure polymorphism ii) it vitrifies in a glassy form abnormally different from the only known crystalline phase at ambient pressure iii) it never crystallises from the melt unless pressure is applied, an intriguing behaviour known as the crystallisation anomaly. Here, by means of ab-initio calculations, we discover the existence of novel B2O3 crystalline polymorphs with structural properties similar to the glass and formation energies comparable to the known ambient crystal. The resulting configurational degeneracy drives the system vitrification at ambient pressure. The degeneracy is lifted under pressure, unveiling the origin of the crystallisation anomaly. This work reconciles the behaviour of B2O3 with that from other glassy systems and reaffirms the role played by polymorphism in a systems ability to vitrify. Some of the predicted crystals are cage-like materials entirely made of three-fold rings, opening new perspectives for the synthesis of boron-based nanoporous materials.

Theory of double-resonant Raman spectra in graphene: intensity and line shape of defect-induced and two-phonon bands

We calculate the double resonant (DR) Raman spectrum of graphene, and determine the lines associated to both phonon-defect processes, and two-phonons ones. Phonon and electronic dispersions reproduce calculations based on density functional theory corrected with GW. Electron-light, -phonon, and -defect scattering matrix elements and the electronic linewidth are explicitly calculated. Defect-induced processes are simulated by considering different kind of idealized defects. For an excitation energy of $epsilon_L=2.4$ eV, the agreement with measurements is very good and calculations reproduce: the relative intensities among phonon-defect or among two-phonon lines; the measured small widths of the D, $D$, 2D and $2D$ lines; the line shapes; the presence of small intensity lines in the 1800, 2000 cm$^{-1}$ range. We determine how the spectra depend on the excitation energy, on the light polarization, on the electronic linewidth, on the kind of defects and on their concentration. According to the present findings, the intensity ratio between the $2D$ and 2D lines can be used to determine experimentally the electronic linewidth. The intensity ratio between the $D$ and $D$ lines depends on the kind of model defect, suggesting that this ratio could possibly be used to identify the kind of defects present in actual samples. Charged impurities outside the graphene plane provide an almost undetectable contribution to the Raman signal.

Temperature dependent phonon renormalization in metallic nanotubes

We measure the temperature dependence of the Raman spectra of metallic and semiconducting nanotubes. We show that the different trend in metallic tubes is due to phonon re-normalization induced by the variation in electronic temperature, which is modeled including non-adiabatic contributions to account for the dynamic, time dependent nature