Function decoherence 

pub fn decoherence(rho: &mut DensityMatrix, rate: f64, dt: f64)

Apply decoherence to a density matrix: off-diagonal elements decay. Models Lindblad-type decoherence: rho_ij *= exp(-rate * dt) for i != j.