Hi
I believe that the examples, both in the tutorials, and the manual, have errors when dealing with dephasing.
In the manual for the 4.6.0 release (page 98), the dephasing operator is defined as: np.sqrt(gamma2) * sz:
c_ops = [np.sqrt(gamma1) * sm, np.sqrt(gamma2) * sz] # relaxation and dephasing
There should be a factor of 2
np.sqrt(gamma2/2) * sz
so that the decay of the off diagonal elements is T2=1/gamma2
Since the Lindbladt operator for dephasing is sqrt(gamma/2)*sz