Showing 1–22 of 22 results for author: Bonini, N

Spin-Flop Ordering from Frustrated Ferro- and Antiferromagnetic Interactions: A Combined Theoretical and Experimental Study of a $\mathrm{Mn}/\mathrm{Fe}(100)$ Monolayer

Authors: C. Grazioli, Dario Alfè, S. R. Krishnakumar, Subhra Sen Gupta, M. Veronese, S. Turchini, Nicola Bonini, Andrea Dal Corso, D. D. Sarma, Stefano Baroni, C. Carbone

Abstract: The occurrence of a noncollinear magnetic structure at a Mn monolayer grown epitaxially on Fe(100) is predicted theoretically, using spinor density-functional theory, and observed experimentally, using x-ray magnetic circular dichroism (XMCD) and linear dichroism (XMLD) spectroscopies. The combined use of XMCD and XMLD at the Mn-absorption edge allows us to assess the existence of ferromagnetic an… ▽ More The occurrence of a noncollinear magnetic structure at a Mn monolayer grown epitaxially on Fe(100) is predicted theoretically, using spinor density-functional theory, and observed experimentally, using x-ray magnetic circular dichroism (XMCD) and linear dichroism (XMLD) spectroscopies. The combined use of XMCD and XMLD at the Mn-absorption edge allows us to assess the existence of ferromagnetic and antiferromagnetic order at the interface, and also to determine the moment orientations with element specificity. The experimental results thus obtained are in excellent agreement with the magnetic structure determined theoretically. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Electron-phonon physics from first principles using the EPW code

Authors: Hyungjun Lee, Samuel Poncé, Kyle Bushick, Samad Hajinazar, Jon Lafuente-Bartolome, Joshua Leveillee, Chao Lian, Francesco Macheda, Hari Paudyal, Weng Hong Sio, Marios Zacharias, Xiao Zhang, Nicola Bonini, Emmanouil Kioupakis, Elena R. Margine, Feliciano Giustino

Abstract: EPW is an open-source software for $\textit{ab initio}$ calculations of electron-phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally-localized Wannier functions to efficiently compute electron-phonon coupling matrix elements on ultra-fine Brillouin zone grids. This data is employed for predictive calculations of temperature-d… ▽ More EPW is an open-source software for $\textit{ab initio}$ calculations of electron-phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally-localized Wannier functions to efficiently compute electron-phonon coupling matrix elements on ultra-fine Brillouin zone grids. This data is employed for predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant new developments in the code that occurred during the period 2016-2022, namely: a transport module for the calculation of charge carrier mobility and conductivity under electric and magnetic fields within the $\textit{ab initio}$ Boltzmann transport equation; a superconductivity module for the calculation of critical temperature and gap structure in phonon-mediated superconductors within the $\textit{ab initio}$ anisotropic multi-band Eliashberg theory; an optics module for calculations of phonon-assisted indirect transitions; a module for the calculation of small and large polarons without supercells using the $\textit{ab initio}$ polaron equations; and a module for calculating electron-phonon couplings, band structure renormalization, and temperature-dependent optical spectra using the special displacement method. For each capability, we outline the methodology and implementation, and provide example calculations. We describe recent code refactoring to prepare EPW for exascale architectures, we discuss efficient parallelization strategies, and report on extreme parallel scaling tests. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 61 pages, 9 figures

Exploring the Effect of the Number of Hydrogen Atoms on the Properties of Lanthanide Hydrides by DMFT

Authors: Yao Wei, Elena Chachkarova, Evgeny Plekhanov, Nicola Bonini, Cedric Weber

Abstract: Lanthanide hydrogen-rich materials have long been considered as one of the candidates with high-temperature superconducting properties in condensed matter physics, and have attracted great interest. Attempts to investigate the effects of different compositions of lanthanide hydrogen-rich materials are ongoing, with predictions and experimental studies in recent years having shown that substances s… ▽ More Lanthanide hydrogen-rich materials have long been considered as one of the candidates with high-temperature superconducting properties in condensed matter physics, and have attracted great interest. Attempts to investigate the effects of different compositions of lanthanide hydrogen-rich materials are ongoing, with predictions and experimental studies in recent years having shown that substances such as LaH 10 , CeH 9 , and LaH 16 exhibit extremely high superconducting temperatures between 150-250 GPa. In particular, researchers have noted that in those materials an increase in the f character at the Fermi level leads to an increase in the superconducting temperature. Here, we further elaborate on the effect of the ratios of lanthanide to hydrogen in these substances with the aim to bring more clarity to the study of superhydrides in these extreme cases by comparing a variety of lanthanide hydrogen-rich materials with different ratios using the DMFT method, and provide ideas for later structural predictions and material property studies. △ Less

