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

Karsten Reuter

Karsten Reuter 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.mes-hallcond-mat.othercond-mat.stat-mechcond-mat.dis-nnphysics.atm-clus

Catalog footprint

What is connected

50works
8topics
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

50 published item(s)

preprint2022arXiv

Machine-learning Based Screening of Lead-free Halide Double Perovskites for Photovoltaic Applications

Lead-free halide double perovskites are promising stable and non-toxic alternatives to methylammonium lead iodide in the field of photovoltaics. In this context, the most commonly used double perovskite is Cs$_2$AgBiBr$_6$, due to its favorable charge transport properties. However, the maximum power conversion efficiency obtained for this material does not exceed 3\%, as a consequence of its wide indirect gap and its intrinsic and extrinsic defects. On the other hand, the materials space that arises from the substitution of different elements in the 4 lattice sites of this structure is large and still mostly unexplored. In this work a neural network is used to predict the band gap of double perovskites from an initial space of 7056 structures and select candidates suitable for visible light absorption. Successive hybrid DFT calculations are used to evaluate the thermodynamic stability, the power conversion efficiency and the effective masses of the selected compounds, and to propose novel potential solar absorbers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predicting binding motifs of complex adsorbates using machine learning with a physics-inspired graph representation

Computational screening in heterogeneous catalysis relies increasingly on machine learning models for predicting key input parameters due to the high cost of computing these directly using first-principles methods. This becomes especially relevant when considering complex materials spaces, e.g. alloys, or complex reaction mechanisms with adsorbates that may exhibit bi- or higher-dentate adsorption motifs. Here we present a data-efficient approach to the prediction of binding motifs and associated adsorption enthalpies of complex adsorbates at transition metals (TMs) and their alloys based on a customized Wasserstein Weisfeiler-Lehman graph kernel and Gaussian Process Regression. The model shows good predictive performance, not only for the elemental TMs on which it was trained, but also for an alloy based on these TMs. Furthermore, incorporation of minimal new training data allows for predicting an out-of-domain TM. We believe the model may be useful in active learning approaches, for which we present an ensemble uncertainty estimation approach.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Selectivity Trends and Role of Adsorbate-Adsorbate Interactions in CO Hydrogenation on Rhodium Catalysts

Predictive-quality computational modeling of heterogeneously catalyzed reactions has emerged as an important tool for the analysis and assessment of activity and activity trends. In contrast, more subtle selectivities and selectivity trends still pose a significant challenge to prevalent microkinetic modeling approaches that typically employ a mean-field approximation (MFA). Here, we focus on CO hydrogenation on Rh catalysts with the possible products methane, acetaldehyde, ethanol and water. This reaction has already been subject to a number of experimental and theoretical studies with conflicting views on the factors controlling activity and selectivity towards the more valuable higher oxygenates. Using accelerated first-principles kinetic Monte Carlo (KMC) simulations and explicitly and systematically accounting for adsorbate-adsorbate interactions through a cluster expansion approach, we model the reaction on the low-index Rh(111) and stepped Rh(211) surfaces. We find that the Rh(111) facet is selective towards methane, while the Rh(211) facet exhibits a similar selectivity towards methane and acetaldehyde. This is consistent with the experimental selectivity observed for larger, predominantly (111)-exposing Rh nanoparticles and resolves the discrepancy to earlier first-principles MFA microkinetic work that found the Rh(111) facet to be selective towards acetaldehyde. While the latter work tried to approximately account for lateral interactions through coverage-dependent rate expressions, our analysis demonstrates that this fails to sufficiently capture concomitant correlations among the adsorbed reaction intermediates that crucially determine the overall selectivity.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Interface between graphene and liquid Cu from molecular dynamics simulations

Controllable synthesis of defect-free graphene is crucial for applications since the properties of graphene are highly sensitive to any deviations from the crystalline lattice. We focus here on the emerging use of liquid Cu catalysts, which has high potential for fast and efficient industrial-scale production of high-quality graphene. The interface between graphene and liquid Cu is studied using force field and ab initio molecular dynamics, revealing a complete or partial embedding of finite-sized flakes. By analyzing flakes of different sizes we find that the size-dependence of the embedding can be rationalized based on the energy cost of embedding versus bending the graphene flake. The embedding itself is driven by the formation of covalent bonds between the under-coordinated edge C atoms and the liquid Cu surface, which is accompanied by a significant charge transfer. In contrast, the central flake atoms are located around or slightly above 3 Å from the liquid Cu surface and exhibit weak vdW-bonding and much lower charge transfer. The structural and electronic properties of the embedded state revealed in our work provides the atomic-scale information needed to develop effective models to explain the special growth observed in experiments where various interesting phenomena such as flake self-assembly and rotational alignment, high growth speeds and low defect densities in the final graphene product have been observed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Atomic structures and orbital energies of 61,489 crystal-forming organic molecules

