Во-вторых, стоит заметить, что поля

являются обобщенными координатой и импульсом соответственно, потому для них справедливо соотношение
![$[\hat {\Psi} (\mathbf x, t), \frac{d}{dt}\hat {\Psi}(\mathbf y, t)] = i\delta (\mathbf x - \mathbf y )$ $[\hat {\Psi} (\mathbf x, t), \frac{d}{dt}\hat {\Psi}(\mathbf y, t)] = i\delta (\mathbf x - \mathbf y )$](https://dxdy.ru/math/fa5717d1517fea9687c7e7566c188e9c82.png)
. Это соотношение использовалось для выведения соотношения

(детали приведены в книгах Вайнберга и Грайнера). Однако далеко не для всех полей производная по времени к полю образует канонический импульс для самого поля. В качестве примера достаточно рассмотреть все поля со спином старше 0. [:)]. Потому для разных полей нужно подбирать соответствующие канонические импульсы, что для произвольного случая есть трудно решаемой (если не нереализуемой) задачей.
Между тем, тот же Вайнберг пишет, что результат для скалярного поля легко обобщается на случай полей произвольного спина. Жаль только, что деталей он не написал. [:)].
P.S. Да, видимо, можно допустить, что ничего страшного в том, что появляются индексы при рассмотрении поля произвольного спина, нет. Вы правы, мне кажется.