of energy we will work with dimensionless quantities: energy $\tilde H = H/J_2$, temperature $\tilde T = 
k_BT/J_2$, and parameter $\tilde K = K_2/J_2$ (in what follows, the tilde will be omitted).

Simulation of faceting {#s3}
======================

We will discus here results of the annealing process investigated by Monte Carlo (MC) simulation in canonical ensemble. Let us start with some details concerning the simulation method. Configurations of BCCSOS model are generated via the classic Metropolis algorithm. A new configuration in Markov chain is generated from the previous configuration by moving one B atom from a site *i* to a site *j*, what is equivalent to the following change of heights: $(h_i, h_j)\rightarrow (h_i-3, 
h_j+3)$. It is important to note that such change of heights is possible only if $h_i$ is the local maximum ( with respect to the nearest neighbors) and $h_j$ is the local minimum because of the restriction on column heights. Let us remember one of the model assumption that adsorbate atoms always stay on the column tops as we study reconstruction of B surface covered by the physical monolayer of A material.

![ Temperature dependence of the surface energy during the warming up cycles (triangles) and the 