Keywords

Summary of all keywords

AIMCriticalPoints
Type

Block

Description

Compute the critical points of the density (Atoms In Molecules). The algorithm starts from a regular mesh of points, and from each of these it walks towards its corresponding critical point.

Enabled
Type

Bool

Default value

No

GUI name

: Critical points and bond paths

Description

Compute the critical points of the density (Atoms In Molecules). The algorithm starts from a regular mesh of points, and from each of these it walks towards its corresponding critical point.

EqvPointsTol
Type

Float

Default value

0.27

Unit

Bohr

Description

If the distance between two critical points is smaller than this value, the two critical points are considered to be the same point.

GridPadding
Type

Float

Default value

0.7

Unit

Bohr

Description

How much extra space is added to the starting guess domain in the search for the critical points

GridSpacing
Type

Float

Default value

0.5

Unit

Bohr

Description

The distance between the initial trial points.

Allow
Type

String

Recurring

True

Description

Debugging feature to let the program continue even when intermediate results seem to be wrong or very inaccurate

ATensor
Type

Block

Description

Hyperfine A-tensor.

Enabled
Type

Bool

Default value

No

GUI name

:A-tensor

Description

Compute the hyperfine A-tensor. Note: Unrestricted calculation is required.

AtomType
Type

Block

Recurring

True

Description

Explicit basis set definition for given atom type.

AutomaticGaussians
Type

Non-standard block

Description

Definition of the automatic gaussians

BasisFunctions
Type

Non-standard block

Description

Definition of the extra Slater-type orbitals

ContractedGaussians
Type

Non-standard block

Description

Definition of the contracted gaussians

Dirac
Type

Non-standard block

Description

Specification of the numerical (‘Herman-Skillman’) free atom, which defines the initial guess for the SCF density, and which also (optionally) supplies Numerical Atomic Orbitals (NOs) as basis functions

FitFunctions
Type

Non-standard block

Description

Slater-type fit functions. Obsolete feature.

BandStructure
Type

Block

Description

Options for the calculation of the band structure.

Automatic
Type

Bool

Default value

Yes

GUI name

Automatic generate path

Description

If True, BAND will automatically generate the standard path through the Brillouin zone. If False BAND will use the user-defined path in BZPath.

DeltaK
Type

Float

Default value

0.1

Unit

1/Bohr

GUI name

Interpolation delta-K

Description

Step (in reciprocal space) for band structure interpolation. Using a smaller number (e.g. 0.03) will result in smoother band curves at the cost of an increased computation time.

Enabled
Type

Bool

Default value

No

GUI name

Calculate band structure

Description

If True, Band will calculate the band structure and save it to file for visualization.

EnergyAboveFermi
Type

Float

Default value

0.75

Unit

Hartree

GUI name

Energy above Fermi level

Description

Bands with minimum energy larger then FermiEnergy + EnergyAboveFermi are not saved to file. Increasing the value of EnergyAboveFermi will result in more unoccupied bands to be saved to file for visualization.

EnergyBelowFermi
Type

Float

Default value

10.0

Unit

Hartree

GUI name

Energy below Fermi level

Description

Bands with maximum energy smaller then FermiEnergy - EnergyBelowFermi are not saved to file. Increasing the value of EnergyBelowFermi will result in more occupied core bands to be saved to file for visualization. Note: EnergyBelowFermi should be a positive number!

FatBands
Type

Bool

Default value

Yes

GUI name

Calculate fatbands

Description

If True, BAND will compute the fat bands (only if BandStructure%Enabled is True). The Fat Bands are the periodic equivalent of the Mulliken population analysis.

UseSymmetry
Type

Bool

Default value

Yes

GUI name

Use symmetry

Description

If True, only the irreducible wedge of the Wigner-Seitz cell is sampled. If False, the whole (inversion-unique) Wigner-Seitz cell is sampled. Note: The Symmetry key does not influence the symmetry of the band structure sampling.

Basis
Type

Block

Description

Definition of the basis set

Core
Type

Multiple Choice

Default value

Large

Options

[None, Small, Medium, Large]

GUI name

Frozen core

Description

Select the size of the frozen core you want to use. Small, Medium, and Large will be interpreted within the basis sets available (of the selected quality), and might refer to the same core in some cases.

Folder
Type

String

Description

Path to a folder containing the basis set files. This can be used for special use-defined basis sets. Cannot be used in combination with ‘Type’

PerAtomType
Type

Block

Recurring

True

Description

Defines the basis set for all atoms of a particular type.

Core
Type

Multiple Choice

Options

[None, Small, Medium, Large]

Description

Size of the frozen core.

File
Type

String

Description

The path to the basis set file. The path can be absolute or relative to $AMSRESOURCES/Band. Specifying the path to the basis file explicitly overrides the automatic basis file selection via the Type and Core subkeys.

Symbol
Type

String

Description

The symbol for which to define the basis set.

Type
Type

Multiple Choice

Options

[SZ, DZ, DZP, TZP, TZ2P, QZ4P]

Description

The basis sets to be used.

PerRegion
Type

Block

Recurring

True

Description

Defines the basis set for all atoms in a region. If specified, this overwrites the values set with the Basis%Type and Basis%PerAtomType keywords for atoms in that region. Note that if this keyword is used multiple times, the chosen regions may not overlap.

Core
Type

Multiple Choice

Default value

Large

Options

[None, Small, Medium, Large]

Description

Size of the frozen core.

Region
Type

String

Description

The identifier of the region for which to define the basis set. Note that this may also be a region expression, e.g. ‘myregion+myotherregion’ (the union of two regions).

Type
Type

Multiple Choice

Default value

DZ

Options

[SZ, DZ, DZP, TZP, TZ2P, QZ4P]

Description

The basis sets to be used.

Type
Type

Multiple Choice

Default value

DZ

Options

[SZ, DZ, DZP, TZP, TZ2P, QZ4P, STO/TZ2P, STO/QZ4P, CORR/QZ6P, CORR/TZ3P, GTO/CC-PV5Z, GTO/CC-PV6Z, GTO/CC-PVQZ, GTO/CC-PVTZ, GTO/CC-PVDZ]

GUI name

Basis set

Description

Select the basis set to use. SZ : Single Z DZ : Double Z DZP : Double Z, 1 polarization function TZP : Triple Z, 1 polarization function TZ2P : Triple Z, 2 polarization functions QZ4P : Quadruple Z, 4 polarization function The basis set chosen will apply to all atoms in your structure. If a matching basis is not found a better type might be used.

BeckeGrid
Type

Block

Description

Options for the numerical integration grid, which is a refined version of the fuzzy cells integration scheme developed by Becke.

Quality
Type

Multiple Choice

Default value

Auto

Options

[Auto, Basic, Normal, Good, VeryGood, Excellent]

Description

Quality of the integration grid. For a description of the various qualities and the associated numerical accuracy see reference. If ‘Auto’, the quality defined in the ‘NumericalQuality’ will be used.

QualityPerRegion
Type

Block

Recurring

True

Description

Sets the grid quality for all atoms in a region. If specified, this overwrites the globally set quality.

Quality
Type

Multiple Choice

Options

[Basic, Normal, Good, VeryGood, Excellent]

Description

The region’s integration grid quality.

Region
Type

String

Description

The identifier of the region for which to set the quality.

RadialGridBoost
Type

Float

Default value

1.0

Description

The number of radial grid points will be boosted by this factor. Some XC functionals require very accurate radial integration grids, so BAND will automatically boost the radial grid by a factor 3 for the following numerically sensitive functionals: LibXC M05, LibXC M05-2X, LibXC M06-2X, LibXC M06-HF, LibXC M06-L, LibXC M08-HX, LibXC M08-SO, LibXC M11-L, LibXC MS0, LibXC MS1, LibXC MS2, LibXC MS2H, LibXC MVS, LibXC MVSH, LibXC N12, LibXC N12-SX, LibXC SOGGA11, LibXC SOGGA11-X, LibXC TH1, LibXC TH2, LibXC WB97, LibXC WB97X, MetaGGA M06L, MetaHybrid M06-2X, MetaHybrid M06-HF, MetaGGA MVS.

BerryPhase
Type

Bool

Default value

No

Description

Boolean that determines whether the dipole as determined through the Berry phase approach should be calculated.

BField
Type

Block

Description

The effect of a magnetic filed can be approximated by the following potential: mu * sigma_i * B, where mu is the Bohr magneton, sigma_i are the Pauli matrices and B is the magnetic field

Bx
Type

Float

Default value

0.0

Unit

Tesla

Description

Value of the x component of the BField

By
Type

Float

Default value

0.0

Unit

Tesla

Description

Value of the y component of the BField

Bz
Type

Float

Default value

0.0

Unit

Tesla

Description

Value of the z component of the BField

Dipole
Type

Bool

Default value

No

GUI name

Bfield is: Atomic dipole

Description

Use an atomic dipole as magnetic field instead of a uniform magnetic field.

DipoleAtom
Type

Integer

Default value

1

GUI name

on atom number

Description

Atom on which the magnetic dipole should be centered (if using the dipole option)

Method
Type

Multiple Choice

Default value

NR_SDOTB

Options

[NR_SDOTB, NR_LDOTB, NR_SDOTB_LDOTB]

Description

There are two terms coupling to an external magnetic field. One is the intrinsic spin of the electron, called S-dot-B, the other one is the orbital momentum call L-dot-B. The L.B is implemented non-relativistically, using GIAOs in the case of a homogeneous magnetic field (not for the dipole case).

Unit
Type

Multiple Choice

Default value

tesla

Options

[tesla, a.u.]

Description

Unit of magnetic filed. The a.u. is the SI version of a.u.

BZPath
Type

Block

Description

Definition of the user-defined path in the Brillouin zone for band structure plotting.

path
Type

Non-standard block

Recurring

True

Description

Definition of the k-points in a path. The vertices of your path should be defined in fractional coordinates (wrt the reciprocal lattice vectors)

Comment
Type

Non-standard block

Description

The content of this block will be copied to the output header as a comment to the calculation.

Convergence
Type

Block

Description

Options and parameters related to the convergence behavior of the SCF procedure.

Criterion
Type

Float

Description

Criterion for termination of the SCF procedure. The default depends on the NumericalQuality and on the number of atoms in the system. Can be used for EngineAutomations

CriterionFactor
Type

Float

Default value

1.0

Description

Multiply Criterion (which depends on system and quality) with this factor. Can be used for EngineAutomations

Degenerate
Type

String

Default value

default

Description

Smooths (slightly) occupation numbers around the Fermi level, so as to insure that nearly-degenerate states get (nearly-) identical occupations. Be aware: In case of problematic SCF convergence the program will turn this key on automatically, unless the key ‘Nodegenerate’ is set in input. The smoothing depends on the argument to this key, which can be considered a ‘degeneration width’. When the argument reads default, the program will use the value 1e-4 a.u. for the energy width.

ElectronicTemperature
Type

Float

Default value

0.0

Unit

Hartree

Description

(KT) Specify this key for a gradient independent electronic temperature

InitialDensity
Type

Multiple Choice

Default value

rho

Options

[rho, psi, frompot]

Description

The SCF is started with a guess of the density. There are the following choices RHO: the sum of atomic density. PSI: construct an initial eigensystem by occupying the atomic orbitals. The guessed eigensystem is orthonormalized, and from this the density is calculated/

LessDegenerate
Type

Bool

Default value

No

Description

If smoothing of occupations over nearly degenerate orbitals is applied (see Degenerate key), then, if this key is set in the input file, the program will limit the smoothing energy range to 1e-4 a.u. as soon as the SCF has converged ‘halfway’, i.e. when the SCF error has decreased to the square root of its convergence criterion.

ModestCriterion
Type

Float

Default value

-1.0

Description

If this is specified band will consider the SCF converged if the error is below this criterion (after using the maximum number of iterations).

NoDegenerate
Type

Bool

Default value

No

Description

This key prevents any internal automatic setting of the key DEGENERATE.

NumBoltz
Type

Integer

Default value

10

Description

The electronic temperature is done with a Riemann Stieltjes numerical integration, between zero and one occupation. This defines the number of points to be used.

SpinFlip
Type

Integer List

GUI name

Flip spin for atoms

Description

List here the atoms for which you want the initial spin polarization to be flipped. This way you can distinguish between ferromagnetic and anti ferromagnetic states. Currently, it is not allowed to give symmetry equivalent atoms a different spin orientation. To achieve that you have to break the symmetry.

SpinFlipEnabled
Type

Bool

Default value

Yes

Description

If set to False, the keys SpinFlip and SpinFlipRegion are ignored. Only useful/convenient when trying to compare in a script the effect of spin flip.

SpinFlipRegion
Type

String

Recurring

True

GUI name

Flip spin for region

Description

Specify here the region for which you want the initial spin polarization to be flipped. This way you can distinguish between ferromagnetic and anti ferromagnetic states. Currently, it is not allowed to give symmetry equivalent atoms a different spin orientation. To achieve that you have to break the symmetry.

