From 44667801f4e387107496d90f7035260861a7b76c Mon Sep 17 00:00:00 2001 From: Valentin Boettcher Date: Sun, 12 Jun 2022 11:03:37 +0200 Subject: [PATCH] better spectral separation --- .../10_antizeno_engine/10_first_anti_zeno.py | 2 +- .../10_antizeno_engine/anti_zeno_engine.org | 28 +++++++++---------- 2 files changed, 15 insertions(+), 15 deletions(-) diff --git a/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py b/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py index 521d7bc..12a45ba 100644 --- a/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py +++ b/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py @@ -82,7 +82,7 @@ def anti_zeno_engine( return model, (τ_mod, τ_c, τ_bath, cycles, model.ω_s) -model, params = anti_zeno_engine() +model, params = anti_zeno_engine(ω_0=6) params diff --git a/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org b/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org index 799aaaf..be2c09b 100644 --- a/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org +++ b/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org @@ -26,7 +26,7 @@ Init ray and silence stocproc. #+end_src #+RESULTS: -: RayContext(dashboard_url='', python_version='3.9.12', ray_version='1.12.1', ray_commit='4863e33856b54ccf8add5cbe75e41558850a1b75', address_info={'node_ip_address': '192.168.100.160', 'raylet_ip_address': '192.168.100.160', 'redis_address': None, 'object_store_address': '/tmp/ray/session_2022-06-11_16-10-06_802932_36816/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-06-11_16-10-06_802932_36816/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-06-11_16-10-06_802932_36816', 'metrics_export_port': 62760, 'gcs_address': '192.168.100.160:42264', 'address': '192.168.100.160:42264', 'node_id': 'a5d1451f14fd9269f1c7c703bd2e23cb84a94170d4ad40b3e7861a51'}) +: RayContext(dashboard_url='', python_version='3.9.12', ray_version='1.12.1', ray_commit='4863e33856b54ccf8add5cbe75e41558850a1b75', address_info={'node_ip_address': '192.168.99.108', 'raylet_ip_address': '192.168.99.108', 'redis_address': None, 'object_store_address': '/tmp/ray/session_2022-06-12_10-55-13_766289_1210189/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-06-12_10-55-13_766289_1210189/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-06-12_10-55-13_766289_1210189', 'metrics_export_port': 43375, 'gcs_address': '192.168.99.108:58576', 'address': '192.168.99.108:58576', 'node_id': 'b177960ddfd37df792a8dacbd7f82440431287248aef27d90193e9fd'}) #+begin_src jupyter-python :results none from hops.util.logging_setup import logging_setup @@ -103,12 +103,12 @@ Init ray and silence stocproc. * Model Definition #+begin_src jupyter-python - model, params = anti_zeno_engine() + model, params = anti_zeno_engine(ω_0=6) params #+end_src #+RESULTS: -| 1.2566370614359172 | 10.053096491487338 | 10.466663926059324 | 8 | (-4 6) | +| 1.2566370614359172 | 10.053096491487338 | 10.466663926059324 | 8 | (0 10) | Let's test the assumptions of the paper. @@ -122,14 +122,14 @@ Let's test the assumptions of the paper. return model.bcf(1)(t) + 2 * (model.thermal_correlations(1)(t).real) plt.plot(model.t, np.abs(thermal_bcf(model.t))/np.abs(thermal_bcf(0))) plt.plot(model.t, model.L[0].operator_norm(model.t)) - plt.plot(model.t, model.H.operator_norm(model.t) - ω_0/2) + #plt.plot(model.t, model.H.operator_norm(model.t) - ω_0/2) plt.plot(model.t, np.exp(- model.t * np.min(np.array(model.bcf_coefficients()[1][0]).real))) #+end_src #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/495a08e3d1d226b8d6a5c7034c56373089b44e01.svg]] +| | +[[file:./.ob-jupyter/6627b6e160d82b9c88aede960db16d75efcbe787.svg]] :END: #+begin_src jupyter-python :tangle nil @@ -140,8 +140,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/c043366e4686530dd5fa8bcddaaf1bbf6ebab95f.svg]] +| | +[[file:./.ob-jupyter/a1a3bf487639295688b6e279e56b9037bc41f08b.svg]] :END: @@ -168,8 +168,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -: -[[file:./.ob-jupyter/bb63ba7429e684f913a8b26a2a9540083c0e30cc.svg]] +: +[[file:./.ob-jupyter/02f5910868ac6afe3e2039db428ac5cccc10491b.svg]] :END: - **too fast decoupling kills it** @@ -200,7 +200,7 @@ Let's test the assumptions of the paper. : pts[1:N+1, 1] = dep1slice : /nix/store/pwhaggpgvxhy410r6vlx448v9lp8n08s-python3-3.9.12-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5224: ComplexWarning: Casting complex values to real discards the imaginary part : pts[N+2:, 1] = dep2slice[::-1] -[[file:./.ob-jupyter/e44159b8feb819c05aa4a795c9acce3ad23cbe67.svg]] +[[file:./.ob-jupyter/45ac729e831799ca145ea562480374c5efea6b5f.svg]] :END: @@ -216,8 +216,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/f3809dd8b346998aea1c4fad127f724b865a932f.svg]] +| | +[[file:./.ob-jupyter/9a0d753dbfd4743e0a0abd43d543885985439bd5.svg]] :END: ** TODO Power and Efficiency @@ -238,7 +238,7 @@ We need the time points where we sample the total energy. #+end_src #+RESULTS: -[[file:./.ob-jupyter/c0bf0862cc10cdbbd0915588568b3cb4344efc07.svg]] +[[file:./.ob-jupyter/d727e5dfb8edb23183c44c78f45d21d3adaf1afd.svg]] #+begin_src jupyter-python with aux.get_data(model) as data: