We investigate low temperature properties of a random Ising model with 1J and 2aJ ͑a fi 1͒ bonds in two dimensions using a cluster heat bath method. It is found that the Binder parameters g L for different sizes of the lattice come together at almost the same temperature, implying the occurrence of the spin glass (SG) phase transition. From results of finite size scaling analyses, we suggest that the SG phase really occurs at low temperatures, which is characterized by a power law decay of spin correlations. [S0031-9007(97) Spin glasses have attracted great challenge for computational physics in these two decades. It is widely believed now in the bond-random Ising model that spin glass (SG) transitions occur at a finite, nonzero temperature T c fi 0 in three dimensions (3D) [1][2][3][4] and at zero temperature T c 0 in two dimensions (2D) [4][5][6][7][8]. Recently, the present authors [9] reexamined the SG phase transition of the 6J Ising model on a square lattice of L 3 L by means of an exchange Monte Carlo method [10] and found that the Binder parameters g L for L # 16 intersect at T fi 0. They also found that better finite-size scaling (FSS) fits of the spin glass susceptibility x SG are obtained when T c fi 0. These results imply the occurrence of the SG phase transition at T c fi 0. If so, it is quite interesting, because it disproves the belief of T c 0. However, there remain two problems which should be considered to suggest T c fi 0 in 2D. One is that g L for a smaller lattice almost saturates below a rather high temperature [7] and its saturation value slightly increases with L [9]. Therefore, it is difficult to see whether the intersection of g L for smaller L suggests the presence of the SG phase at T c fi 0 or is merely due to a finite size effect. The other is that it is still open whether or not the model really exhibits the nature of the SG phase at T , T c , because the estimated transition temperature T c is slightly lower than the lowest temperature which is reached in the simulation. The problems would be solved if we study the same model on bigger lattices at lower temperatures. The saturation of g L at rather high temperatures, however, may be removed if we treat an asymmetric random Ising model with 1J and 2aJ ͑a fi 1͒ bonds, because the energy gap of 2j1 2 ajJ in that model between the ground state and the lowest excitation state is much smaller than that of 4J in the 6J model [11], and, if the lattice is rather small, we may study equilibrium properties at any temperature using a cluster heat bath (CHB) method [12,13].In this Letter, we investigate low temperature properties of the asymmetric random Ising model on the square lattice of L 3 L ͑L # 18͒ using the CHB method. In fact, g L does not saturate down to a very low temperature. We find that as the temperature decreases, g L 's for different L meet at almost the same temperature and then increase together. This property rather resembles that of the 6J model in 3D in which the SG phase transition occurs at T c fi 0. We make the FSS an...