自旋跃迁模拟

输入初态密度矩阵，选择跃迁通道和旋转角度，计算终态密度矩阵 ρ' = UρU†