StartWithMaxSpin
Type

Bool

Default value

Yes

Description

To break the initial perfect symmetry of up and down densities there are two strategies. One is to occupy the numerical orbitals in a maximum spin configuration. The alternative is to add a constant to the potential. See also Vsplit key.

StartWithMaxSpinForSO
Type

Bool

Default value

No

Description

Same as the StartWithMaxSpin option. In case of spin-orbit band always used to split the potential. Now will use maxspin in case of SpinFlip. With this option it will always do that.

CPVector
Type

Integer

Default value

128

GUI name

Vectorlength (blocksize)

Description

The code is vectorized and this key can be used to set the vector length

DensityPlot
Type

Non-standard block

Description

Plots of the density. Goes together with the Restart%DensityPlot and Grid keys.

Dependency
Type

Block

Description

Criteria for linear dependency of the basis and fit set

AllowBasisDependency
Type

Bool

Default value

Yes

Description

Project out the dependent part of the basis set (associated with small eigenvalues of the overlap matrix).

Basis
Type

Float

Default value

1e-08

GUI name

Dependency criterion

Description

Criteria for linear dependency of the basis: smallest eigenvalue of the overlap matrix of normalized Bloch functions.

Core
Type

Float

Default value

0.8

Description

The program verifies that the frozen core approximation is reasonable, by checking the smallest eigen value of the overlap matrix of the core (Bloch) orbitals (which should ideally be one) is bigger than this criterion.

CoreValence
Type

Float

Default value

1e-05

Description

Criterion for dependency of the core functions on the valence basis. The maximum overlap between any two normalized functions in the two respective function spaces should not exceed 1.0-corevalence

Fit
Type

Float

Default value

5e-06

Description

Criterion for dependency of the total set of fit functions. The value monitored is the smallest eigenvalue of the overlap matrix of normalized Bloch sums of symmetrized fit functions.

DIIS
Type

Block

Description

Parameters for the DIIS procedure to obtain the SCF solution

Adaptable
Type

Bool

Default value

Yes

Description

Change automatically the value of dimix during the SCF.

CHuge
Type

Float

Default value

20.0

GUI name

No DIIS (but damping) when coefs >

Description

When the largest coefficient in the DIIS expansion exceeds this value, damping is applied

CLarge
Type

Float

Default value

20.0

GUI name

Reduce DIIS space when coefs >

Description

When the largest DIIS coefficient exceeds this value, the oldest DIIS vector is removed and the procedure re-applied

Condition
Type

Float

Default value

1000000.0

Description

The condition number of the DIIS matrix, the largest eigenvalue divided by the smallest, must not exceed this value. If this value is exceeded, this vector will be removed.

DiMix
Type

Float

Default value

0.2

GUI name

Bias DIIS towards latest vector with

Description

Mixing parameter for the DIIS procedure

DiMixMax
Type

Float

Default value

-1.0

Description

For adaptive diis: A negative value means automatic, see DiMixatnvctrx. If positive it is an absolute upper bound for (adaptive) dimix

DiMixMin
Type

Float

Default value

0.01

Description

An absolute lower bound for adaptive dimix.

NCycleDamp
Type

Integer

Default value

1

GUI name

Do not start DIIS before cycle

Description

Number of initial iterations where damping is applied, before any DIIS is considered

NVctrx
Type

Integer

Default value

20

GUI name

Size of DIIS space

Description

Maximum number of DIIS expansion vectors

Variant
Type

Multiple Choice

Default value

DIIS

Options

[DIIS, LISTi, LISTb, LISTd]

Description

Which variant to use. In case of problematic SCF convergence, first try MultiSecant, and if that does not work the LISTi is the advised method. Note: LIST is computationally more expensive per SCF iteration than DIIS.

DOS
Type

Block

Description

Density-Of-States (DOS) options

CalcDOS
Type

Bool

Default value

Yes

GUI name

Calculate DOS

Description

Whether or not to calculate the density of states.

CalcPDOS
Type

Bool

Default value

No

GUI name

Calculate PDOS

Description

Whether or not to calculate the partial DOS (projections on basis functions). This can be significantly more expensive than calculating the total DOS

CalcPopulationAnalysis
Type

Bool

Default value

Yes

GUI name

Calculate Mulliken charges

Description

Whether or not to calculate the population analysis. Population analysis can become very expensive when there are many symmetry operators, such as in a super cell.

DeltaE
Type

Float

Default value

0.005

Unit

Hartree

GUI name

Delta E

Description

Energy step for the DOS grid. Using a smaller value (e.g. half the default value) will result in a finer sampling of the DOS.

Energies
Type

Integer

Description

Number of equidistant energy-values for the DOS grid. This keyword is superseded by the ‘DeltaE’ keyword.

File
Type

String

Description

Write the DOS (plain text format) to the specified file instead of writing it to the standard output.

IntegrateDeltaE
Type

Bool

Default value

Yes

Description

This subkey handles which algorithm is used to calculate the data-points in the plotted DOS. If true, the data-points represent an integral over the states in an energy interval. Here, the energy interval depends on the number of Energies and the user-defined upper and lower energy for the calculation of the DOS. The result has as unit [number of states / (energy interval * unit cell)]. If false, the data-points do represent the number of states for a specific energy and the resulting plot is equal to the DOS per unit cell (unit: [1/energy]). Since the resulting plot can be a wild function and one might miss features of the DOS due to the step length between the energies, the default is set to the integration algorithm.

Max
Type

Float

Unit

Hartree

Description

User defined upper bound energy (with respect to the Fermi energy)

Min
Type

Float

Unit

Hartree

Description

User defined lower bound energy (with respect to the Fermi energy)

StoreCoopPerBasPair
Type

Bool

Default value

No

GUI name

Calculate COOP

Description

Calculate the COOP (crystal orbital overlap population).

DosBas
Type

Non-standard block

Description

Used to specify the fragment basis for the DOS.

DumpBasisOnly
Type

Bool

Default value

No

Description

Dump basis and fit set files use for each atom.

EffectiveMass
Type

Block

Description

In a semi-conductor, the mobility of electrons and holes is related to the curvature of the bands at the top of the valence band and the bottom of the conduction band. With the effective mass option, this curvature is obtained by numerical differentiation. The estimation is done with the specified step size, and twice the specified step size, and both results are printed to give a hint on the accuracy. The easiest way to use this key is to enabled it without specifying any extra options.

Enabled
Type

Bool

Default value

No

GUI name

Effective mass

Description

Compute the EffectiveMass.

KPointCoord
Type

Float List

Unit

1/Bohr

Recurring

True

GUI name

At K-point

Description

Coordinate of the k-points for which you would like to compute the effective mass.

NumAbove
Type

Integer

Default value

1

GUI name

Include N bands above

Description

Number of bands to take into account above the Fermi level.

NumBelow
Type

Integer

Default value

1

GUI name

Include N bands below

Description

Number of bands to take into account below the Fermi level.

StepSize
Type

Float

Default value

0.001

Description

Size of the step taken in reciprocal space to perform the numerical differentiation

EFG
Type

Block

Description

The electronic charge density causes an electric field, and the gradient of this field couples with the nuclear quadrupole moment, that some (non-spherical) nuclei have and can be measured by several spectroscopic techniques. The EFG tensor is the second derivative of the Coulomb potential at the nuclei. For each atom it is a 3x3 symmetric and traceless matrix. Diagonalization of this matrix gives three eigenvalues, which are usually ordered by their decreasing absolute size and denoted as V_{xx}, V_{yy}, V_{zz}. The result is summarized by the largest eigenvalue and the asymmetry parameter.

Enabled
Type

Bool

Default value

No

GUI name

EFG (electric field gradient): Calculate

Description

Compute the EFG tensor (for nuclear quadrupole interaction).

EigThreshold
Type

Float

Default value

0.01

Description

Threshold for printing the eigenvectors coefficients (Print Eigens)

ElectronHole
Type

Block

Description

Allows one to specify an occupied band which shall be depopulated, where the electrons are then moved to the Fermi level. For a spin-restricted calculation 2 electrons are shifted and for a spin-unrestricted calculation only one electron is shifted.

BandIndex
Type

Integer

Description

Which occupied band shall be depopulated.

SpinIndex
Type

Integer

Description

Defines the spin of the shifted electron (1 or 2).

EmbeddingPotential
Type

Block

Description

An external potential can be read in and will be added to the effective Kohn-Sham potential. It has to be on the becke grid

Filename
Type

String

Default value

Description

Name of the file containing the embedding potential.

PotentialName
Type

String

Default value

Description

Name of variable containing the potential.

EnforcedSpinPolarization
Type

Float

GUI name

Spin polarization

Description

Enforce a specific spin-polarization instead of occupying according to the aufbau principle. The spin-polarization is the difference between the number of alpha and beta electron. Thus, a value of 1 means that there is one more alpha electron than beta electrons. The number may be anything, including zero, which may be of interest when searching for a spin-flipped pair, that may otherwise end up in the (more stable) parallel solution.

ESR
Type

Block

Description

Zeeman g-tensor. The Zeeman g-tensor is implemented using two-component approach of Van Lenthe and co-workers in which the g-tensor is computed from a pair of spinors related to each other by time-reversal symmetry. Note: the following options are necessary for ESR: ‘Relativistic zora spin’ and ‘Kspace 1’

Enabled
Type

Bool

Default value

No

GUI name

ESR: g-tensor

Description

Compute Zeeman g-tensor. The Zeeman g-tensor is implemented using two-component approach of Van Lenthe and co-workers in which the g-tensor is computed from a pair of spinors related to each other by time-reversal symmetry. Note: the following options are necessary for ESR: ‘Relativistic zora spin’ and ‘Kspace 1’

Excitations
Type

Block

Description

Excitation energies: UV/Vis

Fermi
Type

Block

Description

Technical parameter used in determining the Fermi energy, which is carried out at each cycle of the SCF procedure.

Delta
Type

Float

Default value

0.0001

Description

Convergence criterion: upper and lower bounds for the Fermi energy and the corresponding integrated charge volumes must be equal within delta.

Eps
Type

Float

Default value

1e-10

Description

After convergence of the Fermi energy search procedure, a final estimate is defined by interpolation and the corresponding integrated charge volume is tested. It should be exact, to machine precision. Tested is that it deviates not more than eps.

MaxTry
Type

Integer

Default value

15

Description

Maximum number of attempts to locate the Fermi energy. The procedure is iterative in nature, narrowing the energy band in which the Fermi energy must lie, between an upper and a lower bound. If the procedure has not converged sufficiently within MaxTry iterations, the program takes a reasonable value and constructs the charge density by interpolation between the functions corresponding to the last used upper and lower bounds for the Fermi energy.

FermiSurface
Type

Block

Description

Calculation of the Fermi surface for metals

Enabled
Type

Bool

Default value

No

GUI name

Calculate Fermi surface

Description

Calculate the Fermi surface if the system has no band gap (i.e. is a metal). The result can be visualized with amsbands.

KIntegForSymmetricKGrid
Type

Integer

Default value

-1

Description

If the (default) regular k-grid is used, a symmetric one is created to determine the Fermi surface. If this key is not specified an automatic value of kInteg is used. Odd values trigger quadratic interpolation.

NMesh
Type

Integer

Default value

7

Description

Improves the matching of the interpolated quadratic surface. For better results it makes more sense to increaes KIntegForSymmetricKGrid.

FormFactors
Type

Integer

Default value

2

Description

Number of stars of K-vectors for which the form factors are computed

Fragment
Type

Block

Recurring

True

Description

Defines a fragment. You can define several fragments for a calculation.

AtomMapping
Type

Non-standard block

Description

Format ‘indexFragAt indexCurrentAt’. One has to associate the atoms of the fragment to the atoms of the current calculation. So, for each atom of the fragment the indexFragAt has to be associated uniquely to the indexCurrentAt for the current calculation.

Filename
Type

String

Description

Filename of the fragment. Absolute path or path relative to the executing directory.

Labels
Type

Non-standard block

Description

This gives the possibility to introduce labels for the fragment orbitals. See examples.

FuzzyPotential
Type

Non-standard block

Description

Atomic (fuzzy cell) based, external, electric potential. See example.

FuzzyUnitCellGrid
Type

Block

Description

Undocumented.

AtomRadiusLSG
Type

Float

Default value

0.0

Description

Undocumented.

CellPartitionDelta
Type

Float

Default value

4.0

Description

Undocumented.

CellPartitionInterpolationCubic
Type

Bool

Default value

No

Description

Undocumented.

CellPartitionInterpolationMesh
Type

Integer

Default value

100

Description

Undocumented.

CellPartitionVersion
Type

Integer

Default value

2

Description

Undocumented.

CentralizeNaturalLSG
Type

Bool

Default value

No

Description

Undocumented.

InterpolateCellPartition
Type

Bool

Default value

No

Description

Undocumented.

NumIntExtraL
Type

Integer

Default value

0

Description

Undocumented.

NumIntExtraRad
Type

Integer

Default value

0

Description

Undocumented.

