FLEUR

Introduction

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.

Environment variables

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.

FLEUR Calculator

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

Spin-polarized calculation

If the atoms object has non-zero magnetic moments, a spin-polarized calculation will be performed by default.

Utility functions

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.

FLEUR.write_inp(atoms)[source]

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.

FLEUR.initialize_density(atoms)[source]

Creates a new starting density.

FLEUR.calculate(atoms)[source]

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.

FLEUR.relax(atoms)[source]

Currently, user has to manually define relaxation parameters (atoms to relax, relaxation directions, etc.) in inp file before calling this function.

Examples

Lattice constant of fcc Ni

from numpy import linspace

from ase.calculators.fleur import FLEUR
from ase.build 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.