Showing 1–50 of 59 results for author: Rignanese, G

Optical materials discovery and design with federated databases and machine learning

Authors: Victor Trinquet, Matthew L. Evans, Cameron J. Hargreaves, Pierre-Paul De Breuck, Gian-Marco Rignanese

Abstract: Combinatorial and guided screening of materials space with density-functional theory and related approaches has provided a wealth of hypothetical inorganic materials, which are increasingly tabulated in open databases. The OPTIMADE API is a standardised format for representing crystal structures, their measured and computed properties, and the methods for querying and filtering them from remote re… ▽ More Combinatorial and guided screening of materials space with density-functional theory and related approaches has provided a wealth of hypothetical inorganic materials, which are increasingly tabulated in open databases. The OPTIMADE API is a standardised format for representing crystal structures, their measured and computed properties, and the methods for querying and filtering them from remote resources. Currently, the OPTIMADE federation spans over 20 data providers, rendering over 30 million structures accessible in this way, many of which are novel and have only recently been suggested by machine learning-based approaches. In this work, we outline our approach to non-exhaustively screen this dynamic trove of structures for the next-generation of optical materials. By applying MODNet, a neural network-based model for property prediction, within a combined active learning and high-throughput computation framework, we isolate particular structures and chemistries that should be most fruitful for further theoretical calculations and for experimental study as high-refractive-index materials. By making explicit use of automated calculations, federated dataset curation and machine learning, and by releasing these publicly, the workflows presented here can be periodically re-assessed as new databases implement OPTIMADE, and new hypothetical materials are suggested. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Assessing carrier mobility, dopability, and defect tolerance in the chalcogenide perovskite BaZrS$_3$

Authors: Zhenkun Yuan, Diana Dahliah, Romain Claes, Andrew Pike, David P. Fenning, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: The chalcogenide perovskite BaZrS$_3$ has attracted much attention as a promising solar absorber for thin-film photovoltaics. Here, we use first-principles calculations to evaluate its carrier transport and defect properties. We find that BaZrS$_3$ has a phonon-limited electron mobility of 37 cm$^2$/Vs comparable to that in halide perovskites but lower hole mobility of 11 cm$^2$/Vs. The defect com… ▽ More The chalcogenide perovskite BaZrS$_3$ has attracted much attention as a promising solar absorber for thin-film photovoltaics. Here, we use first-principles calculations to evaluate its carrier transport and defect properties. We find that BaZrS$_3$ has a phonon-limited electron mobility of 37 cm$^2$/Vs comparable to that in halide perovskites but lower hole mobility of 11 cm$^2$/Vs. The defect computations indicate that BaZrS$_3$ is intrinsically n-type due to shallow sulfur vacancies, but that strong compensation by sulfur vacancies will prevent attempts to make it p-type. We also establish that BaZrS$_3$ is a defect-tolerant absorber with few low formation energy, deep intrinsic defects. Among the deep defects, sulfur interstitials are the strongest nonradiative recombination centers which in sulfur-rich conditions would limit the carrier lifetime to 10 ns. Our work highlights the material's intrinsic limitations in carrier mobility and suggests suppressing the formation of sulfur interstitials to reach long carrier lifetime. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Facilities and practices for linear response Hubbard parameters U and J in Abinit

Authors: Lórien MacEnulty, Matteo Giantomassi, Bernard Amadon, Gian-Marco Rignanese, David D. O'Regan

Abstract: Members of the DFT+U family of functionals are increasingly prevalent methods of addressing errors intrinsic to (semi-) local exchange-correlation functionals at minimum computational cost, but require their parameters U and J to be calculated in situ for a given system of interest, simulation scheme, and runtime parameters. The SCF linear response approach offers ab initio acquisition of the U an… ▽ More Members of the DFT+U family of functionals are increasingly prevalent methods of addressing errors intrinsic to (semi-) local exchange-correlation functionals at minimum computational cost, but require their parameters U and J to be calculated in situ for a given system of interest, simulation scheme, and runtime parameters. The SCF linear response approach offers ab initio acquisition of the U and has recently been extended to compute the J analogously, which measures localized errors related to exchange-like effects. We introduce a renovated post-processor, the lrUJ utility, together with this detailed best-practices guide, to enable users of the popular, open-source Abinit first-principles simulation suite to engage easily with in situ Hubbard parameters and streamline their incorporation into material simulations of interest. Features of this utility, which may also interest users and developers of other DFT codes, include $n$-degree polynomial regression, error analysis, Python plotting facilities, didactic documentation, and avenues for further developments. In this technical introduction and guide, we place particular emphasis on the intricacies and potential pitfalls introduced by the projector augmented wave (PAW) method, SCF mixing schemes, and non-linear response, several of which are translatable to DFT+U(+J) implementations in other packages. △ Less

Submitted 10 July, 2024; v1 submitted 9 April, 2024; originally announced April 2024.

Comments: 46 pages, 6 figures (+3 figures in appendix), 2 tables (+1 table in appendix)

Systematic assessment of various universal machine-learning interatomic potentials

Authors: Haochen Yu, Matteo Giantomassi, Giuliana Materzanini, Junjie Wang, Gian-Marco Rignanese

Abstract: Machine-learning interatomic potentials have revolutionized materials modeling at the atomic scale. Thanks to these, it is now indeed possible to perform simulations of \abinitio quality over very large time and length scales. More recently, various universal machine-learning models have been proposed as an out-of-box approach avoiding the need to train and validate specific potentials for each pa… ▽ More Machine-learning interatomic potentials have revolutionized materials modeling at the atomic scale. Thanks to these, it is now indeed possible to perform simulations of \abinitio quality over very large time and length scales. More recently, various universal machine-learning models have been proposed as an out-of-box approach avoiding the need to train and validate specific potentials for each particular material of interest. In this paper, we review and evaluate five different universal machine-learning interatomic potentials (uMLIPs), all based on graph neural network architectures which have demonstrated transferability from one chemical system to another. The evaluation procedure relies on data both from a recent verification study of density-functional-theory implementations and from the Materials Project. Through this comprehensive evaluation, we aim to provide guidance to materials scientists in selecting suitable models for their specific research problems, offer recommendations for model selection and optimization, and stimulate discussion on potential areas for improvement in current machine-learning methodologies in materials science. △ Less

Submitted 21 March, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 7 figures

Developments and applications of the OPTIMADE API for materials discovery, design, and data exchange

Authors: Matthew L. Evans, Johan Bergsma, Andrius Merkys, Casper W. Andersen, Oskar B. Andersson, Daniel Beltrán, Evgeny Blokhin, Tara M. Boland, Rubén Castañeda Balderas, Kamal Choudhary, Alberto Díaz Díaz, Rodrigo Domínguez García, Hagen Eckert, Kristjan Eimre, María Elena Fuentes Montero, Adam M. Krajewski, Jens Jørgen Mortensen, José Manuel Nápoles Duarte, Jacob Pietryga, Ji Qi, Felipe de Jesús Trejo Carrillo, Antanas Vaitkus, Jusong Yu, Adam Zettel, Pedro Baptista de Castro , et al. (34 additional authors not shown)

Abstract: The Open Databases Integration for Materials Design (OPTIMADE) application programming interface (API) empowers users with holistic access to a growing federation of databases, enhancing the accessibility and discoverability of materials and chemical data. Since the first release of the OPTIMADE specification (v1.0), the API has undergone significant development, leading to the upcoming v1.2 relea… ▽ More The Open Databases Integration for Materials Design (OPTIMADE) application programming interface (API) empowers users with holistic access to a growing federation of databases, enhancing the accessibility and discoverability of materials and chemical data. Since the first release of the OPTIMADE specification (v1.0), the API has undergone significant development, leading to the upcoming v1.2 release, and has underpinned multiple scientific studies. In this work, we highlight the latest features of the API format, accompanying software tools, and provide an update on the implementation of OPTIMADE in contributing materials databases. We end by providing several use cases that demonstrate the utility of the OPTIMADE API in materials research that continue to drive its ongoing development. △ Less

Submitted 5 April, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

Native point defects in HgCdTe infrared detector material: Identifying deep centers from first principles

Authors: Wei Chen, Gian-Marco Rignanese, Jifeng Liu, Geoffroy Hautier

Abstract: We investigate the native point defects in the long-wavelength infrared (LWIR) detector material Hg$_{0.75}$Cd$_{0.25}$Te using a dielectric-dependent hybrid density functional combined with spin-orbit coupling. Characterizing these point defects is essential as they are responsible for intrinsic doping and nonradiative recombination centers in the detector material. The dielectric-dependent hybri… ▽ More We investigate the native point defects in the long-wavelength infrared (LWIR) detector material Hg$_{0.75}$Cd$_{0.25}$Te using a dielectric-dependent hybrid density functional combined with spin-orbit coupling. Characterizing these point defects is essential as they are responsible for intrinsic doping and nonradiative recombination centers in the detector material. The dielectric-dependent hybrid functional allows for an accurate description of the band gap ($E_g$) for Hg$_{1-x}$Cd$_{x}$Te (MCT) over the entire compositional range, a level of accuracy challenging with standard hybrid functionals. Our comprehensive examination of the native point defects confirms that cation vacancies $V_\text{Hg(Cd)}$ are the primary sources of $p$-type conductivity in the LWIR material given their low defect formation energies and the presence of a shallow acceptor level ($-$/0) near the valence-band maximum (VBM). In addition to the shallow acceptor level, the cation vacancies exhibit a deep charge transition level (2$-$/$-$) situated near the midgap, characteristic of nonradiative recombination centers. Our results indicate that Hg interstitial could also be a deep center in the LWIR MCT through a metastable configuration under the Hg-rich growth conditions. While an isolated Te antisite does not show deep levels, the formation of $V_\text{Hg}$-Te$_\text{Hg}$ defect complex introduces a deep acceptor level within the band gap. △ Less

