-
Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms
Authors:
Sebastian Kokott,
Florian Merz,
Yi Yao,
Christian Carbogno,
Mariana Rossi,
Ville Havu,
Markus Rampp,
Matthias Scheffler,
Volker Blum
Abstract:
Hybrid density functional approximations (DFAs) offer compelling accuracy for ab initio electronic-structure simulations of molecules, nanosystems, and bulk materials, addressing some deficiencies of computationally cheaper, frequently used semilocal DFAs. However, the computational bottleneck of hybrid DFAs is the evaluation of the non-local exact exchange contribution, which is the limiting fact…
▽ More
Hybrid density functional approximations (DFAs) offer compelling accuracy for ab initio electronic-structure simulations of molecules, nanosystems, and bulk materials, addressing some deficiencies of computationally cheaper, frequently used semilocal DFAs. However, the computational bottleneck of hybrid DFAs is the evaluation of the non-local exact exchange contribution, which is the limiting factor for the application of the method for large-scale simulations. In this work, we present a drastically optimized resolution-of-identity-based real-space implementation of the exact exchange evaluation for both non-periodic and periodic boundary conditions in the all-electron code FHI-aims, targeting high-performance CPU compute clusters. The introduction of several new refined Message Passing Interface (MPI) parallelization layers and shared memory arrays according to the MPI-3 standard were the key components of the optimization. We demonstrate significant improvements of memory and performance efficiency, scalability, and workload distribution, extending the reach of hybrid DFAs to simulation sizes beyond ten thousand atoms. As a necessary byproduct of this work, other code parts in FHI-aims have been optimized as well, e.g., the computation of the Hartree potential and the evaluation of the force and stress components. We benchmark the performance and scaling of the hybrid DFA based simulations for a broad range of chemical systems, including hybrid organic-inorganic perovskites, organic crystals and ice crystals with up to 30,576 atoms (101,920 electrons described by 244,608 basis functions).
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Roadmap on Data-Centric Materials Science
Authors:
Stefan Bauer,
Peter Benner,
Tristan Bereau,
Volker Blum,
Mario Boley,
Christian Carbogno,
C. Richard A. Catlow,
Gerhard Dehm,
Sebastian Eibl,
Ralph Ernstorfer,
Ádám Fekete,
Lucas Foppa,
Peter Fratzl,
Christoph Freysoldt,
Baptiste Gault,
Luca M. Ghiringhelli,
Sajal K. Giri,
Anton Gladyshev,
Pawan Goyal,
Jason Hattrick-Simpers,
Lara Kabalan,
Petr Karpov,
Mohammad S. Khorrami,
Christoph Koch,
Sebastian Kokott
, et al. (36 additional authors not shown)
Abstract:
Science is and always has been based on data, but the terms "data-centric" and the "4th paradigm of" materials research indicate a radical change in how information is retrieved, handled and research is performed. It signifies a transformative shift towards managing vast data collections, digital repositories, and innovative data analytics methods. The integration of Artificial Intelligence (AI) a…
▽ More
Science is and always has been based on data, but the terms "data-centric" and the "4th paradigm of" materials research indicate a radical change in how information is retrieved, handled and research is performed. It signifies a transformative shift towards managing vast data collections, digital repositories, and innovative data analytics methods. The integration of Artificial Intelligence (AI) and its subset Machine Learning (ML), has become pivotal in addressing all these challenges. This Roadmap on Data-Centric Materials Science explores fundamental concepts and methodologies, illustrating diverse applications in electronic-structure theory, soft matter theory, microstructure research, and experimental techniques like photoemission, atom probe tomography, and electron microscopy. While the roadmap delves into specific areas within the broad interdisciplinary field of materials science, the provided examples elucidate key concepts applicable to a wider range of topics. The discussed instances offer insights into addressing the multifaceted challenges encountered in contemporary materials research.
△ Less
Submitted 1 May, 2024; v1 submitted 1 February, 2024;
originally announced February 2024.
-
Self-interaction corrected SCAN functional for molecules and solids in the numeric atom-center orbital framework
Authors:
Sheng Bi,
Christian Carbogno,
Igor Ying Zhang,
Matthias Scheffler
Abstract:
Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC metho…
▽ More
Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC method is known for its significant numerical instability. In this study, we introduce a novel constraint that facilitates self-consistent localization of the SIC orbitals in the spirit of Edmiston-Ruedenberg orbitals [Rev. Mod. Phys. 35, 457 (1963)]. Our practical implementation within the all-electron numeric atom-centered orbitals code FHI-aims guarantees efficient and stable convergence of the self-consistent PZ-SIC equations for both molecules and solids. We further demonstrate that our PZ-SIC approach effectively mitigates the SIE in the meta-GGA SCAN functional, significantly improving the accuracy for ionization potentials, charge-transfer energies, and band gaps for a diverse selection of molecules and solids. However, our PZ-SIC method does have its limitations. It can not improve the already accurate SCAN results for properties such as cohesive energies, lattice constants, and bulk modulus in our test sets. This highlights the need for new-generation DFAs with more comprehensive applicability.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
Spontaneous polarization in NaNbO$_{3}$ film on NdGaO$_{3}$ and DyScO$_{3}$ substrates
Authors:
Kisung Kang,
Saud Bin Anooz,
Jutta Schwarzkopf,
Christian Carbogno
Abstract:
Pure NaNbO$_{3}$ is an antiferroelectric material at room temperature that irreversibly transforms to a ferroelectric polar state when subjected to an external electrical field or lattice strain. Experimentally, it has been observed that NaNbO$_{3}$ films grown on NdGaO$_{3}$ exhibit an electrical polarization along the [001]$_{\mathrm{PC}}$ direction, whereas films on DyScO$_{3}$ substrates exhib…
▽ More
Pure NaNbO$_{3}$ is an antiferroelectric material at room temperature that irreversibly transforms to a ferroelectric polar state when subjected to an external electrical field or lattice strain. Experimentally, it has been observed that NaNbO$_{3}$ films grown on NdGaO$_{3}$ exhibit an electrical polarization along the [001]$_{\mathrm{PC}}$ direction, whereas films on DyScO$_{3}$ substrates exhibit a polarization along the [011]$_{\mathrm{PC}}$ direction. These effects have been attributed to the realization of different lattice symmetries in the films due to the incorporation of lattice strain imposed by the use of oxide substrates with different lattice parameters. However, the underlying atomistic mechanisms of the resulting phase symmetry in the films are hardly clear, given that NaNbO$_{3}$ features a diverse and complex phase diagram. In turn, these also impede a straightforward tailoring and optimization of the resulting macroscopic properties on different substrates. To clarify this issue, we perform all-electron first-principles calculations for several potential NaNbO$_{3}$ polymorphs under stress and strain. The computed properties, including the ferroelectric polarization, reveal that an orthorhombic $Pmc2_{1}$ phase is realized on NdGaO$_{3}$ substrates since this is the only phase with an out-of-plane polarization under a compressive strain. Conversely, the monoclinic $Pm$ phase is consistent for the samples grown on DyScO$_{3}$ substrate, since this phase exhibits a spontaneous in-plane polarization along [011]$_{\mathrm{PC}}$ under tensile strain.
△ Less
Submitted 27 November, 2023; v1 submitted 22 November, 2023;
originally announced November 2023.
-
On the Uncertainty Estimates of Equivariant-Neural-Network-Ensembles Interatomic Potentials
Authors:
Shuaihua Lu,
Luca M. Ghiringhelli,
Christian Carbogno,
Jinlan Wang,
Matthias Scheffler
Abstract:
Machine-learning (ML) interatomic potentials (IPs) trained on first-principles datasets are becoming increasingly popular since they promise to treat larger system sizes and longer time scales, compared to the {\em ab initio} techniques producing the training data. Estimating the accuracy of MLIPs and reliably detecting when predictions become inaccurate is key for enabling their unfailing usage.…
▽ More
Machine-learning (ML) interatomic potentials (IPs) trained on first-principles datasets are becoming increasingly popular since they promise to treat larger system sizes and longer time scales, compared to the {\em ab initio} techniques producing the training data. Estimating the accuracy of MLIPs and reliably detecting when predictions become inaccurate is key for enabling their unfailing usage. In this paper, we explore this aspect for a specific class of MLIPs, the equivariant-neural-network (ENN) IPs using the ensemble technique for quantifying their prediction uncertainties. We critically examine the robustness of uncertainties when the ENN ensemble IP (ENNE-IP) is applied to the realistic and physically relevant scenario of predicting local-minima structures in the configurational space. The ENNE-IP is trained on data for liquid silicon, created by density-functional theory (DFT) with the generalized gradient approximation (GGA) for the exchange-correlation functional. Then, the ensemble-derived uncertainties are compared with the actual errors (comparing the results of the ENNE-IP with those of the underlying DFT-GGA theory) for various test sets, including liquid silicon at different temperatures and out-of-training-domain data such as solid phases with and without point defects as well as surfaces. Our study reveals that the predicted uncertainties are generally overconfident and hold little quantitative predictive power for the actual errors.
△ Less
Submitted 31 August, 2023;
originally announced September 2023.
-
Extrapolation to complete basis-set limit in density-functional theory by quantile random-forest models
Authors:
Daniel T. Speckhard,
Christian Carbogno,
Luca Ghiringhelli,
Sven Lubeck,
Matthias Scheffler,
Claudia Draxl
Abstract:
The numerical precision of density-functional-theory (DFT) calculations depends on a variety of computational parameters, one of the most critical being the basis-set size. The ultimate precision is reached with an infinitely large basis set, i.e., in the limit of a complete basis set (CBS). Our aim in this work is to find a machine-learning model that extrapolates finite basis-size calculations t…
▽ More
The numerical precision of density-functional-theory (DFT) calculations depends on a variety of computational parameters, one of the most critical being the basis-set size. The ultimate precision is reached with an infinitely large basis set, i.e., in the limit of a complete basis set (CBS). Our aim in this work is to find a machine-learning model that extrapolates finite basis-size calculations to the CBS limit. We start with a data set of 63 binary solids investigated with two all-electron DFT codes, exciting and FHI-aims, which employ very different types of basis sets. A quantile-random-forest model is used to estimate the total-energy correction with respect to a fully converged calculation as a function of the basis-set size. The random-forest model achieves a symmetric mean absolute percentage error of lower than 25% for both codes and outperforms previous approaches in the literature. Our approach also provides prediction intervals, which quantify the uncertainty of the models' predictions.
△ Less
Submitted 1 June, 2023; v1 submitted 26 March, 2023;
originally announced March 2023.
-
Heat flux for semi-local machine-learning potentials
Authors:
Marcel F. Langer,
Florian Knoop,
Christian Carbogno,
Matthias Scheffler,
Matthias Rupp
Abstract:
The Green-Kubo (GK) method is a rigorous framework for heat transport simulations in materials. However, it requires an accurate description of the potential-energy surface and carefully converged statistics. Machine-learning potentials can achieve the accuracy of first-principles simulations while allowing to reach well beyond their simulation time and length scales at a fraction of the cost. In…
▽ More
The Green-Kubo (GK) method is a rigorous framework for heat transport simulations in materials. However, it requires an accurate description of the potential-energy surface and carefully converged statistics. Machine-learning potentials can achieve the accuracy of first-principles simulations while allowing to reach well beyond their simulation time and length scales at a fraction of the cost. In this paper, we explain how to apply the GK approach to the recent class of message-passing machine-learning potentials, which iteratively consider semi-local interactions beyond the initial interaction cutoff. We derive an adapted heat flux formulation that can be implemented using automatic differentiation without compromising computational efficiency. The approach is demonstrated and validated by calculating the thermal conductivity of zirconium dioxide across temperatures.
△ Less
Submitted 28 March, 2023; v1 submitted 25 March, 2023;
originally announced March 2023.
-
Anharmonicity in Thermal Insulators: An Analysis from First Principles
Authors:
Florian Knoop,
Thomas A. R. Purcell,
Matthias Scheffler,
Christian Carbogno
Abstract:
The anharmonicity of atomic motion limits the thermal conductivity in crystalline solids. However, a microscopic understanding of the mechanisms active in strong thermal insulators is lacking. In this letter, we classify 465 experimentally known materials with respect to their anharmonicity and perform fully anharmonic ab initio Green-Kubo calculations for 58 of them, finding 28 thermal insulators…
▽ More
The anharmonicity of atomic motion limits the thermal conductivity in crystalline solids. However, a microscopic understanding of the mechanisms active in strong thermal insulators is lacking. In this letter, we classify 465 experimentally known materials with respect to their anharmonicity and perform fully anharmonic ab initio Green-Kubo calculations for 58 of them, finding 28 thermal insulators with $κ< 10$ W/mK including 6 with ultralow $κ\lesssim 1$ W/mK. Our analysis reveals that the underlying strong anharmonic dynamics is driven by the exploration of meta-stable intrinsic defect geometries. This is at variance with the frequently applied perturbative approach, in which the dynamics is assumed to evolve around a single stable geometry.
△ Less
Submitted 26 April, 2023; v1 submitted 26 September, 2022;
originally announced September 2022.
-
Ab initio Green-Kubo simulations of heat transport in solids: Method and implementation
Authors:
Florian Knoop,
Matthias Scheffler,
Christian Carbogno
Abstract:
Ab initio Green-Kubo (aiGK) simulations of heat transport in solids allow for assessing lattice thermal conductivity in anharmonic or complex materials from first principles. In this work, we present a detailed account of their practical application and evaluation with an emphasis on noise reduction and finite-size corrections in semiconductors and insulators. To account for such corrections, we p…
▽ More
Ab initio Green-Kubo (aiGK) simulations of heat transport in solids allow for assessing lattice thermal conductivity in anharmonic or complex materials from first principles. In this work, we present a detailed account of their practical application and evaluation with an emphasis on noise reduction and finite-size corrections in semiconductors and insulators. To account for such corrections, we propose strategies in which all necessary numerical parameters are chosen based on the dynamical properties displayed during molecular dynamics simulations in order to minimize manual intervention. This paves the way for applying the aiGK method in semi-automated and high-throughput frameworks. The proposed strategies are presented and demonstrated for computing the lattice thermal conductivity at room temperature in the mildly anharmonic periclase MgO, and for the strongly anharmonic marshite CuI.
△ Less
Submitted 26 April, 2023; v1 submitted 2 September, 2022;
originally announced September 2022.
-
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.
-
Accelerating Materials-Space Exploration for Thermal Insulators by Mapping Materials Properties via Artificial Intelligence
Authors:
Thomas A. R. Purcell,
Matthias Scheffler,
Luca M. Ghiringhelli,
Christian Carbogno
Abstract:
Reliable artificial-intelligence models have the potential to accelerate the discovery of materials with optimal properties for various applications, including superconductivity, catalysis, and thermoelectricity. Advancements in this field are often hindered by the scarcity and quality of available data and the significant effort required to acquire new data. For such applications, reliable surrog…
▽ More
Reliable artificial-intelligence models have the potential to accelerate the discovery of materials with optimal properties for various applications, including superconductivity, catalysis, and thermoelectricity. Advancements in this field are often hindered by the scarcity and quality of available data and the significant effort required to acquire new data. For such applications, reliable surrogate models that help guide materials space exploration using easily accessible materials properties are urgently needed. Here, we present a general, data-driven framework that provides quantitative predictions as well as qualitative rules for steering data creation for all datasets via a combination of symbolic regression and sensitivity analysis. We demonstrate the power of the framework by generating an accurate analytic model for the lattice thermal conductivity using only 75 experimentally measured values. By extracting the most influential material properties from this model, we are then able to hierarchically screen 732 materials and find 80 ultra-insulating materials.
△ Less
Submitted 6 June, 2023; v1 submitted 27 April, 2022;
originally announced April 2022.
-
Electron-phonon coupling in d-electron solids: A temperature dependent study of rutile TiO2 by first-principles theory and two-photon photoemission
Authors:
Honghui Shang,
Adam Argondizzo,
Shijing Tan,
Jin Zhao,
Patrick Rinke,
Christian Carbogno,
Matthias Scheffler,
Hrvoje Petek
Abstract:
Rutile TiO2 is a paradigmatic transition metal oxide with applications in optics, electronics, photocatalysis, etc., that are subject to pervasive electron-phonon interaction. To understand how energies of its electronic bands, and in general semiconductors or metals where the frontier orbitals have a strong d-band character, depend on temperature, we perform a comprehensive theoretical and experi…
▽ More
Rutile TiO2 is a paradigmatic transition metal oxide with applications in optics, electronics, photocatalysis, etc., that are subject to pervasive electron-phonon interaction. To understand how energies of its electronic bands, and in general semiconductors or metals where the frontier orbitals have a strong d-band character, depend on temperature, we perform a comprehensive theoretical and experimental study of the effects of electron-phonon (e-p) interactions. In a two-photon photoemission (2PP) spectroscopy study we observe an unusual temperature dependence of electronic band energies within the conduction band of reduced rutile TiO2, which is contrary to the well understood sp-band semiconductors and points to a so far unexplained dichotomy in how the e-p interactions affect differently the materials where the frontier orbitals are derived from the sp- and d-orbitals. To develop a broadly applicable model, we employ state-of-the-art first-principles calculations that explain how phonons promote interactions between the Ti-3d orbitals of the conduction band within the octahedral crystal field. The characteristics differences in e-p interactions experienced by the Ti 3d-orbitals of rutile TiO2 crystal lattice are contrasted with the more familiar behavior of the Si 2s-orbitals of stishovite SiO2 polymorph, in which the frontier 2s-orbital experience a similar crystal field with the opposite effect...
△ Less
Submitted 8 September, 2020;
originally announced September 2020.
-
Numerical Quality Control for DFT-based Materials Databases
Authors:
Christian Carbogno,
Kristian Sommer Thygesen,
Björn Bieniek,
Claudia Draxl,
Luca M. Ghiringhelli,
Andris Gulans,
Oliver T. Hofmann,
Karsten W. Jacobsen,
Sven Lubeck,
Jens Jørgen Mortensen,
Mikkel Strange,
Elisabeth Wruss,
Matthias Scheffler
Abstract:
Electronic-structure theory is a strong pillar of materials science. Many different computer codes that employ different approaches are used by the community to solve various scientific problems. Still, the precision of different packages has only recently been scrutinized thoroughly, focusing on a specific task, namely selecting a popular density functional, and using unusually high, extremely pr…
▽ More
Electronic-structure theory is a strong pillar of materials science. Many different computer codes that employ different approaches are used by the community to solve various scientific problems. Still, the precision of different packages has only recently been scrutinized thoroughly, focusing on a specific task, namely selecting a popular density functional, and using unusually high, extremely precise numerical settings for investigating 71 monoatomic crystals. Little is known, however, about method- and code-specific uncertainties that arise under numerical settings that are commonly used in practice. We shed light on this issue by investigating the deviations in total and relative energies as a function of computational parameters. Using typical settings for basis sets and k-grids, we compare results for 71 elemental and 63 binary solids obtained by three different electronic-structure codes that employ fundamentally different strategies. On the basis of the observed trends, we propose a simple, analytical model for the estimation of the errors associated with the basis-set incompleteness. We cross-validate this model using ternary systems obtained from the NOMAD Repository and discuss how our approach enables the comparison of the heterogeneous data present in computational materials databases.
△ Less
Submitted 31 January, 2022; v1 submitted 24 August, 2020;
originally announced August 2020.
-
Anharmonicity Measure for Materials
Authors:
Florian Knoop,
Thomas A. R. Purcell,
Matthias Scheffler,
Christian Carbogno
Abstract:
Theoretical frameworks used to qualitatively and quantitatively describe nuclear dynamics in solids are often based on the harmonic approximation. However, this approximation is known to become inaccurate or to break down completely in many modern functional materials. Interestingly, there is no reliable measure to quantify anharmonicity so far. Thus, a systematic classification of materials in te…
▽ More
Theoretical frameworks used to qualitatively and quantitatively describe nuclear dynamics in solids are often based on the harmonic approximation. However, this approximation is known to become inaccurate or to break down completely in many modern functional materials. Interestingly, there is no reliable measure to quantify anharmonicity so far. Thus, a systematic classification of materials in terms of anharmonicity and a benchmark of methodologies that may be appropriate for different strengths of anharmonicity is currently impossible. In this work, we derive and discuss a statistical measure that reliably classifies compounds across temperature regimes and material classes by their "degree of anharmonicity". This enables us to distinguish "harmonic" materials, for which anharmonic effects constitute a small perturbation on top of the harmonic approximation, from strongly "anharmonic" materials, for which anharmonic effects become significant or even dominant and the treatment of anharmonicity in terms of perturbation theory is more than questionable. We show that the analysis of this measure in real and reciprocal space is able to shed light on the underlying microscopic mechanisms, even at conditions close to, e.g., phase transitions or defect formation. Eventually, we demonstrate that the developed approach is computationally efficient and enables rapid high-throughput searches by scanning over a set of several hundred binary solids. The results show that strong anharmonic effects beyond the perturbative limit are not only active in complex materials or close to phase transitions, but already at moderate temperatures in simple binary compounds.
△ Less
Submitted 10 August, 2020; v1 submitted 25 June, 2020;
originally announced June 2020.
-
Accessing the anisotropic non-thermal phonon populations in black phosphorus
Authors:
Hélène Seiler,
Daniela Zahn,
Marios Zacharias,
Patrick Hildebrandt,
Thomas Vasileiadis,
Yoav William Windsor,
Yingpeng Qi,
Christian Carbogno,
Claudia Draxl,
Ralph Ernstorfer,
Fabio Caruso
Abstract:
We combine femtosecond electron diffuse scattering experiments and first-principles calculations of the coupled electron-phonon dynamics to provide a detailed momentum-resolved picture of the ultrafast lattice thermalization in a thin film of black phosphorus. The measurements reveal the emergence of highly anisotropic non-thermal phonon populations which persist for several picoseconds following…
▽ More
We combine femtosecond electron diffuse scattering experiments and first-principles calculations of the coupled electron-phonon dynamics to provide a detailed momentum-resolved picture of the ultrafast lattice thermalization in a thin film of black phosphorus. The measurements reveal the emergence of highly anisotropic non-thermal phonon populations which persist for several picoseconds following excitation of the electrons with a light pulse. Combining ultrafast dynamics simulations based on the time-dependent Boltzmann formalism and calculations of the structure factor, we reproduce the experimental data and identify the vibrational modes primarily responsible for the carrier relaxation via electron-phonon coupling and the subsequent lattice thermalization via phonon-phonon scattering. In particular, we attribute the non-equilibrium lattice dynamics of black phosphorus to highly-anisotropic phonon-assisted scattering processes, which are primarily mediated by high-energy optical phonons. Our approach paves the way towards unravelling and controlling microscopic energy-flow pathways in two-dimensional materials and van der Waals heterostructures, and may also be extended to other non-equilibrium phenomena involving coupled electron-phonon dynamics such as superconductivity, phase transitions or polaron physics.
△ Less
Submitted 27 April, 2021; v1 submitted 23 June, 2020;
originally announced June 2020.
-
A machine learning route between band mapping and band structure
Authors:
Rui Patrick Xian,
Vincent Stimper,
Marios Zacharias,
Maciej Dendzik,
Shuo Dong,
Samuel Beaulieu,
Bernhard Schölkopf,
Martin Wolf,
Laurenz Rettig,
Christian Carbogno,
Stefan Bauer,
Ralph Ernstorfer
Abstract:
Electronic band structure (BS) and crystal structure are the two complementary identifiers of solid state materials. While convenient instruments and reconstruction algorithms have made large, empirical, crystal structure databases possible, extracting quasiparticle dispersion (closely related to BS) from photoemission band mapping data is currently limited by the available computational methods.…
▽ More
Electronic band structure (BS) and crystal structure are the two complementary identifiers of solid state materials. While convenient instruments and reconstruction algorithms have made large, empirical, crystal structure databases possible, extracting quasiparticle dispersion (closely related to BS) from photoemission band mapping data is currently limited by the available computational methods. To cope with the growing size and scale of photoemission data, we develop a pipeline including probabilistic machine learning and the associated data processing, optimization and evaluation methods for band structure reconstruction, leveraging theoretical calculations. The pipeline reconstructs all 14 valence bands of a semiconductor and shows excellent performance on benchmarks and other materials datasets. The reconstruction uncovers previously inaccessible momentum-space structural information on both global and local scales, while realizing a path towards integration with materials science databases. Our approach illustrates the potential of combining machine learning and domain knowledge for scalable feature extraction in multidimensional data.
△ Less
Submitted 15 November, 2022; v1 submitted 20 May, 2020;
originally announced May 2020.
-
Fully Anharmonic, Non-Perturbative Theory of Vibronically Renormalized Electronic Band Structures
Authors:
Marios Zacharias,
Matthias Scheffler,
Christian Carbogno
Abstract:
We develop a first-principles approach for the treatment of vibronic interactions in solids that overcomes the main limitations of state-of-the-art electron-phonon coupling formalisms. In particular, anharmonic effects in the nuclear dynamics are accounted to all orders via ab initio molecular dynamics simulations. This non-perturbative, self-consistent approach evaluates the response of the wave…
▽ More
We develop a first-principles approach for the treatment of vibronic interactions in solids that overcomes the main limitations of state-of-the-art electron-phonon coupling formalisms. In particular, anharmonic effects in the nuclear dynamics are accounted to all orders via ab initio molecular dynamics simulations. This non-perturbative, self-consistent approach evaluates the response of the wave functions along the computed anharmonic trajectory; thus it fully considers the coupling between nuclear and electronic degrees of freedom. We validate and demonstrate the merits of the concept by calculating temperature-dependent, momentum-resolved spectral functions for silicon and the cubic perovskite SrTiO3, a strongly anharmonic material featuring soft modes. In the latter case, our approach reveals that anharmonicity and higher-order vibronic couplings contribute substantially to the electronic-structure at finite-temperatures, noticeably affecting band gaps and effective masses, and hence macroscopic properties such as transport coefficients.
△ Less
Submitted 24 June, 2020; v1 submitted 23 March, 2020;
originally announced March 2020.
-
Parametrically Constrained Geometry Relaxations for High-Throughput Materials Science
Authors:
Maja-Olivia Lenz,
Thomas A. R. Purcell,
David Hicks,
Stefano Curtarolo,
Matthias Scheffler,
Christian Carbogno
Abstract:
Reducing parameter spaces via exploiting symmetries has greatly accelerated and increased the quality of electronic-structure calculations. Unfortunately, many of the traditional methods fail when the global crystal symmetry is broken, even when the distortion is only a slight perturbation (e.g. Jahn-Teller like distortions). Here we introduce a flexible and generalizable parametric relaxation sch…
▽ More
Reducing parameter spaces via exploiting symmetries has greatly accelerated and increased the quality of electronic-structure calculations. Unfortunately, many of the traditional methods fail when the global crystal symmetry is broken, even when the distortion is only a slight perturbation (e.g. Jahn-Teller like distortions). Here we introduce a flexible and generalizable parametric relaxation scheme, and implement it in the all-electron code FHI-aims. This approach utilizes parametric constraints to maintain symmetry at any level. After demonstrating the method's ability to relax metastable structures, we highlight its adaptability and performance over a test set of 359 materials, across thirteen lattice prototypes. Finally we show how these constraints can reduce the number of steps needed to relax local lattice distortions by an order of magnitude. The flexibility of these constraints enables a significant acceleration of the high-throughput searches for novel materials for numerous applications.
△ Less
Submitted 30 October, 2019; v1 submitted 5 August, 2019;
originally announced August 2019.
-
Simultaneous Learning of Several Materials Properties from Incomplete Databases with Multi-Task SISSO
Authors:
Runhai Ouyang,
Emre Ahmetcik,
Christian Carbogno,
Matthias Scheffler,
Luca M. Ghiringhelli
Abstract:
The identification of descriptors of materials properties and functions that capture the underlying physical mechanisms is a critical goal in data-driven materials science. Only such descriptors will enable a trustful and efficient scanning of materials spaces and possibly the discovery of new materials. Recently, the sure-independence screening and sparsifying operator (SISSO) has been introduced…
▽ More
The identification of descriptors of materials properties and functions that capture the underlying physical mechanisms is a critical goal in data-driven materials science. Only such descriptors will enable a trustful and efficient scanning of materials spaces and possibly the discovery of new materials. Recently, the sure-independence screening and sparsifying operator (SISSO) has been introduced and was successfully applied to a number of materials-science problems. SISSO is a compressed-sensing based methodology yielding predictive models that are expressed in form of analytical formulas, built from simple physical properties. These formulas are systematically selected from an immense number (billions or more) of candidates. In this work, we describe a powerful extension of the methodology to a 'multi-task learning' approach, which identifies a single descriptor capturing multiple target materials properties at the same time. This approach is specifically suited for a heterogeneous materials database with scarce or partial data, e.g., in which not all properties are reported for all materials in the training set. As showcase examples, we address the construction of materials-properties maps for the relative stability of octet-binary compounds, considering several crystal phases simultaneously, and the metal/insulator classification of binary materials distributed over many crystal-prototypes.
△ Less
Submitted 3 January, 2019;
originally announced January 2019.
-
Optimizations of the Eigensolvers in the ELPA Library
Authors:
P. Kus,
A. Marek,
S. S. Koecher,
H. -H. Kowalski,
C. Carbogno,
Ch. Scheurer,
K. Reuter,
M. Scheffler,
H. Lederer
Abstract:
The solution of (generalized) eigenvalue problems for symmetric or Hermitian matrices is a common subtask of many numerical calculations in electronic structure theory or materials science. Solving the eigenvalue problem can easily amount to a sizeable fraction of the whole numerical calculation. For researchers in the field of computational materials science, an efficient and scalable solution of…
▽ More
The solution of (generalized) eigenvalue problems for symmetric or Hermitian matrices is a common subtask of many numerical calculations in electronic structure theory or materials science. Solving the eigenvalue problem can easily amount to a sizeable fraction of the whole numerical calculation. For researchers in the field of computational materials science, an efficient and scalable solution of the eigenvalue problem is thus of major importance. The ELPA-library is a well-established dense direct eigenvalue solver library, which has proven to be very efficient and scalable up to very large core counts. In this paper, we describe the latest optimizations of the ELPA-library for new HPC architectures of the Intel Skylake processor family with an AVX-512 SIMD instruction set, or for HPC systems accelerated with recent GPUs. We also describe a complete redesign of the API in a modern modular way, which, apart from a much simpler and more flexible usability, leads to a new path to access system-specific performance optimizations. In order to ensure optimal performance for a particular scientific setting or a specific HPC system, the new API allows the user to influence in straightforward way the internal details of the algorithms and of performance-critical parameters used in the ELPA-library. On top of that, we introduced an autotuning functionality, which allows for finding the best settings in a self-contained automated way. In situations where many eigenvalue problems with similar settings have to be solved consecutively, the autotuning process of the ELPA-library can be done "on-the-fly". Practical applications from materials science which rely on so-called self-consistency iterations can profit from the autotuning. On some examples of scientific interest, simulated with the FHI-aims application, the advantages of the latest optimizations of the ELPA-library are demonstrated.
△ Less
Submitted 3 November, 2018;
originally announced November 2018.
-
Artificial Intelligence for High-Throughput Discovery of Topological Insulators: the Example of Alloyed Tetradymites
Authors:
Guohua Cao,
Runhai Ouyang,
Luca M. Ghiringhelli,
Matthias Scheffler,
Huijun Liu,
Christian Carbogno,
Zhenyu Zhang
Abstract:
Significant advances have been made in predicting new topological materials using high-throughput empirical descriptors or symmetry-based indicators. To date, these approaches have been applied to materials in existing databases, and are severely limited to systems with well-defined symmetries, leaving a much larger materials space unexplored. Using tetradymites as a prototypical class of examples…
▽ More
Significant advances have been made in predicting new topological materials using high-throughput empirical descriptors or symmetry-based indicators. To date, these approaches have been applied to materials in existing databases, and are severely limited to systems with well-defined symmetries, leaving a much larger materials space unexplored. Using tetradymites as a prototypical class of examples, we uncover a novel two-dimensional descriptor by applying an artificial intelligence (AI) based approach for fast and reliable identification of the topological characters of a drastically expanded range of materials, without prior determination of their specific symmetries and detailed band structures. By leveraging this descriptor that contains only the atomic number and electronegativity of the constituent species, we have readily scanned a huge number of alloys in the tetradymite family. Strikingly, nearly half of which are identified to be topological insulators, revealing a much larger territory of the topological materials world. The present work also attests the increasingly important role of such AI-based approaches in modern materials discovery.
△ Less
Submitted 27 February, 2020; v1 submitted 14 August, 2018;
originally announced August 2018.
-
Benefits from using mixed precision computations in the ELPA-AEO and ESSEX-II eigensolver projects
Authors:
Andreas Alvermann,
Achim Basermann,
Hans-Joachim Bungartz,
Christian Carbogno,
Dominik Ernst,
Holger Fehske,
Yasunori Futamura,
Martin Galgon,
Georg Hager,
Sarah Huber,
Thomas Huckle,
Akihiro Ida,
Akira Imakura,
Masatoshi Kawai,
Simone Köcher,
Moritz Kreutzer,
Pavel Kus,
Bruno Lang,
Hermann Lederer,
Valeriy Manin,
Andreas Marek,
Kengo Nakajima,
Lydia Nemec,
Karsten Reuter,
Michael Rippl
, et al. (8 additional authors not shown)
Abstract:
We first briefly report on the status and recent achievements of the ELPA-AEO (Eigenvalue Solvers for Petaflop Applications - Algorithmic Extensions and Optimizations) and ESSEX II (Equipping Sparse Solvers for Exascale) projects. In both collaboratory efforts, scientists from the application areas, mathematicians, and computer scientists work together to develop and make available efficient highl…
▽ More
We first briefly report on the status and recent achievements of the ELPA-AEO (Eigenvalue Solvers for Petaflop Applications - Algorithmic Extensions and Optimizations) and ESSEX II (Equipping Sparse Solvers for Exascale) projects. In both collaboratory efforts, scientists from the application areas, mathematicians, and computer scientists work together to develop and make available efficient highly parallel methods for the solution of eigenvalue problems. Then we focus on a topic addressed in both projects, the use of mixed precision computations to enhance efficiency. We give a more detailed description of our approaches for benefiting from either lower or higher precision in three selected contexts and of the results thus obtained.
△ Less
Submitted 4 June, 2018;
originally announced June 2018.
-
Analysis of Topological Transitions in Two-dimensional Materials by Compressed Sensing
Authors:
Carlos Mera Acosta,
Runhai Ouyang,
Adalberto Fazzio,
Matthias Scheffler,
Luca M. Ghiringhelli,
Christian Carbogno
Abstract:
Quantum spin-Hall insulators (QSHIs), i.e., two-dimensional topological insulators (TIs) with a symmetry-protected band inversion, have attracted considerable scientific interest in recent years. In this work, we have computed the topological Z2 invariant for 220 functionalized honeycomb lattices that are isoelectronic to functionalized graphene. Besides confirming the TI character of well-known m…
▽ More
Quantum spin-Hall insulators (QSHIs), i.e., two-dimensional topological insulators (TIs) with a symmetry-protected band inversion, have attracted considerable scientific interest in recent years. In this work, we have computed the topological Z2 invariant for 220 functionalized honeycomb lattices that are isoelectronic to functionalized graphene. Besides confirming the TI character of well-known materials such as functionalized stanene, our study identifies 45 yet unreported QSHIs. We applied a compressed-sensing approach to identify a physically meaningful descriptor for the Z2 invariant that only depends on the properties of the material's constituent atoms. This enables us to draw a map of materials, in which metals, trivial insulators, and QSHI form distinct regions. This analysis yields fundamental insights in the mechanisms driving topological transitions. The transferability of the identified model is explicitly demonstrated for an additional set of honeycomb lattices with different functionalizations that are not part of the original set of 220 graphene-type materials used to identify the descriptor. In this class, we predict 74 more novel QSHIs that have not been reported in literature yet.
△ Less
Submitted 28 May, 2018;
originally announced May 2018.
-
All-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
Authors:
Honghui Shang,
Nathaniel Raimbault,
Patrick Rinke,
Matthias Scheffler,
Mariana Rossi,
Christian Carbogno
Abstract:
Within density-functional theory, perturbation theory~(PT) is the state-of-the-art formalism for assessing the response to homogeneous electric fields and the associated material properties, e.g., polarizabilities, dielectric constants, and Raman intensities. Here we derive a real-space formulation of PT and present an implementation within the all-electron, numeric atom-centered orbitals electron…
▽ More
Within density-functional theory, perturbation theory~(PT) is the state-of-the-art formalism for assessing the response to homogeneous electric fields and the associated material properties, e.g., polarizabilities, dielectric constants, and Raman intensities. Here we derive a real-space formulation of PT and present an implementation within the all-electron, numeric atom-centered orbitals electronic structure code FHI-aims that allows for massively-parallel calculations. As demonstrated by extensive validation, this allows the rapid computation of accurate response properties of molecules and solids. As an application showcase, we present harmonic and anharmonic Raman spectra, the latter obtained by combining hundreds of thousands of PT calculations with \textit{ab initio} molecular dynamics. By using the PBE exchange-correlation functional with many-body van der Waals corrections, we obtain spectra in good agreement with experiment especially with respect to lineshapes for the isolated paracetamol molecule and two polymorphs of the paracetamol crystal.
△ Less
Submitted 21 June, 2018; v1 submitted 2 March, 2018;
originally announced March 2018.
-
Formation of Vacancies in Si- and Ge-based Clathrates: Role of Electron Localization and Symmetry Breaking
Authors:
Amrita Bhattacharya,
Christian Carbogno,
Bodo Boehme,
Michael Baitinger,
Yuri Grin,
Matthias Scheffler
Abstract:
The formation of framework vacancies in Si- and Ge-based type-I clathrates is studied as function of filling the cages with K and Ba atoms using density-functional theory. Our analysis reveals the relevance of structural disorder, geometric relaxation, electronic saturation, as well as vibrational and configurational entropy. In the Si clathrates we find that vacancies are unstable, but very diffe…
▽ More
The formation of framework vacancies in Si- and Ge-based type-I clathrates is studied as function of filling the cages with K and Ba atoms using density-functional theory. Our analysis reveals the relevance of structural disorder, geometric relaxation, electronic saturation, as well as vibrational and configurational entropy. In the Si clathrates we find that vacancies are unstable, but very differently, in Ge clathrates up to three vacancies per unit cell can be stabilized. This contrasting behavior is largely driven by the different energy gain on populating the electronic vacancy states, which originates from the different degree of localization of the valence orbitals of Si and Ge. This also actuates a qualitatively different atomic relaxation of the framework.
△ Less
Submitted 19 April, 2017;
originally announced April 2017.
-
Lattice Dynamics Calculations based on Density-functional Perturbation Theory in Real Space
Authors:
Honghui Shang,
Christian Carbogno,
Patrick Rinke,
Matthias Scheffler
Abstract:
A real-space formalism for density-functional perturbation theory (DFPT) is derived and applied for the computation of harmonic vibrational properties in molecules and solids. The practical implementation using numeric atom-centered orbitals as basis functions is demonstrated exemplarily for the all-electron Fritz Haber Institute ab initio molecular simulations (FHI-aims) package. The convergence…
▽ More
A real-space formalism for density-functional perturbation theory (DFPT) is derived and applied for the computation of harmonic vibrational properties in molecules and solids. The practical implementation using numeric atom-centered orbitals as basis functions is demonstrated exemplarily for the all-electron Fritz Haber Institute ab initio molecular simulations (FHI-aims) package. The convergence of the calculations with respect to numerical parameters is carefully investigated and a systematic comparison with finite-difference approaches is performed both for finite (molecules) and extended (periodic) systems. Finally, the scaling tests and scalability tests on massively parallel computer systems demonstrate the computational efficiency.
△ Less
Submitted 7 March, 2017; v1 submitted 12 October, 2016;
originally announced October 2016.
-
Ab initio Green-Kubo Approach for the Thermal Conductivity of Solids
Authors:
Christian Carbogno,
Rampi Ramprasad,
Matthias Scheffler
Abstract:
We herein present a first-principles formulation of the Green-Kubo method that allows the accurate assessment of the non-radiative thermal conductivity of solid semiconductors and insulators in equilibrium ab initio molecular dynamics calculations. Using the virial for the nuclei, we propose a unique ab initio definition of the heat flux. Accurate size- and time convergence are achieved within mod…
▽ More
We herein present a first-principles formulation of the Green-Kubo method that allows the accurate assessment of the non-radiative thermal conductivity of solid semiconductors and insulators in equilibrium ab initio molecular dynamics calculations. Using the virial for the nuclei, we propose a unique ab initio definition of the heat flux. Accurate size- and time convergence are achieved within moderate computational effort by a robust, asymptotically exact extrapolation scheme. We demonstrate the capabilities of the technique by investigating the thermal conductivity of extreme high and low heat conducting materials, namely diamond Si and tetragonal ZrO$_2$.
△ Less
Submitted 6 February, 2017; v1 submitted 24 August, 2016;
originally announced August 2016.
-
Towards a Common Format for Computational Material Science Data
Authors:
Luca M. Ghiringhelli,
Christian Carbogno,
Sergey Levchenko,
Fawzi Mohamed,
Georg Huhs,
Martin Lueders,
Micael Oliveira,
Matthias Scheffler
Abstract:
Information and data exchange is an important aspect of scientific progress. In computational materials science, a prerequisite for smooth data exchange is standardization, which means using agreed conventions for, e.g., units, zero base lines, and file formats. There are two main strategies to achieve this goal. One accepts the heterogeneous nature of the community which comprises scientists from…
▽ More
Information and data exchange is an important aspect of scientific progress. In computational materials science, a prerequisite for smooth data exchange is standardization, which means using agreed conventions for, e.g., units, zero base lines, and file formats. There are two main strategies to achieve this goal. One accepts the heterogeneous nature of the community which comprises scientists from physics, chemistry, bio-physics, and materials science, by complying with the diverse ecosystem of computer codes and thus develops "converters" for the input and output files of all important codes. These converters then translate the data of all important codes into a standardized, code-independent format. The other strategy is to provide standardized open libraries that code developers can adopt for shaping their inputs, outputs, and restart files, directly into the same code-independent format. We like to emphasize in this paper that these two strategies can and should be regarded as complementary, if not even synergetic. The main concepts and software developments of both strategies are very much identical, and, obviously, both approaches should give the same final result. In this paper, we present the appropriate format and conventions that were agreed upon by two teams, the Electronic Structure Library (ESL) of CECAM and the NOMAD (NOvel MAterials Discovery) Laboratory, a European Centre of Excellence (CoE). This discussion includes also the definition of hierarchical metadata describing state-of-the-art electronic-structure calculations.
△ Less
Submitted 16 July, 2016;
originally announced July 2016.
-
Signatures of nonadiabatic O2 dissociation at Al(111): First-principles fewest-switches study
Authors:
Christian Carbogno,
Joerg Behler,
Karsten Reuter,
Axel Gross
Abstract:
Recently, spin selection rules have been invoked to explain the discrepancy between measured and calculated adsorption probabilities of molecular oxygen reacting with Al(111). In this work, we inspect the impact of nonadiabatic spin transitions on the dynamics of this system from first principles. For this purpose the motion on two distinct potential-energy surfaces associated to different spin…
▽ More
Recently, spin selection rules have been invoked to explain the discrepancy between measured and calculated adsorption probabilities of molecular oxygen reacting with Al(111). In this work, we inspect the impact of nonadiabatic spin transitions on the dynamics of this system from first principles. For this purpose the motion on two distinct potential-energy surfaces associated to different spin configurations and possible transitions between them are inspected by means of the Fewest Switches algorithm. Within this framework we especially focus on the influence of such spin transitions on observables accessible to molecular beam experiments. On this basis we suggest experimental setups that can validate the occurrence of such transitions and discuss their feasibility.
△ Less
Submitted 2 October, 2009;
originally announced October 2009.
-
Fingerprints for spin-selection rules in the interaction dynamics of O2 at Al(111)
Authors:
Christian Carbogno,
Joerg Behler,
Axel Gross,
Karsten Reuter
Abstract:
We performed mixed quantum-classical molecular dynamics simulations based on first-principles potential-energy surfaces to demonstrate that the scattering of a beam of singlet O2 molecules at Al(111) will enable an unambiguous assessment of the role of spin-selection rules for the adsorption dynamics. At thermal energies we predict a sticking probability that is substantially less than unity, wi…
▽ More
We performed mixed quantum-classical molecular dynamics simulations based on first-principles potential-energy surfaces to demonstrate that the scattering of a beam of singlet O2 molecules at Al(111) will enable an unambiguous assessment of the role of spin-selection rules for the adsorption dynamics. At thermal energies we predict a sticking probability that is substantially less than unity, with the repelled molecules exhibiting characteristic kinetic, vibrational and rotational signatures arising from the non-adiabatic spin transition.
△ Less
Submitted 22 July, 2008;
originally announced July 2008.