Data science and machine learning in materials science require large datasets of technologically relevant molecules or materials. Currently, publicly available molecular datasets with realistic molecular geometries and spectral properties are rare. We here supply a diverse benchmark spectroscopy dataset of 61,489 molecules extracted from organic crystals in the Cambridge Structural Database (CSD), denoted OE62. Molecular equilibrium geometries are reported at the Perdew-Burke-Ernzerhof (PBE) level of density functional theory (DFT) including van der Waals corrections for all 62k molecules. For these geometries, OE62 supplies total energies and orbital eigenvalues at the PBE and the PBE hybrid (PBE0) functional level of DFT for all 62k molecules in vacuum as well as at the PBE0 level for a subset of 30,876 molecules in (implicit) water. For 5,239 molecules in vacuum, the dataset provides quasiparticle energies computed with many-body perturbation theory in the $G_0W_0$ approximation with a PBE0 starting point (denoted GW5000 in analogy to the GW100 benchmark set (M. van Setten et al. J. Chem. Theory Comput. 12, 5076 (2016))).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adsorption structures and energetics of molecules on metal surfaces: Bridging experiment and theory

Adsorption geometry and stability of organic molecules on surfaces are key parameters that determine the observable properties and functions of hybrid inorganic/organic systems (HIOSs). Despite many recent advances in precise experimental characterization and improvements in first-principles electronic structure methods, reliable databases of structures and energetics for large adsorbed molecules are largely amiss. In this review, we present such a database for a range of molecules adsorbed on metal single-crystal surfaces. The systems we analyze include noble-gas atoms, conjugated aromatic molecules, carbon nanostructures, and heteroaromatic compounds adsorbed on five different metal surfaces. The overall objective is to establish a diverse benchmark dataset that enables an assessment of current and future electronic structure methods, and motivates further experimental studies that provide ever more reliable data. Specifically, the benchmark structures and energetics from experiment are here compared with the recently developed van der Waals (vdW) inclusive density-functional theory (DFT) method, DFT+vdW$^{\mathrm{surf}}$. In comparison to 23 adsorption heights and 17 adsorption energies from experiment we find a mean average deviation of 0.06 Å and 0.16 eV, respectively. This confirms the DFT+vdW$^{\mathrm{surf}}$ method as an accurate and efficient approach to treat HIOSs. A detailed discussion identifies remaining challenges to be addressed in future development of electronic structure methods, for which the here presented benchmark database may serve as an important reference.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Communication: Charge-Population Based Dispersion Interactions for Molecules and Materials

We introduce a system-independent method to derive effective atomic C$_6$ coefficients and polarizabilities in molecules and materials purely from charge population analysis. This enables the use of dispersion-correction schemes in electronic structure calculations without recourse to electron-density partitioning schemes and expands their applicability to semi-empirical methods and tight-binding Hamiltonians. We show that the accuracy of our method is en par with established electron-density partitioning based approaches in describing intermolecular C$_6$ coefficients as well as dispersion energies of weakly bound molecular dimers, organic crystals, and supramolecular complexes. We showcase the utility of our approach by incorporation of the recently developed many-body dispersion (MBD) method [Tkatchenko et al., Phys. Rev. Lett. 108, 236402 (2012)] into the semi-empirical Density Functional Tight-Binding (DFTB) method and propose the latter as a viable technique to study hybrid organic-inorganic interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy Dissipation during Diffusion at Metal Surfaces: Disentangling the Role of Phonons vs. Electron-Hole Pairs

Helium spin echo experiments combined with ab initio-based Langevin molecular dynamics simulations are used to quantify the adsorbate-substrate coupling during the thermal diffusion of Na atoms on Cu(111). An analysis of trajectories within the local density friction approximation allows the contribution from electron-hole pair excitations to be separated from the total energy dissipation. Despite the minimal electronic friction coefficient of Na and the relatively small mass mismatch to Cu promoting efficient phononic dissipation, about $(20\pm5)\%$ of the total energy loss is attributable to electronic friction. The results suggest a significant role of electronic non-adiabaticity in the rapid thermalization generally relied upon in adiabatic diffusion theories.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Function-Space Based Solution Scheme for the Size-Modified Poisson-Boltzmann Equation in Full-Potential DFT

The size-modified Poisson-Boltzmann (MPB) equation is an efficient implicit solvation model which also captures electrolytic solvent effects. It combines an account of the dielectric solvent response with a mean-field description of solvated finite-sized ions. We present a general solution scheme for the MPB equation based on a fast function-space oriented Newton method and a Green's function preconditioned iterative linear solver. In contrast to popular multi-grid solvers this approach allows to fully exploit specialized integration grids and optimized integration schemes. We describe a corresponding numerically efficient implementation for the full-potential density-functional theory (DFT) code FHI-aims. We show that together with an additional Stern layer correction the DFT+MPB approach can describe the mean activity coefficient of a KCl aqueous solution over a wide range of concentrations. The high sensitivity of the calculated activity coefficient on the employed ionic parameters thereby suggests to use extensively tabulated experimental activity coefficients of salt solutions for a systematic parametrization protocol.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Global Structure Search for Molecules on Surfaces: Efficient Sampling with Curvilinear Coordinates