Submitted 30 March, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 7pages, 5figures

High temperature superconductivity in the lanthanide hydrides at extreme pressures

Authors: Yao Wei, Francesco Macheda, Zelong Zhao, Terence Tse, Evgeny Plekhanov, Nicola Bonini, Cedric Weber

Abstract: Hydrogen-rich superhydrides are promising high-Tc superconductors, with superconductivity experimentally observed near room temperature, as shown in recently discovered lanthanide superhydrides at very high pressures, e.g. LaH10 at 170 GPa and CeH9 at 150 GPa. Superconductivity is believed to be closely related with the high vibrational modes of the bound hydrogen ions. Here we study the limit of… ▽ More Hydrogen-rich superhydrides are promising high-Tc superconductors, with superconductivity experimentally observed near room temperature, as shown in recently discovered lanthanide superhydrides at very high pressures, e.g. LaH10 at 170 GPa and CeH9 at 150 GPa. Superconductivity is believed to be closely related with the high vibrational modes of the bound hydrogen ions. Here we study the limit of extreme pressures (above 200 GPa) where lanthanide hydrides with large hydrogen content have been reported. We focus on LaH16 and CeH16, two prototype candidates for achieving a large electronic contribution from hydrogen in the electron-phonon coupling. In this work, we propose a first principles calculation platform with the inclusion of many-body corrections to evaluate the detailed physical properties of the Ce-H and La-H systems and to understand the structure, stability and superconductivity of these systems at ultra-high pressure. We provide a practical approach to further investigate conventional superconductivity in hydrogen rich superhydrides. We report that density functional theory provides accurate structure and phonon frequencies, but many-body corrections lead to an increase of the critical temperature, that is associated with spectral weight transfer of the f-states. △ Less

Submitted 13 January, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 8 pages, 5 figures. arXiv admin note: substantial text overlap with arXiv:2107.12316

Theory and computation of Hall scattering factor in graphene

Authors: Francesco Macheda, Samuel Poncé, Feliciano Giustino, Nicola Bonini

Abstract: The Hall scattering factor, $r$, is a key quantity for establishing carrier concentration and drift mobility from Hall measurements; in experiments it is usually assumed to be 1. In this paper we use a combination of analytical and \textit{ab initio} modelling to determine $r$ in graphene. While at high carrier densities $r \approx 1$ in a wide temperature range, at low doping the temperature depe… ▽ More The Hall scattering factor, $r$, is a key quantity for establishing carrier concentration and drift mobility from Hall measurements; in experiments it is usually assumed to be 1. In this paper we use a combination of analytical and \textit{ab initio} modelling to determine $r$ in graphene. While at high carrier densities $r \approx 1$ in a wide temperature range, at low doping the temperature dependence of $r$ is very strong with values as high as 4 below 300~K. These high values are due to the linear bands around the Dirac cone and the carrier scattering rates due to acoustic phonons. At higher temperatures $r$ can instead become as low as $0.5$ due to the contribution of both holes and electrons and the role of optical phonons. Finally, we provide a simple analytical model to compute accurately $r$ in graphene in a wide range of temperatures and carrier densities. △ Less

Submitted 21 November, 2021; originally announced November 2021.

Computational Materials Discovery for Lanthanide Hydrides at high pressure: predicting High Temperature superconductivity

Authors: Evgeny Plekhanov, Zelong Zhao, Francesco Macheda, Yao Wei, Nicola Bonini, Cedric Weber

