| 查看: 2173 | 回复: 26 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
电导率计算求助
|
|||
|
大家好, 我想向大家请教一个用kubo方法计算电导率的问题。我是用下面的这个公式计算电导率的。 σ_xx=e^2 * ℏ /4π 1/A ∑_ktr[v_x * (g^r-g^a ) * v_x * (g^r-g^a )] 其中,v_x是速度算符(dh/dkx);g^r和g^a是延迟和超前格林函数:g^r=1/(e_f-h+i*gamma); g^a=1/(e_f-h+i*gamma)。 假设我的hamiltonian是一个2乘2的矩阵,其对应的本征矢也是两个分别为2乘1的向量(|ϕ_n>)。对于迹部分的计算我想用下面的两种方式计算。 第一种计算方式,设o= v_x * (g^r-g^a ) * v_x * (g^r-g^a ),tr[v_x * (g^r-g^a ) * v_x * (g^r-g^a )]=tr(o)=<ϕ_11|o_11|ϕ_11>+<ϕ_21|o_22|ϕ_21>,其中o_11和o_22分别是这个2乘2矩阵的第1行第1列元素和第2行第2列元素。|ϕ_11>和|ϕ_21>分别是本征矢的第1行第1列元素和第2行第1列元素。 第二种计算方式,是把格林函数中的hamiltonian矩阵转化为hamiltonian的本征值。 σ_xx=e^2 * ℏ /4π*1/a * ∑_ktr[v_x * (g^r-g^a ) * v_x * (g^r-g^a )]=e^2 * ℏ /4π*1/a * ∑_k∑_n<ϕ_n|[v_x * (g^r-g^a ) * v_x * (g^r-g^a )]|ϕ_n>=e^2 * ℏ /4π*1/a * ∑_k∑_n∑_m<ϕ_n|[v_x * (g^r-g^a ) *|ϕ_m><ϕ_m|* v_x * (g^r-g^a )]|ϕ_n>=e^2 * ℏ /4π*1/a * ∑_k∑_n,m<ϕ_n|v_x|ϕ_m> * (g^r_m-g^a_m) * <ϕ_m|v_x|ϕ_n> * (g^r_n-g^a_n)。 其中,g^r_n=1/(e_f-e_n+i*gamma); g^a_n=1/(e_f-e_n+i*gamma);e_n是hamiltonian的本征值。也就是说这里把格林函数由一个矩阵转化成一个数。而迹的部分也就变成列下面的计算组合。 ∑_n,m<ϕ_n|v_x|ϕ_m> * (g^r_m-g^a_m) * <ϕ_m|v_x|ϕ_n> * (g^r_n-g^a_n)=<ϕ_1|v_x|ϕ_1> * (g^r_1-g^a_1) * <ϕ_1|v_x|ϕ_1> * (g^r_1-g^a_1)+<ϕ_1|v_x|ϕ_2> * (g^r_2-g^a_2) * <ϕ_2|v_x|ϕ_1> * (g^r_1-g^a_1)+<ϕ_1|v_x|ϕ_1> * (g^r_1-g^a_1) * <ϕ_1|v_x|ϕ_2> * (g^r_2-g^a_2)+<ϕ_2|v_x|ϕ_2> * (g^r_2-g^a_2) * <ϕ_2|v_x|ϕ_2> * (g^r_2-g^a_2) 这里的|ϕ_1>和|ϕ_2>分别是hamiltonian的第1个和第2个本征矢,g^r_1=1/(e_f-e_1+i*gamma)和g^r_2=1/(e_f-e_2+i*gamma);而e_1和e_2分别是hamiltonian的第1个和第2个本征值。 1. 请问这两种方法是否都正确呢?是否彼此对等呢?我是指最终计算出的结果是否一样呢?如果不正确,要如何修改呢? 2. 这两种方法,我都需要在对k点求和的时候乘以dk_x和dk_y。假设我的第一布里渊区被划分成500乘500的k网格,dk_x=b1/499;dk_x=b2/499,其中b1和b2分别是倒格矢空间中元胞的基矢。因为每个基矢都被这500个k点分割成了499份。 我的这种理解是否正确呢? 3. 因为1/A ∫(dk_x * dk_y)/(4π^2 )=1/N;其中a和n分别是实空间中晶胞的面积和倒格矢空间中的k点总数(N=500*500)。这样一来,帖子一开始的电导率计算公式也可以写成下面的形式。 σ_xx=e^2 * ℏ /4π*1/N * ∑_ktr[v_x (g^r-g^a ) v_x (g^r-g^a )]或者σ_xx=e^2 * ℏ /4π*1/N * ∑_k∑_(n,m)⟨ϕ_n│v_x│ϕ_m ⟩ * ⟨ϕ_m│v_x│ϕ_n ⟩ * (g^r_n -g^a_n ) * (g^r_m - g^a_m) 我的这种理解是否正确呢? 4. 因为周期性的因素,倒格矢空间中每一行的第1个和第500个k格点上的hamiltonian应该是完全相同的。这样一来,我在计算电导率加和的时候,是不是只需要把从第1个k点到第499个k点上的电导率相加就行了,而不需要再加入第500个k点上的电导率。 我的这种理解是否正确呢? 我在用fortran写了两个电导率计算的程序。一个是基于第1种计算方式,另一个是基于第2种计算方式。可是都出了问题。我想可能错误出在上面的4点问题上,能否请这方面的专家给些建议和答案呢?如果可以,我可以把我的代码发给你们看下的。多谢啦。 |
» 猜你喜欢
一志愿西南交大,材料专硕317求调剂
已经有5人回复
304求调剂
已经有7人回复
0703化学一志愿211 总分320求调剂
已经有5人回复
四川大学085601材料工程专硕 初试294求调剂
已经有4人回复
工科,求调剂
已经有3人回复
333求调剂
已经有3人回复
(081700)化学工程与技术-298分求调剂
已经有35人回复
308 085701 四六级已过求调剂
已经有13人回复
315求调剂
已经有10人回复
学硕285求调剂
已经有45人回复
» 抢金币啦!回帖就可以得到:
鲍红丽课题组 研究生招生启事
+1/479
【代人发帖勿站内信】【成都征男友】89年未婚博士漂亮女孩
+1/100
华南师范大学(211)博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/96
海南大学徐月山老师招生第二批博士名额2~3个,2026年9月份入学(高端设备开发方向)
+1/86
郑州大学付永柱课题组诚招2026年入学博士生
+1/83
sciencedirect 网页打不开了,怎么解决啊?
+1/39
南开大学罗景山教授课题组招聘博士后(光/电催化、钙钛矿光伏方向)
+2/38
中科院化学所 宋延林 课题组招聘合成化学方向博士后(开展打印合成化学方向研究)
+1/33
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1名
+1/26
上海交通大学药学院高杨课题组诚聘博士后
+1/16
食品科学与工程083200-297分
+1/8
【第三轮招生,31号截止】澳科大诚招2026秋全奖博士研究生(纳米医学/生物材料方向)
+1/7
哈尔滨理工大学材料与化工学院刘刚老师团队招募“申请-考核”制博士
+2/6
上海第二工业大学-朱大海课题组招生(过线就能调剂)
+1/5
总分291英一、数二-中核论文(导师一作、本人二作)-一志愿211学硕-自我介绍详细
+1/3
欢迎报考浙江工业大学埃克塞特联合学院硕士研究生
+1/3
食品调剂-总分291-有中核论文导师一作、本人二作-一志愿南昌大学学硕-介绍详细
+1/3
333求调剂
+1/2
哈尔滨工业大学(深圳)绿色化学团队招收2026年秋季入学博士生
+1/2
315求调剂
+1/1
26楼2020-11-15 20:34:55
3楼2020-03-09 17:55:19
13楼2020-03-18 10:49:25
14楼2020-04-11 17:25:26
简单回复
tzynew4楼
2020-03-09 17:56
回复
fo3mt2(金币+1): 谢谢参与
8 发自小木虫Android客户端
danghaha10楼
2020-03-12 19:54
回复













回复此楼
电阻率不是直接测出的吗