greater pyramids are formed on the surface. Pyramids reach maximal sizes just below the temperature of the phase transition where dependence of sizes on annealing time is the strongest. It worth noting that dependence of facets size on temperature qualitatively agrees with the experimental results [@mad94].

![ Temperature dependence of the $\delta h^2$ during the heating and the cooling cycles for 3 different annealing times. Each label contains symbols for heating, cooling and the number of MC steps[]{data-label="figure6"}](fig6){width="7cm"}

The value of $\delta h^2$ rapidly decreases above the phase transition where very weak dependence on temperature is observed. In the cooling down process large pyramids are formed just below the phase transition temperature, and further decreasing of temperature practically does not change the size of pyramids.

We check dependence of surface growth on the initial temperature $T_0$ and we find that the behaviour of $\delta h^2$ near the transition point is very similar to that in Fig. \[figure6\] for process which starts at $T_0=0.3$ or $T_0=1.0$. Finally we analyze the influence of the temperature increment, $\Delta T$, on the surface energy and the $\delta h^2$ near the phase transition for $\tau=3\times 10^5$. We found that for $\Delta 
T=0.01$ the hysteresis loop 