| 查看: 643 | 回复: 0 | ||
[求助]
abinit 如何用剪刀差修正静态非线性响应？
|
|
我现在正用abinit教程中的lesson on static non-linear properties计算晶体的非线性响应（非线性光学系数），计算结果和实验值有点差别，看到参考文献上有用剪刀修正的。但是sciss这个参数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计算应该分别计算。请问我这个输入文件中怎样引入剪刀差？是否需要大幅度改动输入文件？ 以下是我的 .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 |
回复此楼» 猜你喜欢
不应该看fileCode
已经有6人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有7人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
关于Filecode分析方法
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
