Фиктивные индексы и ограничения в системе Mathematica

Я хотел бы, чтобы Mathematica упрощала выражение в соответствии с некоторыми ограничениями. Например:

M = p.k p.k

учитывая ограничение:

p(mu) * p(nu) = \delta(mu,nu)

должен дать

M = k^2

Но я не знаю, как указать фиктивные индексы, и я также не совсем уверен, где я должен разместить само ограничение. В аргументе Simplify?