Efficient structure search is a major challenge in computational materials science. We present a modification of the basin hopping global geometry optimization approach that uses a curvilinear coordinate system to describe global trial moves. This approach has recently been shown to be efficient in structure determination of clusters [Nano Letters 15, 8044-8048 (2015)] and is here extended for its application to covalent, complex molecules and large adsorbates on surfaces. The employed, automatically constructed delocalized internal coordinates are similar to molecular vibrations, which enhances the generation of chemically meaningful trial structures. By introducing flexible constraints and local translation and rotation of independent geometrical subunits we enable the use of this method for molecules adsorbed on surfaces and interfaces. For two test systems, trans-$β$-ionylideneacetic acid adsorbed on a Au(111) surface and methane adsorbed on a Ag(111) surface, we obtain superior performance of the method compared to standard optimization moves based on Cartesian coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal and electronic fluctuations of flexible adsorbed molecules: Azobenzene on Ag(111)

We investigate the thermal and electronic collective fluctuations that contribute to the finite-temperature adsorption properties of flexible adsorbates on surfaces on the example of the molecular switch azobenzene C$_{12}$H$_{10}$N$_{2}$ on the Ag(111) surface. Using first-principles molecular dynamics simulations we obtain the free energy of adsorption that accurately accounts for entropic contributions, whereas the inclusion of many-body dispersion interactions accounts for the electronic correlations that govern the adsorbate binding. We find the adsorbate properties to be strongly entropy-driven, as can be judged by a kinetic molecular desorption prefactor of 10$^{24}$ s$^{-1}$ that largely exceeds previously reported estimates. We relate this effect to sizable fluctuations across structural and electronic observables. Comparison of our calculations to temperature-programmed desorption measurements demonstrates that finite-temperature effects play a dominant role for flexible molecules in contact with polarizable surfaces, and that recently developed first-principles methods offer an optimal tool to reveal novel collective behavior in such complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Critical analysis of fragment-orbital DFT schemes for the calculation of electronic coupling values

We present a critical analysis of the popular fragment-orbital density-functional theory (FO-DFT) scheme for the calculation of electronic coupling values. We discuss the characteristics of different possible formulations or 'flavors' of the scheme which differ by the number of electrons in the calculation of the fragments and the construction of the Hamiltonian. In addition to two previously described variants based on neutral fragments, we present a third version taking a different route to the approximate diabatic state by explicitly considering charged fragments. In applying these FO-DFT flavors to the two molecular test sets HAB7 (electron transfer) and HAB11 (hole transfer) we find that our new scheme gives improved electronic couplings for HAB7 (-6.2% decrease in mean relative signed error) and greatly improved electronic couplings for HAB11 (-15.3% decrease in mean relative signed error). A systematic investigation of the influence of exact exchange on the electronic coupling values shows that the use of hybrid functionals in FO-DFT calculations improves the electronic couplings, giving values close to or even better than more sophisticated constrained DFT calculations. Comparing the accuracy and computational cost of each variant we devise simple rules to choose the best possible flavor depending on the task. For accuracy, our new scheme with charged-fragment calculations performs best, while numerically more efficient at reasonable accuracy is the variant with neutral fragments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electronic friction-based vibrational lifetimes of molecular adsorbates: Beyond the independent atom approximation

We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intra-molecular contributions and overcomes the systematic underestimation of the non-adiabatic losses obtained within the prevalent independent atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fingerprints of energy dissipation for exothermic surface chemical reactions: O$_2$ on Pd(100)

We present first-principles calculations of the sticking coefficient of O$_2$ at Pd(100) to assess the effect of phononic energy dissipation on this kinetic parameter. For this we augment dynamical simulations on six-dimensional potential energy surfaces (PESs) representing the molecular degrees of freedom with various effective accounts of surface mobility. In comparison to the prevalent frozen-surface approach energy dissipation is found to qualitatively affect the calculated sticking curves. At the level of a generalized Langevin oscillator (GLO) model we achieve good agreement with experimental data. The agreement is similarly reached for PESs based on two different semi-local density-functional theory functionals. This robustness of the simulated sticking curve does not extend to the underlying adsorption mechanism, which is predominantly directly dissociative for one functional or molecularly trapped for the other. Completely different adsorption mechanisms therewith lead to rather similar sticking curves that agree equally well with the experimental data. This highlights the danger of the prevalent practice to extract corresponding mechanistic details from simple fingerprints of measured sticking data for such exothermic surface reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

First-principles embedded-cluster calculations of the neutral and charged oxygen vacancy at the rutile TiO$_2$(110) surface

We perform full-potential screened-hybrid density-functional theory (DFT) calculations to compare the thermodynamic stability of neutral and charged states of the surface oxygen vacancy at the rutile TiO$_2$(110) surface. Solid-state (QM/MM) embedded-cluster calculations are employed to account for the strong TiO$_2$ polarization response to the charged defect states. Similar to the situation for the bulk O vacancy, the +2 charge state $V_{\rm O}^{2+}$ is found to be energetically by far most stable. Only for Fermi-level positions very close to the conduction band, small polarons may at best be trapped by the charged vacancy. The large decrease of the $V_{\rm O}^{2+}$ formation energy with decreasing Fermi-level position indicates strongly enhanced surface O vacancy concentrations for $p$-doped samples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interfacial charge rearrangement and intermolecular interactions: Density-functional theory study of free-base porphine adsorbed on Ag(111) and Cu(111)

