We include the ηc meson into the η-η′-G mixing formalism
constructed in our previous work, where G represents the pseudoscalar
gluball. The mixing angles in this tetramixing matrix are constrained by
theoretical and experimental implications from relevant hadronic processes.
Especially, the angle between ηc and G is found to be about 11∘
from the measured decay widths of the ηc meson. The pseudoscalar glueball
mass mG, the pseudoscalar densities mqq,ss,cc and the U(1) anomaly
matrix elements associated with the mixed states are solved from the anomalous
Ward identities. The solution mG≈1.4 GeV obtained from the
η-η′-G mixing is confirmed, while mqq grows to above the pion
mass, and thus increases perturbative QCD predictions for the branching ratios
Br(B→η′K). We then analyze the ηc-mixing effects on charmonium
magnetic dipole transitions, and on the B→η(′)KS branching
ratios and CP asymmetries, which further improve the consistency between
theoretical predictions and data. A predominant observation is that the
ηc mixing enhances the perturbative QCD predictions for Br(B→η′K)
by 18%, but does not alter those for Br(B→ηK). The puzzle due to the
large Br(B→η′K) data is then resolved.Comment: 12 pages, version to appear in PR