二维六角晶格体系极化率RPA的计算结果总结

六角形布里渊区不方便进行直接计算，所以需要进行平移操作，将布里渊区平移为如下平行四边形，平移后的均匀取样方法见以前的文章。

能量可以由紧束缚模型计算，也可以从第一性原理计算中读取，然后就可以在平移之后的布里渊区，均匀的画出能带和费米面，Matlib代码如下：

Nx = 200;
Ny = 200;
a = 4.18879028;
b1 = a*[0  ,-1.7321];
b2 = a*[1.5,-0.8660];
e1 = b1/Nx;
e2 = b2/Ny;
KPOS = zeros(Nx,Ny,2);
for i=1:Nxfor j=1:NyKPOS(i,j,:) = i*e1 + j*e2;end
end
Eks = zeros(Nx,Ny,1);
for i = 1:Nxfor j = 1:NyHk = Hamiltonian_TaS2_k(KPOS(i,j,1),KPOS(i,j,2));[~,Ek]=eig(Hk); Eks(i,j) =  Ek(1,1);end
end
Ef=0;
mesh(KPOS(:,:,1),KPOS(:,:,2),Eks);
[C,h] = contour(KPOS(:,:,1),KPOS(:,:,2),Eks,[Ef,0,0]);

然后是极化率的计算，由于已经均匀取样，就可以按照如下方法计算

mu = 0;
T = 0.001;
delta = 0.001;
chi = zeros(Nx,Ny);
for i = 1:Nxfor j = 1:Nyfor k=1:Nxfor l=1:Nyindex_kq_x = mod(i+k,Nx)+1;index_kq_y = mod(j+l,Ny)+1;Ek = Eks(k,l);Ekq = Eks(index_kq_x,index_kq_y);chi(i,j) = chi(i,j) + (Fermi_funtion(Ek,mu,T)-Fermi_funtion(Ekq,mu,T))/(Ek-Ekq+1i*delta);endendend
end
mesh(KPOS(:,:,1),KPOS(:,:,2),-real(chi))

在以上计算中，我们发现，动量k已经被我们用均匀坐标整数化，只在画图时体现出来。到这里，我们基本解决了六角晶格极化率计算的问题。

沿着高对称线计算（Gamma-M-Gamma）

Gamma-K-Gamma

经过RPA计算以后，原来的峰会更尖。

总结，计算的结果，高度依赖于能带在费米面附近的形状。

二维六角晶格体系极化率RPA的计算结果总结相关推荐

用二维数组进行学生的成绩排序和计算
//在控制台中输入学生的学号,语文,数学,英语成绩:然后输出各科成绩,平均成绩//和总成绩Console.WriteLine("总学生数:");int num = Convert. ...
二维静电场有限元matlab,实验2 MATLAB有限元计算
MATLAB电磁场有限元计算实验目的: 1.了解有限元算法的原理,熟练运用MATLAB环境的PDE工具. 2.熟练运用PDE工具分析简单的电磁场边值问题. 实验内容: 一. 有限元简介在电磁场的计 ...
kwant搭建六角晶格矩形区域散射计算LDOS和透射率
主要是熟悉API的调用 import numpy as np import kwant import matplotlib.pyplot as pltrt3 = np.sqrt(3)a1 = np.a ...
qr二维码遇到的一些问题
编码模式,大体定义如下: typedef enum { QR_MODE_NUL = -1, ///< Terminator (NUL character). Internal ...
关于QR二维码的编码模式
编码模式,大体定义如下: typedef enum { QR_MODE_NUL = -1, ///< Terminator (NUL character). Internal ...
单层MoS2和WS2晶体/Zn(Ⅱ)-PDA诱导合成二维连续稳定ZIF-8膜/新型大孔的二维结晶聚酰亚胺COFs(PI-COFs)
硼元素是在元素周期表中的III 族元素,也是碳元素的近邻元素,促使人们寻找硼的二维晶体--硼烯(borophene).三价硼原子的外层电子结构为2s22p1,相比碳原子中典型的sp2或sp3杂化,由于 ...
java后台识别二维码的方法
1.google zxing 缺点:识别准确度一般,简单标准的二维码可以使用这种方式. 优点:使用起来很简单,易上手. package vip.xiaonuo.common.util;import c ...
二维码原理--开胃小知识
二维码是什么原理?不知道从什么时候开始,我们的生活突然之间就充满了二维码,看网页要扫二维码,加好友要扫二维码,现在连楼下卖草莓的大爷都支持扫码支付,那么,你有没有想过,这个长得很奇怪的二维码,到底是怎 ...
Opencv+ZBar识别条形码、二维码
文章目录 Opencv识别条形码.二维码 1.ZBar环境配置 2.一维码(条形码)识别 3. 二维码的识别 4. Opencv识别二维码结束 Opencv识别条形码.二维码最近的一次作业,恰好之 ...

二维六角晶格体系极化率RPA的计算结果总结

二维六角晶格体系极化率RPA的计算结果总结相关推荐

最新文章

热门文章