We employ dispersion-corrected density-functional theory to study the adsorption of tetrapyrrole 2H-porphine (2H-P) at Cu(111) and Ag(111). Various contributions to adsorbate-substrate and adsorbate-adsorbate interactions are systematically extracted to analyze the self-assembly behavior of this basic building block to porphyrin-based metal-organic nanostructures. This analysis reveals a surprising importance of substrate-mediated van der Waals interactions between 2H-P molecules, in contrast to negligible direct dispersive interactions. The resulting net repulsive interactions rationalize the experimentally observed tendency for single molecule adsorption.

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

Spin Manipulation by Creation of Single-Molecule Radical Cations

All-trans-retinoic acid (ReA), a closed-shell organic molecule comprising only C, H, and O atoms, is investigated on a Au(111) substrate using scanning tunneling microscopy and spectroscopy. In dense arrays single ReA molecules are switched to a number of states, three of which carry a localized spin as evidenced by conductance spectroscopy in high magnetic fields. The spin of a single molecule may be reversibly switched on and off without affecting its neighbors. We suggest that ReA on Au is readily converted to a radical by the abstraction of an electron.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structure Sensitivity in Oxide Catalysis: First-Principles Kinetic Monte Carlo Simulations for CO Oxidation at RuO$_2$(111)

We present a density-functional theory based kinetic Monte Carlo study of CO oxidation at the (111) facet of RuO$_2$. We compare the detailed insight into elementary processes, steady-state surface coverages and catalytic activity to equivalent published simulation data for the frequently studied RuO$_2$(110) facet. Qualitative differences are identified in virtually every aspect ranging from binding energetics over lateral interactions to the interplay of elementary processes at the different active sites. Nevertheless, particularly at technologically relevant elevated temperatures, near-ambient pressures and near-stoichiometric feeds both facets exhibit almost identical catalytic activity. These findings challenge the traditional definition of structure sensitivity based on macroscopically observable turnover frequencies and allow to scrutinize the applicability of structure sensitivity classifications developed for metals to oxide catalysis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Embedded-Cluster Calculations in a Numeric Atomic Orbital Density-Functional Theory Framework

We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level DFT calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO2(110).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling Heat Dissipation at the Nanoscale: An Embedding Approach for Chemical Reaction Dynamics on Metal Surfaces

We present an embedding technique for metallic systems that makes it possible to model energy dissipation into substrate phonons during surface chemical reactions from first principles. The separation of chemical and elastic contributions to the interaction potential provides a quantitative description of both electronic and phononic band structure. Application to the dissociation of O$_2$ at Pd(100) predicts translationally "hot" oxygen adsorbates as a consequence of the released adsorption energy (ca. 2.6 eV). This finding questions the instant thermalization of reaction enthalpies generally assumed in models of heterogeneous catalysis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temperature-dependent templated growth of porphine thin films on the (111) facets of copper and silver

The templated growth of the basic porphyrin unit, free-base porphine (2H-P), is characterized by means of X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine-structure (NEXAFS) spectroscopy measurements and density functional theory (DFT). The DFT simulations allow the deconvolution of the complex XPS and NEXAFS signatures into contributions originating from five inequivalent carbon atoms, which can be grouped into C-N and C-C bonded species. Polarization-dependent NEXAFS measurements reveal an intriguing organizational behavior: On both Cu(111) and Ag(111), for coverages up to one monolayer, the molecules adsorb undeformed and parallel to the respective metal surface. Upon increasing the coverage, however, the orientation of the molecules in the thin films depends on the growth conditions. Multilayers deposited at low temperatures (LT) exhibit a similar average tilting angle (30 degree relative to the surface plane) on both substrates. Conversely, for multilayers grown at room temperature a markedly different scenario exists. On Cu(111) the film thickness is self-limited to a coverage of approximately two layers, while on Ag(111) multilayers can be grown easily and, in contrast to the bulk 2H-P crystal, the molecules are oriented perpendicular to the surface. This difference in molecular orientation results in a modified line-shape of the C 1s XPS signatures that is dependent on the incident photon energy, which is explained by comparison with depth-resolved DFT calculations. Simulations of ionization energies for differently stacked molecules show no indication for a packing-induced modification of the multilayer XP spectra, thus indicating that the comparison of single molecule calculations to multilayer data is justified.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermodynamics of Surface Defects at the Aspirin/Water Interface

We present a simulation scheme to calculate defect formation free energies at a molecular crystal/water interface based on force-field molecular dynamics (MD) simulations. To this end we adopt and modify existing approaches to calculate binding free energies of biological ligand/receptor complexes to be applicable to common surface defects, such as step edges and kink sites. We obtain statistically accurate and reliable free energy values for the aspirin/water interface, which can be applied to estimate the distribution of defects using well-established thermodynamic relations. As a show case we calculate the free energy upon dissolving molecules from kink sites at the interface. This free energy can be related to the solubility concentration and we obtain solubility values in excellent agreement with experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

$M$Si$_{20}$H$_{20}$ Aggregates: From Simple Building Blocks to Highly Magnetic Functionalized Materials

