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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf计算超导转变温度
131/2返回列表12下一页
查看: 6819  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 3/150 2026-08-08 03:27 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了，论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