Scalar glueballs: Constraints from the decays into $\eta$ or $\eta'$

We study the mixing of the scalar glueball into the isosinglet mesons f0(1370), f0(1500), and f0(1710) to describe the two-body decays to pseudoscalars. We use an effective Hamiltonian and employ the two-angle mixing scheme for η and η′. In this framework, we analyze existing data and look forward to new data into η and η′ channels. For now, the f0(1710) has the largest glueball component and a sizable branching ratio into ηη′, testable at BESIII.