fix: apply review suggestions from PR #138

Author: NiallCarbery

tutorials-v4/lectures/Lecture-13-Resonance-flourescence.md

@@ -105,11 +105,11 @@ axes[0].set_title("Bloch Vector Components vs Time")
 
 axes[1].plot(result.times, result.expect[5], "b", label=r"$P_e$")
 
-axes[1].set_ylabel(r"$P_e$", fontsize=16)
+axes[1].set_ylabel(r"$P_g$", fontsize=16)
 axes[1].set_xlabel("time", fontsize=16)
 axes[1].legend()
 axes[1].set_ylim(0, 1)
-axes[1].set_title("Excited State Population vs Time");
+axes[1].set_title("Ground State Population vs Time");
 ```
 
 ```python
@@ -120,13 +120,15 @@ for idx, gamma0 in enumerate([0.1 * Omega, 0.5 * Omega, 1.0 * Omega]):
     HL, c_ops = system_spec(Omega, gamma0, N)
     result = mesolve(HL, psi0, tlist, c_ops, e_ops)
 
-    ax.plot(result.times, result.expect[5], "b",
-            label=fr"$P_e$ ($\gamma_0={gamma0:.2f}$)")
+    ax.plot(
+        result.times, result.expect[5], "b",
+        label=fr"$P_g$ ($\gamma_0={gamma0:.2f}$)"
+    )
 
 ax.set_ylim(0, 1)
 ax.set_xlabel("time", fontsize=16)
-ax.set_ylabel(r"$P_e$", fontsize=16)
-ax.set_title("Excited State Population for Different Decay Rates")
+ax.set_ylabel(r"$P_g$", fontsize=16)
+ax.set_title("Ground State Population for Different Decay Rates")
 ax.legend();
 ```
 

tutorials-v4/lectures/Lecture-3A-Dicke-model.md

@@ -241,7 +241,7 @@ fig, axes = plt.subplots(1, 1, figsize=(12, 6))
 
 for NN in N_vec:
 
-    entropy_cavity, entropy_spin = calculcate_entropy(MM, NN, g_vec)
+    entropy_cavity, entropy_spin = calculate_entropy(MM, NN, g_vec)
 
     axes.plot(g_vec, entropy_cavity, "b", label="N = %d" % NN)
     axes.plot(g_vec, entropy_spin, "r--")

tutorials-v5/lectures/Lecture-13-Resonance-flourescence.md

@@ -105,11 +105,11 @@ axes[0].set_title("Bloch Vector Components vs Time")
 
 axes[1].plot(result.times, result.expect[5], "b", label=r"$P_e$")
 
-axes[1].set_ylabel(r"$P_e$", fontsize=16)
+axes[1].set_ylabel(r"$P_g$", fontsize=16)
 axes[1].set_xlabel("time", fontsize=16)
 axes[1].legend()
 axes[1].set_ylim(0, 1)
-axes[1].set_title("Excited State Population vs Time");
+axes[1].set_title("Ground State Population vs Time");
 ```
 
 ```python
@@ -120,13 +120,15 @@ for idx, gamma0 in enumerate([0.1 * Omega, 0.5 * Omega, 1.0 * Omega]):
     HL, c_ops = system_spec(Omega, gamma0, N)
     result = mesolve(HL, psi0, tlist, c_ops, e_ops=e_ops)
 
-    ax.plot(result.times, result.expect[5], "b",
-            label=fr"$P_e$ ($\gamma_0={gamma0:.2f}$)")
+        ax.plot(
+            result.times, result.expect[5], "b",
+            label=fr"$P_g$ ($\gamma_0={gamma0:.2f}$)"
+        )
 
 ax.set_ylim(0, 1)
 ax.set_xlabel("time", fontsize=16)
-ax.set_ylabel(r"$P_e$", fontsize=16)
-ax.set_title("Excited State Population for Different Decay Rates")
+ax.set_ylabel(r"$P_g$", fontsize=16)
+ax.set_title("Ground State Population for Different Decay Rates")
 ax.legend();
 ```
 

tutorials-v5/lectures/Lecture-2B-Single-Atom-Lasing.md

@@ -130,7 +130,7 @@ if rate > 0.0:
 Here we evolve the system with the Lindblad master equation solver, and we request that the expectation values of the operators $a^\dagger a$ and $\sigma_+\sigma_-$ are returned by the solver by passing the list `[a.dag()*a, sm.dag()*sm]` as the fifth argument to the solver.
 
 ```python
-opt = {'nsteps': 2000} # allow extra time-steps
+opt = {'nsteps': 2000}  # allow extra time-steps
 output = mesolve(H, psi0, tlist, c_ops, e_ops=[a.dag() * a, sm.dag() * sm],
                  options=opt)
 ```