diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/09f1f4b3ba44ce70adf73cf8d463082f17ccd1bd38f093bbf70c823832f6cf69_delta_0.01.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/09f1f4b3ba44ce70adf73cf8d463082f17ccd1bd38f093bbf70c823832f6cf69_delta_0.01.npy new file mode 100644 index 0000000..d69a94c Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/09f1f4b3ba44ce70adf73cf8d463082f17ccd1bd38f093bbf70c823832f6cf69_delta_0.01.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/16ef657884259c28c4741600f5ebc514797c16a5222565a82d3cb9fd3a6257f2_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/16ef657884259c28c4741600f5ebc514797c16a5222565a82d3cb9fd3a6257f2_delta_0.4.npy new file mode 100644 index 0000000..c24baed Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/16ef657884259c28c4741600f5ebc514797c16a5222565a82d3cb9fd3a6257f2_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/1d71370526bda3467a28da7b26161957bbaf96c0557d0c2d56da8580d028ba2f_delta_0.01.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/1d71370526bda3467a28da7b26161957bbaf96c0557d0c2d56da8580d028ba2f_delta_0.01.npy new file mode 100644 index 0000000..d69a94c Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/1d71370526bda3467a28da7b26161957bbaf96c0557d0c2d56da8580d028ba2f_delta_0.01.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/2226671cc747bb3e78687eade4ab059069272522c26d57cfd550f3a5303812a0_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/2226671cc747bb3e78687eade4ab059069272522c26d57cfd550f3a5303812a0_delta_0.4.npy new file mode 100644 index 0000000..1b91296 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/2226671cc747bb3e78687eade4ab059069272522c26d57cfd550f3a5303812a0_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta.npy new file mode 100644 index 0000000..4c110dd Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.2.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.2.npy new file mode 100644 index 0000000..669e9ac Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.2.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.3.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.3.npy new file mode 100644 index 0000000..5b6730f Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.3.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.4.npy new file mode 100644 index 0000000..6716907 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/22ea0b849504352255ee6476865d2487742aa5f7caa809a2c1f8d51229d7e1df_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/26c5115438493c72c3318beebf5a22187645914d8099e24dbc5b96c3693275e3_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/26c5115438493c72c3318beebf5a22187645914d8099e24dbc5b96c3693275e3_delta_0.4.npy new file mode 100644 index 0000000..cc18629 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/26c5115438493c72c3318beebf5a22187645914d8099e24dbc5b96c3693275e3_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/3144749732eb55f9020aa4854a56866c13a2a00cd3145cea9693c9c0b9ddb3c9_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/3144749732eb55f9020aa4854a56866c13a2a00cd3145cea9693c9c0b9ddb3c9_delta_0.4.npy new file mode 100644 index 0000000..4093752 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/3144749732eb55f9020aa4854a56866c13a2a00cd3145cea9693c9c0b9ddb3c9_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/358d660806be6914721a1d7be0eda1eaaf3d9d60bce24199f597675195fcbfd0_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/358d660806be6914721a1d7be0eda1eaaf3d9d60bce24199f597675195fcbfd0_delta_0.4.npy new file mode 100644 index 0000000..d94bf23 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/358d660806be6914721a1d7be0eda1eaaf3d9d60bce24199f597675195fcbfd0_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta.npy new file mode 100644 index 0000000..0a8d360 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.2.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.2.npy new file mode 100644 index 0000000..97e12d0 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.2.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.3.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.3.npy new file mode 100644 index 0000000..3933957 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.3.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.4.npy new file mode 100644 index 0000000..1b91296 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/36f4897232ee2a3b8ae25569affeec04ba36ed785fef5ed20641a37dee1441cd_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/37ded4eeccde59ae0ff9509c0b328adddfb6cdb9943271db395a7eace8389e7b_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/37ded4eeccde59ae0ff9509c0b328adddfb6cdb9943271db395a7eace8389e7b_delta_0.4.npy new file mode 100644 index 0000000..eb4c298 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/37ded4eeccde59ae0ff9509c0b328adddfb6cdb9943271db395a7eace8389e7b_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/3f1398a05ca4ee20bf3a5903516d9571b82fa2e1a98d4316a41ba44c59c6b916_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/3f1398a05ca4ee20bf3a5903516d9571b82fa2e1a98d4316a41ba44c59c6b916_delta_0.4.npy new file mode 100644 index 0000000..56799c0 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/3f1398a05ca4ee20bf3a5903516d9571b82fa2e1a98d4316a41ba44c59c6b916_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/46b4c7c1d0ded6f87d5999e2eac1a5effe002c20412dcdd5d55478149352ddc7_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/46b4c7c1d0ded6f87d5999e2eac1a5effe002c20412dcdd5d55478149352ddc7_delta_0.4.npy new file mode 100644 index 0000000..82653ae Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/46b4c7c1d0ded6f87d5999e2eac1a5effe002c20412dcdd5d55478149352ddc7_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/4d41549f37822da9c8f2313f27a0441116c3293a0ce078b87c10fc1724f38e36_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/4d41549f37822da9c8f2313f27a0441116c3293a0ce078b87c10fc1724f38e36_delta_0.4.npy new file mode 100644 index 0000000..1b91296 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/4d41549f37822da9c8f2313f27a0441116c3293a0ce078b87c10fc1724f38e36_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/51910290a72cc7f1fac3b9dda01893fdf777e1844edb8dee9d6b9763830ba938_delta_0.01.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/51910290a72cc7f1fac3b9dda01893fdf777e1844edb8dee9d6b9763830ba938_delta_0.01.npy new file mode 100644 index 0000000..d69a94c Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/51910290a72cc7f1fac3b9dda01893fdf777e1844edb8dee9d6b9763830ba938_delta_0.01.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/55bbab5a873e54423860d605291ded4817365feb76848e7eed882bebbc82ba69_delta_0.01.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/55bbab5a873e54423860d605291ded4817365feb76848e7eed882bebbc82ba69_delta_0.01.npy new file mode 100644 index 0000000..d69a94c Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/55bbab5a873e54423860d605291ded4817365feb76848e7eed882bebbc82ba69_delta_0.01.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/56a2e183907b4eb0b4d2c675f0ee0cc75cd64686db97e5d0978b35785974be96_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/56a2e183907b4eb0b4d2c675f0ee0cc75cd64686db97e5d0978b35785974be96_delta_0.4.npy new file mode 100644 index 0000000..6106748 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/56a2e183907b4eb0b4d2c675f0ee0cc75cd64686db97e5d0978b35785974be96_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/5dbbdd1d17d21e6e845093f2e12c54b4ebea1c0d08171320309953a319491f1e_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/5dbbdd1d17d21e6e845093f2e12c54b4ebea1c0d08171320309953a319491f1e_delta_0.4.npy new file mode 100644 index 0000000..530585b Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/5dbbdd1d17d21e6e845093f2e12c54b4ebea1c0d08171320309953a319491f1e_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/6c5b85b5d369a6a360244a132f2b6e3f9f91ff96c42f0f95fc60b34856a2ec00_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/6c5b85b5d369a6a360244a132f2b6e3f9f91ff96c42f0f95fc60b34856a2ec00_delta_0.4.npy new file mode 100644 index 0000000..39cf2c4 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/6c5b85b5d369a6a360244a132f2b6e3f9f91ff96c42f0f95fc60b34856a2ec00_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/7633774a5d0db77358d3e55d16e8d88214e5c6c5ee8b8673b1737a272be4e0c6_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/7633774a5d0db77358d3e55d16e8d88214e5c6c5ee8b8673b1737a272be4e0c6_delta_0.4.npy new file mode 100644 index 0000000..76e7731 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/7633774a5d0db77358d3e55d16e8d88214e5c6c5ee8b8673b1737a272be4e0c6_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/80225e80ed5aca9197228284d01ffc80b0d8cbd880c27bda280014dded27697b_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/80225e80ed5aca9197228284d01ffc80b0d8cbd880c27bda280014dded27697b_delta_0.4.npy new file mode 100644 index 0000000..2040c4e Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/80225e80ed5aca9197228284d01ffc80b0d8cbd880c27bda280014dded27697b_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/88286708af4677e5be4adad62fb54b02d26ed883652303400b7d5c4fb22d0514_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/88286708af4677e5be4adad62fb54b02d26ed883652303400b7d5c4fb22d0514_delta_0.4.npy new file mode 100644 index 0000000..cbc2014 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/88286708af4677e5be4adad62fb54b02d26ed883652303400b7d5c4fb22d0514_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/9496f7e236720b58aba98fe9ae688be798115479152d57a8d810bdac3304de8d_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/9496f7e236720b58aba98fe9ae688be798115479152d57a8d810bdac3304de8d_delta_0.4.npy new file mode 100644 index 0000000..0000b98 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/9496f7e236720b58aba98fe9ae688be798115479152d57a8d810bdac3304de8d_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/9d068e13160029cd7b10fa3050b4c6e02701b95b012a83fef4dc37efbd8ab123_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/9d068e13160029cd7b10fa3050b4c6e02701b95b012a83fef4dc37efbd8ab123_delta_0.4.npy new file mode 100644 index 0000000..e3765ec Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/9d068e13160029cd7b10fa3050b4c6e02701b95b012a83fef4dc37efbd8ab123_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/a352ae8b1cc052ab6cd89e830f49bd8c26fc9f2ded8f2953b1a5abaa7ec7327f_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/a352ae8b1cc052ab6cd89e830f49bd8c26fc9f2ded8f2953b1a5abaa7ec7327f_delta_0.4.npy new file mode 100644 index 0000000..0e9c096 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/a352ae8b1cc052ab6cd89e830f49bd8c26fc9f2ded8f2953b1a5abaa7ec7327f_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/b1ba97bf213e4b4bd1cac4c26e4c8283b62065cbf8b98c0c9df8ac9f7c8e83c4_delta_0.01.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/b1ba97bf213e4b4bd1cac4c26e4c8283b62065cbf8b98c0c9df8ac9f7c8e83c4_delta_0.01.npy new file mode 100644 index 0000000..d69a94c Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/b1ba97bf213e4b4bd1cac4c26e4c8283b62065cbf8b98c0c9df8ac9f7c8e83c4_delta_0.01.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/c299f299e024e9b2f3f248029e1ecc909f57732a55b6b94d9d669da70773e459_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/c299f299e024e9b2f3f248029e1ecc909f57732a55b6b94d9d669da70773e459_delta_0.4.npy new file mode 100644 index 0000000..c9193f2 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/c299f299e024e9b2f3f248029e1ecc909f57732a55b6b94d9d669da70773e459_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta.npy new file mode 100644 index 0000000..13078dc Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.2.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.2.npy new file mode 100644 index 0000000..64d8ba0 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.2.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.3.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.3.npy new file mode 100644 index 0000000..f1718e6 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.3.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.4.npy new file mode 100644 index 0000000..ea4a567 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/cb8eb9b85dd52a182764d78ea032ffe6859e9e3b84741c1fd3486d84769b4054_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/dd470fcecc6d07bffa7d44f5917825f7432fc94aa5d49c4d44d2c901032fe7df_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/dd470fcecc6d07bffa7d44f5917825f7432fc94aa5d49c4d44d2c901032fe7df_delta_0.4.npy new file mode 100644 index 0000000..803e249 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/dd470fcecc6d07bffa7d44f5917825f7432fc94aa5d49c4d44d2c901032fe7df_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/e8126a0070e4e4ec086966c6c8d9e64a67927acef9f55b4cef26a237ef85927d_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/e8126a0070e4e4ec086966c6c8d9e64a67927acef9f55b4cef26a237ef85927d_delta_0.4.npy new file mode 100644 index 0000000..88687c1 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/e8126a0070e4e4ec086966c6c8d9e64a67927acef9f55b4cef26a237ef85927d_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/fb7425ce3546cabb6dfac91e9318bf936a2ab46bf4e4518e9b9517d482be8f7f_delta_0.4.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/fb7425ce3546cabb6dfac91e9318bf936a2ab46bf4e4518e9b9517d482be8f7f_delta_0.4.npy new file mode 100644 index 0000000..f114950 Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/fb7425ce3546cabb6dfac91e9318bf936a2ab46bf4e4518e9b9517d482be8f7f_delta_0.4.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/.cache/ff083fbe9f769a29d103d91a2c9847e6ab43e1c20c8d5e5eaae0689b4a6534ea_delta_0.01.npy b/python/energy_flow_proper/08_dynamic_one_bath/.cache/ff083fbe9f769a29d103d91a2c9847e6ab43e1c20c8d5e5eaae0689b4a6534ea_delta_0.01.npy new file mode 100644 index 0000000..d69a94c Binary files /dev/null and b/python/energy_flow_proper/08_dynamic_one_bath/.cache/ff083fbe9f769a29d103d91a2c9847e6ab43e1c20c8d5e5eaae0689b4a6534ea_delta_0.01.npy differ diff --git a/python/energy_flow_proper/08_dynamic_one_bath/cutoff_dep.py b/python/energy_flow_proper/08_dynamic_one_bath/cutoff_dep.py index 6da6872..a935571 100644 --- a/python/energy_flow_proper/08_dynamic_one_bath/cutoff_dep.py +++ b/python/energy_flow_proper/08_dynamic_one_bath/cutoff_dep.py @@ -15,6 +15,8 @@ from hops.util.logging_setup import logging_setup import logging logging_setup(logging.INFO, show_stocproc=False) +print(es.models_table(ω_models)) + plt.plot(ωs, [model.full_thermal_bcf(0).real * model.bcf_scale for model in ω_models]) # plt.plot(ωs, [model.bcf_scale for model in ω_models]) @@ -23,14 +25,18 @@ for model in ω_models: plt.legend() for model in ω_models: - vals = abs(model.full_thermal_bcf(model.t)) - plt.plot(model.t, vals / vals.max(), linewidth=1) + vals = abs(model.bcf(model.t)) + plt.plot(model.t, model.bcf_scale * vals, linewidth=1) plt.xlim(-.1, 6) ωs = np.linspace(.01, 30, 1000) for model in ω_models: - plt.plot(ωs, model.full_thermal_spectral_density(ωs), linewidth=1, label=fr"$\omega_c={model.ω_c:.2f}$") + plt.plot(ωs, model.bcf_scale * model.spectral_density(ωs), linewidth=1, label=fr"$\omega_c={model.ω_c:.2f}$") plt.legend() +plt.xlabel(r"$\omega$") +plt.ylabel(r"$J(\omega)$") + +fs.export_fig("omega_sd", tikz=False, x_scaling=.49, y_scaling=.3) bath_scales = [1/(model.bcf_coefficients()[1][0].real.min()) for model in ω_models] therm_scales = [Δ/(model.bcf_scale) for model in ω_models] @@ -38,6 +44,7 @@ plt.plot(ωs, bath_scales) #plt.plot(ωs, therm_scales) fig, ax = plt.subplots() +ax.axhline(0, color="gray") for model, data in aux.model_data_iterator(ω_models): _, _, bar = pu.plot_with_σ( model.t, @@ -45,13 +52,20 @@ for model, data in aux.model_data_iterator(ω_models): ax=ax, label=fr"$\omega_c={model.ω_c:.2f}$", ) +ax.set_xlim(0, model.t[strobe_indices][3]) +ax.set_ylabel(r"$\langle H_\mathrm{I}\rangle$") +ax.set_xlabel(r"$\tau$") +ax.legend() +#fs.export_fig("omega_interactions", fig, tikz=False, x_scaling=.49, y_scaling=.3) fig, ax = plt.subplots() import scipy + for model, data in aux.model_data_iterator(ω_models): + energy = model.total_energy_from_power(data) * (1 / es.ergo(model.T)) _, _, bar = pu.plot_with_σ( model.t, - model.total_energy_from_power(data) * (1/es.ergo(model.T)), + energy, ax=ax, label=fr"$\omega_c={model.ω_c:.2f}$", strobe_frequency=Δ, @@ -60,7 +74,7 @@ for model, data in aux.model_data_iterator(ω_models): pu.plot_with_σ( model.t, - (model.total_energy_from_power(data) * (1/es.ergo(model.T))), + energy, ax=ax, # label=fr"$ω_c={model.ω_c:.2f}$", # strobe_frequency=Δ, @@ -69,27 +83,108 @@ for model, data in aux.model_data_iterator(ω_models): color=bar.lines[0].get_color(), ) - τ_off = 1 / (model.bcf_coefficients()[1][0].real.min()) - a0 = abs(model.bcf(0)) - τ_off = scipy.optimize.fsolve(lambda τ: abs(model.bcf(τ)) - (a0/300), 10) + # τ_off = 1 / (model.bcf_coefficients()[1][0].real.min()) + # # a0 = abs(model.bcf(0)) + # # τ_off = scipy.optimize.fsolve(lambda τ: abs(model.bcf(τ)) - (a0/300), 10) + # ax.axvline( + # τ_off, + # color=bar.lines[0].get_color(), + # linestyle="dotted", + # linewidth=1, + # zorder=-1, + # ) + + ts = model.t[strobe_indices] + e_s = energy.value[strobe_indices] + + + ind_min = e_s.argmin() + t_min = ts[ind_min] + e_min = e_s[ind_min] + ax.axvline( - τ_off, + t_min, color=bar.lines[0].get_color(), + ymin=0, + ymax=abs((e_min + 0.41) / (0.45 + 0.1)), linestyle="dotted", - linewidth=1, - zorder=-1, ) + ax.axhline( + e_min, + color=bar.lines[0].get_color(), + xmin=0, + xmax=t_min / model.t.max(), + linestyle="dotted", + ) + ax.set_xlabel(r"$\tau$") ax.set_ylabel(r"$(\langle H\rangle_\tau -\langle H\rangle_0)/\mathcal{W}_\mathrm{max}$") # ax.axhline(, color="grey", linewidth=1, linestyle="dashed") # ax.plot(model.t, model.L.operator_norm(model.t), linewidth=1) +ax.set_xlim(0, model.t.max()) +ax.set_ylim(-0.41, 0.1) ax.legend() -#fs.export_fig("omegas_total") +fs.export_fig("omegas_total", tikz=False, y_scaling=0.4) + +import scipy + +opt_powers = [] +opt_energies = [] +memories = [] +for model, data in aux.model_data_iterator(ω_models): + energy = model.total_energy_from_power(data) + ts = model.t[strobe_indices] + e_s = energy.value[strobe_indices] + + powers = np.divide(e_s, ts, where=ts > 0) + p_min = powers.min() + ind_min = powers.argmin() + # t_min = ts[ind_min] + # e_min = e_s[ind_min] + + opt_powers.append(abs(p_min)) + opt_energies.append(abs(e_s.min()) * (1 / es.ergo(model.T))) + + # bcf = model.full_thermal_bcf + # α0 = bcf(0).real + # obj = lambda τ: (abs(bcf(τ))/α0 - 1e-1) ** 2 + # τ_off = scipy.optimize.basinhopping( + # obj, 1 + # ).x + memories.append(3 / model.ω_c) + # memories.append(τ_off) + +fig, ax = plt.subplots() +lines = ax.plot(memories, opt_powers, marker="o", markersize=3) +ax.set_xlabel(r"$\tau_\mathrm{B}=3 / \omega_c$") +ax.set_ylabel(r"$P_\mathrm{max}$", color=lines[0].get_color()) + +ax2 = ax.twinx() + +ax2.plot([], []) +lines = ax2.plot(memories, opt_energies, marker="o", markersize=3) +ax2.set_ylabel( + r"$|\Delta \langle H\rangle|_\mathrm{max} / \mathcal{W}_\mathrm{max}$", + color=lines[0].get_color(), +) + +# def f(t, s, h): +# t = t-s +# out = np.zeros_like(t) +# np.log(t, out=out, where=t>1) + +# return out * h + +# opt,_ = scipy.optimize.curve_fit(f, memories, opt_powers) +# smooth_mem = np.linspace(memories[0], memories[-1], 1000) +# ax.plot(smooth_mem, f(smooth_mem, *opt)) # +plt.axvline(1 * np.pi * 2/Δ, color="gray", linestyle="dotted") +fs.export_fig("omega_energies_and_powers", fig, tikz=False, x_scaling=.5, y_scaling=.35) fig, ax = plt.subplots() for model, data in aux.model_data_iterator(ω_models): - fs.plot_with_σ( + pu.plot_with_σ( model.t, model.system_energy(data), ax=ax, @@ -105,16 +200,16 @@ ax.legend() fig, ax = plt.subplots() int_0 = ω_models[0].interaction_energy(aux.get_data(ω_models[0])).for_bath(0) -for model, data in aux.model_data_iterator(ω_models): - fs.plot_with_σ( +for model, data in aux.model_data_iterator(ω_models[1:]): + pu.plot_with_σ( model.t, - (model.interaction_energy(data) - int_0) * (1 / abs(int_0).max.value), + (abs(model.interaction_energy(data)) - abs(int_0)) * (1 / abs(int_0).max.value), ax=ax, bath=0, - label=fr"$ω_c={model.ω_c:.2f}$", + label=fr"$\omega_c={model.ω_c:.2f}$", linewidth=0.5, ) -#ax.legend() +ax.legend() fig, ax = plt.subplots() fig.set_size_inches(fs.get_figsize("poster", 0.49)) @@ -213,9 +308,11 @@ ax.legend() ax.set_xlim((0,10)) fs.export_fig("energy_shovel_preview") +from hops.util.utilities import relative_entropy + fig, ax = plt.subplots() for model,data in aux.model_data_iterator(ω_models): - fs.plot_with_σ(model.t, EnsembleValue(relative_entropy(data, strobe_indices[-1])), ax=ax, label=model.ω_c, strobe_frequency=Δ) + ax.plot(model.t, relative_entropy(data, strobe_indices[-1])[0], ax=ax, label=model.ω_c, strobe_frequency=Δ) ax.legend() fig, ax = plt.subplots()