gstate input variables¶
This document lists and provides the description of the name (keywords) of the gstate input variables to be used in the input file for the abinit executable.
algalch¶
Mnemonics: ALGorithm for generating ALCHemical pseudopotentials
Mentioned in topic(s): topic_AtomTypes
Variable type: integer
Dimensions: (ntypalch)
Default value: ntypalch * 1
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v3: t92.abi
Used for the generation of alchemical pseudopotentials, that is, when ntypalch is nonzero.
Give the algorithm to be used to generate the ntypalch alchemical potentials from the different %npspalch pseudopotentials dedicated to this use.
Presently, algalch can only have the value 1, that is:
 the local potentials are mixed, thanks to the mixalch mixing coefficients
 the form factors of the nonlocal projectors are all preserved, and all considered to generate the alchemical potential
 the scalar coefficients of the nonlocal projectors are multiplied by the proportion of the corresponding type of atom that is present in mixalch
 the characteristic radius for the core charge is a linear combination of the characteristic radii of the core charges, build with the mixalch mixing coefficients
 the core charge function f(r/r_c) is a linear combination of the core charge functions, build with the mixalch mixing coefficients
Later, other algorithms for the mixing might be included.
Important
Note that alchemical mixing cannot be used with PAW.
auxc_ixc¶
Mnemonics: AUxiliary XC functional for hybrid functional, IXC number
Mentioned in topic(s): topic_Hybrids
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v7: t67.abi
Specification of an auxiliary exchangecorrelation functional, thanks to its ixc value, to possibly replace the heavy evaluation of an hybrid functional at specific occasions, e.g. when the Fock operator is frozen during the selfconsistent cycle, thanks to fockoptmix == 11, or when evaluating the correction to forces due to the density residual. This auxiliary exchange correlation functional might be rescaled, thanks to auxc_scal when fockoptmix == 11. If gwcalctyp == 5, 15 or 25, auxc_ixc refers to ixc_sigma instead of ixc.
auxc_scal¶
Mnemonics: AUxiliary XC functional for hybrid functional SCALing factor
Mentioned in topic(s): topic_Hybrids
Variable type: real
Dimensions: scalar
Default value: 1.0
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v7: t67.abi
Possible scaling factor for the auxiliary exchangecorrelation functional defined by auxc_ixc that has the goal to replace the Fock operator or hybrid functional when fockoptmix == 11.
The default value 1.0 corresponds to the unmodified xc functional. When the auxiliary functional is used to replace the hybrid functional in SCF loops, a value of 1.5 has been observed to speed up somehow the convergence.
boxcenter¶
Mnemonics: BOX CENTER
Mentioned in topic(s): topic_TDDFT
Variable type: real
Dimensions: (3)
Default value: [0.5, 0.5, 0.5]
Added in version: before_v9
Test list (click to open). Rarely used, [11/1101] in all abinit tests, [1/151] in abinit tutorials
Defines the center of the box, in reduced coordinates. At present, this information is only used in the case of TimeDependent DFT computation of the oscillator strength. One must take boxcenter such as to be roughly the center of the cluster or molecule. The default is sensible when the vacuum surrounding the cluster or molecule has xred 0 or 1. On the contrary, when the cluster or molecule is close to the origin, it is better to take boxcenter = [0.0, 0.0, 0.0].
boxcutmin¶
Mnemonics: BOX CUToff MINimum
Mentioned in topic(s): topic_Planewaves, topic_TuningSpeedMem
Variable type: real
Dimensions: scalar
Default value: 2.0
Added in version: before_v9
Test list (click to open). Moderately used, [19/1101] in all abinit tests, [8/151] in abinit tutorials
 mpiio: t42.abi
 tutorespfn: teph4mob_5.abi, teph4mob_6.abi, teph4mob_7.abi, teph4isotc_2.abi, teph4isotc_5.abi, teph4zpr_6.abi, teph4zpr_7.abi, teph4zpr_8.abi
 v3: t43.abi
 v5: t75.abi, t76.abi
 v7: t95.abi
 v8: t02.abi, t17.abi
 v9: t29.abi, t57.abi, t60.abi, t61.abi
The box cutoff ratio is the ratio between the radius of the Gsphere that can be inserted in the FFT box (see ngfft) and the Gsphere used to represent the wavefunctions as computed from the input ecut. In order for the density to be exact (in the case of the plane wave part, not the PAW onsite terms), this ratio should be at least two.
Warning
If one uses a smaller ratio (e.g. 1.5), one will gain speed at the price of accuracy. Keep in mind that using a value of boxcutmin too close to 1 may lead to runtime errors in the FFT routines so a value larger than 1.1 is strongly recommended.
It should also be stressed that the quality of the forces is affected by the value of boxcutmin . The default value (2.0) is OK for ground state calculations and structural relaxations. On the contrary, the computations of vibrational frequencies based on finitedifference methods such as those implemented by phonopy turn out to be rather sensitive to the quality of the forces. Tests have shown that, for particular systems, phonopy calculations require a larger boxcutmin (e.g. 3.0) to reproduce the DFPT results. The optimal value of boxcutmin likely depends on the pseudopotentials and the system under investigation, yet this is something worth keeping in mind.
Note
Prior to v8.9, the use of boxcutmin for DFPT calculations was forbidden. However, after testing, it was seen that the deterioration in phonon band structures could be alleviated to a large extent by the imposition of the Acoustic Sum Rule asr.
cellcharge¶
Mnemonics: CELL CHARGE
Mentioned in topic(s): topic_Coulomb
Variable type: real
Dimensions: (nimage)
Default value: 0
*Added in version: v9.3.4
Test list (click to open). Moderately used, [13/1101] in all abinit tests, [7/151] in abinit tutorials
 tutoparal: timages_01.abi, timages_02.abi, timages_03.abi, timages_04.abi, tpsic_01.abi, tpsic_02.abi, tpsic_03.abi
 v2: t42.abi
 v3: t54.abi
 v5: t42.abi
 v6: t67.abi
 v7: t26.abi
 v9: t22.abi
Used to establish charge balance between the number of electrons filling the bands and the nominal charge associated with the atomic cores. The code adds up the number of valence electrons provided by the pseudopotentials of each type (call this “zval”), then add cellcharge , to get the number of electrons per unit cell, %nelect. Then, if iscf is positive, the code adds up the band occupancies (given in array occ) for all bands at each k point, then multiplies by the k point weight wtk at each k point. Call this sum “nelect_occ” (for the number of electrons from occupation numbers). It is then required that: nelect_occ = %nelect. To treat a neutral system, which is desired in nearly all cases, one must use cellcharge = 0. To treat a system missing one electron per unit cell, set cellcharge = +1.
cellcharge superceeds the old charge input variable, whose name was rather unspecific.
When there are several images, cellcharge might depend on the image number, but ONLY when imgmov=6 and occopt=0 or 2. In the checking routine, %nelect is considered separately for each image, while in the remaining of the code, %nelect(1) is propagated, so that %nelect is still a scalar. This is consistent with the pSIC algorithm, see [Sadigh2015] and [Sadigh2015a].
charge¶
Mnemonics: CHARGE
Mentioned in topic(s): topic_Coulomb
Variable type: real
Dimensions: scalar
Default value: 0
Added in version: before_v9, obsolete
Test list (click to open). Moderately used, [64/1101] in all abinit tests, [12/151] in abinit tutorials
 mpiio: t62.abi, t62.abi …
 paral: t07.abi, t07.abi, t07.abi …
 tutoparal: tdmft_2.abi, timages_01.abi, timages_02.abi …
 tutorespfn: telast_2.abi, tffield_1.abi …
 tutorial: tbase1_4.abi, tfold2bloch_1.abi …
 v2: t06.abi, t30.abi, t40.abi …
 v3: t54.abi, t77.abi …
 v4: t02.abi, t08.abi, t46.abi …
 v5: t06.abi, t11.abi, t23.abi …
 v6: t42.abi, t43.abi, t67.abi …
 v7: t09.abi, t10.abi, t17.abi …
 v8: t24.abi, t25.abi, t26.abi …
 v9: t01.abi, t02.abi, t03.abi …
The charge input variable is obsolete, and has been replaced by cellcharge. It is still read during a transitional period, likely up to the end of v9 life cycle.
chkexit¶
Mnemonics: CHecK whether the user want to EXIT
Mentioned in topic(s): topic_Control
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [0/151] in abinit tutorials
If chkexit is 1 or 2, ABINIT will check whether the user wants to interrupt the run (using the keyword “exit” on the top of the input file or creating a file named “abinit.exit”: see the end of section 3.2 of the abinit help file).
 0 → the check is not performed at all
 1 → the check is not performed frequently (after each SCF step)
 2 → the check is performed frequently (after a few bands, at each k point)
In all cases, the check is performed at most every 2 seconds of CPU time.
chkprim¶
Mnemonics: CHecK whether the cell is PRIMitive
Mentioned in topic(s): topic_crystal, topic_UnitCell, topic_SmartSymm
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Moderately used, [65/1101] in all abinit tests, [8/151] in abinit tutorials
 bigdft: t00.abi, t20.abi, t21.abi …
 builtin: testin_bigdft.abi …
 libxc: t19.abi, t73.abi …
 mpiio: t42.abi …
 paral: t09.abi, t09.abi, t09.abi …
 tutoparal: tpsic_01.abi …
 tutorial: tbase4_4.abi, tpositron_2.abi, tpositron_3.abi …
 v1: t93.abi, t96.abi …
 v2: t55.abi, t59.abi, t60.abi …
 v3: t20.abi, t73.abi …
 v4: t95.abi, t96.abi, t99.abi …
 v5: t01.abi, t07.abi, t09.abi …
 v6: t21.abi …
 v7: t11.abi, t35.abi, t76.abi …
 v8: t71.abi, t72.abi, t73.abi …
 v9: t01.abi, t02.abi, t03.abi …
If the symmetry finder is used (see nsym), a nonzero value of chkprim will make the code stop if a nonprimitive cell is used. If chkprim = 0, a warning is issued, but the run does not stop.
If you are generating the atomic and cell geometry using spgroup, you might generate a PRIMITIVE cell using brvltt = 1.
chksymbreak¶
Mnemonics: CHecK SYMmetry BREAKing
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Moderately used, [74/1101] in all abinit tests, [8/151] in abinit tutorials
 bigdft: t07.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 etsf_io: t02.abi, t30.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t19.abi …
 paral: t08.abi, t08.abi, t08.abi …
 seq: tsv2_82.abi, tsv3_05.abi, tsv4_80.abi …
 tutorial: tbs_1.abi, tbs_2.abi, tbs_3.abi …
 v3: t62.abi …
 v4: t09.abi, t42.abi, t72.abi …
 v5: t08.abi, t16.abi, t18.abi …
 v6: t07.abi, t20.abi, t24.abi …
 v67mbpt: t31.abi, t32.abi, t33.abi …
 v7: t09.abi …
 v8: t34.abi, t35.abi, t71.abi …
 v9: t92.abi …
 wannier90: t04.abi, t11.abi …
This variable governs the behaviour of the code when there is a potential source of symmetry breaking related to the k point grid.
When chksymbreak = 1, the code stops if the k point grid is nonsymmetric, in case kptopt =1, 2, or 4. Also, the code stops if nshiftk is not 1, 2 or 4.
Note that the check is disabled when the number of kpoints in the BZ is greater than 40^3.
When chksymbreak = 0, there is no such check.
Explanation: In the groundstate calculation, such breaking of the symmetry is usually harmless. However, if the user is doing a calculation of phonons using DFPT (rfphon = 1), the convergence with respect to the number of k points will be worse with a nonsymmetric grid than with a symmetric one.
So, it was decided to warn the user about such problem already at the level of the ground state calculations, although such warning might be irrelevant.
Concerning the values of nshiftk, usage of values other then 1, 2, or 4 can hardly be understood, as it will yield in many cases a nonhomogeneous or nonsymmetric k point grid. So, it is usually an error of the user.
If you encounter a problem outlined above, you have some choices: change your k point grid, to make it more symmetric, and/or respect nshiftk=1, 2, or 4, or ignore the problem, and set chksymbreak = 0.
chksymtnons¶
Mnemonics: CHecK SYMmetry of TNONS
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_crystal
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: 9.2.0
Test list (click to open). Moderately used, [37/1101] in all abinit tests, [4/151] in abinit tutorials
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t21.abi, t22.abi, t24.abi …
 seq: tsv2_81.abi …
 tutoparal: tdfpt_03.abi, tdfpt_04.abi …
 tutorial: tbase4_7.abi, tbase4_8.abi …
 v3: t10.abi, t11.abi, t12.abi …
 v4: t01.abi, t46.abi …
 v5: t49.abi, t62.abi …
 v6: t23.abi, t34.abi …
 v7: t09.abi …
 v8: t12.abi, t30.abi …
 v9: t17.abi, t23.abi, t25.abi …
This variable governs the behaviour of the code when there is a potential symmetry breaking, related to the presence of nonsymmorphic translations not leaving the FFT exchangecorrelation grid invariant.
When chksymtnons = 1, the code stops if the nonsymmorphic translation part of the symmetry operations has components that are not zero, or simple fractions with 2, 3, 4, 5, 6, 8, 9, 10 or 12 as denominators. Also, suggestions to bypass the problem are made in the output file.
When chksymtnons = 2, the code makes similar check, but does not stop after providing in the output file suggestions to bypass the problem.
When chksymtnons = 0, the code skips the check.
Explanation: In groundstate or DFPT calculations, such breaking of the symmetry is harmless. However, for a GW or BSE calculation, the presence of nonsymmorphic translations that are not coherent with the FFT grid will cause problems (e.g. enormous memory reservation, inducing segfault). In the GW or BSE parts, indeed, one needs to reconstruct the wavefunctions in the full Brillouin zone for calculating both the polarizability and the selfenergy. The wavefunctions in the full Brillouin zone are obtained from the irreducible wedge by applying the symmetry operations of the space group of the crystal. In the present implementation, the symmetrisation of the wavefunctions is done in real space on the FFT mesh that, therefore, has to be coherent both with the rotational part as well as with the fractional translation of each symmetry operation. If the condition above (2, 3, 4, 5, 6, 7, 8, 9, 10, or 12 as denominator) is not met, the GW/BSE code will not be able to find a symmetry preserving FFT mesh.
So, it was decided to warn the user about such problem already at the level of the groundstate calculations, although such warning might be irrelevant.
If you encounter the problem outlined above, you have two choices: change your atomic positions (translate them) such that the origin appears as the most symmetric point; or ignore the problem, and set chksymtnons = 2 or 0. If chksymtnons = 2, ABINIT makes a suggestion of a possible global translation, and corresponding translated atomic positions.
chrgat¶
Mnemonics: CHARGE of the AToms
Mentioned in topic(s): topic_ConstrainedDFT
Variable type: real
Dimensions: [ natrd ] if natrd<natom,
[ natom ] otherwise.
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [7/1101] in all abinit tests, [0/151] in abinit tutorials
Gives the target integrated charge in case of constrained DFT calculations, see constraint_kind. Given in atomic unit of charge (=minus the charge of the electron). Note that this number is the net positive charge inside the sphere: one subtract from the nucleus charge %ziontypat the integrated valence electron density in a sphere defined by ratsph. The latter has indeed a negative value. Note that if the sphere radius ratsph is not sufficiently large, the amount of electrons will be smaller than expected based on chemical intuition. This means that there is in this case a bias toward too positive integrated charges. By contrast, if the sphere radius is too large, the spheres will overlap, and the electrons in the interatomic region will be double counted.
constraint_kind¶
Mnemonics: CONSTRAINT KIND in constrained DFT
Mentioned in topic(s): topic_ConstrainedDFT
Variable type: integer
Dimensions: (ntypat)
Default value: 0
Only relevant if: iscf > 1 and iscf < 10 and ionmov /= 4
Added in version: before_v9
Test list (click to open). Moderately used, [12/1101] in all abinit tests, [0/151] in abinit tutorials
If constraint_kind is nonzero for at least one type of atom, the constrained DFT algorithm is activated. constraint_kind defines, for each type of atom, the kind of constraint(s) imposed by constrained DFT. When constraint_kind is zero for an atom type, there is not constraint applied to this atom type. Otherwise, different constraints can be imposed on the total charge (ion+electronic) and/or magnetization, computed inside a sphere of radius ratsph, possibly smeared within a width ratsm. Such integrated ion+electronic charge might be imposed to be equal to chrgat, while the magnetization might be compared to spinat. The first digit of constraint_kind defines the constraint on the charge, while the second digit defines the constraint on the magnetization.
When constraint_kind is 10 or above, the charge constraint will be imposed.
When constraint_kind =1 or 11, the exact value (vector in the noncollinear case, amplitude and sign in the collinear case) of the magnetization is constrained; When constraint_kind =2 or 12, only the magnetization axis is constrained (only meaningful in the noncollinear case, albeit allowed); When constraint_kind =3 or 13, only the magnetization magnitude is constrained. When constraint_kind =4 or 14, only the magnetization direction is constrained (only meaningful in the noncollinear case, not allowed in the collinear case);
For the algorithm, see topic_ConstrainedDFT. It makes important use of the potential residual, so the algorithm works only with iscf between 2 and 9. The balance between the potential residual, and the density/magnetization constraint is governed by magcon_lambda. The spherical integral is governed by ratsph and ratsm.
Note that while a spherical integral around an atom might reasonably well capture the magnetization of an atom within a solid or within a molecule, so that the sum of such magnetizations might be reasonably close to the total magnetization of the solid, such a procedure hardly gives the total charge of the solid: the space between the spheres is too large when the spheres do not overlap, while overlapping spheres will not deliver the correct total charge of the system.
Note that constraint_kind defines constraints for types of atoms, not for specific atoms. Atoms of the same type are supposed to incur the same constraint. If the user wants to impose different constraints on atoms of the same type (in principle), it is possible (and easy) to pretend that they belong to different types, even if the same pseudopotential file is used for these atoms. There is an example in test v8[24], the hydrogen dimer, where the charge around the first atom is constrained, and the charge around the second atom is left free.
The difference between constraint_kind =4 or 14 and constraint_kind =2 or 12 lies in the fact that constraint_kind =2 or 12 will consider similarly a magnetization vector and its opposite, as this constraint is just alignment on an axis, while constraint_kind =4 or 14 enforces that the scalar product of the target magnetization direction (normalized vector) and the actual optimized magnetization direction (normalized vector) is positive.
Incidentally, ionmov==4 is not allowed in the present implementation of constrained DFT because the motion of atoms and simultaneous computation of constraints would be difficult to handle.
cpuh¶
Mnemonics: CPU time limit in Hours
Characteristics: NO_MULTI, INPUT_ONLY
Mentioned in topic(s): topic_Control
Variable type: real
Dimensions: scalar
Default value: 0.0
The use of this variable forbids the use of: cpum or cpus
Added in version: before_v9
Test list (click to open). Rarely used, [0/1101] in all abinit tests, [0/151] in abinit tutorials
Only one of the three real parameters cpus, cpum and cpuh can be defined in the input file to set up a CPU time limit. When the job reaches that limit, it will try to end smoothly. However, note that this might still take some time. If the user want a firm CPU time limit, the present parameter must be reduced sufficiently. Intuition about the actual margin to be taken into account should come with experience. A zero value has no action of the job.
Tip
One can pass the timelimit to abinit via the command line option:
abinit timelimit hours:minutes:seconds
This approach is much more powerful especially when the job must be submitted to the queue via a submission script e.g. a Slurm script. In this case, indeed, one can define a shell variable for the time limit and use this variable to pass the time limit to Slurm and Abinit at the same time.
Use abinit help
for further information.
cpum¶
Mnemonics: CPU time limit in Minutes
Characteristics: NO_MULTI, INPUT_ONLY
Mentioned in topic(s): topic_Control
Variable type: real
Dimensions: scalar
Default value: 0.0
The use of this variable forbids the use of: cpuh or cpus
Added in version: before_v9
Test list (click to open). Rarely used, [0/1101] in all abinit tests, [0/151] in abinit tutorials
Only one of the three real parameters cpus, cpum and cpuh can be defined in the input file to set up a CPU time limit. When the job reaches that limit, it will try to end smoothly. However, note that this might still take some time. If the user want a firm CPU time limit, the present parameter must be reduced sufficiently. Intuition about the actual margin to be taken into account should come with experience. A zero value has no action of the job.
cpus¶
Mnemonics: CPU time limit in seconds
Characteristics: NO_MULTI, INPUT_ONLY
Mentioned in topic(s): topic_Control
Variable type: real
Dimensions: scalar
Default value: 0.0
The use of this variable forbids the use of: cpuh or cpum
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [1/151] in abinit tutorials
 tutoparal: tgspw_01.abi
 v4: t89.abi
