The article is focused on recovering of two-phase microstructures of effective properties for 2D elasticity (see e.g., ref. [1]) predicted by two, stress-based topology optimization methods: isotropic material design (IMD) and cubic material design (CMD)-see e.g., refs. [2-5]. In these articles, the reconstruction of microstructures in subdomains where the auxetic properties [6-13] of the optimal, nonhomogeneous isotropic or cubic elastic material appear has been undertaken. In the present article, being an extension of the previous work, [14] the homogenization method based on parametric description of the shape of the domain filled by the given isotropic material inside the representative volume elements (RVEs), instead of the combinatorial search among the admissible microstructures, is implemented. Within the 2D setting discussed, the local cubic symmetry of the material reduces to the local symmetry of a square. The main goal is to recover isotropic microstructures or microstructures of symmetry of a square at arbitrary points of the body in a relatively simple way, considering the heterogeneity of the distribution of optimal elastic moduli. In the case of an isotropic material, the results obtained in ref. [14] are used as a starting point to the further search of the auxetic microstructures, realizing the Cherkaev-Gibiansky bounds. [15] In the case of 2D models of materials showing auxetic properties, a series of works devoted to their research at the level of regular atomic structures (e.g., hard cyclic hexamers [HCHs] or hard cyclic tetramers and many others) were particularly important in the development of the theory of such materials-see refs. [16-24]. Especially important in this respect is the ref. [16], in which a constant thermodynamic tension Monte Carlo (MC) simulations revealed the first planar, isotropic (and chiral) thermodynamically stable structure of negative Poisson's ratio in tilted phase formed by HCHs. In ref. [17], the exact solution for negative Poisson's ratio at high densities and in the static limit (i.e., at zero temperature) was obtained for planar, isotropic (and chiral) structure of cyclic hexamers interacting through nearest-neighbouring site-site inversepower potential. In ref. [18], MC simulations performed on hard, cyclic tetramers revealed an interesting feature of such a very simple model of molecules composed of four identical hard disks: depending on its shape, all possible material behavioursregarding negative Poisson's ratio could arise and thermodynamically stable auxetic (or partially auxetic) chiral structures of square symmetry could be spontaneously formed. Analyzing in ref. [19], other 2D cyclic trimers model of triatomic molecules in which "atoms" are distributed on vertices of equilateral triangles, the authors have found the exact solution at zero temperature for the planar auxetic structure which is isotropic but not chiral (soft cyclic trimers). The both chiral and nonchiral molecular models were studied in ref. [20], where the influence of molecular geome...