This commit is contained in:
guanjihuan 2023-12-27 10:42:30 +08:00
parent d599ae945b
commit 118e32aed8
3 changed files with 21 additions and 2 deletions

View File

@ -1,7 +1,7 @@
[metadata] [metadata]
# replace with your username: # replace with your username:
name = guan name = guan
version = 0.1.65 version = 0.1.66
author = guanjihuan author = guanjihuan
author_email = guanjihuan@163.com author_email = guanjihuan@163.com
description = An open source python package description = An open source python package

View File

@ -1,6 +1,6 @@
Metadata-Version: 2.1 Metadata-Version: 2.1
Name: guan Name: guan
Version: 0.1.65 Version: 0.1.66
Summary: An open source python package Summary: An open source python package
Home-page: https://py.guanjihuan.com Home-page: https://py.guanjihuan.com
Author: guanjihuan Author: guanjihuan

View File

@ -551,4 +551,23 @@ def hamiltonian_of_kagome_lattice(kx, ky, t=1):
hamiltonian[1, 2] = 2*math.cos(k3_dot_a3) hamiltonian[1, 2] = 2*math.cos(k3_dot_a3)
hamiltonian = hamiltonian + hamiltonian.transpose().conj() hamiltonian = hamiltonian + hamiltonian.transpose().conj()
hamiltonian = -t*hamiltonian hamiltonian = -t*hamiltonian
return hamiltonian
# 超蜂窝晶格的哈密顿量(倒空间)
@guan.statistics_decorator
def hamiltonian_of_hyperhoneycomb_lattice(kx, ky, kz, t=1, a=1):
import cmath
import numpy as np
hamiltonian = np.zeros((4, 4), dtype=complex)
fx = (t*cmath.exp(1j*np.sqrt(3)/2*kx*a)+t*cmath.exp(-1j*np.sqrt(3)/2*kx*a))*cmath.exp(1j*kz/2*a)
fy = (t*cmath.exp(1j*np.sqrt(3)/2*ky*a)+t*cmath.exp(-1j*np.sqrt(3)/2*ky*a))*cmath.exp(1j*kz/2*a)
fz = cmath.exp(1j*kz*a)
hamiltonian[0, 1] = fx
hamiltonian[1, 0] = np.conj(fx)
hamiltonian[2, 3] = fy
hamiltonian[3, 2] = np.conj(fy)
hamiltonian[1, 2] = fz
hamiltonian[2, 1] = np.conj(fz)
hamiltonian[3, 0] = fz
hamiltonian[0, 3] = np.conj(fz)
return hamiltonian return hamiltonian