Search

Improving charge carrier mobilities in organic semiconductors is a challenging task that has hitherto primarily been tackled by empirical structural tuning of promising core compounds. Knowledge-based methods can greatly accelerate such local exploration, while a systematic analysis of large chemical databases can point towards promising design strategies. Here, we demonstrate such data mining by clustering an in-house database of >64.000 organic molecular crystals for which two charge-transport descriptors, the electronic coupling and the reorganization energy, have been calculated from first principles. The clustering is performed according to the Bemis-Murcko scaffolds of the constituting molecules and according to the sidegroups with which these molecular backbones are functionalized. In both cases, we obtain statistically significant structure-property relationships with certain scaffolds (sidegroups) consistently leading to favorable charge-transport properties. Functionalizing promising scaffolds with favorable sidegroups results in engineered molecular crystals for which we indeed compute improved charge-transport properties.

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the diffculty to describe the entire double layer forming at the electrified solid-liquid interface (SLI) within supercell sizes tractable by commonly employed density-functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread {\em ab initio} thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode's electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Computational screening for new and improved catalyst materials relies on accurate and low-cost predictions of key parameters such as adsorption energies. Here, we use recently developed compressed sensing methods to identify descriptors whose predictive power extends over a wide range of adsorbates, multi-metallic transition metal surfaces and facets. The descriptors are expressed as non-linear functions of intrinsic properties of the clean catalyst surface, e.g. coordination numbers, d-band moments and density of states at the Fermi level. From a single density-functional theory calculation of these properties, we predict adsorption energies at all potential surface sites, and thereby also the most stable geometry. Compared to previous approaches such as scaling relations, we find our approach to be both more general and more accurate for the prediction of adsorption energies on alloys with mixed-metal surfaces, already when based on training data including only pure metals. This accuracy can be systematically improved by adding also alloy adsorption energies to the training data., Comment: This document is the unedited Author's version of a Submitted Work that was subsequently accepted for publication in ACS Catalysis, copyright American Chemical Society after peer review. To access the final edited and published work see https://pubs.acs.org/doi/10.1021/acscatal.8b04478

Widely employed Near-Edge X-Ray Absorption Fine Structure (NEXAFS) spectroscopy probes a system by excitation of core electrons to unoccupied states. A variety of different methodologies are available to simulate corresponding spectra from first-principles. Core-level occupation constraints within ground-state Density-Functional Theory (DFT) represent a numerically most efficient means to this end that provides access to large systems, examples being surface adsorption, proteins, polymers, liquids, and buried, condensed phase interfaces (e.q. solid-liquid and solid-solid). Here, we systematically investigate the performance of different realizations of this approximate approach through the simulation of K-edge NEXAFS-spectra of a set of carbon and nitrogen-containing organic molecules. Variational collapse to the ground state and oscillatory convergence are the major complications of these approximate computational protocols. We present a modified version of the maximum-overlap method to achieve a self-consistent inclusion of electrons in virtual states for systems where convergence is hampered due to degeneracies. Our results demonstrate that reliable spectra allowing for a semi-quantitative analysis of experimental data are already obtained at the semi-local level of density functionals and with standard numeric atomic orbital basis sets., Comment: 12 pages, 10 figures, The following article has been accepted by The Journal of Chemical Physics. After it is published, it will be found at https://dx.doi.org/10.1063/1.5083618

We first briefly report on the status and recent hievements of the ELPA-AEO (Eigen-value Solvers for Petaflop Applications -- Algorithmic Extensions and Optimizations) and ESSEX II (Equipping Sparse Solvers for Exascale) projects. In both collaboratory efforts, scientists from the application areas, mathematicians, and computer scientists work together to develop and make available efficient highly parallel methods for the solution of eigenvalue problems. Then we focus on a topic addressed in both projects, the use of mixed precision computations to enhance efficiency. We give a more detailed description of our approaches for benefiting from either lower or higher precision in three selected contexts and of the results thus obtained.

Realistic oxide materials are often semiconductors, in particular at elevated temperatures, and their surfaces contain undercoordiated atoms at structural defects such as steps and corners. Using hybrid density-functional theory and ab initio atomistic thermodynamics, we investigate the interplay of bond-making, bond-breaking, and charge-carrier trapping at the corner defects at the (100) surface of a p-doped MgO in thermodynamic equilibrium with an O2 atmosphere. We show that by manipulating the coordination of surface atoms one can drastically change and even reverse the order of stability of reduced versus oxidized surface sites., Comment: 5 papges, 4 figures

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., Comment: 16 pages, 8 figures, 3 tables

We present high-pressure x-ray photoelectron spectroscopy (HP-XPS) and first-principles kinetic Monte Carlo study addressing the nature of the active surface in CO oxidation over Pd(100). Simultaneously measuring the chemical composition at the surface and in the near-surface gas phase, we reveal both O-covered pristine Pd(100) and a surface oxide as stable, highly active phases in the near-ambient regime accessible to HP-XPS. Surprisingly, no adsorbed CO can be detected during high CO2 production rates, which can be explained by a combination of a remarkably short residence time of the CO molecule on the surface and mass-transfer limitations in the present setup.

We present a density-functional theory study addressing the energetics and electronic structure properties of isolated oxygen adatoms at the SrTiO3(001) surface. Together with a surface lattice oxygen atom, the adsorbate is found to form a peroxide-type molecular species. This gives rise to a non-trivial topology of the potential energy surface for lateral adatom motion, with the most stable adsorption site not corresponding to the one expected from a continuation of the perovskite lattice. With computed modest diffusion barriers below 1 eV, it is rather the overall too weak binding at both regular SrTiO3(001) terminations that could be a critical factor for oxide film growth applications., Comment: 7 pages including 5 figures; related publications can be found at http://www.fhi-berlin.mpg.de/th/th.html

