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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf计算超导转变温度
51/1返回列表
查看: 6847  |  回复: 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的回帖

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的回帖
查看全部 13 个回答

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的回帖

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的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 12/600 2026-08-12 01:27 by 南京中西医
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +4 Tide man 2026-08-11 4/200 2026-08-12 00:39 by liuqidi
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 29/1450 2026-08-12 00:24 by zhanghaozhu
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