added g_cd and g_3

EQuS · Jul 25, 2024 · dc795af · dc795af
1 parent ff326e6
commit dc795af
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 7 deletions.
diff --git a/tutorials/ats-coupler/0-resonator-ats-resonator.ipynb b/tutorials/ats-coupler/0-resonator-ats-resonator.ipynb
@@ -39,7 +39,9 @@
     {
      "data": {
       "text/plain": [
-       "{'g_ex': Array(-0.00212182+0.j, dtype=complex128),\n",
+       "{'g_3': Array(1.28654004e-05+0.j, dtype=complex128),\n",
+       " 'g_cd': Array(2.57308009e-05+0.j, dtype=complex128),\n",
+       " 'g_ex': Array(-0.00212182+0.j, dtype=complex128),\n",
        " 'ω_ATS': Array(6.02181564, dtype=float64),\n",
        " 'ω_ResonatorA': Array(4.99584201, dtype=float64),\n",
        " 'ω_ResonatorB': Array(7.00151547, dtype=float64)}"
@@ -80,17 +82,17 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "text/plain": [
-       "Array([-2.14061657e-01, -2.14015472e-01, -2.14154028e-01, -2.11546336e-01,\n",
-       "       -4.03521234e+01,  1.18006174e+05], dtype=float64)"
+       "Array([ 9.07343001e+00,  9.07285802e+00,  9.07207109e+00,  9.07084541e+00,\n",
+       "       -1.22035928e+01,  4.27481506e+05], dtype=float64)"
       ]
      },
-     "execution_count": 7,
+     "execution_count": 6,
      "metadata": {},
      "output_type": "execute_result"
     }

diff --git a/tutorials/ats-coupler/rar_helpers.py b/tutorials/ats-coupler/rar_helpers.py
@@ -55,8 +55,8 @@ def create_devices_linear_rar(params):
     Ej_c = params.get("ATS__E_J", 37.0)
     dEj_c = params.get("ATS__dE_J", 0.0)
     Ej2_c = params.get("ATS__E_J_2", 0.0)
-    phi_sum_c = params.get("ATS__phi_sum", 0.25)
-    phi_delta_c = params.get("ATS__phi_sum", 0.0)
+    phi_sum_c = params.get("ATS__phi_sum_ext", 0.25)
+    phi_delta_c = params.get("ATS__phi_delta_ext", 0.25)
 
     ats = qs.ATS.create(
         N_CONS["ats"]["bare"],
@@ -143,6 +143,8 @@ def get_metrics_linear_rar(params):
 
     drive_strength = params.get("ATS__drive_strength", DRIVE_STRENGTH)
     metrics[f"g_ex"] = ϕ["ATS"]["ResonatorA"] * ϕ["ATS"]["ResonatorB"] * drive_strength * ats.params["Ej"]
+    metrics[f"g_3"] = (1/2) * ϕ["ATS"]["ResonatorA"]**2 * ϕ["ATS"]["ResonatorB"] * drive_strength * ats.params["Ej"]
+    metrics[f"g_cd"] = ϕ["ATS"]["ResonatorA"]**2 * ϕ["ATS"]["ResonatorB"] * drive_strength * ats.params["Ej"]
     return ϕ, metrics, system