| 1 | 1/1 | 返回列表 |
| 查看: 585 | 回復(fù): 0 | ||
[求助]
abinit 如何用剪刀差修正靜態(tài)非線性響應(yīng)?
|
|
我現(xiàn)在正用abinit教程中的lesson on static non-linear properties計算晶體的非線性響應(yīng)(非線性光學(xué)系數(shù)),計算結(jié)果和實驗值有點差別,看到參考文獻上有用剪刀修正的。但是sciss這個參數(shù)typical use is for response to electric field (rfelfd=3), but NOT for d/dk (rfelfd=2) and phonon responses,。在FAQ中也說在剪刀修正sciss引入的時候,rfphon 1 和 rfelfd 1計算應(yīng)該分別計算。請問我這個輸入文件中怎樣引入剪刀差?是否需要大幅度改動輸入文件? 以下是我的 .in文件 ixc 3 #Parameters of the SCF cycles #**************************** iscf 5 nstep 100 #Plane wave basis and k-point grid #********************************* ecut 56 ecutsm 0.5 dilatmx 1.05 ngkpt 7 7 7 nshiftk 1 shiftk 0 0 0 ndtset 5 jdtset 1 2 3 4 5 #DATASET1 : scf calculation: GS WF in the BZ #******************************************** prtden1 1 kptopt1 1 toldfe1 1.0d-12 #DATASET2 : non scf calculation: GS WF in the whole BZ #***************************************************** getden2 1 kptopt2 2 iscf2 -2 getwfk2 1 tolwfr2 1.0d-22 nband2 26 #DATASET3 : ddk on the MP grid of k-points #***************************************** getwfk3 2 rfdir3 1 1 1 rfelfd3 2 tolwfr3 1.0d-22 nband3 26 kptopt3 2 #DATASET4 : ddE on the MP grid of k-points #***************************************** prtden4 1 getwfk4 2 getddk4 3 rfdir4 1 1 1 rfphon4 1 rfatpol4 1 8 tolvrs4 1.0d-12 rfelfd4 3 nband4 26 kptopt4 2 prepanl4 1 #DATASET5 : 3DTE calculation #*************************** getden5 1 get1den5 4 getwfk5 2 get1wf5 4 nband5 26 kptopt5 2 optdriver5 5 rf1elfd5 1 rf1phon5 1 rf1atpol5 1 8 rf1dir5 1 1 1 rf2elfd5 1 rf2dir5 1 1 1 rf3elfd5 1 rf3dir5 1 1 1 nbdbuf 0 #define the crystaline ntypat 3 znucl 30.00000 32.00000 15.00000 natom 8 acell 9.8122000262E-01 9.8122000262E-01 9.8122000262E-01 Bohr rprim -5.1636766581E+00 5.1636766581E+00 1.0118349606E+01 5.1636766581E+00 -5.1636766581E+00 1.0118349606E+01 5.1636766581E+00 5.1636766581E+00 -1.0118349606E+01 typat 1 2 3 3 1 2 3 3 scalecart 3*1 xred 5.0000000000E-01 5.0000000000E-01 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 3.7500000000E-01 3.6780769517E-01 4.9280769517E-01 1.1780769517E-01 6.2500000000E-01 9.9280769517E-01 2.5000000000E-01 7.5000000000E-01 5.0000000000E-01 7.5000000000E-01 2.5000000000E-01 5.0000000000E-01 8.7500000000E-01 8.8219230483E-01 5.0719230483E-01 6.3219230483E-01 1.2500000000E-01 7.1923048258E-03 |

| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +5 | 一定上岸aaa 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 303求調(diào)劑 +7 | 安憶靈 2026-03-22 | 8/400 |
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
[考研] 機械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|