-
Ab-initio Van der Waals electrodynamics: polaritons and electron scattering from plasmons and phonons in BN-capped graphene
Authors:
Francesco Macheda,
Francesco Mauri,
Thibault Sohier
Abstract:
Plasmons and polar phonons are elementary electrodynamic excitations of matter. In 2d and at long wavelengths, they couple to light and act as the system polaritons. They also dictate the scattering of charged carriers. Van der Waals heterostructures offer the opportunity to couple excitations from different layers via long-range Coulomb interactions, modifying both their dispersion and their scat…
▽ More
Plasmons and polar phonons are elementary electrodynamic excitations of matter. In 2d and at long wavelengths, they couple to light and act as the system polaritons. They also dictate the scattering of charged carriers. Van der Waals heterostructures offer the opportunity to couple excitations from different layers via long-range Coulomb interactions, modifying both their dispersion and their scattering of electrons. Even when the excitations do not couple, they are still influenced by the screening from all layers, leading to complex dynamical interactions between electrons, plasmons and polar phonons. We develop an efficient ab initio model to solve the dynamical electric response of Van der Waals heterostructures, accompanied by a formalism to extract relevant spectroscopic and transport quantities. Notably, we obtain scattering rates for electrons of the heterostructure coupling remotely with electrodynamic excitations. We apply those developments to BN-capped graphene, in which polar phonons from BN couple to plasmons in graphene. We study the nature of the coupled excitations, their dispersion and their coupling to graphene's electrons. Regimes driven by either phonons or plasmons are identified, as well as a truly hybrid regime at long wavelengths. Those are studied as a function of the graphene's Fermi level and the number of BN layers. In contrast with descriptions in terms of surface-optical phonons, we find that the electron-phonon interaction stems from different modes. Moreover, the dynamical screening of the coupling between BN's LO phonons and graphene's electrons crosses over from inefficient to metal-like depending on the relative value of the phonons' frequency and the energetic onset of interband transitions. While the coupling is significant in general, the associated scattering of graphene's carriers is found to be negligible in the context of electronic transport.
△ Less
Submitted 21 March, 2024;
originally announced March 2024.
-
Theory of infrared double-resonance Raman spectrum in graphene: the role of the zone-boundary electron-phonon enhancement
Authors:
Lorenzo Graziotto,
Francesco Macheda,
Thibault Sohier,
Matteo Calandra,
Francesco Mauri
Abstract:
We theoretically investigate the double-resonance Raman spectrum of monolayer graphene down to infrared laser excitation energies. By using first-principles density functional theory calculations, we improve upon previous theoretical predictions based on conical models or tight-binding approximations, and rigorously justify the evaluation of the electron-phonon enhancement found in Ref. [Venanzi,…
▽ More
We theoretically investigate the double-resonance Raman spectrum of monolayer graphene down to infrared laser excitation energies. By using first-principles density functional theory calculations, we improve upon previous theoretical predictions based on conical models or tight-binding approximations, and rigorously justify the evaluation of the electron-phonon enhancement found in Ref. [Venanzi, T., Graziotto, L. et al., Phys. Rev. Lett. 130, 256901 (2023)]. We proceed to discuss the effects of such enhancement on the room temperature graphene resistivity, hinting towards a possible reconciliation of theoretical and experimental discrepancies.
△ Less
Submitted 19 January, 2024; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Density-functional perturbation theory for one-dimensional systems: implementation and relevance for phonons and electron-phonon interactions
Authors:
Norma Rivano,
Nicola Marzari,
Thibault Sohier
Abstract:
The electronic and vibrational properties and electron-phonon couplings of one-dimensional materials will be key to many prospective applications in nanotechnology. Dimensionality strongly affects these properties and has to be correctly accounted for in first-principles calculations. Here we develop and implement a formulation of density-functional and density-functional perturbation theory that…
▽ More
The electronic and vibrational properties and electron-phonon couplings of one-dimensional materials will be key to many prospective applications in nanotechnology. Dimensionality strongly affects these properties and has to be correctly accounted for in first-principles calculations. Here we develop and implement a formulation of density-functional and density-functional perturbation theory that is tailored for one-dimensional systems. A key ingredient is the inclusion of a Coulomb cutoff, a reciprocal-space technique designed to correct for the spurious interactions between periodic images in periodic-boundary conditions. This restores the proper one-dimensional open-boundary conditions, letting the true response of the isolated one-dimensional system emerge. In addition to total energies, forces and stress tensors, phonons and electron-phonon interactions are also properly accounted for. We demonstrate the relevance of the present method on a portfolio of realistic systems: BN atomic chains, BN armchair nanotubes, and GaAs nanowires. Notably, we highlight the critical role of the Coulomb cutoff by studying previously inaccessible polar-optical phonons and Frohlich electron-phonon couplings. We also develop and apply analytical models to support the physical insights derived from the calculations and we discuss their consequences on electronic lifetimes. The present work unlocks the possibility to accurately simulate the linear response properties of one-dimensional systems, sheds light on the transition between dimensionalities and paves the way for further studies in several fields, including charge transport, optical coupling and polaritronics.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
Electron-phonon interaction and phonons in 2d doped semiconductors
Authors:
Francesco Macheda,
Thibault Sohier,
Paolo Barone,
Francesco Mauri
Abstract:
Electron-phonon interaction and phonon frequencies of doped polar semiconductors are sensitive to long-range Coulomb forces and can be strongly affected by screening effects of free carriers, the latter changing significantly when approaching the two-dimensional limit. We tackle this problem within a linear-response dielectric-matrix formalism, where screening effects can be properly taken into ac…
▽ More
Electron-phonon interaction and phonon frequencies of doped polar semiconductors are sensitive to long-range Coulomb forces and can be strongly affected by screening effects of free carriers, the latter changing significantly when approaching the two-dimensional limit. We tackle this problem within a linear-response dielectric-matrix formalism, where screening effects can be properly taken into account by generalized effective charge functions and the inverse scalar dielectric function, allowing for controlled approximations in relevant limits. We propose complementary computational methods to evaluate from first principles both effective charges -- encompassing all multipolar components beyond dynamical dipoles and quadrupoles -- and the static dielectric function of doped two-dimensional semiconductors, and provide analytical expressions for the long-range part of the dynamical matrix and the electron-phonon interaction in the long-wavelength limit. As a representative example, we apply our approach to study the impact of doping in disproportionated graphene, showing that optical Fröhlich and acoustic piezoelectric couplings, as well as the slope of optical longitudinal modes, are strongly reduced, with a potential impact on the electronic/intrinsic scattering rates and related transport properties.
△ Less
Submitted 23 December, 2022;
originally announced December 2022.
-
Infrared-active phonons in one-dimensional materials and their spectroscopic signatures
Authors:
Norma Rivano,
Nicola Marzari,
Thibault Sohier
Abstract:
Dimensionality provides a clear fingerprint on the dispersion of infrared-active, polar-optical phonons. For these phonons, the local dipoles parametrized by the Born effective charges drive the LO-TO splitting of bulk materials; this splitting actually breaks down in two-dimensional materials. Here, we extend the existing theory to the one-dimensional (1D) case. Combining an analytical model with…
▽ More
Dimensionality provides a clear fingerprint on the dispersion of infrared-active, polar-optical phonons. For these phonons, the local dipoles parametrized by the Born effective charges drive the LO-TO splitting of bulk materials; this splitting actually breaks down in two-dimensional materials. Here, we extend the existing theory to the one-dimensional (1D) case. Combining an analytical model with the implementation of density-functional perturbation theory in 1D boundary conditions, we show that the dielectric splitting in the dispersion relations collapses logarithmically at the zone center. The dielectric properties and the radius of the 1D materials are linked by the present work to these red shifts, opening novel IR and Raman characterization avenues.
△ Less
Submitted 21 August, 2022;
originally announced August 2022.
-
The impact of valley profile on the mobility and Kerr rotation of transition metal dichalcogenides
Authors:
Thibault Sohier,
Pedro M. M. C. de Melo,
Zeila Zanolli,
Matthieu Jean Verstraete
Abstract:
The transport and optical properties of semiconducting transition metal dichalcogenides around room temperature are dictated by electron-phonon scattering mechanisms within a complex, spin-textured and multi-valley electronic landscape. The relative positions of the valleys are critical, yet they are sensitive to external parameters and very difficult to determine directly. We propose a first-prin…
▽ More
The transport and optical properties of semiconducting transition metal dichalcogenides around room temperature are dictated by electron-phonon scattering mechanisms within a complex, spin-textured and multi-valley electronic landscape. The relative positions of the valleys are critical, yet they are sensitive to external parameters and very difficult to determine directly. We propose a first-principle model as a function valley positions to calculate carrier mobility and Kerr rotation angles. The model brings valuable insights, as well as quantitative predictions of macroscopic properties for a wide range of carrier density. The doping-dependant mobility displays a characteristic peak, the height depending on the position of the valleys. The Kerr rotation signal is enhanced when same spin-valleys are aligned, and quenched when opposite spin-valleys are populated. We provide guidelines to optimize these quantities with respect to experimental parameters, as well as the theoretical support for \emph{in situ} characterization of the valley positions.
△ Less
Submitted 24 October, 2022; v1 submitted 1 July, 2022;
originally announced July 2022.
-
Gate control of spin-layer-locking FETs and application to monolayer LuIO
Authors:
Rong Zhang,
Antimo Marrazzo,
Matthieu Verstraete,
Nicola Marzari,
Thibault Sohier
Abstract:
A recent 2D spinFET concept proposes to switch electrostatically between two separate sublayers with strong and opposite intrinsic Rashba effects. This concept exploits the spin-layer locking mechanism present in centrosymmetric materials with local dipole fields, where a weak electric field can easily manipulate just one of the spin channels. Here, we propose a novel monolayer material within thi…
▽ More
A recent 2D spinFET concept proposes to switch electrostatically between two separate sublayers with strong and opposite intrinsic Rashba effects. This concept exploits the spin-layer locking mechanism present in centrosymmetric materials with local dipole fields, where a weak electric field can easily manipulate just one of the spin channels. Here, we propose a novel monolayer material within this family, lutetium oxide iodide (LuIO). It displays one of the largest Rashba effects among 2D materials (up to $k_R = 0.08$ Å$^{-1}$), leading to a $π/2$ rotation of the spins over just 1 nm. The monolayer had been predicted to be exfoliable from its experimentally-known 3D bulk counterpart, with a binding energy even lower than graphene. We characterize and model with first-principles simulations the interplay of the two gate-controlled parameters for such devices: doping and spin channel selection. We show that the ability to split the spin channels in energy diminishes with doping, leading to specific gate-operation guidelines that can apply to all devices based on spin-layer locking.
△ Less
Submitted 15 June, 2021;
originally announced June 2021.
-
Hot-Carrier Cooling in High-Quality Graphene is Intrinsically Limited by Optical Phonons
Authors:
Eva A. A. Pogna,
Xiaoyu Jia,
Alessandro Principi,
Alexander Block,
Luca Banszerus,
Jincan Zhang,
Xiaoting Liu,
Thibault Sohier,
Stiven Forti,
Karuppasamy Soundarapandian,
Bernat Terrés,
Jake D. Mehew,
Chiara Trovatello,
Camilla Coletti,
Frank H. L. Koppens,
Mischa Bonn,
Niek van Hulst,
Matthieu J. Verstraete,
Hailin Peng,
Zhongfan Liu,
Christoph Stampfer,
Giulio Cerullo,
Klaas-Jan Tielrooij
Abstract:
Many promising optoelectronic devices, such as broadband photodetectors, nonlinear frequency converters, and building blocks for data communication systems, exploit photoexcited charge carriers in graphene. For these systems, it is essential to understand, and eventually control, the cooling dynamics of the photoinduced hot-carrier distribution. There is, however, still an active debate on the dif…
▽ More
Many promising optoelectronic devices, such as broadband photodetectors, nonlinear frequency converters, and building blocks for data communication systems, exploit photoexcited charge carriers in graphene. For these systems, it is essential to understand, and eventually control, the cooling dynamics of the photoinduced hot-carrier distribution. There is, however, still an active debate on the different mechanisms that contribute to hot-carrier cooling. In particular, the intrinsic cooling mechanism that ultimately limits the cooling dynamics remains an open question. Here, we address this question by studying two technologically relevant systems, consisting of high-quality graphene with a mobility >10,000 cm$^2$V$^{-1}$s$^{-1}$ and environments that do not efficiently take up electronic heat from graphene: WSe$_2$-encapsulated graphene and suspended graphene. We study the cooling dynamics of these two high-quality graphene systems using ultrafast pump-probe spectroscopy at room temperature. Cooling via disorder-assisted acoustic phonon scattering and out-of-plane heat transfer to the environment is relatively inefficient in these systems, predicting a cooling time of tens of picoseconds. However, we observe much faster cooling, on a timescale of a few picoseconds. We attribute this to an intrinsic cooling mechanism, where carriers in the hot-carrier distribution with enough kinetic energy emit optical phonons. During phonon emission, the electronic system continuously re-thermalizes, re-creating carriers with enough energy to emit optical phonons. We develop an analytical model that explains the observed dynamics, where cooling is eventually limited by optical-to-acoustic phonon coupling. These fundamental insights into the intrinsic cooling mechanism of hot carriers in graphene will play a key role in guiding the development of graphene-based optoelectronic devices.
△ Less
Submitted 5 March, 2021;
originally announced March 2021.
-
Remote free-carrier screening to boost the mobility of Fröhlich-limited 2D semiconductors
Authors:
Thibault Sohier,
Marco Gibertini,
Matthieu Verstraete
Abstract:
Van der Waals heterostructures provide a versatile tool to not only protect or control, but also enhance the properties of a 2D material. We use ab initio calculations and semi-analytical models to find strategies which boost the mobility of a current-carrying 2D semiconductor within an heterostructure. Free-carrier screening from a metallic "screener" layer remotely suppresses electron-phonon int…
▽ More
Van der Waals heterostructures provide a versatile tool to not only protect or control, but also enhance the properties of a 2D material. We use ab initio calculations and semi-analytical models to find strategies which boost the mobility of a current-carrying 2D semiconductor within an heterostructure. Free-carrier screening from a metallic "screener" layer remotely suppresses electron-phonon interactions in the current-carrying layer. This concept is most effective in 2D semiconductors whose scattering is dominated by screenable electron-phonon interactions, and in particular the Fröhlich coupling to polar-optical phonons. Such materials are common and characterised by overall low mobilities in the small doping limit, and much higher ones when the 2D material is doped enough for electron-phonon interactions to be screened by its own free carriers. We use GaSe as a prototype and place it in a heterostructure with doped graphene as the "screener" layer and BN as a separator. We develop an approach to determine the electrostatic response of any heterostructure by combining the responses of the individual layers computed within density-functional perturbation theory. Remote screening from graphene can suppress the long-wavelength Fröhlich interaction, leading to a consistently high mobility around $500$ to $600$ cm$^2$/Vs for carrier densities in GaSe from $10^{11}$ to $10^{13}$ cm$^{-2}$. Notably, the low-doping mobility is enhanced by a factor 2.5. This remote free-carrier screening is more efficient than more conventional manipulation of the dielectric environment, and it is most effective when the separator (BN) is thin.
△ Less
Submitted 15 February, 2021; v1 submitted 10 November, 2020;
originally announced November 2020.
-
Profiling novel high-conductivity 2D semiconductors
Authors:
Thibault Sohier,
Marco Gibertini,
Nicola Marzari
Abstract:
When complex mechanisms are involved, pinpointing high-performance materials within large databases is a major challenge in materials discovery. We focus here on phonon-limited conductivities, and study 2D semiconductors doped by field effects. Using state-of-the-art density-functional perturbation theory and Boltzmann transport equation, we discuss 11 monolayers with outstanding transport propert…
▽ More
When complex mechanisms are involved, pinpointing high-performance materials within large databases is a major challenge in materials discovery. We focus here on phonon-limited conductivities, and study 2D semiconductors doped by field effects. Using state-of-the-art density-functional perturbation theory and Boltzmann transport equation, we discuss 11 monolayers with outstanding transport properties. These materials are selected from a computational database of exfoliable materials providing monolayers that are dynamically stable and that do not have more than 6 atoms per unit cell. We first analyse electron-phonon scattering in two well-known systems: electron-doped InSe and hole-doped phosphorene. Both are single-valley systems with weak electron-phonon interactions, but they represent two distinct pathways to fast transport: a steep and deep isotropic valley for the former and strongly anisotropic electron-phonon physics for the latter. We identify similar features in the database and compute the conductivities of the relevant monolayers. This process yields several high-conductivity materials, some of them only very recently emerging in the literature (GaSe, Bi$_2$SeTe$_2$, Bi$_2$Se$_3$, Sb$_2$SeTe$_2$), others never discussed in this context (AlLiTe$_2$, BiClTe, ClGaTe, AuI). Comparing these 11 monolayers in detail, we discuss how the strength and angular dependency of the electron-phonon scattering drives key differences in the transport performance of materials despite similar valley structure. We also discuss the high conductivity of hole-doped WSe$_2$, and how this case study shows the limitations of a selection process that would be based on band properties alone.
△ Less
Submitted 31 July, 2020;
originally announced July 2020.
-
Extraordinary high room-temperature carrier mobility in graphene-WSe$_2$ heterostructures
Authors:
L. Banszerus,
T. Sohier,
A. Epping,
F. Winkler,
F. Libisch,
F. Haupt,
K. Watanabe,
T. Taniguchi,
K. Müller-Caspary,
N. Marzari,
F. Mauri,
B. Beschoten,
C. Stampfer
Abstract:
High carrier mobilities play a fundamental role for high-frequency electronics, integrated optoelectronics as well as for sensor and spintronic applications, where device performance is directly linked to the magnitude of the carrier mobility. Van der Waals heterostructures formed by graphene and hexagonal boron nitride (hBN) already outperform all known materials in terms of room temperature mobi…
▽ More
High carrier mobilities play a fundamental role for high-frequency electronics, integrated optoelectronics as well as for sensor and spintronic applications, where device performance is directly linked to the magnitude of the carrier mobility. Van der Waals heterostructures formed by graphene and hexagonal boron nitride (hBN) already outperform all known materials in terms of room temperature mobility. Here, we show that the mobility of today's best graphene/hBN devices can be surpassed by more than a factor of three by heterostructures formed by tungsten diselenide (WSe$_2$), graphene and hBN, which can have mobilities as high as 350,000 cm$^2$/(Vs) at room temperature, and resistivities as low as 15 Ohm. The resistivity of these devices shows a much weaker temperature dependence than the one of graphene on any other known substrate. The origin of this behaviour points to modified acoustic phonon bands in graphene and questions our understanding of electron-phonon scattering in van der Waals heterostructures.
△ Less
Submitted 20 September, 2019;
originally announced September 2019.
-
Electric field exfoliation and high-T$_{\text{C}}$ superconductivity in field-effect hole-doped hydrogenated diamond (111)
Authors:
Davide Romanin,
Thibault Sohier,
Dario Daghero,
Francesco Mauri,
Renato S. Gonnelli,
Matteo Calandra
Abstract:
We investigate the possible occurrence of field-effect induced superconductivity in the hydrogenated $(111)$ diamond surface by first-principles calculations. By computing the band alignment between bulk diamond and the hydrogenated surface we show that the electric field exfoliates the sample, separating the electronic states at the valence band top from the bulk projected ones. At the hole dopin…
▽ More
We investigate the possible occurrence of field-effect induced superconductivity in the hydrogenated $(111)$ diamond surface by first-principles calculations. By computing the band alignment between bulk diamond and the hydrogenated surface we show that the electric field exfoliates the sample, separating the electronic states at the valence band top from the bulk projected ones. At the hole doping values considered here, ranging from $n=2.84\times10^{13}$cm$^{-2}$ to $n=6\times 10^{14} $ cm$^{-2}$, the valence band top is composed of up to three electronic bands hosting holes with different effective masses. These bands resemble those of the undoped surface, but they are heavily modified by the electric field and differ substantially from a rigid doping picture. We calculate superconducting properties by including the effects of charging of the slab and of the electric field on the structural properties, electronic structure, phonon dispersion and electron-phonon coupling. We find that at doping as large as $n=6\times 10^{14} $ cm$^{-2}$, the electron-phonon interaction is $λ=0.81$ and superconductivity emerges with $T_{\text{C}}\approx 29-36$K. Superconductivity is mostly supported by in-plane diamond phonon vibrations and to a lesser extent by some out-of-plane vibrations. The relevant electron-phonon scattering processes involve both intra and interband scattering so that superconductivity is multiband in nature.
△ Less
Submitted 5 August, 2019;
originally announced August 2019.
-
Valley-engineering mobilities in two-dimensional materials
Authors:
Thibault Sohier,
Marco Gibertini,
Davide Campi,
Giovanni Pizzi,
Nicola Marzari
Abstract:
Two-dimensional materials are emerging as a promising platform for ultrathin channels in field-effect transistors. To this aim, novel high-mobility semiconductors need to be found or engineered. While extrinsic mechanisms can in general be minimized by improving fabrication processes, the suppression of intrinsic scattering (driven e.g. by electron-phonon interactions) requires to modify the elect…
▽ More
Two-dimensional materials are emerging as a promising platform for ultrathin channels in field-effect transistors. To this aim, novel high-mobility semiconductors need to be found or engineered. While extrinsic mechanisms can in general be minimized by improving fabrication processes, the suppression of intrinsic scattering (driven e.g. by electron-phonon interactions) requires to modify the electronic or vibrational properties of the material. Since intervalley scattering critically affects mobilities, a powerful approach to enhance transport performance relies on engineering the valley structure. We show here the power of this strategy using uniaxial strain to lift degeneracies and suppress scattering into entire valleys, dramatically improving performance. This is shown in detail for arsenene, where a 2% strain stops scattering into 4 of the 6 valleys, and leads to a 600% increase in mobility. The mechanism is general and can be applied to many other materials, including in particular the isostructural antimonene and blue phosphorene.
△ Less
Submitted 29 May, 2019; v1 submitted 28 February, 2019;
originally announced February 2019.
-
Enhanced electron-phonon interaction in multi-valley materials
Authors:
Evgeniy Ponomarev,
Thibault Sohier,
Marco Gibertini,
Helmuth Berger,
Nicola Marzari,
Nicolas Ubrig,
Alberto F. Morpurgo
Abstract:
Through a combined theoretical and experimental effort, we uncover a yet unidentified mechanism that strengthens considerably electron-phonon coupling in materials where electron accumulation leads to population of multiple valleys. Taking atomically-thin transition-metal dichalcogenides as prototypical examples, we establish that the mechanism results from a phonon-induced out-of-phase energy shi…
▽ More
Through a combined theoretical and experimental effort, we uncover a yet unidentified mechanism that strengthens considerably electron-phonon coupling in materials where electron accumulation leads to population of multiple valleys. Taking atomically-thin transition-metal dichalcogenides as prototypical examples, we establish that the mechanism results from a phonon-induced out-of-phase energy shift of the different valleys, which causes inter-valley charge transfer and reduces the effectiveness of electrostatic screening, thus enhancing electron-phonon interactions. The effect is physically robust, it can play a role in many materials and phenomena, as we illustrate by discussing experimental evidence for its relevance in the occurrence of superconductivity.
(short abstract due to size limitations - full abstract in the manuscript)
△ Less
Submitted 23 January, 2019;
originally announced January 2019.
-
Mobility of 2D materials from first principles in an accurate and automated framework
Authors:
Thibault Sohier,
Davide Campi,
Nicola Marzari,
Marco Gibertini
Abstract:
We present a first-principles approach to compute the transport properties of 2D materials in an accurate and automated framework. We use density-functional perturbation theory in the appropriate bidimensional setup with open-boundary conditions in the third direction. The materials are charged by field effect via planar counter-charges. In this approach, we obtain electron-phonon matrix elements…
▽ More
We present a first-principles approach to compute the transport properties of 2D materials in an accurate and automated framework. We use density-functional perturbation theory in the appropriate bidimensional setup with open-boundary conditions in the third direction. The materials are charged by field effect via planar counter-charges. In this approach, we obtain electron-phonon matrix elements in which dimensionality and doping effects are inherently accounted for, without the need for post-processing corrections. This treatment highlights some unexpected consequences, such as an increase of electron-phonon coupling with doping in transition-metal dichalcogenides.We use symmetries extensively and identify pockets of relevant electronic states to minimize the number of electron-phonon interactions to compute; the integrodifferential Boltzmann transport equation is then linearized and solved beyond the relaxation-time approximation. We apply the entire protocol to a set of much studied materials with diverse electronic and vibrational band structures: electron-doped MoS2, WS2, WSe2, phosphorene, arsenene, and hole-doped phosphorene. Among these, hole-doped phosphorene is found to have the highest mobility, with a room temperature value around 600 cm^2/(Vs). Last, we identify the factors that affect most phonon-limited mobilities, such as the number and the anisotropy of electron and hole pockets, to provide a broader understanding of the driving forces behind high mobilities in two-dimensional materials.
△ Less
Submitted 22 November, 2018; v1 submitted 31 August, 2018;
originally announced August 2018.
-
Density functional perturbation theory for gated two-dimensional heterostructures: Theoretical developments and application to flexural phonons in graphene
Authors:
Thibault Sohier,
Matteo Calandra,
Francesco Mauri
Abstract:
The ability to perform first-principles calculations of electronic and vibrational properties of two-dimensional heterostructures in a field-effect setup is crucial for the understanding and design of next-generation devices. We present here an implementation of density functional perturbation theories tailored for the case of two-dimensional heterostructures in field-effect configuration. Key ing…
▽ More
The ability to perform first-principles calculations of electronic and vibrational properties of two-dimensional heterostructures in a field-effect setup is crucial for the understanding and design of next-generation devices. We present here an implementation of density functional perturbation theories tailored for the case of two-dimensional heterostructures in field-effect configuration. Key ingredients are the inclusion of a truncated Coulomb interaction in the direction perpendicular to the slab and the possibility of simulating charging of the slab via field-effects. With this implementation we can access total energies, force and stress tensors, the vibrational properties and the electron-phonon interaction. We demonstrate the relevance of the method by studying flexural acoustic phonons and their coupling to electrons in graphene doped by field-effect. In particular, we show that while the electron-phonon coupling to those phonons can be significant in neutral graphene, it is strongly screened and negligible in doped graphene, in disagreement with other recent first-principles reports. Consequently, the gate-induced coupling with flexural acoustic modes would not be detectable in transport measurements on doped graphene.
△ Less
Submitted 4 September, 2017; v1 submitted 14 May, 2017;
originally announced May 2017.
-
Breakdown of optical phonons' splitting in two-dimensional materials
Authors:
Thibault Sohier,
Marco Gibertini,
Matteo Calandra,
Francesco Mauri,
Nicola Marzari
Abstract:
We investigate the long-wavelength dispersion of longitudinal and transverse optical phonon modes in polar two-dimensional materials, multilayers, and their heterostructures. Using analytical models and density-functional perturbation theory in a two-dimensional framework, we show that, at variance with the three-dimensional case, these modes are degenerate at the zone center but the macroscopic e…
▽ More
We investigate the long-wavelength dispersion of longitudinal and transverse optical phonon modes in polar two-dimensional materials, multilayers, and their heterostructures. Using analytical models and density-functional perturbation theory in a two-dimensional framework, we show that, at variance with the three-dimensional case, these modes are degenerate at the zone center but the macroscopic electric field associated with the longitudinal-optical modes gives rise to a finite slope at the zone center in their corresponding phonon dispersions. This slope increases linearly with the number of layers and it is determined solely by the Born effective charges of the material and the dielectric properties of the surrounding media. Screening from the environment can greatly reduce the slope splitting between the longitudinal and transverse optical modes and can be seen in the experimentally relevant case of boron nitride-graphene heterostructures. As the phonon momentum increases, the intrinsic screening properties of the two-dimensional material dictate the transition to a momentum-independent splitting similar to that of three-dimensional materials. These considerations are essential to understand electrical transport and optical coupling in two-dimensional systems.
△ Less
Submitted 10 May, 2017; v1 submitted 21 December, 2016;
originally announced December 2016.
-
Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds
Authors:
Nicolas Mounet,
Marco Gibertini,
Philippe Schwaller,
Davide Campi,
Andrius Merkys,
Antimo Marrazzo,
Thibault Sohier,
Ivano E. Castelli,
Andrea Cepellotti,
Giovanni Pizzi,
Nicola Marzari
Abstract:
We search for novel two-dimensional materials that can be easily exfoliated from their parent compounds. Starting from 108423 unique, experimentally known three-dimensional compounds we identify a subset of 5619 that appear layered according to robust geometric and bonding criteria. High-throughput calculations using van-der-Waals density-functional theory, validated against experimental structura…
▽ More
We search for novel two-dimensional materials that can be easily exfoliated from their parent compounds. Starting from 108423 unique, experimentally known three-dimensional compounds we identify a subset of 5619 that appear layered according to robust geometric and bonding criteria. High-throughput calculations using van-der-Waals density-functional theory, validated against experimental structural data and calculated random-phase-approximation binding energies, allow to identify 1825 compounds that are either easily or potentially exfoliable, including all that are commonly exfoliated experimentally. In particular, the subset of 1036 easily exfoliable cases---layered materials held together mostly by dispersion interactions and with binding energies up to $30-35$ meV$\cdot\textÅ^{-2}$---provides a wealth of novel structural prototypes and simple ternary compounds, and a large portfolio to search materials for optimal properties. For the 258 compounds with up to 6 atoms per primitive cell we comprehensively explore vibrational, electronic, magnetic, and topological properties, identifying in particular 56 ferromagnetic and antiferromagnetic systems, including half-metals and half-semiconductors.
△ Less
Submitted 22 July, 2020; v1 submitted 16 November, 2016;
originally announced November 2016.
-
Two-dimensional Fröhlich interaction in transition-metal dichalcogenide monolayers: Theoretical modeling and first-principles calculations
Authors:
Thibault Sohier,
Matteo Calandra,
Francesco Mauri
Abstract:
We perform ab initio calculations of the coupling between electrons and small-momentum polar-optical phonons in monolayer transition metal dichalcogenides of the 2H type: MoS2, MoSe2, MoTe2, WS2, and WSe2. The so-called Fröhlich interaction is fundamentally affected by the dimensionality of the system. In a plane-wave framework with periodic boundary conditions, this coupling is affected by the sp…
▽ More
We perform ab initio calculations of the coupling between electrons and small-momentum polar-optical phonons in monolayer transition metal dichalcogenides of the 2H type: MoS2, MoSe2, MoTe2, WS2, and WSe2. The so-called Fröhlich interaction is fundamentally affected by the dimensionality of the system. In a plane-wave framework with periodic boundary conditions, this coupling is affected by the spurious interaction between the 2D material and its periodic images. To overcome this, we perform density functional perturbation theory calculations with a truncated Coulomb interaction in the out-of-plane direction. We show that the 2D Fröhlich interaction is much stronger than assumed in previous ab initio studies. We provide analytical models depending on the effective charges and dielectric properties of the materials to interpret our ab initio calculations. Screening is shown to play a fundamental role in the phonon-momentum dependency of the polar-optical coupling, with a crossover between two regimes depending on the dielectric properties of the material relative to its environment. The Fröhlich interaction is screened by the dielectric environment in the limit of small phonon momenta and sharply decreases due to stronger screening by the monolayer at finite momenta. The small-momentum regime of the ab initio Fröhlich interaction is reproduced by a simple analytical model, for which we provide the necessary parameters. At larger momenta, however, direct ab initio calculations of electron-phonon interactions are necessary to capture band-specific effects. We compute and compare the carrier relaxation times associated to the scattering by both LO and A1 phonon modes. While both modes are capable of relaxing carriers on timescales under the picosecond at room temperature, their absolute and relative importance vary strongly depending on the material, the band, and the substrate.
△ Less
Submitted 4 November, 2016; v1 submitted 26 May, 2016;
originally announced May 2016.
-
Density-functional calculation of static screening in 2D materials: the long-wavelength dielectric function of graphene
Authors:
Thibault Sohier,
Matteo Calandra,
Francesco Mauri
Abstract:
We calculate the long-wavelength static screening properties of both neutral and doped graphene in the framework of density-functional theory. We use a plane-wave approach with periodic images in the third dimension and truncate the Coulomb interactions to eliminate spurious interlayer screening. We carefully address the issue of extracting two dimensional dielectric properties from simulated thre…
▽ More
We calculate the long-wavelength static screening properties of both neutral and doped graphene in the framework of density-functional theory. We use a plane-wave approach with periodic images in the third dimension and truncate the Coulomb interactions to eliminate spurious interlayer screening. We carefully address the issue of extracting two dimensional dielectric properties from simulated three-dimensional potentials. We compare this method with analytical expressions derived for two dimensional massless Dirac fermions in the random phase approximation. We evaluate the contributions of the deviation from conical bands, exchange-correlation and local-fields. For momenta smaller than twice the Fermi wavevector, the static screening of graphene within the density-functional perturbative approach agrees with the results for conical bands within random phase approximation and neglecting local fields. For larger momenta, we find that the analytical model underestimates the static dielectric function by $\approx 10%$, mainly due to the conical band approximation.
△ Less
Submitted 9 April, 2015; v1 submitted 9 March, 2015;
originally announced March 2015.
-
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.
-
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.