Thermal Stability of Organic Monolayers Grafted to Si(111): Insights from ReaxFF Reactive Molecular Dynamics Simulations

We used the ReaxFF reactive molecular dynamics simulations to investigate the chemical mechanisms and kinetics of thermal decomposition processes of silicon surfaces grafted with different organic molecules via Si-C bonds at atomistic level. In this work, we considered the Si(111) surface grafted with n-alkyl (ethyl, propyl, pentyl, and decyl) layers in 50% coverage and, Si-CH₃, Si-CCCH₃ and Si-CHCHCH₃ layers in full coverage. Si radicals primarily formed by the homolytic cleavage of Si-C bonds play a key role in the dehydrogenation processes that lead to the decomposition of the monolayers. Contrary to commonly proposed mechanisms that only involve a single Si atom center, we found that the main decomposition pathways require two Si lattice atoms to proceed. The ability of surface silyl radicals to dehydrogenate the organic molecules depends on the flexibility of the carbon backbones of the organic molecules as well as on the C-H bond strength. The dehydrogenation of n-alkyl chains mainly involves the H atoms of the β-carbon (leading to 1-alkene desorption). However, as the surface coverage decreases, the flexibility of the alkyl chains allows for the dehydrogenation of any methylene group and even the terminal methyl group of the long decyl layer. On the contrary, the rigid carbon backbone of the Si-CCCH₃ and Si-CHCHCH₃ moieties hinders the dehydrogenation of the terminal methyl group, which confers these layers a higher thermal stability. For all layers, the surface ends up mostly hydrogenated as Si-C bonds break and new Si-H bonds are formed during the dehydrogenation reactions.