应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf计算超导转变温度
131/2返回列表12下一页
查看: 6197  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhenzidan

新虫 (初入文坛)

[求助] pwscf计算超导转变温度 已有2人参与

最近用pwscf计算超导转变温度，用的是最新版本(5.1)，按照师兄所言，取最后接近收敛的Tc=20K即可。但有个问题：在最后输出的lambda.out文件（见下文）中，默认输出的degauss变化很大(0.020到0.600)，从而导致N(Ef)的值和Tc的变化也很大，对比了5.0的版本，发现5.0版本的输出的lambda.out中，degauss默认输出只有0.005到0.05，因此5.1算出来的结果就不能取最后接近收敛的Tc=20 K了，不知道我的理解是否正确，恳请各位指点，谢谢！
附5.1版本的lambda.out文件：
     lambda = 0.189642 (  0.189293)  <log w>=  868.577K  N(Ef)= 2.282784 at degauss=0.020
     lambda = 0.225349 (  0.225151)  <log w>=  824.337K  N(Ef)= 2.501259 at degauss=0.040
     lambda = 0.268664 (  0.268490)  <log w>=  826.683K  N(Ef)= 2.807269 at degauss=0.060
     lambda = 0.305953 (  0.305711)  <log w>=  834.678K  N(Ef)= 3.116082 at degauss=0.080
     lambda = 0.337887 (  0.337570)  <log w>=  838.878K  N(Ef)= 3.394924 at degauss=0.100
     lambda = 0.368098 (  0.367726)  <log w>=  836.917K  N(Ef)= 3.646690 at degauss=0.120
     lambda = 0.399433 (  0.399031)  <log w>=  828.882K  N(Ef)= 3.878951 at degauss=0.140
     lambda = 0.432727 (  0.432336)  <log w>=  815.718K  N(Ef)= 4.096595 at degauss=0.160
     lambda = 0.467553 (  0.467219)  <log w>=  799.440K  N(Ef)= 4.301970 at degauss=0.180
     lambda = 0.503118 (  0.502885)  <log w>=  782.463K  N(Ef)= 4.495636 at degauss=0.200
     lambda = 0.538240 (  0.538141)  <log w>=  766.386K  N(Ef)= 4.676941 at degauss=0.220
     lambda = 0.571782 (  0.571836)  <log w>=  752.378K  N(Ef)= 4.844604 at degauss=0.240
     lambda = 0.603029 (  0.603241)  <log w>=  740.900K  N(Ef)= 4.997265 at degauss=0.260
     lambda = 0.631527 (  0.631891)  <log w>=  732.094K  N(Ef)= 5.133888 at degauss=0.280
     lambda = 0.656678 (  0.657184)  <log w>=  725.565K  N(Ef)= 5.253952 at degauss=0.300
     lambda = 0.678616 (  0.679248)  <log w>=  721.060K  N(Ef)= 5.357470 at degauss=0.320
     lambda = 0.697467 (  0.698207)  <log w>=  718.221K  N(Ef)= 5.444908 at degauss=0.340
     lambda = 0.713189 (  0.714022)  <log w>=  716.654K  N(Ef)= 5.517070 at degauss=0.360
     lambda = 0.726153 (  0.727065)  <log w>=  716.244K  N(Ef)= 5.574981 at degauss=0.380
     lambda = 0.736509 (  0.737483)  <log w>=  716.652K  N(Ef)= 5.619785 at degauss=0.400
     lambda = 0.744609 (  0.745634)  <log w>=  717.738K  N(Ef)= 5.652671 at degauss=0.420
     lambda = 0.750644 (  0.751709)  <log w>=  719.388K  N(Ef)= 5.674820 at degauss=0.440
     lambda = 0.754762 (  0.755859)  <log w>=  721.377K  N(Ef)= 5.687369 at degauss=0.460
     lambda = 0.757349 (  0.758468)  <log w>=  723.722K  N(Ef)= 5.691388 at degauss=0.480
     lambda = 0.758516 (  0.759650)  <log w>=  726.340K  N(Ef)= 5.687870 at degauss=0.500
     lambda = 0.758469 (  0.759614)  <log w>=  729.163K  N(Ef)= 5.677722 at degauss=0.520
     lambda = 0.757196 (  0.758346)  <log w>=  732.148K  N(Ef)= 5.661772 at degauss=0.540
     lambda = 0.755147 (  0.756297)  <log w>=  735.273K  N(Ef)= 5.640764 at degauss=0.560
     lambda = 0.752240 (  0.753387)  <log w>=  738.456K  N(Ef)= 5.615367 at degauss=0.580
     lambda = 0.748602 (  0.749743)  <log w>=  741.707K  N(Ef)= 5.586178 at degauss=0.600
