diff --git a/hiro_models/one_qubit_model.py b/hiro_models/one_qubit_model.py
index 7adf92a..c870b14 100644
--- a/hiro_models/one_qubit_model.py
+++ b/hiro_models/one_qubit_model.py
@@ -583,7 +583,7 @@ class QubitModelMutliBath(Model):
 
         return [
             bcf_scale(δ, L, self.t.max(), s, ω)
-            for δ, L, s, ω in zip(self.δ, self.L, self.s, self.ω_c)
+            for δ, L, s, ω in zip(self.δ, self.coupling_operators, self.s, self.ω_c)
         ]
 
     def bcf(