BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthias Scheffler

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phphysics.comp-phcond-mat.othermtrl-thcond-mat.str-elBiological PhysicsBiomoleculesphysics.data-ancond-mat.mes-hallcond-mat.softcond-mat.stat-mechInformation TheoryMachine Learningmath.ITphysics.atm-clus

Catalog footprint

What is connected

53works
16topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

53 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

FAIR data enabling new horizons for materials research

The prosperity and lifestyle of our society are very much governed by achievements in condensed matter physics, chemistry and materials science, because new products for sectors such as energy, the environment, health, mobility and information technology (IT) rely largely on improved or even new materials. Examples include solid-state lighting, touchscreens, batteries, implants, drug delivery and many more. The enormous amount of research data produced every day in these fields represents a gold mine of the twenty-first century. This gold mine is, however, of little value if these data are not comprehensively characterized and made available. How can we refine this feedstock; that is, turn data into knowledge and value? For this, a FAIR (findable, accessible, interoperable and reusable) data infrastructure is a must. Only then can data be readily shared and explored using data analytics and artificial intelligence (AI) methods. Making data 'findable and AI ready' (a forward-looking interpretation of the acronym) will change the way in which science is carried out today. In this Perspective, we discuss how we can prepare to make this happen for the field of materials science.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Interface to high-performance periodic coupled-cluster theory calculations with atom-centered, localized basis functions