PartitionFunctionTol
Type

Float

Default value

1e-08

Description

Undocumented.

PruneLatticeSummedGrid
Type

Bool

Default value

Yes

Description

Undocumented.

ReduceAccuracyLSG
Type

Bool

Default value

No

Description

Undocumented.

SimpleLatticeSummedGrid
Type

Bool

Default value

No

Description

Undocumented.

Grid
Type

Block

Description

Options for the regular grid used for plotting (e.g. density plot). Used ICW the restart option.

ExtendX
Type

Float

Default value

0.0

Unit

Bohr

Description

Extend the default regular grid along the x-direction by the specified amount: [x_min, x_max] => [x_min - ExtendX/2, x_max + ExtendX/2].

ExtendY
Type

Float

Default value

0.0

Unit

Bohr

Description

Extend the default regular grid along the y-direction by the specified amount: [y_min, y_max] => [y_min - ExtendY/2, y_max + ExtendY/2].

ExtendZ
Type

Float

Default value

0.0

Unit

Bohr

Description

Extend the default regular grid along the z-direction by the specified amount: [z_min, z_max] => [z_min - ExtendZ/2, z_max + ExtendZ/2].

FileName
Type

String

Default value

Description

Read in the grid from a file. The file format of the grid is: three numbers per line (defining the x, y and z coordinates of the points).

Type
Type

Multiple Choice

Default value

coarse

Options

[coarse, medium, fine]

Description

The default regular grids.

UserDefined
Type

Non-standard block

Description

Once can define the regular grid specification in this block. See example.

GridBasedAIM
Type

Block

Description

Invoke the ultra fast grid based Bader analysis.

Enabled
Type

Bool

Default value

No

GUI name

Bader (AIM): Atomic properties

Description

Invoke the ultra fast grid based Bader analysis.

Iterations
Type

Integer

Default value

40

Description

The maximum number of steps that may be taken to find the nuclear attractor for a grid point.

SmallDensity
Type

Float

Default value

1e-06

Description

Value below which the density is ignored. This should not be chosen too small because it may lead to unassignable grid points.

UseStartDensity
Type

Bool

Default value

No

Description

Whether the analysis is performed on the startup density (True) or on the final density (False).

GrossPopulations
Type

Non-standard block

Description

Partial DOS (pDOS) are generated for the gross populations listed under this key. See example.

GW
Type

Block

Description

Instruct ADF to perform a G0W0 calculation.

AdaptiveMixing
Type

Float List

Description

Requests to use adaptive mixing instead of DIIS and sets the staring mixing parameter for mixing of Green’s function in case of self-consistency. Adapative mixing is recommenened in case a qsGW calculation does not converge with DIIS. It is ignored in non-selfconsistent calculation and overwritten by DIIS when DIIS is also present.

AnalyticalIntegration
Type

Bool

Default value

No

Description

Use analytical integration to calculate the self-energy. Very slow, unless the system is very small but useful to check the accuracy of the frequency integration

Converge
Type

Block

Description

Sets convergence criteria for the GW calculation in self-consistent case

Density
Type

Float List

Default value

[1e-08, 1e-05]

Description

First Criterion for self-consistency procedure to terminate. Criterion is the trace of the density matrix. Ignored in non-selfconsistent Calculation and in eigenvalue self-consistent GW It is possible to run a qsGW calculation with an inner SCF loop which updates the static part of the elf-energy only. This can be useful to accelerate the convergence in case linear mixing is used. It is not recommended to use linear mixing, so it is also not recommened to use that inner loop as well. The second number in this list specifies the convergence criterion for the inner SCF loop.

HOMO
Type

Float

Default value

0.003

Unit

eV

GUI name

HOMO energy convergence

Description

Criterion for self-consistency procedure to terminate. The self-consistent GW calculation terminates, when the difference between the HOMO QP energies between 2 consecutive iterations is below this number. The LUMO energy converged faster than the HOMO energy so when the HOMO energy is converged according to this criterion, the LUMO energy will be converged as well. In non-selfconsistent Calculation, this criterion is ignored.

DIIS
Type

Integer

Default value

10

Description

Requests to use DIIS. This is the Default. Number of expansion coefficients can be requested as well. Ignored in non-selfconsistent calculation

Enabled
Type

Bool

Default value

No

GUI name

Calculate GW quasi-particle energies

Description

Enable the calculation of the GW quasi-particle energies.

FixedGrids
Type

Bool

Default value

No

Description

In a self-consistent GW calculation, do not recalculate Grids. Can be useful in case of convergence problems. Only relevant for qsGW and qsGW0. In case of evGW and evGW0, the grids are always kept fixed.

LinearMixing
Type

Float List

Description

Requests to use linear mixing instead of DIIS and sets the mixing parameter for linear mixing of Green’s function in case of self-consistency. It is ignored in non-selfconsistent calculation and overwritten by DIIS when DIIS is also present.

LinearizeQPequations
Type

Bool

Default value

No

Description

Instead of solving the non-linear QP equations in a G0W0 (or evGW calculation) by bisection exacly, linearize them by first-order Taylor expansion. This is not recommended since it does not save computational time when used together with analytical continuation (as implemented in AMS). It might however be useful for benchmarking or for validating results. If the results os the linearization differ by a lot (for instance, more than 0.1 eV in frontier QP energies) from the non-linearized results, this might indicate that the GW calculation is not reliable.

OffDiagonalEFermi
Type

Bool

Default value

No

Description

Analytically continue the off-diagonal elements of the KSF2 qsGW Hamiltonian at the Fermi-energy instead of omega=0. Typically leads to slightly lower QP energies, i.e. higher ionization potentials. The HOMO-LUMO gaps are typically not affected.

PrintAllSolutions
Type

Bool

Default value

No

Description

Print out all solutions for all requested states. Detects multiple solutions of the QP equations.

QPHamiltonian
Type

Multiple Choice

Default value

KSF2

Options

[KSF1, KSF2]

Description

The quasi-particle Hamiltonian can be constructed in different ways. KSF1 refers to the original construction by Kotani, Van Schilfgaarde anf Faleev (KSF) which is also implemented in TURBOMOLE. KSF2 refers to an alternative construction by KSF. KSF1 is not recommended since it is numerically less stable than KSF2. The results are typically very similar. The QP energies at which the matrix elements are evaluated can be tweaked further, see the two subsequent keys: However, KSF2 is recommened since it typically leads to QP energies with the best agreement with experiment. Ignored when not a quasi-particle self-consistent GW calculation is performed

ScissorShift
Type

Bool

Default value

No

Description

Only calculate the HOMO and LUMO QP energies and shift the remaining QP energies by the same amount. This is a rather crude approximaiton and not recommended. It might again be useful for benchmarking purposes.

SelfConsistency
Type

Multiple Choice

Default value

G0W0

Options

[G0W0, EVGW0, EVGW, QSGW0, QSGW]

Description

Sets the level of self-consistency in a GW calculation. G0W0 calculates a a on-shot, perturbative correction to the KS eigenvalues. In evGW and evGW0, the quasi-particle energies are updated until self-consistency is reached. evGW0 requests that the Green’s function is evaluated self-consistently but not the screened interaction. In qsGW, the density is updated as well, however, the self-energy is mapped to a static effective potential and the Dyson equation is solved by diagonalization instead of inversion. The results of a qsGW are independent of the choice of the underlying exchange-correlation functional and are usually the most accurate ones. The same is done in qsGW0, but the screened interaction is not updated.

SelfEnergy
Type

Multiple Choice

Default value

GW

Options

[HF, GW, G3W2]

Description

Controls the form of the self-energy. GW is the default and corresponds to the standard GW calculation. G3W2 is a GW calculation plus a perturbative second-order statically screened exchange correction (second order expansion in the self-energy). Note, that there the self-energy is always static.

nIterations
Type

Integer List

Default value

[10]

GUI name

Number of iterations

Description

The maximum number of iterations within the (partially or fully) self-consistent GW calculation has to converge. Ignored when Formalism is set to G0W0

nStates
Type

Integer

Default value

5

GUI name

N states

Description

Number of Quasiparticle States to be printed to output. The default is 5 states which in this case means that min(5, Number of particle states) occupied and min(5, Number of hole states) hole states are printed. The whole list of states can be printed by setting this parameter to -1’

preconditionQSGW
Type

Bool

Default value

No

Description

If true, the QSGW equations are solved but prior to each diagonalization, i.e. a G0W0 calculation is performed to find the optimal QP energies at which to analytically continue the self-energy. This is in principle a more consistent construction than KSF1 or KSF2 since the diagonal elements are consistent with G0W0. In KSF1 and KSF2, the diagonal elements are evaluated at the QP energies from the previous iteration which is equivalent to a zeroth-order Taylor expansion of the diaognal elements around the previous QP energies.Enabling this option typically leads to slightly lower QP energies.

HubbardU
Type

Block

Description

Options for Hubbard-corrected DFT calculations.

Atom
Type

Block

Recurring

True

Description

Specify Hubbard parameters (U,l) for a certain element

Element
Type

String

Description

Name of the element, such as Cu or Zn

LValue
Type

Multiple Choice

Default value

s

Options

[s, p, d, f]

Description

L value of the shell to apply the Hubbard model to

UValue
Type

Float

Default value

0.0

Description

Hubbard U value.

IgnoreForPEDA
Type

Bool

Default value

No

Description

Ignore the hubbard energy term when calculating the energy of psi_0 for the energy decomposition analysis (EDA).

LValue
Type

String

Default value

Description

For each atom type specify the l value (0 - s orbitals, 1 - p orbitals, 2 - d orbitals). A negative value is interpreted as no l-value.

PrintOccupations
Type

Bool

Default value

Yes

Description

Whether or not to print the occupations during the SCF.

Region
Type

Block

Recurring

True

Description

Specify Hubbard parameters (U,l) for all atoms in a certain region

LValue
Type

Multiple Choice

Default value

s

Options

[s, p, d, f]

Description

L value of the shell to apply the Hubbard model to

Name
Type

String

Description

Name of the region

UValue
Type

Float

Default value

0.0

Description

Hubbard U value.

UValue
Type

String

Default value

Description

For each atom type specify the U value (in atomic units). A value of 0.0 is interpreted as no U.

Integration
Type

Block

Description

Options for the Voronoi numerical integration scheme. Deprecated. Use BeckeGrid instead.

AccInt
Type

Float

Default value

3.5

Description

General parameter controlling the accuracy of the Voronoi integration grid. A value of 3 would be basic quality and a value of 7 would be good quality.

IntegrationMethod
Type

Multiple Choice

Default value

Becke

Options

[Becke, Voronoi]

Description

Choose the real-space numerical integration method. Note: the Voronoi integration scheme is deprecated.

KGrpX
Type

Integer

Default value

5

GUI name

Number of K-points at once

Description

Absolute upper bound on the number of k-points processed together. This only affects the computational performance.

KSpace
Type

Block

Description

Options for the k-space integration (i.e. the grid used to sample the Brillouin zone)

Quality
Type

Multiple Choice

Default value

Auto

Options

[Auto, GammaOnly, Basic, Normal, Good, VeryGood, Excellent]

GUI name

K-space

Description

Select the quality of the K-space grid used to sample the Brillouin Zone. If ‘Auto’, the quality defined in the ‘NumericalQuality’ will be used. If ‘GammaOnly’, only one point (the gamma point) will be used. The actual number of K points generated depends on this option and on the size of the unit cell. The larger the real space cell, the fewer K points will be generated. The CPU-time and accuracy strongly depend on this option.

Regular
Type

Block

Description

Options for the regular k-space integration grid.

NumberOfPoints
Type

Integer List

Description

Use a regular grid with the specified number of k-points along each reciprocal lattice vector. For 1D periodic systems you should specify only one number, for 2D systems two numbers, and for 3D systems three numbers.

Symmetric
Type

Block

Description

Options for the symmetric k-space integration grid.

KInteg
Type

Integer

GUI name

Accuracy

Description

Specify the accuracy for the Symmetric method. 1: absolutely minimal (only the G-point is used) 2: linear tetrahedron method, coarsest spacing 3: quadratic tetrahedron method, coarsest spacing 4,6,… (even): linear tetrahedron method 5,7…. (odd): quadratic method The tetrahedron method is usually by far inferior.

Type
Type

Multiple Choice

Default value

Regular

Options

[Regular, Symmetric]

GUI name

K-space grid type

Description

The type of k-space integration grid used to sample the Brillouin zone (BZ) used. ‘Regular’: simple regular grid. ‘Symmetric’: symmetric grid for the irreducible wedge of the first BZ (useful when high-symmetry points in the BZ are needed to capture the correct physics of the system, graphene being a notable example).

LDOS
Type

Block

Description

