A search for exotic Higgs boson decays to light pseudoscalars in the final state of two muons and two $\tau$ leptons is performed using proton-proton collision data recorded by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Masses of the pseudoscalar boson between 15.0 and 62.5 GeV are probed, and no significant excess of data is observed above the prediction of the standard model. Upper limits are set on the branching fraction of the Higgs boson to two light pseudoscalar bosons in different types of two-Higgs-doublet models extended with a complex scalar singlet.
Expected and observed 95% CL upper limits on (sigma(pp->h)/sigma(pp->hSM)) * B(h -> aa -> mumutautau) as a function of m(a), where h(SM) is the Higgs boson of the standard model, h is the observed particle with mass of 125 GeV, and a denotes a light Higgs-like state, as obtained from the 13 TeV data.