@Aaboud:2018zhk] and found them to be also compatible with the SM prediction. A measurement of the couplings of the boson to second generation leptons [@ATLAS-CONF-2019-028; @Sirunyan:2018hbu] and quarks is the next target.

In this paper, we focus on the search for bosons decaying to $\PQc\PAQc$, a charm quark-antiquark pair. The boson to charm quark Yukawa coupling $y_{{\PQc}}$ can be significantly modified by physics beyond the SM [@Ghosh:2015gpa; @Botella:2016krk; @Harnik:2012pb; @Altmannshofer:2016zrn]. Allowed range for unobserved decays constitutes the bound on the charm quark Yukawa coupling. Indirect constraints on $y_{{\PQc}}$ obtained from a global fit to existing boson data result in an upper bound on $\kappa_{\PQc} \equiv y_{\PQc} / y_{\PQc}^{\text{SM}}$ of 6.2 [@Perez:2015aoa] at 95% confidence level (), assuming the absence of non-SM production mechanisms. A direct measurement of this process is extremely challenging at a hadron collider. The branching fraction of this process according to SM computations, $\mathcal{B}({\ensuremath{\PH\to\ccbar}\xspace})=0.0288^{+0.0016}_{-0.0006}$ [@deFlorian:2016spz], is a factor 20 smaller than that of $\PH \to \bbbar$, and there is a very large background from SM processes comprised uniquely of jets produced through the strong interaction, referred to as quantum chromodynamics (QCD) multijet events. Results from direct searches for ${\ensuremath{\PH\to\ccbar}\xspace}$ at the LHC in the ($\PZ\to\ell\ell$, $\ell=\Pe$ or 