Coupled cluster (CC) theory is often considered the gold standard of quantum-chemistry. For solids, however, the available software is scarce. We present CC-aims, which can interface ab initio codes with localized atomic orbitals and the CC for solids (CC4S) code by the group of A. Grüneis. CC4S features a continuously growing selection of wave function-based methods including perturbation and CC theory. The CC-aims interface was developed for the FHI-aims code (https://fhi-aims.org) but is implemented such that other codes may use it as a starting point for corresponding interfaces. As CC4S offers treatment of both molecular and periodic systems, the CC-aims interface is a valuable tool, where DFT is either too inaccurate or too unreliable, in theoretical chemistry and materials science alike.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

All-electron periodic $G_0W_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks

We present an all-electron, periodic {\GnWn} implementation within the numerical atomic orbital (NAO) basis framework. A localized variant of the resolution-of-the-identity (RI) approximation is employed to significantly reduce the computational cost of evaluating and storing the two-electron Coulomb repulsion integrals. We demonstrate that the error arising from localized RI approximation can be reduced to an insignificant level by enhancing the set of auxiliary basis functions, used to expand the products of two single-particle NAOs. An efficient algorithm is introduced to deal with the Coulomb singularity in the Brillouin zone sampling that is suitable for the NAO framework. We perform systematic convergence tests and identify a set of computational parameters, which can serve as the default choice for most practical purposes. Benchmark calculations are carried out for a set of prototypical semiconductors and insulators, and compared to independent reference values obtained from an independent $G_0W_0$ implementation based on linearized augmented plane waves (LAPW) plus high-energy localized orbitals (HLOs) basis set, as well as experimental results. With a moderate (FHI-aims \textit{tier} 2) NAO basis set, our $G_0W_0$ calculations produce band gaps that typically lie in between the standard LAPW and the LAPW+HLO results. Complementing \textit{tier} 2 with highly localized Slater-type orbitals (STOs), we find that the obtained band gaps show an overall convergence towards the LAPW+HLO results. The algorithms and techniques developed in this work pave the way for efficient implementations of correlated methods within the NAO framework.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anharmonicity Measure for Materials

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electron-phonon coupling in d-electron solids: A temperature dependent study of rutile TiO2 by first-principles theory and two-photon photoemission

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...

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fully Anharmonic, Non-Perturbative Theory of Vibronically Renormalized Electronic Band Structures

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pentacene and Tetracene Molecules and Films on H/Si(111): Level Alignment from Hybrid Density Functional Theory

The electronic properties of hybrid organic-inorganic semiconductor interfaces depend strongly on the alignment of the electronic carrier levels in the organic/inorganic components. In the present work, we address this energy level alignment from first principles theory for two paradigmatic organic-inorganic semiconductor interfaces, the singlet fission materials tetracene and pentacene on H/Si(111), using all-electron hybrid density functional theory. For isolated tetracene on H/Si(111), a type I-like heterojunction (lowest-energy electron and hole states on Si) is found. For isolated pentacene, the molecular and semiconductor valence band edges are degenerate. For monolayer films, we show how to construct supercell geometries with up to 1,192 atoms, which minimize the strain between the inorganic surface and an organic monolayer film. Based on these models, we predict the formation of type II heterojunctions (electron states on Si, hole-like states on the organic species) for both acenes, indicating that charge separation at the interface between the organic and inorganic components is favored. The paper discusses the steps needed to find appropriate low-energy interface geometries for weakly bonded organic molecules and films on inorganic substrates from first principles, a necessary prerequisite for any computational level alignment prediction.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Length Dependence of Ionization Potentials of Trans-Acetylenes: Internally-Consistent DFT/GW Approach

We follow the evolution of the Ionization Potential (IP) for the paradigmatic quasi-one-dimensional trans-acetylene family of conjugated molecules, from short to long oligomers and to the infinite polymer trans-poly-acetylene (TPA). Our results for short oligomers are very close to experimental available data. We find that the IP varies with oligomer length and converges to the given value for TPA with a smooth, coupled inverse-length-exponential behavior. Our prediction is based on an "internally-consistent" scheme to adjust the exchange mixing parameter $α$ of the PBEh hybrid density functional, so as to obtain a description of the electronic structure consistent with the quasiparticle approximation for the IP. This is achieved by demanding that the corresponding quasiparticle correction, in the GW@PBEh approximation, vanishes for the IP when evaluated at PBEh($α^{ic}$). We find that $α^{ic}$ is also system-dependent and converges with increasing oligomer length, allowing to capture the dependence of IP and other electronic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-Lattice Kinetic Monte Carlo Simulations from First-Principles: Reduction of the Pd(100) Surface Oxide by CO

We present a multi-lattice kinetic Monte Carlo (kMC) approach that efficiently describes the atomistic dynamics of morphological transitions between commensurate structures at crystal surfaces. As an example we study the reduction of a $(\sqrt{5}\times \sqrt{5})R27^{\circ}$ PdO(101) overlayer on Pd(100) in a CO atmosphere. Extensive density-functional theory calculations are used to establish an atomistic pathway for the oxide reduction process. First-principles multi-lattice kMC simulations on the basis of this pathway fully reproduce the experimental temperature dependence of the reduction rate [Fernandes et al., Surf. Sci. 2014, 621, 31-39] and highlight the crucial role of elementary processes special to the boundary between oxide and metal domains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-Consistent Green Function Embedding for Advanced Electronic Structure Methods Based on a Dynamical Mean-Field Concept

We present an embedding scheme for periodic systems that facilitates the treatment of the physically important part (here the unit cell) with advanced electronic-structure methods, that are computationally too expensive for periodic systems. The rest of the periodic system is treated with computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self- consistent manner. Our scheme is based on the concept of dynamical mean-field theory (DMFT) formulated in terms of Green functions. In contrast to the original DMFT formulation for correlated model Hamiltonians, we here consider the unit cell as local embedded cluster in a first-principles way, that includes all electronic degrees of freedom. Our real-space dynamical mean-field embedding (RDMFE) scheme features two nested Dyson equations, one for the embedded cluster and another for the periodic surrounding. The total energy is computed from the resulting Green functions. The performance of our scheme is demonstrated by treating the embedded region with hybrid functionals and many-body perturbation theory in the GW approach for simple bulk systems. The total energy and the density of states converge rapidly with respect to the computational parameters and approach their bulk limit with increasing cluster (i.e., unit cell) size.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Static correlation and electron localization in molecular dimers from the self-consistent RPA and GW approximation

We investigate static correlation and delocalization errors in the self-consistent GW and random-phase approximation (RPA) by studying molecular dissociation of the H_2 and LiH molecules. Although both approximations contain topologically identical diagrams, the non-locality and frequency dependence of the GW self-energy crucially influence the different energy contributions to the total energy as compared to the use of a static local potential in the RPA. The latter leads to significantly larger correlation energies which allow for a better description of static correlation at intermediate bond distances. The substantial error found in GW is further analyzed by comparing spin-restricted and spin-unrestricted calculations. At large but finite nuclear separation their difference gives an estimate of the so-called fractional spin error normally determined only in the dissociation limit. Furthermore, a calculation of the dipole moment of the LiH molecule at dissociation reveals a large delocalization error in GW making the fractional charge error comparable to the RPA. The analyses are supplemented by explicit formulae for the GW Green's function and total energy of a simplified two-level model providing additional insights into the dissociation limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Why graphene growth is very different on the C face than on the Si face of SiC: Insights from surface equilibria and the (3$\times$3)-3C-SiC($\bar{\text{1}}\bar{\text{1}}\bar{\text{1}}$) reconstruction

We address the stability of the surface phases that occur on the C-side of 3C-SiC($\bar{1} \bar{1} \bar{1}$) at the onset of graphene formation. In this growth range, experimental reports reveal a coexistence of several surface phases. This coexistence can be explained by a Si-rich model for the unknown (3$\times$3) reconstruction, the known (2$\times$2)$_{C}$ adatom phase, and the graphene covered (2$\times$2)$_{C}$ phase. By constructing an $ab$ $initio$ surface phase diagram using a van der Waals corrected density functional, we show that the formation of a well defined interface structure like the "buffer-layer" on the Si side is blocked by Si-rich surface reconstructions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A multi-scale approach to the electronic structure of doped semiconductor surfaces

The inclusion of the global effects of semiconductor doping poses a unique challenge for first-principles simulations, because the typically low concentration of dopants renders an explicit treatment intractable. In systems which do not display long-range band bending, a satisfactory remedy is offered by the use of "pseudoatoms", with a fractional nuclear charge matching the bulk doping concentration. However, this alone is not always sufficient for charged surfaces, where the width of the space-charge region (SCR) often exceeds realistic supercell dimensions. One generalization of the pseudoatom approach which overcomes this difficulty relies on the introduction of an artificially high doping level within a slab calculation, in conjunction with a multi-scale electrostatic energy correction. Here, we present an alternative technique that naturally extends the pseudoatom approach while bypassing the need for calculations with an unrealistically high doping level. It is based on the introduction of excess charge, mimicking free charge carriers from the SCR, along with a fixed sheet of counter-charge mimicking the SCR-related field. Self-consistency is obtained by imposing charge conservation and Fermi level equilibration between the bulk, treated semi-classically, and the electronic states of the slab/surface, which are treated quantum-mechanically. The method, which we call CREST - the Charge-Reservoir Electrostatic Sheet Technique - can be used with standard electronic structure codes. We validate the approach using a simple tight-binding model, which allows for comparison of its results with calculations encompassing the full SCR explicitly. We then employ it with density functional theory, where it is used to obtain insights into the electronic structures of the "clean-cleaved" Si(111) surface and its buckled (2x1) reconstruction, at various doping densities.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-Principles Description of Charge Transfer in Donor-Acceptor Compounds from Self-Consistent Many-Body Perturbation Theory

We investigate charge transfer in prototypical molecular donor-acceptor compounds using hybrid density functional theory (DFT) and the GW approximation at the perturbative level (G0W0) and at full self-consistency (sc-GW). For the systems considered here, no charge transfer should be expected at large intermolecular separation according to photoemission experiment and accurate quantum-chemistry calculations. The capability of hybrid exchange-correlation functionals of reproducing this feature depends critically on the fraction of exact exchange $α$, as for small values of $α$ spurious fractional charge transfer is observed between the donor and the acceptor. G0W0 based on hybrid DFT yields the correct alignment of the frontier orbitals for all values of $α$. However, G0W0 has no capacity to alter the ground-state properties of the system, because of its perturbative nature. The electron density in donor-acceptor compounds thus remains incorrect for small $α$ values. In sc-GW, where the Green's function is obtained from the iterative solution of the Dyson equation, the electron density is updated and reflects the correct description of the level alignment at the GW level, demonstrating the importance of self-consistent many-body approaches for the description of ground- and excited-state properties in donor-acceptor systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Accuracy of van der Waals Inclusive Density-Functional Theory Exchange-Correlation Functionals for Ice at Ambient and High Pressures

Density-functional theory (DFT) has been widely used to study water and ice for at least 20 years. However, the reliability of different DFT exchange-correlation (xc) functionals for water remains a matter of considerable debate. This is particularly true in light of the recent development of DFT based methods that account for van der Waals (vdW) dispersion forces. Here, we report a detailed study with several xc functionals (semi-local, hybrid, and vdW inclusive approaches) on ice Ih and six proton ordered phases of ice. Consistent with our previous study [Phys. Rev. Lett. 107, 185701 (2011)] which showed that vdW forces become increasingly important at high pressures, we find here that all vdW inclusive methods considered improve the relative energies and transition pressures of the high-pressure ice phases compared to those obtained with semi-local or hybrid xc functionals. However, we also find that significant discrepancies between experiment and the vdW inclusive approaches remain in the cohesive properties of the various phases, causing certain phases to be absent from the phase diagram. Therefore, room for improvement in the description of water at ambient and high pressures remains and we suggest that because of the stern test the high pressure ice phases pose they should be used in future benchmark studies of simulation methods for water.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bond Breaking and Bond Formation: How Electron Correlation is Captured in Many-Body Perturbation Theory and Density-Functional Theory

For the paradigmatic case of H2-dissociation we compare state-of-the-art many-body perturbation theory (MBPT) in the GW approximation and density-functional theory (DFT) in the exact-exchange plus random-phase approximation for the correlation energy (EX+cRPA). For an unbiased comparison and to prevent spurious starting point effects both approaches are iterated to full self-consistency (i.e. sc-RPA and sc-GW). The exchange-correlation diagrams in both approaches are topologically identical, but in sc-RPA they are evaluated with non-interacting and in sc-GW with interacting Green functions. This has a profound consequence for the dissociation region, where sc-RPA is superior to sc-GW. We argue that for a given diagrammatic expansion, the DFT framework outperforms the many-body framework when it comes to bond-breaking. We attribute this to the difference in the correlation energy rather than the treatment of the kinetic energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Concentration of Vacancies at Metal Oxide Surfaces: Case Study of MgO (100)

We investigate effects of doping on formation energy and concentration of oxygen vacancies at a metal oxide surface, using MgO (100) as an example. Our approach employs density-functional theory, where the performance of the exchange-correlation functional is carefully analyzed, and the functional is chosen according to a fundamental condition on DFT ionization energies. The approach is further validated by CCSD(T) calculations for embedded clusters. We demonstrate that the concentration of oxygen vacancies at a doped oxide surface is largely determined by formation of a macroscopically extended space charge region.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How Cations Change Peptide Structure

Specific interactions between cations and proteins have a strong impact on peptide and protein structure. We here shed light on the nature of the underlying interactions, especially regarding the effects on the polyamide backbone structure. To do so, we compare the conformational ensembles of model peptides in isolation and in the presence of either Li+ or Na+ cations by state-of-the-art density-functional theory (including van der Waals effects) and gas-phase infrared spectroscopy. These monovalent cations have a drastic effect on the local backbone conformation of turn-forming peptides, by disruption of the H bonding networks and the resulting severe distortion of the backbone conformations. In fact, Li+ and Na+ can even have different conformational effects on the same peptide. We also assess the predictive power of current approximate density functionals for peptide-cation systems and compare to results from established protein force fields as well as to high-level quantum chemistry (CCSD(T)).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interface dipoles of organic molecules on Ag(111) in hybrid density-functional theory

We investigate the molecular acceptors 3,4,9,10-perylene-tetracarboxylic acid dianhydride (PTCDA), 2,3,5,6-tetra uoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ), and 4,5,9,10-pyrenetetraone (PYTON) on Ag(111) using densityfunctional theory. For two groups of the HSE(α, ω) family of exchange-correlation functionals (ω= 0 and ω= 0.2Å) we study the isolated components as well as the combined systems as a function of the amount of exact-exchange (α). We find that hybrid functionals favour electron transfer to the adsorbate. Comparing to experimental work-function data, we report for (α) ca. 0.25 a notable but small improvement over (semi)local functionals for the interface dipole. Although Kohn-Sham eigenvalues are only approximate representations of ionization energies, incidentally, at this value also the density of states agrees well with the photoelectron spectra. However, increasing (α) to values for which the energy of the lowest unoccupied molecular orbital matches the experimental electron affinity in the gas phase worsens both the interface dipole and the density of states. Our results imply that semi-local DFT calculations may often be adequate for conjugated organic molecules on metal surfaces and that the much more computationally demanding hybrid functionals yield only small improvements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-consistent GW: All-electron implementation with localized basis functions