Abstract: Hydrogen-rich superhydrides are believed to be very promising high-T$_c$ superconductors, with experimentally observed critical temperatures near room temperature, as shown in recently discovered lanthanide superhydrides at very high pressures, e.g. LaH$_{10}$ at 170 GPa and CeH$_9$ at 150 GPa. With the motivation of discovering new hydrogen-rich high-T$_c$ superconductors at lowest possible press… ▽ More Hydrogen-rich superhydrides are believed to be very promising high-T$_c$ superconductors, with experimentally observed critical temperatures near room temperature, as shown in recently discovered lanthanide superhydrides at very high pressures, e.g. LaH$_{10}$ at 170 GPa and CeH$_9$ at 150 GPa. With the motivation of discovering new hydrogen-rich high-T$_c$ superconductors at lowest possible pressure, quantitative theoretical predictions are needed. In these promising compounds, superconductivity is mediated by the highly energetic lattice vibrations associated with hydrogen and their interplay with the electronic structure, requiring fine descriptions of the electronic properties, notoriously challenging for correlated $f$ systems. In this work, we propose a first-principles calculation platform with the inclusion of many-body corrections to evaluate the detailed physical properties of the Ce-H system and to understand the structure, stability and superconductivity of CeH$_9$ at high pressure. We report how the prediction of T$_c$ is affected by the hierarchy of many-body corrections, and obtain a compelling increase of T$_c$ at the highest level of theory, which goes in the direction of experimental observations. Our findings shed a significant light on the search for superhydrides in close similarity with atomic hydrogen within a feasible pressure range. We provide a practical platform to further investigate and understand conventional superconductivity in hydrogen rich superhydrides. △ Less

Submitted 7 September, 2021; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 7 pages, 7 figures

First-principles predictions of Hall and drift mobilities in semiconductors

Authors: Samuel Ponce, Francesco Macheda, Elena Roxana Margine, Nicola Marzari, Nicola Bonini, Feliciano Giustino

Abstract: Significant progress on parameter-free calculations of carrier mobilities in real materials has been made during the past decade; however, the role of various approximations remains unclear and a unified methodology is lacking. Here, we present and analyse a comprehensive and efficient approach to compute the intrinsic, phonon-limited drift and Hall carrier mobilities of semiconductors, within the… ▽ More Significant progress on parameter-free calculations of carrier mobilities in real materials has been made during the past decade; however, the role of various approximations remains unclear and a unified methodology is lacking. Here, we present and analyse a comprehensive and efficient approach to compute the intrinsic, phonon-limited drift and Hall carrier mobilities of semiconductors, within the framework of the first-principles Boltzmann transport equation. The methodology exploits a novel approach for estimating quadrupole tensors and including them in the electron-phonon interactions, and capitalises on a rigorous and efficient procedure for numerical convergence. The accuracy reached in this work allows to assess common approximations, including the role of exchange and correlation functionals, spin-orbit coupling, pseudopotentials, Wannier interpolation, Brillouin-zone sampling, dipole and quadrupole corrections, and the relaxation-time approximation. A detailed analysis is showcased on ten prototypical semiconductors, namely diamond, silicon, GaAs, 3C-SiC, AlP, GaP, c-BN, AlAs, AlSb, and SrO. By comparing this extensive dataset with available experimental data, we explore the intrinsic nature of phonon-limited carrier transport and magnetotransport phenomena in these compounds. We find that the most accurate calculations predict Hall mobilities up to a factor of two larger than experimental data; this could point to promising experimental improvements in the samples quality, or to the limitations of density-function theory in predicting the carrier effective masses and overscreening the electron-phonon matrix elements. By setting tight standards for reliability and reproducibility, the present work aims to facilitate validation and verification of data and software towards predictive calculations of transport phenomena in semiconductors. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 40 pages and 26 figures

Journal ref: Phys. Rev. Research 3, 043022 (2021)

Calculating DMFT forces in ab-initio ultrasoft pseudopotential formalism

Authors: Evgeny Plekhanov, Nicola Bonini, Cedric Weber

Abstract: In this paper, we show how to calculate analytical atomic forces within self-consistent density functional theory + dynamical mean-field theory (DFT+DMFT) approach in the case when ultra-soft or norm-conserving pseudopotentials are used. We show how to treat the non-local projection terms arising within the pseudopotential formalism and circumvent the problem of non-orthogonality of the Kohn-Sham… ▽ More In this paper, we show how to calculate analytical atomic forces within self-consistent density functional theory + dynamical mean-field theory (DFT+DMFT) approach in the case when ultra-soft or norm-conserving pseudopotentials are used. We show how to treat the non-local projection terms arising within the pseudopotential formalism and circumvent the problem of non-orthogonality of the Kohn-Sham eigenvectors. Our approach is, in principle, independent of the DMFT solver employed, and here was tested with the Hubbard I solver. We benchmark our formalism by comparing against the forces calculated in Ce$_{2}$O$_{3}$ and PrO$_2$ by numerical differentiation of the total free energy, as well as by comparing the energy profiles against the numerically integrated analytical forces. △ Less