Only one of the three real parameters cpus , cpum and cpuh can be defined in the input file to set up a CPU time limit. When the job reaches that limit, it will try to end smoothly. However, note that this might still take some time. If the user want a firm CPU time limit, the present parameter must be reduced sufficiently. Intuition about the actual margin to be taken into account should come with experience. A zero value has no action of the job.
diecut¶
Mnemonics: DIElectric matrix energy CUToff
Characteristics: ENERGY
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: real
Dimensions: scalar
Default value: 2.2
Added in version: before_v9
Test list (click to open). Moderately used, [13/1101] in all abinit tests, [0/151] in abinit tutorials
Kinetic energy cutoff that controls the number of planewaves used to represent the dielectric matrix: (1/2) [ 2 \pi \GG_{diel,max}]^2 = diecut with \GG_{diel,max} being the maximum length of the reciprocal space planewave wavevectors for the dielectric matrix. Can be specified in Ha (the default), Ry, eV or Kelvin, since diecut has the ENERGY characteristics. (1 Ha = 27.2113845 eV) All planewaves inside this “basis sphere” centered at \GG=0 are included in the basis. This is useful only when iprcel > =21, which means that a preconditioning scheme based on the dielectric matrix is used.
NOTE: a negative diecut will define the same dielectric basis sphere as the corresponding positive value, but the FFT grid will be identical to the one used for the wavefunctions. The much smaller FFT grid, used when diecut is positive, gives exactly the same results.
No meaning for RF calculations yet.
diegap¶
Mnemonics: DIElectric matrix GAP
Characteristics: ENERGY
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: real
Dimensions: scalar
Default value: 0.1
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [0/151] in abinit tutorials
Gives a rough estimation of the dielectric gap between the highest energy level computed in the run, and the set of bands not represented. Used to extrapolate dielectric matrix when iprcel >= 21. Can be specified in Ha (the default), Ry, eV or Kelvin, since diegap has the ENERGY characteristics (1 Ha = 27.2113845 eV).
No meaning for RF calculations yet.
dielam¶
Mnemonics: DIElectric matrix LAMbda
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: real
Dimensions: scalar
Default value: 0.5
Only relevant if: iprcel >= 21
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [0/151] in abinit tutorials
Gives the amount of occupied states with mean energy given by the highest level computed in the run, included in the extrapolation of the dielectric matrix.
No meaning for RF calculations yet.
dielng¶
Mnemonics: model DIElectric screening LeNGth
Characteristics: LENGTH
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: real
Dimensions: scalar
Default value: 1.0774841
Added in version: before_v9
Test list (click to open). Moderately used, [68/1101] in all abinit tests, [28/151] in abinit tutorials
 atompaw: t02.abi …
 builtin: testin_v1.abi …
 libxc: t41.abi, t42.abi, t43.abi …
 paral: t71.abi, t71.abi, t71.abi …
 seq: tsv4_90.abi …
 tutorespfn: teph4mob_1.abi, teph4mob_4.abi, trf1_1.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t38.abi …
 v3: t09.abi, t43.abi, t60.abi …
 v4: t69.abi, t97.abi …
 v5: t54.abi, t60.abi, t66.abi …
 v6: t18.abi, t34.abi, t61.abi …
 v67mbpt: t02.abi, t04.abi, t19.abi …
 v7: t68.abi …
 v8: t21.abi, t22.abi …
 v9: t33.abi, t34.abi, t35.abi …
Used for screening length (in Bohr) of the model dielectric function, diagonal in reciprocal space. By default, given in Bohr atomic units (1 Bohr=0.5291772108 Angstrom), although Angstrom can be specified, if preferred, since dielng has the LENGTH characteristics. This model dielectric function is as follows ({\bf K} being a wavevector):
The inverse of this model dielectric function will be applied to the residual, to give the preconditioned change of potential. Right at {\bf K}=0, \epsilon({\bf K}) is imposed to be 1.
If the preconditioning were perfect, the change of potential would lead to an exceedingly fast solution of the selfconsistency problem (two or three steps). The present model dielectric function is excellent for rather homogeneous unit cells. When {\bf K}>0, it tends to the macroscopic dielectric constant, eventually divided by the mixing factor diemix (or diemixmag for magnetization). For metals, simply put diemac to a very large value (10^6 is OK) The screening length dielng governs the length scale to go from the macroscopic regime to the microscopic regime, where it is known that the dielectric function should tend to 1. It is on the order of 1 Bohr for metals with medium density of states at the Fermi level, like Molybdenum, and for Silicon. For metals with a larger DOS at the Fermi level (like Iron), the screening will be more effective, so that dielng has to be decreased by a factor of 24. This works for GS and RF calculations.
diemac¶
Mnemonics: model DIElectric MACroscopic constant
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: real
Dimensions: scalar
Default value: 1000000.0
Added in version: before_v9
Test list (click to open). Moderately used, [527/1101] in all abinit tests, [54/151] in abinit tutorials
 atompaw: t02.abi …
 bigdft: t00.abi, t20.abi, t21.abi …
 builtin: testin_fast.abi, testin_v5.abi, testin_bigdft.abi …
 etsf_io: t04.abi, t21.abi, t22.abi …
 fast: t00.abi, t01.abi, t02.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t05.abi, t09.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tmbt_1.abi, tpsic_01.abi, tpsic_02.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: tdepes_1.abi, tdepes_2.abi, tdepes_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t08.abi, t09.abi, t10.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t06.abi, t07.abi, t08.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t02.abi, t05.abi …
 v6: t03.abi, t04.abi, t12.abi …
 v67mbpt: t02.abi, t04.abi, t06.abi …
 v7: t01.abi, t03.abi, t09.abi …
 v8: t03.abi, t18.abi, t19.abi …
 v9: t05.abi, t07.abi, t09.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
A rough knowledge of the macroscopic dielectric constant diemac of the system is a useful help to speedup the SCF procedure: a model dielectric function, see the keyword dielng, is used for that purpose. It is especially useful for speeding up the treatment of rather homogeneous unit cells.
Some hint: The value of diemac should usually be bigger than 1.0, on physical grounds.
 For metals, simply put diemac to a very large value (the default 10^6 is OK)
 For silicon, use 12.0. A similar value is likely to work well for other semiconductors
 For wider gap insulators, use 2.0 … 4.0
 For molecules in an otherwise empty big box, try 1.5 … 3.0
Systems that combine a highly polarisable part and some vacuum are rather badly treated by the model dielectric function. One has to use the “extrapolar” technique, activated by the input variable iprcel. In sufficiently homogeneous systems, you might have to experiment a bit to find the best diemac . If you let diemac to its default value, you might even never obtain the selfconsistent convergence! For response function calculations, use the same values as for GS. The improvement in speed can be considerable for small (but nonzero) values of the wavevector.
diemix¶
Mnemonics: model DIElectric MIXing factor
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: real
Dimensions: scalar
Default value: 1.0 if %usepaw == 0 or iprcel !=0,
0.7 if %usepaw == 1 or iprcel == 0,
None otherwise.
Only relevant if: diemix >= 0.0 and diemix <= 1.0
Added in version: before_v9
Test list (click to open). Moderately used, [227/1101] in all abinit tests, [1/151] in abinit tutorials
 bigdft: t00.abi, t14.abi, t20.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v5.abi, testin_bigdft.abi …
 etsf_io: t21.abi, t22.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t03.abi, t05.abi, t05.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 paral: t05.abi, t05.abi, t05.abi …
 seq: tsv2_81.abi …
 tutorial: ttddft_1.abi …
 v1: t08.abi, t09.abi, t10.abi …
 v2: t01.abi, t03.abi, t04.abi …
 v3: t06.abi, t07.abi, t08.abi …
 v4: t06.abi, t07.abi, t35.abi …
 v5: t00.abi, t14.abi, t15.abi …
 v6: t07.abi, t12.abi, t13.abi …
 v7: t11.abi, t24.abi, t25.abi …
 v8: t21.abi, t22.abi, t30.abi …
 wannier90: t04.abi …
Gives overall factor of the preconditioned residual density/potential to be transferred in the SCF cycle. It should be between 0.0 and 1.0.
If the model dielectric function were perfect, diemix should be 1.0. By contrast, if the model dielectric function does nothing (when diemac = 1.0 or dielng is larger than the size of the cell), diemix can be used to damp the amplifying factor inherent to the SCF loop. For molecules, a value on the order 0.5 or 0.33 is rather usual. When mod(iscf,10)=3, 4,5 or 7, diemix is only important at the few first iterations when anharmonic effects are important, since these schemes compute their own mixing factor for selfconsistency. Also note that a different value of diemix can be used for the magnetization (see diemixmag).
diemixmag¶
Mnemonics: model DIElectric MIXing factor for the MAGgnetization
Mentioned in topic(s): topic_spinpolarisation
Variable type: real
Dimensions: scalar
Default value: diemix if 70 < iprcel and iprcel < 80,
diemix if iprcel == 0,
diemix if iscf<10,
diemix otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [14/1101] in all abinit tests, [0/151] in abinit tutorials
Gives overall factor of the preconditioned residual magnetization/magnetic field to be transferred in the SCF cycle (see diemix for further information). For the time being, apply only when the SCF mixing is done on the density (iscf > =10).
A negative value of diemixmag means that magnetization is only preconditioned by ABS( diemixmag ), without the use of any preconditioner.
When SCF cycle has some difficulties to converge, changing the value of diemixmag can have a positive effect. In particular diemixmag = 4 is a good choice (i.e. diemixmag = 4, no other preconditioner on magnetization).
dosdeltae¶
Mnemonics: DOS DELTA in Energy
Characteristics: ENERGY
Mentioned in topic(s): topic_printing, topic_ElecDOS
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
Defines the linear grid resolution (energy increment) to be used for the computation of the DensityOfStates, when prtdos is nonzero. If dosdeltae is set to zero (the default value), the actual increment is 0.001 Ha if prtdos = 1 or 4 (smearing technique), and the much smaller value 0.00005 Ha if prtdos = 2, 3 or 5 (tetrahedron technique). This different default value arises because the smearing technique gives a quite smooth DOS, while the DOS from the tetrahedron method is rapidly varying.
enunit¶
Mnemonics: ENergy UNITs
Mentioned in topic(s): topic_Output, topic_ElecBandStructure
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [230/1101] in all abinit tests, [12/151] in abinit tutorials
 builtin: testin_v1.abi, testin_v5.abi, testin_etsf_io.abi …
 etsf_io: t00.abi, t09.abi, t21.abi …
 fast: t02.abi, t03.abi, t05.abi …
 libxc: t03.abi, t19.abi, t41.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t03.abi, t03.abi, t03.abi …
 seq: tsv2_81.abi …
 tutoplugs: tw90_4.abi …
 tutorespfn: tdepes_1.abi, tdepes_2.abi, tdepes_3.abi …
 tutorial: tbase3_5.abi, tgw2_1.abi, tgw2_2.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t11.abi, t12.abi, t26.abi …
 v3: t01.abi, t09.abi, t18.abi …
 v4: t01.abi, t03.abi, t04.abi …
 v5: t00.abi, t01.abi, t03.abi …
 v6: t04.abi, t34.abi, t37.abi …
 v67mbpt: t01.abi, t03.abi, t04.abi …
 v7: t07.abi, t16.abi, t26.abi …
 v8: t02.abi, t17.abi …
 wannier90: t03.abi …
Governs the units to be used for output of eigenvalues (and if any, phonon frequencies)
 0 → print eigenvalues in hartree;
 1 → print eigenvalues in eV;
 2 → print eigenvalues in both hartree and eV.
If phonon frequencies are to be computed:
 0 → phonon frequencies in Hartree and cm1;
 1 → phonon frequencies in eV and THz;
 2 → phonon frequencies in hartree, eV, cm1, Thz and Kelvin.
expert_user¶
Mnemonics: EXPERTise of the USER
Mentioned in topic(s): topic_UnitCell, topic_crystal, topic_SmartSymm, topic_GeoOpt, topic_kpoints
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: 9.2.2
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v9: t19.abi
If set to 0, the checking provided by ABINIT is maximum (default values of chkprim, chksymbreak, chksymtnons, chkdilatmx). If nonzero (up to three), the abovementioned checking input variables are all disabled (set to zero) although it is still possible to activate particular tests by specifying input variables directly in the input file. In the future, the level three will always be the maximum allowed value, with all checks set to zero, while a more refined behaviour might be implemented for expert_user ==1 or 2).
extfpmd_nbcut¶
Mnemonics: EXTended FPMD: Number of Bands at CUT
Mentioned in topic(s): topic_ExtFPMD
Variable type: integer
Dimensions: scalar
Default value: 25
Added in version: 9.5.2
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v9: t92.abi
Specify the number of bands to use when averaging over last bands to get the energy shift factor when useextfpmd = 2 or 3.
extfpmd_nbcut must be less than nband.
fband¶
Mnemonics: Factor for the number of BANDs
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_BandOcc
Variable type: real
Dimensions: scalar
Default value: 0.125 if occopt == 1,
0.5 if occopt > 2,
0.0 if usewvl == 1,
0.0 otherwise.
Added in version: before_v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [0/151] in abinit tutorials
Governs the number of bands to be used in the code in the case the parameter nband is not defined in the input file (which means that occopt is not equal to 0 or 2).
In case fband is 0.0, the code computes from the pseudopotential files and the geometry data contained in the input file, the number of electrons present in the system. Then, it computes the minimum number of bands that can accommodate them, and use that value for nband.
In case fband differs from zero, other bands will be added, just larger than fband times the number of atoms. This parameter is not echoed in the top of the main output file, but only the parameter nband that it allowed to compute. It is also not present in the dtset array (no internal). The default values are chosen such as to give naturally some conduction bands. This improves the robustness of the code, since this allows one to identify lack of convergence coming from (near)degeneracies at the Fermi level. In the metallic case, the number of bands generated might be too small if the smearing factor is large. The occupation numbers of the higher bands should be small enough such as to neglect higher bands. It is difficult to automate this, so a fixed default value has been chosen.
fock_icutcoul¶
Mnemonics: Integer that governs the CUToff for COULomb interaction
Mentioned in topic(s): topic_Hybrids, topic_Coulomb
Variable type: integer
Dimensions: scalar
Default value: 3
Added in version: before_v9
Test list (click to open). Rarely used, [0/1101] in all abinit tests, [0/151] in abinit tutorials
!!!Under development!!! Electronic structure calculations for isolated systems, 1D and 2D systems present a slow convergence with respect to the size of the supercell due to the long ranged Coulomb interaction and the high degree of nonlocality of the operators involved. Thus, restricting the range of the Coulomb interaction, in order to prevent supercell images to interact can significantly speedup convergence, or even can make convergence happen. Also, even in the groundstate case, a cutoff Coulomb interaction might prove useful.
fock_icutcoul defines the particular expression to be used for the Fock operator in reciprocal space (see icutcoul for the Hartree contributions to ground state calculations, and gw_icutcoul for the corresponding treatment in GW calculations).
The choice of fock_icutcoul depends on the dimensionality and the character of the XC functional used (or otherwise the presence of the exclusive treatment of the shortrange exchange interaction). Possible values of fock_icutcoul are from 0 to 5, but currently are available options 0 and 5. Option 5 is hard coded as the method to be applied to HSE functionals.
Like for icutcoul, for 1dimensional and 2dimensional systems, the geometry of the system has to be specified explicitly. This is done thanks to vcutgeo. For 0, 1 and 2dimensional systems, a cutoff length has to be provided, thanks to rcut.
 0 → Sphere (molecules, but also 3Dcrystals, see below). See rcut.
 1 → (W.I.P.) cylinder (nanowires, nanotubes). See vcutgeo and rcut.
 2 → (W.I.P) Surface. See vcutgeo and rcut.
 3 → (W.I.P) 3D crystal (Coulomb interaction without cutoff).
 4 → (W.I.P.)ERF, longrange only Coulomb interaction.
 5 → ERFC, shortrange only Coulomb interaction (e.g. as used in the HSE functional).
Note that Spencer and Alavi showed that the spherical cutoff can efficiently be used also for 3D systems [Spencer2008]. In the latter case, use a negative value for the cutoff radius of the sphere (rcut<0), which is automatically calculated so that the volume enclosed in the sphere is equal to the volume of the solid.
fockdownsampling¶
Mnemonics: FOCK operator, kgrid DOWNSAMPLING
Mentioned in topic(s): topic_Hybrids
Variable type: integer
Dimensions: (3)
Default value: 3*1
Only relevant if: %usefock == 1
Added in version: before_v9
Test list (click to open). Rarely used, [10/1101] in all abinit tests, [0/151] in abinit tutorials
The k wavevector grid used to compute the Fock operator in the full Brillouin Zone might be a subset of the full Brillouin Zone of the k point grid used for the wavefunctions (see kptopt for its definition). fockdownsampling allows one to define such a reduced set of k wavevectors, yielding a significant speed up at the expense (possibly) of accuracy. The final grid has nkpthf k points, with coordinates %kptns_hf.
In the simplest case, the three values of fockdownsampling are equal, and simply gives the factor by which the initial k point grid will be downsampled along every direction. If fockdownsampling is 3*1, the sampling for the Fock operator is the same as the sampling for the wavefunctions. If fockdownsampling is 3*2 the sampling will have 8 times less k points. Conventionally, if fockdownsampling is 3*0, then the Fock operator is obtained solely from the Gamma point. Also, as soon as fockdownsampling is not 3*1 or 3*0, the k point grid from which a subset will be taken is obtained by imposing nshiftk = 1.
A more accurate description is now given, as one can achieve a better control than described above, with differing values of fockdownsampling and also with negative numbers. One starts from the k point grid defined by kptrlatt, with nshiftk = 1. The absolute value of each of the three numbers of fockdownsampling is used to sample the corresponding axis (in reduced coordinate), as described above. Moreover, the obtained k grid might even be further downsampled by specifying negative numbers: if all three are negative, one further downsamples each 2x2x2 subset of k points by taking its FCC subset (so, 4 points instead of 8); if two are negative, one downsamples each 2x2x2 subset of k points by taking its BCC subset (so, 2 points instead of 8); is only one is negative, then the two other axes are sampled using a facecentered sampling. Finally, some of the values might be taken as 0, in which case the corresponding direction is sampled by only one layer of points (if two are zero, a line of points is obtained).
fockoptmix¶
Mnemonics: FOCK operator: OPTions for MIXing
Mentioned in topic(s): topic_Hybrids
Variable type: integer
Dimensions: scalar
Default value: 0
Only relevant if: %usefock == 1
Added in version: before_v9
Test list (click to open). Rarely used, [5/1101] in all abinit tests, [0/151] in abinit tutorials
Governs the mixing algorithm at the level of the Fock operator, i.e. how to mix it, and how the underlying SCF calculation is to be performed. It is the most relevant when the Fock operator is not updated at each SCF step (nnsclohf/=0).
The last digit of fockoptmix governs what happens at the level of the SCF algorithm, when the Fock operator is updated.
 If fockoptmix == 0: the SCF algorithm is not restarted (it continues to use the previous potential/density pairs without worrying).
 If fockoptmix == 1: the SCF algorithm is restarted (the previous potential/density pairs are discarded).
