mirror of https://github.com/abinit/abinit.git
156 lines
4.1 KiB
Plaintext
156 lines
4.1 KiB
Plaintext
# 'Full' DFPT calculation of third derivatives using Second Order Sternheimer equation
|
|
# Test on AlAs, with norm conserving pseudopotentials
|
|
# (L. Baguet, 05.2018)
|
|
# nsppol = 2
|
|
|
|
# Enable output for nonlinear (full DFPT only)
|
|
#*****************************************************
|
|
# nonlinear_info 1 # print details of 3rd derivatives in .out file (no time consuming)
|
|
# nonlinear_info 2 # nonlinear_info=1 + debug_mode activated in nonlinear (time consuming)
|
|
# nonlinear_info 3 # nonlinear_info=1 + debug_mode activated in rf2_init (time consuming)
|
|
# nonlinear_info 4 # nonlinear_info=1 + debug_mode activated in both nonlinear and rf2_init (time consuming)
|
|
|
|
# Elementary cell
|
|
#*********************************
|
|
acell 3*10.64
|
|
rprim 0.0 0.5 0.5
|
|
0.5 0.0 0.5
|
|
0.5 0.5 0.0
|
|
|
|
# Atoms
|
|
#***********************
|
|
natom 2
|
|
ntypat 2
|
|
znucl 13 33
|
|
typat 1 2
|
|
|
|
xred 0.00 0.00 0.00
|
|
0.25 0.25 0.25
|
|
|
|
# Exchange correlation functional
|
|
#********************************
|
|
ixc 7
|
|
|
|
# Polarization
|
|
#*******************************
|
|
nsppol 2
|
|
spinmagntarget 0
|
|
|
|
# SCF procedure
|
|
#*******************************
|
|
nstep 100
|
|
|
|
# Bands
|
|
#*******************************
|
|
occopt 1
|
|
nband 8
|
|
|
|
# Plane wave basis set
|
|
#**************************************
|
|
ecut 3
|
|
ecutsm 0
|
|
|
|
# K point grid
|
|
#**************************************
|
|
ngkpt 2 2 2
|
|
nshiftk 4
|
|
shiftk 0.5 0.5 0.5
|
|
0.5 0.0 0.0
|
|
0.0 0.5 0.0
|
|
0.0 0.0 0.5
|
|
|
|
# Number of Datasets
|
|
#**************************************
|
|
ndtset 7
|
|
|
|
# For all datasets (except 1 and 2)
|
|
#*******************************
|
|
kptopt 2
|
|
tolwfr 1.0d-26
|
|
getden 2
|
|
getwfk 2
|
|
|
|
# DATASET1 : Ground state (density)
|
|
#*******************************
|
|
getden1 0
|
|
getwfk1 0
|
|
tolvrs1 1.0d-8
|
|
kptopt1 1
|
|
|
|
# DATASET2 : Ground state (highly converged wavefunction)
|
|
#*******************************
|
|
getden2 1
|
|
getwfk2 1
|
|
|
|
# DATASET3 : ddk (SCF cycles are useless)
|
|
#*******************************
|
|
rfddk3 1
|
|
# For a more effective non self-consistent computation:
|
|
nstep3 1
|
|
nline3 100
|
|
tolrde3 1.0d-30 # tolrde is choosen to be much lower than tolwfr.
|
|
# This way the conjugate gradient steps stop at tolwfr, and not tolrde (usually around 1.0d-3).
|
|
# If nline is sufficiently large, the computation converges in one step only.
|
|
|
|
# DATASET4 : Phonons, Electric field
|
|
#*******************************
|
|
rfelfd4 3
|
|
rfphon4 1
|
|
getddk4 3
|
|
prtden4 1
|
|
prepanl4 1
|
|
|
|
# DATASET5 : dkdk
|
|
#*******************************
|
|
rf2_dkdk5 1
|
|
getddk5 3
|
|
prepanl5 1
|
|
# Non self-consistent computation : can converge in one step. If not, increase nline instead of nstep. tolrde is not used.
|
|
nstep5 1
|
|
nline5 100
|
|
|
|
# DATASET6 : dkde
|
|
#*******************************
|
|
rf2_dkde6 1
|
|
getddk6 3
|
|
get1den6 4
|
|
getdelfd6 4
|
|
getdkdk6 5
|
|
prepanl6 1
|
|
# Non self-consistent computation : can converge in one step. If not, increase nline instead of nstep. tolrde is not used.
|
|
nstep6 1
|
|
nline6 100
|
|
|
|
# DATASET7 : 3DTE calculation (full DFPT)
|
|
#*****************************************
|
|
optdriver7 5 # for nonlinear calculation
|
|
usepead7 0
|
|
getddk7 3
|
|
get1den7 4
|
|
get1wf7 4
|
|
getdkde7 6
|
|
d3e_pert1_phon7 1
|
|
d3e_pert1_elfd7 1
|
|
d3e_pert2_elfd7 1
|
|
d3e_pert3_elfd7 1
|
|
|
|
pp_dirpath "$ABI_PSPDIR"
|
|
pseudos "PseudosTM_pwteter/13al.pspnc, PseudosTM_pwteter/33as.pspnc"
|
|
|
|
#%%<BEGIN TEST_INFO>
|
|
#%% [setup]
|
|
#%% executable = abinit
|
|
#%% [files]
|
|
#%% files_to_test =
|
|
#%% t82.abo, tolnlines =1 , tolabs =1.0E-8, tolrel = 1.0E-8, fld_options = -medium;
|
|
#%% [paral_info]
|
|
#%% max_nprocs = 10
|
|
#%% [extra_info]
|
|
#%% authors = L. Baguet
|
|
#%% keywords = NC, DFPT, NONLINEAR
|
|
#%% description =
|
|
#%% 'Full' DFPT computation of third derivatives in Nonlinear (dataset 7).
|
|
#%% Preceded by resolution of Second-order Sternheimer equations (dataset 5 and 6).
|
|
#%% Test the polarized case (with a non-polarized system). Must give same results than t81.in
|
|
#%%<END TEST_INFO>
|