We present a calculation of the thrust distribution in Higgs decays to quarks and gluons, \( H\to b\overline{b} \) , \( H\to c\overline{c} \) , and H → gg, including the resummation of large logarithmic corrections that arise in the two-particle limit at next-to-next-to-leading logarithmic (NNLL) accuracy, and match it to fixed-order results for three-particle decays at next-to-next-to-leading order (NNLO) in the strong coupling. The resummation is performed analytically within the Ares framework and combined with the fixed-order results using the logR matching technique. The fixed-order calculation is carried out numerically with the NNLOjet parton-level event generator, using the antenna subtraction method. We perform detailed cross-validation in the two-particle region, demonstrating that the expansion of the NNLL resummed result correctly reproduces the logarithmic structure of the fixed-order calculation to \( \mathcal{O}\left({\alpha}_{\textrm{s}}^3\right) \) , up to a predictable N3LL term at \( \mathcal{O}\left({\alpha}_{\textrm{s}}^3L\right) \) . In addition to providing the first NNLO+NNLL accurate predictions for the thrust distribution in Higgs decays to quarks and gluons, we analytically extract the \( \mathcal{O}\left({\alpha}_{\textrm{s}}^2\right) \) hard-virtual correction c2 and the \( {\alpha}_{\textrm{s}}^3L \) term G31 in both the \( H\to q\overline{q} \) (q = b, c) and H → gg decay channels.