Density-functional theory based global geometry optimization is used to scrutinize the possibility of using endohedrally-doped hydrogenated Si clusters as building blocks for constructing highly magnetic materials. In contrast to the known clathrate-type facet-sharing, the clusters exhibit a predisposition to aggregation through double Si-Si bridge bonds. For the prototypical CrSi$_{20}$H$_{20}$ cluster we show that reducing the degree of hydrogenation may be used to control the number of reactive sites to which other cages can be attached, while still preserving the structural integrity of the building block itself. This leads to a toolbox of CrSi$_{20}$H$_{20-2n}$ monomers with different number of double "docking sites", that allows building network architectures of any morphology. For (CrSi$_{20}$H$_{18}$)$_{2}$ dimer and [CrSi$_{20}$H$_{16}$](CrSi$_{20}$H$_{18}$)$_{2}$ trimer structures we illustrate that such aggregates conserve the high spin moments of the dopant atoms and are therefore most attractive candidates for cluster-assembled materials with unique magnetic properties. The study suggests that the structural completion of the individual endohedral cages within the doubly-bridge bonded structures and the high thermodynamic stability of the obtained aggregates are crucial for potential synthetic polymerization routes $via$ controlled dehydrogenation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

CO oxidation on Pd(100) vs. PdO(sqrt5xsqrt5)R27^o: First-Principles Kinetic Phase Diagrams and Bistability Conditions

We present first-principles kinetic Monte Carlo (1p-kMC) simulations addressing the CO oxidation reaction at Pd(100) for gas-phase conditions ranging from ultra-high vacuum (UHV) to ambient pressures and elevated temperatures. For the latter technologically relevant regime there is a long-standing debate regarding the nature of the active surface. The pristine metallic surface, an ultra-thin (sqrt5xsqrt5)R27^o PdO(101) surface oxide, and thicker oxide layers have each been suggested as the active state. We investigate these hypotheses with 1p-kMC simulations focusing on either the Pd(100) surface or the PdO(101) surface oxide and intriguingly obtain a range of (T, p)-conditions where both terminations appear metastable. The predicted bistability regime nicely ties in with oscillatory behavior reported experimentally by Hendriksen and coworkers [Catal. Today 105, 234 (2005)]. Within this regime we find that both surface terminations exhibit very similar intrinsic reactivity, which puts doubts on attempts to assign the catalytic function to just one active state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computational screening study towards redox-active metal-organic frameworks

The metal-organic framework (MOF) MFU-4l containing Co(II) centers and Cl- ligands has recently shown promising redox activity. Aiming for further improved MOF catalysts for oxidation processes employing molecular oxygen we present a density-functional theory (DFT) based computational screening approach to identify promising metal center and ligand combinations within the MFU-4l structural family. Using the O2 binding energy as a descriptor for the redox property, we show that relative energetic trends in this descriptor can reliably be obtained at the hybrid functional DFT level and using small cluster (scorpionate-type complex) models. Within this efficient computational protocol we screen a range of metal center / ligand combinations and identify several candidate systems that offer more exothermic O2 binding than the original Co/Cl-based MFU-4l framework.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Control of Intermolecular Bonds by Deposition Rates at Room Temperature: Hydrogen Bonds versus Metal Coordination in Trinitrile Monolayers