Submitted 16 December, 2021; v1 submitted 9 February, 2021; originally announced February 2021.

Comments: 12 pages, 3 figures

Journal ref: Phys. Rev. B 104 235131 (2021)

Thermomechanical properties of honeycomb lattices from internal-coordinates potentials: the case of graphene and hexagonal boron nitride

Authors: Francesco Libbi, Nicola Bonini, Nicola Marzari

Abstract: Lattice dynamics in low-dimensional materials and, in particular, the quadratic behaviour of the flexural acoustic modes play a fundamental role in their thermomechanical properties. A first-principles evaluation of these can be very demanding, and can be affected by numerical noise that breaks translational or rotational invariance. In order to overcome these challenges, we study the Gartstein in… ▽ More Lattice dynamics in low-dimensional materials and, in particular, the quadratic behaviour of the flexural acoustic modes play a fundamental role in their thermomechanical properties. A first-principles evaluation of these can be very demanding, and can be affected by numerical noise that breaks translational or rotational invariance. In order to overcome these challenges, we study the Gartstein internal-coordinate potential and tune its 13 parameters on the first-principles interatomic force constants for graphene. We show that the resulting potential not only reproduces very well the phonon dispersions of graphene, but also those of carbon nanotubes of any diameter and chirality. The addition of a cubic term allows also to reproduce the dominant anharmonic terms, leading to a very good estimate of the lattice thermal conductivity. Finally, this potential form works very well also for boron nitride, provided it is fitted on the short-range (analytical) part of the interatomic force constants, and augmented thereafter with the long-range dielectric contribution. This consideration underscores how potentials based on short-ranged descriptors should be fit, in polar materials, to the short-range part of the first-principles interactions, and complemented by long-range analytical dielectric models parametrized on the same first-principles calculations. △ Less

Submitted 26 October, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

First-principles study of electronic transport and structural properties of Cu$_{12}$Sb$_4$S$_{13}$ in its high-temperature phase

Authors: Cono di Paola, Francesco Macheda, Savio Laricchia, Cedric Weber, Nicola Bonini

Abstract: We present an ab initio study of the structural and electronic transport properties of tetrahedrite, Cu$_{12}$Sb$_4$S$_{13}$, in its high-temperature phase. We show how this complex compound can be seen as the outcome of an ordered arrangement of S-vacancies in a semiconducting fematinite-like structure (Cu3SbS4). Our calculations confirm that the S-vacancies are the natural doping mechanism in th… ▽ More We present an ab initio study of the structural and electronic transport properties of tetrahedrite, Cu$_{12}$Sb$_4$S$_{13}$, in its high-temperature phase. We show how this complex compound can be seen as the outcome of an ordered arrangement of S-vacancies in a semiconducting fematinite-like structure (Cu3SbS4). Our calculations confirm that the S-vacancies are the natural doping mechanism in this thermoelectric compound and reveal a similar local chemical environment around crystallographically inequivalent Cu atoms, shedding light on the debate on XPS measurements in this compound. To access the electrical transport properties as a function of temperature we use the Kubo-Greenwood formula applied to snapshots of first-principles molecular dynamics simulations. This approach is essential to effectively account for the interaction between electrons and lattice vibrations in such a complex crystal structure where a strong anharmonicity plays a key role in stabilising the high-temperature phase. Our results show that the Seebeck coeffcient is in good agreement with experiments and the phonon-limited electrical resistivity displays a temperature trend that compares well with a wide range of experimental data. The predicted lower bound for the resistivity turns out to be remarkably low for a pristine mineral in the Cu-Sb-S system but not too far from the lowest experimental data reported in literature. The Lorenz number turns out to be substantially lower than what expected from the free-electron value in the Wiedemann-Franz law, thus providing an accurate way to estimate the electronic and lattice contributions to the thermal conductivity in experiments, of great significance in this very low thermal conductivity crystalline material. △ Less

Submitted 6 July, 2020; v1 submitted 3 July, 2020; originally announced July 2020.

Magneto-transport phenomena in p-doped diamond from first principles

Authors: Francesco Macheda, Nicola Bonini

