RESULT: position of the low frequency peaks

src/rabifun/plots.py

@@ -53,17 +53,26 @@ def plot_sidebands(ax, params: Params):
     :param params: system parameters
     """
     energy = params.rabi_splitting
-    sidebands = (
+
+    first_sidebands = np.abs(
+        -params.laser_detuning + np.array([1, -1]) * energy / 2 - params.Δ / 2
+    )
+    second_sidebands = (
         params.Ω - params.laser_detuning + np.array([1, -1]) * energy / 2 - params.Δ / 2
     )
 
     ax.axvline(params.Ω - params.Δ, color="black", label="steady state")
 
-    for n, sideband in enumerate(sidebands):
+    for n, sideband in enumerate(first_sidebands):
         ax.axvline(
             sideband,
-            color=f"C{n}",
-            label=f"rabi-sideband {n}",
+            color=f"C1",
+        )
+
+    for n, sideband in enumerate(second_sidebands):
+        ax.axvline(
+            sideband,
+            color=f"C2",
         )
 
     ax.legend()

src/rabifun/system.py

@@ -40,7 +40,7 @@ class Params:
 
     @property
     def rabi_splitting(self):
-        return 1 / 2 * np.sqrt(4 * self.d**2 + 4 * self.Δ**2)
+        return np.sqrt(self.d**2 + self.Δ**2)
 
 
 class RuntimeParams: