Chapter 6
Keyword Reference

The cell and param file keywords used specifically for phonon and related calculations are listed here alphabetically with brief descriptions.

6.1 CELL file keywords

PHONON_FINE_KPOINT_LIST List of phonon wavevectors on the fine grid (Block)
Phonon frequencies are calculated on a coarse set of wavevectors using DFPT and interpolated onto this finer list of points.
PHONON_FINE_KPOINT_MP_GRID Fine MP-grid of phonon wavevectors (Integer Vector)
Phonon frequencies are calculated on a coarse set of wavevectors using DFPT and interpolated onto this finer grid of wavevectors.
PHONON_FINE_KPOINT_MP_OFFSET Origin offset of the fine phonon MP grid (Real Vector)
The offset of the fine MP grid at which the phonons calculated using DFPT are interpolated.
PHONON_FINE_KPOINT_MP_SPACING The spacing of points on the fine MP set for phonons (Physical)
This specifies the minimum spacing between points on a Monkhorst-Pack grid that phonons will be interpolated onto from the coarser phonon grid.
PHONON_FINE_KPOINT_PATH Path of phonon wavevectors on a fine scale (Block)
Phonon frequencies are calculated on a coarse set of wavevectors using DFPT and interpolated onto this finer path.
PHONON_FINE_KPOINT_PATH_SPACING The fine spacing of points on a path at which phonons are calculated (Physical)
The spacing of k-points along a path (specified by PHONON_FINE_KPOINT_PATH) at which phonons will be interpolated from a coarser grid
PHONON_GAMMA_DIRECTIONS Phonon gamma-point LO/TO splitting (Block)
This is a list of directions along which q 0 will be calculated for the non-analytic LO/TO term in a phonon calculation at q = 0. Fractional coordinates must be used. Default value: The k-point before gamma in the k-point list, or the one after or (0.1, 0, 0)
PHONON_KPOINTS_LIST Alias for PHONON_KPOINT_LIST (Block)
Default value: Determined from PHONON_KPOINT_LIST
PHONON_KPOINTS_PATH Alias for PHONON_KPOINT_PATH (Block)
Default value: Determined from PHONON_KPOINT_PATH
PHONON_KPOINTS_PATH_SPACING Alias for PHONON_KPOINT_PATH_SPACING (Physical)
Default value: Determined from PHONON_KPOINT_PATH_SPACING
PHONON_KPOINT_LIST List of phonon k-points (Block)
A list of discrete k-points at which phonon frequencies and eigenvectors will be calculated. Default value: SCF k-points are used if an alternative PHONON_KPOINT_* specifier is not given.
PHONON_KPOINT_MP_GRID Phonon wavevector Monkhorst-Pack grid (Integer Vector)
The phonon wavevectors defined by a Monkhorst-Pack grid. Symmetry (if specified) will be used to generate the wavevector list and weights.
PHONON_KPOINT_MP_OFFSET Phonon wavevector Monkhorst Pack grid offset (Real Vector)
The offset of the origin of the Monkhorst-Pack set for phonons in fractional coordinates Default value: (0, 0, 0).
PHONON_KPOINT_MP_SPACING Phonon wavevector Monkhorst-Pack grid density (Physical)
The density of wavevectors on a a Monkhorst-Pack grid for phonon calculations. Units of inverse length should be specified. Default value: 0.1 A-1.
PHONON_KPOINT_PATH Phonon dispersion k-point path (Block)
The path continuous through the BZ on which phonon dispersion is calculated. This is specified in fractional coordinates. Default value: None.
PHONON_KPOINT_PATH_SPACING Phonon dispersion path spacing (Physical)
The maximum spacing between kpoints along the path specified by PHONON_KPOINT_PATH. Units of inverse length must be specified. Default value: 0.1 A-1.
PHONON_SUPERCELL_MATRIX Supercell matrix for finite difference phonon calculations (Block)
The supercelling matrix for force constant matrix calculations. The supercell matrix is specified by a 3x3 integer matrix which gives the supercell used in finite-difference phonon calculations.
SUPERCELL_KPOINTS_LIST SCF k-points for FD phonon supercell (Block)
A list of k-points in the Brillouin zone (with associated weights) used for BZ integration during a supercell FD phonon calculation. The k-point weights must sum to 1. Default value: Generated from SUPERCELL_KPOINTS_MP_SPACING and the crystal symmetry.
SUPERCELL_KPOINTS_MP_GRID SCF Monkhorst-Pack grid for FD phonon supercell calculation (Integer Vector)
The k-points defined by a Monkhorst-Pack grid when doing a finite displacement phonon calculation. Symmetry (if specified) will be used to generate the k-point list and weights. Default value: Generated from SUPERCELL_KPOINTS_MP_SPACING.
SUPERCELL_KPOINTS_MP_OFFSET SCF Monkhorst Pack grid offset for a FD phonon supercell calculation (Real Vector)
The offset of the origin of the Monkhorst-Pack set in fractional coordinates when performing a finite displacement phonon calculation. Default value: (0, 0, 0).
SUPERCELL_KPOINTS_MP_SPACING SCF Monkhorst-Pack grid density for FD phonon supercell calculation (Physical)
The k-point density of a Monkhorst-Pack grid for a supercell FD phonon calculation. Units of inverse length should be specified. Default value: 0.1 A-1.

