“…Later, Xavier et al [9] considered the helical boundary conditions to calculate the conformal anomaly at the tricritical point with the estimation of the phase boundary. The nature of the phase transitions in the 2D BC model has been extensively examined also with various other numerical methods which include the Monte Carlo renormalization group [10,11], histogram method [12,13], Wang-Landau sampling [14][15][16][17][18][19], simulated * dongheekim@gist.ac.kr tempering [20], tensor renormalization group [21], replicaexchange Monte Carlo [22], and quite recently parallel multicanonical simulations [23]. While we focus on the pure BC model in this work, the effects of randomness have been also extensively studied [16,17,[24][25][26][27][28][29], where, for instance, the phase diagram was found to be shifted in presence of bond randomness at low temperatures [17].…”