This paper describes an all-electron implementation of the self-consistent GW (sc-GW) approach -- i.e. based on the solution of the Dyson equation -- in an all-electron numeric atom-centered orbital (NAO) basis set. We cast Hedin's equations into a matrix form that is suitable for numerical calculations by means of i) the resolution of identity technique to handle 4-center integrals; and ii) a basis representation for the imaginary-frequency dependence of dynamical operators. In contrast to perturbative G0W0, sc-GW provides a consistent framework for ground- and excited-state properties and facilitates an unbiased assessment of the GW approximation. For excited-states, we benchmark sc-GW for five molecules relevant for organic photovoltaic applications: thiophene, benzothiazole, 1,2,5-thiadiazole, naphthalene, and tetrathiafulvalene. At self-consistency, the quasi-particle energies are found to be in good agreement with experiment and, on average, more accurate than G0W0 based on Hartree-Fock (HF) or density-functional theory with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. Based on the Galitskii-Migdal total energy, structural properties are investigated for a set of diatomic molecules. For binding energies, bond lengths, and vibrational frequencies sc-GW and G0W0 achieve a comparable performance, which is, however, not as good as that of exact-exchange plus correlation in the random-phase approximation (EX+cRPA) and its advancement to renormalized second-order perturbation theory (rPT2). Finally, the improved description of dipole moments for a small set of diatomic molecules demonstrates the quality of the sc-GW ground state density.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermodynamic equilibrium conditions of graphene films on SiC

