Thiol-based molecular overlayers adsorbed on C60: Role of the end-group and charge state on the stability of the complexes

Abstract

We present pseudo-potential density functional theory calculations dedicated to analyze the stability and electronic properties of thiol-based molecular overlayers adsorbed on C60. We consider short molecules having a S atom as a headgroup, alkyl chains containing one to three C atoms, and a CH3 species as a terminal group. The thiol molecules are bonded to the carbon surface (through the S atom) with adsorption energies that vary in the range of ~1-2 eV and with S-C bond lengths of ~1.8 Å. For neutral C60(SCH3)n complexes, low energy atomic configurations are obtained when the thiol groups are distributed on the surface forming small molecular domains (e.g., pairs, trimers, or tetramer configurations of neighboring thiol molecules). In contrast, less stable random distributions are defined by orientationally disordered overlayers with highly distorted underlying carbon networks. The inclusion of London dispersion interaction slightly affects the structure of the molecular coating but increases the adsorption energies by values as large as 0.3 eV. Interestingly, the relative stability of the previous adsorbed phases differ from the one obtained when considering single sulfur adsorption on C60, a result that reveals the crucial role played by the terminal CH3 groups on the structure of the molecular coating. The positive (negative) charging of the [C60(SCH3)n](±q) complexes, with q as large as 8e, changes the geometrical structure and the chemical nature of the ligand shell inducing lateral molecular displacements, S-S bonding between neighboring thiols, as well as the partial degradation of the molecular coating. Finally, we consider the stability of two-component mixed overlayers formed by the coadsorption of CH3-, OH-, and NH2-terminated alkanethiols of the same length. In agreement with the results found on Au surfaces, we obtain lowest energy atomic configurations when molecular domains of a single component are stabilized on C60, a result that could be of fundamental importance in biomedical applications.