Source author record

Luca M. Ghiringhelli

Luca M. Ghiringhelli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Ab initio approach for thermodynamic surface phases with full consideration of anharmonic effects -- the example of hydrogen at Si(100)

A reliable description of surfaces structures in a reactive environment is crucial to understand materials functions. We present a first-principles theory of replica-exchange grand-canonical-ensemble molecular dynamics (REGC-MD) and apply it to evaluate phase equilibria of surfaces in reactive gas-phase environment. We identify the different surface phases and locate phase boundaries including triple and critical points. The approach is demonstrated by addressing open questions for the Si(100) surface in contact with a hydrogen atmosphere. In the range from 300 to 1 000 K, we find 25 distinct thermodynamically stable surface phases, for which we also provide microscopic descriptions. Most of the identified phases, including few order-disorder phase transitions, have not yet been observed experimentally. The REGC-MD-derived phase diagram shows significant, qualitative differences to the description by the state-of-the-art "ab initio atomistic thermodynamics" approach.

preprint2022arXiv

Hierarchical symbolic regression for identifying key physical parameters correlated with bulk properties of perovskites

Symbolic regression identifies key physical parameters describing materials properties by uncovering correlations as nonlinear analytical expressions. However, the pool of expressions grows rapidly with complexity, compromising its efficiency. We tackle this challenge by a hierarchical approach: identified expressions are used as input parameters for obtaining more complex expressions. Crucially, this framework can transfer knowledge among properties, highlighting physical relationships. We demonstrate this strategy by using the Sure-Independence-Screening-and-Sparsifying-Operator (SISSO) approach to identify expressions correlated with the lattice constant and cohesive energy, which are then used to model the bulk modulus of ABO3 perovskites.

preprint2022arXiv

Numerical Quality Control for DFT-based Materials Databases

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.

preprint2022arXiv

TCMI: a non-parametric mutual-dependence estimator for multivariate continuous distributions

The identification of relevant features, i.e., the driving variables that determine a process or the properties of a system, is an essential part of the analysis of data sets with a large number of variables. A mathematical rigorous approach to quantifying the relevance of these features is mutual information. Mutual information determines the relevance of features in terms of their joint mutual dependence to the property of interest. However, mutual information requires as input probability distributions, which cannot be reliably estimated from continuous distributions such as physical quantities like lengths or energies. Here, we introduce total cumulative mutual information (TCMI), a measure of the relevance of mutual dependences that extends mutual information to random variables of continuous distribution based on cumulative probability distributions. TCMI is a non-parametric, robust, and deterministic measure that facilitates comparisons and rankings between feature sets with different cardinality. The ranking induced by TCMI allows for feature selection, i.e., the identification of variable sets that are nonlinear statistically related to a property of interest, taking into account the number of data samples as well as the cardinality of the set of variables. We evaluate the performance of our measure with simulated data, compare its performance with similar multivariate-dependence measures, and demonstrate the effectiveness of our feature-selection method on a set of standard data sets and a typical scenario in materials science.

preprint2021arXiv

Materials genes of heterogeneous catalysis from clean experiments and artificial intelligence

Heterogeneous catalysis is an example of a complex materials function, governed by an intricate interplay of several processes, e.g., the different surface chemical reactions, and the dynamic re-structuring of the catalyst material at reaction conditions. Modelling the full catalytic progression via first-principles statistical mechanics is impractical, if not impossible. Instead, we show here how a tailored artificial-intelligence approach can be applied, even to a small number of materials, to model catalysis and determine the key descriptive parameters ("materials genes") reflecting the processes that trigger, facilitate, or hinder catalyst performance. We start from a consistent experimental set of "clean data", containing nine vanadium-based oxidation catalysts. These materials were synthesized, fully characterized, and tested according to standardized protocols. By applying the symbolic-regression SISSO approach, we identify correlations between the few most relevant materials properties and their reactivity. This approach highlights the underlying physicochemical processes, and accelerates catalyst design.

preprint2020arXiv

Artificial Intelligence for High-Throughput Discovery of Topological Insulators: the Example of Alloyed Tetradymites

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.

preprint2016arXiv

Towards a Common Format for Computational Material Science Data

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.

preprint2015arXiv

Big Data of Materials Science - Critical Role of the Descriptor

Statistical learning of materials properties or functions so far starts with a largely silent, non-challenged step: the choice of the set of descriptive parameters (termed descriptor). However, when the scientific connection between the descriptor and the actuating mechanisms is unclear, causality of the learned descriptor-property relation is uncertain. Thus, trustful prediction of new promising materials, identification of anomalies, and scientific advancement are doubtful. We analyse this issue and define requirements for a suited descriptor. For a classical example, the energy difference of zincblende/wurtzite and rocksalt semiconductors, we demonstrate how a meaningful descriptor can be found systematically.

preprint2015arXiv

Computational Design of Nanoclusters by Property-Based Genetic Algorithms: Tuning the Electronic Properties of (TiO$_2$)$_n$ Clusters

In order to design clusters with desired properties, we have implemented a suite of genetic algorithms tailored to optimize for low total energy, high vertical electron affinity (VEA), and low vertical ionization potential (VIP). Applied to (TiO$_2$)$_n$ clusters, the property-based optimization reveals the underlying structure-property relations and the structural features that may serve as active sites for catalysis. High VEA and low VIP are correlated with the presence of several dangling-O atoms and their proximity, respectively. We show that the electronic properties of (TiO$_2$)$_n$ up to n=20 correlate more strongly with the presence of these structural features than with size.

preprint2015arXiv

Strengthening gold-gold bonds by complexing gold clusters with noble gases

