Gedämpfter harmonischer Quantenoszillator - Entwicklung des kohärenten Zustands

Ich versuche, die folgende Hauptgleichung zu lösen (auch ähnlich dem gedämpften harmonischen Quantenoszillator):

$$\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)$$ mit kohärentem Anfangszustand: $\hat{\rho}(0) = |\alpha\rangle\langle\alpha|$. Meine Idee war, die Sudarshan-Funktion und die Gilmore D-Algebra zu verwenden, um eine Differentialgleichung zu schreiben. Der erste Schritt besteht darin, den Dichtematrixoperator in kohärente Zustandsbasis zu zerlegen: $$\hat{\rho}(t) = \int d^2 \beta\ P(\beta, t) |\beta\rangle \langle \beta|$$ und handle mit Operatoren, die in der Anfangsgleichung erscheinen:

$$\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)$$ Schließlich erhalten wir: $$\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)$$

Ich bin mir ziemlich sicher, dass die Ableitung korrekt ist, da die Differentialgleichung die Einheitsspur beibehält, dh

$$\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$$ Meine Idee war nun, die Potenzierung zu verwenden und die fast endgültige Lösung wie folgt zu schreiben: $$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)$$ und dann von der Definition der Delta-Funktion $$\delta^{(2)}(\alpha - \beta) = \frac{1}{\pi^2}\int d^2\eta\ e^{-i\eta^*(\alpha^* - \beta^*)}e^{-i\eta(\alpha - \beta)}$$ ich kann schreiben $$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})$$ Den letzten Schritt finden Sie hier . Am Ende bekomme ich $$\rho(t) = e^{\Gamma t}|\alpha e^{-t\Gamma/2}\rangle\langle \alpha e^{-t\Gamma/2}|$$ Das Problem ist, dass die Spur nicht erhalten bleibt, wenn sie ins Unendliche wächst. Irgendeine Idee, wo ich einen Fehler gemacht habe?