First-principles surface phase diagrams reveal that epitaxial monolayer graphene films on the Si side of 3C-SiC(111) can exist as thermodynamically stable phases in a narrow range of experimentally controllable conditions, defining a path to the highest-quality graphene films. Our calculations are based on a van der Waals corrected density functional. The full, experimentally observed (6 sqrt(3)x 6 sqrt(3))-R30 supercells for zero- to trilayer graphene are essential to describe the correct interface geometries and the relative stability of surface phases and possible defects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Benchmark of GW Methods for Azabenzenes: Is the GW Approximation Good Enough?

Many-body perturbation theory in the GW approximation is a useful method for describing electronic properties associated with charged excitations. A hierarchy of GW methods exists, starting from non-self-consistent G0W0, through partial self-consistency in the eigenvalues (ev-scGW) and in the Green function (scGW0), to fully self-consistent GW (scGW). Here, we assess the performance of these methods for benzene, pyridine, and the diazines. The quasiparticle spectra are compared to photoemission spectroscopy (PES) experiments with respect to all measured particle removal energies and the ordering of the frontier orbitals. We find that the accuracy of the calculated spectra does not match the expectations based on their level of self-consistency. In particular, for certain starting points G0W0 and scGW0 provide spectra in better agreement with the PES than scGW.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assessment of correlation energies based on the random-phase approximation

