mirror of https://github.com/abinit/abinit.git
110 lines
4.1 KiB
Plaintext
110 lines
4.1 KiB
Plaintext
# NaI pSIC-PBE - 8 atom cell. Too small number of atoms, still a polaron forms with an expanded lattice parameter.
|
|
|
|
#Parameters for the image dynamics
|
|
imgmov 6
|
|
nimage 3
|
|
#ntimimage 30 # This value will give a converged result
|
|
ntimimage 6 # This value is simply to start the optimization, and show the reduction of forces of about a factor 5.
|
|
tolimg 5.0d-6 # Note that for the time being, this is the only stopping criterion for the itimimage loop, except for ntimimage.
|
|
# It would be better to implement another stopping criterion ...
|
|
prtvol 1
|
|
imgwfstor 0
|
|
mixesimgf -59.00 +80.0 -20.0 # Must be coherent with the occupation numbers
|
|
|
|
# Bands and Occupations
|
|
nsppol 2
|
|
nband 32
|
|
occopt 0
|
|
|
|
occ_1img 32*1 32*1
|
|
cellcharge_1img 0.0
|
|
|
|
occ_2img 31*1 0.975 32*1
|
|
cellcharge_2img 0.025
|
|
|
|
occ_3img 31*1 0.950 32*1
|
|
cellcharge_3img 0.050
|
|
|
|
# Define the parallelism distribution explicitly, on order to make the user conscious that this is
|
|
# a medium-grained parallelism.
|
|
# autoparal 1
|
|
paral_kgb 1
|
|
npimage 3
|
|
np_spkpt 4
|
|
npband 2
|
|
bandpp 1 # maybe not the most efficient choice
|
|
|
|
###########################################################################
|
|
#Common to all images
|
|
###########################################################################
|
|
|
|
# Optimization parameters
|
|
ionmov 2
|
|
tolmxf 1.0d-5 # This criterion is irrelevant in the present status of the code.
|
|
|
|
# SCF procedure
|
|
nstep 20
|
|
ecut 17
|
|
pawecutdg 40
|
|
toldff 1.0d-6 # For geometry optimization of linear combination of images, the criterion usually preferred, based on tolrff
|
|
# cannot be used straighforwardly, because each image will have non zero forces,
|
|
# and only the linear combination of forces gives
|
|
# Added in v9.11:
|
|
tolwfr_diago 1d-30
|
|
# Default of tolwfr_diago is tolwfr (or, for LOBPCG only, 1d-20 if tolwfr is not defined)
|
|
# To reproduce old behaviour : set stringent value of tolwfr_diago
|
|
|
|
diemac 2.0d0
|
|
|
|
# k-point grid # Not for production. Simply, to yield one k point thanks to symmetries.
|
|
kptopt 1
|
|
ngkpt 2 2 2
|
|
|
|
#Definition of the unit cell
|
|
#acell 3*12.227544 # This is the optimized PAW PBE lattice parameter
|
|
acell 3*13 # Enlarged to stabilize the polaron in the pSICn but not stabilized with cellcharge 1.0
|
|
|
|
#Definition of the atom types and unit cell
|
|
ntypat 2
|
|
znucl 11 53
|
|
pp_dirpath "$ABI_PSPDIR/Psdj_paw_pbe_std"
|
|
pseudos "Na.xml, I.xml"
|
|
|
|
#Definition of atoms and positions
|
|
natom 8
|
|
typat 4*1 4*2
|
|
|
|
xred -3.0863689067E-02 -3.0863689067E-02 2.5197195346E-31
|
|
1.9850978742E-02 4.8014902126E-01 5.0000000000E-01
|
|
4.8014902126E-01 1.9850978742E-02 5.0000000000E-01
|
|
5.3086368907E-01 5.3086368907E-01 -1.5606781514E-30
|
|
4.2000000000E-01 8.0000000000E-02 0.0000000000E+00
|
|
8.0000000000E-02 4.2000000000E-01 0.0000000000E+00
|
|
-2.1313787448E-02 -2.1313787448E-02 5.0000000000E-01
|
|
5.2131378745E-01 5.2131378745E-01 5.0000000000E-01
|
|
|
|
#%%<BEGIN TEST_INFO>
|
|
#%% [setup]
|
|
#%% executable = abinit
|
|
#%% [files]
|
|
#%% [paral_info]
|
|
#%% max_nprocs = 24
|
|
#%% nprocs_to_test = 24
|
|
#%% [NCPU_24]
|
|
#%% files_to_test =
|
|
#%% tpsic_03_MPI24.abo, tolnlines = 500, tolabs = 1.081e-05, tolrel = 0.001
|
|
#%% [extra_info]
|
|
#%% authors = X. Gonze, C. Tantardini
|
|
#%% keywords = pSIC
|
|
#%% description =
|
|
#%% NaI with a sligthly expanded lattice, to allow self-localisation of a hole polaron, using pSIC technique.
|
|
#%% Realistic ecut and reasonable k point sampling, but small cell of only 8 atoms for testing.
|
|
#%% The present test starts from a reasonable I-I interatomic distance, and search for the optimized polaron geometry, without constraint.
|
|
#%%
|
|
#%% Parallelism until 12 procs is quite efficient. However, setting npband 2 does not bring noticeable improvement.
|
|
#%% Although this is still subject for further investigation, the hypothesis is that the density/xc/potential part
|
|
#%% is not improved with npband 2, while it is actually a major contribution, even with 32 bands.
|
|
#%% Timing analysis is to be performed.
|
|
#%% topics = CrossingBarriers
|
|
#%%<END TEST_INFO>
|