mirror of https://github.com/abinit/abinit.git
157 lines
4.1 KiB
Plaintext
157 lines
4.1 KiB
Plaintext
# 'Full' DFPT calculation of third derivatives using Second Order Sternheimer equation
|
|
# Test on AlAs, with PAW pseudopotentials
|
|
# (L. Baguet, 04.2018)
|
|
|
|
# Enable ouput 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
|
|
|
|
# SCF procedure
|
|
#*******************************
|
|
nstep 100
|
|
|
|
# Bands
|
|
#*******************************
|
|
occopt 1
|
|
nband 4
|
|
|
|
# Plane wave basis set
|
|
#**************************************
|
|
ecut 3.5
|
|
ecutsm 0
|
|
pawecutdg 7.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
|
|
ndtset 3
|
|
jdtset 5 6 7
|
|
|
|
# PAW option
|
|
#*******************************
|
|
pawxcdev 0 # non-zero pawxcdev is not allowed for dataset 7, so we use pawxcdev=0 for all
|
|
|
|
# For all datasets (except 1 and 2)
|
|
#*******************************
|
|
kptopt 2
|
|
tolwfr 1.0d-28
|
|
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 "Psdj_paw_pw_std/Al.xml, As.LDA_PW-JTH_sp.xml"
|
|
|
|
#%%<BEGIN TEST_INFO>
|
|
#%% [setup]
|
|
#%% executable = abinit
|
|
#%% test_chain = t83.abi, t84.abi
|
|
#%% [files]
|
|
#%% files_to_test =
|
|
#%% t84.abo, tolnlines =0 , tolabs =1.0E-8, tolrel = 1.0E-8, fld_options = -medium;
|
|
#%% [paral_info]
|
|
#%% max_nprocs = 10
|
|
#%% [extra_info]
|
|
#%% authors = L. Baguet
|
|
#%% keywords = PAW, 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).
|
|
#%%<END TEST_INFO>
|