The secondtolast (tens) digit governs the possible modification of the XC functional inside the SCF loop to take into account the lack of update of the Fock operator. Irrelevant when the unit digit is 0. If the value 1 is used (so, e.g. fockoptmix == 11), an auxiliary xc functional is used inside the SCF loop, for a frozen ACE Fock operator. This auxiliary functional is specified thanks to auxc_ixc and auxc_scal.
The thirdtolast (hundreds) digit governs the mixing of the Fock operator itself with its previous occurrences. Irrelevant when the unit digit is 0.
hyb_mixing¶
Mnemonics: HYBrid MIXING coefficient for unscreened fock operator
Mentioned in topic(s): topic_Hybrids
Variable type: real
Dimensions: scalar
Default value: 999.0
Comment: With the default, hyb_mixing is initialized from ixc.
Only relevant if: %usefock > 0
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [0/151] in abinit tutorials
Mixing coefficient for the unscreened Fock operator in case of hybrid functionals. HartreeFock corresponds to 1.0, PBE0 to 0.25, and B3LYP to 0.2.
ABINIT knows the correct value from ixc. Experts might nevertheless tune this mixing coefficient.
hyb_mixing_sr¶
Mnemonics: HYBrid MIXING coefficient for ShortRange screened fock operator
Mentioned in topic(s): topic_Hybrids
Variable type: real
Dimensions: scalar
Default value: 999.0
Comment: With the default, hyb_mixing_sr is initialized from ixc.
Only relevant if: %usefock > 0
Added in version: before_v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [0/151] in abinit tutorials
Mixing coefficient for the screened Fock operator in case of hybrid functionals. HSE has 0.25.
ABINIT knows the correct value from ixc. Experts might nevertheless tune this mixing coefficient.
hyb_range_dft¶
Mnemonics: HYBrid RANGE for the DFT leftover from the screened fock operator
Mentioned in topic(s): topic_Hybrids
Variable type: real
Dimensions: scalar
Default value: 999.0 or hyb_range_fock if it is defined by the user
Comment: With the default=999.0, hyb_range_dft is initialized from ixc.
Only relevant if: %usefock > 0
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
Range of the DFT leftover from the screened Fock operator in case of hybrid functionals (actually, coefficient of the distance appearing in the erf function, thus it has the dimension of an inverse distance).
As described in the LibXC sources (and copied in the ABINIT doc, see ixc = 428), there is a mess due to an error in the original publication. ABINIT knows the LibXC value from ixc, that might not agree with the definitions from other codes. Usually, hyb_range_dft is the same as hyb_range_fock, see the latter for the different values. However, there is a noticeable exception, the HSE03 from the original paper (not the HSE03 from VASP), for which hyb_range_dft = 0.188988 while hyb_range_fock = 0.106066.
hyb_range_fock¶
Mnemonics: HYBrid RANGE for the screened FOCK operator
Mentioned in topic(s): topic_Hybrids
Variable type: real
Dimensions: scalar
Default value: 999.0 or hyb_range_dft if it is defined by the user
Comment: With the default=999.0, hyb_range_fock is initialized from ixc.
Only relevant if: %usefock > 0
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
Range of the screened Fock operator in case of hybrid functionals (actually, coefficient of the distance appearing in the erf function, thus it has the dimension of an inverse distance).
As described in the LibXC sources (and copied in the ABINIT doc, see ixc = 428), there is a mess due to an error in the original publication. ABINIT knows the LibXC value from ixc, that might not agree with the definitions from other codes. Usually, hyb_range_dft is the same as hyb_range_fock , with one exception explained in hyb_range_dft. The HSE06 value from LibCX is 0.11, the one of Quantum Espresso is 0.106, the one of VASP is 0.105835 (=0.2 \AA^{1}). The HSE03 value from LibCX is 0.106066 (=0.15/\sqrt{2})), the one of VASP is 0.1587531 (=0.3 \AA^{1}).
iatsph¶
Mnemonics: Index for the ATomic SPHeres of the atomprojected densityofstates
Mentioned in topic(s): topic_printing, topic_ElecBandStructure, topic_ElecDOS
Variable type: integer
Dimensions: (natsph)
Default value: [ 1 … natsph ]
Only relevant if: prtdos == 3 or 4 or pawfatbnd in [1,2]
Added in version: before_v9
Test list (click to open). Rarely used, [9/1101] in all abinit tests, [2/151] in abinit tutorials
 mpiio: t26.abi
 tutorial: tpaw1_4.abi, tspin_3.abi
 v4: t35.abi, t38.abi
 v5: t19.abi, t20.abi, t40.abi
 v7: t04.abi
iatsph gives the number of the natsph atoms around which the sphere for atomprojected densityofstates will be build, in the prtdos = 3 or 4 cases. The radius of these spheres is given by ratsph. If pawfatbnd = 1 or 2, it gives the number of the natsph atoms around which atomprojected band structure will be built.
icoulomb¶
Mnemonics: Index for the COULOMB treatment
Mentioned in topic(s): topic_Coulomb
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [33/1101] in all abinit tests, [9/151] in abinit tutorials
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_bigdft.abi …
 tutoparal: tgswvl_1.abi, tgswvl_1.abi, tgswvl_1.abi …
Defines the type of computation (reciprocal space or real space) used for Hartree potential, local part of pseudopotential and ionion interaction:
 icoulomb = 0: usual reciprocal space computation, using icutcoul, gw_icutcoul and fock_icutcoul to define the Hartree potential, and using Ewald correction.
 icoulomb = 1: free boundary conditions are used when the Hartree potential is computed, real space expressions of pseudopotentials are involved (restricted to GTH pseudopotentials) and simple coulomb interaction gives the ionion energy. The wavelet Coulomb solver is used in this case.
icutcoul¶
Mnemonics: Integer that governs the CUToff for COULomb interaction
Mentioned in topic(s): topic_Coulomb
Variable type: integer
Dimensions: scalar
Default value: 3
Added in version: before_v9
Test list (click to open). Moderately used, [68/1101] in all abinit tests, [4/151] in abinit tutorials
 libxc: t44.abi, t45.abi …
 paral: t71.abi, t71.abi, t71.abi …
 tutoparal: tucrpa_2.abi …
 tutorial: tbs_2.abi, tbs_3.abi, tbs_4.abi …
 v3: t30.abi, t31.abi …
 v4: t84.abi, t85.abi, t86.abi …
 v5: t63.abi, t64.abi, t65.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t23.abi, t24.abi, t25.abi …
 v9: t31.abi, t90.abi …
 wannier90: t03.abi …
Electronic structure calculations for isolated systems, 1D and 2D systems present a slow convergence with respect to the size of the supercell due to the long ranged Coulomb interaction and the high degree of nonlocality of the operators involved. Thus, restricting the range of the Coulomb interaction, in order to prevent supercell images to interact can significantly speedup convergence, or even can make convergence happen. Also, even in the groundstate case, a cutoff Coulomb interaction might prove useful.
icutcoul defines the particular expression to be used for the Coulomblike terms in reciprocal space in groundstate calculations. See gw_icutcoul for GW calculationts, and fock_icutcoul for the Focklike terms in groundstate calculations e.g. using hybrid functionals. . The choice of icutcoul depends on the dimensionality of the system. Possible values of icutcoul are from 0 to 5. For 1dimensional and 2dimensional systems, the geometry of the system has to be specified explicitly. This is done thanks to vcutgeo. For 0, 1 and 2dimensional systems, a cutoff length has to be provided, thanks to rcut.
 0 → Sphere (molecules, but also 3Dcrystals, see below). See rcut.
 1 → (W.I.P.) cylinder (nanowires, nanotubes). See vcutgeo and rcut.
 2 → Surface. See vcutgeo and rcut.
 3 → 3D crystal (Coulomb interaction without cutoff).
 4 → ERF, longrange only Coulomb interaction.
 5 → ERFC, shortrange only Coulomb interaction (e.g. as used in the HSE functional). (W.I.P.)
Note that Spencer and Alavi showed that the spherical cutoff can efficiently be used also for 3D systems [Spencer2008]. In the latter case, use a negative value for the cutoff radius of the sphere (rcut<0), which is automatically calculated so that the volume enclosed in the sphere is equal to the volume of the solid.
iprcel¶
Mnemonics: Integer for PReConditioning of ELectron response
Mentioned in topic(s): topic_SCFAlgorithms
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [27/1101] in all abinit tests, [3/151] in abinit tutorials
Used when iscf > 0, to define the SCF preconditioning scheme. Potential based preconditioning schemes for the SCF loop (electronic part) are still a subject of active research. The present parameter (electronic part) describes the way the change of potential is derived from the residual. The possible values of iprcel correspond to:
 0 → model dielectric function described by diemac, dielng and diemix.
 larger or equal to 21 → will compute the dielectric matrix according to diecut, dielam, diegap. This methodology is described in [Anglade2008].
 Between 21 and 29 → for the first few steps uses the same as option 0 then compute RPA dielectric function, and use it as such.
 Between 31 and 39 → for the first few steps uses the same as option 0 then compute RPA dielectric function, and use it, with the mixing factor diemix.
 Between 41 and 49 → compute the RPA dielectric matrix at the first step, and recompute it at a later step, and take into account the mixing factor diemix.
 Between 51 and 59 → same as between 41 and 49, but compute the RPA dielectric matrix by another mean
 Between 61 and 69 → same as between 41 and 49, but compute the electronic dielectric matrix instead of the RPA one.
 Between 71 and 78 → STILL UNDER DEVELOPMENT – NOT USABLE; Use the modified Kerker preconditioner with a realspace formulation (basic formulation is shown at dielng). The dielectric matrix is approximated thanks to diemac and dielng. Note that diemix is also used.
 79 → STILL UNDER DEVELOPMENT – NOT USABLE; same as previous but with an alternate algorithm.
 141 to 169 → same as Between 41 and 69 (but, the dielectric matrix is also recomputed every iprcel modulo 10 step).
The computation of the dielectric matrix (for 0 [100]< iprcel < 70 [100]) is based on the extrapolar approximation, see [Anglade2008]. This approximation can be tuned with diecut, dielam, and diegap. Yet its accuracy mainly depends on the number of conduction bands included in the system. Having 2 to 10 empty bands in the calculation is usually enough (use nband).
NOTES:
 The step at which the dielectric matrix is computed or recomputed is determined by modulo( iprcel ,10). The recomputation happens just once in the calculation for iprcel < 100.
 For nonhomogeneous relatively large cells iprcel = 45 will likely give a large improvement over iprcel = 0.
 In case of PAW and iprcel > 0, see pawsushat input variable. By default, an approximation (which can be suppressed) is done for the computation of susceptibility matrix.
 For extremely large inhomogeneous cells where computation of the full dielectric matrix takes too many weeks, 70 < iprcel < 80 is advised.
 For nsppol = 2 or nspinor = 2 with metallic occopt, only mod(iprcel,100) <50 is allowed.
 No meaning for RF calculations yet.
 The exchange term in the full dielectric matrix diverges for vanishing densities. Therefore the values of iprcel beyond 60 must not be used for cells containing vacuum, unless ones computes this matrix for every step ( iprcel = 161).
iqpt¶
Mnemonics: Index for QPoinT generation
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [11/1101] in all abinit tests, [3/151] in abinit tutorials
 paral: t60.abi, t60.abi, t60.abi
 tutorespfn: tdepes_3.abi, tdepes_4.abi, teph_1.abi
 v6: t35.abi, t78.abi
 v7: t50.abi, t51.abi, t55.abi
Only used if nqpt = 1, and qptopt = 1 to 4.
Defines the index of the Q point to be selected in the list of q points generated by ngqpt, qptrlatt, nshiftq, and shiftq.
If iqpt = 0, then the q point is Gamma (0 0 0).
The usual working mode is to define a series of values for iqpt , starting with iqpt = 0 or 1 (so through the definition of iqpt: ), and increasing it by one for each dataset (thanks to iqpt+ ).
ivalence¶
Mnemonics: Index of the highest VALENCE band
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_BandOcc, topic_DeltaSCF
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [0/151] in abinit tutorials
Only used if [[occopt] = 9.
When occopt==9, the lowest ivalence bands are considered to be valence bands and nqfd holes are constrained to exist in them. Accordingly, higher energy bands (with index >= ivalence + 1) are considered to be conduction bands, and nqfd electrons are constrained there. See [Paillard2019] for more details about the general method implemented.
ixcpositron¶
Mnemonics: Integer for the eXchangeCorrelation applied to the electronPOSITRON interaction
Mentioned in topic(s): topic_positron
Variable type: integer
Dimensions: scalar
Default value: 1
Comment: (Teter parameterization). However, if all the pseudopotentials have the same value of pspxc, the initial value of ixc will be that common value
Only relevant if: positron/=0
Added in version: before_v9
Test list (click to open). Rarely used, [10/1101] in all abinit tests, [7/151] in abinit tutorials
 tutorial: tpositron_1.abi, tpositron_2.abi, tpositron_3.abi, tpositron_4.abi, tpositron_5.abi, tpositron_6.abi, tpositron_7.abi
 v4: t96.abi
 v5: t07.abi
 v7: t35.abi
Relevant only when positron/=0. Define the type of electronpositron correlation that is used in case of a electronpositron twocomponent DFT calculation. Define also the analytical formula of the enhancement factor used to compute the electronpositron annihilation rate:
Electronpositron correlation functional:
 ixcpositron=1: LDA zero positron density limit parametrized by Arponen & Pajanne and provided by Boronski & Nieminen [Arponen1979a],[Boronski1986]
 ixcpositron=11: LDA zero positron density limit parametrized by Arponen & Pajanne and fitted by Sterne & Kaiser [Arponen1979a],[Sterne1991]
 ixcpositron=2: LDA electronpositron correlation provided by Puska, Seitsonen, and Nieminen [Arponen1979a],[Puska1995]
 ixcpositron=3: GGA zero positron density limit parametrized by Arponen & Pajanne and provided by Boronski & Nieminen [Arponen1979a],[Boronski1986],[Barbiellini1995]
 ixcpositron=31: GGA zero positron density limit parametrized by Arponen & Pajanne and fitted by Sterne & Kaiser [Arponen1979a],[Sterne1991],[Barbiellini1995]
Annihilation rate enhancement factor:
 ixcpositron=1: Boronski and Nieminen full modelisation and RPA limit [Arponen1979a]
 ixcpositron=11: Sterne and Kaiser [Boronski1986]
 ixcpositron=2: Puska, Seitsonen and Nieminen [Sterne1991]
 ixcpositron=3: Boronski and Nieminen full modelisation and RPA limit [Arponen1979a], with GGA corrections
 ixcpositron=31: Sterne and Kaiser [Boronski1986], with GGA corrections
jellslab¶
Mnemonics: include a JELLium SLAB in the cell
Mentioned in topic(s): topic_Artificial
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [0/151] in abinit tutorials
If set to 1, a slab of uniform positive background charge density, that is, a jellium slab, is included in the calculation cell. A portion of the unit cell is filled with such positive charge density distribution which is equal to a bulkmean value n_{bulk} between two edges and zero in the vacuum region if present. For the sake of convenience the unit cell is supposed to have the third crystal primitive lattice vector orthogonal to the other ones so that the portion of the cell filled by the jellium slab can be defined through its edges along z. The bulkmean positive charge density is fixed by the input variable slabwsrad, while the position of the slab edges along z is defined through the input variables slabzbeg and slabzend.
kptbounds¶
Mnemonics: K PoinT BOUNDarieS
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_kpoints, topic_ElecBandStructure
Variable type: real
Dimensions: (3,abs(kptopt)+1))
Default value: None
Added in version: before_v9
Test list (click to open). Moderately used, [15/1101] in all abinit tests, [5/151] in abinit tutorials
 libxc: t42.abi
 tutoparal: tdmft_1.abi, tucrpa_1.abi
 tutorespfn: tdepes_4.abi
 tutorial: tbase3_5.abi, tfold2bloch_1.abi
 v2: t43.abi
 v4: t89.abi
 v5: t20.abi, t54.abi
 v6: t14.abi, t34.abi
 v7: t76.abi, t77.abi
 v8: t04.abi
It is used to generate the circuit to be followed by the band structure, when kptopt is negative (it is not read if kptopt is zero or positive).
There are abs(kptopt) segments to be defined, each of which starting from the end point of the preceeding one. Thus, the number of points to be input is abs(kptopt)+1. They form a circuit starting at kptbounds (1:3,1)/kptnrm and ending at kptbounds (1:3,abs(kptopt)+1)/kptnrm. The number of divisions of each segment can be defined either using the array ndivk or the variable ndivsm (preferred) that just defines the number of divisions for the smallest segment.
As for kpt, kptbounds is specified using the primitive vectors in reciprocal space. If your Bravais lattice is simple, then it should be quite easy to find the coordinates of the end points. On the other hand, for centered, bodycentered, facecentered, hexagonal, and rhombohedral Bravais lattice, the conversion might be more difficult. See the description of kpt for an explanation of how to convert data from the “conventional” cartesian coordinates to the primitive vectors in the reciprocal space. In order to help a bit, we list below a series of typical values, for the FCC, BCC, hexagonal and rhombohedral Bravais lattices. Note: all the data below are given in dimensionless units; they have to be rescaled by the actual lengths defined by the acell values. However, kptbounds values can be used as such, if the values of rprim given below are adopted.
A. FCC lattice
Suppose the primitive vectors in real space are given by
rprim 0 1 1 1 0 1 1 1 0
or
rprim 0 1/2 1/2 1/2 0 1/2 1/2 1/2 0
(these two possibilities only differ by a scaling factor, irrelevant for the definition of the k points in the primitive vectors in reciprocal space). Then, the reciprocal primitive vectors (in conventional cartesian coordinates) are
(1/2 1/2 1/2), (1/2 1/2 1/2), (1/2 1/2 1/2)
or
(1 1 1), (1 1 1), (1 1 1)
and, in both cases, the coordinates of several special points with respect to primitive vectors in reciprocal space are
X (0 1/2 1/2) (conventional cartesian coordinate 1/2 0 0)
X'(1/2 1/2 1 ) (conventional cartesian coordinate 1/2 1/2 0) (an other instance of X, in another Brillouin zone)
L (1/2 1/2 1/2) (conventional cartesian coordinate 1/4 1/4 1/4)
L'(1/2 0 0 ) (conventional cartesian coordinate 1/4 1/4 1/4) (an other instance of L, on another face of the BZ)
W (1/4 1/2 3/4) (conventional cartesian coordinate 1/2 1/4 0)
U (1/4 5/8 5/8) (conventional cartesian coordinate 1/2 1/8 1/8)
K (3/8 3/8 3/4) (conventional cartesian coordinate 3/8 3/8 0)
Note that K is actually equivalent to U, by spatial and translational symmetry. So, if you want to specify a typical circuit, the following might do the work: LGammaXWK,ULWXK,UGamma with
kptbounds 1/2 0 0 0 0 0 0 1/2 1/2 1/4 1/2 3/4 3/8 3/8 3/4 1/2 1/2 1/2 1/4 1/2 3/4 1/2 1/2 1 3/8 3/8 3/4 0 0 0
The lengths of segments (this information is useful to draw the band structure, with the correct relative scale between special points) can be found using the conventional cartesian coordinates:
l(LGamma)=\sqrt{3}/4=0.433… \n l(GammaX)=1/2=0.5 \n l(XW)=1/4=0.25 \n l(WK)=\sqrt{2}/8=0.177… \n l(KL)=\sqrt{6}/8=0.306… \n l(LW)=\sqrt{2}/4=0.354… \n l(WX)=1/4=0.25 \n l(XK)=\sqrt{2}/8=0.177… \n l(KGamma)=3\sqrt{2}/8=0.530… \n
B. BCC lattice
Suppose the primitive vectors in real space are given by
rprim 1 1 1 1 1 1 1 1 1
(as for the FCC lattice, there is a scale invariance). Then, the reciprocal primitive vectors (in conventional cartesian coordinates) are (0 ½ ½), (½ 0 ½), and (½ ½ 0) and the coordinates of several special points with respect to primitive vectors in reciprocal space are
H (1/2 1/2 1/2) (conventional cartesian coordinate 1/2 0 0)
N ( 0 0 1/2) (conventional cartesian coordinate 1/4 1/4 0)
P ( 1/4 1/4 1/4) (conventional cartesian coordinate 1/4 1/4 1/4)
So, if you want to specify a typical circuit, the following might do the work: GammaHNGammaPNPH
kptbounds 0 0 0 1/2 1/2 1/2 0 0 1/2 0 0 0 1/4 1/4 1/4 0 0 1/2 1/4 1/4 1/4 1/2 1/2 1/2
The lengths of segments (this information is useful to draw the band structure, with the correct relative scale between special points) can be found using the conventional cartesian coordinates:
l(GammaH)=1/2=0.5 \n l(HN)=\sqrt{2}/4=0.354… \n l(NGamma)=\sqrt{2}/4=0.354… \n l(GammaP)=\sqrt{3}/4=0.433… \n l(PN)=1/4=0.25 \n l(NP)=1/4=0.25 \n l(PH)=\sqrt{3}/4=0.433… \n
C. Hexagonal lattices
Suppose the primitive vectors in real space are given by
rprim 1 0 0 1/2 sqrt(0.75) 0 0 0 1
The coordinates of several special points with respect to primitive vectors in reciprocal space are
M (1/2 0 0) or (0 1/2 0) or (1/2 1/2 0)
L (1/2 0 1/2) or (0 1/2 1/2) or (1/2 1/2 1/2)
K (1/3 1/3 0) or (2/3 1/3 0) or (1/3 2/3 0)
H (1/3 1/3 1/2) or (2/3 1/3 1/2) or (1/3 2/3 1/2)
A (0 0 1/2)
So, if you want to specify a typical circuit, the following might do the work: KGammaMKHALHLMGammaA
kptbounds 1/3 1/3 0 0 0 0 1/2 0 0 1/3 1/3 0 1/3 1/3 1/2 0 0 1/2 1/2 0 1/2 1/3 1/3 1/2 1/2 0 1/2 1/2 0 0 0 0 0 0 0 1/2
In order to find the lengths of segments (this information is useful to draw the band structure, with the correct relative scale between special points) one needs to know the a and c lattice parameters. Also, in what follows, we omit the 2\pi factor sometimes present in the definition of the reciprocal space vectors. The reciprocal vectors are (1/a\: 1/(\sqrt{3}a)\: 0), (0\: 2/(\sqrt{3}a)\: 0), (0\: 0\: 1/c). The lengths of the abovementioned segments can be computed as:
l(KGamma)=2/(3a)=0.666…/a \n l(GammaM)=1/(\sqrt{3}a)=0.577…/a \n l(MK)=1/(3a)=0.333…/a \n l(KH)=1/(2c)=0.5…/c \n l(HA)=2/(3a)=0.666…/a \n l(AL)=1/(\sqrt{3}a)=0.577…/a \n l(LH)=1/(3a)=0.333…/a \n l(HL)=1/(3a)=0.333…/a \n l(LM)=1/(2c)=0.5…/c \n l(MGamma)=1/(\sqrt{3}a)=0.577…/a \n l(GammaA)=1/(2c)=0.5…/c \n
D. Rhombohedral lattices
Rhombohedral lattices are characterised by two parameters, the length of the primitive vectors, that we will denote a0, and the angle they form, \gamma. These can be directly input of ABINIT, as acell and angdeg
This will generate the primitive vectors in real space, with
acell a0 a0 a0 and rprim a 0 c a/2 a*sqrt(0.75) c a/2 a*sqrt(0.75) c
with,
 a^2+c^2=1,
 a^2=2/3(1\cos(\gamma)),
 c^2=1/3(1+2\cos(\gamma)),
 (a/c)^2=2(1\cos(\gamma))/(1+2\cos(\gamma)), and also
 \cos(\gamma)=(1(a/c)^2/2)/(1+(a/c)^2).
