FLEUR is a density-functional theory code which uses the full potential linearized augmented plane-wave (FLAPW) method. FLEUR can be applied to any element in the periodic table, and the code is well suited especially to surfaces and magnetic materials.
In order to use FLEUR through ASE, two environment variables have to be set. FLEUR_INPGEN should point to the simple input generator of FLEUR, and FLEUR to the actual executable. Note that FLEUR has different executables e.g. for cases with and without inversion symmetry, so the environment variable has to be set accordingly. As an example, the variables could be set like:
$ export FLEUR_INPGEN=$HOME/fleur/v25/inpgen.x
$ export FLEUR=$HOME/fleur/v25/fleur.x
or:
$ export FLEUR="mpirun -np 32 $HOME/fleur/v25/fleur.x"
for parallel calculations.
Currently, limited number of FLEUR parameters can be set via the ASE interface Below follows a list of supported parameters
keyword | type | default value | description |
---|---|---|---|
xc | str | 'LDA' | XC-functional. Must be one of ‘LDA’, ‘PBE’, ‘RPBE’ |
kpts | seq | \(\Gamma\)-point | k-point sampling |
convergence | dict | {'energy': 0.0001} | Convergence criteria (meV) |
width | float | Width of Fermi smearing (eV) | |
kmax | float | Plane-wave cut-off (a.u.) | |
mixer | dict | Mixing parameters ‘imix’, ‘alpha’, and ‘spinf’ | |
maxiter | int | 40 | Maximum number of SCF steps |
maxrelax | int | 20 | Maximum number of relaxation steps |
workdir | str | Current dir | Working directory for the calculation |
seq: A sequence of three int‘s. dict: A dictionary
If the atoms object has non-zero magnetic moments, a spin-polarized calculation will be performed by default.
As only a subset of FLEUR parameters can currently be specified through ASE interface, the interface defines some utility functions for cases where manual editing of the FLEUR input file inp is necessary.
Write the inp input file of FLEUR.
First, the information from Atoms is written to the simple input file and the actual input file inp is then generated with the FLEUR input generator. The location of input generator is specified in the environment variable FLEUR_INPGEN.
Finally, the inp file is modified according to the arguments of the FLEUR calculator object.
Converge a FLEUR calculation to self-consistency.
Input files should be generated before calling this function FLEUR performs always fixed number of SCF steps. This function reduces the number of iterations gradually, however, a minimum of five SCF steps is always performed.
Lattice constant of fcc Ni
from numpy import linspace
from ase.calculators.fleur import FLEUR
from ase.lattice import bulk
from ase.io.trajectory import Trajectory
atoms = bulk('Ni', a=3.52)
calc = FLEUR(xc='PBE', kmax=3.6, kpts=(10, 10, 10), workdir='lat_const')
atoms.set_calculator(calc)
traj = Trajectory('Ni.traj','w', atoms)
cell0 = atoms.get_cell()
for s in linspace(0.95, 1.05, 7):
cell = cell0 * s
atoms.set_cell((cell))
ene = atoms.get_potential_energy()
traj.write()
See the equation of states tutorial for analysis of the results.