24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2202  |  回复: 26

fo3mt2

金虫 (小有名气)


[交流] 电导率计算求助

大家好,

我想向大家请教一个用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点问题上,能否请这方面的专家给些建议和答案呢?如果可以,我可以把我的代码发给你们看下的。多谢啦。
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

danghaha

新虫 (文坛精英)



fo3mt2(金币+1): 谢谢参与
3楼2020-03-09 17:55:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

霸气侧漏的鱼

新虫 (小有名气)



fo3mt2(金币+1): 谢谢参与
电阻率不是直接测出的吗
13楼2020-03-18 10:49:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)



fo3mt2(金币+1): 谢谢参与
电导率是手动算吗?

发自小木虫Android客户端
14楼2020-04-11 17:25:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)



fo3mt2(金币+1): 谢谢参与
能自己写代码非常好!基于tb模型和Kubo公式计算电导,有地方可参考的。至于第4点,你理解的对,从1到499就可以了,最后一点,周期性边界条件决定了它与第1个点的等价,就不需要重复计算了。

发自小木虫Android客户端
15楼2020-11-15 02:41:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


你的哈密顿量如何体现出对k点的相关性?没有说清楚和想清楚对v速度的计算。

发自小木虫Android客户端
16楼2020-11-15 02:54:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


26楼2020-11-15 20:34:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)


引用回帖:
26楼: Originally posted by valenhou001 at 2020-11-15 20:34:55
参考这篇文献:
https://sci-hub.se/https://doi.org/10.1016/S0010-4655(99)00206-4
及相关的程序:
http://www.cpc.cs.qub.ac.uk/summaries/ADKI_v1_0.html

非常感谢你的建议和相关链接。我已经解决了电导率计算的问题了。再次非常感谢你的帮助。
27楼2020-11-17 03:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
syhorchid2楼
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): 谢谢参与
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   回复  
2020-11-15 14:21   回复  
相关版块跳转 我要订阅楼主 fo3mt2 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 329求调剂,一志愿西北工业大学,材料工程(085601) +7 小小机灵虫 2026-03-29 13/650 2026-04-03 18:38 by ls刘帅
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +5 redcom 2026-04-01 6/300 2026-04-03 15:15 by Howard28
[考研] 求调剂机会 +5 意染ivy 2026-04-03 5/250 2026-04-03 15:13 by qoooooo614
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 280求调剂 +17 咕噜晓晓 2026-04-02 18/900 2026-04-03 15:03 by hypershenger
[考研] 复试调剂 +5 春日来信- 2026-04-03 5/250 2026-04-03 15:01 by buqi613
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +4 哇呼哼呼哼 2026-04-02 4/200 2026-04-02 13:10 by yulian1987
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 化学0703 调剂 306分 一志愿211 +12 26要上岸 2026-03-28 12/600 2026-04-01 11:10 by chemdavid
[考研] 本2一志愿C9-333分,材料科学与工程,求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 一志愿华中师范化学332分求调剂 +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
[考研] 105500药学求调剂,一志愿山东大学药学,348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 291求调剂 +5 Y-cap 2026-03-29 6/300 2026-03-29 13:18 by mumin1990
信息提示
请填处理意见