try to eek out the "schwebungen"

python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno_cc.py

@@ -28,37 +28,46 @@ from anti_zeno_engine import *
     params
 ) = anti_zeno_engine(
     Δ=11,
-    ε=.1,
+    ε=1,#.1,
     ω_c=1,
-    ε_couple=0.7,
-    n=10,
+    ε_couple=.68,
+    n=40,
     detune=.5,
     ω_0=20,
-    T_c=1e3,
-    T_h=1e4,
-    δ=[3.2*.01/4, 1*.01/4],
+    T_c=5e3,
+    T_h=5e4,
+    δ=[3.2*.01/10, 1.08*.01/10],
     γ=.2,
     switch_cycles=1,
     therm_initial_state=False,
-    ε_init=.00001/2,
+    ε_init=.000001/2,
     terms=7,
-    dt=0.01,
-    sp_tol=1e-3,
+    dt=0.01/2,
+    sp_tol=1e-4,
     init_time_steps=10,
+    interpolation_multiplier=10,
 )
 model.k_max = 3
 
-aux.integrate(model, 1)
+ts = model.t # np.linspace(0,10,1000)
+proc = model.thermal_process(1)
+import hops
+z=hops.core.utility.uni_to_gauss(np.random.rand(proc.get_num_y() * 2))
+proc.new_process(z)
+pu.plot_complex(ts, proc(ts) * model.bcf_scales[1])
 
 params.cycles
 
+aux.integrate(model, 10)
+
 plot_az_coupling_diagram(model, params)
 
 plot_az_sd_overview(model, params)
-fs.export_fig("with_gap_coupling_diagram", tikz=False, y_scaling=.4)
 
 plot_total_power(model, params)
 
 plot_excited_state(model, params)
 
-plot_power_output(model, params)
+plot_coherences(model, params)
+
+plot_power_output(model, params, steady_index=0)

python/energy_flow_proper/10_antizeno_engine/continously_couple.org

@@ -48,13 +48,13 @@ Init ray and silence stocproc.
       Δ=11,
       ε=1,#.1,
       ω_c=1,
-      ε_couple=.0001,
-      n=1,
+      ε_couple=.68,
+      n=40,
       detune=.5,
       ω_0=20,
-      T_c=1e3,
-      T_h=1e4,
-      δ=[3.2*.01/4, 1.08*.01/4],
+      T_c=5e3,
+      T_h=5e4,
+      δ=[3.2*.01/10, 1.08*.01/10],
       γ=.2,
       switch_cycles=1,
       therm_initial_state=False,
@@ -79,19 +79,13 @@ Init ray and silence stocproc.
   pu.plot_complex(ts, proc(ts) * model.bcf_scales[1])
 #+end_src
 
-#+RESULTS:
-:RESULTS:
-| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:> |
-[[file:./.ob-jupyter/9e7c77f1ba1686b8d5c068ea4463029eaf871a52.svg]]
-:END:
-
 
 #+begin_src jupyter-python
   params.cycles
 #+end_src
 
 #+RESULTS:
-: 122
+: 10
 
 #+begin_src jupyter-python
   aux.integrate(model, 10)
@@ -111,7 +105,7 @@ Init ray and silence stocproc.
 #+RESULTS:
 :RESULTS:
 | <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:xlabel= | $\tau$ | > |
-[[file:./.ob-jupyter/b21178f0e184aef19d53f25ad13a0514a71f1356.svg]]
+[[file:./.ob-jupyter/2a8f8502ac9e11101c59c75c15da03cab48dbc86.svg]]
 :END:
 
 #+begin_src jupyter-python
@@ -119,7 +113,10 @@ Init ray and silence stocproc.
 #+end_src
 
 #+RESULTS:
-: 7ef94938-dccc-41cf-a788-0729c098f789
+:RESULTS:
+| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:xlabel= | $\omega$ | ylabel= | Spectral Density | > |
+[[file:./.ob-jupyter/afa180fb667b63a127dbecc86f43980e65eabb6a.svg]]
+:END:
 
 
 #+begin_src jupyter-python
@@ -129,17 +126,53 @@ Init ray and silence stocproc.
 #+RESULTS:
 :RESULTS:
 | <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:> |
-[[file:./.ob-jupyter/f26abd3f142832f7b78b92c34f41832fad9462ed.svg]]
+[[file:./.ob-jupyter/a2ee618f4377def0d3f99b20f4d3fc914cbb1198.svg]]
 :END:
 
 #+begin_src jupyter-python
