We test in $(A_{n-1},A_{m-1})$ Argyres-Douglas theories with $mathrm{gcd}(n,m)=1$ the proposal of Songs in arXiv:1612.08956 that the Macdonald index gives a refined character of the dual chiral algebra. In particular, we extend the analysis to higher rank theories and Macdonald indices with surface operator, via the TQFT picture and Gaiotto-Rastelli-Razamats Higgsing method. We establish the prescription for refined characters in higher rank minimal models from the dual $(A_{n-1},A_{m-1})$ theories in the large $m$ limit, and then provide evidence for Songs proposal to hold (at least) in some simple modules (including the vacuum module) at finite $m$. We also discuss some observed mismatch in our approach.