Submitted 13 March, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 12 pages, 7 figures

Surface Enhanced Infrared Absorption mechanism and modification of the plasmonic response

Authors: Tanguy Colleu, Adam Fekete, Xavier Gonze, Alexandre Cloots, Vincent Liégeois, Gian-Marco Rignanese, Luc Henrard

Abstract: Surface Enhanced Infrared Absorption (SEIRA) is an experimental method where trace amount of a compound can be detected with high sensibility. This high detection sensibility is the result of the interaction of the molecules with a localized plasmon, usually from a metallic nano-particle. In this study we numerically investigate by discrete dipole approximation the origin of the Fano-like response… ▽ More Surface Enhanced Infrared Absorption (SEIRA) is an experimental method where trace amount of a compound can be detected with high sensibility. This high detection sensibility is the result of the interaction of the molecules with a localized plasmon, usually from a metallic nano-particle. In this study we numerically investigate by discrete dipole approximation the origin of the Fano-like response of the system, including the induced transparency when the plasmon resonance and the molecular vibrational mode coincide. The detailed analysis of the localization of the absorption show that the modification of the absorption cross-section when the molecule is present comes from a change of the plasmonic resonance, not from the direct molecular response which is negligible. This sheds a new light on the SEIRA mechanism. In particular, it demonstrates that the sensibility is associated with the influence of the molecule on the plasmon resonance rather than with the local field enhancement itself. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 12 pages, 5 figures

Discovery of the Zintl-phosphide BaCd$_{2}$P$_{2}$ as a long carrier lifetime and stable solar absorber

Authors: Zhenkun Yuan, Diana Dahliah, Muhammad Rubaiat Hasan, Gideon Kassa, Andrew Pike, Shaham Quadir, Romain Claes, Cierra Chandler, Yihuang Xiong, Victoria Kyveryga, Philip Yox, Gian-Marco Rignanese, Ismaila Dabo, Andriy Zakutayev, David P. Fenning, Obadiah G. Reid, Sage Bauers, Jifeng Liu, Kirill Kovnir, Geoffroy Hautier

Abstract: Thin-film photovoltaics offers a path to significantly decarbonize our energy production. Unfortunately, current materials commercialized or under development as thin-film solar cell absorbers are far from optimal as they show either low power conversion efficiency or issues with earth-abundance and stability. Entirely new and disruptive materials platforms are rarely discovered as the search for… ▽ More Thin-film photovoltaics offers a path to significantly decarbonize our energy production. Unfortunately, current materials commercialized or under development as thin-film solar cell absorbers are far from optimal as they show either low power conversion efficiency or issues with earth-abundance and stability. Entirely new and disruptive materials platforms are rarely discovered as the search for new solar absorbers is traditionally slow and serendipitous. Here, we use first principles high-throughput screening to accelerate this process. We identify new solar absorbers among known inorganic compounds using considerations on band gap, carrier transport, optical absorption but also on intrinsic defects which can strongly limit the carrier lifetime and ultimately the solar cell efficiency. Screening about 40,000 materials, we discover the Zintl-phosphide BaCd$_{2}$P$_{2}$ as a potential high-efficiency solar absorber. Follow-up experimental work confirms the predicted promises of BaCd$_{2}$P$_{2}$ highlighting an optimal band gap for visible absorption, bright photoluminescence, and long carrier lifetime of up to 30 ns even for unoptimized powder samples. Importantly, BaCd$_{2}$P$_{2}$ does not contain any critical elements and is highly stable in air and water. Our work opens an avenue for a new family of stable, earth-abundant, high-performance Zintl-based solar absorbers. It also demonstrates how recent advances in first principles computation can accelerate the search of photovoltaic materials by combining high-throughput screening with experiment. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Journal ref: Joule 8 (2024) 1412-1429

A substitutional quantum defect in WS$_2$ discovered by high-throughput computational screening and fabricated by site-selective STM manipulation

Authors: John C. Thomas, Wei Chen, Yihuang Xiong, Bradford A. Barker, Junze Zhou, Weiru Chen, Antonio Rossi, Nolan Kelly, Zhuohang Yu, Da Zhou, Shalini Kumari, Edward S. Barnard, Joshua A. Robinson, Mauricio Terrones, Adam Schwartzberg, D. Frank Ogletree, Eli Rotenberg, Marcus M. Noack, Sinéad Griffin, Archana Raja, David A. Strubbe, Gian-Marco Rignanese, Alexander Weber-Bargioni, Geoffroy Hautier

Abstract: Point defects in two-dimensional materials are of key interest for quantum information science. However, the space of possible defects is immense, making the identification of high-performance quantum defects extremely challenging. Here, we perform high-throughput (HT) first-principles computational screening to search for promising quantum defects within WS$_2$, which present localized levels in… ▽ More Point defects in two-dimensional materials are of key interest for quantum information science. However, the space of possible defects is immense, making the identification of high-performance quantum defects extremely challenging. Here, we perform high-throughput (HT) first-principles computational screening to search for promising quantum defects within WS$_2$, which present localized levels in the band gap that can lead to bright optical transitions in the visible or telecom regime. Our computed database spans more than 700 charged defects formed through substitution on the tungsten or sulfur site. We found that sulfur substitutions enable the most promising quantum defects. We computationally identify the neutral cobalt substitution to sulfur (Co$_{\rm S}^{0}$) as very promising and fabricate it with scanning tunneling microscopy (STM). The Co$_{\rm S}^{0}$ electronic structure measured by STM agrees with first principles and showcases an attractive new quantum defect. Our work shows how HT computational screening and novel defect synthesis routes can be combined to design new quantum defects. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 38 pages, 19 figures

How to verify the precision of density-functional-theory implementations via reproducible and universal workflows

Authors: Emanuele Bosoni, Louis Beal, Marnik Bercx, Peter Blaha, Stefan Blügel, Jens Bröder, Martin Callsen, Stefaan Cottenier, Augustin Degomme, Vladimir Dikan, Kristjan Eimre, Espen Flage-Larsen, Marco Fornari, Alberto Garcia, Luigi Genovese, Matteo Giantomassi, Sebastiaan P. Huber, Henning Janssen, Georg Kastlunger, Matthias Krack, Georg Kresse, Thomas D. Kühne, Kurt Lejaeghere, Georg K. H. Madsen, Martijn Marsman , et al. (20 additional authors not shown)

Abstract: In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a firs… ▽ More In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies). △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: Main text: 23 pages, 4 figures. Supplementary: 68 pages. Nature Review Physics 2023

Journal ref: Nat. Rev. Phys. 6, 45 (2024)

Nonunique fraction of Fock exchange for defects in two-dimensional materials

Authors: Wei Chen, Sinéad M. Griffin, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: By investigating the vacancy and substitutional defects in monolayer WS$_2$ with hybrid functionals, we find that there is no unique amount of Fock exchange that concurrently satisfies the generalized Koopmans' condition and reproduces the band gap and band-edge positions. Fixing the mixing parameter of Fock exchange based upon the band gap can lead to qualitatively incorrect defect physics in two… ▽ More By investigating the vacancy and substitutional defects in monolayer WS$_2$ with hybrid functionals, we find that there is no unique amount of Fock exchange that concurrently satisfies the generalized Koopmans' condition and reproduces the band gap and band-edge positions. Fixing the mixing parameter of Fock exchange based upon the band gap can lead to qualitatively incorrect defect physics in two-dimensional materials. Instead, excellent agreement is achieved with experiment and many-body perturbation theory within $GW$ approximation once the mixing parameters are tuned individually for the defects and the band edges. We show the departure from a unique optimized mixing parameter is inherent to two-dimensional systems as the band edges experience a reduced screening whilst the localized defects are subject to bulklike screening. △ Less

Submitted 9 October, 2022; v1 submitted 1 June, 2022; originally announced June 2022.

Comments: Main manuscript (7 pages, 3 figures) with supplemental material (3 pages, 5 figures)

Shared Metadata for Data-Centric Materials Science

Authors: Luca M. Ghiringhelli, Carsten Baldauf, Tristan Bereau, Sandor Brockhauser, Christian Carbogno, Javad Chamanara, Stefano Cozzini, Stefano Curtarolo, Claudia Draxl, Shyam Dwaraknath, Ádám Fekete, James Kermode, Christoph T. Koch, Markus Kühbach, Alvin Noe Ladines, Patrick Lambrix, Maja-Olivia Lenz-Himmer, Sergey Levchenko, Micael Oliveira, Adam Michalchuk, Ron Miller, Berk Onat, Pasquale Pavone, Giovanni Pizzi, Benjamin Regler , et al. (10 additional authors not shown)

Abstract: The expansive production of data in materials science, their widespread sharing and repurposing requires educated support and stewardship. In order to ensure that this need helps rather than hinders scientific work, the implementation of the FAIR-data principles (Findable, Accessible, Interoperable, and Reusable) must not be too narrow. Besides, the wider materials-science community ought to agree… ▽ More The expansive production of data in materials science, their widespread sharing and repurposing requires educated support and stewardship. In order to ensure that this need helps rather than hinders scientific work, the implementation of the FAIR-data principles (Findable, Accessible, Interoperable, and Reusable) must not be too narrow. Besides, the wider materials-science community ought to agree on the strategies to tackle the challenges that are specific to its data, both from computations and experiments. In this paper, we present the result of the discussions held at the workshop on "Shared Metadata and Data Formats for Big-Data Driven Materials Science". We start from an operative definition of metadata, and what features a FAIR-compliant metadata schema should have. We will mainly focus on computational materials-science data and propose a constructive approach for the FAIRification of the (meta)data related to ground-state and excited-states calculations, potential-energy sampling, and generalized workflows. Finally, challenges with the FAIRification of experimental (meta)data and materials-science ontologies are presented together with an outlook of how to meet them. △ Less

