A search for long-lived particles (LLPs) decaying in the CMS muon detectors is presented. A data sample of proton-proton collisions at √𝑠 =13  TeV corresponding to an integrated luminosity of 138  fb−1, recorded at the LHC in 2016–2018, is used. The decays of LLPs are reconstructed as high multiplicity clusters of hits in the muon detectors. In the context of twin Higgs models, the search is sensitive to LLP masses from 0.4 to 55 GeV and a broad range of LLP decay modes, including decays to hadrons, 𝜏 leptons, electrons, or photons. No excess of events above the standard model background is observed. The most stringent limits to date from LHC data are set on the branching fraction of the Higgs boson decay to a pair of LLPs with masses below 10 GeV. This search also provides the best limits for various intervals of LLP proper decay length and mass. Finally, this search sets the first limits at the LHC on a dark quantum chromodynamic sector whose particles couple to the Higgs boson through gluon, Higgs boson, photon, vector, and dark-photon portals, and is sensitive to branching fractions of the Higgs boson to dark quarks as low as 2 ×10−3.