Local Density-Of-States information. This can be used to generate STM images in the Tersoff-Hamann approximation (see https://doi.org/10.1103/PhysRevB.31.805)

DeltaNeg
Type

Float

Default value

0.0001

Unit

Hartree

Description

Lower bound energy (Shift-DeltaNeg)

DeltaPos
Type

Float

Default value

0.0001

Unit

Hartree

Description

Upper bound energy (Shift+DeltaPos)

Shift
Type

Float

Default value

0.0

Unit

Hartree

Description

The energy bias with respect to the Fermi level.

MBPT
Type

Block

Description

Technical aspects of the MP2 algorithm.

Dependency
Type

Bool

Default value

Yes

Description

If true, to improve numerical stability, almost linearly-dependent combination of basis functions are removed from the Green’s function that are used in the MBPT equations. Disabeling this key is stringly discouraged. Its value can however be changed. The key to adjust this value is RiHartreeFock%DependencyThreshold

ExcludeCore
Type

Bool

Description

If active, excludes core states from the calculation of the optiomal imaginary time and frequency grids. The core states are still included in all parts of the calculations. In case a frozen care calculation is performed, this key is ignored. For MP2 and double hybrid calculation, it defauls to false. For RPA and GW calculations, it defaults to true.

FitSetQuality
Type

Multiple Choice

Default value

Auto

Options

[Auto, VeryBasic, Basic, Normal, Good, VeryGood]

Description

Specifies the fit set to be used in the MBPT calculation. ‘Normal’ quality is generally sufficient for basis sets up to and including TZ2P. For larger basis sets (or for benchmarking purposes) a ‘VeryGood’ fit set is recommended. Note that the FitSetQuality heavily influences the computational cost of the calculation. If not specified or ‘Auto’, the RIHartreeFock%FitSetQuality is used.

Formalism
Type

Multiple Choice

Default value

Auto

Options

[Auto, RI, LT, All]

Description

Specifies the formalism for the calculation of the MP2 correlation energy. ‘LT’ means Laplace Transformed MP2 (also referred to as AO-PARI-MP2), ‘RI’ means that a conventional RI-MP2 is carried out. If ‘Auto’, LT will be used in case of DOD double hybrids and SOS MP2, and RI will be used in all other cases. ‘All’ means that both RI and LT formalisms are used in the calculation. For a RPA or GW calculation, the formalism is always LT, irrespective of the formalism specified with this key.

FrequencyGridType
Type

Multiple Choice

Default value

LeastSquare

Options

[LeastSquare, GaussLegendre]

Description

Use Gauss-legendre grid for imaginary frequency integration in RPA and GW calculations instead of the usually used Least-Square optimized ones. Has the advantage that it can be systematically converged and an arbitrary number of grid points can be used. Typically more grid points will be needed to get the same level of accuracy. However, the convergence of the results with the size of the grid can be more systematic. These grids can only be used when Formalism is set to RI.

IntegrationQuality
Type

Multiple Choice

Options

[VeryBasic, Basic, Normal, Good, VeryGood]

Description

Specifies the integration quality to be used in the MBPT calculation. If not specified, the RIHartreeFock%IntegrationQuality is used.

SigmaFunctionalParametrization
Type

Multiple Choice

Default value

W1

Options

[W1, W2, S1, S2]

Description

Only relevant if a sigma-functional calculation is performed. Possible choices for the parametrization of the sigma-functional. Not all options are supported for all functionals.

ThresholdQuality
Type

Multiple Choice

Options

[VeryBasic, Basic, Normal, Good, VeryGood, Excellent]

Description

Controls the distances between atomic centers for which the product of two basis functions is not fitted any more. Especially for spatially extended, large systems, ‘VERYBASIC’ and ‘BASIC’ can lead to large computational savings, but the fit is also more approximate. If not specified, the RIHartreeFock%ThresholdQuality is used.

UseScaledZORA
Type

Bool

Default value

Yes

Description

If true, use the scaled ZORA orbital energies instead of the ZORA orbital energies in the MBPT equations.

nFrequency
Type

Integer

Default value

12

GUI name

Number of frequency points

Description

Number of imaginary frequency points. This key is only relevant for RPA and GW and will be ignored if used in an AO-PARI-MP2 calculation. 12 Points is the default for a RPA calculation. It is technically possible to use a different number of imaginary frequency points than for imaginary time. The maximum number of points which can be requested for imaginary frequency integration is 42. Important note: The computation time and memory requirements roughyl scale linearly with the number of imaginary frequency points. However, memory can be an issue for RPA and GW when the number of imaginary frequency points is high. In case a job crashes, it is advised to increase the number of nodes since the necessary memory distributes over all nodes.

nLambda
Type

Integer

Default value

1

GUI name

Number of lambda points

Description

Size of coupling constant integration grid for SOSEX variants in RPA. Default is 4 points

nTime
Type

Integer

GUI name

Number of time points

Description

Number of imaginary time points (only relevant in case the Laplace Transformed (LT) formalism is used). In the many-body-perturbation theory module in ADF, the polarizability (or Kohn-Sham density response function) is evaluated in imaginary time to exploit sparsity in the AO basis. For MP2, this is often referred to as a Laplace transform. For MP2, 9 points are the default. This is a safe choice, guaranteeing accuracies higher than 1 Kj/mol for most systems (For many simple organic systems, 6 points are sufficient for good accuracy). Only for systems with a very small HOMO-LUMO gap or low-lying core states (heavy elements starting from the 4th row of the periodic table) more points might be necessary. In principle, the same considerations apply for RPA and GW as well, however, the accuracy requirements are somewhat higher and 12 point are the default for RPA. In a GW calculation, the number of points is adjusted according to the numerical quality. Using less than 9 points is strongly discouraged except for the simplest molecules. In ADF2019, it can happen that the algorithm determining the imaginary time grid does not converge. In this case, the usual reason is that the number of points is too small and more points need to be specified. Starting from AMS2020, this does not happen any more. In case the imaginary time grid does not converge, the number of points is automatically adjusted until it does. The computation time of AO-PARI-MP2, RPA, and GW scales linearly with the number of imaginary time points.

MolecularNMR
Type

Block

Description

Options for the calculations of the NMR shielding tensor for molecules, excluding periodic systems. Implements the Schreckenbach method like ADF.

Enabled
Type

Bool

Default value

No

Description

Compute NMR shielding.

MultiSecantConfig
Type

Block

Description

Parameters for the Multi-secant SCF convergence method.

CMax
Type

Float

Default value

20.0

GUI name

Max coeff

Description

Maximum coefficient allowed in expansion

InitialSigmaN
Type

Float

Default value

0.1

GUI name

Initial

Description

This is a lot like a mix factor: bigger means bolder

MaxSigmaN
Type

Float

Default value

0.3

GUI name

Max

Description

Upper bound for the SigmaN parameter

MaxVectors
Type

Integer

Default value

20

GUI name

Number of cycles to use

Description

Maximum number of previous cycles to be used

MinSigmaN
Type

Float

Default value

0.01

GUI name

Min

Description

Lower bound for the SigmaN parameter

NEGF
Type

Block

Description

Options for the NEGF (non-equilibrium green function) transport calculation.

AlignChargeTol
Type

Float

Default value

0.1

Description

In an alignment run you want to get the number of electrons in the center right. This number specifies the criterion for that.

AlignmentFile
Type

String

Default value

Description

Band result file (.rkf) corresponding to the alignment calculation.

Alpha
Type

Float

Default value

1e-05

Description

A charge error needs to be translated in a potential shift. DeltaV = alpha * DeltaQ

ApplyShift1
Type

Bool

Default value

Yes

Description

Apply the main shift, obtained from comparing matrix elements in the leads with those from the tight-binding run. Strongly recommended.

ApplyShift2
Type

Bool

Default value

Yes

Description

Apply the smaller alignment shift. This requires an extra alignment run. Usually this shift is smaller.

AutoContour
Type

Bool

Default value

Yes

Description

Use automatic contour integral.

BiasPotential
Type

Float

Default value

0.0

Description

Apply a bias potential (atomic units). Can be negative. One has to specify the ramp potential with the FuzzyPotential key. This is mostly conveniently done with the GUI.

BoundOccupationMethod
Type

Integer

Default value

1

Description

See text. Only relevant with NonEqDensityMethod equal 2 or 3.

CDIIS
Type

Bool

Default value

No

Description

Make the normal DIIS procedure aware of the align charge error

CheckOverlapTol
Type

Float

Default value

0.01

Description

BAND checks how well the TB overlap matrix S(R=0) represents the overlap matrix in the lead region. Elements corresponding to the outer layer are neglected, because when using a frozen core they have bigger errors.

ContourQuality
Type

Multiple Choice

Default value

good

Options

[basic, normal, good, verygood]

Description

The density matrix is calculated numerically via a contour integral. Changing the quality influences the number of points. This influences a lot the performance.

DEContourInt
Type

Float

Default value

-1.0

Description

The energy interval for the contour grid. Defaults depends on the contour quality

DERealAxisInt
Type

Float

Default value

-1.0

Description

The energy interval for the real axis grid. Defaults depends on the contour quality.

DeltaPhi0
Type

Float

Default value

0.0

Description

Undocumented.

DeltaPhi1
Type

Float

Default value

0.0

Description

Undocumented.

DoAlignment
Type

Bool

Default value

No

Description

Set this to True if you want to do an align run. Between the leads there should be lead material. The GUI can be of help here.

EMax
Type

Float

Default value

5.0

Unit

eV

Description

The maximum energy for the transmission grid (with respect to the Fermi level of the lead)

EMin
Type

Float

Default value

-5.0

Unit

eV

Description

The minimum energy for the transmission grid (with respect to the Fermi level of the lead)

Eta
Type

Float

Default value

1e-05

Description

Small value used for the contour integral: stay at least this much above the real axis. This value is also used for the evaluation of the Transmission and dos.

IgnoreOuterLayer
Type

Bool

Default value

Yes

Description

Whether or not to ignore the outer layer.

KT
Type

Float

Default value

0.001

Description

k-Boltzman times temperature.

LeadFile
Type

String

Default value

Description

File containing the tight binding representation of the lead.

NE
Type

Integer

Default value

100

Description

The number of energies for the transmission energy grid.

NonEqDensityMethod
Type

Integer

Default value

1

Description

See text.

SGFFile
Type

String

Default value

Description

The result from the SGF program. Contains the Fermi energy of the lead.

YContourInt
Type

Float

Default value

0.3

Description

The density is calculated via a contour integral. This value specifies how far above the real axis the (horizontal part of the) contour runs. The value is rounded in such a way that it goes exactly halfway between two Fermi poles. There is a trade off: making it bigger makes the integrand more smooth, but the number of enclosed poles increases. For low temperatures it makes sense to lower this value, and use a smaller deContourInt.

YRealaxisInt
Type

Float

Default value

1e-05

Description

The non-Equilibrium density is calculated near the real axis.

NeutralizingDensity
Type

Multiple Choice

Default value

None

Options

[None, rho(atoms), rho(valence/atoms), rho(neutralizing/atoms), rho(homogeneous)]

Description

For charged systems an artificial compensating density can be used to make it neutral again. This fictitious density only affects the Coulomb potential. For charged periodic systems neutralization is required, as otherwise the Coulomb potential diverges.

NeutralizingDensityDetails
Type

Block

Description

DiffuseFactor
Type

Float

Default value

1.0

Description

The bigger this number, the more diffuse (extended) the neutralizing density becomes. Works only for rho(neutralizing/atoms)

HomogeneousDensity
Type

Block

Description

xxx

Origin
Type

Float List

Default value

[0.0, 0.0, 0.0]

Unit

Bohr

Description

Range
Type

Float

Default value

10.0

Unit

Bohr

Description

Width
Type

Float

Default value

1.0

Unit

Bohr

Description

NewResponse
Type

Block

Description

The TD-CDFT calculation to obtain the dielectric function is computed when this block is present in the input. Several important settings can be defined here.

ActiveESpace
Type

Float

Default value

5.0

Unit

eV

GUI name

Active energy space

Description

Modifies the energy threshold (DeltaE^{max}_{thresh} = omega_{high} + ActiveESpace) for which single orbital transitions (DeltaEpsilon_{ia} = Epsilon_{a}^{virtual} - Epsilon_{i}^{occupied}) are taken into account.

ActiveXYZ
Type

String

Default value

t

Description

Expects a string consisting of three letters of either ‘T’ (for true) or ‘F’ (for false) where the first is for the X-, the second for the Y- and the third for the Z-component of the response properties. If true, then the response properties for this component will be evaluated.

DensityCutOff
Type

Float

Default value

0.001

GUI name

Volume cutoff

Description

For 1D and 2D systems the unit cell volume is undefined. Here, the volume is calculated as the volume bordered by the isosurface for the value DensityCutoff of the total density.

EShift
Type

Float

Default value

0.0

Unit

eV

GUI name

Shift

Description

Energy shift of the virtual crystal orbitals.

FreqHigh
Type

Float

Default value

3.0

Unit

eV

Description

Upper limit of the frequency range for which response properties are calculated (omega_{high}).

FreqLow
Type

Float

Default value

1.0

Unit

eV

Description

Lower limit of the frequency range for which response properties are calculated. (omega_{low})

NFreq
Type

Integer

Default value

5

Description

Number of frequencies for which a linear response TD-CDFT calculation is performed.

NewResponseKSpace
Type

Block

Description

Modify the details for the integration weights evaluation in reciprocal space for each single-particle transition. Only influencing the NewResponse code.

Eta
Type

Float

Default value

1e-05

Description

Defines the small, finite imaginary number i*eta which is necessary in the context of integration weights for single-particle transitions in reciprocal space.

SubSimp
Type

Integer

Default value

3

Description

determines into how many sub-integrals each integration around a k point is split. This is only true for so-called quadratic integration grids. The larger the number the better the convergence behavior for the sampling in reciprocal space. Note: the computing time for the weights is linear for 1D, quadratic for 2D and cubic for 3D!

NewResponseSCF
Type

Block

Description

Details for the linear-response self-consistent optimization cycle. Only influencing the NewResponse code.

Bootstrap
Type

Integer

Default value

0

Description

defines if the Berger2015 kernel (Bootstrap 1) is used or not (Bootstrap 0). If you chose the Berger2015 kernel, you have to set NewResponseSCF%XC to ‘0’. Since it shall be used in combination with the bare Coulomb response only. Note: The evaluation of response properties using the Berger2015 is recommend for 3D systems only!

COApproach
Type

Bool

Default value

Yes

Description

The program automatically decides to calculate the integrals and induced densities via the Bloch expanded atomic orbitals (AO approach) or via the cyrstal orbitals (CO approach). The option COApproach overrules this decision.

COApproachBoost
Type

Bool

Default value

No

GUI name

CO Approach Boost

Description

Keeps the grid data of the Crystal Orbitals in memory. Requires significantly more memory for a speedup of the calculation. One might have to use multiple computing nodes to not run into memory problems.

Criterion
Type

Float

Default value

0.001

Description

For the SCF convergence the RMS of the induced density change is tested. If this value is below the Criterion the SCF is finished. Furthermore, one can find the calculated electric susceptibility for each SCF step in the output and can therefore decide if the default value is too loose or too strict.

DIIS
Type

Block

Description

Parameters influencing the DIIS self-consistency method

Enabled
Type

Bool

Default value

Yes

Description

If not enabled simple mixing without DIIS acceleration will be used.

MaxSamples
Type

Integer

Default value

10

Description

Specifies the maximum number of samples considered during the direct inversion of iteration of subspace (DIIS) extrapolation of the atomic charges during the SCC iterations. A smaller number of samples potentially leads to a more aggressive convergence acceleration, while a larger number often guarantees a more stable iteration. Due to often occurring linear dependencies within the set of sample vectors, the maximum number of samples is reached only in very rare cases.

MaximumCoefficient
Type

Float

Default value

10.0

Description

When the diis expansion coefficients exceed this threshold, the solution is rejected. The vector space is too crowded. The oldest vector is discarded, and the expansion is re-evaluated.

MinSamples
Type

Integer

Default value

-1

Description

When bigger than one, this affects the shrinking of the DIIS space on linear depence. It will not reduce to a smaller space than MinSamples unless there is extreme dependency.

MixingFactor
Type

Float

Default value

0.2

Description

The parameter used to mix the DIIS linear combination of previously sampled atomic charge vectors with an analogous linear combination of charge vectors resulting from population analysis combination. It can assume real values between 0 and 1.

LowFreqAlgo
Type

Bool

Default value

Yes

GUI name

Low Frequency Algorithm

Description

Numerically more stable results for frequencies lower than 1.0 eV. Note: for a graphene monolayer the conical intersection results in a very small band gap (zero band gap semi-conductor). This leads ta a failing low frequency algorithm. One can then chose to use the algoritm as originally proposed by Kootstra by setting the input value to *false*. But, this can result in unreliable results for frequencies lower than 1.0 eV!

NCycle
Type

Integer

Default value

20

GUI name

Cycles

Description

Number of SCF cycles for each frequency to be evaluated.

XC
Type

Integer

Default value

1

Description

Influences if the bare induced Coulomb response (XC 0) is used for the effective, induced potential or the induced potential derived from the ALDA kernel as well (XC 1).

NMR
Type

Block

Description

Options for the calculations of the NMR shielding tensor.

Correction_r
Type

Bool

Default value

Yes

Description

Undocumented.

Enabled
Type

Bool

Default value

No

Description

Compute NMR shielding.

MS0
Type

Float

Default value

0.01

Description

Undocumented.

NMRAtom
Type

Integer

Default value

0

Description

The index of the atom atom (in input order) for which NMR should be computed.

Numeric
Type

Bool

Default value

No

Description

Undocumented.

Original
Type

Bool

Default value

No

Description

Undocumented.

Print_jp
Type

Bool

Description

Print paramagnetic current.

SuperCell
Type

Bool

Default value

Yes

Description

This is the switch between the two methods, either the super cell (true), or the single-dipole method (false)

Test
Type

Bool

Description

Key for printing all intrinsic tensors.

Test_E
Type

Bool

Description

Test of energy levels.

Test_S
Type

Bool

Description

Test of overlap matrix.

UseSharedMemory
Type

Bool

Default value

Yes

Description

Whether or not to use shared memory in the NMR calculation.

NOCVdRhoPlot
Type

Non-standard block

Description

Goes together with the Restart%NOCVdRhoPlot and Grid keys. See example.

NOCVOrbitalPlot
Type

Non-standard block

Description

Goes together with the Restart%NOCVOrbitalPlot and Grid keys. See example.

NuclearModel
Type

Multiple Choice

Default value

PointCharge

Options

[PointCharge, Gaussian, Uniform]

Description

Specify what model to use for the nucleus. For the Gaussian model the nuclear radius is calculated according to the work of Visscher and Dyall (L. Visscher, and K.G. Dyall, Dirac-Fock atomic electronic structure calculations using different nuclear charge distributions, Atomic Data and Nuclear Data Tables 67, 207 (1997))

NUElstat
Type

Integer

Default value

50

Description

Number of outward (parabolic) integration points (for elliptical integration of the electrostatic interaction)

NumericalQuality
Type

Multiple Choice

Default value

Normal

Options

[Basic, Normal, Good, VeryGood, Excellent]

Description

Set the quality of several important technical aspects of a BAND calculation (with the notable exception of the basis set). It sets the quality of: BeckeGrid (numerical integration), ZlmFit (density fitting), KSpace (reciprocal space integration), and SoftConfinement (basis set confinement). Note: the quality defined in the block of a specific technical aspects supersedes the value defined in NumericalQuality (e.g. if I specify ‘NumericalQuality Basic’ and ‘BeckeGrid%Quality Good’, the quality of the BeckeGrid will be ‘Good’)

NVElstat
Type

Integer

Default value

80

Description

Number of angular (elliptic) integration points (for elliptical integration of the electrostatic interaction)

Occupations
Type

Non-standard block

Description

Allows one to input specific occupations numbers. Applies only for calculations that use only one k-point (i.e. pseudo-molecule calculations). See example.

OldResponse
Type

Block

Description

Options for the old TD-CDFT implementation.

Berger2015
Type

Bool

Default value

No

Description

Use the parameter-free polarization functional by A. Berger (Phys. Rev. Lett. 115, 137402). This is possible for 3D insulators and metals. Note: The evaluation of response properties using the Berger2015 is recommend for 3D systems only!

CNT
Type

Bool

Description

Use the CNT parametrization for the longitudinal and transverse kernels of the XC kernel of the homogeneous electron gas. Use this in conjunction with the NewVK option.

CNVI
Type

Float

Default value

0.001

Description

The first convergence criterion for the change in the fit coefficients for the fit functions, when fitting the density.

CNVJ
Type

Float

Default value

0.001

Description

the second convergence criterion for the change in the fit coefficients for the fit functions, when fitting the density.

Ebndtl
Type

Float

Default value

0.001

Unit

Hartree

Description

the energy band tolerance, for determination which routines to use for calculating the numerical integration weights, when the energy band posses no or to less dispersion.

Enabled
Type

Bool

Default value

No

Description

If true, the response function will be calculated using the old TD-CDFT implementation

Endfr
Type

Float

Default value

3.0

Unit

eV

Description

The upper bound frequency of the frequency range over which the dielectric function is calculated

Isz
Type

Integer

Default value

0

Description

Integer indicating whether or not scalar zeroth order relativistic effects are included in the TDCDFT calculation. 0 = relativistic effects are not included, 1 = relativistic effects are included. The current implementation does NOT work with the option XC%SpinOrbitMagnetization equal NonCollinear

Iyxc
Type

Integer

Default value

0

Description

integer for printing yxc-tensor (see http://aip.scitation.org/doi/10.1063/1.1385370). 0 = not printed, 1 = printed.

NewVK
Type

Bool

Description

Use the slightly modified version of the VK kernel (see https://aip.scitation.org/doi/10.1063/1.1385370). When using this option one uses effectively the static option, even for metals, so one should check carefully the convergence with the KSPACE parameter.

Nfreq
Type

Integer

Default value

5

Description

the number of frequencies for which a linear response TD-CDFT calculation is performed.

QV
Type

Bool

Description

Use the QV parametrization for the longitudinal and transverse kernels of the XC kernel of the homogeneous electron gas. Use this in conjunction with the NewVK option. (see reference).

Shift
Type

Float

Default value

0.0

Unit

eV

Description

energy shift for the virtual crystal orbitals.

Static
Type

Bool

Description

An alternative method that allows an analytic evaluation of the static response (normally the static response is approximated by a finite small frequency value). This option should only be used for non-relativistic calculations on insulators, and it has no effect on metals. Note: experience shows that KSPACE convergence can be slower.

Strtfr
Type

Float

Default value

1.0

Unit

eV

Description

is the lower bound frequency of the frequency range over which the dielectric function is calculated.

OrbitalPlot
Type

Non-standard block

Description

Goes together with the Restart%OrbitalPlot and Grid keys. See Example.

Output
Type

Block

Description

Control the output.

Print
Type

Block

Recurring

True

Description

Level
Type

Multiple Choice

Options

[None, Error, Warning, Minimal, Normal, Detail, TooMuchDetail]

Description

Section
Type

Multiple Choice

Options

[Prepare, SCF, Properties]

Description

OverlapPopulations
Type

Non-standard block

Description

Overlap population weighted DOS (OPWDOS), also known as the crystal orbital overlap population (COOP).

PEDA
Type

Bool

Default value

No

Description

If present in combination with the fragment block, the decomposition of the interaction energy between fragments is invoked.

PEDANOCV
Type

Block

Description

Options for the decomposition of the orbital relaxation (pEDA).

EigvalThresh
Type

Float

Default value

0.001

GUI name

Use NOCVs with ev larger than

Description

The threshold controls that for all NOCV deformation densities with NOCV eigenvalues larger than EigvalThresh the energy contribution will be calculated and the respective pEDA-NOCV results will be printed in the output

Enabled
Type

Bool

Default value

No

GUI name

Perform PEDA-NOCV analysis

Description

If true in combination with the fragment blocks and the pEDA key, the decomposition of the orbital relaxation term is performed.

PeriodicSolvation
Type

Block

Description

Additional options for simulations of periodic structures with solvation.

NStar
Type

Integer

Default value

4

Description

This option, expecting an integer number (>2), handles the accuracy for the construction of the COMSO surface. The larger the given number the more accurate the construction.

RemovePointsWithNegativeZ
Type

Bool

Default value

No

GUI name

Only above slab

Description

Whether the COSMO surface is constructed on both sides of a surface. If one is only interested in the solvation effect on the upper side of a surface (in the Z direction), then this option should be set to ‘True’

SymmetrizeSurfacePoints
Type

Bool

Default value

Yes

Description

Whether or not the COSMO point should be symmetrized

PopThreshold
Type

Float

Default value

0.01

Description

Threshold for printing Mulliken population terms. Works with ‘Print orbpop’

PotentialNoise
Type

Float

Default value

0.0001

Description

The initial potential for the SCF procedure is constructed from a sum-of-atoms density. Added to this is some small noise in the numerical values of the potential in the points of the integration grid. The purpose of the noise is to help the program break the initial symmetry, if that would lower the energy, by effectively inducing small differences between (initially) degenerate orbitals.

Print
Type

String

Recurring

True

Description

One or more strings (separated by blanks) from a pre-defined set may be typed after the key. This induces printing of various kinds of information, usually only used for debugging and checking. The set of recognized strings frequently changes (mainly expands) in the course of software-developments. Useful arguments may be symmetry, and fit.

Programmer
Type

Block

Description

Miscellaneous technical options.

SharedMemorySandwichingThreshold
Type

Integer

Default value

5000

Description

When the nr. of basis functions exceeds this threshold shared memory will be used to calculate matrix elements. Unless UseSharedMemoryForSandwiching is explicitly set in the input.

StoreDOSPerBas
Type

Bool

Default value

Yes

Description

Whether or not to store the parial DOS per basis function. This allows you to view any partial DOS with amsspectra and amsbands. Requires the CalcPDOS option to be on.

StoreOrbitals
Type

Bool

Default value

Yes

Description

Copy information on band.rkf needed for orbital plotting and restarts. This can be a lot of information. DOS and BandStructure require StoreOrbitals=true.

UseSharedMemoryForSandwiching
Type

Bool

Default value

Yes

GUI name

Use shared memory

Description

When calculating matrix elements the array will be shared. This saves memory at the cost of locking overhead. If not specified this will depend on the threshold SharedMemorySandwichingThreshold

UseTurnoverRuleForXcMatrix
Type

Bool

Default value

No

Description

Undocumented.

Usesharedmemory
Type

Bool

Default value

Yes

GUI name

Use shared memory

Description

When running more then one task, share memory between those tasks. This saves a lot of memory. Only disable it in case of problems.

PropertiesAtNuclei
Type

Non-standard block

Description

A number of properties can be obtained near the nucleus. An average is taken over a tiny sphere around the nucleus. The following properties are available: vxc[rho(fit)], rho(fit), rho(scf), v(coulomb/scf), rho(deformation/fit), rho(deformation/scf).

RadialDefaults
Type

Block

Description

Options for the logarithmic radial grid of the basis functions used in the subprogram Dirac

NR
Type

Integer

Default value

3000

Description

Number of radial points. With very high values (like 30000) the Dirac subprogram may not converge.

NRPerType
Type

Integer List

Description

If present overrides NR. The list needs to be as long as there are atom types

RMax
Type

Float

Default value

100.0

Unit

Bohr

Description

Upper bound of the logarithmic radial grid

RMin
Type

Float

Default value

1e-06

Unit

Bohr

Description

Lower bound of the logarithmic radial grid

RMinPerType
Type

Float List

Unit

Bohr

Description

If specified overrides RMin. The list needs to be as long as there are atom types (different elements)

Relativity
Type

Block

Description

Options for relativistic effects.

Level
Type

Multiple Choice

Default value

Scalar

Options

[None, Scalar, Spin-Orbit]

GUI name

Relativity (ZORA)

Description

None: No relativistic effects. Scalar: Scalar relativistic ZORA. This option comes at very little cost. SpinOrbit: Spin-orbit coupled ZORA. This is the best level of theory, but it is (4-8 times) more expensive than a normal calculation. Spin-orbit effects are generally quite small, unless there are very heavy atoms in your system, especially with p valence electrons (like Pb). See also the SpinOrbitMagnetization key.

ResponseInducedDensityPlot
Type

Non-standard block

Description

Goes together with Restart%ResponseInducedDensityPlot and Grid.

Restart
Type

Block

Description

Tells the program that it should restart with the restart file, and what to restart.

CheckAtomicPositions
Type

Bool

Default value

Yes

Description

If set to True: For restarting the SCF the atomic positions will be checked, and may not deviate too much.

DensityPlot
Type

Bool

Default value

No

Description

Goes together with the DensityPlot block and Grid blocks

File
Type

String

Default value

Description

Name of the restart file.

LoadEigenSystem
Type

Bool

Default value

No

GUI name

Load: eigen system

Description

At each step of the SCF load the section eigensystem from the restart file, forcing constant eigenvalues and vectors.

NOCVOrbitalPlot
Type

Bool

Default value

No

Description

Goes together with the NOCVOrbitalPlot and Grid blocks.

NOCVdRhoPlot
Type

Bool

Default value

No

Description

Goes together with the NOCVdRhoPlot and Grid blocks.

OrbitalPlot
Type

Bool

Default value

No

Description

Goes together with the OrbitalPlot and Grid

ResponseInducedDensityPlot
Type

Bool

Default value

No

Description

Goes together with the ResponseInducedDensityPlot and Grid blocks.

SCF
Type

Bool

Default value

No

GUI name

Restart: SCF

Description

Continue the SCF procedure using the orbital coefficients and occupations from the restart file.

UseDensityMatrix
Type

Bool

Default value

No

Description

If set to True: For restarting the SCF the density matrix will be used. Requires you to set ‘Save DensityMatrix’ in the previous run.

VTKFile
Type

String

Default value

Description

If specified a vtk file with be created with this name. If the extesion is ‘.txt’, a text file is created. Setting it to ‘CUBE’ one or more (one for each component) files in the cube format are generated with an automatic naming scheme.

VoronoiGrid
Type

Bool

Default value

No

Description

Copy the section Num In Params to the current file.

RIHartreeFock
Type

Block

Description

DependencyCoreRange
Type

Float

Description

Basis functions may be given a core character based on the range. For now only active in Band and only if present in the input

DependencyThreshold
Type

Float

Default value

0.001

Description

To improve numerical stability, almost linearly-dependent combination of basis functions are removed from the Hartree-Fock exchange matrix. If you obtain unphysically large bond energy in an Hybrid calculation, or an unphysically low correlation energy in an RPA, MP2, or double hybrid calculation, you might try setting the DependencyThreshold to a larger value (e.g. 3.0E-3) Note, that in GW calculations and GW-BSE calculations the default for this key is 5.0e-3.

FitGenerationDetails
Type

Block

Description

Technical details about how the RI Hartree-Fock fit functions are generated.

BoostL
Type

Bool

Default value

No

Description

Add extra max(l)+1 diffuse function When l denotes the highest angular momentum present in the primary basis, FromBasisProducts will generate auxiliary fit functions with up to 2l angular momentum. When this key is set to true, the maximum angular momentum in the auxiliary fit set becomes 2l+2. Typically, this option is not needed and when precision issues arise, it is rather advised to adjust the OneCenterDependencyThreshold key to a smaller value.

LapackWorkAround
Type

Bool

Default value

No

Description

GetFitFunctionsForAtomType diagonalization done with Lapack instead of Scalapack

Method
Type

Multiple Choice

Default value

Auto

Options

[Auto, FromBasisProducts]

Description

The way in which fit functions are generated. The main distinction is whether it depends on the basis functions used. When FromBasisProducts is used, the auxiliary basis is generated directly from the products of primary basis functions. This has the advantage that the auxiliary fit adapts automatcally to the basis set size. Especially for basis sets of QZ quality or larger, this is often necessary to otbain highly precise correlation energies using RPA or double hybrids FromBasisProducts option is also useful for GW or BSE calculations with basis sets of QZ quality or larger.

OneCenterDependencyThreshold
Type

Float

Default value

1e-08

Description

This key is only active when FromBasisProducts is chosen as method to generate the auxiliary basis. This threshold controls the size, and at the samw time, the precision of the auxiliary basis set. A smaller number leads to a larger auxiliary fit set. The default value of 1e-8 is typically sufficient to converge correlation energies and QP energies to a very high precision. It corresponds to an auxiliary basis which is typically 8-9 times larger than the primary basis.

UseBandRadialGrid
Type

Bool

Default value

Yes

Description

Only applies to band. The band logarithmic grid ranges (by default) from 1e-6 to 100 with 3000 points. Otherwise 300 points will be used. For 0-periodicity (molecules) it is advaisable to set this key to false since lots of memory is needed to evaluate all necessary integrals.

FitSetQuality
Type

Multiple Choice

Default value

Auto

Options

[Auto, VeryBasic, Basic, Normal, Good, VeryGood, Excellent, FromBasisProducts]

Description

The quality of auxiliary fit set employed in the RI scheme. If ‘Auto’, the value of the RIHartreeFock Quality option will be used. Normal quality is generally sufficient for basis sets up to and including TZ2P. For larger basis sets (or for benchmarking purposes) a VeryGood fit set is recommended. Note that the FitSetQuality heavily influences the computational cost of the calculation.

IntegrationQuality
Type

Multiple Choice

Options

[VeryBasic, Basic, Normal, Good, VeryGood, Excellent]

Description

Quality of the numerical integration for evaluating the integrals between basis functions and fit functions. If IntegrationQuality is not defined in input, the value defined in RIHartreeFock%Quality will be used.

Quality
Type

Multiple Choice

Default value

Auto

Options

[Auto, VeryBasic, Basic, Normal, Good, VeryGood, Excellent]

Description

Numerical accuracy of the RI procedure. If ‘Auto’, the quality specified in the ‘NumericalQuality’ will be used.

QualityPerRegion
Type

Block

Recurring

True

Description

Sets the fit-set quality for all atoms in a region. If specified, this overwrites the globally set quality.

Quality
Type

Multiple Choice

Options

[VeryBasic, Basic, Normal, Good, VeryGood, Excellent]

Description

This region’s quality of the auxiliary fit set employed in the RI scheme.

Region
Type

String

Description

The identifier of the region for which to set the quality.

ResponseQuality
Type

Multiple Choice

Options

[VeryBasic, Basic, Normal, Good, VeryGood, Excellent]

Description

Numerical accuracy of the RI procedure for the Response module.

ThresholdQuality
Type

Multiple Choice

Options

[VeryBasic, Basic, Normal, Good, VeryGood, Excellent]

Description

Linear scaling thresholds (also used for determining at what range the multiple approximation is used). To disable all linear scaling thresholds set this to Excellent.

UseMe
Type

Bool

Default value

Yes

Description

Set to False if you want to use the old RI scheme (ADF only)

Save
Type

String

Recurring

True

Description

Save scratch files or extra data that would be otherwise deleted at the end of the calculation. e.g. ‘TAPE10’ (containing the integration grid) or ‘DensityMatrix’

SCF
Type

Block

Description

Controls technical SCF parameters.

Eigenstates
Type

Bool

Description

The program knows two alternative ways to evaluate the charge density iteratively in the SCF procedure: from the P-matrix, and directly from the squared occupied eigenstates. By default the program actually uses both at least one time and tries to take the most efficient. If present, Eigenstates turns off this comparison and lets the program stick to one method (from the eigenstates).

Iterations
Type

Integer

Default value

300

GUI name

Maximum number of cycles

Description

The maximum number of SCF iterations to be performed.

Method
Type

Multiple Choice

Default value

MultiStepper

Options

[DIIS, MultiSecant, MultiStepper]

Description

Choose the general scheme used to converge the density in the SCF. In case of scf problems one can try the MultiSecant alternative at no extra cost per SCF cycle. For more details see the DIIS and MultiSecantConfig block.

Mixing
Type

Float

Default value

0.075

Description

Initial ‘damping’ parameter in the SCF procedure, for the iterative update of the potential: new potential = old potential + mix (computed potential-old potential). Note: the program automatically adapts Mixing during the SCF iterations, in an attempt to find the optimal mixing value.

MultiStepperPresetPath
Type

String

Default value

DFTB/default2023.inc

Description

Name of file containing a SCFMultiStepper key block. This will be used if no Explicit SCFMultiStepper block is in the input, and Method=MultiStepper. If the path is not absolute, it is relative to $AMSHOME/data/presets/multi_stepper’

PMatrix
Type

Bool

Description

If present, evaluate the charge density from the P-matrix. See also the key Eigenstates.

PrintAllOccupiedBands
Type

Bool

Default value

No

Description

When printing the ranges of the bands, include all occupied ones.

PrintAllVirtualBands
Type

Bool

Default value

No

Description

When printing the ranges of the bands, include all virtual ones.

PrintAlwaysBandRanges
Type

Bool

Default value

No

Description

Normally the ranges of the bands are only printed at the last SCF cycle

Rate
Type

Float

Default value

0.99

Description

Minimum rate of convergence for the SCF procedure. If progress is too slow the program will take measures (such as smearing out occupations around the Fermi level, see key Degenerate of block Convergence) or, if everything seems to fail, it will stop

SCFMultiStepper
Type

Block

Description

To solve the self-consistent problem multiple steppers can be tried during stints using the ones that give the best progress.

AlwaysChangeStepper
Type

Bool

Default value

No

Description

When the progress is fine there is no reason to change the stepper. In practice this is always set to true, because also the Stepper%ExpectedSlope can be used to achieve similar behavior.

ErrorGrowthAbortFactor
Type

Float

Default value

1000.0

Description

Abort stint when the error grows too much, compared to the error at the start of the stint.

FractionalStepFactor
Type

Float

Default value

-1.0

Description

Multiply the step by this factor. If smaller than zero this is not used.

MinStintCyclesForAbort
Type

Integer

Default value

0

Description

Look at ErrorGrowthAbortFactor only when a number of steps has been completed since the start of the stint. A value of 0 means always.

Stepper
Type

Block

Recurring

True

Description

??

AbortSlope
Type

Float

Default value

100.0

Description

If the slope (at the end of a stint) is larger than this: abort the stepper

DIISStepper
Type

Block

Description

DIIS stepper

EDIISAlpha
Type

Float

Default value

0.01

Description

The extra energy vector is weighed by this factor. .

MaxCoefficient
Type

Float

Default value

20.0

Description

The largest allowed value of the expansion coefficients. If exceed the number of vectors is reduces until the criterion is met.

MaxVectors
Type

Integer

Default value

10

Description

Maximum number of previous densities to be used (size of the history).

MinVectors
Type

Integer

Default value

-1

Description

Try to prevent to make nVectors shrink below this value, by allowing for significantly larger coefficents.

Mix
Type

Float

Default value

0.2

Description

Also known as greed. It determines the amount of output density to be used. May be changed by the MixAdapter.

ErrorGrowthAbortFactor
Type

Float

Default value

-1.0

Description

Abort stint when the error grows too much, compared to the error at the start of the stint. Overides global ErrorGrowthAbortFactor when set to a value > 0

ExpectedSlope
Type

Float

Default value

-100.0

Description

If the slope of the total SCF is better than this keep on going.

FractionalStepFactor
Type

Float

Default value

-1.0

Description

Multiply the step by this factor. If smaller than zero this is not used.

MaxInitialError
Type

Float

Description

Only use the stepper when error is smaller than this.

MaxIterationNumber
Type

Integer

Default value

-1

Description

Stepper will only be active for iterations smaller than this number. (Negative value means: Ignore this option)

MaxStintNumber
Type

Integer

Default value

-1

Description

Stepper will only be active for stints smaller than this number. (Negative value means: Ignore this option)

MinInitialError
Type

Float

Description

Only use the stepper when error is larger than this.

MinIterationNumber
Type

Integer

Default value

-1

Description

Stepper will only be active for iterations larger than this number.

MinStintCyclesForAbort
Type

Integer

Default value

0

Description

Look at ErrorGrowthAbortFactor only when a number of steps has been completed since the start of the stint. A value of 0 means always. Overides global value.

MinStintNumber
Type

Integer

Default value

-1

Description

Stepper will only be active for stints larger than this number.

MixAdapter
Type

Block

Description

Generic mix adapter

ErrorGrowthPanicFactor
Type

Float

Default value

10.0

Description

When the error increases more than this factor, this mix is reduced a lot.

GrowthFactor
Type

Float

Default value

1.1

Description

When the mix is considered too low it is multiplied by this factor. Otherwise it is divided by it.

MaxMix
Type

Float

Default value

0.3

Description

Do not grow the mix above this value.

MinMix
Type

Float

Default value

0.1

Description

Do not shrink the mix below this value.

NTrialMixFactors
Type

Integer

Default value

3

Description

Only used with Type=Trial. Must be an odd number.

TrialMode
Type

Multiple Choice

Default value

CurrentMixCentered

Options

[CurrentMixCentered, FullRange]

Description

How are the NTrialMixFactors chosen?

Type
Type

Multiple Choice

Default value

Error

Options

[Error, Energy, UnpredictedStep, Trial]

Description

Adapt the mix factor based on the observed progress (slope).

MixStepper
Type

Block

Description

Simple mixing stepper, only using the previous (in/out) denstity.

Mix
Type

Float

Default value

0.1

Description

???.

MultiSecantStepper
Type

Block

Description

Multi secant stepper.

MaxCoefficient
Type

Float

Default value

20.0

Description

???.

MaxVectors
Type

Integer

Default value

10

Description

???.

Mix
Type

Float

Default value

0.2

Description

???.

Variant
Type

Multiple Choice

Default value

MSB2

Options

[MSB1, MSB2, MSR1, MSR1s]

Description

There are several version of the Multi secant method.

StintLength
Type

Integer

Description

Override global StintLength.

StintLength
Type

Integer

Default value

10

Description

A stepper is active during a number of SCF cycles, called a stint.

UsePreviousStintForErrorGrowthAbort
Type

Bool

Default value

No

Description

The error is normally checked against the first error of the stint. With this option that will be the one from the previous stint, if performed with the same stepper.

VSplit
Type

Float

Default value

0.05

Description

To disturb degeneracy of alpha and beta spin MOs the value of this key is added to the beta spin potential at the startup.

Screening
Type

Block

Description

For the periodic solvation potential and for the old (not default anymore) fitting method, BAND performs lattice summations which are in practice truncated. The precision of the lattice summations is controlled by the options in this block.

CutOff
Type

Float

Description

Criterion for negligibility of tails in the construction of Bloch sums. Default depends on Accuracy.

DMadel
Type

Float

Description

One of the parameters that define the screening of Coulomb-potentials in lattice sums. Depends by default on Accuracy, rmadel, and rcelx. One should consult the literature for more information

NoDirectionalScreening
Type

Bool

Description

Real space lattice sums of slowly (or non-) convergent terms, such as the Coulomb potential, are computed by a screening technique. In previous releases, the screening was applied to all (long-range) Coulomb expressions. Screening is only applied in the periodicity directions. This key restores the original situation: screening in all directions

RCelx
Type

Float

Description

Max. distance of lattice site from which tails of atomic functions will be taken into account for the Bloch sums. Default depends on Accuracy.

RMadel
Type

Float

Description

One of the parameters that define screening of the Coulomb potentials in lattice summations. Depends by default on Accuracy, dmadel, rcelx. One should consult the literature for more information.

SelectedAtoms
Type

Integer List

Description

With this key you can select atoms. This has an effect on a few of options, like NMR and EFG.

Skip
Type

String

Recurring

True

Description

Skip the specified part of the Band calculation (expert/debug option).

SoftConfinement
Type

Block

Description

In order to make the basis functions more compact, the radial part of the basis functions is multiplied by a Fermi-Dirac (FD) function (this ‘confinement’ is done for efficiency and numerical stability reasons). A FD function goes from one to zero, controlled by two parameters. It has a value 0.5 at Radius, and the decay width is Delta.

Delta
Type

Float

Unit

Bohr

Description

Explicitely specify the delta parameter of the Fermi-Dirac function (if not specified, it will be 0.1*Radius).

Quality
Type

Multiple Choice

Default value

Auto

Options

[Auto, Basic, Normal, Good, VeryGood, Excellent]

GUI name

Confinement

Description

In order to make the basis functions more compact, the radial part of the basis functions is multiplied by a Fermi-Dirac (FD) function (this ‘confinement’ is done for efficiency and numerical stability reasons). A FD function goes from one to zero, controlled by two parameters. It has a value 0.5 at Radius, and the decay width is Delta. This key sets the two parameters ‘Radius’ and ‘Delta’. Basic: Radius=7.0, Delta=0.7; Normal: Radius=10.0, Delta=1.0; Good: Radius=20.0, Delta=2.0; VeryGood and Excellent: no confinement at all. If ‘Auto’, the quality defined in the ‘NumericalQuality’ will be used.

Radius
Type

Float

Unit

Bohr

Description

Explicitely specify the radius parameter of the Fermi-Dirac function.

Solvation
Type

Block

Description

Options for the COSMO (Conductor like Screening Model) solvation model.

CVec
Type

Multiple Choice

Default value

EXACT

Options

[EXACT, FITPOT]

GUI name

Calculate Coulomb interaction

Description

Choose how to calculate the Coulomb interaction matrix between the molecule and the point charges on the surface: - EXACT: use exact density, and integrate against the potential of the point charges. This may have inaccuracies when integration points are close to the point charges. - FITPOT: evaluate the molecular potential at the positions of the point charges, and multiply with these charges.

Charge
Type

Block

Description

Select the algorithm to determine the charges.

Conv
Type

Float

Default value

1e-08

Description

Charge convergence threshold in iterative COSMO solution.

Corr
Type

Bool

Default value

Yes

GUI name

Correct for outlying charge

Description

Correct for outlying charge.

Iter
Type

Integer

Default value

1000

Description

Maximum number of iterations to solve COSMO equations.

Method
Type

Multiple Choice

Default value

CONJ

Options

[CONJ, INVER]

GUI name

Charge determination method

Description

INVER: matrix inversion, CONJ: biconjugate gradient method. The CONJ method is guaranteed to converge with small memory requirements and is normally the preferred method.

Enabled
Type

Bool

Default value

No

GUI name

Include COSMO solvation

Description

Use the Conductor like Screening Model (COSMO) to include solvent effects.

Radii
Type

Non-standard block

Description

The values are the radii of the atomic spheres. If not specified the default values are those by Allinge. Format: ‘AtomType value’. e.g.: ‘H 0.7’

SCF
Type

Multiple Choice

Default value

VAR

Options

[VAR, PERT, NONE]

GUI name

Handle charges

Description

Determine the point charges either Variational (VAR) or after the SCF as a Perturbation (PERT).

Solvent
Type

Block

Description

Solvent details

Del
Type

Float

Description

Del is the value of Klamt’s delta_sol parameter, only relevant in case of Klamt surface.

Emp
Type

Float

Description

Emp is the empirical scaling factor x for the energy scaling.

Eps
Type

Float

Description

User-defined dielectric constant of the solvent (overrides the Eps value of the solvent defined in ‘Name’)

Name
Type

Multiple Choice

Default value

Water

Options

[AceticAcid, Acetone, Acetonitrile, Ammonia, Aniline, Benzene, BenzylAlcohol, Bromoform, Butanol, isoButanol, tertButanol, CarbonDisulfide, CarbonTetrachloride, Chloroform, Cyclohexane, Cyclohexanone, Dichlorobenzene, DiethylEther, Dioxane, DMFA, DMSO, Ethanol, EthylAcetate, Dichloroethane, EthyleneGlycol, Formamide, FormicAcid, Glycerol, HexamethylPhosphoramide, Hexane, Hydrazine, Methanol, MethylEthylKetone, Dichloromethane, Methylformamide, Methypyrrolidinone, Nitrobenzene, Nitrogen, Nitromethane, PhosphorylChloride, IsoPropanol, Pyridine, Sulfolane, Tetrahydrofuran, Toluene, Triethylamine, TrifluoroaceticAcid, Water]

GUI name

Solvent

Description

Name of a pre-defined solvent. A solvent is characterized by the dielectric constant (Eps) and the solvent radius (Rad).

Rad
Type

Float

Unit

Angstrom

Description

User-defined radius of the solvent molecule (overrides the Rad value of the solvent defined in ‘Name’).

Surf
Type

Multiple Choice

Default value

Delley

Options

[Delley, Wsurf, Asurf, Esurf, Klamt]

GUI name

Surface type

Description

Within the COSMO model the molecule is contained in a molecule shaped cavity. Select one of the following surfaces to define the cavity: - Wsurf: Van der Waals surface - Asurf: solvent accessible surface - Esurf: solvent excluding surface - Klamt: Klamt surface - Delley: Delley surface.

SolvationSM12
Type

Block

Description

Options for Solvation Model 12 (SM12).

ARO
Type

Float

Default value

0.0

Description

Square of the fraction of non-hydrogen atoms in the solvent that are aromatic carbon atoms (carbon aromaticity)

Acid
Type

Float

Default value

0.82

Description

Abraham hydrogen bond acidity parameter

Base
Type

Float

Default value

0.35

Description

Abraham hydrogen bond bacicity parameter

BornC
Type

Float

Default value

3.7

Description

Coulomb constant for General Born Approximation

BornRadiusConfig
Type

Block

Description

MaxCellDistance
Type

Float

Default value

30.0

Unit

Bohr

Description

Max distance from the centra cell used when computing the Born radii for periodic systems

PointsPerBohr
Type

Integer

Default value

10

Description

UseLegendreGrid
Type

Bool

Default value

Yes

Description

Chgal
Type

Float

Default value

2.474

Description

Exponential of Pauli’s bond order

Cust
Type

String

Description

Custom solvent input

Debug
Type

String

Description

Prints a lot of information about every pass on CDS and ENP code, keywords: ENP, CDS

EPS
Type

Float

Default value

78.36

Description

The dielectric constant

Enabled
Type

Bool

Default value

No

GUI name

Include SM12 solvation

Description

Whether to use the Solvation Model 12 (SM12) in the calculation.

HALO
Type

Float

Default value

0.0

Description

Square of the fraction of non-hydrogen atoms in the solvent molecule that are F, Cl, or Br (electronegative halogenicity)

Kappa
Type

Float

Default value

0.0

Description

Factor for Debye screening

PostSCF
Type

Bool

Default value

No

Description

Whether to apply the solvation potential during the SCF or only calculate the solvation energy after the SCF.

PrintSM12
Type

Bool

Default value

No

Description

Prints out an in-depth breakdown of solvation energies

RadSolv
Type

Float

Default value

0.4

Description

The radius distance between the solute and solvent

Ref
Type

Float

Default value

1.3328

Description

Refractive index of solvent

Solv
Type

Multiple Choice

Default value

WATER

Options

[ACETICACID, ACETONITRILE, ACETOPHENONE, ANILINE, ANISOLE, BENZENE, BENZONITRILE, BENZYLALCOHOL, BROMOBENZENE, BROMOETHANE, BROMOFORM, BROMOOCTANE, N-BUTANOL, SEC-BUTANOL, BUTANONE, BUTYLACETATE, N-BUTYLBENZENE, SEC-BUTYLBENZENE, T-BUTYLBENZENE, CARBONDISULFIDE, CARBONTETRACHLORIDE, CHLOROBENZENE, CHLOROFORM, CHLOROHEXANE, M-CRESOL, CYCLOHEXANE, CYCLOHEXANONE, DECALIN, DECANE, DECANOL, 1-2-DIBROMOETHANE, DIBUTYLETHER, O-DICHLOROBENZENE, 1-2-DICHLOROETHANE, DIETHYLETHER, DIISOPROPYLETHER, N-N-DIMETHYLACETAMIDE, N-N-DIMETHYLFORMAMIDE, 2-6-DIMETHYLPYRIDINE, DIMETHYLSULFOXIDE, DODECANE, ETHANOL, ETHOXYBENZENE, ETHYLACETATE, ETHYLBENZENE, FLUOROBENZENE, 1-FLUORO-N-OCTANE, HEPTANE, HEPTANOL, HEXADECANE, HEXADECYLIODIDE, HEXANE, HEXANOL, IODOBENZENE, ISOBUTANOL, ISOOCTANE, ISOPROPANOL, ISOPROPYLBENZENE, P-ISOPROPYLTOLUENE, MESITYLENE, METHANOL, METHOXYETHANOL, METHYLENECHLORIDE, N-METHYLFORMAMIDE, 2-METHYLPYRIDINE, 4-METHYL-2-PENTANONE, NITROBENZENE, NITROETHANE, NITROMETHANE, O-NITROTOLUENE, NONANE, NONANOL, OCTANE, OCTANOL, PENTADECANE, PENTANE, PENTANOL, PERFLUOROBENZENE, PHENYLETHER, PROPANOL, PYRIDINE, TETRACHLOROETHENE, TETRAHYDROFURAN, TETRAHYDROTHIOPHENEDIOXIDE, TETRALIN, TOLUENE, TRIBUTYLPHOSPHATE, TRIETHYLAMINE, 1-2-4-TRIMETHYLBENZENE, UNDECANE, WATER, XYLENE, 1-2-DIBROMOETHANE_WATER, 1-2-DICHLOROETHANE_WATER, BENZENE_WATER, CARBONTETRACHLORIDE_WATER, CHLOROBENZENE_WATER, CHLOROFORM_WATER, CYCLOHEXANE_WATER, DIBUTYLETHER_WATER, DIETHYLETHER_WATER, ETHYLACETATE_WATER, HEPTANE_WATER, HEXANE_WATER, NITROBENZENE_WATER, OCTANOL_WATER]

GUI name

Solvent

Description

List of predefined solvents

Tens
Type

Float

Default value

103.62

Description

Macroscopic surface tension of the solvent at the air/solvent interface at 298K (cal*mol^-1*Ang^-2)

TopologicalExtrapolation
Type

Block

Description

Method to extrapolate the long range Coulomb potential, needed for periodic calculations

FirstCell
Type

Integer

Default value

5

Description

First cell for the topological extrapolation of the long range part of the Coulomb Potential.

LastCell
Type

Integer

Default value

10

Description

Last cell for the topological extrapolation of the long range part of the Coulomb Potential.

Order
Type

Integer

Default value

3

Description

Order of the topological extrapolation of the long range part of the Coulomb Potential.

StopAfter
Type

String

Default value

BAND

Description

Specifies that the program is stopped after execution of a specified program-part (subroutine).

StoreHamAsMol
Type

Bool

Default value

No

Description

Undocumented, used for (at least) NEGF.

StoreHamiltonian
Type

Bool

Description

Undocumented.

StoreHamiltonian2
Type

Bool

Default value

No

Description

determine the tight-binding representation of the overlap an fock matrix. Used for (at least) NEGF.

StrainDerivatives
Type

Block

Description

Undocumented.

Analytical
Type

Bool

Description

Whether or not to use analytical strain derivatives. By default this is determined automatically, and used if possible.

AnalyticalElectrostatic
Type

Bool

Default value

No

Description

Undocumented.

Analyticalkinetic
Type

Bool

Default value

No

Description

Undocumented.

Analyticalpulay
Type

Bool

Default value

No

Description

Undocumented.

Analyticalxc
Type

Bool

Default value

No

Description

Undocumented.

Celltopoorder
Type

Integer

Default value

20

Description

Undocumented.

Coreorthoption
Type

Integer

Default value

2

Description

Undocumented.

Fitrho0numintextrarad
Type

Integer

Default value

0

Description

Undocumented.

Fitrho0prune
Type

Bool

Default value

Yes

Description

Undocumented.

Kinviadagger
Type

Bool

Default value

No

Description

Undocumented.

Lmaxmultipoleexpansion
Type

Integer

Default value

4

Description

Undocumented.

Naiveelstat
Type

Bool

Default value

No

Description

Undocumented.

Numericaldefdef
Type

Bool

Default value

Yes

Description

Undocumented.

Numericaldefdeflong
Type

Bool

Default value

No

Description

Undocumented.

Pairgridlowerangularorder
Type

Integer

Default value

5

Description

Undocumented.

Pairgridradpointsincrease
Type

Integer

Default value

0

Description

Undocumented.

Renormalizechargefitrho0
Type

Bool

Default value

No

Description

Undocumented.

Shiftmultipoleorigin
Type

Bool

Default value

Yes

Description

Undocumented.

Skipinlgwsmodule
Type

Bool

Default value

Yes

Description

Undocumented.

Subtractatomicxc
Type

Bool

Default value

No

Description

Undocumented.

Usesymmetry
Type

Bool

Default value

No

Description

Undocumented.

Usevstrainderrho
Type

Bool

Default value

No

Description

Undocumented.

fitrho0numintextral
Type

Integer

Default value

0

Description

Undocumented.

SubSymmetry
Type

Integer List

Description

The indices of the symmetry operators to maintain.

Tails
Type

Block

Description

Ignore function tails.

Bas
Type

Float

Default value

1e-06

GUI name

Basis functions

Description

Cut off the basis functions when smaller than the specified threshold.

Title
Type

String

Default value

Description

Title of the calculation, which will be printed in the output file.

Unrestricted
Type

Bool

Default value

No

Description

Controls whether Band should perform a spin-unrestricted calculation. Spin-unrestricted calculations are computationally roughly twice as expensive as spin-restricted.

UnrestrictedOnlyReference
Type

Bool

Default value

No

Description

Undocumented.

UnrestrictedReference
Type

Bool

Default value

No

Description

Undocumented.

UnrestrictedStartup
Type

Bool

Default value

No

Description

Undocumented.

UseInversionSymmetryInReciprocalSpace
Type

Bool

Default value

Yes

Description

Whether to use inversion symmetry in reciprocal space. This is almost always a valid assumption.

UseSymmetry
Type

Bool

Default value

Yes

Description

Whether or not to exploit symmetry during the calculation.

XC
Type

Block

Description

Exchange Correlation functionals

DFTHalf
Type

Block

Description

DFT-1/2 method for band gaps. See PRB vol 78,125116 2008. This method can be used in combination with any functional. For each active atom type (see ActiveAtomType) Band will perform SCF calculations at different screening cut-off values (see ScreeningCutOffs) and pick the cut-off value that maximizes the band gap. If multiple atom types are active, the screening cut-off optimizations are done one type at the time (in the same order as the ActiveAtomType blocks appear in the input).

ActiveAtomType
Type

Block

Recurring

True

Description

Use the DFT-1/2 method for the atom-type specified in this block.

AtomType
Type

String

Description

Atom-type to use. You can activate all atom-types by specifying ‘All’.

IonicCharge
Type

Float

Default value

0.5

Description

The amount of charge to be removed from the atomic HOMO.

ScreeningCutOffs
Type

Float List

Default value

[0.0, 1.0, 2.0, 3.0, 4.0, 5.0]

Unit

Bohr

Description

List of screening cut-offs (to screen the asymptotic IonicCharge/r potential). Band will loop over these values and find the cut-off that maximizes the band-gap. If only one number is provided, Band will simply use that value.

Enabled
Type

Bool

Default value

No

GUI name

Use method

Description

Whether the DFT-1/2 method will be used.

Prepare
Type

Bool

Default value

No

Description

Analyze the band structure to determine reasonable settings for an DFT-1/2 calculation. If this is possible the list of active atom types is written to the output. This can be used in a next run as the values for ActiveAtomType. The DFTHalf%Enabled key should be set to false

SelfConsistent
Type

Bool

Default value

Yes

Description

Apply the extra potential during the SCF, or only afterwards. Applying DFT-1/2 only post SCF increases the band gap, compared to the self-consistent one.

DoubleHybrid
Type

String

Description

Specifies the double hybrid functional that should be used during the SCF.

EmpiricalScaling
Type

Multiple Choice

Default value

None

Options

[None, SOS, SCS, SCSMI]

Description

Calculate the (SOS/SCS/SCSMI)-MP2 correlation energy.

GLLBKParameter
Type

Float

Default value

0.382

Description

K parameter for the GLLB functionals. See equation (20) of the paper.

HartreeFock
Type

Bool

Default value

No

Description

Stand alone HF calculation.

MP2
Type

Bool

Default value

No

Description

Calculate the MP2 correlation energy after the HF SCF is completed.

RPA
Type

Multiple Choice

Default value

None

Options

[None, Direct, Sigma, SOSEX, SOSSX]

Description

Specifies that RPA is used an possibly also a post-RPA method. By default, direct RPA is used

diracgga
Type

String

Default value

Description

GGA for the dirac .

dispersion
Type

String

Default value

DEFAULT

Description

The dispersion correction model to be used.

gga
Type

String

Default value

NONE

Description

GGA XC functional.

lda
Type

String

Default value

VWN

Description

LDA XC functional.

libxc
Type

String

Default value

NONE

Description

Functional using the LicXC library.

libxcdensitythreshold
Type

Float

Default value

1e-10

Description

Density threshold for LibXC functionals.

metagga
Type

String

Default value

NONE

Description

MetaGG XC functional.

model
Type

String

Default value

LB94

Description

Model potential. The possible choices are LB94, GLLB-SC, BGLLB-VWN, and BGLLB-LYP

spinorbitmagnetization
Type

String

Default value

collinearz

Description

Type of Spin-Orbit magnetization.

tb_mbjafactor
Type

Float

Default value

-1.23456789

Description

a parameter for the TB-MBJ model potential.

tb_mbjbfactor
Type

Float

Default value

-1.23456789

Description

b parameter for the TB-MBJ model potential..

tb_mbjcfactor
Type

Float

Default value

-1.23456789

Description

c parameter for the TB-MBJ model potential..

tb_mbjefactor
Type

Float

Default value

-1.23456789

Description

e parameter for the TB-MBJ model potential..

usexcfun
Type

Bool

Default value

No

Description

Whether ot not the XCFun library should be used.

xcfun
Type

Bool

Default value

No

Description

Functional for the XCFun library.

ZlmFit
Type

Block

Description

Options for the density fitting scheme ‘ZlmFit’.

AllowBoost
Type

Bool

Default value

Yes

Description

Allow automatic atom-dependent tuning of maximum l of spherical harmonics expansion. Whether or not this boost is needed for a given atom is based on an heuristic estimate of how complex the density around that atom is.

DensityThreshold
Type

Float

Description

Threshold below which the electron density is considered to be negligible. Depends on Quality and is normally 1.0e-7

FGaussianW
Type

Float

Default value

1.0

Description

Only for 3D periodic systems. Width of the Gaussian functions replacing the S and P Zlms for Fourier transform.

FGridSpacing
Type

Float

Description

Only for 3D periodic systems. Spacing for the Fourier grid. By default, this depends on the quality.

FKSpaceCutOff
Type

Float

Description

Only for 3D periodic systems. Cut-off of the grid in k-space for the Fourier transform.

FirstTopoCell
Type

Integer

Default value

5

Description

First cell for the topological extrapolation of the long range part of the Coulomb Potential.

LMargin
Type

Integer

Description

User-defined l-margin, i.e., l_max for fitting is max(lMargin + l_max_basis_function, 2*l_max_basis_function). Depends on Quality and normally is 4

LastTopoCell
Type

Integer

Default value

10

Description

Last cell for the topological extrapolation of the long range part of the Coulomb Potential.

NumStarsPartitionFun
Type

Integer

Default value

5

Description

Number of cell stars to consider when computing the partition function.

OrderTopoTrick
Type

Integer

Default value

3

Description

Order of the topological extrapolation of the long range part of the Coulomb Potential.

PartitionFunThreshold
Type

Float

Default value

0.0

Description

Threshold for the partition functions: if an integration point has a partition function weight smaller than this threshold, it will be discarded.

Quality
Type

Multiple Choice

Default value

Auto

Options

[Auto, Basic, Normal, Good, VeryGood, Excellent]

GUI name

Spline Zlm fit

Description

Quality of the density-fitting approximation. For a description of the various qualities and the associated numerical accuracy see reference. If ‘Auto’, the quality defined in the ‘NumericalQuality’ will be used.

QualityPerRegion
Type

Block

Recurring

True

Description

Sets the ZlmFit quality for all atoms in a region. If specified, this overwrites the globally set quality.

Quality
Type

Multiple Choice

Options

[Basic, Normal, Good, VeryGood, Excellent]

Description

The region’s quality of the ZlmFit.

Region
Type

String

Description

The identifier of the region for which to set the quality.