Alternatively, these values of rprim might directly be the input of ABINIT (then, the balance of the scaling factor might be adjusted between acell and rprim).
Unlike for the simple cubic, FCC, BCC, hexagonal (and some other) Bravais lattice, the topology of the Brillouin zone will depend on the \gamma (or a/c) value. We give below information concerning the case when \cos(\gamma) is positive, that is, (a/c)^2 lower than 2.
The coordinates of several special points with respect to primitive vectors in reciprocal space will not depend on the a/c ratio, but some others will depend on it. So, some care has to be exercised. Notations for the Brillouin Zone special points are the same as in [Gonze1990].
L (1/2 0 0) or (0 1/2 0) or (0 0 1/2) (or with negative signs)
T (1/2 1/2 1/2)
X (1/2 1/2 0) or (1/2 0 1/2) or (0 1/2 1/2) (or with separate negative signs)
W (5/6  (a/c)^2/6, 1/2, 1/6 + (a/c)^2/6 ) = (1 0 1)*(1(a/c)^2/2)/3 + (1 1 1)/2
U ( (1+(a/c)^2)/6, (8(a/c)^2)/12, (8(a/c)^2)/12 ) = (1 1/2 1/2)*(1(a/c)^2/2)/3 + (1 1 1)/2
K (1 0 1)*(1+(a/c)^2/4)/3
So, if you want to specify a typical circuit, the following might do the work (the representative points on lines of symmetry are indicated  there are sometimes more than one way to go from one point to another): XVKSigma GammaLambdaTQWYLsigmaGammasigmaX. The suggestion is to sample this path with the following coordinates for the special points X, Gamma, T, L, Gamma, X:
kptbounds 1/2 0 1/2 0 0 0 1/2 1/2 1/2 1 1/2 0 1 0 0 1 1/2 1/2
In order to find the lengths of segments (this information is useful to draw the band structure, with the correct relative scale between special points) one needs to know the a and c lattice parameters. Also, in what follows, we omit the 2\pi factor sometimes present in the definition of the reciprocal space vectors. The reciprocal vectors are ( 2/(3a)\: 0\: 1/(3c) ), ( 1/(3a)\: 1/(\sqrt{3}a)\: 1/(3c) ), ( 1/(3a)\: 1/(\sqrt{3}a)\: 1/(3c) ). The lengths of the abovementioned segments can be computed as:
l(XGamma)=2/(\sqrt{3}a)=1.155…/a \n l(KGamma)=4(1+(a/c)^2/4)/(3\sqrt{3}a) \n l(GammaT)=1/(2c) \n l(TL)=2/(\sqrt{3}a)=1.155…/a \n l(TW)=4(1(a/c)^2/2)/(3\sqrt{3}a) \n l(LGamma)=\sqrt{4/(a^2)+1/(c^2)}/3 \n l(GammaX)=2\sqrt{1/(a^2)+1/(c^2)}/3 \n
kptrlatt¶
Mnemonics: K  PoinTs grid: Real space LATTice
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: (3,3)
Default value: 0
*The use of this variable forbids the use of: ngkpt
Added in version: before_v9
Test list (click to open). Moderately used, [64/1101] in all abinit tests, [7/151] in abinit tutorials
 gpu: t03.abi, t05.abi, t05.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t21.abi, t22.abi, t24.abi …
 tutorespfn: trf1_1.abi, trf1_2.abi, trf1_3.abi …
 v2: t12.abi, t74.abi, t99.abi …
 v3: t22.abi …
 v4: t69.abi …
 v5: t08.abi, t16.abi, t18.abi …
 v6: t07.abi, t28.abi, t36.abi …
 v67mbpt: t22.abi …
 v7: t03.abi, t14.abi, t21.abi …
 v8: t01.abi, t01_triqs2_0.abi, t03.abi …
 v9: t66.abi …
This input variable is used only when kptopt is positive. It partially defines the k point grid. The other piece of information is contained in shiftk. kptrlatt cannot be used together with ngkpt.
The values kptrlatt (1:3,1), kptrlatt (1:3,2), kptrlatt (1:3,3) are the coordinates of three vectors in real space, expressed in the %rprimd coordinate system (reduced coordinates). They defines a superlattice in real space. The k point lattice is the reciprocal of this superlattice, possibly shifted (see shiftk).
If neither ngkpt nor kptrlatt are defined, ABINIT will automatically generate a set of k point grids, and select the best combination of kptrlatt and shiftk that allows one to reach a sufficient value of kptrlen. See this latter variable for a complete description of this procedure.
kptrlen¶
Mnemonics: K  PoinTs grid: Real space LENgth
Mentioned in topic(s): topic_kpoints
Variable type: real
Dimensions: scalar
Default value: 30.0
Added in version: before_v9
Test list (click to open). Moderately used, [18/1101] in all abinit tests, [1/151] in abinit tutorials
This input variable is used only when kptopt is positive and nonzero.
Preliminary explanation: The k point lattice defined by ngkpt or kptrlatt is used to perform integrations of periodic quantities in the Brillouin Zone, like the density or the kinetic energy. One can relate the error made by replacing the continuous integral by a sum over k point lattice to the Fourier transform of the periodic quantity. Erroneous contributions will appear only for the vectors in real space that belong to the reciprocal of the k point lattice, except the origin. Moreover, the expected size of these contributions usually decreases exponentially with the distance. So, the length of the smallest of these real space vectors is a measure of the accuracy of the k point grid.
When either ngkpt or kptrlatt is defined, kptrlen is not used as an input variable, but the length of the smallest vector will be placed in this variable, and echoed in the output file.
On the other hand, when neither ngkpt nor kptrlatt are defined, ABINIT will automatically generate a large set of possible k point grids, and select among this set, the grids that give a length of smallest vector LARGER than kptrlen , and among these grids, the one that, when used with kptopt = 1, reduces to the smallest number of k points. Note that this procedure can be timeconsuming. It is worth doing it once for a given unit cell and set of symmetries, but not use this procedure by default. The best is then to set prtkpt = 1, in order to get a detailed analysis of the set of grids.
If some layer of vacuum is detected in the unit cell (see the input variable vacuum), the computation of kptrlen will ignore the dimension related to the direction perpendicular to the vacuum layer, and generate a bi dimensional k point grid. If the system is confined in a tube, a one dimensional k point grid will be generated. For a cluster, this procedure will only generate the Gamma point.
magcon_lambda¶
Mnemonics: MAGnetization CONstraint LAMBDA parameter
Mentioned in topic(s): topic_MagMom
Variable type: real
Dimensions: scalar
Default value: 0.01
Added in version: before_v9
Test list (click to open). Moderately used, [13/1101] in all abinit tests, [0/151] in abinit tutorials
This variable gives the amplitude of the penalty function imposed on the magnetization vectors on each atom (turned on with flag variable magconon=1 to 3). Typical values for magcon_lambda are 0.001 to 0.1. The SCF convergence will be difficult if magcon_lambda is too large. If magcon_lambda is too small, the penalty will not be very effective and it will give magnetization not close to the desired spinat target. In case of convergence problem, it can help to start with a small value of magcon_lambda and to increase it by reading the wavefunction obtained with a lower magcon_lambda value. See variable magconon for more details.
magconon¶
Mnemonics: turn MAGnetization CONstraint ON
Mentioned in topic(s): topic_MagMom
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v7: t05.abi
Turns on the imposition of a constraint on the magnetization, using a penalty function. For each atom, the magnetization is calculated in a sphere (radius ratsph) and a penalty function is applied to bring it to the input values of spinat. The constraint can be either on the direction/axis only ( magconon = 1) or on the full vector ( magconon = 2). The penalty function has an amplitude magcon_lambda that should be neither too big (bad or impossible convergence) nor too small (no effect). The penalty function is documented in [Ma2015] as being a Lagrange approach, which is a misnomer for the algorithm that they describe. It has the drawback of being unable to deliver the exact sought value for the magnetization. So, the true Lagrange approach has to be preferred, except for testing purposes. This is provided by the algorithm governed by the input variable constraint_kind, which is actually also much more flexible than the implementation corresponding to magconon .
Final subtlety: when magconon = 1, if nspden==2 (collinear case), then the direction of magnetization is constraint (positive or negative along z), while if nspden==4, then the axis of magnetization is constraint (the actual direction is not imposed, both directions are equivalent). This might be confusing.
mixalch¶
Mnemonics: MIXing coefficients for ALCHemical potentials
Characteristics: EVOLVING
Mentioned in topic(s): topic_AtomTypes
Variable type: real
Dimensions: (%npspalch,ntypalch)
Default value: None
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [0/151] in abinit tutorials
Used for the generation of alchemical pseudoatoms, that is, when ntypalch is nonzero.
This array gives, for each type of alchemical pseudoatom (there are ntypalch such pseudoatoms), the mixing coefficients of the basic %npspalch pseudopotentials for alchemical use. For each type of alchemical pseudoatom, the sum of the mixing coefficients must equal 1.
The actual use of the mixing coefficients is defined by the input variable algalch. Note that the masses of the atoms, amu are also mixed according to the value of mixalch , by default.
Example 1. Suppose that we want to describe Ba(0.25) Sr(0.75) Ti O_3. The input variables related to the construction of the alchemical Ba(0.25) Sr(0.75) potential will be:
npsp 4 ! 4 pseudopotentials should be read.
znucl 8 40 56 38 ! The nuclear charges. Note that the two
! atoms whose pseudopotentials are to be mixed
! are mentioned at the end of the series.
ntypat 3 ! There will be three types of atoms.
ntypalch 1 ! One pseudoatom will be alchemical.
! Hence, there will be ntyppure=2 pure pseudoatoms,
! with znucl 8 (O) and 40 (Ti), corresponding to
! the two first pseudopotentials. Out of the
! four pseudopotentials, npspalch=2 are left
! for alchemical purposes, with znucl 56 (Ba)
! and 38 (Sr).
mixalch 0.25 0.75 ! For that unique pseudoatom to be
! generated, here are the mixing coefficients,
! to be used to combine the Ba and Sr pseudopotentials.
Example 2. More complicated, and illustrate some minor drawback of the design of input variables. Suppose that one wants to generate Al(0.25)Ga(0.75) As(0.10)Sb(0.90). The input variables will be:
npsp 4 ! 4 pseudopotentials should be read
znucl 13 31 33 51 ! The atomic numbers. All pseudopotentials
! will be used for some alchemical purpose
ntypat 2 ! There will be two types of atoms.
ntypalch 2 ! None of the atoms will be "pure".
! Hence, there will be npspalch=4 pseudopotentials
! to be used for alchemical purposes.
mixalch 0.25 0.75 0.0 0.0 ! This array is a (4,2) array, arranged in the
0.0 0.0 0.1 0.9 ! usual Fortran order.
Minor drawback: one should not forget to fill mixalch with the needed zero’s, in this later case.
In most cases, the use of mixalch will be as a static (nonevolving) variable. However, the possibility to have different values of mixalch for different images has been coded. A population of cells with different atomic characteristics can thus be considered, and can be made to evolve, e.g. with a genetic algorithm (not coded in v7.0.0 though). There is one restriction to this possibility: the value of %ziontypat for the atoms that are mixed should be identical.
natsph¶
Mnemonics: Number of ATomic SPHeres for the atomprojected densityofstates
Mentioned in topic(s): topic_printing, topic_ElecBandStructure, topic_ElecDOS
Variable type: integer
Dimensions: scalar
Default value: natom
Only relevant if: prtdos == 3 or 4 or pawfatbnd in [1,2]
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [2/151] in abinit tutorials
 mpiio: t26.abi
 tutorial: tpaw1_4.abi, tspin_3.abi
 v4: t38.abi
 v5: t19.abi, t20.abi, t40.abi
 v7: t04.abi
natsph gives the number of atoms around which the sphere for atomprojected densityofstates will be built, in the prtdos = 3 or 4 case. The indices of these atoms are given by iatsph. The radius of these spheres is given by ratsph. If pawfatbnd = 1 or 2, it gives the number of atoms around which atomprojected band structure will be built (the indices of these atoms are given by iatsph).
natsph_extra¶
Mnemonics: Number of ATomic SPHeres for the lprojected densityofstates in EXTRA set
Mentioned in topic(s): topic_printing
Variable type: integer
Dimensions: scalar
Default value: 0
Only relevant if: prtdos == 3 or 4 or pawfatbnd in [1,2]
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v7: t04.abi
natsph_extra gives the number of extra spheres for which the angular momentumprojected densityofstates will be built, in the prtdos = 3 or 4 case. The radius of these spheres is given by ratsph_extra. This simulates the STS signal for an STM tip atom placed at the sphere position, according to the chemical nature of the tip (s p d wave etc…). If pawfatbnd = 1 or 2, it gives the number of spheres in which lprojected band structure will be built. The position of the spheres is given by the xredsph_extra variable.
nbdbuf¶
Mnemonics: Number of BanDs for the BUFfer
Mentioned in topic(s): topic_SCFControl, topic_BandOcc
Variable type: integer
Dimensions: scalar
Default value: 2*nspinor if optdriver == 0 and iscf<0,
2*nspinor if optdriver == 1 and 3<=occopt and occopt<= 8,
0 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [205/1101] in all abinit tests, [20/151] in abinit tutorials
 bigdft: t22.abi …
 fast: t09.abi, t11.abi, t12.abi …
 libxc: t10.abi, t67.abi, t68.abi …
 mpiio: t69.abi, t69.abi …
 paral: t05.abi, t05.abi, t05.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tmbt_1.abi, tucrpa_2.abi …
 tutorespfn: tdepes_1.abi, tdepes_2.abi, tdepes_3.abi …
 tutorial: tbs_1.abi, tgw1_1.abi, tgw1_2.abi …
 v2: t07.abi, t11.abi, t26.abi …
 v3: t18.abi, t30.abi, t31.abi …
 v4: t02.abi, t52.abi, t56.abi …
 v5: t23.abi, t39.abi, t41.abi …
 v6: t08.abi, t09.abi, t10.abi …
 v67mbpt: t01.abi, t02.abi, t07.abi …
 v7: t41.abi, t47.abi, t55.abi …
 v8: t36.abi, t41.abi, t47.abi …
 v9: t05.abi, t07.abi, t09.abi …
nbdbuf gives the number of bands, the highest in energy, that, among the nband bands, are to be considered as part of a buffer. A negative value is interpreted as percentage of nband (added in v9).
This concept is useful in three situations: in nonselfconsistent calculations, for the determination of the convergence tolerance; for response functions of metals, to avoid instabilities, and also when finite electric fields or nonlinear responses (with electric field perturbations) are considered. For the two first, the need of a buffer is a natural requirement of the problem, so that the default value is changed to 2 automatically, as explained in the following. The third case is only for implementation convenience.
In nonselfconsistent GS calculations (iscf<0), the highest levels might be difficult to converge, if they are degenerate with another level, that does not belong to the set of bands treated. Then, it might take extremely long to reach tolwfr, although the other bands are already extremely wellconverged, and the energy of the highest bands (whose residual are not yet good enough), is also rather well converged.
In response to this problem, for nonzero nbdbuf , the largest residual (residm), to be later compared with tolwfr, will be computed only in the set of nonbuffer bands (this modification applies for nonselfconsistent as well as selfconsistent calculation, for GS as well as RF calculations). For a GS calculation, with iscf<0, supposing nbdbuf is not initialized in the input file, then ABINIT will overcome the default nbdbuf value, and automatically set nbdbuf to 2.
In metallic RF calculations, in the conjugate gradient optimisation of firstorder wavefunctions, there is an instability situation when the q wavevector of the perturbation brings the eigenenergy of the highest treated band at some kpoint higher than the lowest untreated eigenenergy at some k+q point. If one accepts a buffer of frozen states, this instability can be made to disappear. Frozen states receive automatically a residual value of 0.1 For a RF calculation, with 3 <= occopt <= 7, supposing nbdbuf is not initialized in the input file, then ABINIT will overcome the default nbdbuf value, and automatically set nbdbuf to 2. This value might be too low in some cases.
Also, the number of active bands, in all cases, is imposed to be at least 1, irrespective of the value of nbdbuf .
ndivk¶
Mnemonics: Number of DIVisions of K lines
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_kpoints, topic_ElecBandStructure
Variable type: integer
Dimensions: (abs(kptopt))
Default value: None
Comment: Will be generated automatically from ndivsm if the latter is defined, which is much more convenient.
Only relevant if: kptopt < 0
The use of this variable forbids the use of: ndivsm
Added in version: before_v9
Test list (click to open). Rarely used, [10/1101] in all abinit tests, [3/151] in abinit tutorials
 tutoparal: tdmft_1.abi, tucrpa_1.abi
 tutorial: tfold2bloch_1.abi
 v2: t43.abi
 v4: t89.abi
 v5: t20.abi, t54.abi
 v6: t34.abi
 v7: t76.abi, t77.abi
