Narayana numbers in invariants for interacting particle systems of Itoh-Mallows-Shepp

Abstract We consider an interacting particle system on star graphs. As in the case of the KdV equation, we have infinitely many invariants (here, martingale invariants). It enables us to obtain the limiting distribution of the Markov chain. Each of the martingale invariants is a homogeneous polynomial with coefficients of Narayana numbers. By a change of variables, the identity for the enumeration of plane unlabeled trees, which gives Narayana numbers, becomes the key identity to obtain the probability of death states.