Error running master equation in QuantumOptics.jl

[apkille]

Hi @DominicDams, it looks like the issue is that ρ₀ is a DenseOperator wrapped around an array of type Matrix{Int64} rather than Matrix{ComplexF64}. Try setting

ρ₀ = DenseOperator(b, ComplexF64[1 0; 0 1])

and see if that works.

@DominicDams

σ₋ = sigmam(b)#DenseOperator(b,[ 0 1; 0 0])
σ₊ =sigmap(b)# DenseOperator(b, [ 0 0; 1 0])
σz = sigmaz(b)#DenseOperator(b, [1 0 ; 0 -1])
σx = sigmax(b)#DenseOperator(b,[0 1;1 0])

On a related and slightly pedantic note: the comments here are incorrect. If you run any of these lines of code in the REPL, you’ll see that you are creating an object of type Operator wrapped around a matrix with complex entries.

julia> σ₋ = sigmam(b)
Operator(dim=2x2)
  basis: Spin(1/2)sparse([2], [1], ComplexF64[1.0 + 0.0im], 2, 2)