qutip

To simulate open quantum system dynamics, define your Hamiltonian and collapse operators, then use the mesolve solver to numerically integrate the Lindblad master equation over a time grid and extract expectation values.

The mesolve solver directly integrates the density matrix master equation for mixed states, while mcsolve uses Monte Carlo quantum trajectories to simulate individual stochastic jump events for measurement-driven dissipation behavior.

Yes, modeling decoherence requires manually constructing collapse operators representing system losses, which you pass alongside the Hamiltonian to solvers like mesolve to predict photon number decay and steady-state observables.

Compute steady states and correlation functions by configuring your Hamiltonian and collapse operators, running the appropriate QuTiP dynamics solver, and then analyzing the resulting state objects to extract spectra and observables.

The best way to visualize Wigner functions and Bloch spheres is to run time evolution solvers like sesolve or mesolve, then use built-in analysis functions to plot the resulting quantum states directly.

Use quantum trajectories via mcsolve when you need to simulate individual stochastic measurement-driven jump events, whereas master equation solvers like mesolve are better for tracking continuous density matrix evolution.