Obtaining all the density matrix elements as a function of time from mesolve

[ssau...@gmail.com]

Hi. All.

I am using qutip's mesolve for the open system dynamics of a system. I need to save all the density matrix elements as a function of time. mesolve only gives the option of returning the expectation value of a given operator given to it as an input. I am wondering if there is a way of having each element of the density matrix as a function of time without having to write an operator for each element and getting the expectation values individually. Please help if you have any suggestions.

Saumya Biswas

[Nathan Shammah]

Hi Saumya, 

Set the optional argument `store_states = True` in the options. 

For example with this snippet

```

opts = Options(store_states = True) 

result = mesolve(H, rho0, tlist, c_ops, e_ops, options = opts)

rho_t = result.states

```

Best wishes, 

Nathan 

[Saumya biswas]

Thanks Nathan!!

--
You received this message because you are subscribed to the Google Groups "QuTiP: Quantum Toolbox in Python" group.
To unsubscribe from this group and stop receiving emails from it, send an email to qutip+un...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.