Submitted 23 August, 2023; v1 submitted 29 May, 2022; originally announced May 2022.

A simple denoising approach to exploit multi-fidelity data for machine learning materials properties

Authors: Xiaotong Liu, Pierre-Paul De Breuck, Linghui Wang, Gian-Marco Rignanese

Abstract: Machine-learning models have recently encountered enormous success for predicting the properties of materials. These are often trained based on data that present various levels of accuracy, with typically much less high- than low-fidelity data. In order to extract as much information as possible from all available data, we here introduce an approach which aims to improve the quality of the data th… ▽ More Machine-learning models have recently encountered enormous success for predicting the properties of materials. These are often trained based on data that present various levels of accuracy, with typically much less high- than low-fidelity data. In order to extract as much information as possible from all available data, we here introduce an approach which aims to improve the quality of the data through denoising. We investigate the possibilities that it offers in the case of the prediction of the band gap relying on both limited experimental data and density-functional theory relying different exchange-correlation functionals (with an increasing amount of data as the accuracy of the functional decreases). We explore different ways to combine the data into training sequences and analyze the effect of the chosen denoiser. Finally, we analyze the effect of applying the denoising procedure several times until convergence. Our approach provides an improvement over existing methods to exploit multi-fidelity data. △ Less

Submitted 21 April, 2022; originally announced April 2022.

SnP$_2$S$_6$: A Promising Infrared Nonlinear Optical Crystal with Strong Non-Resonant Second Harmonic Generation and Phase-matchability

Authors: Jingyang He, Seng Huat Lee, Francesco Naccarato, Guillaume Brunin, Rui Zu, Yuanxi Wang, Leixin Miao, Huaiyu Wang, Nasim Alem, Geoffroy Hautier, Gian-Marco Rignanese, Zhiqiang Mao, Venkatraman Gopalan

Abstract: High-power infrared laser systems with broadband tunability are of great importance due to their wide range of applications in spectroscopy and free-space communications. These systems require nonlinear optical (NLO) crystals for wavelength up/down conversion using sum/difference frequency generation, respectively. NLO crystals need to satisfy many competing criteria, including large nonlinear opt… ▽ More High-power infrared laser systems with broadband tunability are of great importance due to their wide range of applications in spectroscopy and free-space communications. These systems require nonlinear optical (NLO) crystals for wavelength up/down conversion using sum/difference frequency generation, respectively. NLO crystals need to satisfy many competing criteria, including large nonlinear optical susceptibility, large laser induced damage threshold (LIDT), wide transparency range and phase-matchability. Here, we report bulk single crystals of SnP_2S_6 with a large non-resonant SHG coefficient of d33= 53 pm/V at 1550nm and a large LIDT of 350 GW/cm^2 for femtosecond laser pulses. It also exhibits a broad transparency range from 0.54 μみゅーm to 8.5μみゅーm (bandgap of ~2.3 eV) and can be both Type I and Type II phase-matched. The complete linear and SHG tensors are measured as well as predicted by first principles calculations, and they are in excellent agreement. A proximate double-resonance condition in the electronic band structure for both the fundamental and the SHG light is shown to enhance the non-resonant SHG response. Therefore, SnP2S6 is an outstanding candidate for infrared laser applications. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Assessing the quality of relaxation-time approximations with fully-automated computations of phonon-limited mobilities

Authors: Romain Claes, Guillaume Brunin, Matteo Giantomassi, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: The mobility of carriers, as limited by their scattering with phonons, can now routinely be obtained from first-principles electron-phonon coupling calculations. However, so far, most computations have relied on some form of simplification of the linearized Boltzmann transport equation based on either the self-energy, the momentum- or constant relaxation time approximations. Here, we develop a hig… ▽ More The mobility of carriers, as limited by their scattering with phonons, can now routinely be obtained from first-principles electron-phonon coupling calculations. However, so far, most computations have relied on some form of simplification of the linearized Boltzmann transport equation based on either the self-energy, the momentum- or constant relaxation time approximations. Here, we develop a high-throughput infrastructure and an automatic workflow and we compute 69 phonon-limited mobilities in semiconductors. We compare the results resorting to the approximations with the exact iterative solution. We conclude that the approximate values may deviate significantly from the exact ones and are thus not reliable. Given the minimal computational overhead, our work encourages to rely on this exact iterative solution and warns on the possible inaccuracy of earlier results reported using relaxation time approximations. △ Less

Submitted 25 February, 2022; originally announced February 2022.

Comments: 6 pages, 2 figures

Giant thermoelectric figure of merit in multivalley high-complexity-factor LaSO

Authors: Roberta Farris, Francesco Ricci, Giulio Casu, Diana Dahliah, Geoffroy Hautier, Gian-Marco Rignanese, Vincenzo Fiorentini

Abstract: We report a giant thermoelectric figure of merit $ZT$ (up to 6 at 1100 K) in $n$-doped lanthanum oxysulphate LaSO. Thermoelectric coefficients are computed from ab initio bands within Bloch-Boltzmann theory in an energy-, chemical potential- and temperature-dependent relaxation time approximation. The lattice thermal conductivity is estimated from a model employing the ab initio phonon and Grüneis… ▽ More We report a giant thermoelectric figure of merit $ZT$ (up to 6 at 1100 K) in $n$-doped lanthanum oxysulphate LaSO. Thermoelectric coefficients are computed from ab initio bands within Bloch-Boltzmann theory in an energy-, chemical potential- and temperature-dependent relaxation time approximation. The lattice thermal conductivity is estimated from a model employing the ab initio phonon and Grüneisen-parameter spectrum. The main source of the large $ZT$ is the significant power factor which correlates with a large band complexity factor. We also suggest a possible $n$-type dopant for the material based on ab initio calculations. △ Less

Submitted 27 November, 2021; originally announced November 2021.

Journal ref: Physical Review Materials 5, 125406 (2021)

Indirect light absorption model for highly strained silicon infrared sensors

Authors: Nicolas Roisin, Guillaume Brunin, Gian-Marco Rignanese, Denis Flandre, Jean-Pierre Raskin

Abstract: The optical properties of silicon can be greatly tuned by applying strain and opening new perspectives, particularly in applications where infrared is key. In this work, we use a recent model for the indirect light absorption of silicon and include the effects of tensile and compressive uniaxial strains. The model is based on material properties such as the bandgap, the conduction and valence band… ▽ More The optical properties of silicon can be greatly tuned by applying strain and opening new perspectives, particularly in applications where infrared is key. In this work, we use a recent model for the indirect light absorption of silicon and include the effects of tensile and compressive uniaxial strains. The model is based on material properties such as the bandgap, the conduction and valence band density-of-states effective masses, and the phonon frequencies, which are obtained from first principles including strain up to +2% along the [110] and [111] directions. We show that the limit of absorption can increase from 1.14 (1.09) to 1.35 $μみゅー$m (0.92 eV) under 2% strain and that the absorption increases by a factor of 55 for the zero-strain cutoff wavelength of 1.14 $μみゅー$m when a 2% compressive strain is applied in the [110] direction. We demonstrate that this effect is mainly due to the impact of strain on the electronic bandgaps of silicon, directly followed by the valence band density-of-states effective mass. △ Less

Submitted 2 November, 2021; v1 submitted 17 August, 2021; originally announced August 2021.

Comments: Published in Journal of Applied Physics (4 August 2021)

Journal ref: J. Appl. Phys. 130, 055105 (2021)

Common workflows for computing material properties using different quantum engines

Authors: Sebastiaan P. Huber, Emanuele Bosoni, Marnik Bercx, Jens Bröder, Augustin Degomme, Vladimir Dikan, Kristjan Eimre, Espen Flage-Larsen, Alberto Garcia, Luigi Genovese, Dominik Gresch, Conrad Johnston, Guido Petretto, Samuel Poncé, Gian-Marco Rignanese, Christopher J. Sewell, Berend Smit, Vasily Tseplyaev, Martin Uhrin, Daniel Wortmann, Aliaksandr V. Yakutovich, Austin Zadoks, Pezhman Zarabadi-Poor, Bonan Zhu, Nicola Marzari , et al. (1 additional authors not shown)

Abstract: The prediction of material properties through electronic-structure simulations based on density-functional theory has become routinely common, thanks, in part, to the steady increase in the number and robustness of available simulation packages. This plurality of codes and methods aiming to solve similar problems is both a boon and a burden. While providing great opportunities for cross-verificati… ▽ More The prediction of material properties through electronic-structure simulations based on density-functional theory has become routinely common, thanks, in part, to the steady increase in the number and robustness of available simulation packages. This plurality of codes and methods aiming to solve similar problems is both a boon and a burden. While providing great opportunities for cross-verification, these packages adopt different methods, algorithms, and paradigms, making it challenging to choose, master, and efficiently use any one for a given task. Leveraging recent advances in managing reproducible scientific workflows, we demonstrate how developing common interfaces for workflows that automatically compute material properties can tackle the challenge mentioned above, greatly simplifying interoperability and cross-verification. We introduce design rules for reproducible and reusable code-agnostic workflow interfaces to compute well-defined material properties, which we implement for eleven different quantum engines and use to compute three different material properties. Each implementation encodes carefully selected simulation parameters and workflow logic, making the implementer's expertise of the quantum engine directly available to non-experts. Full provenance and reproducibility of the workflows is guaranteed through the use of the AiiDA infrastructure. All workflows are made available as open-source and come pre-installed with the Quantum Mobile virtual machine, making their use straightforward. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Journal ref: npj Comput Mater 7, 136 (2021)

