| 查看: 608 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
论文撤稿了
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
河北省自然科学基金
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
10