working on timing comparisons

EQuS · May 20, 2024 · f6e1e9b · f6e1e9b
1 parent ce67f8c
commit f6e1e9b
Show file tree

Hide file tree

Showing 4 changed files with 537 additions and 9 deletions.
diff --git a/qcsys/devices/resonator.py b/qcsys/devices/resonator.py
@@ -23,12 +23,18 @@ def common_ops(self):
         ops["a"] = jqt.destroy(N)
         ops["a_dag"] = jqt.create(N)
         ops["phi"] = self.phi_zpf()*(ops["a"] + ops["a_dag"])  
+        ops["n"] = 1j * self.n_zpf() * (ops["a_dag"] - ops["a"])
+
         return ops
 
     def phi_zpf(self):
         """Return Phase ZPF."""
         return (2*self.params["Ec"]/self.params["El"])**(.25)
 
+    def n_zpf(self):
+        n_zpf = (self.params["El"] / (32.0 * self.params["Ec"])) ** (0.25)
+        return n_zpf
+
     def get_linear_ω(self):
         """Get frequency of linear terms."""
         return jnp.sqrt(8*self.params["El"]*self.params["Ec"])

diff --git a/tutorials/1-demo.ipynb b/tutorials/1-demo.ipynb
diff --git a/tutorials/3-fluxonium.ipynb b/tutorials/3-fluxonium.ipynb
@@ -77,8 +77,8 @@
     "axs = qubit.plot_wavefunctions(phi_vals, max_n=6, mode=\"real\")\n",
     "axs.set_ylim(-6,11)\n",
     "axs.set_xlim(-2,2)\n",
-    "axs.set_xlabel(\"$\\\\theta/\\Phi_0$\")\n",
-    "axs.set_title(\"$\\\\theta_{\\\\rm ext} = 0.5\\Phi_0$\")"
+    "axs.set_xlabel(\"$\\\\theta/\\\\Phi_0$\")\n",
+    "axs.set_title(\"$\\\\theta_{\\\\rm ext} = 0.5\\\\Phi_0$\")"
    ]
   },
   {
@@ -394,7 +394,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.9.18"
+   "version": "3.12.2"
   }
  },
  "nbformat": 4,

diff --git a/tutorials/5-time-comparisons.ipynb b/tutorials/5-time-comparisons.ipynb