OPTIMADE, an API for exchanging materials data

Authors: Casper W. Andersen, Rickard Armiento, Evgeny Blokhin, Gareth J. Conduit, Shyam Dwaraknath, Matthew L. Evans, Ádám Fekete, Abhijith Gopakumar, Saulius Gražulis, Andrius Merkys, Fawzi Mohamed, Corey Oses, Giovanni Pizzi, Gian-Marco Rignanese, Markus Scheidgen, Leopold Talirz, Cormac Toher, Donald Winston, Rossella Aversa, Kamal Choudhary, Pauline Colinet, Stefano Curtarolo, Davide Di Stefano, Claudia Draxl, Suleyman Er , et al. (31 additional authors not shown)

Abstract: The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API throug… ▽ More The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API through worked examples on each of the public materials databases that support the full API specification. △ Less

Submitted 25 August, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 11 pages, 1 table

Journal ref: Sci Data 8, 217 (2021)

Robust model benchmarking and bias-imbalance in data-driven materials science: a case study on MODNet

Authors: Pierre-Paul De Breuck, Matthew L. Evans, Gian-Marco Rignanese

Abstract: As the number of novel data-driven approaches to material science continues to grow, it is crucial to perform consistent quality, reliability and applicability assessments of model performance. In this paper, we benchmark the Materials Optimal Descriptor Network (MODNet) method and architecture against the recently released MatBench v0.1, a curated test suite of materials datasets. MODNet is shown… ▽ More As the number of novel data-driven approaches to material science continues to grow, it is crucial to perform consistent quality, reliability and applicability assessments of model performance. In this paper, we benchmark the Materials Optimal Descriptor Network (MODNet) method and architecture against the recently released MatBench v0.1, a curated test suite of materials datasets. MODNet is shown to outperform current leaders on 6 of the 13 tasks, whilst closely matching the current leaders on a further 2 tasks; MODNet performs particularly well when the number of samples is below 10,000. Attention is paid to two topics of concern when benchmarking models. First, we encourage the reporting of a more diverse set of metrics as it leads to a more comprehensive and holistic comparison of model performance. Second, an equally important task is the uncertainty assessment of a model towards a target domain. Significant variations in validation errors can be observed, depending on the imbalance and bias in the training set (i.e., similarity between training and application space). By using an ensemble MODNet model, confidence intervals can be built and the uncertainty on individual predictions can be quantified. Imbalance and bias issues are often overlooked, and yet are important for successful real-world applications of machine learning in materials science and condensed matter. △ Less

Submitted 7 July, 2021; v1 submitted 3 February, 2021; originally announced February 2021.

Journal ref: J. Phys.: Condens. Matter 33 404002 (2021)

CZTS Raman spectra beyond kesterite: a first-principles study

Authors: S. P. Ramkumar, H. P. C. Miranda, X. Gonze, G. -M. Rignanese

Abstract: Cu$_2$ZnSnS$_4$ is an earth-abundant photovoltaic absorber material predicted to provide a sustainable solution for commercial solar applications. One of the main limitations restricting its commercialization is the issue of cation disorder. Raman spectroscopy has been a sought after technique to characterize disorder in CZTS while a clear consensus between theoretical and experimental results is… ▽ More Cu$_2$ZnSnS$_4$ is an earth-abundant photovoltaic absorber material predicted to provide a sustainable solution for commercial solar applications. One of the main limitations restricting its commercialization is the issue of cation disorder. Raman spectroscopy has been a sought after technique to characterize disorder in CZTS while a clear consensus between theoretical and experimental results is yet to be achieved. In the present study, via the virtual crystal approximation, we take into account the progressive nature of Cu-Zn disorder in CZTS: we obtain the phonon frequencies at zone-center within the density functional perturbation theory formalism, and further compute the Raman spectra for the disordered phases, achieving a consensus between theory and experiment. These calculations confirm the presence of complete disorder in Cu-Zn 2$a$, 2$c$ and 2$d$ Wyckoff sites. They also show that the Raman intensities of two prominent $A$ phonon modes characterized by motion of S atoms, also known to be experimentally significant, play a key role in understanding the nature of disorder in CZTS. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 16 pages, 5 figures (+7 figures in the appendix)

Ferroelectricity and multiferroicity in anti-Ruddlesden-Popper structures

Authors: Maxime Markov, Louis Alaerts, Henrique P. C. Miranda, Guido Petretto, Wei Chen, Janine George, Eric Bousquet, Philippe Ghosez, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: Combining ferroelectricity with other properties such as visible light absorption or long-range magnetic order requires the discovery of new families of ferroelectric materials. Here, through the analysis of a high-throughput database of phonon band structures, we identify a new structural family of anti-Ruddlesden-Popper phases A$_4$X$_2$O (A=Ca, Sr, Ba, Eu, X=Sb, P, As, Bi) showing ferroelectric… ▽ More Combining ferroelectricity with other properties such as visible light absorption or long-range magnetic order requires the discovery of new families of ferroelectric materials. Here, through the analysis of a high-throughput database of phonon band structures, we identify a new structural family of anti-Ruddlesden-Popper phases A$_4$X$_2$O (A=Ca, Sr, Ba, Eu, X=Sb, P, As, Bi) showing ferroelectric and anti-ferroelectric behaviors. The discovered ferroelectrics belong to the new class of hyperferroelectrics which polarize even under open-circuit boundary conditions. The polar distortion involves the movement of O anions against apical A cations and is driven by geometric effects resulting from internal chemical strains. Within this new structural family, we show that Eu$_4$Sb$_2$O combines coupled ferromagnetic and ferroelectric order at the same atomic site, a very rare occurrence in materials physics. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: 21 pages, 4 figures with supplemental information attached

MODNet -- accurate and interpretable property predictions for limited materials datasets by feature selection and joint-learning

Authors: Pierre-Paul De Breuck, Geoffroy Hautier, Gian-Marco Rignanese

Abstract: In order to make accurate predictions of material properties, current machine-learning approaches generally require large amounts of data, which are often not available in practice. In this work, an all-round framework is presented which relies on a feedforward neural network, the selection of physically-meaningful features and, when applicable, joint-learning. Next to being faster in terms of tra… ▽ More In order to make accurate predictions of material properties, current machine-learning approaches generally require large amounts of data, which are often not available in practice. In this work, an all-round framework is presented which relies on a feedforward neural network, the selection of physically-meaningful features and, when applicable, joint-learning. Next to being faster in terms of training time, this approach is shown to outperform current graph-network models on small datasets. In particular, the vibrational entropy at 305 K of crystals is predicted with a mean absolute test error of 0.009 meV/K/atom (four times lower than previous studies). Furthermore, joint-learning reduces the test error compared to single-target learning and enables the prediction of multiple properties at once, such as temperature functions. Finally, the selection algorithm highlights the most important features and thus helps understanding the underlying physics. △ Less

Submitted 7 July, 2021; v1 submitted 30 April, 2020; originally announced April 2020.

Comments: 5 pages, 2 figures

Boron phosphide as a \emph{p}-type transparent conductor: optical absorption and transport through electron-phonon coupling

Authors: Viet-Anh Ha, Bora Karasulu, Ryo Maezono, Guillaume Brunin, Joel Basile Varley, Gian-Marco Rignanese, Bartomeu Monserrat, Geoffroy Hautier

Abstract: Boron phosphide has recently been identified as a potential high hole mobility transparent conducting material. This promise arises from its low hole effective masses. However, BP has a relatively small 2 eV indirect band gap which will affect its transparency. In this work, we computationally study both optical absorption across the indirect gap and phonon-limited electronic transport to quantify… ▽ More Boron phosphide has recently been identified as a potential high hole mobility transparent conducting material. This promise arises from its low hole effective masses. However, BP has a relatively small 2 eV indirect band gap which will affect its transparency. In this work, we computationally study both optical absorption across the indirect gap and phonon-limited electronic transport to quantify the potential of boron phosphide as a \emph{p}-type transparent conductor. We find that phonon-mediated indirect optical absorption is weak in the visible spectrum and that the phonon-limited hole mobility is very high (around 900 cm$^2$/Vs) at room temperature. This exceptional mobility comes from a combination of low hole effective mass and very weak scattering by polar phonon modes. We rationalize the weak scattering by the less ionic bonding in boron phosphide compared to oxides. We suggest this could be a general advantage of non-oxides for \emph{p}-type transparent conducting applications. Using our computed properties, we assess the transparent conductor figure of merit of boron phosphide and shows that it exceeds by one order of magnitude that of established \emph{p}-type transparent conductors, confirming the potential of this material. △ Less

Submitted 11 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. Materials 4, 065401 (2020)

Phonon-limited electron mobility in Si, GaAs and GaP with exact treatment of dynamical quadrupoles

Authors: Guillaume Brunin, Henrique Pereira Coutada Miranda, Matteo Giantomassi, Miquel Royo, Massimiliano Stengel, Matthieu J. Verstraete, Xavier Gonze, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense $\mathbf{q}$ meshes using Fourier transforms and ab initio models to… ▽ More We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense $\mathbf{q}$ meshes using Fourier transforms and ab initio models to describe the long-range potentials generated by dipoles and quadrupoles. To reduce significantly the computational cost, we take advantage of crystal symmetries and employ the linear tetrahedron method and double-grid integration schemes, in conjunction with filtering techniques in the Brillouin zone. We report results for the electron mobility in Si, GaAs, and GaP obtained with this new methodology. △ Less