Gives the number of divisions of each of the segments of the band structure, whose path is determined by kptopt and kptbounds. In this case, the absolute value of kptopt is the number of such segments.
For example, suppose that the number of segment is just one (kptopt = 1), a value ndivk = 4 will lead to the computation of points with relative coordinates 0.0, 0.25, 0.5, 0.75 and 1.0, along the segment in consideration.
Now, suppose that there are two segments (kptopt = 2), with ndivk (1)=4 and ndivk (2)=2, the computation of the eigenvalues will be done at 7 points, 5 belonging to the first segment, with relative coordinates 0.0, 0.25, 0.5, 0.75 and 1.0, the last one being also the starting point of the next segment, for which two other points must be computed, with relative coordinates 0.5 and 1.0.
It is easy to compute disconnected circuits (nonchained segments), by separating the circuits with the value ndivk = 1 for the intermediate segment connecting the end of one circuit with the beginning of the next one (in which case no intermediate point is computed along this segment).
Alternatively it is possible to generate automatically the array ndivk by just specifying the number of divisions for the smallest segment. See the related input variable ndivsm. This is much more convenient than using ndivk .
ndivsm¶
Mnemonics: Number of DIVisions for the SMallest segment
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_kpoints, topic_ElecBandStructure
Variable type: integer
Dimensions: scalar
Default value: None
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [3/151] in abinit tutorials
 libxc: t42.abi
 tutorespfn: tdepes_4.abi, teph4zpr_4.abi
 tutorial: tbase3_5.abi
 v6: t14.abi
 v7: t88.abi
 v8: t04.abi, t44.abi
This variable defines the number of divisions used to sample the smallest segment of the circuit employed in a band structure calculations (see related input variables kptopt and kptbounds). If ndivsm is given in the input file, there is no need to specify the number of divisions to be used for the other segments. Indeed ndivk is automatically calculated inside the code in order to generate a path where the number of divisions in each segment is proportional to the length of the segment itself. This option is activated only when kptopt is negative. In this case, the absolute value of kptopt is the number of such segments.
ngfft¶
Mnemonics: Number of Grid points for Fast Fourier Transform
Mentioned in topic(s): topic_Planewaves
Variable type: integer
Dimensions: (3)
Default value: [0, 0, 0]
Comment: (automatic selection of optimal values)
Added in version: before_v9
Test list (click to open). Moderately used, [66/1101] in all abinit tests, [0/151] in abinit tutorials
Gives the size of fast Fourier transform (FFT) grid in three dimensions. Each number must be composed of the factors 2, 3, and 5 to be consistent with the radices available in our FFT.
If no ngfft is provided or if ngfft is set to 0 0 0, the code will automatically provide an optimal set of ngfft values, based on acell, rprim and ecut (see also boxcutmin for speed/accuracy concerns). This is the recommended procedure, of course. The total number of FFT points is the product: ngfft (1) x ngfft (2) x ngfft (3)=%nfft.
When ngfft is made smaller than recommended values (e.g. by setting boxcutmin to a value smaller than 2.0 or by setting ngfft manually), the code runs faster and the equations in effect are approximated by a low pass Fourier filter. The code reports to standard output (unit 06) a parameter “boxcut” which is the smallest ratio of the FFT box side to the \GG vector basis sphere diameter. When boxcut is less than 2 the Fourier filter approximation is being used. When boxcut gets less than about 1.5 the approximation may be too severe for realistic results and should be tested against larger values of ngfft . When boxcut is larger than 2, ngfft could be reduced without loss of accuracy. In this case, the small variations that are observed are solely due to the xc quadrature, that may be handled with intxc = 1 to even reduce this effect.
Internally, ngfft is an array of size 18. The present components are stored in ngfft (1:3), while
 ngfft (4:6) contains slightly different (larger) values, modified for efficiency of the FFT
 ngfft (7) is fftalg
 ngfft (8) is fftcache
 ngfft (9) is set to 0 if the parallelization of the FFT is not activated, while it is set to 1 if it is activated.
 ngfft (10) is the number of processors of the FFT group
 ngfft (11) is the index of the processor in the group of processors
 ngfft (12) is n2proc, the number of xz planes, in reciprocal space, treated by the processor
 ngfft (13) is n3proc, the number of xy planes, in real space, treated by the processor
 ngfft (14) is mpi_comm_fft, the handle on the MPI communicator in charge of the FFT parallelisation
 ngfft (15:18) are not yet used
The number of points stored by this processor in real space is n1*n2*n3proc, while in reciprocal space, it is n1*n2proc*n3.
ngqpt¶
Mnemonics: Number of Grid points for Q PoinTs generation
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: integer
Dimensions: (3)
Default value: [0, 0, 0]
Only relevant if: nqpt == 1 and kptopt >= 0
The use of this variable forbids the use of: qptrlatt
Added in version: before_v9
Test list (click to open). Moderately used, [31/1101] in all abinit tests, [13/151] in abinit tutorials
 tutorespfn: tdepes_3.abi, tdepes_4.abi, teph_1.abi …
 v6: t78.abi …
 v7: t50.abi, t51.abi, t57.abi …
 v8: t44.abi, t57.abi, t65.abi …
 v9: t53.abi, t54.abi, t55.abi …