Self-assembled monolayers of 1,3,5-tris(4'-biphenyl-4"-carbonitrile)benzene, a large functional trinitrile molecule, on the (111) surfaces of copper and silver under ultrahigh vacuum conditions were studied by scanning tunneling microscopy and low-energy electron diffraction. A densely packed hydrogen-bonded polymorph was equally observed on both surfaces. Additionally, deposition onto Cu(111) yielded a well-ordered metal-coordinated porous polymorph that coexisted with the hydrogen-bonded structure. The required coordination centers were supplied by the adatom gas of the Cu(111) surface. On Ag(111), however, the well-ordered metal-coordinated network was not observed. Differences between the adatom reactivities on copper and silver and the resulting bond strengths of the respective coordination bonds are held responsible for this substrate dependence. By utilizing ultralow deposition rates, we demonstrate that on Cu(111) the adatom kinetics plays a decisive role in the expression of intermolecular bonds and hence structure selection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Excited-state potential-energy surfaces of metal-adsorbed organic molecules from Linear Expansion Δ-Self-Consistent Field Density-Functional Theory (ΔSCF-DFT)

Accurate and efficient simulation of excited state properties is an important and much aspired cornerstone in the study of adsorbate dynamics on metal surfaces. To this end, the recently proposed linear expansion ΔSelf-Consistent Field (leΔSCF) method by Gavnholt et al. [Phys. Rev. B 78, 075441 (2008)] presents an efficient alternative to time consuming quasi-particle calculations. In this method the standard Kohn-Sham equations of Density-Functional Theory are solved with the constraint of a non-equilibrium occupation in a region of Hilbert-space resembling gas-phase orbitals of the adsorbate. In this work we discuss the applicability of this method for the excited-state dynamics of metal-surface mounted organic adsorbates, specifically in the context of molecular switching. We present necessary advancements to allow for a consistent quality description of excited-state potential-energy surfaces (PESs), and illustrate the concept with the application to Azobenzene adsorbed on Ag(111) and Au(111) surfaces. We find that the explicit inclusion of substrate electronic states modifies the topologies of intra-molecular excited-state PESs of the molecule due to image charge and hybridization effects. While the molecule in gas phase shows a clear energetic separation of resonances that induce isomerization and backreaction, the surface-adsorbed molecule does not. The concomitant possibly simultaneous induction of both processes would lead to a significantly reduced switching efficiency of such a mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exploring Morphology-Activity Relationships: Ab Initio Wulff Construction for RuO2 Nanoparticles under Oxidizing Conditions

We present a density-functional theory based Wulff construction of the equilibrium shape of RuO2 particles in an oxygen environment. The obtained intricate variations of the crystal habit with the oxygen chemical potential allow for a detailed discussion of the dependence on the oxidizing pretreatment observed in recent powder catalyst studies. The analysis specifically indicates an incomplete particle shape equilibration in previously employed low temperature calcination. Equilibrated particles could be active CO oxidation catalysts with long-term stability in oxidizing feed and then represent an interesting alternative to the previously suggested core-shell concept.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First-principles thermodynamic screening approach to photo-catalytic water splitting with co-catalysts

We adapt the computational hydrogen electrode approach to explicitly account for photo-generated charges and use it to computationally screen for viable catalyst/co-catalyst combinations for photo-catalytic water splitting. The hole energy necessary to thermodynamically drive the reaction is employed as descriptor for the screening process. Using this protocol and hybrid-level density-functional theory we show that water oxidation on bare TiO2 surfaces is thermodynamically more complex than previously thought. This motivates a screening for suitable co-catalysts for this half-reaction, which we carry out for Au particles down to the non-scalable size regime. We find that almost all small Au clusters studied are better suited for water photo-oxidation than an extended Au(111) surface or bare TiO2 facets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multi-Doping of Si Cages: High Spin States beyond the Single-Dopant Septet Limit

Density-functional theory based global geometry optimization is employed to systematically scrutinize the possibility of multi-doping of hydrogenated Si clusters in order to achieve high spin states beyond the septet limit of a single-atom dopant. While our unbiased configurational search reveals that the previously suggested Si18H12 double hexagonal prism structure is generally too small to accommodate two dopants in magnetized state, the larger Si24H24 cage turns out to be suitable for such applications. For dimer dopants M2+ = Cr2+, Mn2+ and CrMn+, the structural integrity of the host cage is conserved in the ground-state structure of corresponding M2+@Si24H24 aggregates, as is the unusually high spin state of the guest dopant, which in case of Cr2+ already exceeds the single-atom dopant septet limit by almost a factor of two. Moreover, the possibility of further increasing the cluster spin moment by encapsulating an even larger number of dopants into a suitably sized hydrogenated Si cage is illustrated for the example of a (CrMn+)2@Si28H28 aggregate with a total number of 18 unpaired electrons. These results strongly suggest multi-doping of Si clusters as a viable route to novel cluster-based materials for magneto-optic applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ready, set and no action: A static perspective on potential energy surfaces commonly used in gas-surface dynamics

In honoring the seminal contribution of Henry Eyring and Michael Polanyi who first introduced the concept of potential energy surfaces (PESs) to describe chemical reactions in gas-phase [Z. Phys. Chem. 12, 279-311, (1931)], this work comes to review and assess state-of-the-art approaches towards first-principle based modeling in the field of gas-surface dynamics. Within the Born-Oppenheimer and frozen surface approximations, the O2-Ag(100) interaction energetics are used as a showcase system to accentuate the complex landscape exhibited by the PESs employed to describe the impingement of diatomics on metal substrates and draw attention to the far-from-trivial task of continuously representing them within all six molecular degrees of freedom. To this end, the same set of ab initio reference data obtained within Density Functional Theory (DFT) are continuously represented by two different state-of-the-art high-dimensional approaches, namely the Corrugation-Reducing Procedure and Neural Networks. Exploiting the numerically undemanding nature of the resulting representations, a detailed static evaluation is performed on both PESs based on an extensive global minima search. The latter proved particularly illuminating in revealing representation deficiencies which affect the dynamical picture yet go otherwise unnoticed within the so-called "divide-and-conquer" approach.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assessing computationally efficient isomerization dynamics: Delta-SCF density-functional theory study of azobenzene molecular switching

We present a detailed comparison of the S0, S1 (n -> π*) and S2 (π-> π*) potential energy surfaces (PESs) of the prototypical molecular switch azobenzene as obtained by Delta-self-consistent-field (Delta-SCF) Density-Functional Theory (DFT), time-dependent DFT (TD-DFT) and approximate Coupled Cluster Singles and Doubles (RI-CC2). All three methods unanimously agree in terms of the PES topologies, which are furthermore fully consistent with existing experimental data concerning the photo-isomerization mechanism. In particular, sum-method corrected Delta-SCF and TD-DFT yield very similar results for S1 and S2, when based on the same ground-state exchange-correlation (xc) functional. While these techniques yield the correct PES topology already on the level of semi-local xc functionals, reliable absolute excitation energies as compared to RI-CC2 or experiment require an xc treatment on the level of long-range corrected hybrids. Nevertheless, particularly the robustness of Delta-SCF with respect to state crossings as well as its numerical efficiency suggest this approach as a promising route to dynamical studies of larger azobenzene-containing systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evaluation of endohedral doping of hydrogenated Si fullerenes as a route to magnetic Si building blocks

Density-functional theory based global geometry optimization is used to scrutinize the possibility of endohedral doping of hydrogenated Si fullerenes as a route to Si nanostructures with high magnetic moments. In contrast to previous suggestions, our unbiased sampling finds the smallest Si16H16 endohedral cage generally too small to encapsulate 3d metal dopant atoms. For the next larger fullerene-like cage though, we identify perfectly symmetric MSi20H20 (M = Co, Ti, V, Cr) cage structures as ground states. These structures conserve the high spin moment of the dopant atom and therewith underscore the potential of this Si nanoform for novel cluster-based materials with unique magnetic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Molecular-level understanding of the WGS and reverse WGS reactions on Rh through hierarchical multiscale approach

Hierarchically combining semi-empirical methods and first-principles calculations we gain a novel and noteworthy picture of the molecular-level mechanisms that govern the water-gas-shift (WGS) and reverse water-gas-shift (r-WGS) reactions on Rh catalysts. Central to this picture is that the WGS and r-WGS follow two different dominant reaction mechanisms: WGS proceeds according to a carboxyl (COOH) mechanism, whereas r-WGS proceeds according to a redox (CO2 {\to} CO + O) mechanism. The obtained results furthermore underscore the danger of common first-principles analyses that focus on a priori selected dominant paths. Not restricted to such bias, our herein proposed hierarchical approach thus constitutes a promising avenue to properly transport and incorporate the ab initio predictive-quality to a new level of system complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-adiabatic effects during the dissociative adsorption of O2 at Ag(111)? A first-principles divide and conquer study

We study the gas-surface dynamics of O2 at Ag(111) with the particular objective to unravel whether electronic non-adiabatic effects are contributing to the experimentally established inertness of the surface with respect to oxygen uptake. We employ a first-principles divide and conquer approach based on an extensive density-functional theory mapping of the adiabatic potential energy surface (PES) along the six O2 molecular degrees of freedom. Neural networks are subsequently used to interpolate this grid data to a continuous representation. The low computational cost with which forces are available from this PES representation allows then for a sufficiently large number of molecular dynamics trajectories to quantitatively determine the very low initial dissociative sticking coefficient at this surface. Already these adiabatic calculations yield dissociation probabilities close to the scattered experimental data. Our analysis shows that this low reactivity is governed by large energy barriers in excess of 1.1 eV very close to the surface. Unfortunately, these adiabatic PES characteristics render the dissociative sticking a rather insensitive quantity with respect to a potential spin or charge non-adiabaticity in the O2-Ag(111) interaction. We correspondingly attribute the remaining deviations between the computed and measured dissociation probabilities primarily to unresolved experimental issues with respect to surface imperfections.

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

When atomic-scale resolution is not enough: Spatial effects in in situ model catalyst studies

We investigate transport effects in in situ studies of defined model catalysts using a multi-scale modeling approach integrating first-principles kinetic Monte Carlo simulations into a fluid dynamical treatment. We specifically address two isothermal flow setups: i) a channel flow with the gas-stream approaching the single crystal from the side, as is representative for reactor scanning tunneling microscopy experiments; and ii) a stagnation flow with perpendicular impingement. Using the CO oxidation at RuO2 (110) as showcase we obtain substantial variations in the gas-phase pressures between the inlet and the catalyst surface. In the channel geometry the mass transfer limitations lead furthermore to pronounced lateral changes in surface composition across the catalyst surface. This prevents the aspired direct relation between activity and catalyst structure. For the stagnation flow the lateral variations are restricted to the edges of the catalyst. This allows to access the desired structure-activity relation using a simple model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electron-hole pairs during the adsorption dynamics of O2 on Pd(100) - Exciting or not?

