Pseudoscalar decays into lepton pairs from rational approximants

The pseudoscalar decays into lepton pairs $P\rightarrow\overline{\ell}\ell$ are analyzed with the machinery of Canterbury approximants, an extension of Pad\'e approximants to bivariate functions. This framework provides an ideal model-independent approach to implement all our knowledge of the pseudoscalar transition form factors driving these decays, can be used for data analysis, and allows to include experimental data and theoretical constraints in an easy way, and determine a systematic error. We find that previous theoretical estimates for these branching ratios have underestimated their theoretical uncertainties. From our updated results, the existing experimental discrepancies for $\pi^0\rightarrow e^+e^-$ and $\eta\rightarrow \mu^+\mu^-$ channels cannot be explained unless the doubly-virtual transition form factors behavior -not yet measured- are out of theoretical expectations, which is an interesting result both for anomalous magnetic moments of leptons, and for physics beyond the standard model.