lambda        omega_log          T_c
   0.18964       868.577              0.000
   0.22535       824.337              0.000
   0.26866       826.683              0.001
   0.30595       834.678              0.016
   0.33789       838.878              0.096
   0.36810       836.917              0.304
   0.39943       828.882              0.728
   0.43273       815.718              1.458
   0.46755       799.440              2.542
   0.50312       782.463              3.963
   0.53824       766.386              5.634
   0.57178       752.378              7.435
   0.60303       740.900              9.259
   0.63153       732.094             11.027
   0.65668       725.565             12.656
   0.67862       721.060             14.128
   0.69747       718.221             15.432
   0.71319       716.654             16.549
   0.72615       716.244             17.499
   0.73651       716.652             18.282
   0.74461       717.738             18.919
   0.75064       719.388             19.419
   0.75476       721.377             19.786
   0.75735       723.722             20.048
   0.75852       726.340             20.210
   0.75847       729.163             20.285
   0.75720       732.148             20.269
   0.75515       735.273             20.197
   0.75224       738.456             20.058
   0.74860       741.707             19.862
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-01-13 11:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Local_Spin

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指导 2015-01-13 15:35:49
zhenzidan: 金币+5, ★★★很有帮助 2015-01-13 18:33:50
ph的输入文件中， el_ph_sigma控制degauss的步长，el_ph_nsigma 控制degauss的步数。设置el_ph_sigma=0.005就是以前版本的默认情况了。
赞一下 回复此楼
2楼2015-01-13 15:28:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenzidan

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by Local_Spin at 2015-01-13 15:28:11
ph的输入文件中， el_ph_sigma控制degauss的步长，el_ph_nsigma 控制degauss的步数。设置el_ph_sigma=0.005就是以前版本的默认情况了。

谢谢前辈指点，那请问 el_ph_nsigma 设为多长合适呢？像我的这个例子中，取Tc=20K有问题吗？再次感谢！
赞一下 回复此楼
3楼2015-01-13 18:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenzidan

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by Local_Spin at 2015-01-13 15:28:11
ph的输入文件中， el_ph_sigma控制degauss的步长，el_ph_nsigma 控制degauss的步数。设置el_ph_sigma=0.005就是以前版本的默认情况了。

前辈，还请赐教啊，实在不知道如何取舍，谢谢！
赞一下 回复此楼
4楼2015-01-15 21:07:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Local_Spin

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by zhenzidan at 2015-01-15 21:07:56
前辈，还请赐教啊，实在不知道如何取舍，谢谢！...

你好，我一般取0.0025。你这个有没有实验值做比较啊。感觉20 K对应的smear过大了些，可能不对。你取0.0025看看，如果不收敛的话，你也许应该加大计算的k点数。
赞一下 回复此楼
5楼2015-01-19 09:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenzidan

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by Local_Spin at 2015-01-19 09:43:13
你好，我一般取0.0025。你这个有没有实验值做比较啊。感觉20 K对应的smear过大了些，可能不对。你取0.0025看看，如果不收敛的话，你也许应该加大计算的k点数。...

谢谢您的回复，这个没有实验值！您的意思是让我设置el_ph_sigma更小一些，然后看lambda.out中的degauss=0.0025时是否已经收敛，如果收敛就取该degauss对应的Tc，如果不收敛，则加大计算K点看看。再次感谢！
赞一下 回复此楼
6楼2015-01-19 16:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)
请问楼主，你的那个超导转变温度的文献是哪篇？谢谢

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
做好自己
7楼2015-01-19 20:57:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

du33lang

铁虫 (初入文坛)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-03-24 07:46:05
楼主，你好。请问进行自恰运算的时候是用的哪种方法呢？因为我看一些资料，都说LDA+U 还不能在计算声子中实现。但是LDA+U 也是计算超导的最简单的方法了。
赞一下 回复此楼
8楼2015-03-07 11:21:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pangrui1985

铜虫 (小有名气)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-03-24 07:46:01
引用回帖:
8楼: Originally posted by du33lang at 2015-03-07 11:21:46
楼主，你好。请问进行自恰运算的时候是用的哪种方法呢？因为我看一些资料，都说LDA+U 还不能在计算声子中实现。但是LDA+U 也是计算超导的最简单的方法了。

LDAU一般给的晶格常数都不对吧，都是用它算电子的性质
赞一下 回复此楼
9楼2015-03-23 17:18:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
6楼: Originally posted by zhenzidan at 2015-01-19 16:12:59
谢谢您的回复，这个没有实验值！您的意思是让我设置el_ph_sigma更小一些，然后看lambda.out中的degauss=0.0025时是否已经收敛，如果收敛就取该degauss对应的Tc，如果不收敛，则加大计算K点看看。再次感谢！...

楼主弄清楚了吗，是这个意思吗？
赞一下 回复此楼
10楼2016-05-09 16:21:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhenzidan 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