Molecular modeling of alkyl monolayers on the Si(100)-2 x 1 surface

Molecular modeling was used to simulate various surfaces derived from the addition of 1-alkenes and 1-alkynes to Si=Si dimers on the Si(100)-2 x 1 surface. The primary aim was to better understand the interactions between adsorbates on the surface and distortions of the underlying silicon crystal due to functionalization. Random addition of ethylene and acetylene was used to determine how the addition of an adduct molecule affects subsequent additions for coverages up to one molecule per silicon dimer, that is, 100% coverage. Randomization subdues the effect that the relative positions of the adsorbates have on the enthalpy of the system. For ethylene and acetylene, the enthalpy of reaction changes less than 3 and 5 kcal/mol, respectively, from the first reacted species up to 100% coverage. As a result, a (near-)complete coverage is predicted, which is in line with experimental data. When 1-alkenes and 1-alkynes add by [2 + 2] addition, the hydrocarbon chains interact differently depending on the direction they project from the surface. These effects were investigated for four-carbon chains: 1-butene and 1-butyne. As expected, the chains that would otherwise intersect bend to avoid each other, raising the enthalpy of the system. For alkyl chains longer than four carbons, the chains are able to reorient themselves in a favorable manner, thus, resulting in a steady reduction in reaction enthalpy of about 2 kcal/mol for each additional methylene unit.