6.2 PARAM file keywords

NUM_BACKUP_ITER md/geom iterations between backups (Integer)
Specifies the number of iterations between backups of all data for restarts, for a geometry optimization or molecular dynamics run. Allowed values: (any integer) > 0 Default value : 5
BACKUP_INTERVAL seconds between backups (Integer)
Specifies the interval, in seconds, between backups of all data for restarts, for a geometry optimization/molecular dynamics/phonon run - if less than or equal to zero then no timed backups. Allowed values: (any) Default value : 0
BORN_CHARGE_SUM_RULE enforce Born charge sum rule (Logical)
Selects whether to explicitly correct the Born effective charge tensor to enforce the sum rule that effective charges sum to zero.
CALCULATE_BORN_CHARGES calculate Born effective charges (Logical)
Selects whether to compute Born effective charge tensors as part of a phonon or efield linear-response calculation. Allowed values: TRUE or FALSE Default value : TRUE
CALCULATE_BORN_CHARGES calculate Born effective charges (Logical)
Selects whether to compute Born effective charge tensors as part of a phonon or efield linear-response calculation. Allowed values: TRUE or FALSE Default value : TRUE
CALCULATE_BORN_CHARGES calculate Born effective charges (Logical)
Selects whether to compute Born effective charge tensors as part of a phonon or efield linear-response calculation. Allowed values: TRUE or FALSE Default value : TRUE
CALCULATE_RAMAN calculate Raman intensities (Logical)
Selects whether to compute Raman intensities as part of a phonon or efield linear-response calculation. Allowed values: TRUE or FALSE Default value : FALSE
CALCULATE_RAMAN calculate Raman intensities (Logical)
Selects whether to compute Raman intensities as part of a phonon or efield linear-response calculation. Allowed values: TRUE or FALSE Default value : FALSE
EFIELD_CALC_ION_PERMITTIVITY calculate zero-frequency permittivity (Logical)
Specifies whether or not to compute the zero-frequency dielectric permittivity based on the ionic response to electric fields. This requires a gamma-point phonon calculation in addition to the EFIELD linear response one. Allowed values: TRUE or FALSE Default value : TRUE
EFIELD_CALC_ION_PERMITTIVITY calculate zero-frequency permittivity (Logical)
Specifies whether or not to compute the zero-frequency dielectric permittivity based on the ionic response to electric fields. This requires a gamma-point phonon calculation in addition to the EFIELD linear response one. Allowed values: TRUE or FALSE Default value : TRUE
EFIELD_CONVERGENCE_WIN convergence tolerance window in EFIELD (Integer)
The LR convergence criteria must be met for EFIELD_CONVERGENCE_WIN iterations before acceptance. Allowed values: (any integer) 2 Default value : 2
EFIELD_ENERGY_TOL E(2) convergence tolerance in EFIELD (Physical)
Tolerance for accepting convergence of the field constants during PHONON calculation. The difference between max and min E(2) values over EFIELD_CONVERGENCE_WIN iterations must be less than this. NB This is an INTENSIVE parameter and has units of volume. Allowed values: (any) > 0.0 Default value : 10-5 A3
EFIELD_FREQ_SPACING Spacing of frequencies in permittivity calculation (Physical)
Spacing of frequencies in calculation of frequency-dependent permittivity. Allowed values: (any) > 0.0 Default value : 1.0 cm-1
EFIELD_IGNORE_MOLEC_MODES Ignore lowest modes in permittivity calculation (String)
Ignore the lowest lying (3,5,6) modes when computing the ionic contribution to the permittivity and polarizability. Allowed values: CRYSTAL, MOLECULE, LINEAR_MOLECULE Default value : CRYSTAL
EFIELD_MAX_CG_STEPS max. number of cg steps in EFIELD (Integer)
The maximum number of conjugate gradient steps in EFIELD calculation before performing a SD reset. Allowed values: (any integer) 0 Default value : 0
EFIELD_MAX_CYCLES maximum cycles in EFIELD (Integer)
The maximum number of SCF cycles in EFIELD calculation regardless of convergence. Allowed values: (any integer) > 0 Default value : 50
EFIELD_OSCILLATOR_Q Q-factor for line-shape broadening (Real)
Oscillator Q-factor for line-shape broadening in calculation of frequency- dependent permittivity. Allowed values: (any) > 0.0 Default value : 50.0
EFIELD_UNIT unit of electric field in output (String)
Controls the units used for all electric field in output - many different units are supported. Default value : eV/A/E
ELEC_METHOD treatment of metals or finite temperature insulator (String)
The treatment of metals or finite temperature insulator to be used. An alias for METALS_METHOD. Allowed values: NONE (=ALLBANDS), DM, EDFT Default value : DM
EXCITED_STATE_SCISSORS “scissors” operator band-gap correction (Physical)
Effectively adds an offset to conduction-band eigenvalues as empirical correction for LDA/GGA underestimation of band-gaps. Allowed values: (any) Default value : 0.0
FIX_OCCUPANCY treat system as an insulator (Logical)
Determines if the system is treated as an insulator or a metal. Allowed values: TRUE or FALSE Default value : FALSE
GEOM_FORCE_TOL geometry optimization force convergence tolerance (Physical)
Tolerance for accepting convergence of the maximum —ionic force— during geometry optimization. Allowed values: (any) > 0.0 Default value : 0.05 eV/A
PHONON_CALC_LO_TO_SPLITTING gamma-point phonon LO/TO correction (Logical)
Selects whether to compute non-analytic contribution to dynamical matrix from long-ranged electric field effects responsible for LO/TO splitting. This requires calculation of the dielectric permittivity by efield linear-response and the Born effective charges. Allowed values: TRUE or FALSE Default value : TRUE
PHONON_CALC_LO_TO_SPLITTING gamma-point phonon LO/TO correction (Logical)
Selects whether to compute non-analytic contribution to dynamical matrix from long-ranged electric field effects responsible for LO/TO splitting. This requires calculation of the dielectric permittivity by efield linear-response and the Born effective charges. Allowed values: TRUE or FALSE Default value : TRUE
PHONON_CONVERGENCE_WIN convergence tolerance window in LR (Integer)
The LR convergence criteria must be met for PHONON_CONVERGENCE_WIN iterations before acceptance. Allowed values: (any integer) 2 Default value : 2
PHONON_ENERGY_TOL E(2) convergence tolerance in LR (Physical)
Tolerance for accepting convergence of the force constants during PHONON calculation. The difference between max and min E(2) values over PHONON_CONVERGENCE_WIN iterations must be less than this. Allowed values: (any) > 0.0 Default value : 10-5 eV/A2
PHONON_FINE_METHOD fine phonon calculation method (String)
Selects which calculation method to use for phonon calculation on a fine grid. Allowed values: NONE, SUPERCELL, INTERPOLATE Default value : SUPERCELL if TASK=THERMODYNAMICS else NONE
PHONON_FINITE_DISP finite displacement amplitude (Physical)
The amplitude of the ionic perturbation to be used in a finite displacement phonon calculation. Allowed values: (any) > 0.0 Default value : 0.01 a0
PHONON_FORCE_CONSTANT_CUTOFF Cutoff for force constant matrix (Physical)
The cutoff for the force constant matrix in a phonon calculation on a fine grid with supercell method. Allowed values: (any) 0.0 Default value : 0.0
PHONON_FINE_CUTOFF_METHOD Selects which method to use to extract non-periodic force constant matrix from periodic supercell. (String)
With the CUMULANT method, all contributions from the periodic supercell are summed with a suitable weighting factor to avoid double counting of image contributions.
The SPHERICAL method, uses a minimum image convention with a spherical cutoff given by PHONON_FORCE_CONSTANT_CUTOFF.
Allowed values: CUMULANT and SPHERICAL. Default value : CUMULANT.
PHONON_FORCE_CONSTANT_CUT_SCALE Scaling factor for aspherical force constant matrix cutoff (Real)
The range of force constant terms included is up to s times halfway to the Wigner Seitz cell boundary. This parameter supplies the value of s. Allowed values: 0.0 (any) 1.0 Default value : 0.0
PHONON_FORCE_CONSTANT_ELLIPSOID Ellipsoid size for force constant matrix (Real)
Alias for PHONON_FORCE_CONSTANT_CUT_SCALE (deprecated).
PHONON_MAX_CG_STEPS max. number of cg steps in LR (Integer)
The maximum number of conjugate gradient steps in PHONON calculation before performing a SD reset. Allowed values: (any integer) 0 Default value : 0
PHONON_MAX_CYCLES maximum cycles in LR (Integer)
The maximum number of SCF cycles in PHONON calculation regardless of convergence. Allowed values: (any integer) 0 or if TASK=PHONON etc PHONON_CONVERGENCE_WIN Default value : 50
PHONON_METHOD phonon calculation method (String)
Selects which calculation method to use for phonons. Allowed values: DFPT, LINEARRESPONSE, FINITEDISPLACEMENT Default value : set by PHONON_FINE_METHOD
PHONON_PRECONDITIONER scheme to use in LR (String)
The preconditioning scheme used by the CG minimiser in LR. Allowed values: RTPA, TPA, PS, NONE Default value : TPA
PHONON_SUM_RULE enforce phonon sum rule (Logical)
Selects whether to explicitly correct the dynamical matrix to enforce the acoustic q=0 phonon sum rule, i.e. that 3 modes have zero frequency at q=0. Allowed values: TRUE or FALSE Default value : FALSE
PHONON_SUM_RULE_METHOD select method to enforce phonon sum rule (String)
Selects a method to use when enforcing acoustic phonon sum rule. Allowed values: NONE : No sum-rule correction will be applied. RECIPROCAL : Correct dynamical matrix D(q) at each q using D(q=0). REALSPACE : Correct the real-space force constant matrix C(R). REAL-RECIP : Correct C(R) in realspace followed by D(q) in reciprocal space. MOLECULAR : Correct D(0) using rotational as well as translational sum-rule. MOLECULAR-1D : Correct D(0) for a linear molecule using rotational as well as translational sum-rule. Default value : RECIPROCAL
PHONON_USE_KPOINT_SYMMETRY reduced or full kpoint set in LR (Logical)
Selects which k-point set to use For each phonon q-vector in LR: T = use the irreducible k-point set of the (reduced) symmetry, F = use the complete fully symmetric k-point set. Allowed values: TRUE or FALSE Default value : TRUE
THERMO_T_NPOINTS Number of points in temperature interval (Integer)
The number of points in the temperature interval for the thermodynamics calculation. Allowed values: (any integer) 1 Default value : 2 if THERMO_T_STOP /= THERMO_T_START or 1 otherwise
THERMO_T_SPACING Temperature spacing (Physical)
The spacing between temperature values for the thermodynamics calculation. Allowed values: (any) > (-epsilon) Default value : THERMO_T_STOP - THERMO_T_START
THERMO_T_START Starting temperature (Physical)
The desired starting temperature for the thermodynamics calculation. Allowed values: (any) > 0.0 Default value : 298 K
THERMO_T_STOP Final temperature (Physical)
The desired final temperature for the thermodynamics calculation. Allowed values: (any) THERMO_T_START Default value : 298