Randomness and frustration are considered to be the key ingredients for the existence of spin glass (SG) phase. In a canonical system, these ingredients are realized by the random mixture of ferromagnetic (FM) and antiferromagnetic (AF) couplings. Th
e study by Bartolozzi {it et al.} [Phys. Rev. B{bf 73}, 224419 (2006)] who observed the presence of SG phase on the AF Ising model on scale free network (SFN) is stimulating. It is a new type of SG system where randomness and frustration are not caused by the presence of FM and AF couplings. To further elaborate this type of system, here we study Heisenberg model on AF SFN and search for the SG phase. The canonical SG Heisenberg model is not observed in $d$-dimensional regular lattices for ($d leq 3$). We can make an analogy for the connectivity density ($m$) of SFN with the dimensionality of the regular lattice. It should be plausible to find the critical value of $m$ for the existence of SG behaviour, analogous to the lower critical dimension ($d_l$) for the canonical SG systems. Here we study system with $m=2,3,4$ and $5$. We used Replica Exchange algorithm of Monte Carlo Method and calculated the SG order parameter. We observed SG phase for each value of $m$ and estimated its corersponding critical temperature.
The edge-cubic spin model on square lattice is studied via Monte Carlo simulation with cluster algorithm. By cooling the system, we found two successive symmetry breakings, i.e., the breakdown of $O_h$ into the group of $C_{3h}$ which then freezes in
to ground state configuration. To characterize the existing phase transitions, we consider the magnetization and the population number as order parameters. We observe that the magnetization is good at probing the high temperature transition but fails in the analysis of the low temperature transition. In contrast the population number performs well in probing the low- and the high-$T$ transitions. We plot the temperature dependence of the moment and correlation ratios of the order parameters and obtain the high- and low-$T$ transitions at $T_h = 0.602(1)$ and $T_l=0.5422(2)$ respectively, with the corresponding exponents of correlation length $ u_h=1.50(1)$ and $ u_l=0.833(1)$. By using correlation ratio and size dependence of correlation function we estimate the decay exponent for the high-$T$ transition as $eta_h=0.260(1)$. For the low-$T$ transition, $eta_l = 0.267(1)$ is extracted from the finite size scaling of susceptibility. The universality class of the low-$T$ critical point is the same as the 3-state Potts model.
