diff --git a/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno_cc.py b/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno_cc.py index ffd8e34..9fd6740 100644 --- a/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno_cc.py +++ b/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) diff --git a/python/energy_flow_proper/10_antizeno_engine/continously_couple.org b/python/energy_flow_proper/10_antizeno_engine/continously_couple.org index 480a66b..621d086 100644 --- a/python/energy_flow_proper/10_antizeno_engine/continously_couple.org +++ b/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: -|
| | -[[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: |
| | -[[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: +|
| | +[[file:./.ob-jupyter/afa180fb667b63a127dbecc86f43980e65eabb6a.svg]] +:END: #+begin_src jupyter-python @@ -129,17 +126,53 @@ Init ray and silence stocproc. #+RESULTS: :RESULTS: |
| | -[[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: -|
| | -[[file:./.ob-jupyter/e01a22e101c610a13ff055f6fb32c0ec6cf57c7c.svg]] +# [goto error] +#+begin_example + --------------------------------------------------------------------------- + RuntimeError Traceback (most recent call last) + Input In [57], in () + ----> 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..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: |
| | -[[file:./.ob-jupyter/92ad75e633d380edb09712b501789c8151445708.svg]] +[[file:./.ob-jupyter/39b065c6a65911b3349c41e28138c9e15ff080da.svg]] :END: @@ -159,6 +192,6 @@ Init ray and silence stocproc. #+RESULTS: :RESULTS: -|
| | EnsembleValue | (((10000 -1.238968933350463e-05 0.00034012433086271946))) | -[[file:./.ob-jupyter/40029fe276ad30c0cf7704976bf9f6ec980ad25a.svg]] +|
| | EnsembleValue | (((40000 -0.000918483945087149 0.00014770367528033136))) | +[[file:./.ob-jupyter/833c93cdbf0589c458fc058a3db8852a1144950e.svg]] :END: diff --git a/python/energy_flow_proper/10_antizeno_engine/flake.lock b/python/energy_flow_proper/10_antizeno_engine/flake.lock index 56b0cb6..e68d171 100644 --- a/python/energy_flow_proper/10_antizeno_engine/flake.lock +++ b/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": {