Finite-temperature thermodynamic theory is developed for equilibrium arrays of two-dimensional monolayer-high islands in heteroepitaxial systems at submonolayer coverage. It is shown that the entropy contribution to the Helmholtz free energy of the system favors formation of smaller islands at higher temperatures which results in a decrease of the average number of atoms in the islands (the island volume) with temperature. The characteristic temperature T char , at which the island volume is significantly decreased compared to its value at T ¼ 0 K, is found to be far below the characteristic energy of the island formation and to lie in a region of several hundreds of K. Such a temperature dependence can be the basis for decisive experiments aimed at distinguishing between thermodynamic and kinetic effects in the formation of arrays of 2D islands. Results of high resolution electron microscopy and photoluminescence spectroscopy of a submonolayer InAs/GaAs(001) system are in agreement with the theory. It confirms that the formation of a submonolayer array of InAs/GaAs(001) islands is predominantly influenced by thermodynamics.The recent breakthrough in quantum dot (QD) technology, allowing a new generation of optoelectronic devices such as semiconductor diode lasers, relies on effects of selforganization at semiconductor surfaces [1]. During the initial years of QD research its main objective was to study basic effects of self-organization, to investigate the optical properties of QD arrays, and to demonstrate the possibility of fabricating QD lasers. To ensure further progress in the fabrication of QD heterostructures one has to master the way of controlling the process of self-organized growth in order to intentionally tune the morphology of QD arrays. To achieve this goal it is necessary to understand the formation mechanisms of QD arrays for any heteroepitaxial system.There exist two basic classes of theoretical models of self-organized formation of QD arrays: thermodynamic models, which describe an array of strained islands as the equilibrium state of a heteroepitaxial system [2,3], and kinetic models, which emphasize strain-driven barriers for adatoms to attach to islands [4,5] or for a new atomic layer to nucleate on facets [6]. The barriers eventually lead to the self-limited growth of the islands. In this connection, the relative importance of thermodynamic and kinetic effects in the self-organized formation of arrays of strained islands still remains a highly debated issue [7]. To elucidate the relative roles of thermodynamic and kinetic effects in