We test in (A n−1 , A m−1 ) Argyres-Douglas theories with gcd(n, m) = 1 the proposal of Song's in [1] 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-Razamat's 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 Song's 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.