# 'Full' DFPT calculation of third derivatives using Second Order Sternheimer equation
# Test on AlAs, with norm conserving pseudopotentials
# (L. Baguet, 04.2018)

# 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

# SCF procedure
#*******************************
   nstep    100

# Bands
#*******************************
  occopt    1
   nband    4

# 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-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
    rfdir3   1 1 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
  rfatpol4   1 2
    rfdir4   1 1 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_atpol7   1 2
  d3e_pert1_elfd7   1
   d3e_pert1_dir7   1 1 1
  d3e_pert2_elfd7   1
   d3e_pert2_dir7   1 1 1
  d3e_pert3_elfd7   1
   d3e_pert3_dir7   1 1 1

 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosTM_pwteter/13al.pspnc, PseudosTM_pwteter/33as.pspnc"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%    t81.out, 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).
#%%<END TEST_INFO>