Abstract: We present a first-principles study of the magnetotransport phenomena in p-doped diamond via the exact solution of the linearized Boltzmann transport equation, in which the materials' parameters, including electron-phonon and phonon-phonon interactions, are obtained from density functional theory. This approach gives results in very good agreement with experimental data for Hall and drift mobiliti… ▽ More We present a first-principles study of the magnetotransport phenomena in p-doped diamond via the exact solution of the linearized Boltzmann transport equation, in which the materials' parameters, including electron-phonon and phonon-phonon interactions, are obtained from density functional theory. This approach gives results in very good agreement with experimental data for Hall and drift mobilities, low- and high-field magnetoresistance and Seebeck coefficient, including the phonon-drag effect, in a range of temperatures and carrier concentrations. In particular, our results provide a detailed characterisation of the exceptionally high values for mobility and Seebeck coefficient, and predict a large magnetic field driven enhancement of the Seebeck coefficient, of up to 30% in a magnetic field of 40 kOe already at room temperature. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Journal ref: Phys. Rev. B 98, 201201 (2018)

Phonon-limited resistivity of graphene by first-principle calculations: electron-phonon interactions, strain-induced gauge field and Boltzmann equation

Authors: Thibault Sohier, Matteo Calandra, Cheol-Hwan Park, Nicola Bonini, Nicola Marzari, Francesco Mauri

Abstract: Electron-phonon coupling in graphene is extensively modeled and simulated from first principles. We find that using an accurate model for the polarizations of the acoustic phonon modes is crucial to obtain correct numerical results. The interactions between electrons and acoustic phonon modes, the gauge field and deformation potential, are calculated at the DFT level in the framework of linear res… ▽ More Electron-phonon coupling in graphene is extensively modeled and simulated from first principles. We find that using an accurate model for the polarizations of the acoustic phonon modes is crucial to obtain correct numerical results. The interactions between electrons and acoustic phonon modes, the gauge field and deformation potential, are calculated at the DFT level in the framework of linear response. The zero-momentum limit of acoustic phonons is interpreted as a strain pattern, allowing the calculation of the acoustic gauge field parameter in the GW approximation. The role of electronic screening on the electron-phonon matrix elements is investigated. We then solve the Boltzmann equation semi-analytically in graphene, including both acoustic and optical phonon scattering. We show that, in the Bloch-Grüneisen and equipartition regimes, the electronic transport is mainly ruled by the unscreened acoustic gauge field, while the contribution due to the deformation potential is negligible and strongly screened. By comparing with experimental data, we show that the contribution of acoustic phonons to resistivity is doping- and substrate-independent. The DFT+GW approach underestimates this contribution to resistivity by about 30 %. Above 270K, the calculated resistivity underestimates the experimental one more severely, the underestimation being larger at lower doping. We show that, beside remote phonon scattering, a possible explanation for this disagreement is the electron-electron interaction that strongly renormalizes the coupling to intrinsic optical-phonon modes. Finally, after discussing the validity of the Matthiessen rule in graphene, we derive simplified analytical solutions of the Boltzmann equation to extract the coupling to acoustic phonons, related to the strain-induced gauge field, directly from experimental data. △ Less

Submitted 10 September, 2014; v1 submitted 3 July, 2014; originally announced July 2014.

Journal ref: Physical Review B 90, 125414 (2014)

Electron-Phonon Interactions and the Intrinsic Electrical Resistivity of Graphene

Authors: Cheol-Hwan Park, Nicola Bonini, Thibault Sohier, Georgy Samsonidze, Boris Kozinsky, Matteo Calandra, Francesco Mauri, Nicola Marzari

Abstract: We present a first-principles study of the temperature- and density-dependent intrinsic electrical resistivity of graphene. We use density-functional theory and density-functional perturbation theory together with very accurate Wannier interpolations to compute all electronic and vibrational properties and electron-phonon coupling matrix elements; the phonon-limited resistivity is then calculated… ▽ More We present a first-principles study of the temperature- and density-dependent intrinsic electrical resistivity of graphene. We use density-functional theory and density-functional perturbation theory together with very accurate Wannier interpolations to compute all electronic and vibrational properties and electron-phonon coupling matrix elements; the phonon-limited resistivity is then calculated within a Boltzmann-transport approach. An effective tight-binding model, validated against first-principles results, is also used to study the role of electron-electron interactions at the level of many-body perturbation theory. The results found are in excellent agreement with recent experimental data on graphene samples at high carrier densities and elucidate the role of the different phonon modes in limiting electron mobility. Moreover, we find that the resistivity arising from scattering with transverse acoustic phonons is 2.5 times higher than that from longitudinal acoustic phonons. Last, high-energy, optical, and zone-boundary phonons contribute as much as acoustic phonons to the intrinsic electrical resistivity even at room temperature and become dominant at higher temperatures. △ Less

