We study the multicritical behavior for the semimetal-insulator transitions on graphenes honeycomb lattice using the Gross-Neveu-Yukawa effective theory with two order parameters: the SO(3) (Heisenberg) order parameter describes the antiferromagnetic
transition, and the $mathbb{Z}_2$ (Ising) order parameter describes the transition to a staggered density state. Their coupling induces multicritical behavior which determines the structure of the phase diagram close to the multicritical point. Depending on the number of fermion flavors $N_f$ and working in the perturbative regime in vicinity of three (spatial) dimensions, we observe first order or continuous phase transitions at the multicritical point. For the graphene case of $N_f=2$ and within our low order approximation, the phase diagram displays a tetracritical structure.
In numerical simulations, spontaneously broken symmetry is often detected by computing two-point correlation functions of the appropriate local order parameter. This approach, however, computes the square of the local order parameter, and so when it
is {it small}, very large system sizes at high precisions are required to obtain reliable results. Alternatively, one can pin the order by introducing a local symmetry breaking field, and then measure the induced local order parameter infinitely far from the pinning center. The method is tested here at length for the Hubbard model on honeycomb lattice, within the realm of the projective auxiliary field quantum Monte Carlo algorithm. With our enhanced resolution we find a direct and continuous quantum phase transition between the semi-metallic and the insulating antiferromagnetic states with increase of the interaction. The single particle gap in units of the Hubbard $U$ tracks the staggered magnetization. An excellent data collapse is obtained by finite size scaling, with the values of the critical exponents in accord with the Gross-Neveu universality class of the transition.