Submitted 21 September, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. B 102, 094308 (2020)

Electron-Phonon Beyond Fröhlich: Dynamical Quadrupoles in Polar and Covalent Solids

Authors: Guillaume Brunin, Henrique Pereira Coutada Miranda, Matteo Giantomassi, Miquel Royo, Massimiliano Stengel, Matthieu J. Verstraete, Xavier Gonze, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: We include the treatment of quadrupolar fields beyond the Fröhlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical… ▽ More We include the treatment of quadrupolar fields beyond the Fröhlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical quadrupoles are not properly treated. △ Less

Submitted 21 September, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. Lett. 125, 136601 (2020)

Photorealistic modelling of metals from first principles

Authors: Gianluca Prandini, Gian-Marco Rignanese, Nicola Marzari

Abstract: The colours of metals have attracted the attention of humanity since ancient times, and coloured metals, in particular gold compounds, have been employed for tools and objects symbolizing the aesthetics of power. In this work we develop a comprehensive framework to obtain the reflectivity and colour of metals, and show that the trends in optical properties and the colours can be predicted by strai… ▽ More The colours of metals have attracted the attention of humanity since ancient times, and coloured metals, in particular gold compounds, have been employed for tools and objects symbolizing the aesthetics of power. In this work we develop a comprehensive framework to obtain the reflectivity and colour of metals, and show that the trends in optical properties and the colours can be predicted by straightforward first-principles techniques based on standard approximations. We apply this to predict reflectivity and colour of several elemental metals and of different types of metallic compounds (intermetallics, solid solutions and heterogeneous alloys), considering mainly binary alloys based on noble metals. We validate the numerical approach through an extensive comparison with experimental data and the photorealistic rendering of known coloured metals. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Journal ref: npj Computational Materials 5, 129 (2019)

Reproducibility in $G_0W_0$ Calculations for Solids

Authors: Tonatiuh Rangel, Mauro Del Ben, Daniele Varsano, Gabriel Antonius, Fabien Bruneval, Felipe H. da Jornada, Michiel J. van Setten, Okan K. Orhan, David D. O'Regan, Andrew Canning, Andrea Ferretti, Andrea Marini, Gian-Marco Rignanese, Jack Deslippe, Steven G. Louie, Jeffrey B. Neaton

Abstract: Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ a… ▽ More Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Journal ref: Comput. Phys. Commun. 255, 107242 (2020)

Computationally-driven, high throughput identification of CaTe and Li$_\textrm{3}$Sb as promising candidates for high mobility $p$-type transparent conducting materials

Authors: Viet-Anh Ha, Guodong Yu, Francesco Ricci, Diana Dahliah, Michiel J. van Setten, Matteo Giantomassi, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: High-performance $p$-type transparent conducting materials (TCMs) must exhibit a rare combination of properties including high mobility, transparency and $p$-type dopability. The development of high-mobility/conductivity $p$-type TCMs is necessary for many applications such as solar cells, or transparent electronic devices. Oxides have been traditionally considered as the most promising chemical s… ▽ More High-performance $p$-type transparent conducting materials (TCMs) must exhibit a rare combination of properties including high mobility, transparency and $p$-type dopability. The development of high-mobility/conductivity $p$-type TCMs is necessary for many applications such as solar cells, or transparent electronic devices. Oxides have been traditionally considered as the most promising chemical space to dig out novel $p$-type TCMs. However, non-oxides might perform better than traditional $p$-type TCMs (oxides) in terms of mobility. We report on a high-throughput (HT) computational search for non-oxide $p$-type TCMs from a large dataset of more than 30,000 compounds which identified CaTe and Li$_\textrm{3}$Sb as very good candidates for high-mobility $p$-type TCMs. From our calculations, both compounds are expected to be $p$-type dopable: intrinsically for Li$_\textrm{3}$Sb while CaTe would require extrinsic doping. Using electron-phonon computations, we estimate hole mobilities at room-temperature to be about 20 and 70 cm$^2$/Vs for CaTe and Li$_\textrm{3}$Sb, respectively. The computed hole mobility for Li$_\textrm{3}$Sb is quite exceptional and comparable with the electron mobility in the best $n$-type TCMs. △ Less

Submitted 13 November, 2018; originally announced November 2018.

Comments: 10 pages, 5 figures

Journal ref: Phys. Rev. Materials 3, 034601 (2019)

Are Small Polarons Always Detrimental to Transparent Conducting Oxides ?

Authors: Guillaume Brunin, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: Transparent conducting oxides (TCOs) are essential to many technologies including solar cells and transparent electronics. The search for high performance n- or p-type TCOs has mainly focused on materials offering transport through band carriers instead of small polarons. In this work, we break this paradigm and demonstrate using well-known physical models that, in certain circumstances, TCOs exhi… ▽ More Transparent conducting oxides (TCOs) are essential to many technologies including solar cells and transparent electronics. The search for high performance n- or p-type TCOs has mainly focused on materials offering transport through band carriers instead of small polarons. In this work, we break this paradigm and demonstrate using well-known physical models that, in certain circumstances, TCOs exhibiting transport by small polarons offer a better combination of transparency and conductivity than materials conducting through band transport. We link this surprising finding to the fundamentally different physics of optical absorption for band and polaronic carriers. Our work rationalizes the good performances of recently emerging small-polaronic Cr-based p-type TCOs such as Sr-doped LaCrO$_3$ and outlines design principles for the development of high-performance TCOs based on transport by small polarons. This opens new avenues for the discovery of high-performance TCOs especially p-type. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. Materials 3, 064602 (2019)

Searching for Materials with High Refractive Index and Wide Band Gap: A First-Principles High-Throughput Study

Authors: Francesco Naccarato, Francesco Ricci, Jin Suntivich, Geoffroy Hautier, Ludger Wirtz, Gian-Marco Rignanese