Submitted 18 March, 2014; originally announced March 2014.

Comments: 7 pages 5 figures

Journal ref: Nano Lett. 14, 1113-1119 (2014)

Temperature evolution of infrared- and Raman-active phonons in graphite

Authors: P. Giura, N. Bonini, G. Creff, J. B. Brubach, P. Roy, M. Lazzeri

Abstract: We perform a comparative experimental and theoretical study of the temperature dependence up to 700 K of the frequency and linewidths of the graphite E1u and E2g optical phonons (~1590 and 1580 cm-1) by infra-red (IR) and Raman spectroscopy. Despite their similar character, the temperature dependence of the two modes is quite different, being, e.g., the frequency shift of the IR-active E1u mode is… ▽ More We perform a comparative experimental and theoretical study of the temperature dependence up to 700 K of the frequency and linewidths of the graphite E1u and E2g optical phonons (~1590 and 1580 cm-1) by infra-red (IR) and Raman spectroscopy. Despite their similar character, the temperature dependence of the two modes is quite different, being, e.g., the frequency shift of the IR-active E1u mode is almost twice as big as that of the Raman active E2g mode. Ab initio calculations of the anharmonic properties are in remarkable agreement with measurements and explain the observed behavior. △ Less

Submitted 20 September, 2012; v1 submitted 26 July, 2012; originally announced July 2012.

Comments: 5 pages, 4 figures

Journal ref: Physical Review B 86, 121404(R) (2012)

The Shear Mode of Multi-Layer Graphene

Authors: P. H. Tan, W. P. Han, W. J. Zhao, Z. H. Wu, K. Chang, H. Wang, Y. F. Wang, N. Bonini, N. Marzari, G. Savini, A. Lombardo, A. C. Ferrari

Abstract: We uncover the interlayer shear mode of multi-layer graphene samples, ranging from bilayer-graphene (BLG) to bulk graphite, and show that the corresponding Raman peak measures the interlayer coupling. This peak scales from~43cm-1 in bulk graphite to~31cm-1 in BLG. Its low energy makes it a probe of near-Dirac point quasi-particles, with a Breit-Wigner-Fano lineshape due to resonance with electroni… ▽ More We uncover the interlayer shear mode of multi-layer graphene samples, ranging from bilayer-graphene (BLG) to bulk graphite, and show that the corresponding Raman peak measures the interlayer coupling. This peak scales from~43cm-1 in bulk graphite to~31cm-1 in BLG. Its low energy makes it a probe of near-Dirac point quasi-particles, with a Breit-Wigner-Fano lineshape due to resonance with electronic transitions. Similar shear modes are expected in all layered materials, providing a direct probe of interlayer interactions △ Less

Submitted 6 June, 2011; originally announced June 2011.

Journal ref: Nature Materials 11, 294-300 (2012)

Dominant Phonon Wavevectors and Strain-induced Splitting of the 2D Graphene Raman Mode

Authors: Rohit Narula, Nicola Bonini, Nicola Marzari, Stephanie Reich

Abstract: The dominant phonon wavevectors $q^{*}$ probed by the 2D Raman mode of graphene are highly anisotropic and rotate with the orientation of the polarizer:analyzer direction relative to the lattice. The corresponding electronic transitions connect the electronic equibandgap contours where the product of the ingoing and outgoing optical matrix elements is strongest, showing a finite component along… ▽ More The dominant phonon wavevectors $q^{*}$ probed by the 2D Raman mode of graphene are highly anisotropic and rotate with the orientation of the polarizer:analyzer direction relative to the lattice. The corresponding electronic transitions connect the electronic equibandgap contours where the product of the ingoing and outgoing optical matrix elements is strongest, showing a finite component along $\bm{K}-\bmΓ$ that sensitively determines $q^{*}$. We revoke the notion of 'inner' and 'outer' processes. Our findings explain the splitting of the 2D mode of graphene under uniaxial tensile strain. The splitting originates from a strain-induced distortion of the phonon dispersion; changes in the electronic band structure and resonance conditions are negligeable for the 2D Raman spectrum. △ Less

Submitted 2 December, 2011; v1 submitted 25 February, 2011; originally announced February 2011.

Comments: 4 pages, 6 figures

Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials

Authors: P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, m. Cococcioni, I. Dabo, A. Dal Corso, S. Fabris, G. Fratesi, S. de Gironcoli, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello , et al. (8 additional authors not shown)

