The Bloch wave functions have been playing a crucial role in the diagnosis of topological phases in noninteracting systems. However, the Bloch waves are no longer applicable in the presence of finite Coulomb interaction and alternative approaches are needed to identify the topological indices. In this paper, we focus on three-dimensional higher-order topological insulators protected by ๐ถ 4 ๐ symmetry and show that the topological index can be computed through eigenstates of inverse Green's function at zero frequency. If there is an additional ๐ 4 rotoinversion symmetry, the topological index ๐ 3 can be determined by eigenvalues of ๐ 4 at high symmetry momenta, similar to the Fu-Kane parity criterion. We also discuss the robustness of the hinge states under deformation of boundary and show that these hinge states exist even when boundary is smooth and without a sharp hinge. Finally, we discuss the realization of this higher-order topological phase, which includes the higherorder topological insulators in tetragonal lattice structure with ๐ถ 4 ๐-preserving magnetic order, and higher-order topological superconductors with ๐ + ๐๐-wave pairing.