mirror of https://github.com/abinit/abinit.git
128 lines
3.3 KiB
Plaintext
128 lines
3.3 KiB
Plaintext
# 'Full' DFPT calculation of third derivatives using Second Order Sternheimer equation
|
|
# Test on AlAs, with PAW pseudopotentials, using OpenMP GPU.
|
|
# (M. Sarraute, 01.2024, variant of v8/t83)
|
|
|
|
# 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)
|
|
|
|
# Enable OpenMP GPU
|
|
#*********************************
|
|
gpu_option 2
|
|
|
|
# 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
|
|
bandpp 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 4
|
|
jdtset 1 2 3 4
|
|
|
|
# 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-20
|
|
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
|
|
tolwfr2 1.0d-22
|
|
|
|
# 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
|
|
|
|
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 = t21.abi
|
|
#%% [files]
|
|
#%% files_to_test =
|
|
#%% t21.abo, tolnlines =0 , tolabs =1.0E-8, tolrel = 1.0E-8, fld_options = -easy;
|
|
#%% [paral_info]
|
|
#%% max_nprocs = 2
|
|
#%% [extra_info]
|
|
#%% authors = M. Sarraute, L. Baguet
|
|
#%% keywords = PAW, DFPT
|
|
#%% description =
|
|
#%% Preparation for a nonlinear calculation : GS + linear response (ddk+phonon+efield), using GEMM nonlop.
|
|
#%% Related tests: v8 83 and gpu_omp t25
|
|
#%%<END TEST_INFO>
|