 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4896 | 回復(fù): 3 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
MgB2超導(dǎo)轉(zhuǎn)變溫度
|
||
|
最近在學(xué)習(xí)用quantum espresso計算超導(dǎo)轉(zhuǎn)變溫度，所以就用比較成熟的MgB2體系作為例子進行實際的計算。但是計算結(jié)果與實驗值差距很大，實驗值是39.3 K，而我最高只能算到4.1 K，想盡各種辦法也沒能有效提高，不知道是參數(shù)設(shè)置的不對還是贗勢選擇有問題�？吹接行┪墨I用QE計算的MgB2的轉(zhuǎn)變溫度與實驗值復(fù)合的很好。但是自己計算的就是不對，有這方面經(jīng)驗的蟲友請幫忙看一下輸入文件，有什么地方是可以改進或是改正的。謝謝了！ 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ù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有6人回復(fù)
-
280求調(diào)劑
已經(jīng)有6人回復(fù)
-
302求調(diào)劑一志愿華中師范大學(xué)
已經(jīng)有5人回復(fù)
-
346求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專碩322分
已經(jīng)有9人回復(fù)
-
372分材料與化工（085600）一志愿湖南大學(xué)求調(diào)劑
已經(jīng)有5人回復(fù)
-
318求調(diào)劑，計算材料方向
已經(jīng)有10人回復(fù)
-
310求調(diào)劑
已經(jīng)有14人回復(fù)
-
315求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有13人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 超導(dǎo)元素臨界溫度與電場，磁場的關(guān)系（即公式）
已經(jīng)有7人回復(fù)
- 超導(dǎo)元素（單質(zhì)）中的同位素有幾種，臨界溫度是多少？
已經(jīng)有3人回復(fù)
- 【求助】超導(dǎo)的原理是什么？
已經(jīng)有4人回復(fù)
Robin122R
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 22
- 帖子: 7
- 在線: 8小時
- 蟲號: 10696870
- 注冊: 2018-10-15
- 專業(yè): 固體力學(xué)
4樓2020-09-25 17:02:08
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 071000生物學(xué)調(diào)劑 +6 | 知昭蔓 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 總分328生物與醫(yī)藥考數(shù)學(xué)求調(diào)劑 +3 | aaadim 2026-04-02 | 3/150 |
|
|
|
[考研] 求調(diào)劑推薦 +3 | 南山南@ 2026-04-01 | 3/150 |
|
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +8 | pipiver 2026-03-30 | 8/400 |
|
|
[考研] 324分 085600材料與化工 +20 | 呆鵝oor 2026-03-27 | 20/1000 |
|
|
|
[考研] 290求調(diào)劑 +5 | dfffsar 2026-03-29 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +25 | lllllos 2026-03-30 | 26/1300 |
|
|
|
[考研] 求調(diào)劑：一志愿：南京大學(xué) 專業(yè)：0705 總分320 ，本科985，四六級已過 +3 | lfy760306 2026-03-31 | 3/150 |
|
|
|
[考研] 318求調(diào)劑 +3 | 篤行致遠. 2026-03-31 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 雙馬尾痞老板2 2026-03-31 | 7/350 |
|
|
|
[考研] 材料工程085601數(shù)二英一335求調(diào)劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|