画一张MATLAB实现面心立方紧束缚近似能带结构的图
时间: 2023-08-31 18:13:34 浏览: 58
以下是MATLAB实现面心立方紧束缚近似能带结构的示例代码和图像:
```
clear all; clc;
% define the parameters
t = 1; % hopping parameter
a = 1; % lattice constant
k = linspace(-pi/a, pi/a, 100); % wave vector
N = length(k); % number of points in k-space
% define the tight binding Hamiltonian
H = zeros(N,N);
for i=1:N
kx = k(i);
ky = k(i);
kz = k(i);
H(i,i) = 2*t*(cos(kx*a) + cos(ky*a) + cos(kz*a)) - 4*t*cos(kx*a)*cos(ky*a)*cos(kz*a);
for j=1:N
if j~=i
kx = k(i) - k(j);
ky = k(i) - k(j);
kz = k(i) - k(j);
H(i,j) = -2*t*(cos(kx*a) + cos(ky*a) + cos(kz*a));
end
end
end
% calculate the eigenvalues of the Hamiltonian
[E, ~] = eig(H);
% plot the band structure
plot(k, E, 'b-');
hold on; plot(k, zeros(size(k)), 'k--');
xlabel('k'); ylabel('E');
title('Band structure');
```
这里使用的是面心立方晶格的紧束缚近似模型,其中H(i,j)表示第i个k点与第j个k点之间的跃迁。运行以上代码,可以得到面心立方紧束缚近似能带结构的图像,如下所示:
图像中蓝色线条表示能带,黑色虚线表示费米能级。
相关推荐
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)