From 2f0690926204c9f35a6d45635ba7165c09974d69 Mon Sep 17 00:00:00 2001 From: Alberto Mercurio Date: Sat, 7 Oct 2023 19:36:30 +0200 Subject: [PATCH] Fixed problem in the new liouvillian_generalized --- src/time_evolution/time_evolution.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/time_evolution/time_evolution.jl b/src/time_evolution/time_evolution.jl index 162815f1..83080ee2 100644 --- a/src/time_evolution/time_evolution.jl +++ b/src/time_evolution/time_evolution.jl @@ -735,8 +735,8 @@ function liouvillian_generalized(H::QuantumObject{<:AbstractArray, OperatorQuant end # Filter in the Hilbert space - σ = σ_filter == nothing ? 20 * sum(abs2, X_op) / length(X_op) : σ_filter - F1 = QuantumObject(gaussian.(Ω, 0, 2*κc), dims=dims) + σ = σ_filter == nothing ? 10 * sum(abs2, X_op) / length(X_op) : σ_filter + F1 = QuantumObject(gaussian.(Ω, 0, σ), dims=dims) F1 = dense_to_sparse(F1, tol) # Filter in the Liouville space @@ -744,7 +744,7 @@ function liouvillian_generalized(H::QuantumObject{<:AbstractArray, OperatorQuant Ω1 = kron(Ω, M1) Ω2 = kron(M1, Ω) Ωdiff = Ω1 .- Ω2 - F2 = QuantumObject(gaussian.(Ωdiff, 0, 2*κc), SuperOperatorQuantumObject, dims) + F2 = QuantumObject(gaussian.(Ωdiff, 0, σ), SuperOperatorQuantumObject, dims) F2 = dense_to_sparse(F2, tol) # Ohmic reservoir