SCF

CalcEigenvalues
Section: SCF
Type: logical
Default: yes

(Experimental) When this variable is set to no, Octopus will not calculate the eigenvalues or eigenvectors of the Hamiltonian. Instead, Octopus will obtain the occupied subspace. The advantage that calculation can be made faster by avoiding subspace diagonalization and other calculations.

This mode cannot be used with unoccupied states.


SCFCalculateDipole
Section: SCF
Type: logical

This variable controls whether the dipole is calculated at the end of a self-consistent iteration. For finite systems the default is yes. For periodic systems the default is no, unless an electric field is being applied in a periodic direction. The single-point Berry`s phase approximation is used for periodic directions. Ref: E Yaschenko, L Fu, L Resca, and R Resta, Phys. Rev. B 58, 1222-1229 (1998).


SCFCalculateForces
Section: SCF
Type: logical

This variable controls whether the forces on the ions are calculated at the end of a self-consistent iteration. The default is yes, unless the system only has user-defined species.


SCFCalculatePartialCharges
Section: SCF
Type: logical
Default: no

(Experimental) This variable controls whether partial charges are calculated at the end of a self-consistent iteration.


SCFCalculateStress
Section: SCF
Type: logical

This variable controls whether the stress on the lattice is calculated at the end of a self-consistent iteration. The default is no.


SCFinLCAO
Section: SCF
Type: logical
Default: no

Performs the SCF cycle with the calculation restricted to the LCAO subspace. This may be useful for systems with convergence problems (first do a calculation within the LCAO subspace, then restart from that point for an unrestricted calculation).


SCF::Convergence

ConvAbsDens
Section: SCF::Convergence
Type: float
Default: 0.0

Absolute convergence of the density:

\(\varepsilon = \int {\rm d}^3r \left| \rho^{out}(\bf r) -\rho^{inp}(\bf r) \right|\).

A zero value (the default) means do not use this criterion.


ConvAbsEv
Section: SCF::Convergence
Type: float
Default: 0.0

Absolute convergence of the sum of the eigenvalues:

\( \varepsilon = \left| \sum_{j=1}^{N_{occ}} \varepsilon_j^{out} - \sum_{j=1}^{N_{occ}} \varepsilon_j^{inp} \right| \)

A zero value (the default) means do not use this criterion.


ConvEigenError
Section: SCF::Convergence
Type: logical
Default: false

If true, the calculation will not be considered converged unless all states have individual errors less than EigensolverTolerance.


ConvEnergy
Section: SCF::Convergence
Type: float
Default: 0.0

Stop the SCF when the magnitude of change in energy during at one SCF iteration is smaller than this value.


ConvForce
Section: SCF::Convergence
Type: float

Absolute convergence of the forces: maximum variation of any component of the ionic forces in consecutive iterations. A zero value means do not use this criterion. The default is zero, except for geometry optimization, which sets a default of 1e-8 H/b.


ConvRelDens
Section: SCF::Convergence
Type: float
Default: 1e-5

Relative convergence of the density:

\(\varepsilon = \frac{1}{N} \mathrm{ConvAbsDens}\).

N is the total number of electrons in the problem. A zero value means do not use this criterion.

If you reduce this value, you should also reduce EigensolverTolerance to a value of roughly 1/10 of ConvRelDens to avoid convergence problems.


ConvRelEv
Section: SCF::Convergence
Type: float
Default: 0.0

Relative convergence of the sum of the eigenvalues:

\(\varepsilon = \frac{ \left| \sum_{j=1}^{N_{occ}} ( \varepsilon_j^{out} - \varepsilon_j^{inp} ) \right|} {\left| \sum_{j=1}^{N_{occ}} \varepsilon_j^{out} \right|} \)


MaximumIterBerry
Section: SCF::Convergence
Type: integer
Default: 10

Maximum number of iterations for the Berry potential, within each SCF iteration. Only applies if a StaticElectricField is applied in a periodic direction. The code will move on to the next SCF iteration even if convergence has not been achieved. -1 means unlimited.


SCF::Eigensolver

CGAdditionalTerms
Section: SCF::Eigensolver
Type: logical
Default: no

Used by the cg solver only. Add additional terms during the line minimization, see PTA92, eq. 5.31ff. These terms can improve convergence for some systems, but they are quite costly. If you experience convergence problems, you might try out this option. This feature is still experimental.


CGDirection
Section: SCF::Eigensolver
Type: integer

Used by the cg solver only. The conjugate direction is updated using a certain coefficient to the previous direction. This coeffiction can be computed in different ways. The default is to use Fletcher-Reeves (FR), an alternative is Polak-Ribiere (PR).
Options:


CGOrthogonalizeAll
Section: SCF::Eigensolver
Type: logical
Default: no

Used by the cg solver only. During the cg iterations, the current band can be orthogonalized against all other bands or only against the lower bands. Orthogonalizing against all other bands can improve convergence properties, whereas orthogonalizing against lower bands needs less operations.


Eigensolver
Section: SCF::Eigensolver
Type: integer

Which eigensolver to use to obtain the lowest eigenvalues and eigenfunctions of the Kohn-Sham Hamiltonian. The default is conjugate gradients (cg), except that when parallelization in states is enabled, the default is lobpcg.
Options:


EigensolverImaginaryTime
Section: SCF::Eigensolver
Type: float
Default: 10.0

The imaginary-time step that is used in the imaginary-time evolution method (Eigensolver = evolution) to obtain the lowest eigenvalues/eigenvectors. It must satisfy EigensolverImaginaryTime > 0.


EigensolverMaxIter
Section: SCF::Eigensolver
Type: integer

Determines the maximum number of iterations that the eigensolver will perform if the desired tolerance is not achieved. The default is 25 iterations for all eigensolvers except for rmdiis, which performs only 3 iterations (only increase it if you know what you are doing).


EigensolverMinimizationIter
Section: SCF::Eigensolver
Type: integer
Default: 5

During the first iterations, the RMMDIIS eigensolver requires some steepest-descent minimizations to improve convergence. This variable determines the number of those minimizations.


EigensolverSkipKpoints
Section: SCF::Eigensolver
Type: logical

Only solve Hamiltonian for k-points with zero weight


EigensolverTolerance
Section: SCF::Eigensolver
Type: float

This is the tolerance for the eigenvectors. The default is 1e-6, except for the ARPACK solver for which it is 0.


Preconditioner
Section: SCF::Eigensolver
Type: integer

Which preconditioner to use in order to solve the Kohn-Sham equations or the linear-response equations. The default is pre_filter, except for curvilinear coordinates, where no preconditioner is applied by default.
Options:


PreconditionerFilterFactor
Section: SCF::Eigensolver
Type: float

This variable controls how much filter preconditioner is applied. A value of 1.0 means no preconditioning, 0.5 is the standard.

The default is 0.5, except for periodic systems where the default is 0.6.

If you observe that the first eigenvectors are not converging properly, especially for periodic systems, you should increment this value.


PreconditionerIterationsMiddle
Section: SCF::Eigensolver
Type: integer

This variable is the number of smoothing iterations on the coarsest grid for the multigrid preconditioner. The default is 1.


PreconditionerIterationsPost
Section: SCF::Eigensolver
Type: integer

This variable is the number of post-smoothing iterations for the multigrid preconditioner. The default is 2.


PreconditionerIterationsPre
Section: SCF::Eigensolver
Type: integer

This variable is the number of pre-smoothing iterations for the multigrid preconditioner. The default is 1.


StatesOrthogonalization
Section: SCF::Eigensolver
Type: integer

The full orthogonalization method used by some eigensolvers. The default is cholesky_serial, except with state parallelization, the default is cholesky_parallel.
Options:


SubspaceDiagonalization
Section: SCF::Eigensolver
Type: integer
Default: standard

Selects the method to perform subspace diagonalization. The default is standard, unless states parallelization is used, when the default is scalapack.
Options:


SCF::LCAO

AtomsMagnetDirection
Section: SCF::LCAO
Type: block

This option is only used when GuessMagnetDensity is set to user_defined. It provides a direction for the magnetization vector of each atom when building the guess density. In order to do that, the user should specify the coordinates of a vector that has the desired direction and norm. Note that it is necessary to maintain the ordering in which the species were defined in the coordinates specifications.

For spin-polarized calculations, the vectors should have only one component; for non-collinear-spin calculations, they should have three components. If the norm of the vector is greater than the number of valence electrons in the atom, it will be rescaled to this number, which is the maximum possible magnetization.


GuessMagnetDensity
Section: SCF::LCAO
Type: integer
Default: ferromagnetic

The guess density for the SCF cycle is just the sum of all the atomic densities. When performing spin-polarized or non-collinear-spin calculations this option sets the guess magnetization density.

For anti-ferromagnetic configurations, the user_defined option should be used.

Note that if the paramagnetic option is used, the final ground state will also be paramagnetic, but the same is not true for the other options.
Options:


LCAOAlternative
Section: SCF::LCAO
Type: logical
Default: false

If this variable is set, the LCAO procedure will use an alternative (and experimental) implementation. It is faster for large systems and parallel in states. It is not working for spinors, however.


LCAOComplexYlms
Section: SCF::LCAO
Type: logical
Default: false

If set to true, and using complex states, complex spherical harmonics will be used, i.e. with \(e^{\pm i m \phi}\). If false, real spherical harmonics with \(\sin(m \phi)\) or \(\cos(m \phi)\) are used. This variable will make it more likely to get states that are eigenvectors of the \(L_z\) operator, with a definite angular momentum.


LCAODiagTol
Section: SCF::LCAO
Type: float
Default: 1e-10

Only applies if LCAOAlternative = true. The tolerance for the diagonalization of the LCAO Hamiltonian.


LCAODimension
Section: SCF::LCAO
Type: integer

(Only applies if LCAOAlternative = no.) Before starting the SCF cycle, an initial LCAO calculation can be performed in order to obtain reasonable initial guesses for spin-orbitals and densities. For this purpose, the code calculates a number of atomic orbitals. The number available for a species described by a pseudopotential is all the orbitals up the maximum angular momentum in the pseudopotential, minus any orbitals that are found to be unbound. For non-pseudopotential species, the number is equal to twice the valence charge. The default dimension for the LCAO basis set will be the sum of all these numbers, or twice the number of required orbitals for the full calculation, whichever is less.

This dimension however can be changed by making use of this variable. Note that LCAODimension cannot be smaller than the number of orbitals needed in the full calculation -- if LCAODimension is smaller, it will be silently increased to meet this requirement. In the same way, if LCAODimension is larger than the available number of atomic orbitals, it will be reduced. If you want to use the largest possible number, set LCAODimension to a negative number.


LCAOExtraOrbitals
Section: SCF::LCAO
Type: logical
Default: false

Only applies if LCAOAlternative = true, and all species are pseudopotentials. (experimental) If this variable is set to yes, the LCAO procedure will add an extra set of numerical orbitals (by using the derivative of the radial part of the original orbitals). Note that this corresponds roughly to adding orbitals with higher principal quantum numbers, but the same angular momentum. This option may cause problems for unoccupied states since you may miss some lower-lying states which correspond to higher angular momenta instead of higher principal quantum number.


LCAOKeepOrbitals
Section: SCF::LCAO
Type: logical
Default: yes

Only applies if LCAOAlternative = true. If set to yes (the default) Octopus keeps atomic orbitals in memory during the LCAO procedure. If set to no, the orbitals are generated each time that they are needed, increasing computational time but saving memory.

When set to yes, Octopus prints the amount of memory per node that is required to store the orbitals.


LCAOMaximumOrbitalRadius
Section: SCF::LCAO
Type: float
Default: 20.0 a.u.

The LCAO procedure will ignore orbitals that have an extent greater that this value.


LCAOScaleFactor
Section: SCF::LCAO
Type: float
Default: 1.0

The coordinates of the atomic orbitals used by the LCAO procedure will be rescaled by the value of this variable. 1.0 means no rescaling.


LCAOStart
Section: SCF::LCAO
Type: integer

Before starting a SCF calculation, Octopus can perform a linear combination of atomic orbitals (LCAO) calculation. These can provide Octopus with a good set of initial wavefunctions and with a new guess for the density. (Up to the current version, only a minimal basis set is used.) The default is lcao_full if at least one species representing an atom is present. The default is lcao_none if all species are species_user_defined, species_charge_density, species_from_file, or species_jellium_slab.

The initial guess densities for LCAO are taken from the atomic orbitals for pseudopotential species; from the natural charge density for species_charge_density, species_point, species_jellium, and species_jellium_slab; or uniform for species_full_delta, species_full_gaussian, species_user_defined, or species_from_file. Pseudopotential species use the pseudo-wavefunctions as orbitals, full-potential atomic species (species_full_delta and species_full_gaussian) use hydrogenic wavefunctions, and others use harmonic-oscillator wavefunctions.

Note: Some pseudopotential files (CPI, FHI for example) do not contain full information about the orbitals. In this case, Octopus generates the starting density from the normalized square root of the local potential. If no orbitals are available at all from the pseudopotential files, Octopus will not be able to perform an LCAO and the initial states will be randomized.


Options:


SCF::Mixing

MixField
Section: SCF::Mixing
Type: integer

Selects what should be mixed during the SCF cycle. Note that currently the exact-exchange part of hybrid functionals is not mixed at all, which would require wavefunction-mixing, not yet implemented. This may lead to instabilities in the SCF cycle, so starting from a converged LDA/GGA calculation is recommended for hybrid functionals. The default depends on the TheoryLevel and the exchange-correlation potential used.
Options:


MixInterval
Section: SCF::Mixing
Type: integer
Default: 1

When this variable is set to a value different than 1 (the defaul) a combined mixing scheme will be used, with MixInterval - 1 steps of linear mixing followed by 1 step of the selected mixing. For the moment this variable only works with DIIS mixing.


MixNumberSteps
Section: SCF::Mixing
Type: integer
Default: 3

In the Broyden and Bowler_Gillan schemes, the new input density or potential is constructed from the values of the densities/potentials of a given number of previous iterations. This number is set by this variable. Must be greater than 1.


Mixing
Section: SCF::Mixing
Type: float
Default: 0.3

The linear, Broyden and DIIS scheme depend on a "mixing parameter", set by this variable. Must be 0 < Mixing <= 1.


MixingPreconditioner
Section: SCF::Mixing
Type: logical
Default: false

(Experimental) If set to yes, Octopus will use a preconditioner for the mixing operator.


MixingResidual
Section: SCF::Mixing
Type: float
Default: 0.05

In the DIIS mixing it is benefitial to include a bit of residual into the mixing. This parameter controls this amount.


MixingScheme
Section: SCF::Mixing
Type: integer
Default: broyden

The scheme used to produce, at each iteration in the self-consistent cycle that attempts to solve the Kohn-Sham equations, the input density from the value of the input and output densities of previous iterations.
Options:


SCF::RDMFT

RDMBasis
Section: SCF::RDMFT
Type: logical
Default: yes

If true, all the energy terms and corresponding derivatives involved in RDMFT will not be calculated on the grid but on the basis of the initial orbitals


RDMConvEner
Section: SCF::RDMFT
Type: float
Default: 1e-6 Ha

Convergence criterion for stopping the overall minimization of the energy with respect to occupation numbers and the orbitals. The minimization of the energy stops when the total energy difference between two subsequent minimizations of the energy with respect to the occupation numbers and the orbitals is smaller than this criterion. It is also used to exit the orbital minimization.


RDMTolerance
Section: SCF::RDMFT
Type: float
Default: 1e-7 Ha

Convergence criterion for stopping the occupation numbers minimization. Minimization is stopped when all derivatives of the energy wrt. each occupation number are smaller than this criterion. The bisection for finding the correct mu that is needed for the occupation number minimization also stops according to this criterion.