The random-phase approximation to the ground state correlation energy (RPA) in combination with exact exchange (EX) has brought Kohn-Sham (KS) density functional theory one step closer towards a universal, "general purpose first principles method". In an effort to systematically assess the influence of several correlation energy contributions beyond RPA, this work presents dissociation energies of small molecules and solids, activation energies for hydrogen transfer and non-hydrogen transfer reactions, as well as reaction energies for a number of common test sets. We benchmark EX+RPA and several flavors of energy functionals going beyond it: second-order screened exchange (SOSEX), single excitation (SE) corrections, renormalized single excitation (rSE) corrections, as well as their combinations. Both the single excitation correction as well as the SOSEX contribution to the correlation energy significantly improve upon the notorious tendency of EX+RPA to underbind. Surprisingly, activation energies obtained using EX+RPA based on a KS reference alone are remarkably accurate. RPA+SOSEX+rSE provides an equal level of accuracy for reaction as well as activation energies and overall gives the most balanced performance, which makes it applicable to a wide range of systems and chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Benzene Adsorbed on Metals: Concerted Effect of Covalency and van der Waals Bonding

The adsorption of aromatic molecules on metal surfaces plays a key role in condensed matter physics and functional materials. Depending on the strength of the interaction between the molecule and the surface, the binding is typically classified as either physisorption or chemisorption. Van der Waals (vdW) interactions contribute significantly to the binding in physisorbed systems, but the role of the vdW energy in chemisorbed systems remains unclear. Here we study the interaction of benzene with the (111) surface of transition metals, ranging from weak adsorption (Ag and Au) to strong adsorption (Pt, Pd, Ir, and Rh). When vdW interactions are accurately accounted for, the barrier to adsorption predicted by standard density functional theory (DFT) calculations essentially vanishes, producing a metastable precursor state on Pt and Ir surfaces. Notably, vdW forces contribute more to the binding of covalently bonded benzene than they do when benzene is physisorbed. Comparison to experimental data demonstrates that some of the recently developed methods for including vdW interactions in DFT allow quantitative treatment of both weakly and strongly adsorbed aromatic molecules on metal surfaces, extending the already excellent performance found for gas-phase molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Density-functional Theory for f electron Systems: the α-γ Phase Transition in Cerium

