版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

zhenzidan

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
帖子: 43
在線: 28.8小時(shí)
蟲(chóng)號(hào): 2560014
注冊(cè): 2013-07-23
專(zhuān)業(yè): 基礎(chǔ)物理學(xué)

[求助] pwscf計(jì)算超導(dǎo)轉(zhuǎn)變溫度已有2人參與

最近用pwscf計(jì)算超導(dǎo)轉(zhuǎn)變溫度，用的是最新版本(5.1)，按照師兄所言，取最后接近收斂的Tc=20K即可。但有個(gè)問(wèn)題：在最后輸出的lambda.out文件（見(jiàn)下文）中，默認(rèn)輸出的degauss變化很大(0.020到0.600)，從而導(dǎo)致N(Ef)的值和Tc的變化也很大，對(duì)比了5.0的版本，發(fā)現(xiàn)5.0版本的輸出的lambda.out中，degauss默認(rèn)輸出只有0.005到0.05，因此5.1算出來(lái)的結(jié)果就不能取最后接近收斂的Tc=20 K了，不知道我的理解是否正確，懇請(qǐng)各位指點(diǎn)，謝謝！
附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

回復(fù)此樓

1樓 2015-01-13 11:36:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Local_Spin

銅蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 201.5
散金: 15
帖子: 60
在線: 60.8小時(shí)
蟲(chóng)號(hào): 1121662
注冊(cè): 2010-10-13
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+2, 謝謝指導(dǎo) 2015-01-13 15:35:49
zhenzidan: 金幣+5, ★★★很有幫助 2015-01-13 18:33:50

ph的輸入文件中， el_ph_sigma控制degauss的步長(zhǎng)，el_ph_nsigma 控制degauss的步數(shù)。設(shè)置el_ph_sigma=0.005就是以前版本的默認(rèn)情況了。

贊一下

回復(fù)此樓

2樓2015-01-13 15:28:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhenzidan

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
帖子: 43
在線: 28.8小時(shí)
蟲(chóng)號(hào): 2560014
注冊(cè): 2013-07-23
專(zhuān)業(yè): 基礎(chǔ)物理學(xué)

引用回帖:

2樓: Originally posted by Local_Spin at 2015-01-13 15:28:11
ph的輸入文件中， el_ph_sigma控制degauss的步長(zhǎng)，el_ph_nsigma 控制degauss的步數(shù)。設(shè)置el_ph_sigma=0.005就是以前版本的默認(rèn)情況了。

謝謝前輩指點(diǎn)，那請(qǐng)問(wèn) el_ph_nsigma 設(shè)為多長(zhǎng)合適呢？像我的這個(gè)例子中，取Tc=20K有問(wèn)題嗎？再次感謝！

贊一下

回復(fù)此樓

3樓2015-01-13 18:38:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhenzidan

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
帖子: 43
在線: 28.8小時(shí)
蟲(chóng)號(hào): 2560014
注冊(cè): 2013-07-23
專(zhuān)業(yè): 基礎(chǔ)物理學(xué)

引用回帖:

前輩，還請(qǐng)賜教啊，實(shí)在不知道如何取舍，謝謝！

贊一下

回復(fù)此樓

4樓2015-01-15 21:07:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Local_Spin

銅蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 201.5
散金: 15
帖子: 60
在線: 60.8小時(shí)
蟲(chóng)號(hào): 1121662
注冊(cè): 2010-10-13
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

4樓: Originally posted by zhenzidan at 2015-01-15 21:07:56
前輩，還請(qǐng)賜教啊，實(shí)在不知道如何取舍，謝謝！...

你好，我一般取0.0025。你這個(gè)有沒(méi)有實(shí)驗(yàn)值做比較啊。感覺(jué)20 K對(duì)應(yīng)的smear過(guò)大了些，可能不對(duì)。你取0.0025看看，如果不收斂的話，你也許應(yīng)該加大計(jì)算的k點(diǎn)數(shù)。

贊一下

回復(fù)此樓

5樓2015-01-19 09:43:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhenzidan

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
帖子: 43
在線: 28.8小時(shí)
蟲(chóng)號(hào): 2560014
注冊(cè): 2013-07-23
專(zhuān)業(yè): 基礎(chǔ)物理學(xué)

引用回帖:

5樓: Originally posted by Local_Spin at 2015-01-19 09:43:13
你好，我一般取0.0025。你這個(gè)有沒(méi)有實(shí)驗(yàn)值做比較啊。感覺(jué)20 K對(duì)應(yīng)的smear過(guò)大了些，可能不對(duì)。你取0.0025看看，如果不收斂的話，你也許應(yīng)該加大計(jì)算的k點(diǎn)數(shù)。...

