fix errors from new version of QuantumToolbox.jl

Author: ytdHuang

QuantumToolbox.jl/time_evolution/kerr.qmd

@@ -1,7 +1,7 @@
 ---
 title: "Kerr nonlinearities"
 author: Li-Xun Cai
-date: 2025-02-08  # last update (keep this comment as a reminder)
+date: 2025-05-17  # last update (keep this comment as a reminder)
 
 engine: julia
 ---
@@ -126,10 +126,10 @@ Since we are considering unitary dynamics, i.e., no dissipation, the dynamics ar
 
 Note that if the keyword argument `e_ops` is not supplied to [`mesolve`](https://qutip.org/QuantumToolbox.jl/stable/resources/api#QuantumToolbox.sesolve), the returned `result` contains the state at each time point at field `states`. Consult the [user guide to `TimeEvolutionSol`](https://qutip.org/QuantumToolbox.jl/stable/users_guide/time_evolution/solution) for more details.
 ```{julia}
-ψ0 = coherent(N, 2.0)
+ρ0 = coherent_dm(N, 2.0)
 tlist = 0:0.01:(2*π / χ)
 
-result = mesolve(H, ψ0, tlist)
+result = mesolve(H, ρ0, tlist)
 ```
 
 We first check the expectation value dynamics of the number operator `n` with the two visualization functions we defined previously.

QuantumToolbox.jl/time_evolution/lowrank.qmd

@@ -1,7 +1,7 @@
 ---
 title: Low rank master equation
 author: Luca Gravina
-date: 2025-04-14  # last update (keep this comment as a reminder)
+date: 2025-05-17  # last update (keep this comment as a reminder)
 
 engine: julia
 ---
@@ -90,7 +90,7 @@ M = latt.N + 1; # Number of states in the LR basis
 Since we will take as initial state for our dynamics the pure state with all spins pointing up, the initial low-rank basis must include at least this state.
 
 ```{julia}
-ϕ = Vector{QuantumObject{KetQuantumObject,Dimensions{M - 1,NTuple{M - 1,Space}},Vector{ComplexF64}}}(undef, M)
+ϕ = Vector{QuantumObject{Ket,Dimensions{M - 1,NTuple{M - 1,Space}},Vector{ComplexF64}}}(undef, M)
 ϕ[1] = kron(fill(basis(2, 1), N_modes)...)
 ```
 
@@ -182,7 +182,7 @@ function f_entropy(p, z, B)
 
     mul!(C, z, sqrt(B))
     mul!(σ, C', C)
-    return entropy_vn(Qobj(Hermitian(σ), type=Operator), base=2)
+    return entropy_vn(Qobj(Hermitian(σ), type=Operator()), base=2)
 end;
 ```