Structure and energetics of liquid water-hydroxyl layers on Pt(111).

The interactions between liquid water and hydroxyl species on Pt(111) surfaces have been intensely investigated due to their importance to fuel cell electrocatalysis. Here we present a molecular dynamics study of their structure and energetics using an ensemble of neural network potentials, which allow us to obtain unprecedented statistical sampling. We first study the energetics of hydroxyl formation, where we find a near-linear adsorption energy profile, which exhibits a soft and gradual increase in the differential adsorption energy at high hydroxyl coverages. This is strikingly different from the predictions of the conventional bilayer model, which displays a kink at 1/3ML OH coverage indicating a sizeable jump in differential adsorption energy, but within the statistical uncertainty of previously reported ab initio molecular dynamics studies. We then analyze the structure of the interface, where we provide evidence for the water-OH/Pt(111) interface being hydrophobic at high hydroxyl coverages. We furthermore explain the observed adsorption energetics by analyzing the hydrogen bonding in the water-hydroxyl adlayers, where we argue that the increase in differential adsorption energy at high OH coverage can be explained by a reduction in the number of hydrogen bonds from the adsorbed water molecules to the hydroxyls.

Is the water/Pt(111) interface ordered at room temperature?

The structure of the water/Pt(111) interface has been a subject of debate over the past decades. Here, we report the results of a room temperature molecular dynamics study based on neural network potentials, which allow us to access long time scale simulations while retaining ab initio accuracy. We find that the water/Pt(111) interface is characterized by a double layer composed of a primary, strongly bound adsorption layer with a coverage of ∼0.15 ML, which is coupled to a secondary, weakly bound adsorption layer with a coverage of ∼0.58 ML. By studying the order of the primary adsorption layer, we find that there is an effective repulsion between the adsorbed water molecules, which gives rise to a dynamically changing, semi-ordered interfacial structure, where the water molecules in the primary adsorption layer are distributed homogeneously across the interface, forming frequent hydrogen bonds to water molecules in the secondary adsorption layer. We further show that these conclusions are beyond the time scales accessible to ab initio molecular dynamics.

Exploration versus Exploitation in Global Atomistic Structure Optimization.

The ability to navigate vast energy landscapes of molecules, clusters, and solids is a necessity for discovering novel compounds in computational chemistry and materials science. For high-dimensional systems, it is only computationally feasible to search a small portion of the landscape, and hence, the search strategy is of critical importance. Introducing Bayesian optimization concepts in an evolutionary algorithm framework, we quantify the concepts of exploration and exploitation in global minimum searches. The method allows us to control the balance between probing unknown regions of the landscape (exploration) and investigating further regions of the landscape known to have low-energy structures (exploitation). The search for global minima structures proves significantly faster with the optimal balance for three test systems (molecular compounds) and to a lesser extent also for a crystalline surface reconstruction. In addition, global search behaviors are analyzed to provide reasonable grounds for an optimal balance for different problems.

Dynamic breaking of a single gold bond.

While one might assume that the force to break a chemical bond gives a measure of the bond strength, this intuition is misleading. If the force is loaded slowly, thermal fluctuations may break the bond before it is maximally stretched, and the breaking force will be less than the bond can sustain. Conversely, if the force is loaded rapidly it is more likely that the maximum breaking force is measured. Paradoxically, no clear differences in breaking force were observed in experiments on gold nanowires, despite being conducted under very different conditions. Here we explore the breaking behaviour of a single Au-Au bond and show that the breaking force is dependent on the loading rate. We probe the temperature and structural dependencies of breaking and suggest that the paradox can be explained by fast breaking of atomic wires and slow breaking of point contacts giving very similar breaking forces.

The atomic simulation environment-a Python library for working with atoms.

The atomic simulation environment (ASE) is a software package written in the Python programming language with the aim of setting up, steering, and analyzing atomistic simulations. In ASE, tasks are fully scripted in Python. The powerful syntax of Python combined with the NumPy array library make it possible to perform very complex simulation tasks. For example, a sequence of calculations may be performed with the use of a simple 'for-loop' construction. Calculations of energy, forces, stresses and other quantities are performed through interfaces to many external electronic structure codes or force fields using a uniform interface. On top of this calculator interface, ASE provides modules for performing many standard simulation tasks such as structure optimization, molecular dynamics, handling of constraints and performing nudged elastic band calculations.

Band Gap Tuning and Defect Tolerance of Atomically Thin Two-Dimensional Organic-Inorganic Halide Perovskites.

Organic-inorganic halide perovskites have proven highly successful for photovoltaics but suffer from low stability, which deteriorates their performance over time. Recent experiments have demonstrated that low dimensional phases of the hybrid perovskites may exhibit improved stability. Here we report first-principles calculations for isolated monolayers of the organometallic halide perovskites (C4H9NH3)2MX2Y2, where M = Pb, Ge, Sn and X,Y = Cl, Br, I. The band gaps computed using the GLLB-SC functional are found to be in excellent agreement with experimental photoluminescence data for the already synthesized perovskites. Finally, we study the effect of different defects on the band structure. We find that the most common defects only introduce shallow or no states in the band gap, indicating that these atomically thin 2D perovskites are likely to be defect tolerant.

