| 查看: 2044 | 回复: 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点问题上，能否请这方面的专家给些建议和答案呢？如果可以，我可以把我的代码发给你们看下的。多谢啦。 |
回复此楼» 猜你喜欢
心脉受损
已经有5人回复
-
博士读完未来一定会好吗
已经有15人回复
-
Springer期刊投稿求助
已经有4人回复
-
读博
已经有3人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
» 抢金币啦！回帖就可以得到:
中国地质大学（武汉）杨明教授组招收2026年博士
+2/414
-
寻找承接企业/科研机构/高校：启明人才计划（上海及长三角方向）
+1/88
-
加拿大/英属哥伦比亚大学曹彦凯课题组招收全奖博士/博后 [机器学习/优化/控制方向]
+1/81
-
加拿大/英属哥伦比亚大学曹彦凯课题组招收全奖博士/博后 [机器学习/优化/控制方向]
+1/80
-
天津大学精密测试技术及仪器全国重点实验室柔性电子技术实验室专任副研究员招聘启事
+1/75
-
QS TOP100英国南安普顿大学数字健康与生医工招博后，博士，Fellowship，访问学者
+1/30
-
澳门科技大学2026年数学博士招生—计算物理与数学课题组:计算流体与相场方法
+1/22
-
浙江工业大学国家优青朱艺涵团队在固态电池解构与设计方向招收2026年博士生2名
+1/22
-
天津医科大学基础医学院张恒课题组博士后招聘
+1/21
-
招聘2026年入学博士生
+1/7
-
哈尔滨工业大学（深圳）赵怡潞课题组诚招博士后
+1/6
-
招若干有分子生物，细胞培养，动物实验背景的人员(中山大学)
+1/5
-
东北师范大学光电转换材料与器件团队招收2026级博士研究生
+1/5
-
北理工柔性电子国家杰青团队招【博士后】【博士】【科研助理】
+1/4
-
招收26年秋季入学博士生（北科大高精尖学院 力学超材料/机器学习/增材制造相关方向）
+1/4
-
澳大利亚南昆士兰大学（UniSQ）量子点课题组 招收CSC全奖博士生
+1/3
-
深圳先进院三院院士成会明团队诚聘液流电池/高分子合成等方向博后、科研助理、工程师
+1/3
-
上海大学长江学者钟云波教授团队招收外场冶金或材料加工方向2026年博士研究生
+1/2
-
浙江师范大学夏永姚/黄健航教授团队招收2026级博士
+1/1
-
东北师范大学国家杰青汤庆鑫教授团队招收博士研究生
+1/1
3楼2020-03-09 17:55:19
电阻率不是直接测出的吗
13楼2020-03-18 10:49:25
14楼2020-04-11 17:25:26
★
fo3mt2(金币+1): 谢谢参与
fo3mt2(金币+1): 谢谢参与
|
能自己写代码非常好！基于tb模型和Kubo公式计算电导，有地方可参考的。至于第4点，你理解的对，从1到499就可以了，最后一点，周期性边界条件决定了它与第1个点的等价，就不需要重复计算了。 发自小木虫Android客户端 |
15楼2020-11-15 02:41:56
16楼2020-11-15 02:54:42
26楼2020-11-15 20:34:55
27楼2020-11-17 03:37:13
简单回复
2020-03-09 17:20
回复
fo3mt2(金币+1): 谢谢参与
tzynew4楼
2020-03-09 17:56
回复
fo3mt2(金币+1): 谢谢参与
8 发自小木虫Android客户端
2020-03-09 20:12
回复
fo3mt2(金币+1): 谢谢参与
2020-03-10 10:20
回复
fo3mt2(金币+1): 谢谢参与
2020-03-10 20:41
回复
fo3mt2(金币+1): 谢谢参与
2020-03-11 22:36
回复
fo3mt2(金币+1): 谢谢参与
暧暧远人村9楼
2020-03-12 18:32
回复
fo3mt2(金币+1): 谢谢参与
1
danghaha10楼
2020-03-12 19:54
回复
psylhh11楼
2020-03-13 20:26
回复
fo3mt2(金币+1): 谢谢参与
nono200912楼
2020-03-15 11:04
回复
fo3mt2(金币+1): 谢谢参与
·
icanspll17楼
2020-11-15 11:34
回复
shehuifei18楼
2020-11-15 11:35
回复
zzhello19楼
2020-11-15 11:35
回复
gofast11220楼
2020-11-15 11:35
回复
likecitywo21楼
2020-11-15 11:37
回复
Haibo_22楼
2020-11-15 11:38
回复
howtoge23楼
2020-11-15 11:40
回复
yaqianye24楼
2020-11-15 11:40
回复
蓉蓉番薯25楼
2020-11-15 14:21
回复