“…Classical HMFT computations with N =4 and N =16 clusters provide a picture consistent with Landau theory, where the Néel order melts through a second order transition into a C 4 symmetric plaquette-VBS characterized by 4-spin strongly resonating plaquettes through a second order transition, before entering the CAF through a first order transition [35]. Interestingly, upon introducing a third nearest neighbor Heisenberg interaction, the plaquette-VBS order has been argued to become stronger [50] and, more recently, to be a higher-order symmetryprotected topological phase [72]. Specifically, we perform ideal classical simulations of Q-HMFT on the J 1 -J 2 model (7), i.e.…”