Defect-Tolerant Monolayer Transition Metal Dichalcogenides.

Localized electronic states formed inside the band gap of a semiconductor due to crystal defects can be detrimental to the material's optoelectronic properties. Semiconductors with a lower tendency to form defect induced deep gap states are termed defect-tolerant. Here we provide a systematic first-principles investigation of defect tolerance in 29 monolayer transition metal dichalcogenides (TMDs) of interest for nanoscale optoelectronics. We find that the TMDs based on group VI and X metals form deep gap states upon creation of a chalcogen (S, Se, Te) vacancy, while the TMDs based on group IV metals form only shallow defect levels and are thus predicted to be defect-tolerant. Interestingly, all the defect sensitive TMDs have valence and conduction bands with a very similar orbital composition. This indicates a bonding/antibonding nature of the gap, which in turn suggests that dangling bonds will fall inside the gap. These ideas are made quantitative by introducing a descriptor that measures the degree of similarity of the conduction and valence band manifolds. Finally, the study is generalized to nonpolar nanoribbons of the TMDs where we find that only the defect sensitive materials form edge states within the band gap.

Two-Dimensional Metal Dichalcogenides and Oxides for Hydrogen Evolution: A Computational Screening Approach.

We explore the possibilities of hydrogen evolution by basal planes of 2D metal dichalcogenides and oxides in the 2H and 1T class of structures using the hydrogen binding energy as a computational activity descriptor. For some groups of systems like the Ti, Zr, and Hf dichalcogenides the hydrogen bonding to the 2H structure is stronger than that to the 1T structure, while for the Cr, Mo, and W dichalcogenides the behavior is opposite. This is rationalized by investigating shifts in the chalcogenide p levels comparing the two structures. We find that usually for a given material only at most one of the two phases will be active for the hydrogen evolution reaction; however, in most cases the two phases are very close in formation energy, opening up the possibility for stabilizing the active phase. The study points to many new possible 2D HER materials beyond the few that are already known.

Catalysis. Assessing the reliability of calculated catalytic ammonia synthesis rates.

We introduce a general method for estimating the uncertainty in calculated materials properties based on density functional theory calculations. We illustrate the approach for a calculation of the catalytic rate of ammonia synthesis over a range of transition-metal catalysts. The correlation between errors in density functional theory calculations is shown to play an important role in reducing the predicted error on calculated rates. Uncertainties depend strongly on reaction conditions and catalyst material, and the relative rates between different catalysts are considerably better described than the absolute rates. We introduce an approach for incorporating uncertainty when searching for improved catalysts by evaluating the probability that a given catalyst is better than a known standard.

mBEEF: an accurate semi-local Bayesian error estimation density functional.

We present a general-purpose meta-generalized gradient approximation (MGGA) exchange-correlation functional generated within the Bayesian error estimation functional framework [J. Wellendorff, K. T. Lundgaard, A. Møgelhøj, V. Petzold, D. D. Landis, J. K. Nørskov, T. Bligaard, and K. W. Jacobsen, Phys. Rev. B 85, 235149 (2012)]. The functional is designed to give reasonably accurate density functional theory (DFT) predictions of a broad range of properties in materials physics and chemistry, while exhibiting a high degree of transferability. Particularly, it improves upon solid cohesive energies and lattice constants over the BEEF-vdW functional without compromising high performance on adsorption and reaction energies. We thus expect it to be particularly well-suited for studies in surface science and catalysis. An ensemble of functionals for error estimation in DFT is an intrinsic feature of exchange-correlation models designed this way, and we show how the Bayesian ensemble may provide a systematic analysis of the reliability of DFT based simulations.

Direct measurement and modulation of single-molecule coordinative bonding forces in a transition metal complex.

Coordination chemistry has been a consistently active branch of chemistry since Werner's seminal theory of coordination compounds inaugurated in 1893, with the central focus on transition metal complexes. However, control and measurement of metal-ligand interactions at the single-molecule level remain a daunting challenge. Here we demonstrate an interdisciplinary and systematic approach that enables measurement and modulation of the coordinative bonding forces in a transition metal complex. Terpyridine is derived with a thiol linker, facilitating covalent attachment of this ligand on both gold substrate surfaces and gold-coated atomic force microscopy tips. The coordination and bond breaking between terpyridine and osmium are followed in situ by electrochemically controlled atomic force microscopy at the single-molecule level. The redox state of the central metal atom is found to have a significant impact on the metal-ligand interactions. The present approach represents a major advancement in unravelling the nature of metal-ligand interactions and could have broad implications in coordination chemistry.