Abstract: Materials combining both a high refractive index and a wide band gap are of great interest for optoelectronic and sensor applications. However, these two properties are typically described by an inverse correlation with high refractive index appearing in small gap materials and vice-versa. Here, we conduct a first-principles high-throughput study on more than 4000 semiconductors (with a special fo… ▽ More Materials combining both a high refractive index and a wide band gap are of great interest for optoelectronic and sensor applications. However, these two properties are typically described by an inverse correlation with high refractive index appearing in small gap materials and vice-versa. Here, we conduct a first-principles high-throughput study on more than 4000 semiconductors (with a special focus on oxides). Our data confirm the general inverse trend between refractive index and band gap but interesting outliers are also identified. The data are then analyzed through a simple model involving two main descriptors: the average optical gap and the effective frequency. The former can be determined directly from the electronic structure of the compounds, but the latter cannot. This calls for further analysis in order to obtain a predictive model. Nonetheless, it turns out that the negative effect of a large band gap on the refractive index can counterbalanced in two ways: (i) by limiting the difference between the direct band gap and the average optical gap which can be realized by a narrow distribution in energy of the optical transitions and (ii) by increasing the effective frequency which can be achieved through either a high number of transitions from the top of the valence band to the bottom of the conduction or a high average probability for these transitions. Focusing on oxides, we use our data to investigate how the chemistry influences this inverse relationship and rationalize why certain classes of materials would perform better. Our findings can be used to search for new compounds in many optical applications both in the linear and non-linear regime (waveguides, optical modulators, laser, frequency converter, etc.). △ Less

Submitted 25 March, 2019; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: The manuscript is made of 10 pages with 6 figures and 1 table. At the end of the main text, there is an appendix with 2 figures. Supplemental material includes 7 figures and 6 tables. In the db.csv file there are the materials properties we computed in this work

Journal ref: Phys. Rev. Materials 3, 044602 (2019)

An unlikely route to low lattice thermal conductivity: small atoms in a simple layered structure

Authors: Wanyue Peng, Guido Petretto, Gian-Marco Rignanese, Geoffroy Hautier, Alexandra Zevalkink

Abstract: In the design of materials with low lattice thermal conductivity, compounds with high density, low speed of sound, and complexity at either the atomic, nano- or microstructural level are preferred. The layered compound Mg$_3$Sb$_2$ defies these prevailing paradigms, exhibiting lattice thermal conductivity comparable to PbTe and Bi$_2$Te$_3$, despite its low density and simple structure. The excell… ▽ More In the design of materials with low lattice thermal conductivity, compounds with high density, low speed of sound, and complexity at either the atomic, nano- or microstructural level are preferred. The layered compound Mg$_3$Sb$_2$ defies these prevailing paradigms, exhibiting lattice thermal conductivity comparable to PbTe and Bi$_2$Te$_3$, despite its low density and simple structure. The excellent thermoelectric performance ($zT$ $\sim$ 1.5) in $n$-type Mg$_3$Sb$_2$ has thus far been attributed to its multi-valley conduction band, while its anomalous thermal properties have been largely overlooked. To explain the origin of the low lattice thermal conductivity of Mg$_3$Sb$_2$, we have used both experimental methods and ab initio phonon calculations to investigate trends in the elasticity, thermal expansion and anharmonicity of $A$Mg$_2Pn_2$ Zintl compounds with $A$ = Mg, Ca, Yb, and $Pn$ = Sb and Bi. Phonon calculations within the quasi-harmonic approximation reveal large mode Grüneisen parameters in Mg$_3$Sb$_2$ compared with isostructural compounds, in particular in transverse acoustic modes involving shearing of adjacent anionic layers. Measurements of the elastic moduli and sound velocity as a function of temperature using resonant ultrasound spectroscopy provide a window into the softening of the acoustic branches at high temperature, confirming their exceptionally high anharmonicity. We attribute the anomalous thermal behavior of Mg$_3$Sb$_2$ to the diminutive size of Mg, which may be too small for the octahedrally-coordinated site, leading to weak, unstable interlayer Mg-Sb bonding. This suggests more broadly that soft shear modes resulting from undersized cations provide a potential route to achieving low lattice thermal conductivity low-density, earth-abundant materials. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: 11 pages, 9 figures

First-principle study of paraelectric and ferroelectric CsH$_2$PO$_4$ including dispersion forces: stability and related vibrational, dielectric and elastic properties

Authors: Benoit Van Troeye, Michiel Jan van Setten, Matteo Giantomassi, Marc Torrent, Gian-Marco Rignanese, Xavier Gonze

Abstract: Using density functional theory (DFT) and density functional perturbation theory (DFPT), we investigate the stability and response functions of CsH$_2$PO$_4$, a ferroelectric material at low temperature. This material cannot be described properly by the usual (semi-)local approximations within DFT. The long-range e$^-$-e$^-$ correlation needs to be properly taken into account, using, for instance,… ▽ More Using density functional theory (DFT) and density functional perturbation theory (DFPT), we investigate the stability and response functions of CsH$_2$PO$_4$, a ferroelectric material at low temperature. This material cannot be described properly by the usual (semi-)local approximations within DFT. The long-range e$^-$-e$^-$ correlation needs to be properly taken into account, using, for instance, Grimme's DFT-D methods, as investigated in this work. We find that DFT-D3(BJ) performs the best for the members of the dihydrogenated alkali phosphate family (KH$_2$PO$_4$, RbH$_2$PO$_4$, CsH$_2$PO$_4$), leading to experimental lattice parameters reproduced with an average deviation of 0.5 %. With these DFT-D methods, the structural, dielectric, vibrational and mechanical properties of CsH$_2$PO$_4$ are globally in excellent agreement with the available experiments ($<$ 2% MAPE for Raman-active phonons). Our study suggests the possible existence of a new low-temperature phase for CsH$_2$PO$_4$, not yet reported experimentally. Finally, we report the implementation of DFT-D contributions to elastic constants within DFPT. △ Less

Submitted 26 January, 2018; originally announced January 2018.

Comments: This paper was published in Physical Review B the 25 January 2017 (21 pages, 4 figures)

Journal ref: B. Van Troeye, M. J. van Setten, M. Giantomassi, M. Torrent, G.-M. Rignanese, and X. Gonze, Phys. Rev. B 95, 024112 (2017)

High-Throughput Identification of Electrides from all Known Inorganic Materials

Authors: Lee A. Burton, Francesco Ricci, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: In this paper, we present the results of a large-scale, high-throughput computational search for electrides among all known inorganic materials. Analyzing a database of density functional theory results on more than 60,000 compounds, we identify 69 new electride candidates. We report on all these candidates and discuss the structural and chemical factors leading to electride formation. Among these… ▽ More In this paper, we present the results of a large-scale, high-throughput computational search for electrides among all known inorganic materials. Analyzing a database of density functional theory results on more than 60,000 compounds, we identify 69 new electride candidates. We report on all these candidates and discuss the structural and chemical factors leading to electride formation. Among these candidates, our work identifies the first partially-filled 3d transition metal containing electrides Ba3CrN3 and Sr3CrN3; an unexpected finding that contravenes conventional chemistry. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: 5 page manuscript in letter format, 27 page Supplementary Information

The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table

Authors: M. J. van Setten, M. Giantomassi, E. Bousquet, M. J. Verstraete, D. R. Hamann, X. Gonze, G. -M. Rignanese

Abstract: First-principles calculations in crystalline structures are often performed with a planewave basis set. To make the number of basis functions tractable two approximations are usually introduced: core electrons are frozen and the diverging Coulomb potential near the nucleus is replaced by a smoother expression. The norm-conserving pseudopotential was the first successful method to apply these appro… ▽ More First-principles calculations in crystalline structures are often performed with a planewave basis set. To make the number of basis functions tractable two approximations are usually introduced: core electrons are frozen and the diverging Coulomb potential near the nucleus is replaced by a smoother expression. The norm-conserving pseudopotential was the first successful method to apply these approximations in a fully ab initio way. Later on, more efficient and more exact approaches were developed based on the ultrasoft and the projector augmented wave formalisms. These formalisms are however more complex and developing new features in these frameworks is usually more difficult than in the norm-conserving framework. Most of the existing tables of norm- conserving pseudopotentials, generated long ago, do not include the latest developments, are not systematically tested or are not designed primarily for high accuracy. In this paper, we present our PseudoDojo framework for developing and testing full tables of pseudopotentials, and demonstrate it with a new table generated with the ONCVPSP approach. The PseudoDojo is an open source project, building on the AbiPy package, for developing and systematically testing pseudopotentials. At present it contains 7 different batteries of tests executed with ABINIT, which are performed as a function of the energy cutoff. The results of these tests are then used to provide hints for the energy cutoff for actual production calculations. Our final set contains 141 pseudopotentials split into a standard and a stringent accuracy table. In total around 70.000 calculations were performed to test the pseudopotentials. The process of developing the final table led to new insights into the effects of both the core-valence partitioning and the non-linear core corrections on the stability, convergence, and transferability of norm-conserving pseudopotentials. ... △ Less

Submitted 30 October, 2017; v1 submitted 27 October, 2017; originally announced October 2017.

Comments: abstract truncated, 17 pages, 25 figures, 8 tables

Journal ref: Computer Physics Communications 226, 39-54, 2018

Convergence and pitfalls of density functional perturbation theory phonons calculations from a high-throughput perspective

Authors: Guido Petretto, Xavier Gonze, Geoffroy Hautier, Gian-Marco Rignanese

Abstract: The diffusion of large databases collecting different kind of material properties from high-throughput density functional theory calculations has opened new paths in the study of materials science thanks to data mining and machine learning techniques. Phonon calculations have already been employed successfully to predict materials properties and interpret experimental data, e.g. phase stability, f… ▽ More The diffusion of large databases collecting different kind of material properties from high-throughput density functional theory calculations has opened new paths in the study of materials science thanks to data mining and machine learning techniques. Phonon calculations have already been employed successfully to predict materials properties and interpret experimental data, e.g. phase stability, ferroelectricity and Raman spectra, so their availability for a large set of materials will further increase the analytical and predictive power at hand. Moving to a larger scale with density functional perturbation calculations, however, requires the presence of a robust framework to handle this challenging task. In light of this, we automatized the phonon calculation and applied the result to the analysis of the convergence trends for several materials. This allowed to identify and tackle some common problems emerging in this kind of simulations and to lay out the basis to obtain reliable phonon band structures from high-throughput calculations, as well as optimizing the approach to standard phonon simulations. △ Less

Submitted 3 January, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Journal ref: Computational Materials Science Volume 144, March 2018, Pages 331-337

Automation methodologies and large-scale validation for $GW$, towards high-throughput $GW$ calculations

Authors: M. J. van Setten, M. Giantomassi, X. Gonze, G. -M. Rignanese, G. Hautier

Abstract: The search for new materials, based on computational screening, relies on methods that accurately predict, in an automatic manner, total energy, atomic-scale geometries, and other fundamental characteristics of materials. Many technologically important material properties directly stem from the electronic structure of a material, but the usual workhorse for total energies, namely density-functiona… ▽ More The search for new materials, based on computational screening, relies on methods that accurately predict, in an automatic manner, total energy, atomic-scale geometries, and other fundamental characteristics of materials. Many technologically important material properties directly stem from the electronic structure of a material, but the usual workhorse for total energies, namely density-functional theory, is plagued by fundamental shortcomings and errors from approximate exchange-correlation functionals in its prediction of the electronic structure. At variance, the $GW$ method is currently the state-of-the-art {\em ab initio} approach for accurate electronic structure. It is mostly used to perturbatively correct density-functional theory results, but is however computationally demanding and also requires expert knowledge to give accurate results. Accordingly, it is not presently used in high-throughput screening: fully automatized algorithms for setting up the calculations and determining convergence are lacking. In this work we develop such a method and, as a first application, use it to validate the accuracy of $G_0W_0$ using the PBE starting point, and the Godby-Needs plasmon pole model ($G_0W_0^\textrm{GN}$@PBE), on a set of about 80 solids. The results of the automatic convergence study utilized provides valuable insights. Indeed, we find correlations between computational parameters that can be used to further improve the automatization of $GW$ calculations. Moreover, we find that $G_0W_0^\textrm{GN}$@PBE shows a correlation between the PBE and the $G_0W_0^\textrm{GN}$@PBE gaps that is much stronger than that between $GW$ and experimental gaps. However, the $G_0W_0^\textrm{GN}$@PBE gaps still describe the experimental gaps more accurately than a linear model based on the PBE gaps. △ Less

Submitted 25 September, 2017; v1 submitted 21 September, 2017; originally announced September 2017.

Comments: 12 pages, 11 figures

Journal ref: Phys. Rev. B 96, 155207 (2017)

A parallel orbital-updating based plane-wave basis method for electronic structure calculations

Authors: Yan Pan, Xiaoying Dai, Stefano de Gironcoli, Xin-Gao Gong, Gian-Marco Rignanese, Aihui Zhou

Abstract: Motivated by the recently proposed parallel orbital-updating approach in real space method, we propose a parallel orbital-updating based plane-wave basis method for electronic structure calculations, for solving the corresponding eigenvalue problems. In addition, we propose two new modified parallel orbital-updating methods. Compared to the traditional plane-wave methods, our methods allow for two… ▽ More Motivated by the recently proposed parallel orbital-updating approach in real space method, we propose a parallel orbital-updating based plane-wave basis method for electronic structure calculations, for solving the corresponding eigenvalue problems. In addition, we propose two new modified parallel orbital-updating methods. Compared to the traditional plane-wave methods, our methods allow for two-level parallelization, which is particularly interesting for large scale parallelization. Numerical experiments show that these new methods are more reliable and efficient for large scale calculations on modern supercomputers △ Less

Submitted 13 February, 2017; originally announced February 2017.

Many-body correlations and coupling in benzene-dithiol junctions

Authors: Tonatiuh Rangel, Andrea Ferretti, Valerio Olevano, Gian-Marco Rignanese

Abstract: Most theoretical studies of nanoscale transport in molecular junctions rely on the combination of the Landauer formalism with Kohn-Sham density functional theory (DFT) using standard local and semilocal functionals to approximate exchange and correlation effects. In many cases, the resulting conductance is overestimated with respect to experiments. Recent works have demonstrated that this discrepa… ▽ More Most theoretical studies of nanoscale transport in molecular junctions rely on the combination of the Landauer formalism with Kohn-Sham density functional theory (DFT) using standard local and semilocal functionals to approximate exchange and correlation effects. In many cases, the resulting conductance is overestimated with respect to experiments. Recent works have demonstrated that this discrepancy may be reduced when including many-body corrections on top of DFT. Here we study benzene-dithiol (BDT) gold junctions and analyze the effect of many-body perturbation theory (MBPT) on the calculation of the conductance with respect to different bonding geometries. We find that the many-body corrections to the conductance strongly depend on the metal-molecule coupling strength. In the BDT junction with the lowest coupling, many-body corrections reduce the overestimation on the conductance to a factor two, improving the agreement with experiments. In contrast, in the strongest coupling cases, many-body corrections on the conductance are found to be sensibly smaller and standard DFT reveals a valid approach. △ Less

Submitted 16 March, 2017; v1 submitted 20 December, 2016; originally announced December 2016.

Comments: 9 pages, 4 figures

Journal ref: Phys. Rev. B 95, 115137 (2017)

First-principles investigation of the structural, dynamical and dielectric properties of kesterite, stannite and PMCA phases of Cu2ZnSnS4

Authors: S. Poyyapakkam Ramkumar, Y. Gillet, A. Miglio, M. J. van Setten, X. Gonze, G. -M. Rignanese

Abstract: Cu2ZnSnS4 (CZTS) is a promising material as an absorber in photovoltaic applications. The measured efficiency, however, is far from the theoretically predicted value for the known CZTS phases. To improve the understanding of this discrepancy we investigate the structural, dynamical, and dielectric of the three main phases of CZTS (kesterite, stannite, and PMCA) using density functional perturbatio… ▽ More Cu2ZnSnS4 (CZTS) is a promising material as an absorber in photovoltaic applications. The measured efficiency, however, is far from the theoretically predicted value for the known CZTS phases. To improve the understanding of this discrepancy we investigate the structural, dynamical, and dielectric of the three main phases of CZTS (kesterite, stannite, and PMCA) using density functional perturbation theory (DFPT). The effect of the exchange-correlation functional on the computed properties is analyzed. A qualitative agreement of the theoretical Raman spectrum with measurements is observed. However, none of the phases correspond to the experimental spectrum within the error bar that is usually to be expected for DFPT. This corroborates the need to consider cation disorder and other lattice defects extensively in this material. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Journal ref: Phys. Rev. B 94, 224302 (2016)

Influence of the "second gap" on the transparency-conductivity compromise in transparent conducting oxides: an ab initio study

Authors: Viet-Anh Ha, David Waroquiers, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: Transparent conducting oxides (TCOs) are essential to many technologies. These materials are doped (\emph{n}- or \emph{p}-type) oxides with a large enough band gap (ideally $>$3~eV) to ensure transparency. However, the high carrier concentration present in TCOs lead additionally to the possibility for optical transitions from the occupied conduction bands to higher states for \emph{n}-type materia… ▽ More Transparent conducting oxides (TCOs) are essential to many technologies. These materials are doped (\emph{n}- or \emph{p}-type) oxides with a large enough band gap (ideally $>$3~eV) to ensure transparency. However, the high carrier concentration present in TCOs lead additionally to the possibility for optical transitions from the occupied conduction bands to higher states for \emph{n}-type materials and from lower states to the unoccupied valence bands for \emph{p}-type TCOs. The "second gap" formed by these transitions might limit transparency and a large second gap has been sometimes proposed as a design criteria for high performance TCOs. Here, we study the influence of this second gap on optical absorption using \emph{ab initio} computations for several well-known \emph{n}- and \emph{p}-type TCOs. Our work demonstrates that most known \emph{n}-type TCOs do not suffer from second gap absorption in the visible even at very high carrier concentrations. On the contrary, \emph{p}-type oxides show lowering of their optical transmission for high carrier concentrations due to second gap effects. We link this dissimilarity to the different chemistries involved in \emph{n}- versus typical \emph{p}-type TCOs. Quantitatively, we show that second gap effects lead to only moderate loss of transmission (even in p-type TCOs) and suggest that a wide second gap, while beneficial, should not be considered as a needed criteria for a working TCO. △ Less

Submitted 13 March, 2016; originally announced March 2016.

Comments: 6 pages, 4 figures, APS March Meeting

Can molecular projected density-of-states (PDOS) be systematically used in electronic conductance analysis?

Authors: Tonatiuh Rangel, Gian-Marco Rignanese, Valerio Olevano

Abstract: Using benzene-diamine and benzene-dithiol molecular junctions as benchmarks, we investigate the widespread analysis of the quantum transport conductance $\mathcal{G}(εいぷしろん)$ in terms of the projected density of states (PDOS) onto molecular orbitals (MOs). We first consider two different methods for identifying the relevant MOs: 1) diagonalization of the Hamiltonian of the isolated molecule, and 2) dia… ▽ More Using benzene-diamine and benzene-dithiol molecular junctions as benchmarks, we investigate the widespread analysis of the quantum transport conductance $\mathcal{G}(εいぷしろん)$ in terms of the projected density of states (PDOS) onto molecular orbitals (MOs). We first consider two different methods for identifying the relevant MOs: 1) diagonalization of the Hamiltonian of the isolated molecule, and 2) diagonalization of a submatrix of the junction Hamiltonian constructed by considering only basis elements localized on the molecule. We find that these two methods can lead to substantially different MOs and hence PDOS. Furthermore, within Method 1, the PDOS can differ depending on the isolated molecule chosen to represent the molecular junction (e.g. benzene-dithiol or -dithiolate); and, within Method 2, the PDOS depends on the chosen basis set. We show that these differences can be critical when the PDOS is used to provide a physical interpretation of the conductance (especially, when it has small values as it happens typically at zero bias). In this work, we propose a new approach trying to reconcile the two traditional methods. Though some improvements are achieved, the main problems are still unsolved. Our results raise more general questions and doubts on a PDOS-based analysis of the conductance. △ Less

