 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 4899 | 回復(fù): 3 | |||
[求助]
MgB2超導(dǎo)轉(zhuǎn)變溫度
|
|
最近在學(xué)習(xí)用quantum espresso計算超導(dǎo)轉(zhuǎn)變溫度，所以就用比較成熟的MgB2體系作為例子進(jìn)行實際的計算。但是計算結(jié)果與實驗值差距很大，實驗值是39.3 K，而我最高只能算到4.1 K，想盡各種辦法也沒能有效提高，不知道是參數(shù)設(shè)置的不對還是贗勢選擇有問題�？吹接行┪墨I(xiàn)用QE計算的MgB2的轉(zhuǎn)變溫度與實驗值復(fù)合的很好。但是自己計算的就是不對，有這方面經(jīng)驗的蟲友請幫忙看一下輸入文件，有什么地方是可以改進(jìn)或是改正的。謝謝了！ cat > mgb2.scf.fit.in << EOF &control calculation='scf' restart_mode='from_scratch', prefix='mgb2', etot_conv_thr = 1.0d-8 forc_conv_thr = 1.0d-7 pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 4, celldm(1) =5.8124, celldm(3)=1.1483, nat= 3, ntyp= 2, ecutwfc =60.0, ecutrho =600, occupations='smearing', smearing='methfessel-paxton', degauss=0.04, la2F = .true., / &electrons conv_thr = 1.0d-8 mixing_beta = 0.7 / &IONS / &CELL press = 0.0 cell_factor = 1.5d0 / ATOMIC_SPECIES Mg 24.305 Mg.pw91-np-van.UPF B 10.811 B.pw91-n-van_ak.UPF ATOMIC_POSITIONS {crystal} Mg 0.000000000 0.000000000 0.000000000 B 0.333333333 0.666666667 0.500000000 B 0.666666667 0.333333333 0.500000000 K_POINTS {automatic} 32 32 32 0 0 0 EOF $ECHO " running the scf calculation with dense k-point grid...\c" $PW_COMMAND < mgb2.scf.fit.in > mgb2.scf.fit.out check_failure $? $ECHO " done" # # SCF at k-mesh good enough for phonons # cat > mgb2.scf.in << EOF &control calculation='scf' restart_mode='from_scratch', prefix='mgb2', etot_conv_thr = 1.0d-8 forc_conv_thr = 1.0d-7 pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 4, celldm(1) =5.8124, celldm(3)=1.1483, nat= 3, ntyp= 2, ecutwfc = 60.0,ecutrho = 600, occupations='smearing', smearing='methfessel-paxton', degauss=0.04 / &electrons conv_thr = 1.0d-8 mixing_beta = 0.7 / &IONS / &CELL press = 0.0 cell_factor = 1.5d0 / ATOMIC_SPECIES Mg 24.305 Mg.pw91-np-van.UPF B 10.811 B.pw91-n-van_ak.UPF ATOMIC_POSITIONS {crystal} Mg 0.000000000 0.000000000 0.000000000 B 0.333333333 0.666666667 0.500000000 B 0.666666667 0.333333333 0.500000000 K_POINTS {automatic} 16 16 16 0 0 0 EOF $ECHO " running the scf calculation...\c" $PW_COMMAND < mgb2.scf.in > mgb2.scf.out check_failure $? $ECHO " done" # cat > mgb2.elph.in << EOF Electron-phonon coefficients for LiCs &inputph iverbosity=1 tr2_ph=1.0d-12, prefix='mgb2', fildvscf='mgb2dv', amass(1)=24.305, amass(2)=10.811, outdir='$TMP_DIR/', fildyn='mgb2.dyn', electron_phonon='interpolated', trans=.true., ldisp=.true. nq1=8, nq2=8, nq3=8 / EOF $ECHO " running the el-ph calculation...\c" $PH_COMMAND < mgb2.elph.in > mgb2.elph.out check_failure $? $ECHO " done" # # q2r and matdyn # cat > q2r.in << EOF &input zasr='simple', fildyn='mgb2.dyn', flfrc='mgb2.fc', la2F=.true. / EOF $ECHO " running q2r...\c" $Q2R_COMMAND < q2r.in > q2r.out check_failure $? $ECHO " done" # # # cat > matdyn.in.freq << EOF &input asr='simple', amass(1)=24.305, amass(2)=10.811, flfrc='mgb2.fc', flfrq='mgb2.freq', la2F=.true., dos=.false. / 19 0.000 0.0 0.0 0.0 0.125 0.0 0.0 0.0 0.250 0.0 0.0 0.0 0.375 0.0 0.0 0.0 0.500 0.0 0.0 0.0 0.750 0.0 0.0 0.0 1.000 0.0 0.0 0.0 0.825 0.125 0.125 0.0 0.750 0.250 0.250 0.0 0.625 0.375 0.375 0.0 0.500 0.500 0.500 0.0 0.325 0.325 0.325 0.0 0.250 0.250 0.250 0.0 0.125 0.125 0.125 0.0 0.000 0.000 0.000 0.0 0.125 0.125 0.000 0.0 0.250 0.250 0.000 0.0 0.325 0.325 0.000 0.0 0.500 0.500 0.000 0.0 EOF $ECHO " running matdyn for frequency calculation...\c" $MATDYN_COMMAND < matdyn.in.freq > matdyn.out.freq check_failure $? $ECHO " done" # # # cat > matdyn.in.dos << EOF &input asr='simple', amass(1)=24.305, amass(2)=10.811, flfrc='mgb2.fc', flfrq='mgb2.freq', la2F=.true., dos=.true. fldos='phonon.dos', nk1=10, nk2=10, nk3=10, ndos=50 / EOF $ECHO " running matdyn for a2F(omega) calculation...\c" $MATDYN_COMMAND < matdyn.in.dos > matdyn.out.dos check_failure $? $ECHO " done" # # # cat > lambda.in << EOF 50 0.05 1 ! emax (something more than highest phonon mode in THz), degauss, smearing method 1 ! Number of q-points for which EPC is calculated, /其余的q點會在計算Tc是加上。 0.000000 0.000000 0.000000 1 elph. 0.000000. 0.000000. 0.000000 0.10 ! \mu the Coloumb coefficient in the modified ! Allen-Dynes formula for T_c (via \omega_log) EOF $ECHO " running lambda.x for lambda calculation...\c" $LAMBDA_COMMAND < lambda.in > lambda.out check_failure $? $ECHO " done" |
回復(fù)此樓» 猜你喜歡
315求調(diào)劑
已經(jīng)有8人回復(fù)
-
312求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿鄭州大學(xué)材料與化工085600，求調(diào)劑
已經(jīng)有4人回復(fù)
-
315分 085602 求調(diào)劑
已經(jīng)有7人回復(fù)
-
298分 070300求調(diào)劑
已經(jīng)有5人回復(fù)
-
349求調(diào)劑
已經(jīng)有9人回復(fù)
-
302求調(diào)劑
已經(jīng)有4人回復(fù)
-
化工調(diào)劑303分，過四級
已經(jīng)有10人回復(fù)
-
321求調(diào)劑
已經(jīng)有14人回復(fù)
-
296材料專碩求調(diào)劑
已經(jīng)有14人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 超導(dǎo)元素臨界溫度與電場，磁場的關(guān)系（即公式）
已經(jīng)有7人回復(fù)
- 超導(dǎo)元素（單質(zhì)）中的同位素有幾種，臨界溫度是多少？
已經(jīng)有3人回復(fù)
- 【求助】超導(dǎo)的原理是什么？
已經(jīng)有4人回復(fù)
ccmin4
銀蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 631.4
- 散金: 30
- 沙發(fā): 1
- 帖子: 60
- 在線: 53.6小時
- 蟲號: 2974959
- 注冊: 2014-02-18
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2014-05-28 13:47:00
3樓2014-06-24 07:32:39
Robin122R
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 22
- 帖子: 7
- 在線: 8小時
- 蟲號: 10696870
- 注冊: 2018-10-15
- 專業(yè): 固體力學(xué)
4樓2020-09-25 17:02:08
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 flying_lw 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +15 | c297914 2026-03-28 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +15 | 太容易1018 2026-04-01 | 15/750 |
|
|
|
[考研] 292分，材料與化工，申請調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
|
[考研] 0805求調(diào)劑 +8 | 是水分 2026-03-31 | 8/400 |
|
|
|
[基金申請] esi高被引論文是不是能對中標(biāo)有所加分和幫助呢 +4 | redcom 2026-04-01 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +12 | Amber00 2026-03-31 | 12/600 |
|
|
|
[考博] 26年申博 +3 | staryer 2026-03-30 | 4/200 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
|
[考研] 省雙一流重點一本大學(xué)招收調(diào)劑 +4 | wwwwffffff 2026-03-31 | 7/350 |
|
|
|
[考研] 262求調(diào)劑 +9 | 勵志一定發(fā)文章 2026-03-31 | 10/500 |
|
|
|
[考研] 0855機械初試280求調(diào)劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 雙馬尾痞老板2 2026-03-31 | 7/350 |
|
|
|
[考研] 080500-315分復(fù)試調(diào)劑 +9 | 上岸3821 2026-03-31 | 9/450 |
|
|
|
[考研] 英一數(shù)一總分334求調(diào)劑 +4 | 陳陽坤 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿浙江大學(xué)工科動力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|