Pybinding纳米带缺陷模型应用

最新推荐文章于 2024-08-28 07:35:17 发布

wwxy261

最新推荐文章于 2024-08-28 07:35:17 发布

阅读量471

点赞数

分类专栏：算法文章标签：数学

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

版权

本文介绍了如何利用Python库Pybinding构建二维材料的紧束缚模型，包括石墨烯的单层模型，并添加缺陷。通过设置不同参数，可以计算能带结构、绘制格点示意图以及计算局域密度态。该库允许在代码中直接生成模型，简化了复杂系统的建模过程。

摘要由CSDN通过智能技术生成

该软件包可以任意在实空间搭建模型，然后设置相应的hopping，同时可以增加缺陷，应力等，该软件包的模型是在Python代码中生成，不是输入文件类型。

个人感觉，目前该软件包的有以下几种应用

1. 计算能带，搭建TB模型后，就可以直接计算能带，避免重复造一些轮子工作。

2. 画出格点模型示意图，感觉对于复杂格点，比较方便。

3. 计算LDOS，或者DOS。

import pybinding as pb
import numpy as np
import matplotlib.pyplot as plt
from math import sqrt, pi
pb.pltutils.use_style()

def monolayer_graphene():
    t = 1      # [eV] nearest neighbour hopping
 
    lat = pb.Lattice(a1=[2*sqrt(3), 0],
                     a2=[sqrt(3)/2, -1.5])
            
                     
    lat.add_sublattices(('A', [0, 0]),
                        ('B', [sqrt(3)/2, -0.5]),
                        ('C', [sqrt(3), 0]),
                        ('D', [sqrt(3)*3/2, -0.5]))
    
    lat.add_hoppings(
        # inside the main cell
        ([0,  0], 'A', 'B', t),