Investigation of Catalytic Finite-Size-Effects of Platinum Metal Clusters.

In this paper, we use density functional theory (DFT) calculations on highly parallel computing resources to study size-dependent changes in the chemical and electronic properties of platinum (Pt) for a number of fixed freestanding clusters ranging from 13 to 1415 atoms, or 0.7-3.5 nm in diameter. We find that the surface catalytic properties of the clusters converge to the single crystal limit for clusters with as few as 147 atoms (1.6 nm). Recently published results for gold (Au) clusters showed analogous convergence with size. However, this convergence happened at larger sizes, because the Au d-states do not contribute to the density of states around the Fermi-level, and the observed level fluctuations were not significantly damped until the cluster reached ca. 560 atoms (2.7 nm) in size.

Ab initio van der waals interactions in simulations of water alter structure from mainly tetrahedral to high-density-like.

The structure of liquid water at ambient conditions is studied in ab initio molecular dynamics simulations in the NVE ensemble using van der Waals (vdW) density-functional theory, i.e., using the new exchange-correlation functionals optPBE-vdW and vdW-DF2, where the latter has softer nonlocal correlation terms. Inclusion of the more isotropic vdW interactions counteracts highly directional hydrogen bonds, which are enhanced by standard functionals. This brings about a softening of the microscopic structure of water, as seen from the broadening of angular distribution functions and, in particular, from the much lower and broader first peak in the oxygen-oxygen pair-correlation function (PCF) and loss of structure in the outer hydration shells. Inclusion of vdW interactions is shown to shift the balance of resulting structures from open tetrahedral to more close-packed. The resulting O-O PCF shows some resemblance with experiment for high-density water (Soper, A. K. and Ricci, M. A. Phys. Rev. Lett. 2000, 84, 2881), but not directly with experiment for ambient water. Considering the accuracy of the new functionals for interaction energies, we investigate whether the simulation protocol could cause the deviation. An O-O PCF consisting of a linear combination of 70% from vdW-DF2 and 30% from low-density liquid water, as extrapolated from experiments, reproduces near-quantitatively the experimental O-O PCF for ambient water. This suggests the possibility that the new functionals may be reliable and that instead larger-scale simulations in the NPT ensemble, where the density is allowed to fluctuate in accordance with proposals for supercooled water, could resolve the apparent discrepancy with the measured PCF.

Nonlocal screening of plasmons in graphene by semiconducting and metallic substrates: first-principles calculations.

We investigate the role of substrates on the collective excitations of graphene by using a first-principles implementation of the density response function within the random-phase approximation. Specifically, we consider graphene adsorbed on SiC(0001) and Al(111) as representative examples of a semiconducting and metallic substrate. On SiC(0001), the long wavelength π plasmons are significantly damped although their energies remain almost unaltered. On Al(111), the long wavelength π plasmons are completely quenched due to the coupling to the metal surface plasmon. The strong damping of the plasmon excitations occurs despite the fact that the single-particle band structure of graphene is completely unaffected by the substrates illustrating the nonlocal nature of the effect.

Anomalous conductance oscillations and half-metallicity in atomic Ag-O chains.

Using spin density functional theory, we study the electronic and magnetic properties of atomically thin, suspended chains containing silver and oxygen atoms in an alternating sequence. Chains longer than 4 atoms develop a half-metallic ground state implying fully spin-polarized charge carriers. The conductances of the chains exhibit weak even-odd oscillations around an anomalously low value of 0.1G0 (G0=2e2/h) which coincide with the averaged experimental conductance in the long chain limit. The unusual conductance properties are explained in terms of a resonating-chain model, which takes the reflection probability and phase shift of a single bulk-chain interface as the only input. The model also explains the conductance oscillations for other metallic chains.

Avalanche size scaling in sheared three-dimensional amorphous solid.

We study the statistics of plastic rearrangement events in a simulated amorphous solid at T=0. Events are characterized by the energy release and the "slip volume", the product of plastic strain and system volume. Their distributions for a given system size L appear to be exponential, but a characteristic event size cannot be inferred, because the mean values of these quantities increase as Lalpha with alpha approximately 3/2. In contrast with results obtained in 2D models, we do not see simply connected avalanches. The exponent suggests a fractal shape of the avalanches, which is also evidenced by the mean fractal dimension and participation ratio.

Bayesian ensemble approach to error estimation of interatomic potentials.

Using a Bayesian approach a general method is developed to assess error bars on predictions made by models fitted to data. The error bars are estimated from fluctuations in ensembles of models sampling the model-parameter space with a probability density set by the minimum cost. The method is applied to the development of interatomic potentials for molybdenum using various potential forms and databases based on atomic forces. The calculated error bars on elastic constants, gamma-surface energies, structural energies, and dislocation properties are shown to provide realistic estimates of the actual errors for the potentials.