No abstract

Using atomic-scale molecular dynamics and energy minimisation techniques with semi-empirical molecular mechanics 3 potential energy functions, we consider the adsorption of a C 60 molecule on pentacene in the known thin film and bulk phases and a series of hypothetical, variably angled, pentacene structures. The thin film phase has a more energetically isotropic surface than the bulk, and exhibits diffusion coefficients that are twice as high, potentially leading to preferable characteristics for ordered film growth. For the variably angled pentacene structures, where the long axis is parallel to the substrate, a relationship was discovered between the angle that the pentacene short axis makes with the surface normal (controlled by the underlying substrate) and the adsorption characteristics of C 60 . There is a transition of the dominant energy minima from between the pentacene rows at low values of w 1 to within the rows at high values of w 1 , where w 1 is the angle the bottom pentacene short axis makes with the surface normal. This shift implies that the likelihood of forming C 60 nanowires on pentacene is greater at extreme values of w 1 for which there would be a clear preference for C 60 to be located between a row or within a row, rather than at intermediate values of w 1 , where there is no clear preference.

The tendency for C(60) nanowires to persist on two monolayers of recumbent pentacene is studied using molecular dynamics (MD) simulations. A review of existing experimental literature for the tilt angle adopted by pentacene on noble metal surfaces shows that studies cover a limited range from 55° to 90°, motivating simulation studies of essentially the entire range of tilt angles (10°-90°) to predict the optimum surface tilt angle for C(60) nanowire formation. The persistence of a 1D nanowire depends sensitively on this tilt angle, the amount of initial tensile strain, and the presence of surface step edges. At room temperature, C(60) nanowires oriented along the pentacene short axes persist for several nanoseconds and are more likely to occur if they reside between, or within, pentacene rows for ϕ ≤ ∼60°. The likelihood of this persistence increases the smaller the tilt angle. Nanowires oriented along the long axes of pentacene molecules are unlikely to form. The limit of stability of nanowires was tested by raising the temperature to 400 K. Nanowires located between pentacene rows survived this temperature rise, but those located initially within pentacene rows are only stable in the range ϕ(1) = 30°-50°. Flatter pentacene surfaces, that is, tilt angles above about 60°, are subject to disorder caused by C(60) molecules "burrowing" into the pentacene surface. An initial strain of 5% applied to the C(60) nanowires significantly decreases the likelihood of nanowire persistence. In contrast, any appreciable surface roughness, even by half a monolayer in height of a third pentacene monolayer, strongly enhances the likelihood of nanowire formation due to the strong binding energy of C(60) molecules to step edges.

Using atomic-scale Molecular Dynamics (MD) and energy minimization techniques in conjunction with semi-empirical MM3 potential energy functions, we consider the adsorption of a C60 molecule on a series of hypothetical pentacene structures that vary only in the tilt of the angle that the short axis of the pentacene molecules makes with the underlying surface (the long axis lying essentially flat, as on a metal substrate). Important relationships were discovered between the angle adopted by the short axis of pentacene on the surface, φ1, and the adsorption and diffusion characteristics of C60. Static energy calculations show that there is a transition of the deepest energy minima from between the pentacene rows at low values of φ1 to within the rows at high values of φ1, where φ1 is the angle the pentacene short axis makes with the surface. MD confirms this trend by the predominant residence locations at the extreme φ1 values. Furthermore, MD results suggest that the C60 traverses the pentacene surface in the east-west direction for lower φ1 values (φ1 ≤ 40°) and in the north-south direction for higher φ1 values (φ1 ≥ 70°). Taking both static and dynamic results together, the most favorable tilt angles for mono-directional nanowire growth should occur between 70° and 80° off-normal.

A novel multilattice kinetic Monte Carlo algorithm is developed for heteroepitactical growth of a hexagonal lattice material (C60) on an oblique lattice material (pentacene). This algorithm captures the behavior of single molecule and small clusters of C60 molecules diffusing, clustering, and reorganizing as monolayers on the surface and switching lattices depending on their local environment. An extensive catalog of jump rates and energy barriers created from molecular dynamics simulations and molecular mechanics is used as the sole input to the method in order to follow the evolution of C60 growth on pentacene for time scales approaching a millisecond that are unattainable using molecular dynamics alone.