24小时热门版块排行榜

返回列表

肖瑞春

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 293.1
红花: 1
帖子: 53
在线: 45.5小时
虫号: 2407422
注册: 2013-04-08
专业: 凝聚态物性 II ：电子结构

[求助] 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

回复此楼

» 猜你喜欢

快乐生活每一天！

1楼 2014-04-10 11:23:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主肖瑞春的主题更新

返回列表