From c7ebe963fc83f7396722a1a7e7943815299e602a Mon Sep 17 00:00:00 2001
From: Valentin Boettcher <hiro@protagon.space>
Date: Mon, 27 May 2024 14:21:59 -0400
Subject: [PATCH] rename d->g_0

---
 scripts/experiments/002_rabi_detuning_scan.py | 10 +++++-----
 src/rabifun/system.py                         |  6 +++---
 2 files changed, 8 insertions(+), 8 deletions(-)

diff --git a/scripts/experiments/002_rabi_detuning_scan.py b/scripts/experiments/002_rabi_detuning_scan.py
index f7e2f51..c384b89 100644
--- a/scripts/experiments/002_rabi_detuning_scan.py
+++ b/scripts/experiments/002_rabi_detuning_scan.py
@@ -13,7 +13,7 @@ def transient_rabi():
         η=0.7,
         Ω=13,
         δ=1 / 4,
-        d=0.02,
+        g_0=0.02,
         laser_detuning=0,
         Δ=0,
         N=3,
@@ -40,7 +40,7 @@ def steady_rabi():
         η=0.7,
         Ω=13,
         δ=1 / 4,
-        d=0.02,
+        g_0=0.02,
         laser_detuning=0,
         Δ=0,
         N=3,
@@ -67,7 +67,7 @@ def noisy_transient_rabi():
         η=0.7,
         Ω=13,
         δ=1 / 4,
-        d=0.02,
+        g_0=0.02,
         laser_detuning=0,
         Δ=0,
         N=3,
@@ -102,7 +102,7 @@ def ringdown_after_rabi():
         η=0.7,
         Ω=13,
         δ=1 / 4,
-        d=0.02,
+        g_0=0.02,
         laser_detuning=laser_detuning,
         Δ=0,
         N=3,
@@ -137,7 +137,7 @@ def sweep():
         η=1,
         Ω=1,
         δ=1 / 4,
-        d=0.0,
+        g_0=0.0,
         laser_detuning=-2,
         Δ=0,
         N=3,
diff --git a/src/rabifun/system.py b/src/rabifun/system.py
index 6df22bf..61e4fae 100644
--- a/src/rabifun/system.py
+++ b/src/rabifun/system.py
@@ -21,7 +21,7 @@ class Params:
     """Decay rate :math:`\eta/2` of the system in frequency units (no
     :math:`2 \pi`)."""
 
-    d: float = 0.01
+    g_0: float = 0.01
     """Drive amplitude in units of :any:`Ω`."""
 
     Δ: float = 0.0
@@ -65,7 +65,7 @@ class Params:
     @property
     def rabi_splitting(self):
         """The Rabi splitting of the system in *frequency units*."""
-        return np.sqrt((self.Ω * self.d) ** 2 + self.Δ**2)
+        return np.sqrt((self.Ω * self.g_0) ** 2 + self.Δ**2)
 
     @property
     def ω_eom(self):
@@ -141,7 +141,7 @@ def make_righthand_side(runtime_params: RuntimeParams, params: Params):
 
         if (params.drive_off_time is None) or (t < params.drive_off_time):
             differential += eom_drive(
-                t, x, 2 * np.pi * params.Ω * params.d, params.ω_eom
+                t, x, 2 * np.pi * params.Ω * params.g_0, params.ω_eom
             )
 
         if (params.laser_off_time is None) or (t < params.laser_off_time):