diff --git a/src/models.jl b/src/models.jl index 329f792..d1fe7a9 100644 --- a/src/models.jl +++ b/src/models.jl @@ -511,8 +511,8 @@ function chaincoeffs_ohmic(nummodes, α, s; ωc=1, soft=false) else if s==0 c0 = sqrt(2α)*ωc - e = fill(0,nummodes) - t = [ωc*(n+1)/(2n+1) for n in 0:(nummodes-2)] + e = fill(0.5*ωc,nummodes) + t = [ωc*(n+1)/(4n+2) for n in 0:(nummodes-2)] else c0 = sqrt(2α/(s+1))*ωc e = [(ωc/2)*(1 + (s^2)/((s+2n)*(2+s+2n))) for n in 0:(nummodes-1)]