We report an unexpectedly strong and complex chemical bonding of rare-gas atoms to neutral gold clusters. The bonding features are consistently reproduced at different levels of approximation within density-functional theory and beyond: from GGA, through hybrid and double-hybrid functionals, up to renormalized second-order perturbation theory. The main finding is that the adsorption of Ar, Kr, and Xe reduces electron-electron repulsion within gold dimer, causing strengthening of the Au-Au bond. Differently from the dimer, the rare-gas adsorption effects on the gold trimer's geometry and vibrational frequencies are mainly due to electron occupation of the trimer's lowest unoccupied molecular orbital. For the trimer, the theoretical results are also consistent with far-infrared multiple photon dissociation experiments.

preprint2014arXiv

A quantum reactive scattering perspective on electronic nonadiabaticity

Based on quantum reactive-scattering theory, we propose a method for studying the electronic nonadiabaticity in collision processes involving electron-ion rearrangements. We investigate the state-to-state transition probability for electron-ion rearrangements with two comparable approaches. In the first approach the information of the electron is only contained in the ground-state Born-Oppenheimer potential-energy surface, which is the starting point of common reactive-scattering calculations. In the second approach, the electron is explicitly taken into account and included in the calculations at the same level as the ions. Hence, the deviation in the results between the two approaches directly reflects the electronic nonadiabaticity during the collision process. To illustrate the method, we apply it to the well-known proton-transfer model of Shin and Metiu (one electron and three ions), generalized by us in order to allow for reactive scattering channels. It is shown that our explicit electron approach is able to capture electronic nonadiabaticity and the renormalization of the reaction barrier near the classical turning points of the potential in nuclear configuration space. In contrast, system properties near the equilibrium geometry of the asymptotic scattering channels are hardly affected by electronic nonadiabatic effects. We also present an analytical expression for the transition amplitude of the asymmetric proton-transfer model based on the direct evaluation of integrals over the involved Airy functions.

preprint2014arXiv

Efficient ab initio schemes for finding thermodynamically stable and metastable atomic structures: Benchmark of cascade genetic algorithms

A first-principles based methodology for efficiently and accurately finding thermodynamically stable and metastable atomic structures is introduced and benchmarked. The approach is demonstrated for gas-phase metal-oxide clusters in thermodynamic equilibrium with a reactive (oxygen) atmosphere at finite pressure and temperature. It consists of two steps. At first, the potential-energy surface is scanned by means of a global-optimization technique, i.e., a massive-parallel first-principles cascade genetic algorithm for which the choice of all parameters is validated against higher-level methods. In particular, we validate a) the criteria for selection and combination of structures used for the assemblage of new candidate structures, and b) the choice of the exchange-correlation functional. The selection criteria are validated against a fully unbiased method: replica-exchange molecular dynamics. Our choice of the exchange-correlation functional, the van-der-Waals-corrected PBE0 hybrid functional, is justified by comparisons up to highest level currently achievable within density-functional theory, i.e., the renormalized second-order perturbation theory, rPT2. In the second step, the low-energy structures are analyzed by means of ab initio atomistic thermodynamics in order to determine compositions and structures that minimize the Gibbs free energy at given temperature and pressure of the reactive atmosphere.

preprint2013arXiv

Stability and metastability of clusters in a reactive atmosphere: theoretical evidence for unexpected stoichiometries of MgMOx

By applying a genetic algorithm in a cascade approach of increasing accuracy, we calculate the composition and structure of MgMOx clusters at realistic temperatures and oxygen pressures. The stable and metastable systems are identified by ab initio atomistic thermodynamics. We find that small clusters (M <= 5) are in thermodynamic equilibrium when x > M. The non-stoichiometric clusters exhibit peculiar magnetic behavior, suggesting the possibility of tuning magnetic properties by changing environmental pressure and temperature conditions. Furthermore, we show that density-functional theory (DFT) with a hybrid exchange-correlation (xc) functional is needed for predicting accurate phase diagrams of metal-oxide clusters. Neither a (sophisticated) force field nor DFT with (semi)local xc functionals are sufficient for even a qualitative prediction.

preprint2012arXiv

Autocatalytic and cooperatively-stabilized dissociation of water on a stepped platinum surface

Water-metal interfaces are ubiquitous and play a key role in many chemical processes, from catalysis to corrosion. Whereas water adlayers on atomically flat transition metal surfaces have been investigated in depth, little is known about the chemistry of water on stepped surfaces, commonly occurring in realistic situations. Using first-principles simulations we study the adsorption of water on a stepped platinum surface. We find that water adsorbs preferentially at the step edge, forming linear clusters or chains, stabilized by the cooperative effect of chemical bonds with the substrate and hydrogen bonds. In contrast with flat Pt, at steps water molecules dissociate forming mixed hydroxyl/water structures, through an autocatalytic mechanism promoted by hydrogen bonding. Nuclear quantum effects contribute to stabilize partially dissociated cluster and chains. Together with the recently demonstrated attitude of water chains adsorbed on stepped Pt surfaces to transfer protons via thermally activated hopping, these findings candidate these systems as viable proton wires.

preprint2012arXiv

Electronic Energy Functionals: Levy-Lieb principle within the Ground State Path Integral Quantum Monte Carlo

We propose a theoretical/computational protocol based on the use of the Ground State (GS) Path Integral (PI) Quantum Monte Carlo (QMC) for the calculation of the kinetic and Coulomb energy density for a system of $N$ interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the $N-1$-conditional probability density within the framework of the Levy-Lieb constrained search principle. The consequences for the development of energy functionals within the context of Density Functional Theory (DFT) are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the $N-1$-conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an {\it on-the-fly} determination of the energy functional for any system.