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

24小时热门版块排行榜    

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

tsglss

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

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
做好自己
7楼2015-01-19 20:57:33
已阅   回复此楼   关注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的回帖

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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +73 AI8RaGaPaSCR 2026-08-07 76/3800 2026-08-14 08:45 by IUa8yhHGR7hP
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
关闭
请填处理意见
关闭