The isostructural α-γ phase transition in cerium is analyzed using density-functional theory with different exchange-correlation functionals, in particular the PBE0 hybrid functional and the exact- exchange plus correlation in the random-phase approximation [(EX+cRPA)@PBE0] approach. We show that the Hartree-Fock exchange part of the hybrid functional actuates two distinct solutions at zero temperature that can be associated with the α and γ phases of cerium. However, despite the relatively good structural and magnetic properties, PBE0 predicts the γ phase to be the stable phase at ambient pressure and zero temperature, in contradiction with low temperature experiments. EX+cRPA reverses the energetic ordering, which emphasizes the importance of correlation for rare- earth systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Impact of vibrational entropy on the stability of unsolvated peptide helices with increasing length

Helices are a key folding motif in protein structure. The question which factors determine helix stability for a given polypeptide or protein is an ongoing challenge. Here we use van der Waals corrected density-functional theory to address a part of this question in a bottom-up approach. We show how intrinsic helical structure is stabilized with length and temperature for a series of experimentally well studied unsolvated alanine based polypeptides, Ac-Alan-LysH+. By exploring extensively the conformational space of these molecules, we find that helices emerge as the preferred structure in the length range n=4-8 not just due to enthalpic factors (hydrogen bonds and their cooperativity, van der Waals dispersion interactions, electrostatics), but importantly also by a vibrational entropic stabilization over competing conformers at room temperature. The stabilization is shown to be due to softer low-frequency vibrational modes in helical conformers than in more compact ones. This observation is corroborated by including anharmonic effects explicitly through \emph{ab initio} molecular dynamics, and generalized by testing different terminations and considering larger helical peptide models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions

Efficient implementations of electronic structure methods are essential for first-principles modeling of molecules and solids. We here present a particularly efficient common framework for methods beyond semilocal density-functional theory, including Hartree-Fock (HF), hybrid density functionals, random-phase approximation (RPA), second-order Møller-Plesset perturbation theory (MP2), and the $GW$ method. This computational framework allows us to use compact and accurate numeric atom-centered orbitals (popular in many implementations of semilocal density-functional theory) as basis functions. The essence of our framework is to employ the "resolution of identity (RI)" technique to facilitate the treatment of both the two-electron Coulomb repulsion integrals (required in all these approaches) as well as the linear density-response function (required for RPA and $GW$). This is possible because these quantities can be expressed in terms of products of single-particle basis functions, which can in turn be expanded in a set of auxiliary basis functions (ABFs). The construction of ABFs lies at the heart of the RI technique, and here we propose a simple prescription for constructing the ABFs which can be applied regardless of whether the underlying radial functions have a specific analytical shape (e.g., Gaussian) or are numerically tabulated. We demonstrate the accuracy of our RI implementation for Gaussian and NAO basis functions, as well as the convergence behavior of our NAO basis sets for the above-mentioned methods. Benchmark results are presented for the ionization energies of 50 selected atoms and molecules from the G2 ion test set as obtained with $GW$ and MP2 self-energy methods, and the G2-I atomization energies as well as the S22 molecular interaction energies as obtained with the RPA method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unified description of ground and excited states of finite systems: the self-consistent GW approach

GW calculations with fully self-consistent G and W -- based on the iterative solution of the Dyson equation -- provide an approach for consistently describing ground and excited states on the same quantum mechanical level. We show that for the systems considered here self-consistent GW reaches the same final Green function regardless of the initial reference state. Self-consistency systematically improves ionization energies and total energies of closed shell systems compared to G_0W_0 based on Hartree-Fock and (semi)local density-functional theory. These improvements also translate to the electron density as exemplified by an improved description of dipole moments and permit us to assess the quality of ground state properties such as bond lengths and vibrational frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Beyond the Random Phase Approximation for the Electron Correlation Energy: The Importance of Single Excitations

The random phase approximation (RPA) for the electron correlation energy, combined with the exact-exchange energy, represents the state-of-the-art exchange-correlation functional within density-functional theory (DFT). However, the standard RPA practice -- evaluating both the exact-exchange and the RPA correlation energy using local or semilocal Kohn-Sham (KS) orbitals -- leads to a systematic underbinding of molecules and solids. Here we demonstrate that this behavior is largely corrected by adding a "single excitation" (SE) contribution, so far not included in the standard RPA scheme. A similar improvement can also be achieved by replacing the non-self-consistent exact-exchange total energy by the corresponding self-consistent Hartree-Fock total energy, while retaining the RPA correlation energy evaluated using Kohn-Sham orbitals. Both schemes achieve chemical accuracy for a standard benchmark set of non-covalent intermolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrogen bonds and van der Waals forces in ice at ambient and high pressures

The first principles approaches, density functional theory (DFT) and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen (H) bonding in ambient and high pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from H bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from DFT functionals which neglect vdW forces are greatly overestimated.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large-scale surface reconstruction energetics of Pt(100) and Au(100) by all-electron DFT