Submitted 27 April, 2015; originally announced April 2015.

Comments: 12 pages, 9 figures

Journal ref: Beilstein J. Nanotechnol. 6, 1247 (2015)

High-Mobility Bismuth-based Transparent P-Type Oxide from High-throughput Material Screening

Authors: Amit Bathia, Geoffroy Hautier, Tan Nilgianskul, Anna Miglio, Gian-Marco Rignanese, Xavier Gonze, Jin Suntivich

Abstract: Transparent oxides are essential building blocks to many technologies, ranging from components in transparent electronics, transparent conductors, to absorbers and protection layers in photovoltaics and photoelectrochemical devices. However, thus far, it has been difficult to develop p-type oxides with wide band gap and high hole mobility; current state-of-art transparent p-type oxides have hole m… ▽ More Transparent oxides are essential building blocks to many technologies, ranging from components in transparent electronics, transparent conductors, to absorbers and protection layers in photovoltaics and photoelectrochemical devices. However, thus far, it has been difficult to develop p-type oxides with wide band gap and high hole mobility; current state-of-art transparent p-type oxides have hole mobility in the range of < 10 cm$^2$/Vs, much lower than their n-type counterparts. Using high-throughput computational screening to guide the discovery of novel oxides with wide band gap and high hole mobility, we report the computational identification and the experimental verification of a bismuth-based double-perovskite oxide that meets these requirements. Our identified candidate, Ba$_2$BiTaO$_6$, has an optical band gap larger than 4 eV and a Hall hole mobility above 30 cm$^2$/Vs. We rationalize this finding with molecular orbital intuitions; Bi$^{3+}$ with filled s-orbitals strongly overlap with the oxygen p, increasing the extent of the metal-oxygen covalency and effectively reducing the valence effective mass, while Ta$^{5+}$ forms a conduction band with low electronegativity, leading to a high band gap beyond the visible range. Our concerted theory-experiment effort points to the growing utility of a data-driven materials discovery and the combination of both informatics and chemical intuitions as a way to discover future technological materials. △ Less

Submitted 14 December, 2014; originally announced December 2014.

Journal ref: Chem. Mater., 2016, 28, pp 30-34

Low-Dimensional Transport and Large Thermoelectric Power Factors in Bulk Semiconductors by Band Engineering of Highly Directional Electronic States

Authors: Daniel I. Bilc, Geoffroy Hautier, David Waroquiers, Gian-Marco Rignanese, Philippe Ghosez

