admin管理员组

文章数量:1332359

I want to get the eigenvectors of the Floquet operator and then derive them with respect to 'k' to be able to do the inner product that is in the Zak phase formula and then integrate over 'k' from -pi to pi. As I will get three eigenvectors, I expect to get three Zak phases.

I am writing code, and I expect help.

import numpy as np
from scipy.integrate import quad
from scipy.linalg import expm

# Define parameters
J1 = 1
J2 = 2
a = 0.25
T = 2

def floquet_operator(k):
    A1 = np.array([[0, J1 + a, 0],
                          [J1 + a, 0, J1],
                          [0, J1, 0]])
    A2 = np.array([
        [0, J2 * np.exp(-1j * k), 0],
        [J2 * np.exp(1j * k), 0, J2 * np.exp(-1j * k)],
        [0, J2 * np.exp(1j * k), 0]
    ], dtype=complex)
    
    exp_A1 = expm(-1j * A1 * (T / 2))
    exp_A2 = expm(-1j * A2 * (T / 2))
    return exp_A1 @ exp_A2

本文标签：

版权声明：本文标题：python - If I have Floquet operator, how to derive its eigenvectors w.r.t k as I want to calculate the Zak phase? - Stack Overfl 内容由网友自发贡献，该文观点仅代表作者本人， 转载请联系作者并注明出处：http://www.betaflare.com/web/1742316744a2451966.html， 本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容，一经查实，本站将立刻删除。