Temperature dependence of the specific heat C at various fields for the model with J_ab = −0.2, ε = 0.25. The system size is 24 × 24 × 24. The peak in C indicates a thermal transition to the magnetically ordered phase, which is used to determine the phase boundary in Fig. 1(b).