P

Functional : GGA X Becke C Lee/Yang/Parr

Projection : Mulliken dipole

Calculation the bulk modulus obtained by fitting to the Murnaghan equation. The energy is plotted relative to the minimum energy for each case, -177.58 eV/atom in Fireball and -179.099 eV/atom for Quantum espresso using the pseudopotentials P.pbe-van_ak.UPF . a0 and B0 are the equilibrium lattice constant and bulk modulus obtained by fitting to the Murnaghan equation of state.

bulk.png

Calculation the band dispersion for the bcc structure with the lattice constant of 4.5 Å

Fireball

bandas for Fireball

Quamtum espresso

bandas

The cut off radio of the orbitals in Fireball are rs = 5.5 (bohr), rp = 5.5 (bohr), rd = 5.5 (bohr), and was generate like qs = 2.0, qp = 2.2, qd = 0.0, . The orbital d had been generated from the ion with charge qs = 2.0, qp = 2.6, qd = 2.0.

wf.png
Other basis
Base Functional Projection method Crystal structure Lattice parameter (Å) Bulk modulus (GPa)
Ps200800p100800d000350 LDA Weighted LWeighted Lowdinouml;wdin bcc 3.056 132
Ps200800p100800d000350 LDA Löwdin bcc 3.060 132
Ps200800p100800d000350 LDA Mulliken dipole bcc 3.059 132
Ps050800p250800d000350 GGA Weighted LWeighted Lowdinouml;wdin bcc 3.057 127
Ps200550p220550d000550Id200260000 GGA Löwdin bcc 3.048 108
Ps000800p300800d000350 GGA Mulliken dipole bcc 3.049 135

The LDA calculation is done with the ncpp file to obtain 015.pp pseudopotential :

ncpp.psot ncpp.psot

The GGA calculation is done with the ncpp file to obtain 015.pp pseudopotential :

ncpp.psot ncpp.psot