  We have performed Kinetic Monte Carlo simulation to study the effect of
diffusion anisotropy and bonding anisotropy on island formation at different
temperatures during the sub-monolayer film growth in Molecular Beam Epitaxy. We
use simple cubic solid on solid model and event based Bortz, Kalos and
Labowitch (BKL) algorithm on Kinetic Monte Carlo method to simulate the
physical phenomena. We have found that surface anisotropy has no significant
role on island elongation however it influences on the island morphology,
growth exponent and island size distribution. Elongated islands were obtained
when bonding anisotropy was included.
