A study of the Haken, Strobl, Reineker model.
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
Main function is hsrmodel()
as found in hsrmodel.py
.
First we create the coherent part based on .
Which results to:
for n in range(N):
for m in range(N):
idx1 = fidx(n,m,N)
for q in range(N):
idx2 = fidx(q,m,N)
L[idx1][idx2] = L[idx1][idx2] -1j*J[abs(n-q)]
idx2 = fidx(n,q,N)
L[idx1][idx2] = L[idx1][idx2] +1j*J[abs(q-m)]