The structure, energetics, and nature of the chemical bonding of phenylthiol adsorbed on the Au(111) surface: implications for density-functional calculations of molecular-electronic conduction

The adsorption of phenylthiol on the Au(111) surface is modeled using Perdew and Wang density-functional calculations. Both direct molecular physisorption and dissociative chemisorption via S-H bond cleavage are considered as well as dimerization to form disulfides. For the major observed product, the chemisorbed thiol, an extensive potential-energy surface is produced as a function of both the azimuthal orientation of the adsorbate and the linear translation of the adsorbate through the key fcc, hcp, bridge, and top binding sites. Key structures are characterized, the lowest-energy one being a broad minimum of tilted orientation ranging from the bridge structure halfway towards the fcc one. The vertically oriented threefold binding sites, often assumed to dominate molecular electronics measurements, are identified as transition states at low coverage but become favored in dense monolayers. A similar surface is also produced for chemisorption of phenylthiol on Ag(111); this displays significant qualitative differences, consistent with the qualitatively different observed structures for thiol chemisorption on Ag and Au. Full contours of the minimum potential energy as a function of sulfur translation over the crystal face are described, from which the barrier to diffusion is deduced to be 5.8 kcal mol(-1), indicating that the potential-energy surface has low corrugation. The calculated bond lengths, adsorbate charge and spin density, and the density of electronic states all indicate that, at all sulfur locations, the adsorbate can be regarded as a thiyl species that forms a net single covalent bond to the surface of strength 31 kcal mol(-1). No detectable thiolate character is predicted, however, contrary to experimental results for alkyl thiols that indicate up to 20%-30% thiolate involvement. This effect is attributed to the asymptotic-potential error of all modern density functionals that becomes manifest through a 3-4 eV error in the lineup of the adsorbate and substrate bands. Significant implications are described for density-functional calculations of through-molecule electron transport in molecular electronics.