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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf计算超导转变温度
131/2返回列表12下一页
查看: 6508  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 297，工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
信息提示
关闭
请填处理意见
关闭