The low-index surfaces of Au and Pt all tend to reconstruct, a fact that is of key importance in many nanostructure, catalytic, and electrochemical applications. Remarkably, some significant questions regarding their structural energies remain even today, in particular for the large-scale quasihexagonal reconstructed (100) surfaces: Rather dissimilar reconstruction energies for Au and Pt in available experiments, and experiment and theory do not match for Pt. We here show by all-electron density-functional theory that only large enough "(5 x N)" approximant supercells capture the qualitative reconstruction energy trend between Au(100) and Pt(100), in contrast to what is often done in the theoretical literature. Their magnitudes are then in fact similar, and closer to the measured value for Pt(100); our calculations achieve excellent agreement with known geometric characteristics and provide direct evidence for the electronic reconstruction driving force.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetism in C or N-doped MgO and ZnO: density-functional study of impurity pairs

It is shown that substitution of C or N for O recently proposed as a way to create ferromagnetism in otherwise nonmagnetic oxide insulators is curtailed by formation of impurity pairs, and the resultant C2 spin=1 dimers as well as the isoelectronic N2^{2+} interact antiferromagneticallly in p-type MgO. For C-doped ZnO, however, we demonstrate using the HSE hybrid functional that a resonance of the spin-polarized C2 ppπ* states with the host conduction band results in a long-range ferromagnetic interaction. Magnetism of open-shell impurity molecules is proposed as a possible route to d0-ferromagnetism in oxide spintronic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Ag-Cu alloy surfaces in an oxidizing environment: a first-principles study

Recent experiments on model catalysts have shown that Ag-Cu alloys have improved selectivity with respect to pure silver for ethylene epoxidation. In this paper we review our first-principles investigations on the (111) surface of this alloy and present new findings on other low index surfaces. We find that, for every surface orientation, the presence of oxygen leads to copper segregation to the surface. Considering the alloy to be in equilibrium with an oxygen atmosphere and accounting for the effect of temperature and pressure, we compute the surface free energy and study the stability of several surface structures. Investigating the dependence of the surface free energy on the surface composition, we construct the phase diagram of the alloy for every surface orientation. Around the temperature, pressure and composition of interest for practical applications, we find that a limited number of structures can be present, including a thin layer of copper oxide on top of the silver surface and copper-free structures. Different surface orientations show a very similar behavior and in particular a single layer with CuO stoichiometry, significantly distorted when compared to a layer of bulk CuO, has a wide range of stability for all orientations. Our results are consistent with, and help explain, recent experimental measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Examination of the concept of degree of rate control by first-principles kinetic Monte Carlo simulations

The conceptual idea of degree of rate control (DRC) approaches is to identify the "rate limiting step" in a complex reaction network by evaluating how the overall rate of product formation changes when a small change is made in one of the kinetic parameters. We examine two definitions of this concept by applying it to first-principles kinetic Monte Carlo simulations of the CO oxidation at RuO2(110). Instead of studying experimental data we examine simulations, because in them we know the surface structure, reaction mechanism, the rate constants, the coverage of the surface and the turn-over frequency at steady state. We can test whether the insights provided by the DRC are in agreement with the results of the simulations thus avoiding the uncertainties inherent in a comparison with experiment. We find that the information provided by using the DRC is non-trivial: It could not have been obtained from the knowledge of the reaction mechanism and of the magnitude of the rate constants alone. For the simulations the DRC provides furthermore guidance as to which aspects of the reaction mechanism should be treated accurately and which can be studied by less accurate and more efficient methods. We therefore conclude that a sensitivity analysis based on the DRC is a useful tool for understanding the propagation of errors from the electronic structure calculations to the statistical simulations in first-principles kinetic Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Non-Adiabatic Potential-Energy Surfaces by Constrained Density-Functional Theory

Non-adiabatic effects play an important role in many chemical processes. In order to study the underlying non-adiabatic potential-energy surfaces (PESs), we present a locally-constrained density-functional theory approach, which enables us to confine electrons to sub-spaces of the Hilbert space, e.g. to selected atoms or groups of atoms. This allows to calculate non-adiabatic PESs for defined charge and spin states of the chosen subsystems. The capability of the method is demonstrated by calculating non-adiabatic PESs for the scattering of a sodium and a chlorine atom, for the interaction of a chlorine molecule with a small metal cluster, and for the dissociation of an oxygen molecule at the Al(111) surface.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Ab initio atomistic thermodynamics and statistical mechanics of surface properties and functions

