A search for a light pseudoscalar Higgs boson (a) decaying from the 125 GeV (or a heavier) scalar Higgs boson (H) is performed using the 2016 LHC proton-proton collision data at root s = 13 TeV, corresponding to an integrated luminosity of 35.9 fb(-1), collected by the CMS experiment. The analysis considers gluon fusion and vector boson fusion production of the H, followed by the decay H -> aa -> mu mu tau tau, and considers pseudoscalar masses in the range 3.6 aa -> mu mu tau tau, down to 1.5 (2.0)x10(-4) for m(H) = 125 (300) GeV.
Model-dependent limits on B(H -> aa) are set within the context of two Higgs doublets plus singlet models, with the most stringent results obtained for Type-III models. These results extend current LHC searches for heavier a bosons that decay to resolved lepton pairs and provide the first such bounds for an H boson with a mass above 125 GeV.