A fresh look into mbarc,b (mbarc,b) and precise f D(s) ,B(s) from heavy-light QCD spectral sum rules

Using recent values of the QCD (non-)perturbative parameters given in Table 1 and an estimate of the N3LO QCD perturbative contributions based on the geometric growth of the PT series, we re-use QCD spectral sum rules (QSSR) known to N2LO PT series and including all dimension-six NP condensate contributions in the full QCD theory, for improving the existing estimates of mbar_c,b and _{f D(s), B(s)} from the open charm and beauty systems. We especially study the effects of the subtraction point on "different QSSR data" and use (for the first time) the Renormalization Group Invariant (RGI) scale-independent quark masses in the analysis. The estimates [rigourous model-independent upper bounds within the SVZ framework] reported in Table 8: f_D /f_π = 1.56 (5) [ ⩽ 1.68 (1) ], f_B /f_π = 1.58 (5) [ ⩽ 1.80 (3) ] and f_Ds /f_K = 1.58 (4) [ ⩽ 1.63 (1) ], f_Bs /f_K = 1.50 (3) [ ⩽ 1.61 (3.5) ], which improve previous QSSR estimates, are in perfect agreement (in values and precisions) with some of the experimental data on f_D,Ds and on recent lattice simulations within dynamical quarks. These remarkable agreements confirm both the success of the QSSR semi-approximate approach based on the OPE in terms of the quark and gluon condensates and of the Minimal Duality Ansatz (MDA) for parametrizing the hadronic spectral function which we have tested from the complete data of the J / ψ and ϒ systems. The values of the running quark masses mbar_c (m_c) = 1286 (66) MeV and mbar_b (m_b) = 4236 (69) MeV from M_D,B are in good agreement though less accurate than the ones from recent J / ψ and ϒ sum rules.