Previous and present "academic" research aiming at atomic scale understanding is mainly concerned with the study of individual molecular processes possibly underlying materials science applications. Appealing properties of an individual process are then frequently discussed in terms of their direct importance for the envisioned material function, or reciprocally, the function of materials is somehow believed to be understandable by essentially one prominent elementary process only. What is often overlooked in this approach is that in macroscopic systems of technological relevance typically a large number of distinct atomic scale processes take place. Which of them are decisive for observable system properties and functions is then not only determined by the detailed individual properties of each process alone, but in many, if not most cases also the interplay of all processes, i.e. how they act together, plays a crucial role. For a "predictive materials science modeling with microscopic understanding", a description that treats the statistical interplay of a large number of microscopically well-described elementary processes must therefore be applied. Modern electronic structure theory methods such as DFT have become a standard tool for the accurate description of individual molecular processes. Here, we discuss the present status of emerging methodologies which attempt to achieve a (hopefully seamless) match of DFT with concepts from statistical mechanics or thermodynamics, in order to also address the interplay of the various molecular processes. The new quality of, and the novel insights that can be gained by, such techniques is illustrated by how they allow the description of crystal surfaces in contact with realistic gas-phase environments.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Island morphology and adatom self-diffusion on Pt(111)

The results of a density-functional-theory study of the formation energies of (100)- and (111)-faceted steps on the Pt(111) surface, as well as of the barrier for diffusion of an adatom on the flat surface, are presented. The step formation energies are found to be in a ratio of 0.88 in favour of the (111)-faceted step, in excellent agreement with experiment; the equilibrium shape of islands should therefore clearly be non-hexagonal. The origin of the difference between the two steps is discussed in terms of the release of stress at the surface through relaxation. For the diffusion barrier, we also find relaxation to be important, leading to a 20% decrease of its energy. The value we obtain, 0.33 eV, however remains higher than available experimental data; possible reasons for this discrepancy are discussed. We find the ratio of step formation energies and the diffusion barrier to be the same whether using the local-density approximation or the generalized-gradient approximation for the exchange-and-correlation energy.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

A Model for the Thermal Expansion of Ag(111) and other Metal Surfaces

We develop a model to study the thermal expansion of surfaces, wherein phonon frequencies are obtained from ab initio total energy calculations. Anharmonic effects are treated exactly in the direction normal to the surface, and within a quasiharmonic approximation in the plane of the surface. We apply this model to the Ag(111) and Al(111) surfaces, and find that our calculations reproduce the experimental observation of a large and anomalous increase in the surface thermal expansion of Ag(111) at high temperatures [P. Statiris, H.C. Lu and T. Gustafsson, Phys. Rev. Lett. 72, 3574 (1994)]. Surprisingly, we find that this increase can be attributed to a rapid softening of the in-plane phonon frequencies, rather than due to the anharmonicity of the out-of-plane surface phonon modes. This provides evidence for a new mechanism for the enhancement of surface anharmonicity. A comparison with Al(111) shows that the two surfaces behave quite differently, with no evidence for such anomalous behavior on Al(111).

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Steering and isotope effects in the dissociative adsorption of H_2/Pd(100)

The interaction of hydrogen with many transition metal surfaces is characterized by a coexistence of activated with non-activated paths to adsorption with a broad distribution of barrier heights. By performing six-dimensional quantum dynamical and classical molecular dynamics calculations using the same potential energy surface derived from ``ab initio'' calculations for the system H_2/Pd(100) we show that these features of the potential energy surface lead to strong steering effects in the dissociative adsorption dynamics. The adsorption dynamics shows only a small isotope effect which is purely due to the quantum nature of hydrogen.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Localized Excitons and Breaking of Chemical Bonds at III-V (110) Surfaces

Electron-hole excitations in the surface bands of GaAs(110) are analyzed using constrained density-functional theory calculations. The results show that Frenkel-type autolocalized excitons are formed. The excitons induce a local surface unrelaxation which results in a strong exciton-exciton attraction and makes complexes of two or three electron-hole pairs more favorable than separate excitons. In such microscopic exciton "droplets" the electron density is mainly concentrated in the dangling orbital of a surface Ga atom whereas the holes are distributed over the bonds of this atom to its As neighbors thus weakening the bonding to the substrate. This finding suggests the microscopic mechanism of a laser-induced emission of neutral Ga atoms from GaAs and GaP (110) surfaces.

0 citations0 reviewsNo code linkedOpen access