hiro/HOPSFlow-Paper
1
0
Fork 0
mirror of https://github.com/vale981/HOPSFlow-Paper synced 2025-03-05 09:41:40 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

increase weight of off-axis element

Browse source
This commit is contained in:
valentin.boettcher@mailbox.tu-dresden.de 2023-07-31 16:50:06 -04:00
parent 8e03fbd401
commit aa3ae06cd8
No known key found for this signature in database
GPG key ID: E034E12B7AF56ACE
3 changed files with 10 additions and 10 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
python/otto_motor/subprojects/cycle_shift/cycle_shift.org
View file

@ -988,7 +988,7 @@ hamiltonian.
#+begin_src jupyter-python :tangle tangle/rot.py
rot_models = []
weights = [.3]
weights = [.5]
for weight in weights:
off_ax = sc.make_model(0, 0)
off_ax.H_bias = ConstantMatrix(weight / 2 * qt.sigmax().full())
@ -1001,17 +1001,17 @@ hamiltonian.
#+begin_src jupyter-python
τs = rot_models[0].t
#plt.plot(τs, np.einsum('tij,ij->t', rot_models[0].H(τs), qt.sigmay().full()).real)
plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
# plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
# plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
# plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
H = rot_models[0].H
plt.plot(τs, list(map(lambda t: get_energy_gap(H(t)), τs)), color="black")
#+end_src
#+RESULTS:
:RESULTS:
| <matplotlib.lines.Line2D | at | 0x7f36244f3340> |
[[file:./.ob-jupyter/e2600fe61c160903c60b7159b388b82b5ce667d0.svg]]
| <matplotlib.lines.Line2D | at | 0x7f36470cdd00> |
[[file:./.ob-jupyter/d5d7c229833877657c2abc95a30a085919c05d40.svg]]
:END:
** Integration

6
python/otto_motor/subprojects/cycle_shift/tangle/cycle_shift.py
View file

@ -373,9 +373,9 @@ fs.export_fig(f"energy_change_off_axis", x_scaling=2, y_scaling=0.7)
τs = rot_models[0].t
#plt.plot(τs, np.einsum('tij,ij->t', rot_models[0].H(τs), qt.sigmay().full()).real)
plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
# plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
# plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
# plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
H = rot_models[0].H
plt.plot(τs, list(map(lambda t: get_energy_gap(H(t)), τs)), color="black")

2
python/otto_motor/subprojects/cycle_shift/tangle/rot.py
View file

@ -22,7 +22,7 @@ import itertools
from hops.util.dynamic_matrix import ConstantMatrix
rot_models = []
weights = [.3]
weights = [.5]
for weight in weights:
off_ax = sc.make_model(0, 0)
off_ax.H_bias = ConstantMatrix(weight / 2 * qt.sigmax().full())