WARNING : Only used if nqpt = 1. If nqpt=0 (which is the default value of nqpt), no reading of ngqpt is done, and the default value of ngqpt is kept.
Similar to ngkpt, but for the qwavevector. At variance with ngkpt, note that only one q point is selected per dataset (see iqpt). Its three positive components give the number of q points of MonkhorstPack grids (defined with respect to primitive axis in reciprocal space) in each of the three dimensions. The use of nshiftq and shiftq, allows one to generate shifted grids, or MonkhorstPack grids defined with respect to conventional unit cells.
For more information on MonkhorstPack grids, see ngkpt.
nline¶
Mnemonics: Number of LINE minimisations
Mentioned in topic(s): topic_SCFControl
Variable type: integer
Dimensions: scalar
Default value: 4
Added in version: before_v9
Test list (click to open). Very frequently used, [1098/1101] in all abinit tests, [151/151] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: tdepes_1.abi, tdepes_2.abi, tdepes_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Gives maximum number of line minimizations allowed in preconditioned conjugate gradient minimization for each band. The default, 4, is fine. Special cases, with degeneracies or neardegeneracies of levels at the Fermi energy may require a larger value of nline (5 or 6 ?). Line minimizations will be stopped anyway when improvement gets small (governed by tolrde). With the input variable nnsclo, governs the convergence of the wavefunctions for fixed potential. Note that nline = 0 can be used to diagonalize the Hamiltonian matrix in the subspace spanned by the input wavefunctions.
npsp¶
Mnemonics: Number of PSeudoPotentials to be read
Characteristics: NO_MULTI
Mentioned in topic(s): topic_AtomTypes, topic_PseudosPAW
Variable type: integer
Dimensions: scalar
Default value: ntypat
Added in version: before_v9
Test list (click to open). Moderately used, [23/1101] in all abinit tests, [3/151] in abinit tutorials
Usually, the number of pseudopotentials to be read, npsp , is equal to the number of type of atoms, ntypat. However, in the case an alchemical mixing of pseudopotential is used, npsp will be bigger than ntypat.
Alchemical pseudopotentials will be present when ntypalch is nonzero. See ntypalch to understand how to use alchemical potentials in ABINIT. The input variables (ntypalch, algalch, mixalch) are active, and generate alchemical potentials from the available pseudopotentials. Also, the inner variables (%ntyppure, %npspalch) become active. See these input variables, especially mixalch, to understand how to use alchemical potentials in ABINIT.
npspalch¶
Mnemonics: Number of PSeudoPotentials that are “ALCHemical”
Characteristics: INTERNAL_ONLY
Mentioned in topic(s): topic_AtomTypes
Variable type: integer
Dimensions: scalar
Default value: npsp%ntyppure
Only relevant if: ntypalch>0
Added in version: before_v9
Gives the number of pseudopotentials that are used for alchemical mixing (when ntypalch is nonzero):
When alchemical mixing of potentials is used (that is, when ntypalch>0), then npspalch must be greater than 0.
nqfd¶
Mnemonics: Number of Quasi FermiDirac excited carriers
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_BandOcc, topic_DeltaSCF
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v9: t91.abi
Only used if [[occopt] = 9.
Controls the numbers of electrons per cell constrained in conduction bands with index strictly greater than ivalence, and the number of holes per cell constrained in valence bands with index between 1 and ivalence. Internal variables ne_qFD (number of electrons) and nh_qFD (number of holes) are presently initialized to nqfd , which is NOT INTERNAL.
nqpt¶
Mnemonics: Number of Q  POINTs
Mentioned in topic(s): topic_qpoints
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [204/1101] in all abinit tests, [19/151] in abinit tutorials
 gpu: t01.abi …
 libxc: t81.abi, t82.abi …
 mpiio: t51.abi, t62.abi, t62.abi …
 paral: t52.abi, t53.abi, t52.abi …
 seq: tsv4_80.abi …
 tutoparal: tdfpt_02.abi, tdfpt_04.abi …
 tutorespfn: tdepes_1.abi, tdepes_2.abi, tdepes_3.abi …
 v2: t01.abi, t03.abi, t04.abi …
 v3: t02.abi, t06.abi, t07.abi …
 v4: t02.abi, t56.abi, t58.abi …
 v5: t05.abi, t21.abi, t23.abi …
 v6: t35.abi, t36.abi, t37.abi …
 v67mbpt: t52.abi …
 v7: t41.abi, t43.abi, t45.abi …
 v8: t07.abi, t36.abi, t41.abi …
 v9: t05.abi, t07.abi, t09.abi …
Determines whether one q point must be read (See the variable %qptn). Can be either 0 or 1.
If 1 and used in groundstate calculation, a global shift of all the kpoints is applied, to give calculation at k+q. In this case, the output wavefunction will be appended by _WFQ instead of _WFK Also, if 1 and a RF calculation is done, defines the wavevector of the perturbation.
For further information about the files file, consult the abinit help file.
nshiftq¶
Mnemonics: Number of SHIFTs for Q point grids
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Rarely used, [6/1101] in all abinit tests, [2/151] in abinit tutorials
 tutorespfn: tdepes_3.abi, tdepes_4.abi
 v6: t35.abi, t78.abi
 v7: t50.abi, t51.abi
WARNING : Only used if nqpt = 1. If nqpt=0 (which is the default value of nqpt), no reading of nshiftq is done, and the default value of nshiftq is kept.
This parameter gives the number of shifted grids to be used concurrently to generate the full grid of q points. It can be used with primitive grids defined either from ngqpt or qptrlatt. The maximum allowed value of nshiftq is 8. The values of the shifts are given by shiftq.
nspden¶
Mnemonics: Number of SPinDENsity components
Mentioned in topic(s): topic_spinpolarisation
Variable type: integer
Dimensions: scalar
Default value: nsppol
Added in version: before_v9
Test list (click to open). Moderately used, [97/1101] in all abinit tests, [5/151] in abinit tutorials
 bigdft: t12.abi …
 builtin: testin_libxc.abi …
 gpu: t03.abi, t05.abi, t05.abi …
 libxc: t00.abi, t02.abi, t05.abi …
 mpiio: t28.abi …
 paral: t28.abi, t32.abi, t83.abi …
 seq: tsv3_03.abi …
 tutorial: tdftu_1.abi, tdftu_2.abi, tdftu_3.abi …
 v1: t39.abi …
 v3: t12.abi, t20.abi, t22.abi …
 v4: t99.abi …
 v5: t08.abi, t14.abi, t15.abi …
 v6: t07.abi, t17.abi, t31.abi …
 v7: t05.abi, t07.abi, t12.abi …
 v8: t01.abi, t01_triqs2_0.abi, t25.abi …
 v9: t03.abi, t29.abi …
 wannier90: t04.abi …
If nspden = 1, no spinmagnetization the density matrix is diagonal, with same values spinup and spindown (compatible with nsppol = 1 only, for both nspinor = 1 or 2)
If nspden = 2, scalar magnetization (the axis is arbitrarily fixed in the z direction) the density matrix is diagonal, with different values for spinup and spindown (compatible with nspinor = 1, either with nsppol = 2 general collinear magnetization or nsppol = 1 antiferromagnetism)
If nspden = 4, vector magnetization: the density matrix is full, with allowed x, y and z magnetization (useful only with nspinor = 2 and nsppol = 1, either because there is spinorbit without timereversal symmetry  and thus spontaneous magnetization, or with spinorbit, if one allows for spontaneous noncollinear magnetism). Available for response functions [Ricci2019]. Also note that, with nspden = 4, timereversal symmetry is not taken into account (at present; this has to be checked) and thus kptopt has to be different from 1 or 2.
The default ( nspden = nsppol) does not suit the case of vector magnetization. Note that the choice of nspden has an influence on the treatment of symmetries. See symafm.
nspinor¶
Mnemonics: Number of SPINORial components of the wavefunctions
Mentioned in topic(s): topic_spinpolarisation
Variable type: integer
Dimensions: scalar
Default value: 2 if pawspnorb == 1,
1 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [85/1101] in all abinit tests, [5/151] in abinit tutorials
 builtin: testin_libxc.abi …
 gpu: t03.abi …
 libxc: t00.abi, t02.abi, t05.abi …
 mpiio: t28.abi …
 paral: t28.abi, t32.abi, t83.abi …
 tutorial: tdftu_1.abi, tdftu_2.abi, tdftu_3.abi …
 v1: t01.abi, t39.abi …
 v2: t74.abi, t75.abi, t76.abi …
 v3: t01.abi, t02.abi, t13.abi …
 v4: t36.abi …
 v5: t16.abi, t17.abi, t20.abi …
 v6: t07.abi, t17.abi …
 v7: t03.abi, t05.abi, t06.abi …
 v8: t01.abi, t01_triqs2_0.abi, t26.abi …
 v9: t03.abi, t29.abi …
If nspinor = 1, usual case: scalar wavefunction (compatible with (nsppol = 1, nspden = 1) as well as (nsppol = 2, nspden = 2) )
If nspinor = 2, the wavefunction is a spinor (compatible with nsppol = 1, with nspden = 1 or 4, but not with nsppol = 2)
When nspinor is 2, the values of istwfk are automatically set to 1. Also, the number of bands, for each kpoint, should be even.
ntypalch¶
Mnemonics: Number of TYPe of atoms that are “ALCHemical”
Mentioned in topic(s): topic_AtomTypes
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [0/151] in abinit tutorials
Used for the generation of alchemical pseudopotentials: when ntypalch is nonzero, alchemical mixing will be used.
Among the ntypat types of atoms, the last ntypalch will be “alchemical” pseudoatoms, while only the first %ntyppure will be uniquely associated with a pseudopotential (the %ntyppure first of these, actually). The ntypalch types of alchemical pseudoatoms are to be made from the remaining %npspalch pseudopotentials.
In this case, the input variables algalch, mixalch are active, and generate alchemical potentials from the available pseudopotentials. See these input variables, especially mixalch, to understand how to use alchemical potentials in ABINIT.
ntyppure¶
Mnemonics: Number of TYPe of atoms that are “PURE”
Characteristics: INTERNAL_ONLY
Mentioned in topic(s): topic_AtomTypes
Variable type: integer
Dimensions: scalar
Default value: ntypatntypalch
Added in version: before_v9
Gives the number of type of atoms that are “pure” when alchemical mixing is used (ntypalch /= 0):
nucdipmom¶
Mnemonics: NUClear DIPole MOMents
Mentioned in topic(s): topic_MagField
Variable type: real
Dimensions: (3,natom)
Default value: 0.0
Only relevant if: %usepaw = 1; pawcpxocc = 2; kptopt > 2
Added in version: before_v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [0/151] in abinit tutorials
Places an array of nuclear magnetic dipole moments on the atomic positions, useful for computing the magnetization in the presence of nuclear dipoles and thus the chemical shielding by the converse method [Thonhauser2009]. The presence of these dipoles breaks time reversal symmetry and lowers the overall spatial symmetry. The dipole moment values are entered in atomic units. For reference, note that one Bohr magneton has value 1/2 in atomic units, while one nuclear Bohr magneton has value 2.7321\times 10^{4} in atomic units.
nwfshist¶
Mnemonics: Number of WaveFunctionS HISTory
Mentioned in topic(s): topic_Wavelets
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [16/1101] in all abinit tests, [9/151] in abinit tutorials
 bigdft: t02.abi, t06.abi, t07.abi, t10.abi, t11.abi, t12.abi, t14.abi
 tutoparal: tgswvl_1.abi, tgswvl_1.abi, tgswvl_1.abi, tgswvl_1.abi, tgswvl_1.abi, tgswvl_1.abi, tgswvl_2.abi, tgswvl_2.abi, tgswvl_2.abi
In the wavelet basis set, the ground state is found by direct minimisation. The algorithm used can be either the steepest descent or the DIIS (Direct Inversion of Iteration Space). When nwfshist = 0, the steepest descent is used ( i.e. there is no history storage of the previous iterations).
If nwfshist is strictly positive, a DIIS is used. A typical value is 6. Using a DIIS increases the memory required by the program since N previous wavefunctions are stored during the electronic minimisation.
occ¶
Mnemonics: OCCupation numbers
Characteristics: EVOLVING
Mentioned in topic(s): topic_BandOcc
Variable type: real
Dimensions: (nband,nsppol)
Commentdims: in case occopt=2, dimensions are (%mband,nkpt,nsppol)
Default value: 0
*Added in version: before_v9
Test list (click to open). Very frequently used, [616/1101] in all abinit tests, [75/151] in abinit tutorials
 atompaw: t04.abi …
 bigdft: t11.abi, t12.abi, t14.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t09.abi, t21.abi …
 fast: t00.abi, t02.abi, t03.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t02.abi, t03.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t03.abi, t03.abi, t03.abi …
 psml: t02.abi, t03.abi, t04.abi …
 seq: tsv2_81.abi, tsv2_82.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_04.abi …
 tutorespfn: telast_6.abi, teph_1.abi, teph4mob_5.abi …
 tutorial: tbase1_2.abi, tbase1_3.abi, tbase1_5.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t06.abi, t09.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t03.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t04.abi …
 v8: t01.abi, t01_triqs2_0.abi, t04.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t04.abi …
Gives occupation numbers for all bands in the problem. Needed if occopt == 0 or occopt == 2. Ignored otherwise (automatically computed). Also ignored when iscf = 2. Typical band occupancy is either 2.0 or 0.0, but will usually be 1.0 or 0.0 for nsppol=2, or nspinor=2, or halfoccupied band, or other choices in special circumstances.
If occopt is not 2, then the occupancies must be the same for each k point. If nsppol=1, the total number of arrays which must be provided is nband, in order of increasing energy. If nsppol=2, the total number of arrays which must be provided is nband*nsppol, first spin up, in order of increasing electronic eigenenergy, then spin down, in order of increasing electronic eigenenergy.
If occopt = 2, then the band occupancies must be provided explicitly for each band, EACH k POINT, and EACH SPINPOLARIZATION, in an array which runs over all bands, k points, and spinpolarizations. The order of entries in the array would correspond to all bands at the first k point (spin up), then all bands at the second k point (spin up), etc, then all kpoints spin down. The total number of array elements which must be provided is ( nband(1)+nband(2)+…+ nband(nkpt) ) * nsppol.
The occupation numbers evolve only for metallic occupations, that is, occopt â‰¥ 3.
When there are several images, occ might depend on the image number, see the description in nimage.
optdriver¶
Mnemonics: OPTions for the DRIVER
Mentioned in topic(s): topic_nonlinear, topic_GWls, topic_ElPhonInt, topic_GW, topic_BSE, topic_DFPT, topic_Susceptibility, topic_SelfEnergy
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [165/1101] in all abinit tests, [31/151] in abinit tutorials
 libxc: t41.abi, t42.abi, t43.abi …
 paral: t54.abi, t54.abi, t54.abi …
 tutoparal: tmbt_2.abi, tmbt_3.abi, tucrpa_2.abi …
 tutorespfn: teph4mob_5.abi, teph4mob_6.abi, teph4mob_7.abi …
 tutorial: tbs_1.abi, tbs_2.abi, tbs_3.abi …
 v3: t30.abi, t31.abi, t83.abi …
 v4: t52.abi, t84.abi, t85.abi …
 v5: t63.abi, t64.abi, t65.abi …
 v6: t63.abi, t64.abi, t65.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t12.abi, t16.abi, t23.abi …
 v8: t44.abi, t47.abi, t57.abi …
 v9: t31.abi, t32.abi, t33.abi …
 wannier90: t03.abi …
For each dataset, choose the task to be done, at the level of the “driver” routine.
The choice is among:
 0 → groundstate calculation (GS), routine gstate
 1 → responsefunction calculation (RF), routine respfn
 2 → susceptibility calculation (SUS), routine suscep
 3 → susceptibility and dielectric matrix calculation (SCR), routine screening
 4 → selfenergy calculation (SIG), routine sigma.
 5 → nonlinear response functions (NONLINEAR), using the 2n+1 theorem, routine nonlinear.
 7 → electronphonon coupling (EPH), see also eph_task input variable.
 8 → Postprocessing of WFK file, routine wfk_analyze. See also wfk_task input variable.
 10 → longwave response functions (LONGWAVE), routine longwave. See also lw_flexo or lw_qdrpl input variables.
 66 → GW using LanczosSternheimer, see input variables whose name start with
gwls_*
.  99 → BetheSalpeter calculation (BSE), routine bethe_salpeter
If one of rfphon, rfddk, rfelfd, or rfstrs is nonzero, while optdriver is not defined in the input file, ABINIT will set optdriver to 1 automatically. These input variables (rfphon, rfddk, rfelfd, and rfstrs) must be zero if optdriver is not set to 1.
optstress¶
Mnemonics: OPTion for the computation of STRESS
Mentioned in topic(s): topic_ForcesStresses
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Moderately used, [45/1101] in all abinit tests, [16/151] in abinit tutorials
 bigdft: t01.abi, t02.abi, t03.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 libxc: t70.abi …
 seq: tsv7_70.abi …
 tutoparal: tgswvl_1.abi, tgswvl_1.abi, tgswvl_1.abi …
 tutorial: tpositron_1.abi, tpositron_2.abi, tpositron_3.abi …
 v5: t07.abi …
 v7: t26.abi, t32.abi, t35.abi …
 v9: t44.abi, t45.abi, t73.abi …
If set to 1, the computation of stresses is done, in the SCF case (under the conditions iscf > 0, prtstm == 0, positron == 0, and either nstep >0, or %usepaw == 0 or irdwfk == 1). Otherwise, to save CPU time, if no optimization of the cell is required, one can skip the computation of stresses. The CPU time saving might be interesting for some PAW calculations.
orbmag¶
Mnemonics: ORBital MAGnetization
Characteristics: DEVELOP
Mentioned in topic(s): topic_MagField
Variable type: integer
Dimensions: scalar
Default value: 0
Only relevant if: %usepaw == 1;
usexcnhat == 0;
nspinor == 1;
paral_atom == 0;
paral_kgb == 0;
(kptopt == 3 or kptopt == 0)
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
Compute quantities related to orbital magnetic moment. The implementation assumes an insulator, so no empty or partially filled bands, and currently restricted to nspinor 1. Such insulators have orbital magnetization zero, except in the presence of nonzero nuclear dipole moments, see nucdipmom. orbmag is parallelized over k points only. The implementation follows the theory outlined in [Gonze2011a] extended to the PAW case; see also [Ceresoli2006]. The computed results are returned in the standard output file, search for “Orbital magnetic moment”. This calculation requires both the ground state and DDK wavefunctions, and is triggered at the end of a DDK calculation.
 orbmag = 1: Compute orbital magnetization and integral of the Berry curvature (Chern number) over the Brillouin zone.
 orbmag = 2: Same as orbmag 1 but also print out values of each term making up total orbital magnetic moment.
 orbmag = 3: Same as orbmag 2 but print out values of each term for each band.
posdoppler¶
Mnemonics: POSitron computation of DOPPLER broadening
Mentioned in topic(s): topic_positron
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [2/151] in abinit tutorials
 tutorial: tpositron_5.abi, tpositron_7.abi
 v7: t35.abi
Relevant only when positron/=0. This input parameter activates the calculation of the Doppler broadening of the electronpositron annihilation radiation. An output file containing the momentum distributions of annihilating electronpositron pairs is created. Such a computation needs a core wavefunction file (per atom type) to be provided. This core WF file should be named ‘psp_file_name.corewf’ (where pspfile_name is the name of the pseudopotential (or PAW) file) or ‘corewf.abinit**ityp**’ (where ityp is the index of the atom type). Core WF files can be obtained with the atompaw tool by the use of prtcorewf keyword.
positron¶
Mnemonics: POSITRON calculation
Mentioned in topic(s): topic_positron
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [11/1101] in all abinit tests, [7/151] in abinit tutorials
 tutorial: tpositron_1.abi, tpositron_2.abi, tpositron_3.abi, tpositron_4.abi, tpositron_5.abi, tpositron_6.abi, tpositron_7.abi
 v4: t95.abi, t96.abi
 v5: t07.abi
 v7: t35.abi
This input parameter can be positive or negative. Negative values for positron are only relevant for PAW calculations. Electronpositron correlation functional is defined by ixcpositron. Other relevant input parameter: posocc (occupation number for the positron).
Positive values for positron : For positron = 1 or 2, will perform the calculation of positron lifetime (and annihilation rate).
 positron = 1:\n Starting from a previous electronic GS density (with positron = 0), a positronic groundstate calculation is performed, considering that the electrons are not perturbed by the presence of the positron. This is almost correct for a positron in a perfect bulk material. But this approximation fails when defects exist in the material (for instance: the positron might be trapped by a vacancy). The electronic density will be automatically read from a _DEN file (with or without irdden keyword). At the end of the SCF cycle, the positron lifetime and annihilation rate are printed out.
Additional information for the use of pseudopotentials:
PAW datasets: nothing to do; simply use usual electronic PAW datasets
Normconserving pseudopotentials: One has to use specific pseudopotentials for the positron calculation. They must be of the FHI type (pspcod=6), and must contain at their end, the allelectrons core density generated with FHI98PP. They must have lmax=lloc=0 (check that this works for the electronic GS !! No ghost, etc...). Otherwise, their are similar to an usual FHI pseudopotential.
 positron = 2:\n Starting from a previous positronic GS density (with positron = 1 ), an electronic groundstate calculation is performed, keeping the positronic density constant. The positronic density will be automatically read from a _DEN file (with or without getden/irdden keyword). At the end of the SCF cycle, the positron lifetime and annihilation rate are printed out.
Additional information for the use of pseudopotentials:
PAW datasets: nothing to do; simply use usual electronic PAW datasets
Normconserving pseudopotentials: One has to use specific pseudopotentials for the electron calculation. They must be of the FHI type (pspcod=6), and must contain at their end, the allelectrons core density generated with FHI98PP.
 Typical use:\n The calculation is done in several steps: The first one is a normal GS calculation for the electrons, with positron = 0. The only specific thing to do is to set prtden = 1 (this is the default for ABINIT v6.x+). This will create the associated _DEN file which will be used as input file for the positronic GS calculation. The second step is the GS calculation of the positron and subsequently its lifetime, with positron =1. One has to define also ixcpositron. Then, it is possible to perform an additional step, computing the GS electronic density in presence of the positron, with positron = 2 and so on… This procedure can be automated (for PAW only) by the use of a negative value for positron . At the end, a converged value of the positron lifetime (decomposed in several contributions) is printed. See also posdoppler keyword for the calculation of Doppler broadening.
Negative values for positron :\n For positron <0, will perform an automatic calculation of electrons and positron densities in the twocomponent DFT context; then will compute positron lifetime (and annihilation rate).

positron = 1:\n Starting from scratch, will first perform a usual electronic groundstate calculation until convergence (controlled by the use of one of the tolerance keywords). Then will perform a positronic ground state calculation in presence of the electrons and ions; then an electronic ground state calculation in presence of the positron and the ions and so on until the total energy is converged. The convergence of the total energy of the ions+electrons+positron system is controlled by the use of the postoldfe, postoldff and posnstep input keywords. With positron = 1, at the beginning of each new electronic/positronic step, the wave functions are unknown.

positron = 10:\n Same as positron = 1 except that the electronic/positronic wave functions are stored in memory. Consequently, the total number of iterations to reach the convergence (\DeltaEtotal<postoldfe or \DeltaForces<postoldff) is smaller. But, this can increase the total amount of memory needed by the code.

positron = 2:\n Same as positron = 1 except that the twocomponent DFT cycle is forced to stop at the end of an electronic step.

positron = 20:\n Same as positron = 10 except that the twocomponent DFT cycle is forced to stop at the end of an electronic step.
Advice for use: There are two typical cases which have to be differently treated:

A positron in a perfect bulk system:\n In that case, the positron is delocalized in the whole crystal. Its density is almost zero. Thus, the “zero density positron limit” has to be used. ixcpositron has to be chosen accordingly. In order to have the zero density positron limit it is advised to follow these points:

1 Put a small positronic charge (by setting a posocc to a small value) OR use a big supercell.
 2 Use only k=\Gamma wave vector for the positronic calculation.
 3 Use the manual procedure in 2 steps: first positron = 0 and then positron = 1; avoid the positron = 2 step and the automatic procedure ( positron <0).
In principle, the positron lifetime should converge with the value of posocc or the size of the supercell.
 A positron trapped in a default (vacancy):\n In that case, the positron is localized in the default. Its density can be localized in the simulation cell (provided that the cell is sufficiently large) and influences the electronic density. So, it is advised to use the automatic procedure ( positron <0) or the manual procedure with several positron = 0,1,2,1,… steps. Kpoints can be used as in usual electronic calculations. Also note that it is possible to use forces and stresses to perform structural minimization.
References: [Arponen1979a], [Boronski1986], [Sterne1991], [Puska1995], [Barbiellini1995]
posnstep¶
Mnemonics: POSitron calculation: max. Number of STEPs for the twocomponent DFT
Mentioned in topic(s): topic_positron
Variable type: integer
Dimensions: scalar
Default value: 50
Added in version: before_v9
Test list (click to open). Rarely used, [6/1101] in all abinit tests, [4/151] in abinit tutorials
 tutorial: tpositron_3.abi, tpositron_4.abi, tpositron_5.abi, tpositron_7.abi
 v5: t07.abi
 v7: t35.abi
Relevant only when positron<0. Sets the maximum number of electronic/positronic iterations that, when reached, will cause the twocomponent DFT SCF cycle to stop. The code will first compute the electronic groundstate, then the positronic ground state in the electronic density, then the electronic ground state in the positronic density until \DeltaEtotal<postoldfe or \DeltaForces<postoldff or the number of electronic/positronic steps is posnstep .
posocc¶
Mnemonics: POSitron calculation: OCCupation number for the positron
Mentioned in topic(s): topic_positron
Variable type: real
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Rarely used, [9/1101] in all abinit tests, [7/151] in abinit tutorials
 tutorial: tpositron_1.abi, tpositron_2.abi, tpositron_3.abi, tpositron_4.abi, tpositron_5.abi, tpositron_6.abi, tpositron_7.abi
 v5: t07.abi
 v7: t35.abi
Relevant only when positron/=0. Sets the occupation number for the positron. Has to be <=1. Changing posocc is only useful for bulk calculation when one wants to perform lifetime computations using a small simulation cell (can avoid the use of a supercell). It simulates the dispersion of the positron in the whole crystal.
postoldfe¶
Mnemonics: POSitron calculation: TOLerance on the DiFference of total Energy
Characteristics: ENERGY
Mentioned in topic(s): topic_positron
Variable type: real
Dimensions: scalar
Default value: 1e06 if postoldff = 0,
0.0 otherwise.
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [2/151] in abinit tutorials
 tutorial: tpositron_3.abi, tpositron_4.abi
 v5: t07.abi
 v7: t35.abi
Relevant only when positron<0. Sets a tolerance for absolute difference of total energy (of ions+electrons+positron system) that, when reached, will cause the SCF cycle to stop before the number of steps is nstep or the number of electronic/positronic steps is posnstep.
Can be specified in Ha (the default), Ry, eV or Kelvin, since postoldfe has the ENERGY characteristics. One and only one of postoldfe or postoldff can be set.
postoldff¶
Mnemonics: POSitron calculation: TOLerance on the DiFference of Forces
Mentioned in topic(s): topic_positron
Variable type: real
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v5: t07.abi
Relevant only when positron < 0. Sets a tolerance for absolute difference of maximum force (in hartree/Bohr) acting on ions (due to ions+electrons+positron system) that, when reached, will cause the SCF cycle to stop before the number of SCF steps is nstep or the number of electronic/positronic steps is posnstep. One and only one of postoldfe or postoldff can be set.
prtdensph¶
Mnemonics: PRinT integral of DENsity inside atomic SPHeres
Mentioned in topic(s): topic_printing
Variable type: integer
Dimensions: scalar
Default value: 1 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [16/1101] in all abinit tests, [0/151] in abinit tutorials
When this flag is activated, values of integral(s) of total density inside sphere(s) around each atom are printed in output file (for each spin component). Spheres around atoms are defined by a radius given by ratsph keyword. Note: integral of density inside a sphere around an atom can be used to determine a rough approximation of the local magnetic moment; this is particularly useful for antiferromagnetic systems. The algorithm to compute this integral is particularly primitive: the points on the FFT grids, belonging to the interior of the sphere are determined, and the value of the functions on these points are summed, taking into account a fixed volume attributed to each point. In particular, the integral as a function of the radius will be a constant, except when a new point enters the sphere, in which case a sudden jump occurs. However, since the purpose of this output is to get a rough idea of the repartition of the density, this is not a real problem. If you are interested in a more accurate estimation of the density within a sphere, you should use the cut3d postprocessor (cut3d help file).
prtebands¶
Mnemonics: PRinT Electron BANDS
Mentioned in topic(s): topic_printing
Variable type: integer
Dimensions: scalar
Default value: 0 if nimage > 1,
1 otherwise.
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [0/151] in abinit tutorials
This option activates the output of the electron eigenvalues. Possible values:
 0 Disable the output of the band energies.
 1 Write eigenvalues in xmgrace format. A file with extension
EBANDS.agr
is produced at the end of the run. Usexmgrace file_EBANDS.agr
to visualize the band energies  2 Write eigenvalues in gnuplot format. The code produces a
EBANDS.dat
file with the eigenvalues and afile_EBANDS.gnuplot
script. Usegnuplot file_EBANDS.gnuplot
to visualize the band energies.
qpt¶
Mnemonics: Q PoinT
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: real
Dimensions: (3)
Default value: [0, 0, 0]
Added in version: before_v9
Test list (click to open). Moderately used, [224/1101] in all abinit tests, [29/151] in abinit tutorials
 gpu: t01.abi …
 libxc: t81.abi, t82.abi …
 mpiio: t51.abi, t62.abi, t62.abi …
 paral: t52.abi, t53.abi, t52.abi …
 seq: tsv4_80.abi …
 tutoparal: tdfpt_02.abi, tdfpt_04.abi …
 tutorespfn: tdepes_1.abi, tdepes_2.abi, tdepes_3.abi …
 v2: t01.abi, t03.abi, t04.abi …
 v3: t02.abi, t06.abi, t07.abi …
 v4: t02.abi, t56.abi, t58.abi …
 v5: t05.abi, t21.abi, t23.abi …
 v6: t35.abi, t36.abi, t37.abi …
 v67mbpt: t52.abi …
 v7: t41.abi, t43.abi, t45.abi …
 v8: t07.abi, t36.abi, t41.abi …
 v9: t05.abi, t07.abi, t09.abi …
WARNING : Only used if nqpt = 1. If nqpt=0 (which is the default value of nqpt), no reading of qpt is done, and the default value of qpt , namely the null vector, is kept.
Combined with qptnrm, define the q vector %qptn(1:3) in the case qptopt = 0.
This input variable is not internal (%qptn(1:3) is used instead), but is used to echo the value of %qptn(1:3), with renormalisation factor one.
qptnrm¶
Mnemonics: Q PoinTs NoRMalization
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: real
Dimensions: scalar
Default value: 1.0
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
Only used if nqpt = 1 and qptopt = 0
Provides renormalization of qpt. Must be positive, nonzero. The actual q vector (renormalized) is %qptn(1:3)= qpt(1:3)/ qptnrm .
qptopt¶
Mnemonics: QPoinTs OPTion
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [13/1101] in all abinit tests, [3/151] in abinit tutorials
 paral: t60.abi, t60.abi, t60.abi
 tutorespfn: tdepes_3.abi, tdepes_4.abi, teph_1.abi
 v6: t35.abi, t78.abi
 v7: t50.abi, t51.abi, t55.abi, t57.abi
 v8: t65.abi
WARNING : Only used if nqpt = 1. If nqpt=0 (which is the default value of nqpt), no reading of qptopt is done, and the default value of qptopt is kept.
Controls the set up to generate the Q point %qptn(1:3) to be used for the specific dataset, either as a shift of kpoint grid in groundstate calculations, or as a standalone phonon wavevector.
There are two basic techniques to generate the Q point: either by specifying it directly, possibly with a renormalisation factor ( qptopt = 0), or extracting it from a grid a Q points ( qptopt = 1 to 4), using the index iqpt. At variance with the similar generation of k points, only ONE q point can be used per dataset.
With qptopt = 1 to 4, rely on ngqpt or qptrlatt, as well as on nshiftq and shiftq to set up a q point grid, from which the q point with number iqpt will be selected. The values qptopt = 1 to 4 differ by the treatment of symmetries. Note that the symmetries are recomputed starting from the values of %rprimd, xred and spinat. So, the explicit value of symrel are not used. This is to allow doing calculations with nsym = 1, sometimes needed for Tdependent electronic structure, still decreasing the number of q points in the case qptopt = 1 or qptopt = 3.

0 → read directly qpt, and its (possible) renormalisation factor qptnrm.

1 → Take fully into account the symmetry to generate the grid of q points in the Irreducible Brillouin Zone only. (This is the usual mode for RF calculations)

2 → Take into account only the timereversal symmetry: q points will be generated in half the Brillouin zone.

3 → Do not take into account any symmetry: q points will be generated in the full Brillouin zone.

4 → Take into account all the symmetries EXCEPT the timereversal symmetry to generate the k points in the Irreducible Brillouin Zone.
In the case of a grid of q points, the auxiliary variables kptrlen, ngkpt and prtkpt might help you to select the optimal grid, similarly to the case of the K point grid.
qptrlatt¶
Mnemonics: Q  PoinTs grid: Real space LATTice
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: integer
Dimensions: (3,3)
Default value: 0
*The use of this variable forbids the use of: ngqpt
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v6: t35.abi
WARNING : Only used if nqpt = 1. If nqpt=0 (which is the default value of nqpt), no reading of qptrlatt is done, and the default value of qptrlatt is kept.
This input variable is used only when qptopt is positive. It partially defines the q point grid. The other piece of information is contained in shiftq. qptrlatt cannot be used together with ngqpt.
The values qptrlatt (1:3,1), qptrlatt (1:3,2), qptrlatt (1:3,3) are the coordinates of three vectors in real space, expressed in the %rprimd coordinate system (reduced coordinates). They defines a superlattice in real space. The k point lattice is the reciprocal of this superlattice, possibly shifted (see shiftq).
If neither ngqpt nor qptrlatt are defined, ABINIT will automatically generate a set of k point grids, and select the best combination of qptrlatt and shiftq that allows one to reach a sufficient value of kptrlen. See this latter variable for a complete description of this procedure.
ratsm¶
Mnemonics: Radii of the ATomic spheres SMearing
Mentioned in topic(s): topic_printing, topic_MagMom, topic_ConstrainedDFT
Variable type: real
Dimensions: scalar
Default value: 0.05 if any(constraint_kind > 1),
0.0 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [12/1101] in all abinit tests, [0/151] in abinit tutorials
Smearing width for the atomic spheres whose radius is determined by ratsph. For each spherical zone around each atom, the integrating function goes from 1.0 to 0.0 in an interval from ratsph ratsm to ratsph. The function is the same as the one used to smear the kinetic energy, see ecutsm.
ratsph¶
Mnemonics: Radii of the ATomic SPHere(s)
Mentioned in topic(s): topic_printing, topic_MagMom, topic_ElecBandStructure, topic_ElecDOS, topic_ConstrainedDFT
Variable type: real
Dimensions: (ntypat)
Default value: AUTO_FROM_PSP if %usepaw == 1,
2.0 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [19/1101] in all abinit tests, [2/151] in abinit tutorials
Relevant only when prtdensph = 1, or magconon/=0, or any(constraint_kind(:)/=0) (that is, constrained DFT), or prtdos = 3. In most cases (see later for prtdos = 3), provides the radius of the spheres around each atom in which the total charge density or magnetization will be integrated. The integral within the sphere is obtained by a sum over real space FFT points inside the sphere, multiplied by a function that is one inside the sphere, except in a small boundary zone determined by ratsm, where this fonction goes smoothly from 1 to 0. In case of PAW, ratsph radius has to be greater or equal to PAW radius of considered atom type (which is read from the PAW dataset file; see rc_sph or r_paw). In case of constrained DFT, note that the sphere for different atoms are not allowed to overlap.
When prtdos = 3 or 4 :
Provides the radius of the spheres around the natsph atoms of indices iatsph, in which the local DOS and its angularmomentum projections will be analysed.
The choice of this radius is quite arbitrary. In a planewave basis set, there is no natural definition of an atomic sphere. However, it might be wise to use the following welldefined and physically motivated procedure: from the Bader analysis, one can define the radius of the sphere that contains the same charge as the Bader volume. This “Equivalent Bader charge atomic radius” might then be used to perform the present analysis. See the aim help file for more explanations. Another physically motivated choice would be to rely on another charge partitioning, like the Hirshfeld one (see the cut3d utility cut3d help file). The advantage of using charge partitioning schemes comes from the fact that the sum of atomic DOS, for all angular momenta and atoms, integrated on the energy range of the occupied states, gives back the total charge. If this is not an issue, one could rely on the half of the nearest neighbour distances, or any scheme that allows one to define an atomic radius. Note that the choice of this radius is however critical for the balance between the s, p and d components. Indeed, the integrated charge within a given radius, behave as a different power of the radius, for the different channels s, p, d. At the limit of very small radii, the s component dominates the charge contained in the sphere.
ratsph_extra¶
Mnemonics: Radii of the ATomic SPHere(s) in the EXTRA set
Characteristics: LENGTH
Mentioned in topic(s): topic_printing
Variable type: real
Dimensions: scalar
Default value: 2.0 Bohr
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v7: t04.abi
Radius for extra spheres the DOS is projected into. See natsph_extra and xredsph_extra for the number and positions of the spheres.
rcut¶
Mnemonics: Radius of the CUToff for coulomb interaction
Mentioned in topic(s): topic_Coulomb, topic_GWls, topic_Susceptibility, topic_SelfEnergy
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Moderately used, [16/1101] in all abinit tests, [0/151] in abinit tutorials
Truncation of the Coulomb interaction in real space. The meaning of rcut is governed by the cutoff shape options icutcoul, gw_icutcoul and/or fock_icutcoul.
If rcut is negative, the cutoff is automatically calculated so to enclose the same volume inside the cutoff as the volume of the primitive cell.
rmm_diis¶
Mnemonics: Activate the RMMDIIS eigensolver for GS calculations.
Mentioned in topic(s): topic_TuningSpeedMem, topic_SCFAlgorithms
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: 9.3.0
Test list (click to open). Rarely used, [5/1101] in all abinit tests, [0/151] in abinit tutorials
This variable activates the Residual Minimization Method, Direct Inversion in the Iterative Subspace (RMMDIIS) eigensolver to accelerate GS computations, structural relaxations and moleculardynamics simulations. The algorithm is inspired to [Kresse1996] although the ABINIT implementation is not completely equivalent to the original formulation. RMMDIIS can be used both with NC and PAW pseudos and is fully compatible with the paral_kgb distribution. This variable has no effect when optdriver > 0.
It is worth noting that RMMDIIS is usually employed in conjunction with another eigenvalue solver that is used to perform the initial SCF iterations. RMMDIIS, indeed, is not guaranteed to converge to the ground state as the algorithm will find the eigenvector/eigenvalue pair closest to the input trial state. The reliability of RMMDIIS thus strongly depends on the initialisation of the SCF cycle and on the quality of the input wavefunctions that are supposed to be reasonably close to the final solution.
Both the CG and the LOBPCG solver can be used in conjunction with RMMDIIS although it is strongly suggested to use LOBPCG with paral_kgb = 1 to take advantage of its better efficiency and improved parallel scalability. To activate RMMDIIS with e.g. LOBPCG, it is sufficient to use:
paral_kgb 1
rmm_diis 1
in the input file and then select the value of npkpt, npband, npfft, npspinor according to the system.
Tip
Using use_gemm_nonlop = 1 with bandpp > 1 can lead to a significant speedup at the price of increased memory requirements.
In the case of standard SCF calculations, Abinit activates the RMMDIIS solver after 3 + rmm_diis SCF iterations. When performing structural relaxations, RMMDIIS is activated after rmm_diis SCF iterations once the first GS calculation is completed. This means that using rmm_diis = 1 for a structural relaxation leads to:
 4 SCF iterations with the CG/LOBPCG eigensolver followed by RMMDIIS when are performing the **first GS calculation**.
 1 SCF iterations with CG/LOBPCG followed by RMMDIIS for all the subsequent relaxation steps.
A negative value rmm_diis (e.g. 3) can be used to bypass the initial CG/LOBPCG iterations but this option should be used with extreme care and it is not recommended in general.
RMMDIIS usually requires less walltime per iteration when compared to other approaches since there is no explicit orthogonalization while optimizing the trial states. Only a single fullband Cholesky orthogonalization is performed per SCF iteration before recomputing the new density. As a consequence, one RMMDIIS iteration is usually faster (sometimes even by a factor two) than one CG/LOBPCG iteration, especially in systems with relatively large %mpw. However, the additional steps of the algorithm (subspace rotation and Cholesky orthogonalization) present poor MPIscalability hence this part will start to dominate the walltime in systems with large nband.
The algorithm can also be used for NSCF band structure calculations although one should not expect RMMDIIS to provide highenergy states of the same quality as the one obtained with other eigenvalue solvers. Although RMMDIIS can be used for computing band structures and electron DOS with iscf = 2, we do not recommend using this solver to produce WFK files with many empty states as required by manybody calculations.
From the above discussion, it should be clear that RMMDIIS is less stable than the other eigenvalue solvers and you should be aware (and accept) that RMMDIIS calculations may fail. In the best case scenario, one clearly sees the problem as the SCF cycle does not convergence and/or the code aborts with a nonzero exit status returned by one of the LAPACK routines. In the worst case scenario, RMMDIIS may seamlessly converge to the wrong solution. This is especially true if the first LOBPCG/CG iterations fail to provide a reasonablygood initial starting point.
In problematic cases, you may want to increase the value of rmm_diis so that more SCF iterations are done with LOBPCG/CG before activating RMMDIIS. RMMDIIS is indeed very sensitive to density mixing and increasing rmm_diis is a possible solution if you don’t want to play with the variables defining the mixing algorithm.
As a general rule, we strongly suggest to include in the calculation more nband bands that what is strictly needed to accommodate all the valence electrons (let’s say > 10% of the occupied bands). Increasing nband, indeed, increases the number of vectors used for the RayleighRitz subspace rotation and this step is crucial for finding the most accurate variational approximation to the eigenvectors before starting the DIIS optimization.
For a given precision, RMMDIIS usually requires more iterations than the other eigensolvers. For performance reasons, one should avoid using tight tolerances. Something of the order of tolvrs = 1e8 or toldfe = 1e10 to stop the SCF cycle should be fine. Avoid using (tolwfr) (criterion on the residuals) as converge criterion for SCF cycles Use tolwfr only if you are using RMMDIIS for NSCF band structure calculations (as this is the only converge criterion available for NSCF calculations).
Last but not least, note that the default implementation is quite aggressive at the level of memory allocations. A less memoryintensive version of the algorithm can be activated via rmm_diis_savemem = 1.
rmm_diis_savemem¶
Mnemonics: Save memory in the RMMDIIS eigensolver.
Mentioned in topic(s): topic_TuningSpeedMem
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: 9.3.0
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
This variable allows one to save memory in the RMMDIIS eigensolver when rmm_diis /= 0. By default, the RMMDIIS routine allocates two extra arrays to reduce the number of applications of the Hamiltonian. The size of these arrays depends on the number of planewaves treated by each processor and nband. The amount of memory scales with npband and [[npfft] yet this extra memory is not negligible and the code may go out of memory for large systems. In this case, one can use rmm_diis_savemem = 1 to activate a version of RMMDIIS that avoids these extra allocations.
scphon_supercell¶
Mnemonics: Self Consistent PHONon SUPERCELL
Mentioned in topic(s): topic_DFPT
Variable type: integer
Dimensions: (3)
Default value: [1, 1, 1]
Added in version: before_v9
Test list (click to open). Rarely used, [0/1101] in all abinit tests, [0/151] in abinit tutorials
Give extent, in number of primitive unit cells, of the supercell being used for a selfconsistent phonon calculation. Presumes the phonon frequencies and eigenvectors have been calculated in the original primitive unit cell, on a grid of qpoints which corresponds to the supercell in the present calculation. Experimental.
scphon_temp¶
Mnemonics: Self Consistent PHONon TEMPerature
Characteristics: ENERGY
Mentioned in topic(s): topic_DFPT
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [0/1101] in all abinit tests, [0/151] in abinit tutorials
Temperature which is imposed on phonon distribution, in the selfconsistent scheme of [Souvatzis2008]. Determines the extent of the finite displacements used, and consequent anharmonic effects. Experimental.
shiftq¶
Mnemonics: SHIFT for Q points
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_qpoints
Variable type: real
Dimensions: (3,nshiftq)
Default value: None if nshiftq>1,
[0.5, 0.5, 0.5] otherwise.
Added in version: before_v9
Test list (click to open). Rarely used, [8/1101] in all abinit tests, [2/151] in abinit tutorials
 tutorespfn: tdepes_3.abi, tdepes_4.abi
 v6: t35.abi, t78.abi
 v7: t50.abi, t51.abi, t89.abi
 v8: t65.abi
WARNING : Only used if nqpt = 1. If nqpt=0 (which is the default value of nqpt), no reading of shiftq is done, and the default value of shiftq is kept.
It is used only when qptopt >= 0, and must be defined if nshiftq is larger than 1. shiftq (1:3,1:nshiftq) defines nshiftq shifts of the homogeneous grid of q points based on ngqpt or qptrlatt.
See shiftk for more information on the definition, use, and suitable values for these shifts.
slabwsrad¶
Mnemonics: jellium SLAB WignerSeitz RADius
Characteristics: LENGTH
Mentioned in topic(s): topic_Artificial
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [0/151] in abinit tutorials
Fix the bulkmean positive charge density n_{bulk} of a jellium slab (if the latter is employed, e.g. jellslab/=0). Often called r_s (see for example [Lang1970]), slabwsrad is the radius of a sphere which has the same volume as the average volume per particle in a homogeneous electron gas with density n_{bulk}, so:
For example, the bulk aluminum fcc lattice constant is a=4.0495 Angstroms WebElements, each cubic centered cell includes 4 Al atoms and each atom has 3 valence electrons, so the average volume per electron is a^3/12=37.34 Bohr^3 which has to be equal to \frac{4 \pi}{3} r_s^3. Consequently Al has approximately r_s=2.07 Bohr, while for example magnesium has r_s=2.65 Bohr, sodium 3.99 Bohr. By default, given in Bohr atomic units (1 Bohr=0.5291772108 Angstroms).
slabzbeg¶
Mnemonics: jellium SLAB BEGinning edge along the zdirection
Mentioned in topic(s): topic_Artificial
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [0/151] in abinit tutorials
Define the edges of the jellium slab (if used, so if jellslab/=0) along z, namely the slab starts at a point along z which is expressed in Bohr by slabzbeg and it ends at a point expressed in Bohr by slabzend. The zdirection is parallel to the third crystal primitive lattice vector which has to be orthogonal to the other ones, so the length of the cell along z is %rprimd(3,3). In addition slabzbeg and slabzend have to be such that:
0 â‰¤ **slabzbeg** < [[slabzend]] â‰¤ [[rprimd]](3,3)
Together with slabwsrad they define the jellium positive charge density distribution n_{+}(x,y,z) in this way:
so the positive charge density is invariant along the xy plane as well as the electrostatic potential generated by it.
slabzend¶
Mnemonics: jellium SLAB ENDing edge along the zdirection
Mentioned in topic(s): topic_Artificial
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [4/1101] in all abinit tests, [0/151] in abinit tutorials
Define the edges of the jellium slab (if used, so if jellslab/=0) along z, namely the slab starts at a point along z which is expressed in Bohr by slabzbeg and it ends at a point expressed in Bohr by slabzend . The zdirection is parallel to the third crystal primitive lattice vector which has to be orthogonal to the other ones, so the length of the cell along z is %rprimd(3,3). In addition slabzbeg and slabzend have to be such that:
0 â‰¤ [[slabzbeg]] < **slabzend** â‰¤ [[rprimd]](3,3)
Together with slabwsrad they define the jellium positive charge density distribution n_{+}(x,y,z) in this way:
so the positive charge density is invariant along the xy plane as well as the electrostatic potential generated by it.
slk_rankpp¶
Mnemonics: ScaLapacK matrix RANK Per Process
Mentioned in topic(s): topic_parallelism
Variable type: integer
Dimensions: scalar
Default value: [1000]
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 paral: t31.abi
This variable controls how the number of processes to be used in Scalapack diagonalization algorithm: np_slk will be calculated according to this value. This value is the matrix rank that each process will hold for the diagonalization. For a 1000x1000 matrix with default value, scalapack won’t be used (Lapack will be used). For a 2000x2000 matrix with default value, scalapack will be used with 2000/1000=2 MPI processes. For a 2000x2000 matrix with a slk_rank=500, scalapack will be used with 2000/500=4 MPI processes. In case of hybrid MPI+OpenMP, the number of thread is also taken into account.
WARNING None of the available scalapack library are threadsafe (2018). Therefore using both scalapack and OpenMP is highly unpredictable. Furthermore, using multithreaded linear algebra library (MKL ACML…) is more efficient than pure MPI scalapack.
Usually it is better to define this variable and let the code do the rest.
so_psp¶
Mnemonics: SpinOrbit treatment for each PSeudoPotential
Mentioned in topic(s): topic_spinpolarisation
Variable type: integer
Dimensions: (npsp)
Default value: npsp * 1
Only relevant if: nspinor == 2 and %usepaw == 0
Added in version: before_v9
Test list (click to open). Moderately used, [49/1101] in all abinit tests, [0/151] in abinit tutorials
For each type of atom (each pseudopotential), specify the treatment of spinorbit interaction (if nspinor == 2 and Normconserving pseudopotentials i.e. %usepaw == 0) For PAW calculations with SOC, please refer to pawspnorb.
 If 0: no spinorbit interaction, even if nspinor = 2
 If 1: treat spinorbit as specified in the pseudopotential file.
 If 2: treat spinorbit in the HGH form (usual form, although not allowed for all pseudopotentials)
 If 3: treat spinorbit in the HFN form (HemstreetFongNelson) (actually, not implemented).
For typical usage, the default value is OK. If the spinorbit needs to be turned off for one atom, 0 might be relevant. Note however, that the code will stop if nspinor == 2 is used and one of the pseudopotential does not contain the information about the spinorbit interaction (this is the case for some old pseudopotentials). Indeed, for spinorial calculations, turning off the spinorbit interaction is unphysical, and also does not save CPU time. It should only be done for test purposes
Note that if nspinor == 1, the spinorbit cannot be treated anyhow, so the value of so_psp is irrelevant.
Prior to v5.4, the input variable so_typat was used, in place of so_psp . Because the values 0 and 1 have been switched between so_psp and so_typat, it was dangerous to continue to allow the use of so_typat.
spinat¶
Mnemonics: SPIN for AToms
Mentioned in topic(s): topic_spinpolarisation, topic_crystal, topic_MagMom, topic_ConstrainedDFT
Variable type: real
Dimensions: [3, natrd ] if natrd<natom,
[3, natom ] otherwise.
Default value: 0.0
Added in version: before_v9
Test list (click to open). Moderately used, [135/1101] in all abinit tests, [17/151] in abinit tutorials
 atompaw: t04.abi …
 builtin: testin_libxc.abi …
 gpu: t03.abi, t05.abi, t05.abi …
 libxc: t00.abi, t02.abi, t05.abi …
 paral: t28.abi, t56.abi, t56.abi …
 psml: t02.abi, t06.abi, t10.abi …
 tutorial: tbase1_5.abi, tbase2_1.abi, tbase2_2.abi …
 v1: t08.abi, t09.abi, t39.abi …
 v2: t50.abi, t51.abi, t54.abi …
 v3: t12.abi, t15.abi, t19.abi …
 v4: t08.abi, t21.abi, t22.abi …
 v5: t06.abi, t08.abi, t16.abi …
 v6: t07.abi, t17.abi, t45.abi …
 v7: t05.abi, t06.abi, t10.abi …
 v8: t25.abi, t26.abi, t27.abi …
 v9: t02.abi, t03.abi, t41.abi …
 wannier90: t04.abi …
Gives the initial electronic spinmagnetization for each atom, in unit of \hbar/2, as well as, in case of fixed magnetization calculations (see constraint_kind and magconon), the target value of the magnetization.
Note that if nspden = 2, the zcomponent must be given for each atom, in triplets (0 0 zcomponent). For example, the electron of an hydrogen atom can be spin up (0 0 1.0) or spin down (0 0 1.0).
This value is only used to create the first exchange and correlation potential. It is not checked against the initial occupation numbers occ for each spin channel. It is meant to give an easy way to break the spin symmetry, and to allow to find stable local spin fluctuations, for example: antiferromagnetism, or the spontaneous spatial spin separation of elongated H_2 molecule.

If the atom manipulator is used, spinat will be related to the preprocessed set of atoms, generated by the atom manipulator. The user must thus foresee the effect of this atom manipulator (see objarf).

If the atom manipulator is not used, and the symmetries are not specified by the user (nsym = 0), spinat will be used, if present, to determine the antiferromagnetic characteristics of the symmetry operations, see symafm. In case of collinear antiferromagnetism (nsppol = 1, nspinor = 1, nspden = 2), these symmetries are used to symmetrize the density. In case of noncollinear magnetism (nsppol = 1, nspinor = 2, nspden = 4), they are also used to symmetrize the density. In the latter case, this strongly constrains the magnetization (imposing its direction). If the user want to let all degrees of freedom of the magnetization evolve, it is then recommended to put nsym = 1.

If the symmetries are specified, and the irreducible set of atoms is specified, the antiferromagnetic characteristics of the symmetry operations symafm will be used to generate spinat for all the nonirreducible atoms.

In the case of PAW+U calculations using the dmatpawu initial occupation matrix, and if nspden = 4, spinat is also used to determine the direction of the integrated magnetization matrix.
stmbias¶
Mnemonics: Scanning Tunneling Microscopy BIAS voltage
Characteristics: ENERGY
Mentioned in topic(s): topic_STM
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v4: t46.abi
Gives, in Hartree, the bias of the STM tip, with respect to the sample, in order to generate the STM density map. Used with positive iscf, occopt = 7 (metallic, gaussian), nstep = 1, and positive prtstm, this value is used to generate a charge density map from electrons close to the Fermi energy, in a (positive or negative) energy range. Positive stmbias will lead to the inclusion of occupied (valence) states only, while negative stmbias will lead to the inclusion of unoccupied (conduction) states only. Can be specified in Ha (the default), Ry, eV or Kelvin, since stmbias has the ENERGY characteristics (0.001 Ha = 27.2113845 meV = 315.773 Kelvin). With occopt = 7, one has also to specify an independent broadening tsmear.
supercell_latt¶
Mnemonics: SUPERCELL LATTice
Mentioned in topic(s): topic_UnitCell
Variable type: integer
Dimensions: (3)
Default value: [1, 1, 1]
Added in version: before v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v8: t94.abi
Factor that multiplies the three primitive vectors to generate a supercell. The starting nuclei types typat and coordinates xred and xcart are simply replicated in the different images of the original primitive cell. If present, spinat and chrgat are also replicated without change. Note that supercell_latt is not echoed. It is immediately used to change natom, typat, and all input variables that depend on natom and typat.
This input variable cannot be used jointly with nondefault values of natrd and nobj. Also, nsym is set to 1. Also note that vel, if defined, must have the dimensions of the number of atoms in the supercell, not the initial natom, because it is not replicated. supercell_latt is typically used to initialize molecular dynamics runs.
symafm¶
Mnemonics: SYMmetries, AntiFerroMagnetic characteristics
Mentioned in topic(s): topic_spinpolarisation
Variable type: integer
Dimensions: (nsym)
Default value: nsym * 1
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
In case the material is magnetic, nspden=2 or 4, additional symmetry operations might appear, that change the sign of the magnetization (spinflip). They have been introduced by Shubnikov in 1951 [Bradley1972]. They can be used by ABINIT to decrease the CPU time, either by decreasing the number of kpoints or by suppressing the explicit treatment of one spin channel, or by decreasing the number of perturbations in DFPT.
symafm should be set to +1 for all the usual symmetry operations, that do not change the sign of the magnetization, while it should be set to 1 for the magnetizationchanging operations (spinflip). If the symmetry operations are not specified by the user in the input file, that is, if nsym = 0, then ABINIT will use the values of spinat to determine the content of symafm .
The symmetries that can act on the magnetization can yield decreased CPU time (and usually also memory decrease) in the following cases:
 antiferromagnetism (nsppol = 1, nspinor = 1, nspden = 2)
 noncollinear magnetism (nsppol = 1, nspinor = 2, nspden = 4)
Also in the case nsppol = 2, nspinor = 1, nspden = 2 they might simply yield better accuracy (or faster convergence), but there is no automatic gain of CPU time or memory, although it is not as clear cut as in the above cases.
IMPORTANT : The meaning of symafm is different in the nspden = 2 case (collinear magnetism), and in the nspden = 4 case (noncollinear magnetism, with explicit treatment of magnetization as a vector). Indeed in the first case, it is supposed that the magnetization vector is not affected by the real space symmetry operations (socalled black and white symmetry groups). By contrast, in the second case, the real space symmetry operations act on the magnetization vector. The rationale for such different treatment comes from the fact that the treatment of spinorbit coupling is incompatible with collinear magnetism nspden=2, so there is no need to worry about it in this case. On the contrary, many calculations with nspden=2 will include spinorbit coupling. The symmetry operations should thus act coherently on the spinorbit coupling, which implies that the real space operations should act also on the magnetization vector in the nspden=4 case. So, with nspden=4, even with symafm =1, symmetry operations might change the magnetization vector, e.g. possibly reverse it from one atom to another atom. Still, when real space operations also act on the magnetization vector, nothing prevents to have ADDITIONAL “spinflip” operations, which is indeed then the meaning of symafm =1 in the nspden=4 case.
Let’s illustrate this with an example. Take an H_2 system, with the two H atoms quite distant from each other. The electron on the first H atom might be 1s spin up, and the electron on the second atom might be 1s spin down. With nspden=2, the inversion symmetry centered in the middle of the segment joining the two atoms will NOT act on the spin, so that the actual symmetry operation that leaves the system invariant is an inversion (symrel= 1 0 0 0 1 0 0 0 1) accompanied by a spinflip with symafm =1. By contrast, with nspden=4, the inversion symmetry centered in the middle of the segment joining the two atoms will reverse the spin direction as well, so that the proper symmetry operation is symrel= 1 0 0 0 1 0 0 0 1 but no additional spinflip is needed to obtain a symmetry operation that leaves the system invariant, so that symafm =1.
Although this might seem confusing, ABINIT is able to recognise the correct symmetry operations from the available atomic coordinates, from spinat, and from nspden, so that the user should hardly be affected by such different conventions. However, the use of %ptgroupma and spgroupma to define the antiferromagnetic operations of symmetry should be done carefully.
Related variables: symrel, tnons, %ptgroupma, spgroupma.
timopt¶
Mnemonics: TIMing OPTion
Characteristics: NO_MULTI
Mentioned in topic(s): topic_Control
Variable type: integer
Dimensions: scalar
Default value: 1 if SEQUENTIAL,
0 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [92/1101] in all abinit tests, [18/151] in abinit tutorials
 fast: t02.abi, t03.abi, t21.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 tutoparal: tdfpt_02.abi, tdfpt_04.abi, tgspw_01.abi …
 v1: t21.abi …
 v2: t59.abi …
 v4: t97.abi …
 v5: t63.abi, t66.abi, t75.abi …
 v6: t03.abi, t30.abi, t34.abi …
 v67mbpt: t06.abi, t14.abi …
 v8: t21.abi, t22.abi …
 v9: t57.abi …
This input variable allows one to modulate the use of the timing routines.
 If 0 → as soon as possible, suppresses all calls to timing routines
 If 1 → usual timing behaviour, with short analysis, appropriate for sequential execution
 If 2 → close to timopt = 1, except that the analysis routine does not time the timer, appropriate for parallel execution.
 If 3 → close to timopt = 1, except that the different parts of the lobpcg routine are timed in detail.
 If 4 → close to timopt = 1, except that the different parts of the lobpcg routine are timed in detail. A different splitting of lobpcg than for timopt = 3 is provided.
 If 10 → relevant only when the wavelet basis set is used (i.e. usewvl = 1).
When activated, a file named
wvl_timings.yaml
, in YAML format, is created. It contains a time analysis of the BigDFT wavelet routines. See the tutorial on parallelism using wavelets  If 1 → a full analysis of timings is delivered
 If 2 → a full analysis of timings is delivered, except timing the timer
 If 3 → a full analysis of timings is delivered, including the detailed timing of the different parts of the lobpcg routine. (this takes time, and is discouraged for too small runs  the timing would take more time than the run !). The timer is timed.
 If 4 → a full analysis of timings is delivered, including the detailed timing of the different parts of the lobpcg routine. A different splitting of lobpcg than for timopt = 3 is provided (this takes time, and is discouraged for too small runs  the timing would take more time than the run !). The timer is timed. The sum of the independent parts is closer to 100% than for timopt = 3.
tl_nprccg¶
Mnemonics: TaiL maximum Number of PReConditioner Conjugate Gradient iterations
Mentioned in topic(s): topic_Wavelets
Variable type: integer
Dimensions: scalar
Default value: 30
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 bigdft: t09.abi
This variable is similar to wvl_nprccg but for the preconditioner iterations during the tail corrections (see tl_radius).
tl_radius¶
Mnemonics: TaiL expansion RADIUS
Characteristics: LENGTH
Mentioned in topic(s): topic_Wavelets
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 bigdft: t09.abi
In the wavelet computation case, the linkage between the grid and the free boundary conditions can be smoothed using an exponential decay. This means a correction on the energy at the end on each wavefunction optimisation run. If this parameter is set to zero, no tail computation is done. On the contrary, put it to a positive value makes the tail correction available. The value correspond to a length in atomic units being the spacial expansion with the exponential decay around the grid.
tphysel¶
Mnemonics: Temperature (PHYSical) of the ELectrons
Characteristics: ENERGY
Mentioned in topic(s): topic_BandOcc
Variable type: real
Dimensions: scalar
Default value: 0.0
Added in version: before_v9
Test list (click to open). Rarely used, [5/1101] in all abinit tests, [0/151] in abinit tutorials
Gives, in Hartree, the physical temperature of the system, in case occopt = 4, 5, 6, or 7.
Can be specified in Ha (the default), Ry, eV or Kelvin, since tphysel has the ENERGY characteristics (0.001 Ha = 27.2113845 meV = 315.773 Kelvin). One has to specify an independent broadening tsmear. The combination of the two parameters tphysel and tsmear is described in [Verstraete2002]. Note that the signification of the entropy is modified with respect to the usual entropy. The choice has been made to use tsmear as a prefactor of the entropy, to define the entropy contribution to the free energy.
tsmear¶
Mnemonics: Temperature of SMEARing
Characteristics: ENERGY
Mentioned in topic(s): topic_BandOcc, topic_STM
Variable type: real
Dimensions: scalar
Default value: 0.01
Added in version: before_v9
Test list (click to open). Moderately used, [292/1101] in all abinit tests, [42/151] in abinit tutorials
 atompaw: t04.abi …
 bigdft: t14.abi …
 builtin: testin_etsf_io.abi …
 etsf_io: t00.abi, t09.abi …
 gpu: t02.abi, t03.abi, t05.abi …
 libxc: t03.abi, t70.abi, t71.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t03.abi, t03.abi, t03.abi …
 psml: t02.abi, t03.abi, t04.abi …
 seq: tsv2_81.abi …
 tutoparal: tdmft_1.abi, tdmft_2.abi, tgspw_01.abi …
 tutorespfn: telast_6.abi, teph_1.abi, teph4isotc_2.abi …
 tutorial: tbase4_1.abi, tbase4_2.abi, tbase4_3.abi …
 v1: t01.abi, t02.abi, t03.abi …
 v2: t11.abi, t26.abi, t36.abi …
 v3: t01.abi, t18.abi, t22.abi …
 v4: t01.abi, t04.abi, t05.abi …
 v5: t01.abi, t03.abi, t06.abi …
 v6: t07.abi, t08.abi, t09.abi …
 v7: t04.abi, t05.abi, t06.abi …
 v8: t01.abi, t01_triqs2_0.abi, t04.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t04.abi …
Gives the broadening of occupation numbers occ, in the metallic cases (occopt = 3, 4, 5, 6 and 7). Can be specified in Ha (the default), eV, Ry, or Kelvin, since tsmear has the ENERGY characteristics (0.001 Ha = 27.2113845 meV = 315.773 Kelvin). Default is 0.01 Ha. This should be OK using gaussian like smearings (occopt = 4,5,6,7) for a freeelectron metal like Al. For dband metals, you may need to use less. Always check the convergence of the calculation with respect to this parameter, and simultaneously, with respect to the sampling of kpoints (see nkpt) If occopt = 3, tsmear is the physical temperature, as the broadening is based on FermiDirac statistics. However, if occopt = 4, 5, 6, or 7, the broadening is not based on FermiDirac statistics, and tsmear is only a convergence parameter. It is still possible to define a physical temperature, thanks to the input variable tphysel (See also [Verstraete2002]).
useextfpmd¶
Mnemonics: USE EXTended FPMD model
Mentioned in topic(s): topic_ExtFPMD
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: 9.5.2
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v9: t92.abi
Enables the calculation of contributions to the energy, entropy, stresses, number of electrons and chemical potential using the extended first principle molecular dynamics model for high temperature simulations.

useextfpmd = 1 (Recommanded), the energy shift factor will be evaluated by making an integration of the trial potential over the real space and the contributions will be computed with integrals over the band number.

useextfpmd = 2, the energy shift factor will be evaluated by making the average between the eigenvalues and the Fermi gas energy over the last extfpmd_nbcut bands, and the contributions will be computed with integrals over the band number.

useextfpmd = 3, the energy shift factor will be evaluated by making the average between the eigenvalues and the kinetic energies over the last extfpmd_nbcut bands, and the contributions will be computed using the density of states of the Fermi gas.
usekden¶
Mnemonics: USE Kinetic energy DENsity
Mentioned in topic(s): topic_xc
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [13/1101] in all abinit tests, [0/151] in abinit tutorials
If usekden = 1 the kinetic energy density will be computed during the selfconsistent loop, in a way similar to the computation of the density. This is needed if a metaGGA is to be used as XC functional. By default ( usekden = 0), the kinetic energy density is not computed during the self consistent loop.
vacuum¶
Mnemonics: VACUUM identification
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: (3)
Default value: None
Added in version: before_v9
Test list (click to open). Moderately used, [22/1101] in all abinit tests, [4/151] in abinit tutorials
 paral: t03.abi, t03.abi, t03.abi …
 tutorial: tbase4_5.abi, tbase4_6.abi, tbase4_7.abi …
 v1: t41.abi, t55.abi, t56.abi …
 v2: t58.abi …
 v4: t46.abi …
 v5: t57.abi, t62.abi, t73.abi …
 v7: t04.abi …
Establishes the presence (if vacuum = 1) or absence (if vacuum = 0) of a vacuum layer, along the three possible directions normal to the primitive axes.
This information might be used to generate kpoint grids, if kptopt = 0 and neither ngkpt nor kptrlatt are defined (see explanations with the input variable prtkpt). It will allow to select a zero, one, two or threedimensional grid of k points. The coordinate of the k points along vacuum directions is automatically set to zero.
If vacuum is not defined, the input variable vacwidth will be used to determine automatically whether the distance between atoms is sufficient to have the presence or absence of vacuum.
vacwidth¶
Mnemonics: VACuum WIDTH
Characteristics: INPUT_ONLY, LENGTH
Mentioned in topic(s): topic_kpoints
Variable type: real
Dimensions: scalar
Default value: 10.0
Added in version: before_v9
Test list (click to open). Rarely used, [1/1101] in all abinit tests, [0/151] in abinit tutorials
 v2: t58.abi
Give a minimum “projected” distance between atoms to be found in order to declare that there is some vacuum present for each of the three directions. By default, given in Bohr atomic units (1 Bohr=0.5291772108 Angstroms), although Angstrom can be specified, if preferred, since vacwidth has the LENGTH characteristics. The precise requirement is that a slab of width vacwidth , delimited by two planes of constant reduced coordinates in the investigated direction, must be empty of atoms.
vcutgeo¶
Mnemonics: V (potential) CUToff GEOmetry
Mentioned in topic(s): topic_Coulomb, topic_GWls, topic_Susceptibility, topic_SelfEnergy
Variable type: real
Dimensions: (3)
Default value: 3 * 0.0
Only relevant if: icutcoul in [1,2]
Added in version: before_v9
Test list (click to open). Rarely used, [2/1101] in all abinit tests, [0/151] in abinit tutorials
vcutgeo is used in conjunction with icutcoul, fock_icutcoul and/or gw_icutcoul to specify the geometry used to truncate the Coulomb interaction, as well as the particular approach to be used. It has a meaning either for a periodic onedimensional system, typically a nanowire, nanotube or polymer surrounded by vacuum separating the system from images in neighbouring cells (icutcoul = 1) or in the case of periodic twodimensional system, typically a slab with vacuum separating it from images in neighbouring cells (icutcoul = 2). For each geometry, two different definitions of the cutoff region are available (see [IsmailBeigi2006] and [Rozzi2006] for a complete description of the methods)
In the method of IsmailBeigi [IsmailBeigi2006] for onedimensional systems, the cutoff region is given by the WignerSeitz cell centered on the axis of the cylinder. The cutoff region is thus automatically defined by the unit cell and there is no need to specify the value of rcut.
To define a cylinder along the zaxis use the following lines:
icutcoul 1
vcutgeo 0 0 1
Please note that the method of IsmailBeigi is implemented only in the case of an orthorhombic Bravais lattice. For hexagonal lattices, one has to use the method of Rozzi [Rozzi2006]. In this case, the interaction is truncated in a finite cylinder. Contrarily to the first approach, the user has to specify both the radius of the cylinder with rcut as well as the length of the cylinder along the periodic dimension, that should always be smaller than the extension of the Born von Karman box. The length of the cylinder is given in terms of a multiple of the primitive vector along the periodic direction. Another option provided by Rozzi [Rozzi2006] is the infinite length cylinder. In order to activate it in ABINIT, one needs to use a very large negative vcutgeo value on the third direction (i.e. vcutgeo(3) <= 999).
For example, in order to define a finite cylinder along z of radius 2.5 Bohr and length 3*R3,
icutcoul 1
vcutgeo 0 0 3.0 # note the minus sign
rcut 2.5
For twodimensional systems (icutcoul = 2), vcutgeo is used to define the two periodic directions. Also in this case two different techniques are available. In the method of IsmailBeigi, the (positive) nonzero components of vcutgeo define the periodic directions of the infinite surface. The interaction is truncated within a slab of width L where L is the length of the primitive vector of the lattice along the nonperiodic dimension. For example:
icutcoul 2
vcutgeo 1 1 0
In Rozzi’s method, it is also possible to define a finite range for the Coulomb interaction in the periodic directions by employing negative values. For example:
icutcoul 2
vcutgeo 3 2 0
Note that not all kpoint grids are allowed in these 1D and 2D cases: the k point vector component(s) along nonperiodic direction(s) must vanish. So, if the 2D IsmailBeigi technique is used, the zcomponent of the k points must vanish.
wfk_task¶
Mnemonics: WFK TASK
Mentioned in topic(s): topic_ElecBandStructure
Variable type: string
Dimensions: scalar
Default value: 0
Only relevant if: optdriver == 8
Added in version: 9.0.0
Test list (click to open). Rarely used, [6/1101] in all abinit tests, [1/151] in abinit tutorials
This variable defines the quantity to compute starting from a previously generated WFK file. Possible values are:

“wfk_fullbz” → Read WFK file and produce new WFK file with kpoints in the full BZ. Wavefunctions with istwfk > 2 are automatically converted into the full Gsphere representation. This option can be used to interface Abinit with external tools requiring kpoints in the full BZ.

“wfk_einterp” → Read energies from WFK file and interpolate band structure using the parameters specified by einterp.

“wfk_ddk” → Compute DDK matrix elements for all bands and kpoints in the WFK file. The contribution due to the nonlocal part of the pseudopotential can be ignored by setting inclvkb = 0 (not recommended unless you know what you are doing).

“wfk_kpts_erange” → Read WFK file, use starfunction and einterp parameters to interpolate electron energies onto fine kmesh defined by sigma_ngkpt and sigma_shiftk. Find kpoints inside (electron/hole) pockets according to the values specified in sigma_erange. Write KERANGE.nc file with the tables required by the code to automate NSCF band structure calculations inside the pocket(s) and electron lifetime computation in the EPH code when eph_task = 4.
wfmix¶
Mnemonics: WaveFunctions MIXing factor
Mentioned in topic(s): topic_Hybrids
Variable type: real
Dimensions: scalar
Default value: 1.0
Only relevant if: %usefock > 0 and nnsclohf >0 and fockoptmix/100 > 0
Added in version: before_v9
Test list (click to open). Rarely used, [3/1101] in all abinit tests, [0/151] in abinit tutorials
When the wavefunctions are determined using a SCF double loop (hybrid functionals), wfmix provides the mixing factor to obtain the new input wavefunctions by the combination of the earlier input wavefunctions and corresponding (DFTpreconditioned) output wavefunctions at the level of the outer loop, according to the algorithm specified by fockoptmix/100. If wfmix is 1.0, the output wavefunctions only will determine the new input wavefunctions. This might possibly lead to instabilities. If wfmix is smaller than 1.0, the whole iteration procedure is da