I'm analysing an open quantum system where the optical cavity is interacting with a quantum dot. I modelled the system using cavity QED and used Lindblad master equation to model the system.
Now I have obtained an equation for the cavity field operators $\hat{a}^\dagger$ and $\hat{a}$. Now I want to calculate first-order correlation function of field operators.
$$g^{(1)}= \langle \hat{a}^\dagger(t)\hat{a}^\dagger(t+\tau) \hat{a}(t+\tau)\hat{a}(t) \rangle$$
Firstly, I am new to quantum optics and can't figure out how to calculate above analytically or plot it numerically. Secondly, I don't know how I can get the time dependence of Schrodinger operators $\hat{a}^\dagger$ and $\hat{a}$
If you can provide me with useful tools or some hints it will be useful.