Dataset Information

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects.

ABSTRACT: Core-level energies are frequently calculated to explain the X-ray photoelectron spectra of metal-organic hybrid interfaces. The current paper describes how such simulations can be flawed when modeling interfaces between physisorbed organic molecules and metals. The problem occurs when applying periodic boundary conditions to correctly describe extended interfaces and simultaneously considering core hole excitations in the framework of a final-state approach to account for screening effects. Since the core hole is generated in every unit cell, an artificial dipole layer is formed. In this work, we study methane on an Al(100) surface as a deliberately chosen model system for hybrid interfaces to evaluate the impact of this computational artifact. We show that changing the supercell size leads to artificial shifts in the calculated core-level energies that can be well beyond 1 eV for small cells. The same applies to atoms at comparably large distances from the substrate, encountered, for example, in extended, upright-standing adsorbate molecules. We also argue that the calculated work function change due to a core-level excitation can serve as an indication for the occurrence of such an artifact and discuss possible remedies for the problem.

SUBMITTER: Taucher TC

PROVIDER: S-EPMC7557941 | biostudies-literature | 2020 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects.

Taucher Thomas C TC Hofmann Oliver T OT Zojer Egbert E

ACS omega 20200929 40

Core-level energies are frequently calculated to explain the X-ray photoelectron spectra of metal-organic hybrid interfaces. The current paper describes how such simulations can be flawed when modeling interfaces between physisorbed organic molecules and metals. The problem occurs when applying periodic boundary conditions to correctly describe extended interfaces and simultaneously considering core hole excitations in the framework of a final-state approach to account for screening effects. Sin ...[more]

PMID: 33073112

Similar Datasets

Project description:Accurate yet efficient computational models of solvent environment are central for most calculations that rely on atomistic modeling, such as prediction of protein-ligand binding affinities. In this study, we evaluate the accuracy of a recently developed generalized Born implicit solvent model, GBNSR6 (Aguilar et al. J. Chem. Theory Comput. 2010, 6, 3613-3639), in estimating the electrostatic solvation free energies (ΔG(pol)) and binding free energies (ΔΔG(pol)) for small protein-ligand complexes. We also compare estimates based on three different explicit solvent models (TIP3P, TIP4PEw, and OPC). The two main findings are as follows. First, the deviation (RMSD = 7.04 kcal/mol) of GBNSR6 binding affinities from commonly used TIP3P reference values is comparable to the deviations between explicit models themselves, e.g. TIP4PEw vs TIP3P (RMSD = 5.30 kcal/mol). A simple uniform adjustment of the atomic radii by a single scaling factor reduces the RMS deviation of GBNSR6 from TIP3P to within the above "error margin" - differences between ΔΔG(pol) estimated by different common explicit solvent models. The simple radii scaling virtually eliminates the systematic deviation (ΔΔG(pol)) between GBNSR6 and two out of the three explicit water models and significantly reduces the deviation from the third explicit model. Second, the differences between electrostatic binding energy estimates from different explicit models is disturbingly large; for example, the deviation between TIP4PEw and TIP3P estimates of ΔΔG(pol) values can be up to ∼50% or ∼9 kcal/mol, which is significantly larger than the "chemical accuracy" goal of ∼1 kcal/mol. The absolute ΔG(pol) calculated with different explicit models could differ by tens of kcal/mol. These discrepancies point to unacceptably high sensitivity of binding affinity estimates to the choice of common explicit water models. The absence of a clear "gold standard" among these models strengthens the case for the use of accurate implicit solvation models for binding energetics, which may be orders of magnitude faster.

Dataset Information

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects.

Publications

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets