在SCTK的后處理中，calculation = “lambda_mu_k”可以計算出lambda和mu的值

PATH/sctk.x < sctk.in > lambda.out

lambda.out文件

##### Integrate gapeq #####mu : 0.36011E+00lambda : 0.18219E+01##### Write mu.frmsf and lambda_frmsf #####Max : 0.53735334 Min : 0.12076953Max : 2.86872787 Min : 0.23379861##### Clock summary #####

以及生成兩個文件：lambda.frmsf, mu.frmsf
其計算過程如下：
sctk_gapeq_rhs.f90/SUBROUTINE make_lambda_mu_f

lambda = SUM(ZF(1,b_low:b_high,1:nks)*dosd(b_low:b_high,1:nks))mu = SUM( dltf(1,b_low:b_high,1:nks)*dosd(b_low:b_high,1:nks))weight = SUM( dosd(b_low:b_high,1:nks))lambda = lambda / weightmu = mu / weight!WRITE(stdout,'(7x," weight : ",e15.5)') weightWRITE(stdout,'(7x," mu : ",e15.5)') muWRITE(stdout,'(7x,"lambda : ",e15.5)') lambda

其中：lambda為輸出的，我們計算的lambda
mu為輸出的，我們計算的mu
weight為權重，正常沒有輸出，我加了一行進去

ZF為重正化因子，其維度為（2，b_low:b_high,nks）2是自旋，b_low是計算2.2步我們在輸入文件中給出的elph_nbnd_min和elph_nbnd_max，是費米面附近的最小能帶和最大能帶。nks是總K點數，nk1nk2nk3。它是lambda.frmsf文件的后（b_low-b_high+1）*nks個數據。

dosd是對于每個K點的態密度，它沒有輸出文件。它的維度是(b_low:b_high,nks)

McMillan formula中的量

一共四個變量
tc：輸出文件tc.out中的tc就是公式中的tc，且單位一致
輸出文件lambda.out中的lambda就是公式中的lambda，且單位一致
wln：尚未找到 平均聲子頻率，在代碼中為omg
mustar:尚未找到

ambda.out中的mu更像是muC和muS的和。有待考證。

總結

以上是生活随笔為你收集整理的SCTK计算lambda的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。