From d1fae818eb859c427285a3475ba0e36715563229 Mon Sep 17 00:00:00 2001
From: Valentin Boettcher <hiro@protagon.space>
Date: Thu, 24 Feb 2022 17:11:32 +0100
Subject: [PATCH] add a phenomentological bcf coupling strength estimator

---
 two_qubit_model.py | 13 +++++++++++++
 1 file changed, 13 insertions(+)

diff --git a/two_qubit_model.py b/two_qubit_model.py
index 63cfc08..f47e9c4 100644
--- a/two_qubit_model.py
+++ b/two_qubit_model.py
@@ -280,6 +280,19 @@ class TwoQubitModel:
 
         return float(self.δ[i]) ** 2
 
+    def η(self, i: int):
+        """The BCF pre-factor :math:`η` of the ``i``th bath."""
+        ω_c = float(self.ω_c[i])
+        s = float(self.s[i])
+        T = float(self.T[i])
+
+        return (
+            1
+            / (ω_c * s) ** s
+            * np.exp(s)
+            * (max([1, np.exp(ω_c * s * 1 / T) - 1]) if T > 0 else 1)
+        )
+
     def bcf(self, i: int) -> hops.util.bcf.OhmicBCF_zeroTemp:
         """
         The normalized zero temperature BCF of the  ``i``th bath.