Abstract: Quantum ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). Quantum ESPRESSO stands for "opEn Source Package for Research in Electronic Structure, Simulation, and Optimization". It is freely available to rese… ▽ More Quantum ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). Quantum ESPRESSO stands for "opEn Source Package for Research in Electronic Structure, Simulation, and Optimization". It is freely available to researchers around the world under the terms of the GNU General Public License. Quantum ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively-parallel architectures, and a great effort being devoted to user friendliness. Quantum ESPRESSO is evolving towards a distribution of independent and inter-operable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes. △ Less

Submitted 26 July, 2009; v1 submitted 14 June, 2009; originally announced June 2009.

Comments: 36 pages, 5 figures, resubmitted to J.Phys.: Condens. Matter

Journal ref: J. Phys.: Condens. Matter 21, 395502 (2009)

Uniaxial Strain in Graphene by Raman Spectroscopy: G peak splitting, Gruneisen Parameters and Sample Orientation

Authors: T. M. G. Mohiuddin, A. Lombardo, R. R. Nair, A. Bonetti, G. Savini, R. Jalil, N. Bonini, D. M. Basko, C. Galiotis, N. Marzari, K. S. Novoselov, A. K. Geim, A. C. Ferrari

Abstract: Graphene is the two-dimensional building block for carbon allotropes of every other dimensionality. Since its experimental discovery, graphene continues to attract enormous interest, in particular as a new kind of matter, in which electron transport is governed by a Dirac-like wave equation, and as a model system for studying electronic and phonon properties of other, more complex, graphitic mat… ▽ More Graphene is the two-dimensional building block for carbon allotropes of every other dimensionality. Since its experimental discovery, graphene continues to attract enormous interest, in particular as a new kind of matter, in which electron transport is governed by a Dirac-like wave equation, and as a model system for studying electronic and phonon properties of other, more complex, graphitic materials[1-4]. Here, we uncover the constitutive relation of graphene and probe new physics of its optical phonons, by studying its Raman spectrum as a function of uniaxial strain. We find that the doubly degenerate E2g optical mode splits in two components, one polarized along the strain and the other perpendicular to it. This leads to the splitting of the G peak into two bands, which we call G+ and G-, by analogy with the effect of curvature on the nanotube G peak[5-7]. Both peaks red shift with increasing strain, and their splitting increases, in excellent agreement with first-principles calculations. Their relative intensities are found to depend on light polarization, which provides a useful tool to probe the graphene crystallographic orientation with respect to the strain. The singly degenerate 2D and 2D' bands also red shift, but do not split for small strains. We study the Gruneisen parameters for the phonons responsible for the G, D and D' peaks. These can be used to measure the amount of uniaxial or biaxial strain, providing a fundamental tool for nanoelectronics, where strain monitoring is of paramount importance[8, 9] △ Less

Submitted 8 December, 2008; originally announced December 2008.

Journal ref: Phys. Rev. B, 79, 205433 (2009)

Bulk Aluminum at High Pressure: A First-Principles Study

Authors: Michael J. Tambe, Nicola Bonini, Nicola Marzari

Abstract: The behavior of metals at high pressure is of great importance to the fields of shock physics, geophysics, astrophysics, and nuclear materials. In order to further understand the properties of metals at high pressures we studied the equation of state of aluminum using first-principles techniques up to 2500 GPa, pressures within reach of the planned L.L.N.L. National Ignition Facility. Our simula… ▽ More The behavior of metals at high pressure is of great importance to the fields of shock physics, geophysics, astrophysics, and nuclear materials. In order to further understand the properties of metals at high pressures we studied the equation of state of aluminum using first-principles techniques up to 2500 GPa, pressures within reach of the planned L.L.N.L. National Ignition Facility. Our simulations use density-functional theory and density-functional perturbation theory in the generalized gradient approximation at 0K. We found core overlaps to become relevant beyond pressures of 1200 GPa. The equations of state for three phases (fcc, bcc, and hcp) were calculated predicting the fcc-hcp, fcc-bcc, and hcp-bcc transitions to occur at 215 GPa, 307 GPa, and 435 GPa respectively. From the phonon dispersions at increasing pressure, we predict a softening of the lowest transverse acoustic vibrational mode along the [110] direction, which corresponds to a Born instability of the fcc phase at 725 GPa. △ Less

Submitted 17 April, 2008; v1 submitted 26 October, 2007; originally announced October 2007.