During the exothermic adsorption of molecules at solid surfaces dissipation of the released energy occurs via the excitation of electronic and phononic degrees of freedom. For metallic substrates the role of the nonadiabatic electronic excitation channel has been controversially discussed, as the absence of a band gap could favour an easy coupling to a manifold of electronhole pairs of arbitrarily low energies. We analyse this situation for the highly exothermic showcase system of molecular oxygen dissociating at Pd(100), using time-dependent perturbation theory applied to first-principles electronic-structure calculations. For a range of different trajectories of impinging O2 molecules we compute largely varying electron-hole pair spectra, which underlines the necessity to consider the high-dimensionality of the surface dynamical process when assessing the total energy loss into this dissipation channel. Despite the high Pd density of states at the Fermi level, the concomitant non-adiabatic energy losses nevertheless never exceed about 5% of the available chemisorption energy. While this supports an electronically adiabatic description of the predominant heat dissipation into the phononic system, we critically discuss the non-adiabatic excitations in the context of the O2 spin transition during the dissociation process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the stability of "non-magic" endohedrally doped Si clusters: A first-principles sampling study of MSi16^+ (M =Ti,V,Cr)

Density-functional theory is used to study the geometric and electronic structure of cationic Si16^+ clusters with a Ti, V or Cr dopant atom. Through unbiased global geometry optimization based on the basin-hopping approach we confirm that a Frank- Kasper polyhedron with the metal atom at the center represents the ground-state isomer for all three systems. The endohedral cage geometry is thus stabilized even though only VSi16^+ achieves electronic shell closure within the prevalent spherical potential model. Our analysis of the electronic structure traces this diminished role of shell closure for the stabilization back to the adaptive capability of the metal- Si bonding, which is more the result of a complex hybridization than the orginally proposed mere formal charge transfer. The resulting flexibility of the metal-Si bond can help to stabilize also "non-magic" cage-dopant combinations, which suggests that a wider range of materials may eventually be cast into this useful geometry for cluster-assembled materials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Azobenzene versus 3,3',5,5'-tetra-tert-butyl-azobenzene (TBA) at Au(111): Characterizing the role of spacer groups

