Явное определение оператора энергии в модели Изинга

Я смоделировал несколько двухмерных моделей Изинга при критической температуре на треугольной решетке и теперь пытаюсь проверить правильность корреляционных функций. Я уже сделал это для оператора спина ($\sigma$) и я хотел бы сделать то же самое для оператора локальной энергии, который должен удовлетворять следующей пространственной зависимости:

$$\langle\epsilon (x) \epsilon (y) \rangle \propto \frac{1}{(x-y)^2}$$

(масштабирование веса 1)

но я не уверен в определении этого оператора, я бы предположил что-то вроде:

$$\epsilon(x)=-J\,\sigma(x)\sum_{\text{neighbours}}\sigma$$

где J — критическая связь ($=\frac{1}{T_{\text{crit}}}$)

Это правильно?

Редактировать: новое исследование заставило меня отказаться от моего первого предложения : этот оператор энергии должен быть странным под$\text{Z}_2$симметрия системы$\sigma \leftrightarrow - \sigma$что явно не относится к моему предыдущему предположению.