Затухающий квантовый гармонический осциллятор - эволюция когерентного состояния

Я пытаюсь решить следующее основное уравнение (также похожее на затухающий квантовый гармонический осциллятор):

$$\frac{d\hat{\rho}}{dt} = \frac{\Gamma}{2}\left(2\hat{a}\hat{\rho}\hat{a}^{\dagger} - \hat{a}^{\dagger}\hat{a} \hat{\rho} - \hat{\rho}\hat{a}^{\dagger}\hat{a} \right)$$ с начальным когерентным состоянием: $\hat{\rho}(0) = |\alpha\rangle\langle\alpha|$. Моя идея заключалась в том, чтобы использовать функцию Сударшана и D-алгебру Гилмора для записи дифференциального уравнения. Первым шагом является декомпозиция оператора матрицы плотности в когерентном базисе состояния: $$\hat{\rho}(t) = \int d^2 \beta\ P(\beta, t) |\beta\rangle \langle \beta|$$ и действовать с операторами, входящими в исходное уравнение:

$$\hat{a}\hat{\rho}\hat{a}^{\dagger} = \int d^2 \beta\ |\beta|^2 P(\beta, t) |\beta\rangle \langle \beta|$$ $$\hat{a}^{\dagger}\hat{a} \hat{\rho} = \int d^2 \beta\ P(\beta,t)\beta\left(\beta^* + \frac{\partial}{\partial\beta} \right)|\beta\rangle \langle \beta| = \int d^2 \beta\ |\beta\rangle\langle\beta| \left(\beta^* - \frac{\partial}{\partial\beta} \right) \beta P(\beta,t)$$ $$\hat{\rho}\hat{a}^{\dagger}\hat{a} = \int d^2 \beta\ P(\beta,t) \beta^*\left(\beta + \frac{\partial}{\partial\beta^*} \right) |\beta\rangle \langle \beta| = \int d^2 \beta\ |\beta\rangle\langle\beta| \left(\beta - \frac{\partial}{\partial\beta^*} \right)\beta^*P(\beta,t)$$ Наконец мы получаем: $$\frac{\partial P(\beta, t)}{\partial t} = \frac{\Gamma}{2}\left(\beta \frac{\partial}{\partial \beta} + \beta^* \frac{\partial}{\partial \beta^* } +2\right) P(\beta, t)$$

Я почти уверен, что вывод правильный, потому что дифференциальное уравнение сохраняет единичный след, т.е.

$$\frac{d}{dt}\text{Tr}\{\hat{\rho}(t)\} = \int d^2\beta\ \frac{\partial P(\beta,t)}{\partial t} = \frac{\Gamma}{2}\int d^2\beta\ \left(\frac{\partial}{\partial \beta}\beta P(\beta,t) + \frac{\partial}{\partial \beta^*}\beta^* P(\beta,t) \right) = 0$$ Теперь моя идея заключалась в том, чтобы использовать возведение в степень и записать почти окончательное решение в виде: $$P(\beta, t) = \exp\left[t\frac{\Gamma}{2}\left(\beta \frac{\partial}{\partial \beta} + \beta^* \frac{\partial}{\partial \beta^*} +2\right)\right] \delta^{(2)}(\alpha - \beta)$$ а затем из определения дельта-функции $$\delta^{(2)}(\alpha - \beta) = \frac{1}{\pi^2}\int d^2\eta\ e^{-i\eta^*(\alpha^* - \beta^*)}e^{-i\eta(\alpha - \beta)}$$ я могу написать $$P(\beta,t) = e^{\Gamma t} \frac{1}{\pi^2} \int d^2\eta\ e^{-i\eta^*\alpha^2} e^{t\frac{\Gamma}{2} \beta^* \partial_{\beta^*}} e^{i\eta^* \beta^*} \times e^{-i\eta \alpha} e^{t\frac{\Gamma}{2} \beta \partial_{\beta}} e^{i\eta \beta} = e^{\Gamma t}\delta^{2}(\alpha - \beta e^{t\Gamma/2})$$ Последний шаг можно найти здесь . В конце я получаю $$\rho(t) = e^{\Gamma t}|\alpha e^{-t\Gamma/2}\rangle\langle \alpha e^{-t\Gamma/2}|$$ Проблема в том, что он не сохраняет след, когда он растет до бесконечности. Любая идея, где я сделал ошибку?