We present large-scale density-functional theory (DFT) calculations and temperature programmed desorption measurements to characterize the structural, energetic and vibrational properties of the functionalized molecular switch 3,3',5,5'-tetra-tert-butyl-azobenzene (TBA) adsorbed at Au(111). Particular emphasis is placed on exploring the accuracy of the semi-empirical dispersion correction approach to semi-local DFT (DFT-D) in accounting for the substantial van der Waals component in the surface chemical bond. In line with previous findings for benzene and pure azobenzene at coinage metal surfaces, DFT-D significantly overbinds the molecule, but seems to yield an accurate adsorption geometry as far as can be judged from the experimental data. Comparing the trans adsorption geometry of TBA and azobenzene at Au(111) reveals a remarkable insensitivity of the structural and vibrational properties of the -N=N- moiety. This questions the established view of the role of the bulky tert-butyl-spacer groups for the switching of TBA in terms of a mere geometric decoupling of the photochemically active diazo-bridge from the gold substrate.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structural Metastability of Endohedral Silicon Fullerenes

Endohedrally doped Si20 fullerenes appear as appealing building blocks for nanoscale materials. We investigate their structural stability with an unbiased and systematic global geometry optimization method within density-functional theory. For a wide range of metal doping atoms, it was sufficient to explore the Born Oppenheimer surface for only a moderate number of local minima to find structures that clearly differ from the initial endohedral cages, but are considerably more favorable in terms of energy. Previously proposed structures are thus all metastable.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Water adsorption and dissociation on SrTiO3(001) revisited: A density-functional theory study

We present a comprehensive density-functional theory study addressing the adsorption, dissociation and successive diffusion of water molecules on the two regular terminations of SrTiO3(001). Combining the obtained supercell-geometry converged energetics within a first-principles thermodynamics framework we are able to reproduce the experimentally observed hydroxilation of the SrO-termination already at lowest background humidity, whereas the TiO2-termination stays free of water molecules in the regime of low water partial pressures. This different behavior is traced back to the effortless formation of energetically very favorable hydroxyl-pairs on the prior termination. Contrary to the prevalent understanding our calculations indicate that at low coverages also the less water-affine TiO2-termination can readily decompose water, with the often described molecular state only stabilized towards higher coverages.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Azobenzene at Coinage Metal Surfaces: The Role of Dispersive van der Waals Interactions

We use different semi-empirical dispersion correction schemes to assess the role of long-range van der Waals interactions in the adsorption of the prototypical molecular switch azobenzene (C6H5-N2-C6H5) at the coinage metal surfaces Cu(111), Ag(111) and Au(111). Compared to preceding density-functional theory results employing a semi-local exchange and correlation functional we obtain partly sizable changes of the computed adsorption geometry and energetics. The discomforting scatter in the results provided by the different schemes is largely attributed to the unknown form of the damping function in the semi-empirical correction expression. Using the congeneric problem of the adsorption of benzene as a vehicle to connection with experiment, we cautiously conclude that the account of dispersive interactions at the metal surfaces provided by the various schemes is in the right ballpark, with the more recent, general schemes likely to overbind.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Robustness of "cut and splice" genetic algorithms in the structural optimization of atomic clusters

We return to the geometry optimization problem of Lennard-Jones clusters to analyze the performance dependence of "cut and splice" genetic algorithms (GAs) on the employed population size. We generally find that admixing twinning mutation moves leads to an improved robustness of the algorithm efficiency with respect to this a priori unknown technical parameter. The resulting very stable performance of the corresponding mutation+mating GA implementation over a wide range of population sizes is an important feature when addressing unknown systems with computationally involved first-principles based GA sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stabilizing a Molecular Switch at Solid Surfaces: A Density-Functional Theory Study of Azobenzene at Cu(111), Ag(111), and Au(111)

We present a density-functional theory trend study addressing the binding of the trans-cis conformational switch azobenzene (C6H5-N=N-C6H5) at three coinage metal surfaces. From the reported detailed energetic, geometric, and electronic structure data we conclude that the governing factor for the molecule-surface interaction is a competition between covalent bonding of the central azo (-N=N-) bridge on the one hand and the surface interaction of the two closed-shell phenyl (-C6H5) rings on the other. With respect to this factor the cis conformer exhibits a more favorable gas-phase geometric structure and is thus more stabilized at the studied surfaces. With the overall binding still rather weak the relative stability of the two isomers is thereby reduced at Ag(111) and Au(111). This is significantly different at Cu(111), where the cis bonding is strong enough to even reverse the gas-phase energetic order at the level of the employed semi-local electronic exchange and correlation (xc) functional. While this actual reversal may well be affected by the deficiencies due to the approximate xc treatment, we critically discuss that the rationalization of the general effect of the surface on the meta-stable molecular states is quite robust. This should equally hold for the presented analysis of recent tip-manipulation and photo-excitation isomerization experiments from the view point of the derived bonding mechanism.

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