The underpotential deposition of transition metal ions is a critical step in many electrosynthetic approaches. While underpotential deposition has been intensively studied at the atomic level, first-principles calculations in vacuum can strongly underestimate the stability of underpotentially deposited metals. It has been shown recently that the consideration of co-adsorbed anions can deliver more reliable descriptions of underpotential deposition reactions; however, the influence of additional key environmental factors such as the electrification of the interface under applied voltage and the activities of the ions in solution have yet to be investigated. In this work, copper underpotential deposition on gold is studied under realistic electrochemical conditions using a quantum-continuum model of the electrochemical interface. We report here on the influence of surface electrification, concentration effects, and anion co-adsorption on the stability of the copper underpotential deposition layer on the gold (100) surface.npj Computational Materials (2017) 3:1 ; doi:10.1038/s41524-016-0004-9 INTRODUCTION Underpotential deposition (UPD) has played an increasingly important role in the electrochemical preparation of nanomaterials with atomically thin metal film coatings for catalysis, imaging, and sensing applications.1-6 The UPD process is characterized by the formation of a (sub)monolayer of metal ions on a more noble metal substrate in a voltage range more positive than the reversible reduction potential of the adsorbing ion. The voltage at which the adlayer desorbs from the surface during an anodic scan is typically referenced to the bulk stripping potential of the adsorbed metal film and is termed the underpotential shift (ΔΦ upd ). Kolb and co-workers correlated underpotential shifts with differences in work functions of the substrate and the depositing metal (Fig. 1), suggesting that a charge transfer between the adlayer and the substrate may account for the larger adsorption energy of the adatom on the foreign surface.7 Since their seminal work, numerous studies have been performed in an effort to characterize a wide variety of UPD couples with the aim of understanding the voltammetric dependence of the formation and stability of the UPD layer in addition to its composition and surface structure. [8][9][10][11][12][13][14] Simultaneously, the theoretical aspects of UPD have been studied from first principles in an effort to connect calculated adsorption energies and model surface structures to experimentally measured underpotential shifts.10,15-23 For a thorough overview of the modeling of UPD phenomena, we direct the reader to the review by Sudha and Sangaranarayanan as well as a more recent review by Oviedo and co-workers.