The oxidized Ag(111) surface has been studied by a combination of experimental and theoretical methods, scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), and density functional theory (DFT). A large variety of different surface structures is found, depending on the detailed preparation conditions. The observed structures fall into four classes: (a) individually chemisorbed atomic oxygen atoms, (b) three different oxygen overlayer structures, including the well-known p(4 x 4) phase, formed from the same Ag$_6$ and Ag$_{10}$ building blocks, (c) a c(4 x 8) structure not previously observed, and (d) at higher oxygen coverages structures characterized by stripes along the high-symmetry directions of the Ag(111) substrate. Our analysis provides a detailed explanation of the atomic-scale geometry of the Ag$_6$/Ag$_{10}$ building block structures, and the c(4 x 8) and stripe structures are discussed in detail. The observation of many different and co-existing structures implies that the O/Ag(111) system is characterized by a significantly larger degree of complexity than previously anticipated, and this will impact our understanding of oxidation catalysis processes on Ag catalysts., Comment: 12 pages, 9 figures; peer-reviewed and published version

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., Comment: 10 pages including 4 figures; related publications can be found at http://www.fhi-berlin.mpg.de/th/th.html

Recently, spin selection rules have been invoked to explain the discrepancy between measured and calculated adsorption probabilities of molecular oxygen reacting with Al(111). In this work, we inspect the impact of nonadiabatic spin transitions on the dynamics of this system from first principles. For this purpose the motion on two distinct potential-energy surfaces associated to different spin configurations and possible transitions between them are inspected by means of the Fewest Switches algorithm. Within this framework we especially focus on the influence of such spin transitions on observables accessible to molecular beam experiments. On this basis we suggest experimental setups that can validate the occurrence of such transitions and discuss their feasibility., Comment: 13 pages, 7 figures

The accurate description of molecule-surface interactions requires a detailed knowledge of the underlying potential-energy surface (PES). Recently, neural networks (NNs) have been shown to be an efficient technique to accurately interpolate the PES information provided for a set of molecular configurations, e.g. by first-principles calculations. Here, we further develop this approach by building the NN on a new type of symmetry functions, which allows to take the symmetry of the surface exactly into account. The accuracy and efficiency of such symmetry-adapted NNs is illustrated by the application to a six-dimensional PES describing the interaction of oxygen molecules with the Al(111) surface., Comment: 13 pages including 8 figures; related publications can be found at http://www.fhi-berlin.mpg.de/th/th.html

We employ a first-principles lattice-gas Hamiltonian (LGH) approach to determine the lateral interactions between O atoms adsorbed on the Pd(100) surface. With these interactions we obtain an ordering behavior at low coverage that is in quantitative agreement with experimental data. Uncertainties in the approach arise from the finite LGH expansion and from the approximate exchange-correlation (xc) functional underlying the employed density-functional theory energetics. We carefully scrutinize these uncertainties and conclude that they primarily affect the on-site energy, which rationalizes the agreement with the experimental critical temperatures for the order-disorder transition. We also investigate the validity of the frequently applied assumption that the ordering energies can be represented by a sum of pair terms. Restricting our LGH expansion to just pairwise lateral interactions, we find that this results in effective interactions which contain spurious contributions that are of equal size, if not larger than any of the uncertainties e.g. due to the approximate xc functional., Comment: 14 pages including 7 figures; related publications can be found at http://www.fhi-berlin.mpg.de/th/th.html

Using density-functional theory, we predict that the oxidation of the Ru(0001) surface proceeds via the accumulation of sub-surface oxygen in two-dimensional islands between the first and second substrate layer. This leads locally to a decoupling of an O-Ru-O trilayer from the underlying metal. Continued oxidation results in the formation and stacking of more of these trilayers, which unfold into the RuO_2(110) rutile structure once a critical film thickness is exceeded. Along this oxidation pathway, we identify various metastable configurations. These are found to be rather close in energy, indicating a likely lively dynamics between them at elevated temperatures, which will affect the surface chemical and mechanical properties of the material., Comment: 11 pages including 9 figures. Submitted to Phys. Rev. B. Related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.html

We present density-functional theory calculations of Ru 3d and O 1s surface core-level shifts (SCLSs) at an oxygen-rich Ru(0001) surface, namely for the O(1x1)/Ru(0001) chemisorption phase and for two surface terminations of fully oxidized RuO_2(110). Including final-state effects, the computed SCLSs can be employed for the analysis of experimental X-ray photoelectron spectroscopy (XPS) data enabling a detailed study of the oxidation behaviour of the Ru(0001) surface. We show that certain peaks can be used as a fingerprint for the existence of the various phases and propose that the long disputed satellite peak in RuO_2(110) XPS data originates from a hitherto unaccounted surface termination., Comment: 7 pages including 3 figures. Surf. Sci. in press. Related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.html

Using a Green's function approach, we investigate band structure effects in the BEEM current distribution in reciprocal space. In the elastic limit, this formalism provides a 'parameter free' solution to the BEEM problem. At low temperatures, and for thin metallic layers, the elastic approximation is enough to explain the experimental I(V) curves at low voltages. At higher voltages inelastic effects are approximately taken into account by introducing an effective RPA-electron lifetime, much in similarity with LEED theory. For thick films, however, additional damping mechanisms are required to obtain agreement with experiment., Comment: 4 pages, 3 postscript figures, revtex