謝謝您的回復(fù)，這個(gè)沒(méi)有實(shí)驗(yàn)值！您的意思是讓我設(shè)置el_ph_sigma更小一些，然后看lambda.out中的degauss=0.0025時(shí)是否已經(jīng)收斂，如果收斂就取該degauss對(duì)應(yīng)的Tc，如果不收斂，則加大計(jì)算K點(diǎn)看看。再次感謝！

贊一下

回復(fù)此樓

6樓2015-01-19 16:12:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tsglss

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 41 (小學(xué)生)
金幣: 6014.2
紅花: 19
帖子: 1875
在線: 940小時(shí)
蟲(chóng)號(hào): 321641
注冊(cè): 2007-03-11
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

請(qǐng)問(wèn)樓主，你的那個(gè)超導(dǎo)轉(zhuǎn)變溫度的文獻(xiàn)是哪篇？謝謝

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

做好自己

7樓2015-01-19 20:57:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

du33lang

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10.6
散金: 5
紅花: 1
帖子: 22
在線: 9.6小時(shí)
蟲(chóng)號(hào): 2177287
注冊(cè): 2012-12-10
性別: MM
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★
liliangfang: 金幣+1, 謝謝交流 2015-03-24 07:46:05

樓主，你好。請(qǐng)問(wèn)進(jìn)行自恰運(yùn)算的時(shí)候是用的哪種方法呢？因?yàn)槲铱匆恍┵Y料，都說(shuō)LDA+U 還不能在計(jì)算聲子中實(shí)現(xiàn)。但是LDA+U 也是計(jì)算超導(dǎo)的最簡(jiǎn)單的方法了。

贊一下

回復(fù)此樓

8樓2015-03-07 11:21:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

pangrui1985

銅蟲(chóng) (小有名氣)

應(yīng)助: 25 (小學(xué)生)
金幣: 23.2
紅花: 6
帖子: 142
在線: 89.2小時(shí)
蟲(chóng)號(hào): 2808500
注冊(cè): 2013-11-18
專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★
liliangfang: 金幣+1, 謝謝交流 2015-03-24 07:46:01

引用回帖:

8樓: Originally posted by du33lang at 2015-03-07 11:21:46
樓主，你好。請(qǐng)問(wèn)進(jìn)行自恰運(yùn)算的時(shí)候是用的哪種方法呢？因?yàn)槲铱匆恍┵Y料，都說(shuō)LDA+U 還不能在計(jì)算聲子中實(shí)現(xiàn)。但是LDA+U 也是計(jì)算超導(dǎo)的最簡(jiǎn)單的方法了。

LDAU一般給的晶格常數(shù)都不對(duì)吧，都是用它算電子的性質(zhì)

贊一下

回復(fù)此樓

9樓2015-03-23 17:18:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

菲頂頂

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 32.6
散金: 45
紅花: 10
帖子: 631
在線: 250.2小時(shí)
蟲(chóng)號(hào): 1462090
注冊(cè): 2011-10-26
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

6樓: Originally posted by zhenzidan at 2015-01-19 16:12:59
謝謝您的回復(fù)，這個(gè)沒(méi)有實(shí)驗(yàn)值！您的意思是讓我設(shè)置el_ph_sigma更小一些，然后看lambda.out中的degauss=0.0025時(shí)是否已經(jīng)收斂，如果收斂就取該degauss對(duì)應(yīng)的Tc，如果不收斂，則加大計(jì)算K點(diǎn)看看。再次感謝！...

樓主弄清楚了嗎，是這個(gè)意思嗎？

贊一下

回復(fù)此樓

10樓2016-05-09 16:21:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 zhenzidan 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 335分 \| 材料與化工專(zhuān)碩 \| GPA 4.07 \| 有科研經(jīng)歷 +4	cccchenso 2026-03-23	4/200	2026-03-23 23:00 by 徐ckkk
[考研] 材料/農(nóng)業(yè)專(zhuān)業(yè)，07/08開(kāi)頭均可，過(guò)線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪！？！
[考研] 工科0856求調(diào)劑 +5	沐析汀汀 2026-03-21	5/250	2026-03-23 17:56 by 海瑟薇-
[考研] 306求調(diào)劑 +9	chuanzhu川燭 2026-03-18	9/450	2026-03-23 13:17 by luoyongfeng
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚(yú).. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 306求調(diào)劑 +5	來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 269專(zhuān)碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 354求調(diào)劑 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人來(lái)盛
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.

24小時(shí)熱門(mén)版塊排行榜

[求助] pwscf計(jì)算超導(dǎo)轉(zhuǎn)變溫度 已有2人參與

» 收錄本帖的淘帖專(zhuān)輯推薦

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] pwscf計(jì)算超導(dǎo)轉(zhuǎn)變溫度已有2人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: