基于紧束缚模型,用MATLAB计算二维SSH模型,得到投影能带和原胞能带的重复结果

本文详细描述了如何使用MATLAB对二维SSH模型进行紧束缚计算,展示了计算结果中的投影能带和原胞能带,探讨了模型特性及其在材料科学和器件设计中的应用潜力。
摘要由CSDN通过智能技术生成

基于紧束缚模型,使用matlab计算二维SSH模型,结果是投影能带和原胞能带
注:这个是对文章的重复结果

ID:24200731781977720

胜的科研日记


基于紧束缚模型的计算在固体物理领域中具有广泛的应用。在本文中,我们将使用MATLAB软件对二维SSH模型进行计算,并展示其投影能带和原胞能带的结果。

紧束缚模型是一种常用的量子力学模型,用于描述电子的运动行为和能带结构。它通过将晶格的哈密顿量近似为相邻原子间的耦合来研究电子在晶体中的行为。二维SSH模型是紧束缚模型的一种特殊情况,它描述了一维电子链在二维平面上的运动行为。

在本文中,我们首先介绍了紧束缚模型的基本原理和二维SSH模型的数学表达式。然后,我们使用MATLAB软件编写了计算二维SSH模型的程序。程序中,我们采用了适当的边界条件和参数设置,确保了计算的准确性和可靠性。

接下来,我们展示了计算得到的投影能带和原胞能带的结果。投影能带是描述电子在能级上的分布情况的一种方法,它可以帮助我们理解材料的导电性和电子结构。原胞能带则描述了电子在晶体中的动力学行为,反映了晶格结构对电子行为的影响。

通过对投影能带和原胞能带的分析,我们发现了二维SSH模型的一些有趣特性。例如,我们观察到能带结构中的能隙和拓扑性质,这对于设计新的材料和开发新的器件具有重要的意义。此外,我们还研究了参数对能带结构的影响,并讨论了一些与能带结构相关的物理量。

总之,本文使用MATLAB软件对二维SSH模型进行了计算,并展示了其投影能带和原胞能带的结果。通过对计算结果的分析和讨论,我们揭示了二维SSH模型的一些有趣特性,并对其在材料科学和器件设计中的应用进行了展望。希望本文对读者对紧束缚模型和二维SSH模型有着更深入的理解,并为相关领域的研究提供一定的参考价值。

以上相关代码,程序地址:http://wekup.cn/731781977720.html

  • 6
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
束缚模型是一种离散化的模型,用于描述电子在晶格中的行为。一维SSH模型是指具有两种不同类型的原子相间排列的一维晶格模型。下面给出用MATLAB计算SSH模型本征值的代码: ```matlab % 定义模型参数 t1 = 1; t2 = 0.5; Delta = 0.1; % 定义晶格大小和边界条件 N = 50; L = N/2; bc = 'periodic'; % 构造哈密顿矩阵 H = zeros(N,N); for i = 1:N if mod(i,2) == 1 H(i,i) = Delta/2; H(i+1,i) = t1; H(i,i+1) = t1; else H(i,i) = -Delta/2; H(i+1,i) = t2; H(i,i+1) = t2; end end % 处理周期边界条件 if strcmp(bc,'periodic') H(1,N) = t1; H(N,1) = t1; H(N-1,N) = t1; H(N,N-1) = t1; end % 计算本征值和本征态 [V,D] = eig(H); % 绘制能图 k = linspace(-pi,pi,100); E = zeros(size(k)); for i = 1:length(k) Hk = zeros(N,N); for j = 1:N if mod(j,2) == 1 Hk(j,j) = Delta/2; Hk(j+1,j) = t1*exp(-1i*k(i)); Hk(j,j+1) = t1*exp(1i*k(i)); else Hk(j,j) = -Delta/2; Hk(j+1,j) = t2*exp(-1i*k(i)); Hk(j,j+1) = t2*exp(1i*k(i)); end end if strcmp(bc,'periodic') Hk(1,N) = t1*exp(-1i*k(i)); Hk(N,1) = t1*exp(1i*k(i)); Hk(N-1,N) = t1*exp(-1i*k(i)); Hk(N,N-1) = t1*exp(1i*k(i)); end E(i) = min(eig(Hk)); end plot(k,E,'k') xlabel('k') ylabel('E') ``` 在此代码中,我们首先定义了模型参数,包括两种不同的跃迁能级和离子之间的交互作用强度。然后,我们定义了晶格大小和边界条件,并使用这些参数构造了哈密顿矩阵。接下来,我们使用MATLAB内置的函数“eig”计算哈密顿矩阵的本征值和本征态。最后,我们使用计算出的本征值绘制了能图。 运行此代码,即可得到SSH模型的能图,从而得到模型的本征值。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值