-  plot_excited_state(model, params, with_init=True)
+  plot_excited_state(model, params)
 #+end_src
 
 #+RESULTS:
 :RESULTS:
-| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:> |
-[[file:./.ob-jupyter/e01a22e101c610a13ff055f6fb32c0ec6cf57c7c.svg]]
+# [goto error]
+#+begin_example
+  ---------------------------------------------------------------------------
+  RuntimeError                              Traceback (most recent call last)
+  Input In [57], in <cell line: 1>()
+  ----> 1 plot_excited_state(model, params)
+
+  File ~/Documents/Projects/UNI/master/masterarb/python/energy_flow_proper/10_antizeno_engine/plot_utils.py:40, in wrap_plot.<locals>.wrapped(ax, setup_function, *args, **kwargs)
+       37 if not ax:
+       38     fig, ax = setup_function()
+  ---> 40 ret_val = f(*args, ax=ax, **kwargs)
+       41 return (fig, ax, ret_val) if ret_val else (fig, ax)
+
+  File ~/Documents/Projects/UNI/master/masterarb/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.py:285, in plot_excited_state(model, params, ax, with_init)
+      283 @wrap_plot
+      284 def plot_excited_state(model, params, ax=None, with_init=False):
+  --> 285     with aux.get_data(model) as data:
+      286         mask = model.t >= 0 if with_init else model.t > params.τ_init
+      288         plot_with_σ(
+      289             model.t[mask],
+      290             abs(
+     (...)
+      298             ax=ax,
+      299         )
+
+  File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:146, in get_data(model, data_path, read_only, **kwargs)
+      135         return HIData(
+      136             path,
+      137             hi_key=model.hops_config,
+     (...)
+      142             **kwargs,
+      143         )
+      145 else:
+  --> 146     raise RuntimeError(f"No data found for model with hash '{hexhash}'.")
+
+  RuntimeError: No data found for model with hash '043441e996b658af99b232ace67ecc5c8d9ff76b8190a6707bf2556908770a5c'.
+#+end_example
+[[file:./.ob-jupyter/6140cc8a5adedec7ec2beecd22a919c202fd1f03.svg]]
 :END:
 
 #+begin_src jupyter-python
@@ -149,7 +182,7 @@ Init ray and silence stocproc.
 #+RESULTS:
 :RESULTS:
 | <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:> |
-[[file:./.ob-jupyter/92ad75e633d380edb09712b501789c8151445708.svg]]
+[[file:./.ob-jupyter/39b065c6a65911b3349c41e28138c9e15ff080da.svg]]
 :END:
 
 
@@ -159,6 +192,6 @@ Init ray and silence stocproc.
 
 #+RESULTS:
 :RESULTS:
-| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:title= | (center : \(\bar{P}=-0.00001\pm 0.00034\)) | xlabel= | $\tau$ | ylabel= | $\langle{P}\rangle$ | > | EnsembleValue | (((10000 -1.238968933350463e-05 0.00034012433086271946))) |
-[[file:./.ob-jupyter/40029fe276ad30c0cf7704976bf9f6ec980ad25a.svg]]
+| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:title= | (center : \(\bar{P}=-0.00092\pm 0.00015\)) | xlabel= | $\tau$ | ylabel= | $\langle{P}\rangle$ | > | EnsembleValue | (((40000 -0.000918483945087149 0.00014770367528033136))) |
+[[file:./.ob-jupyter/833c93cdbf0589c458fc058a3db8852a1144950e.svg]]
 :END:

python/energy_flow_proper/10_antizeno_engine/flake.lock

@@ -32,11 +32,11 @@
     },
     "nixpkgs": {
       "locked": {
-        "lastModified": 1660396586,
-        "narHash": "sha256-ePuWn7z/J5p2lO7YokOG1o01M0pDDVL3VrStaPpS5Ig=",
+        "lastModified": 1660485612,
+        "narHash": "sha256-sSLW1KaB1adKTJn9+Ja3h3AaS7QCZyhUKiSUStcLg80=",
         "owner": "NixOS",
         "repo": "nixpkgs",
-        "rev": "e105167e98817ba9fe079c6c3c544c6ef188e276",
+        "rev": "6512b21eabb4d52e87ea2edcf31a288e67b2e4f8",
         "type": "github"
       },
       "original": {