python kwant计算简单方格体系透射率

最新推荐文章于 2022-01-21 23:55:49 发布

wwxy261

最新推荐文章于 2022-01-21 23:55:49 发布

阅读量1.7k

点赞数

分类专栏：算法

本文链接：https://blog.csdn.net/wwxy1995/article/details/79811101

版权

import kwant
import matplotlib.pyplot as plt


def make_system(a=1, t=1.0, W=10, L=30):
    """搭建体系TB模型哈密顿量"""

    # 搭建device的哈密顿量
    lat = kwant.lattice.square(a)  # a为晶格常数
    syst = kwant.Builder()

    # 设置体系座位能on_site_energy
    syst[(lat(x,y) for x in range(L) for y in range(W))] = 4 * t
    # 设置体系的最近邻hooping
    syst[lat.neighbors()] = -t

    # 搭建左端的left_lead的哈密顿量
    lead = kwant.Builder(kwant.TranslationalSymmetry((-a, 0)))
    # 左端on_site_energy
    lead[(lat(0, j) for j in range(W))] = 4 * t
    # 左端最近邻hooping
    lead[lat.neighbors()] = -t

    # device和lead相连接
    # 与左端相连接
    syst.attach_lead(lead)
    syst.attach_lead(lead.reversed())
    # 返回体系的哈密顿矩阵
    return syst


def plot_conductance(syst, energies):
    """计算透射率"""
    data = []
    for energy in energies:
        # 生成体系的S矩阵
        smatrix = kwant.smatrix(syst, energy)