Comments: 4 pages, 5 figures, accepted to Phys. Rev. B as a Brief Report. This version has update many figures. Moreover we provided updated and more accurate numbers based on further in-depth analyses of potential computational errors

Phonon anharmonicities in graphite and graphene

Authors: Nicola Bonini, Michele Lazzeri, Nicola Marzari, Francesco Mauri

Abstract: We determine from first-principles the finite-temperature properties--linewidths, line shifts, and lifetimes--of the key vibrational modes that dominate inelastic losses in graphitic materials. In graphite, the phonon linewidth of the Raman-active E2g mode is found to decrease with temperature; such anomalous behavior is driven entirely by electron-phonon interactions, and does not appear in the… ▽ More We determine from first-principles the finite-temperature properties--linewidths, line shifts, and lifetimes--of the key vibrational modes that dominate inelastic losses in graphitic materials. In graphite, the phonon linewidth of the Raman-active E2g mode is found to decrease with temperature; such anomalous behavior is driven entirely by electron-phonon interactions, and does not appear in the nearly-degenerate infrared-active E1u mode. In graphene, the phonon anharmonic lifetimes and decay channels of the A'1 mode at K dominate over E2g at G and couple strongly with acoustic phonons, highlighting how ballistic transport in carbon-based interconnects requires careful engineering of phonon decays and thermalization. △ Less

Submitted 25 October, 2007; v1 submitted 30 August, 2007; originally announced August 2007.

Comments: 5 pages, 4 figures; typos corrected and reference added

Journal ref: Phys. Rev. Lett. 99, 176802 (2007)

Structure and dynamics of oxygen adsorbed on Ag(100) vicinal surfaces

Authors: N. Bonini, A. Kokalj, A. Dal Corso, S. de Gironcoli, S. Baroni

Abstract: The structure and dynamics of atomic oxygen adsorbed on Ag(410) and Ag(210) surfaces have been investigated using density functional theory. Our results show that the adsorption configuration in which O adatoms decorate the upper side of the (110) steps forming O--Ag--O rows is particularly stable for both surfaces. On Ag(210), this arrangement is more stable than other configurations at all the… ▽ More The structure and dynamics of atomic oxygen adsorbed on Ag(410) and Ag(210) surfaces have been investigated using density functional theory. Our results show that the adsorption configuration in which O adatoms decorate the upper side of the (110) steps forming O--Ag--O rows is particularly stable for both surfaces. On Ag(210), this arrangement is more stable than other configurations at all the investigated coverages. On Ag(410), adsorption on the terrace and at the step edge are almost degenerate, the former being slightly preferred at low coverage while the latter is stabilized by increasing the coverage. These findings are substantiated by a comparison between the vibrational modes, calculated within density-functional perturbation theory, and the HREEL spectrum which has been recently measured in these systems. △ Less

Submitted 15 March, 2004; originally announced March 2004.

Comments: 17 pages, 7 tables, 12 figures. To appear in Phys. Rev. B

On-surface and Subsurface Adsorption of Oxygen on Stepped Ag(210) and Ag(410) Surfaces

Authors: A. Kokalj, N. Bonini, A. Dal Corso, S. de Gironcoli, S. Baroni

Abstract: The adsorption of atomic oxygen and its inclusion into subsurface sites on Ag(210) and Ag(410) surfaces have been investigated using density functional theory. We find that--in the absence of adatoms on the first metal layer--subsurface adsorption results in strong lattice distortion which makes it energetically unfavoured. However subsurface sites are significantly stabilised when a sufficient… ▽ More The adsorption of atomic oxygen and its inclusion into subsurface sites on Ag(210) and Ag(410) surfaces have been investigated using density functional theory. We find that--in the absence of adatoms on the first metal layer--subsurface adsorption results in strong lattice distortion which makes it energetically unfavoured. However subsurface sites are significantly stabilised when a sufficient amount of O adatoms is present on the surface. At high enough O coverage on the Ag(210) surface the mixed on-surface + subsurface O adsorption is energetically favoured with respect to the on-surface only adsorption. Instead, on the Ag(410) surface, at the coverage we have considered (3/8 ML), the existence of stable terrace sites makes the subsurface O incorporation less favourable. These findings are compatible with the results of recent HREEL experiments which have actually motivated this work. △ Less

Submitted 5 November, 2003; originally announced November 2003.

Comments: 8 pages, 4 figures and 1 table