Correct s = 0 case

shareloqs · Jun 19, 2024 · 0885a47 · 0885a47
1 parent 2d44f7a
commit 0885a47
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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)]