Abstract: Thermoelectrics are promising to address energy issues but their exploitation is still hampered by low efficiencies. So far, much improvement has been achieved by reducing the thermal conductivity but less by maximizing the power factor. The latter imposes apparently conflicting requirements on the band structure: a narrow energy distribution and a low effective mass. Quantum confinement in nanost… ▽ More Thermoelectrics are promising to address energy issues but their exploitation is still hampered by low efficiencies. So far, much improvement has been achieved by reducing the thermal conductivity but less by maximizing the power factor. The latter imposes apparently conflicting requirements on the band structure: a narrow energy distribution and a low effective mass. Quantum confinement in nanostructures or the introduction of resonant states were suggested as possible solutions to this paradox but with limited success. Here, we propose an original approach to fulfill both requirements in bulk semiconductors. It exploits the highly-directional character of some orbitals to engineer the band-structure and produce a type of low-dimensional transport similar to that targeted in nanostructures, while retaining isotropic properties. Using first-principles calculations, the theoretical concept is demonstrated in Fe$_2$YZ Heusler compounds, yielding power factors 4-5 times larger than in classical thermoelectrics at room temperature. Our findings are totally generic and rationalize the search of alternative compounds with a similar behavior. Beyond thermoelectricity, these might be relevant also in the context of electronic, superconducting or photovoltaic applications. △ Less

Submitted 2 April, 2015; v1 submitted 19 May, 2014; originally announced May 2014.

Comments: 6 pages, 2 figures

Journal ref: Phys. Rev. Lett. 114, 136601 (2015)

Structural, electronic, vibrational and dielectric properties of LaBGeO$_5$ from first principles

Authors: R. Shaltaf, H. K. Juwhari, B. Hamad, J. Khalifeh, G. -M. Rignanese, X. Gonze

Abstract: Structural, electronic, vibrational and dielectric properties of LaBGeO$_5$ with the stillwellite structure are determined based on \textit{ab initio} density functional theory. The theoretically relaxed structure is found to agree well with the existing experimental data with a deviation of less than $0.2\%$. Both the density of states and the electronic band structure are calculated, showing fiv… ▽ More Structural, electronic, vibrational and dielectric properties of LaBGeO$_5$ with the stillwellite structure are determined based on \textit{ab initio} density functional theory. The theoretically relaxed structure is found to agree well with the existing experimental data with a deviation of less than $0.2\%$. Both the density of states and the electronic band structure are calculated, showing five distinct groups of valence bands. Furthermore, the Born effective charge, the dielectric permittivity tensors, and the vibrational frequencies at the center of the Brillouin zone are all obtained. Compared to existing model calculations, the vibrational frequencies are found in much better agreement with the published experimental infrared and Raman data, with absolute and relative rms values of 6.04 cm$^{-1}$, and $1.81\%$, respectively. Consequently, numerical values for both the parallel and perpendicular components of the permittivity tensor are established as 3.55 and 3.71 (10.34 and 12.28), respectively, for the high-(low-)frequency limit. △ Less

Submitted 29 November, 2013; originally announced November 2013.

Accuracy of generalized gradient approximation functionals for density functional perturbation theory calculations

Authors: Lianhua He, Fang Liu, Geoffroy Hautier, Micael J. T. Oliveira, Miguel A. L. Marques, Fernando D. Vila, J. J. Rehr, G. -M. Rignanese, Aihui Zhou

Abstract: We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We… ▽ More We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $αあるふぁ$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $αあるふぁ$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant). △ Less

Submitted 18 September, 2013; originally announced September 2013.

Low Hole Effective Mass p-type Transparent Conducting Oxides: Identification and Design Principles

Authors: Geoffroy Hautier, Anna Miglio, Gerbrand Ceder, Gian-Marco Rignanese, Xavier Gonze

Abstract: The development of high performance transparent conducting oxides (TCOs) is critical to many technologies from transparent electronics to solar cells. While n-type TCOs are present in many devices, current p-type TCOs are not largely commercialized as they exhibit much lower carrier mobilities, due to the large hole effective masses of most oxides. Here, we conduct a high-throughput computational… ▽ More The development of high performance transparent conducting oxides (TCOs) is critical to many technologies from transparent electronics to solar cells. While n-type TCOs are present in many devices, current p-type TCOs are not largely commercialized as they exhibit much lower carrier mobilities, due to the large hole effective masses of most oxides. Here, we conduct a high-throughput computational search on thousands of binary and ternary oxides and identify several highly promising compounds displaying exceptionally low hole effective masses (up to an order of magnitude lower than state of the art p-type TCOs) as well as wide band gaps. In addition to the discovery of specific compounds, the chemical rationalization of our findings opens new directions, beyond current Cu-based chemistries, for the design and development of future p-type TCOs. △ Less

Submitted 22 April, 2013; originally announced April 2013.

Journal ref: Nat. Commun., 4, 2292 (2013)

Origin of magnetism and quasiparticles properties in Cr-doped TiO$_2$

Authors: F. Da Pieve, S. Di Matteo, T. Rangel, M. Giantomassi, D. Lamoen, G. -M. Rignanese, X. Gonze

Abstract: Combining LSDA+$U$ and an analysis of superexchange interactions beyond DFT, we describe the magnetic ground states in rutile and anatase Cr-doped TiO$_2$. In parallel, we correct our LSDA+$U$ ground state through GW corrections ($GW$@LSDA+$U$) that reproduce the position of impurity states and the band gaps in satisfying agreement with experiments. Because of the different topological coordinatio… ▽ More Combining LSDA+$U$ and an analysis of superexchange interactions beyond DFT, we describe the magnetic ground states in rutile and anatase Cr-doped TiO$_2$. In parallel, we correct our LSDA+$U$ ground state through GW corrections ($GW$@LSDA+$U$) that reproduce the position of impurity states and the band gaps in satisfying agreement with experiments. Because of the different topological coordinations of Cr-Cr bonds in the ground states of rutile and anatase, superexchange interactions induce either ferromagnetic or antiferromagnetic couplings of Cr ions. In Cr-doped anatase, this interaction leads to a new mechanism which stabilizes a ferromagnetic ground state, in keeping with experimental evidence, without the need to invoke F-center exchange. △ Less

Submitted 28 February, 2013; originally announced February 2013.

Comments: 5<pages, 4 figures

Band widths and gaps from the Tran-Blaha functional : Comparison with many-body perturbation theory

Authors: David Waroquiers, Aurélien Lherbier, Anna Miglio, Martin Stankovski, Samuel Poncé, Micael J. T. Oliveira, Matteo Giantomassi, Gian-Marco Rignanese, Xavier Gonze

Abstract: For a set of ten crystalline materials (oxides and semiconductors), we compute the electronic band structures using the Tran-Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TB09) functional. The band widths and gaps are compared with those from the local-density approximation (LDA) functional, many-body perturbation theory (MBPT), and experiments. At the density-functional theory (DFT) level, TB09 le… ▽ More For a set of ten crystalline materials (oxides and semiconductors), we compute the electronic band structures using the Tran-Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TB09) functional. The band widths and gaps are compared with those from the local-density approximation (LDA) functional, many-body perturbation theory (MBPT), and experiments. At the density-functional theory (DFT) level, TB09 leads to band gaps in much better agreement with experiments than LDA. However, we observe that it globally underestimates, often strongly, the valence (and conduction) band widths (more than LDA). MBPT corrections are calculated starting from both LDA and TB09 eigenenergies and wavefunctions. They lead to a much better agreement with experimental data for band widths. The band gaps obtained starting from TB09 are close to those from quasi-particle self-consistent GW calculations, at a much reduced cost. Finally, we explore the possibility to tune one of the semi-empirical parameters of the TB09 functional in order to obtain simultaneously better band gaps and widths. We find that these requirements are conflicting. △ Less

Submitted 27 February, 2013; originally announced February 2013.

Comments: 18 pages, 16 figures

Journal ref: Phys. Rev. B, American Physical Society, 2013, 87, 075121

The bandstructure of gold from many-body perturbation theory

Authors: T. Rangel, D. Kecik, P. E. Trevisanutto, G. -M. Rignanese, H. Van Swygenhoven, V. Olevano

Abstract: The bandstructure of gold is calculated using many-body perturbation theory (MBPT). Different approximations within the GW approach are considered. Standard single shot G0W0 corrections shift the unoccupied bands up by ~0.2 eV and the first sp-like occupied band down by ~0.4 eV, while leaving unchanged the 5d occupied bands. Beyond G0W0, quasiparticle self-consistency on the wavefunctions lowers t… ▽ More The bandstructure of gold is calculated using many-body perturbation theory (MBPT). Different approximations within the GW approach are considered. Standard single shot G0W0 corrections shift the unoccupied bands up by ~0.2 eV and the first sp-like occupied band down by ~0.4 eV, while leaving unchanged the 5d occupied bands. Beyond G0W0, quasiparticle self-consistency on the wavefunctions lowers the occupied 5d bands by 0.35 eV. Globally, many-body effects achieve an opening of the interband gap (5d-6sp gap) of 0.35 to 0.75 eV approaching the experimental results. Finally, the quasiparticle bandstructure is compared to the one obtained by the widely used HSE (Heyd, Scuseria, and Ernzerhof) hybrid functional. △ Less

Submitted 23 October, 2012; v1 submitted 20 March, 2012; originally announced March 2012.

Journal ref: Rangel, T., Kecik, D., Trevisanutto, P. E., Rignanese, G.- M., Van Swygenhoven, H., and Olevano, V. (2012). Band structure of gold from many-body perturbation theory. Phys. Rev. B, 86(12):125125