From 0885a47d216af30348fa0ffcd6cffa9805a7b904 Mon Sep 17 00:00:00 2001
From: Thibaut Lacroix <57836508+tfmlaX@users.noreply.github.com>
Date: Wed, 19 Jun 2024 19:52:24 +0200
Subject: [PATCH] Correct s = 0 case

---
 src/models.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

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)]