The quadrupole insulator, a high-order topological insulator, with on-site Hubbard interaction is numerically studied by large-scale projector quantum Monte Carlo (PQMC) simulations. The Green's function formalism is successfully used to characterize topological properties in interacting quadrupole insulators for the first time. We find that the topological quadrupole insulator is stable against weak interactions and turns into a trivial antiferromagnetic (AFM) insulator by a continuous topological phase transition (TPT) for strong interactions. The critical exponents related to the TPT are estimated to be ν = 0.67(4), β = 0.40(2), which are distinct from those of the known AFM transitions and suggest a new universality class.The bulk-edge correspondence [1] is one of the most fundamental concepts for understanding topological insulators (TIs) in d-dimensional systems which possess gapless modes in their (d − 1)-dimensional boundaries [2][3][4]. Recently, topological crystalline insulators generalize this bulk-edge correspondence to high-order topological insulators (HOTIs) [5][6][7][8]. In addition, the concept of HOTIs is also generalized to high-order topological superconductors [9-13]. The predicted Majorana corner/hinge states in such systems are promising building blocks (qubits) for realizing quantum computation [14,15].Attracted by the novel properties and potential applications of HOTIs, many theoretical works have been done to predict HOTIs in realistic systems, such as quasicrystals, sonic crystals, non-hermitian systems, and solid state systems [16][17][18][19][20][21][22][23]. And some experimental works have even shown the direct observation of HOTIs in photonic/phononic crystals, electrical circuits as well as crystalline solids, such as bismuth [24][25][26][27][28][29].Unlike bosonic systems, fermionic systems with spin degrees of freedom are subjected to the Coulumb repulsion which may give rise to many exotic phenomena [30][31][32][33][34][35][36][37]. However, there are only a few studies which discuss correlation effects in HOTIs [38][39][40]. Unbiased numerical research on interacting HOTIs is almost a blank and is the main motivation of this work.