We have implemented lattice-based kMC simulations of Ge/Si(001) heteroepitaxy. Surface reconstructions and mismatch strain energy are explicitly considered. Stranski–Krastanov growth with wetting layer and 3D island formation is reproduced. We observe shape transition from rectangular to square with increasing coverage. Island size evolution from bimodal to unimodal distribution is reproduced.